Open access peer-reviewed chapter

Application of Spin-Orbit Coupling in Exotic Graphene Structures and Biology

Written By

Richard Pinčák and Erik Bartoš

Submitted: February 19th, 2019 Reviewed: July 9th, 2019 Published: February 13th, 2020

DOI: 10.5772/intechopen.88486

Chapter metrics overview

493 Chapter Downloads

View Full Metrics


An important measurable quantity in the carbon nanostructures, including the nanotubular part of the graphitic wormhole, is the spin-orbit coupling. We will present in this chapter spin-orbit coupling for the fermions located in exotic graphene structures as is graphene wormhole and also in biological systems. Considering this influence, the two-component Dirac equation is changed into the usual four-component form. As a consequence, the chiral fermions should be detected close to the wormhole bridge. We will show that the smaller is the radius of the wormhole bridge, the stronger this effect should be. Finally, we will describe the role of spinor fields in the time series of genetic code. The reversed transcription process of the gene expression could be defined by a moduli state space model of a coupling spinor field between the gene of a viral particle and the host cell. As a general result, all states of codon can be computed by the Chern-Simons 3-forms.


  • spinor network structure
  • spin orbit coupling
  • Chern-Simons fields
  • graphene wormhole
  • genetic code

1. Graphitic wormhole

The investigation of unique chemical and mechanical properties of nanostructures, e.g., fullerene, graphene, and nanotubes, promises a wide application in many technical areas. The electronic properties of the nanostructures are basically defined by their hexagonal carbon lattice structure and its variations. New promising results are expected with the preparation of more complicated forms as a wormhole. The wormhole is usually composed of two different kinds of nanostructure: two graphene sheets are connected together with the help of a connecting nanotube [1] (Figure 1). This is achieved by a supply of two sets of six heptagonal defects onto both sides of the given nanotube. There exists the restrictions on the form of the nanotube—the chirality must be 6 n 6 n armchair or 6 n 0 zigzag and a radius of the nanotube is larger than its length.

Figure 1.

Schematic representation of graphitic wormhole consisting from two graphene sheets connected together with the help of a nanotube.

The metric tensor of the wormhole is given by

g μν = Λ 2 r ± 1 0 0 r ± 2 , Λ r ± = a / r ± 2 θ a r ± + θ r ± a , E1

where θ is the Heaviside step function, r and r + are the polar coordinates of lower and upper graphene sheets, and a = r r + is the radius of the wormhole.

1.1 Electronic structure

We consider the continuum gauge field theory, i.e., at each point of a molecular surface, we take into account an influence of different gauge fields that enter the Dirac-like equation for an electron

i v F σ μ μ + Ω μ ia μ ia μ W iA μ ψ = , E2

with σ α matrices as the Pauli matrices, the Fermi velocity v F , the spin connection

Ω μ = 1 8 ω μ αβ σ α σ β , E3

and the covariant derivative μ = μ + Ω μ . The gauge fields a μ , a μ W are caused by the presence of the defects, and by rotational symmetry, the gauge field A μ characterizes the possible magnetic field.

In the case of the wormhole with the metric Eq. (1), the effective flux caused by the presence of the defects is included in the gauge field a μ , and for the particular polar components, it has the values a φ = 3 2 , a r = 0 for two possibilities: the first corresponds to the chiral vector with the form 6 n 6 n , the second corresponds to the chiral vector with the form 6 n 0 and n divisible by 3. In the case of chiral vector of the form 6 n 0 , where n is not divisible by 3, the components of the gauge field are a φ = 1 2 , a r = 0 . Regarding that the components of the spin connection are

Ω φ = i 2 σ 3 r Λ r Λ r + 1 , Ω r = 0 , E4

and after the substitution into Eq. (2), we get the equation

i v F σ μ μ + Ω μ i a μ ψ ± = εψ ± , E5

where each sign corresponds to a different Dirac point

i v F r + 1 r i θ a φ r + 1 2 r ψ B ± = εψ A ± , i v F r 1 r i θ ± a φ r + 1 2 r ψ A ± = εψ B ± E6

for r a and

i v F r a 2 r 1 r i θ ± a φ r 1 2 r ψ B ± = εψ A ± , i v F r a 2 r + 1 r i θ a φ r 1 2 r ψ A ± = εψ B ± E7

for 0 < r a . For r a , the solution is

ψ ± = ψ A ± r φ ψ B ± r φ = c 1 J n a φ 1 / 2 kr i sgn ε J n a φ + 1 / 2 kr + c 2 Y n a φ 1 / 2 kr i sgn ε Y n a φ + 1 / 2 kr , E8

where the energy ε = ± v F k , J n x and Y n x are the Bessel functions.

The zero modes solve the Dirac equation for zero energy. If one choose the component ψ A ± of the solution to be equal to zero, one get from (6) and (7)

r 1 r i θ a φ r + 1 2 r ψ B ± = 0 E9

for r a and

r 1 r i θ ± a φ r 1 2 r ψ B ± = 0 E10

for 0 < r a . For ψ B and the value a φ = 3 2 , the solution is

ψ B r φ r n 2 e i E11

for r a and

ψ B r φ r n + 2 e i E12

for 0 < r a . For both cases, it is normalizable only for n = 0 , and so this is the only solution. In a similar way, we can calculate the zero modes for the component ψ B + . For the value a φ = 1 2 , possible solutions are not strictly normalizable, and the zero modes exist only for the case of the connecting nanotube being armchair or zigzag with the chiral vector 6 n 0 , n divisible by 3 . In other cases the zero modes do not exist.

Recently, in work [2] some peculiarities in the bilayer graphene were analytically predicted. A possible indication of the wormhole could be found in [3, 4], where a new type of zero modes is investigated. These zero modes could be the zero modes studied in this subsection applied to the case of the smallest wormhole.

