UDC 538.91
Chiral Model of Graphene Yu. P. Rybakov
Department of Theoretical Physics Peoples' Friendship University of Russia 6, Mikluho-Maklay str., Moscow, 117198, Russia
The chiral model of graphene based on the SU (2) order parameter is suggested in the long-wave approximation, the ideal graphene plane being determined by the kink-like solution. Corrugation of the graphene surface is described in the form of ripple and rings. The approximate solution corresponding to an infinite carbon nanotube is found.
Key words and phrases: graphene, order parameter, C-nanotube, chiral model.
1. Introduction. Structure of the Model
Since the very discovery of mono-atomic carbon layers called graphenes this material attracted deep interest of researchers due to its extraordinary properties concerning magnetism and high electric conductivity. The concept behind this research is the following. As is well known, the carbon atom possesses of four valence electrons in the so-called hybridized sp-states, the one of them being "free" in graphene lattice and all others forming sp-bonds with the neighbours.
It appears natural to introduce scalar a0 and 3-vector a fields corresponding to the s-orbital and the p-orbital states of the "free" electron respectively. These two fields can be combined into the unitary matrix U £ SU(2) considered as the order parameter of the model in question, the long-wave approximation being adopted, i. e.
U = a0 to + i a ■ t, (1)
where r0 is the unit 2 x 2-matrix and t are the three Pauli matrices, with the SU(2)-condition
a20 + a2 = 1 (2)
being imposed. It is convenient to construct via the differentiation of the chiral field (1) the so-called left chiral current
= U+d^U, (3)
the index ^ running 0,1, 2, 3 and denoting the derivatives with respect to the time x° = ct and the space coordinates xz, i = 1,2,3. Then the simplest Lagrangian density will be given by
C = -11 Sp(l^ 1») - 1 A2a2 (4)
and corresponds to the sigma-model approach in the field theory with the mass term. Here the constant model parameters I and A were introduced. Comparing the Lagrangian density (4) with that of the Landau—Lifshits theory corresponding to the quasiclassical long-wave approximation to the Heisenberg magnetic model [1-4], one can interpret the parameter I in (4) as the exchange energy between the atoms (per spacing).
Inserting (1) into (3) and (4) and taking into account the condition (2), one easily finds the following Lagrangian density:
£ = 11 (d^ao d»ao + d^a • ^a) - 1 A2a2. (5)
Received 15th December, 2010.
For the case of small a-excitations the equations of motion generated by (5) read as
□ a - (A2//)a = 0 and correspond to the dispersion law
w = koc, k2 = k2 + A 2/I,
which in the high-frequency approximation has the linear photon-like form.
First we begin with the static 1D configuration corresponding to the ideal graphene plane, the normal being oriented along the z-axis. In this case the order parameter has the form
U = exp(î p t3), p = p(z), with the Lagrangian density being
c = -1 ip'2 - 2a2 sin2 p. (6)
The Lagrangian (6) yields the equations of motion
2 hp" - A2 sin 2p = 0. (7)
The solution to (7) satisfying the natural boundary conditions
p(-ro) = 0, p(+ro) = n has the well-known kink-like form
p0(z) = 2arctan exp(z/1) (8)
with the characteristic thickness (length parameter)
I = VI / A (9)
and the energy per unit area
E = 2 Jdz (jpp'02 + A2 sin2 Po) = 2AV7.
2. Ripple on graphene surface
Let us now consider small static perturbations to the solution (8) in the vicinity of the ideal graphene plane, i.e. for small z. Since ^0(0) = ^/2, one finds for the perturbations £ = 5a3 and a+ = a\ + ia2 the following equations:
A£ = 0, (A- l-2)a+ = 0. (10)
The Descartes coordinates x, y being the coordinates of the ideal graphene plane z = 0, one easily finds the excitations of the periodic form:
£ = £o ekz cos hx, a+ = ekz cos Kx, K2 = h2 -A2//, (11)
where fcl > 1. The exponential increasing in z of the solution (11) signifies the instability of the ideal graphene plane first mentioned by N.D. Mermin and H. Wagner in 1966 for the case of magnetics.
Rybakov Yu. P. Chiral Model of Graphene
167
There exist also the ring excitations of the axially-symmetric form:
£ = ioekz Jo(kp), a+ = A+ e^ehz Jm(Kp), (12)
where Jm is the Bessel function of the rn-th order, m = 0,1, 2,..., and p, are the polar coordinates in the graphene plane. Thus, one concludes that the graphene plane has the tendency of bending. The corrugation of the graphene plane was observed experimentally [5-7].
In view of (10) it should also be underlined that the dispersion curve reveals the anisotropic character and has the two branches. The first one concerns the transverse a3-perturbations and has the photon-like behaviour. The second one concerns the longitudinal a\- and a2-perturbations and has the "massive" behaviour mentioned above.
3. O-nanotubes
Let us now search for the static axially-symmetric configuration of the form
U = exp(npa), (13)
where
= ^(p), (7 = T\ cos + T2 sin = n<£, n =1, 2, . . . .
The configuration (13) describes the infinite C-nanotube with the hedgehog structure in the transverse section. Substituting (13) into (4), one gets
£ = -1 iU'2 + ^ sin2 ^ - 2A2 sin2 (14)
The corresponding equations of motion for the chiral angle ip(p) read
2 p(p4>')' = (n2 + p2/l2) sin 2^. (15)
After changing the variable ^ = log(p/£) one finds in the limit n ^ p/l the solution to the equation (15) of the kink-like type:
= 2 arctan[exp(n^0 — wq)] (16)
satisfying the boundary conditions
= 0, ^(—ro) = k,
where the parameter R = l exp rq0 plays the role of the tube radius. Inserting (16) into (14), it is possible to calculate the energy of the C-nanotube per unit length:
E = 2nln2 Jd-q sin2 ^ = 4nln. (17)
It is worthwhile to remark that the integer number n in (17) is the so-called topological charge of the degree type (the winding number)
Q = y d^d^ sin ^ = n. (18)
1
s2
As follows from (17) and (18), the energy of the tube per unit length is proportional to the topological charge Q = n.
4. Discussion
The proposed chiral model of the graphene permits one to describe the ripple structure of the real graphene surface illustrating the Mermin-Wagner instability of the 2 D configurations. It also contains very simple graphene and tube solutions, the latter one confirming the existence of carbon nanotubes. In future it is desirable to include in the model the interaction with the electromagnetic field for the description of conductivity and magnetic properties.
References
1. Kosevich A. M., Ivanov B. A., Kovalev A. S. Nonlinear Magnetization Waves. Dynamical and Topological Solitons. — Kiev: Naukova Dumka, 1988. — 190 p.
2. Rybakov Yu. P., Jimenez R. Ochoa. Topological Excitations in Non-Heisenberg Ferromagnetic with Biquadratic Spin Exchange Interaction // Applied Nonlinear Dynamics. — 2001. — Vol. 9, No 4/5. — Pp. 155-161.
3. Izyumov Yu. A., Katsnelson M. I., Skryabin Yu. N. Itinerant Electron Magnetism. — M.: Nauka, 1994. — 368 p.
4. Tyablikov S. V. Methods of Quantum Theory of Magnetism. — M.: Nauka, 1975. — 528 p.
5. Fasolino A., Los J. H., Katsnelson M. I. Intrinsic Ripples in Graphene // Nature Materials. — 2007. — Vol. 6. — Pp. 858-861.
6. Controlled Ripple Texturing of Suspended Graphene and Ultra Thin Graphite Membranes / W. Bao, F. Miao, Z. Chen et al. // Nature Nanotechnology. — 2009. — Vol. 4. — Pp. 562-566.
7. The Structure of Suspended Graphene Sheets / J. C. Meyer, A. K. Geim, M. I. Katsnelson et al. // Nature. — 2007. — Vol. 446. — Pp. 60-63.
УДК 538.91
Киральная модель графена Ю. П. Рыбаков
Кафедра теоретической физики Российский университет дружбы народов 117198 Москва, ул. Миклухо-Маклая, 6, Россия
В длинноволновом приближении предложена киральная модель графена, построенная на основе матричного параметра порядка из группы Яи(2). При этом идеальная графеновая плоскость определяется кинковым решением. Описывается возмущение гра-феновой поверхности в виде ряби и колец. Находится приближённое решение, соответствующее бесконечной углеродной нанотрубке.
Ключевые слова: графен, параметр порядка, С-нанотрубка, киральная модель.