1.2 Case of massive fermions

Up to now we supposed that the fermions appearing in the Dirac equation have the zero mass or that the mass is very small in comparison with their energy, but in [5, 6] it was shown that the Fermi velocity needs to be renormalized due to the elasticity and deformations in a graphene. In our case of the graphitic wormhole, including big deformations, the velocity of fermions close to the wormhole bridge could achieve such values that the relativistic effects can appear or break off the symmetry [7] and the mass of fermions would be non-negligible. The radius of the wormhole and its bridge is very small in comparison with the size of the upper and the lower graphene sheet (Figure 2) and by folding the sheet into a tube they acquire nonzero effective mass as they move along the tube axis. This change of the space topology of graphene from 2D to 1D space compactification is similar to the string theory compactification, and we can imagine a wormhole connecting nanotubes as 1D object.

Figure 2.

The simplest realization of smallest graphitic wormholes.

To include the mass into the Dirac Eq. (2), one can transform the system of equations [8] into the differential equation of the second order

ξξ 1 2 g ξξ ξ g ξξ + j ˜ 2 g ξξ g φφ 3 ξ g φφ j ˜ 2 g ξξ g φφ + E 2 g ξξ u j = 0 . E13

One can suppose cylindrical geometry in order to simplify the equation into

ξξ + E 2 j ˜ 2 R 2 u j = 0 , E14

if a radius vector of the point at the surface will have the form

R = R cos φ R sin φ ξ , E15

with R as the radius of the cylinder. The solution of this equation has the form

u j ξ = Ae + Be , k = j ˜ 2 R 2 E 2 . E16

A similar form has the dispersion relation associated with the massive 1D Dirac equation

k = M 2 E 2 , E17

where M is the mass of corresponding fermion. It is proven [9] that for a suitable choice of the parameters, 2D massless case is in analogy with 1D massive case, and one can rewrite Eq. (13) in the form

ξξ 1 2 g ξξ ξ g ξξ + j ˜ 2 g ξξ g φφ 3 ξ g φφ j ˜ 2 g ξξ g φφ + E 2 M 2 g ξξ u j = 0 , E18

where M is the mass of the corresponding fermion. For different values of M , one can find the corrections of local density of states for the graphitic wormhole. It seems that these massive particles arising in the wormhole nanotubes could create energy bulks on and near the wormhole bridge which should be experimentally measured by the STM or Raman spectroscopy [10]. Another possibly identification of wormhole structure comes from the fact of creation of strain solitons and topological defects by massive particles on the bridge of bilayer graphene which should propagate throughout the graphene sheet. These are almost macroscopic effects and should be caught by the experimental physicists [11].

1.3 Spin-orbit coupling in the wormhole connecting nanotube

An important measurable quantity in carbon nanostructures, which includes a nanotubular part of a graphitic wormhole, is a spin-orbit coupling (SOC) [12, 13]. If one considers this influence, two-component Dirac equation could be changed into the usual four-component form, and as a consequence chiral fermions should be detected close to the wormhole bridge.

One can reflect on two sources of SOC: (1) the interatomic one that preserves the z -component of a spin and (2) so-called Rashba-type coming from the external electric field, which conserves the z -component of an angular momentum J z . In both cases, the strength of SOC is influenced by the nonzero curvature. In the next we will be interested in the first source of the SOC.

Considering the SOC we can write the Dirac equation for the nanotube in the form

H ̂ F A K F B K = 0 f ̂ f ̂ 0 F A K F B K = E F A K F B K , E19


F A K = F A , K F A , K , F B K = F B , K F B , K . E20

The expression f ̂ has the following form

f ̂ = γ k ̂ x i k ̂ y + i δ γ 4 R σ ̂ x r 2 δγp R σ ̂ y , E21


k ̂ x = i R θ , k ̂ y = i y , σ ̂ x r = σ ̂ x cos θ σ ̂ z sin θ . E22

Next one can take

γ = 3 2 aV pp π , γ = 3 2 a V pp σ V pp π , p = 1 3 γ 8 γ , E23

where a is the length of the atomic bond and V pp σ , V pp π are the hopping integrals for the σ and π bond, respectively.

For the interatomic source of the SOC, one has

δ = Δ 3 ε πσ , Δ = i 3 4 m 2 c 2 x l V x p ̂ y V y p ̂ x y l E24

with the difference of energies of the relevant π and σ orbitals

ϵ πσ = ϵ 2 p π ϵ 2 p σ , E25

x l , and y l being the local coordinates. By applying the transformation

H ̂ ' = U ̂ H ̂ U ̂ 1 , U ̂ = exp i σ ̂ y θ 2 0 0 exp i σ ̂ y θ 2 E26

the transformed Hamiltonian H ̂ will have the form with two terms, including the H ̂ SOC term which corresponds to the spin-orbit coupling

H ̂ = H ̂ kin + H ̂ SOC , H ̂ kin = y Id 2 s ̂ y + 1 R θ Id 2 s ̂ x , H ̂ SOC = λ y σ ̂ x s ̂ y λ x σ ̂ y s ̂ x . E27

The operators s ̂ x , y , z are the Pauli matrices, which transform the wave function of the A sublattice into the wave function of the B sublattice and vice versa.

In our model, the SOC is induced by the curvature, and it is described with the help of two strength parameters, namely, λ x and λ y in the form

λ x = γ R 1 2 + 2 δp , λ y = δ γ 4 R , E28

for the case of single-wall carbon nanotube with different magnitude. Here, λ y λ x and for R 0 , both strengths go to infinity. So reminding the previous results, the chiral massive fermions should be detected around the wormhole bridge. For more complicated forms as perturbed nanotube in the wormhole center, the geometry of the corresponding graphene sheets will be curved, and this brings a significant change of the physical properties.

1.4 Graphene black hole

The effects connected with the deformation of a graphene and a consequent change of the distance of the carbon atoms in the layer are described in [14]. It causes the rotation of the p z orbitals and rehybridization of the π and σ orbitals. The procedure leads to the creation of the p n junctions similarly to the case of a transistor. This effect changes the Fermi level which is rising in the far areas from the wormhole center. The electron flux is directed from these areas to the middle where the electric charge is accumulated, and in the case of the deformed wormhole, one can speak about so-called graphene black hole. The form of a middle part of the nanotube plays a big role for this purpose. It cannot be unperturbed because in such a case the effect of the black hole would be disrupted. It can be ensured only in the case when the nanotubular neck is tapering in the direction to its center, because this ensures the decrease of the Fermi level [15]. The related effects which appear on the nanostructures are also described in [16], where the special relativistic-like properties of the Beltrami pseudosphere naturally point to quantum field theory in curved space. In the work the finite temperature local density of states is predicted that is a realization of the Hawking-Unruh effect. Mentioned effect of the graphene black hole could eventually disappear in the presence of external magnetic (electric) field which would cause the transfer of the charge from one wormhole sheet to another one through a nanotube center. This serves as an important model for further investigations of the electron flux in the presence of the defects with the applications in cosmological models.


2. Spinor fields in biological systems

One of the present problems in genetic engineering is the prediction of biological gene variation and the representation of corresponding genetic code. This issue emerges in the plotting graphs related to the connection curvature of a docking processes. The docking process is important in the genes of the protein structure and could be adopted instead of using a very long alphabet notation as the string sequence and the comparison of the sequences of docking. From this point of view, methods of quantum field theory, general relativity, and related tools can be of high interest. The equilibrium between the supersymmetry and the mirror symmetry of the left-handed and right-handed DNA, RNA, nucleic and amino acid molecules can be explained by anti-de Sitter (AdS) correspondence in the Yang-Mills theory and the Chern-Simon currents in biology as the curvature of the spectrum in genetic code of the protein curvature.

Today, a genetical structure is studied by standard alphabet codes A , T , C , G , and U as a sequence of strings for the representation of genetic code for various organisms without any exact definition of a new time series of genetic code [17] in contrast to standard time series modeling. With this representation [18], it is very difficult to calculate the genetic variation [19] and to perform calculations within a framework of self-consistent mathematical theory [20], namely, in the context of string theory and M- and G-theories [21, 22].

There are still attempts to perform empirical data analysis of the genetic variation [23] and to detect the pattern matching over the gene sequence by using algorithm over a standard alphabet code as their time series representation. It seems that one main problem in this field is how can we predict the genetic variation and the gene structure in the viral particle and other organisms, or in the context of new representation, the question is how we can explain the intuition behind a definition of new time series data of gene, e.g., involved in the Batalin-Vilkovisky cohomology of DNA and the viral gene structure. The Chern-Simons current and the anomaly over a superspace of cell membrane can be applied to diagnose new gene diseases, the cloning technology or the gene therapy in medicine. Moreover, a presented method can be improved also in view of describing a useless trash area of DNA, which is considered as unknown part of human genome.

On the other side, another approach based on the usage of a spinor field in the Kolmogorov space of the time series data [24] over the genetic code can represent the gene structure as the ghost and anti-ghost fields of the codon and anticodon. This can be achieved in the frameworks of supersymmetry [25] and G-theory [22]. Results of the works show that all calculations over the codon can be assumed as a new superspace of the time series representation of the gene structure [24].

In [26, 27, 28], we introduced a new representation of the genetic code in the time series using a modeling by strings and D-branes. By applying a spinor field to a superspace in time series data [29], the method allows us to develop supersymmetry for living organisms. In particular, it is possible to control the anomalies in the codon and anticodon ghost fields and construct an algebraic approach for the trash DNA.

The “gravitational” analogy of the Chern-Simons currents in a gravitational physics, emanating from a system of DNA-RNA transcriptions, could have interesting counterparts also in biology. A representation of codons in human genome, derived from the Chern-Simons currents, can be useful in biology to explain the source of connections over protein-docking states. In this perspective, adopting cohomology in biology can be useful as a new modeling tool for plotting genes with spinor field in time series data. Especially, the junk area of DNA, with repeated inactive genes, can be represented by the Chern-Simons currents with extended structures of knot states in a Laurent polynomial of knots.

Further we discuss the role of spinors in the time series of genetic code. We can denote H P 1 as a quaternionic projective space and H 0 X t as a pointed space of DNA alphabet sequence with A , T , C , and G as an equivalent class of A , T , C , G H 0 X t , DNA Φ i X t a ghost field with parity two with H 0 X t = H 0 X t where = is a pointed space. We define an equivalent class of DNA–RNA translation processes by using the notation of a master equation for an interaction between the viral RNA and the host cell DNA by x t y t = DNA RNA . The whole state space model of the viral replication cycle, embedded in the host cell, is denoted by X t , DNA / Y t , RNA = Z t , GENE = H / H = H P 1 as a moduli state space model with the definition of genetic code as an equivalent class of the map α t : X t , DNA A T C G H ; the host cell gene alphabet is defined by a hidden state space X t with the gene β i

A DNA e i πβ i 2 + 0 i + 0 j + 0 k , T DNA 0 + e i πβ i 2 i + 0 j + 0 k , C DNA 0 + 0 i + e i πβ i j + 0 k , G DNA 0 + 0 i + 0 j + e i 2 πβ i k . E29

In the retroviral RNA of the observed state space Y t , RNA is a span by gene α i with the anti-ghost field Φ i , + Y t , RNA of viral particle. We define a pair of ghost and anti-ghost field genes by a middle hidden state in mRNA and ribosomal EPA state in codon and anticodon state as the ghost and anti-ghost fields in the genetic code. One can define a mutual genetic code as passive or dual hidden states s 1 , , s 8 and active eight states s 1 , , s 8 for the spinor field in the genetic code by

A tRNA NU mRNA = e j πα i 2 + 0 i + 0 j + 0 k , U tRNA NA mRNA 0 + e j πα i 2 i + 0 j + 0 k , C tRNA NG mRNA 0 + 0 i + e j πα i j + 0 k , G tRNA NC mRNA 0 + 0 i + 0 j + e j 2 πα i k . E30

The reversed transcription process of the gene expression is defined by a moduli state space model of a coupling spinor field between the gene of a viral particle and the host cell (Figure 3)

Figure 3.

(a) The site of docking between CCR5 Δ 32 co-receptor in host T-cell membrane and V3 loop in HIV viral gp120. The cell membrane has a mirror symmetry structure of D-brane for outer layer. The anti-D-brane structure of the cell membrane is an inner layer of phospholipid. (b) The knot 4 1 model of a short exact sequence of hidden eight states in genetic code.

H P 1 = X t / Y t 1 e 2 π 4 niα e 2 π 4 mjβ m , n = 1,2,3,4 = 1 q DNA q RNA = q DNA q RNA 1 . E31

One can define Sp 1 S 7 H P 1 as a Hopf fibration of eight states of the genetic code s 1 , s 2 , s 8 S 7 = T p M , denoted by s 1 , s 2 , s 8 T p M states of the genetic code of the space of a viral RNA X t and a space of host cell DNA, Y t .

If U A α H P 1 is a chart of local coordinate in a manifold of genetic code over X t / Y t , where A α is defined over the right-hand isomer genetic code A T C G (for the simplicity we use a symbol G also for U) with their dual A α , with the mirror symmetry of a genetic code NA NT NC NG . We have a cycle and a cocycle of an orbifold as a trivialization over the tangent of the living organism manifold, so-called codon and anticodon U i U j U k . Let M g be a living organism manifold with M = H P 1 for a living organism with the Riemannian metric tensor g ij = < T A A α M , T A α M > over a tangent manifold and a cotangent manifold T A α M = S 7 with g ij defined as a tensor behavior field transformation between 64 states in a codon as the distance in a space of the genetic code with i , j = 1,2,3 , 8 . The active states in genetic code are denoted as s 1 , s 2 , s 8 , and the hidden states inside eight states or passive states are denoted as s 1 , s 1 , s 8 . The dual part of a genotype s i is defined by a tangent of a manifold of the genetic code with the Jacobian J = g ij , where the metrics is g ij = < [ s i , s j > . A smallest state in gene is defined as a pair of a genetic code by a classical notation A T , C G with the coordinates A T , C G . One can define a superstate in a pair of genes as ghost and anti-ghost fields in the genetic code with the supersymmetry of D-isomer to L-isomer from the right-hand D-state in the light of the polarization in a nucleic acid as s i Φ i to a left-hand light L-state of an isomer of a light polarization, denoted as s i Φ i +

s 1 = A T T p M , p 1 = s 1 = s 11 = A T T p M s 2 = A NA T p M , s 2 = A NA T p M s 3 = C G T p M , s 3 = C G T p M s 4 = C NC T p M , s 4 = C NC T p M s 5 = T T T p M , s 5 = T T T p M s 6 = T NA T p M , s 6 = T NA T p M s 7 = G G T p M , s 7 = G G T p M s 8 = G NC T p M , s 8 = G NC T p M E32

As known from biochemistry, there exist only s 1 = A T and s 3 = C G observational states of the living organism in nature. However, in the theory of supersymmetry of S 7 Hopf fibration, there exist eight states of the ghost fields with six hidden states in the mirror symmetry. In each state, an orbifold of the living organism, other eight states denoted p 1 p 8 exist. One can define all 64 states with the help of a notation for the Riemann tensor field g ij = < s i , s j > , e.g., g 11 = < s 1 , s 1 > , and we denote a pair of states in a gene as s 1 p 1 . It is a pair of the genetic code g 11 = < A T , A T > ; for the case of codon, we replace T with U . Finally, we have for g 11 = < ( A , T ) , A T > = < ( A , A T > since the states T are hidden states, so one has a codon AAU for g 11 .

If Γ ij k A k is a connection over a tangent of a manifold of X t / Y t with M = H P 1 of the genetic code k A T C G , one can denote Γ μν k F μν as the behavioral Yang-Mills field with its dual in AdS theory of supersymmetry with F μν = F μν . It is a behavior of a protein folding inducing a curvature between a viral glycoprotein and a host cell receptor. The behavior of a Yang-Mills field is an interaction field between the behavior of virus and host cells which can survive by a change of curvature of the protein during the evolution. It is a connection in the sense of an evolutional field not in the sense of traditional gravitational field as usual

Γ ij k = 1 2 g kl j g jk + i g jk k g jk , E33

where g ij = < s i , s j > .

For an equilibrium state of evolution of an organism, we have no change of covariant derivative for the tensor field g

g g = 0 . E34

We can define a F μν A α μ ν as a connection on a principle bundle Sp 1 S 7 H P 1 of a genetic code over a supermanifold of a viral particle A . We can use an exact sequence of the sheave cohomology O X t with the chart over the supermanifold defined by homogeneous coordinates of H P 1 for viral gene X t along the host cell gene O Y t , while the virus attachment to the host cell is defined by the coordinate in sheave O X t / Y t . We apply a supersymmetry AdS theory over a Yang-Mills field of a behavioral field of the genetic code as a connection over A α { A , T , C , G } H P 1 with the Hopf fibration of viral RNA gene F μν A α ν μ = Γ αν μ with an anti self dual field over the gene of the host cell DNA , F μν . A current J of the connection between the fields is defined by the Chern-Pontryagin density for the interaction of behavioral fields. The current varies from the curvature of docking between the behavior of the curvature over amino acid k in X t and its dual curvature in Y t while docking < F μν F μν > , where < > is an average or an expectation operator.

The connection between genes is

R ναβ μ = α A β ν μ β A α ν μ + A α A β ν μ , E35

where R ναβ μ is the curvature over a tangent space of a genetic code A α A T C G H P 1 and A α ν μ Γ αν μ = F μν is the connection of the coupling between two alphabets of two organisms, i.e., from DNA and viral RNA. The above curvature is also useful in other situations like t-RNA docking with DNA in interaction between two D-branes of DNA and RNA. In a gauge field theory of DNA and RNA genetic codes of translation process, it is the group action of Lie-algebra one form. We have an adjoint representation of the genetic code as a translation process over the codon and anticodon of t RNA . One needs to define a new value to measure the curvature in amino acids of protein structure, not in the tangent space of a genetic code, and one needs to define a curvature over the codon and anticodon to represent a curvature of proteins while they are docking. The new value needs to be unique for all states in the codon and to have a meaning of curvature of genetic code with connection over manifold of genetic code. To represent a spectrum of genetic variation as curvature in protein structure while docking, one can introduce a new quantity, so-called Chern-Simons current for biology.

One can translate a genetic code in a codon in three steps. The translation operator of group is given by a behavior matrix in Lie group, a group of supermanifold of living organism with action in three times. It generates a codon representation as an adjoint representation over gene expression, and it is a precise definition of genetic code with parity two of ghost field and anti-ghost field in the Chern-Simons current for the representation of a gene A i with the current density J A i = t 1 t n dJ A i .

Let a knot serve as a representation of anticodon in t-RNA topological structure for amino acid μ with J μ as a representation R of gauge group G of gene geometric translation as group action of transcription process; then the genetic code is an average expectation value of Wilson loop operator of coupling between hidden state of x t and y t twist D-brane and anti-D-brane over superspace of cell membrane, i.e.,

W K R = Tr R P exp K A . E36

The above term gives the asymmetric property of chiral molecule of DNA and RNA, twisted from the left hand to the right hand in a supersymmetry breaking as knot polynomial related to the connection A . By such a way, we can represent the genetic code as Laurent polynomials in variable q with integer coefficients, and for any knot K we have

J K q = i = 1 n a n q n . E37

By using the new parameter of knot q

q = e 2 π k + h E38

one can induce a spinor field for representation of genetic code, where h is the dual coexter number for group action of supersymmetry of gene expression G . It might be the source of evolution from the adaptive behavior derived from the environment. In the next, we set h = 0 in our definition of the Chern-Simons current for biology for the simpler derivation of formulas.

A Chern-Simons current J μ for anomaly quantum system of codons can be also defined as the spectrum of curvature in the genetic code for gene evolution detection. Under the definition we mean a differential 3-forms in cohomology of spin fiber S 3 over the homotopy class S 3 X t / Y t in the codon of t-RNA molecule. A path integral of gene expression is defined by the Chern-Simons theory over knots of codon and anticodon: it is defined by the interaction between codon A i and anticodon between DNA and RNA in the form of integral A i + S CS = D A i exp i S CS

S CS = k 4 π W Tr A dA + 2 3 A A A E39


J q K A i R i = < W K i R i > = < Tr R i P K i A > = D A i exp i S CS Π i W K i R i D A i exp i S CS . E40

The explicit definition of curvature over the connection of genetic code has also new meaning of the genetic spectrum current J μ , μ = 1 , 2 20 of the Chern-Simons current; it is generated from the representation of Lie group over manifold of a host cell.

An example of our approach can serve as a case of phenylalanine ( Phe ), where UUU and UUC definition is

J Phe = ϵ μαβα < 1 2 A α β A γ + 1 3 A α A β A γ > = ϵ μαβα < 1 2 A 2 ν μ d A 2 ν μ + 1 3 A 2 ν μ A 2 ν μ A 2 ν μ > E41

where we explicitly define the differential form of genetic code for Phe by dA 2 = A 2 A 2 A 2 A 3 UU UC , so we have AdA = A 2 A 2 A 2 A 2 A 3 = A 2 A 2 A 2 A 2 A 2 A 3 . The minus sign represents a linear combination of basis for codon. Therefore it follows

J Phe = ϵ μαβα < 1 2 A α β A γ + 1 3 A α A β A γ > = Tr 5 6 U U U 1 2 U U C = Tr H 3 M . E42

For a translation in reversed direction of antigen shift and drift in gene evolution theory, we can use the definition of the group action of reversed direction of time by the CPT theory for anti-ghost field in field of time series of antibody gene as

x t y t = Tr H 3 M = i = 1 3 g i x t / y t = α t y t , g 3 x t / y t β t x t / α t y t ϵ t = Tr H 3 M g F μν . E43

Then a numerical representation for spinor field of curvature in the gene expression by the Chern-Simons action is defined as follows:

S CS = k 4 π Tr A dA + 2 3 A A A E44

where k = 1,2,3 n are

J Amino = J k 2 k + 2 sin π k + 2 = D A S CS . E45

The derivation of the Chern-Simons current can be done by a simple algorithm [30], i.e., the Chern-Simons current maps the string of genetic code into numerical values by explicit formulas. It can be used to plot the time series data directly into the superspace of gene expression. We transform the alphabet string values, which cannot be computed in the classical standard definition of genetic code, into the Chern-Simons current of time series data of genetic code with k = 1,2,3 , 64 over spinor field with ground field of real values. We think that the approach is more suitable for computational programs used in data analyses.


3. Circular Artin braid group representation for spinor field in genetic code

In each cell division, the telomeres are shortened [31], and total length of DNA is changing. As the result of shorter biological clock from cell division, the living things die. In order to understand cell cocycle and division mechanism of telomerase aging, one can explain the source of cancer as a source of age acceleration and its relationship to telomere shortening mechanism. It is a source of braid group operation [32] so-called self-diffeomorphism in the genetic code. The age acceleration is a relative measurement between the chronical clock and the biological clock in telomere. Up to now, scientists understand that a telomere and telomerase are the locations of ancient viruses that rely on DNA in the chromosomes of living organisms. Telomere is composed of the repeated sequence of TTAGGG dt where 1000 dt 2000 . The size of the duplicate sequence at the end of this open chromosome is amplified by six braids caused by six superspaces in time series data of organisms. The G alphabet might be suitable to be chosen as hidden time scale in the biological clock.

Here we assume that all genetic code cannot be completely separated and biological clock in telomere length is parametrized by a hidden state of the number of dt G alphabet in TTAGGG n = dt repeated pattern of the telomere. The element is Grothendieck topology over an adjoint cofunctor; it is a self-diffeomorphism ξ : B 3 c B 3 c . The loop braid generator for B 3 c is a quaternionic field in genetic code. We define their explicit forms and their permutations over the symmetric group by a chosen basis in Clifford algebra as

σ D 1 2 1 2 1 2 , σ R 1 2 1 2 1 2 , σ P 1 2 1 2 1 2 , E46

and σ D = σ D 1 , σ R = σ R 1 , and σ P = σ P 1 . We have

σ D 1 2 1 2 1 2 , σ R 1 2 1 2 1 2 , σ P 1 2 1 2 1 2 , E47

therefore one can write eight bases for spinor field in the genetic code in braid form as follows

σ G 0,0,0 σ R 1 σ D = σ R 1 σ D , σ A 0,0,1 σ R 1 σ D = σ D σ R σ R 1 σ D , σ U 0,1,0 σ R 1 σ D = σ P σ D σ R 1 σ D , σ C 1,0,0 σ R 1 σ D = σ R σ P σ R 1 σ D , σ NA 0,1,1 σ R 1 σ D = σ R σ P σ P σ D σ R 1 σ D , σ NU 1,1,0 σ R 1 σ D = σ D σ R σ R σ P σ R 1 σ D , σ NC 1,0,1 σ R 1 σ D = σ P σ D σ D σ R σ R 1 σ D , σ NG 1,1,1 σ R 1 σ D = σ R σ P σ P σ D σ D σ R σ R 1 σ D . E48

We may also use θ = 2 πs with spin quantum number s being an integer for retrotransposon and half-integer for geneon, so that

e i θ = e 2 iπs = 1 2 s ψ 2 ψ 1 . E49

In a 3 -dimensional position space, the geneon and retrotransposon statistics operators are 1 and + 1 , respectively. By the same way, in two-dimensional position space, the abelian anyonic statistics operators e i θ are 1 -dimensional representations of eight loop braid elements σ 1 , σ 1 2 , σ 1 3 , , σ 1 8 in circular Artin braid group B 2 3 acting on the space of wave functions (Figure 4).

Figure 4.

The left picture shows biological Artin braid element σ 1 B 2 in complex plane C . There are eight equivalent classes span by eight orders of σ 1 . It is a representation of eight states in three genetic codes in codon as braid element in σ 1 A μ σ 1 A μ σ 1 A μ B 2 3 . The red color line represents the curvature from the physiology of biological time series data, and the blue color represents the active and passive behavior field layers. The right picture shows a member of loop braid group. Three circles S 1 represent the sources of closed 3 -balls B 3 c , the structure group of affine transform of 3 behavior fields in genetic code σ i A μ B 2 3 action on affine fibre bundle of behavior field in the genetic code σ ω ± , σ im ± , σ m ± B 3 c . We have affine group as loop braid group in genetic code by B 2 3 B 3 c .

3.1 Classification of loop braid group in genetic code

We classify three types of loop braid group operations; it is a representation of an anyon for protein folding. For two-dimensional representation of D-brane in loop braid group for the genetic code, we define abelian anyon for biology in (2 + 1) dimensions, the extra dimensions used to represent the homotopy path of protein folding.

The loop braid group, LB n BIO , for a genetic code has three types of generators, σ = σ D σ R σ P , ρ = σ A σ U σ C σ G , and τ = σ ω σ im σ m . The generators σ i ρ i B 2 3 i = 1 n 1 and τ i B 3 c i = 1 n fulfill the relations

τ i τ j = τ j τ i , i j τ i 2 = 1 , i = 1 , , n σ i τ j = τ j σ i , i j > 1 ρ i τ j = τ i ρ j , i j > 1 τ i ρ i = ρ i τ i + 1 , i = 1 , , n 1 τ i σ i = σ i τ i + 1 , i = 1 , , n 1 τ i + 1 σ i = ρ i σ i 1 ρ i τ i , i = 1 , , n 1 . E50

For a group operation of the genotype G and a representation ρ for gene translation as an anyon, we have ρ : G U 1 = S 2 as a representation in Ω n BIO and in LB n BIO . In order to visualize 3D folding structure of the protein structure, we define three types of loop braid group operations in biology. All loop elements of the representation of amino acids arose from group operations over the superspace of time series data. These tree types are a translation, reflection, and rotation, ρ Type I , ρ Type II , and ρ Type III . For the translation as a string of amino acids, we have the anyon Type I for biology,

ρ Type I σ i ρ Type I σ i + 1 ρ Type I σ i = ρ Type I σ i + 1 ρ Type I σ i ρ Type I σ i + 1 ,
ρ Type I σ i + 1 ρ Type I σ i ρ Type I σ i + 1 = ρ Type I σ i + 1 ρ Type I σ i + 1 ρ Type I σ i , E51
ρ Type I σ i = ρ Type I σ i + 1 .

It is isomorphic to the mapping class group of the infinitely punctured disk, a discrete set of punctures limiting to the boundary of the disk.

By an analogy with the action of the symmetric group by permutations in various mathematical settings, there exists a natural action of the braid group on n -tuples of objects or on the n -folded tensor product that involves some twistors. Let us consider an arbitrary group G , and let X be the set of all n -tuples of elements of G whose product is the identity element of G . The kernel of the homomorphism LB n BIO Ω n BIO is a subgroup of LB n BIO called pure loop braid group for biology on n strands and denoted as LP n BIO . In the pure braid, the beginning and end of each strand are in the same position. Pure braid groups fit into a short exact sequence

1 LF n 1 LP n BIO LP n 1 BIO 1 . E52

This sequence splits, and therefore pure braid groups are realized as iterated semi-direct products of free groups.

The braid group B 3 is the universal central extension of the modular group PSL 2 , with these sitting as lattices inside the universal covering group. If we define O D , O R , and O P as active layers over the superspace of DNA, RNA, and protein, O D , O R , and O P , as passive layers, we can define braid group in genetic code by a curvature inside DNA, RNA, and protein folding structure. It is a source of an acceleration of biological clock in the epigenetic code. We let σ D , σ R , σ P , σ D , σ R , σ P , and σ A μ and A μ = A , U , C , G , NA , NU , NC , NG be a loop braid group elements in the genetic code. They are the circular Artin braid groups for the genetic code. One can then define

Ψ R = σ D σ R σ D , Ψ P = σ D σ R , Ψ P = σ R σ D σ R . E53

The braid group operation gives Ψ P

σ D Ψ P σ D 1 = σ R Ψ P σ R 1 = Ψ P E54

implying that Ψ P is in the center of B 3 . It is a wave function of protein transition anyon state. If G = LB n BIO acts on X t , we get

σ i D Ψ i d 1 d i 1 d i d i + 1 d n = Ψ i d 1 d i 1 d i + 1 d i + 1 1 d i d i + 1 d i + 2 d n . E55

If G = LB n BIO acts on Y t , we get

σ i R Ψ i r 1 r i 1 r i r i + 1 r n = Ψ i r 1 r i 1 r i + 1 r i + 1 1 r i r i + 1 r i + 2 r n . E56

If G = LB n BIO acts on P t = X t / Y t , we get

σ i P Ψ i p 1 p i 1 p i p i + 1 p n = Ψ i p 1 p i 1 p i + 1 p i + 1 1 p i p i + 1 p i + 2 p n . E57

Thus the elements d i and d i + 1 exchange places in DNA strand by an analogy with genetic variation. If d i is twisted by the inner automorphism corresponding to d i + 1 at the position i = dt , the product of the d components remains the identity element. It may be checked that the braid group relations are satisfied and this formula indeed defines a group action of LB n BIO on X t .


4. Conclusions

The spin-orbit coupling is an important quantity, which is measurable in the carbon nanostructures, including the graphitic wormhole (or its nanotubular part); it can also help to identify the wormhole structure in details. SOC in a graphene could be induced by the nonzero curvature; in the particular case of the wormhole with negative curvature, the chiral fermions penetrating through the connecting nanotube in the wormhole structure could be created. The two-component Dirac equation is changed into the usual four-component form. As a consequence, the chiral fermions should be detected close to the wormhole bridge; the effect is stronger if the radius of the wormhole bridge is smaller. Moreover, one can detect permanently oriented flow when the chiral fermions prefer only one direction of the massive or massless fermionic current from the upper graphene sheet to the lower one, depending on the wormhole curvature.

We also describe the role of spinor fields in the time series of genetic code. The reversed transcription process of the gene expression could be defined by a moduli state space model of a coupling spinor field between the gene of a viral particle and the host cell. As a general result, all states of codon can be computed by the Chern-Simons 3-forms. The Chern-Simons current, coming from ghost and anti-ghost fields of supersymmetry theory, can be used to define a spectrum of gene expression in new time series data where a spinor field, as alternative representation of a gene, is adopted instead of using the alphabet sequence of standard bases A , T , C , U , and G . Effort is also directed toward the explanation of the adaptive behavior of immunosystem and to find the source of cancer from the physiology of telomere malfunction in DNA repairing state. Similar examination of Holo-Hilbert spectral analysis of the Chern-Simons current in V3 loop genotypes was performed recently in [33]. A genetic variation in V3 loop genotypes was forecast by using the imaging generated from tensor correlation network with an autoregressive integrated moving average model, support spinor model, and convolutional neural network algorithms.

The reported results of the work have promissory perspective for their extension to interdisciplinary areas as machine learning, econophysics, or biological sciences.



The work is partly supported by Scientific Grant Agency VEGA Grant No. 2/0009/19 and No. 2/0153/17. R. Pinčák would like to thank the TH division in CERN for hospitality.


  1. 1. González J, Herrero J. Graphene wormholes: A condensed matter illustration of dirac fermions in curved space. Nuclear Physics B. 2010;825(3):426-443
  2. 2. Pudlak M, Pincak R. Edge states of graphene bilayer strip. European Physical Journal B. 2013;86(3):1-4. Article number: 107
  3. 3. Castro EV, López-Sancho MP, Vozmediano MAH. New type of vacancy-induced localized states in multilayer graphene. Physical Review Letters. 2010;104(3):1-4. Article number: 036802
  4. 4. Castro EV, López-Sancho MP, Vozmediano MAH. Vacancy induced zero energy modes in graphene stacks: The case of abc trilayer. Solid State Communications. 2012;152(15):1483-1488
  5. 5. Kochetov EA, Osipov VA, Pincak R. Electronic properties of disclinated flexible membrane beyond the inextensional limit: Application to graphene. Journal of Physics. Condensed Matter. 2010;22(39):1-11. Article number: 395502
  6. 6. Rostami H, Asgari R. Electronic ground-state properties of strained graphene. Physical Review B: Condensed Matter and Materials Physics. 2012;86(15):1-7. Article number: 155435
  7. 7. Jiang C-W, Zhou X, Xie R-H, Li F-L. Semiclassical molecular dynamics simulations for ultrafast processes in molecules. Quantum Matter. 2013;2(6):353-363
  8. 8. Pincak R, Smotlacha J. Analogies in electronic properties of graphene wormhole and perturbed nanocylinder. European Physical Journal B. 2013;86(11):1-7. Article number: 480
  9. 9. Alhaidari AD, Jellal A, Choubabi EB, Bahlouli H. Dynamical mass generation via space compactification. Quantum Matter. 2013;2(2):140-143
  10. 10. Tan PH, Han WP, Zhao WJ, Wu ZH, Chang K, Wang H. The shear mode of multilayer graphene. Nature Materials, 11(4):294–300, 2012
  11. 11. Alden JS, Tsen AW, Huang PY, Hovden R, Brown L, Park J, et al. Strain solitons and topological defects in bilayer graphene. Proceedings of the National Academy of Sciences. 2013;110(28):11256-11260
  12. 12. Kane CL, Mele EJ. Quantum spin hall effect in graphene. Physical Review Letters. 2005;95(22):1-4. Article number: 226801
  13. 13. Ando T. Spin-orbit interaction in carbon nanotubes. Journal of the Physical Society of Japan. 2000;69(6):1757-1763
  14. 14. Atanasov V, Saxena A. Tuning the electronic properties of corrugated graphene: Confinement, curvature, and band-gap opening. Physical Review B: Condensed Matter and Materials Physics. 2010;81(20):1-8. Article number: 205409
  15. 15. Pincak R, Pudlak M, Smotlacha J. Electronic Properties of Single and Double Wall Carbon Nanotubes. New York: Nova Science Publishers, Inc; 2012
  16. 16. Iorio A, Lambiase G. The hawking-unruh phenomenon on graphene. Physics Letters, Section B: Nuclear, Elementary Particle and High-Energy Physics. 2012;716(2):334-337
  17. 17. José MV, Zamudio GS, Morgado ER. A unified model of the standard genetic code. Royal Society Open Science. 2017;4(3):1-13. Article number: 160908
  18. 18. Hoffman M. Straightening out the protein folding puzzle. Science. 1991;253(5026):1357-1358
  19. 19. Jagadesh A, Salam AAA, Zadeh VR, Arunkumar G. Genetic analysis of neuraminidase gene of influenza a(h1n1)pdm09 virus circulating in Southwest India from 2009 to 2012. Journal of Medical Virology. 2017;89(2):202-212
  20. 20. Beyl LM. Super-lagrangians. Physical Review D. 1979;19(6):1732-1745
  21. 21. Sepehri A, Pincak R. G-theory and its reduction to f (R)-gravity in FRW universe. International Journal of Geometric Methods in Modern Physics. 2018;15(9):1-15. Article number: 1850144
  22. 22. Capozziello S, Saridakis EN, Bamba K, Sepehri A, Rahaman F, Ali AF, et al. Cosmic space and pauli exclusion principle in a system of m 0-branes. International Journal of Geometric Methods in Modern Physics. 2017;14(6):1-32. Article number: 1750095
  23. 23. Antoneli F, Forger M, Gaviria PA, Hornos JEM. On amino acid and codon assignment in algebraic models for the genetic code. International Journal of Modern Physics B. 2010;24(4):435-463
  24. 24. Kanjamapornkul K, Pinčák R. Kolmogorov space in time series data. Mathematical Methods in the Applied Sciences. 2016;39(15):4463-4483
  25. 25. Capozziello S, Pincak R, Kanjamapornkul K. Anomaly on superspace of time series data. Zeitschrift fur Naturforschung: Section A Journal of Physical Sciences. 2017;72(12):1077-1091
  26. 26. Bartoš E, Pinčák R. Identification of market trends with string and D2-brane maps. Physica. 2017;A479:57-70
  27. 27. Pinčák R. The string prediction models as invariants of time series in the forex market. Physica. 2013;A392:6414-6426
  28. 28. Pinčák R, Bartoš E. With string model to time series forecasting. Physica. 2015;A436:135-146
  29. 29. Kanjamapornkul K, Pinčák R, Bartoš E. The study of Thai stock market across the 2008 financial crisis. Physica. 2016;A462:117-133
  30. 30. Capozziello S, Pincak R, Kanjamapornkul K, Saridakis EN. The chern-simons current in systems of dna-rna transcriptions. Annalen der Physik. 2018;530(4):1-35. Article number: 1700271
  31. 31. Harley CB, Futcher AB, Greider CW. Telomeres shorten during ageing of human fibroblasts. Nature. 1990;345(6274):458-460
  32. 32. Bilson-Thompson S. Braided topology and the emergence of matter. Journal of Physics: Conference Series. 2012;360:1-4. Article number: 012056
  33. 33. Pinčák R, Kanjamapornkul K, Bartoš E. Forecasting laurent polynomial in the chern–simons current of v3 loop time series. Annals of Physics. 2019;531:1-23. Article number: 1800375

Written By

Richard Pinčák and Erik Bartoš

Submitted: February 19th, 2019 Reviewed: July 9th, 2019 Published: February 13th, 2020