# Dislocations and torsion in graphene and related systems

## Abstract

A continuum model to study the influence of dislocations on the electronic properties of condensed matter systems is described and analyzed. The model is based on a geometrical formalism that associates a density of dislocations with the torsion tensor and uses the technique of quantum field theory in curved space. When applied to two-dimensional systems with Dirac points like graphene we find that dislocations couple in the form of vector gauge fields similar to these arising from curvature or elastic strain. We also describe the ways to couple dislocations to normal metals with a Fermi surface.

###### pacs:

81.05.Uw, 75.10.Jm, 75.10.Lp, 75.30.Ds## I Introduction

Since the experimental synthesis of graphene the interest in this material has mostly focussed on its electronic properties Neto et al. (2009). As time goes by and as the experiments become more accurate the interest is being displaced to the morphological aspects and their possible influence on the electronic properties Gass et al. (2008). Part of the interest on these issues is prompted by the necessity to find the mechanisms limiting the mobility of the graphene samples Katsnelson and Geim (2008), an important issue for the potential technological applications. Another source of interest lies on the special mechanical and elastic properties of graphene recently tested in experiments Lee et al. (2009); Deshpande et al. (2009); Bao et al. (2009); Teague et al. (2009); Sun et al. (2009) and on the fascinating possibility to explore the physics of real two dimensional crystals.

Local disorder in graphene has been thoroughly studied and we refer to the review article Neto et al. (2009) for a fairly complete list of references. A different type of disorder is provided by the observation of ripples in suspended graphene Meyer et al. (2007a, b) and in graphene grown on a substrate Stolyarova et al. (2007); Ishigami et al. (2007). The elastic and mechanical properties of graphitic structures have been studied intensely in the past, mostly in the context of understanding the formation of fullerenes and nanotubes. Very little work has been done for the flat graphene sheet Fasolino et al. (2007); Castro-Neto and Kim (2008) and topological defects have been often excluded in these studies. In the fullerene literature it was established that the formation of topological defects (substitution of a hexagonal ring by other polygons) is the natural way in which the graphitic net heals vacancies and other damages produced for instance by irradiation Lee et al. (2005). Among those, disclinations (isolated pentagon or heptagon rings), dislocations (pentagon-heptagon pairs) and Stone-Wales (SW) defects (special dislocation dipoles) were found to have the least formation energy and activation barriers. Dislocations and SW defects have been observed in carbon structures Hashimoto et al. (2004) and are known to have a strong influence on the electronic properties of nanotubes. Dislocations in irradiated graphitic structures have been described in Ewels et al. (2002) and experimental observations have been reported in graphene grown on Ir in Coraux et al. (2008).

In previous works we used a general relativity formalism to study the influence of curved portions of the lattice on the electronic properties of graphene by coupling the Dirac fermions to the corresponding curved space. We found that curvature generates a fictitious gauge field - a result that is also obtained in more conventional descriptions - and, in addition, we predicted a space-dependent parameter in the kinetic term (similar to a space-dependent Fermi velocity) that depends on the intrinsic curvature of the sample and has non-trivial effects on the density of states Cortijo and Vozmediano (2007a, b); de Juan et al. (2007); Vozmediano et al. (2008) and on the conductivity Cortijo and Vozmediano (2009). In the present work we extend this formalism to include a finite density of dislocations. Dislocations are modelled by adding torsion to the graphene sheet either curved or flat.

Torsion was first considered by Cartan and then by Sciama and Kibble in order to deal with spin in General Relativity Hehl et al. (1976). Within this scheme the spin of a particle turns out to be related to the torsion just as its mass is related to the curvature. Torsion is generally ignored in the context of GR, because the torsion field is algebraic in the Einstein-Cartan theory (its equations of motion are constraints, it has no dynamics) and as a consequence there is no torsion in the absence of matter. When matter generating torsion is added, its observable effects are suppressed by the smallness of the gravitational coupling and they are negligible Shapiro (2002). These two drawbacks are overcome when torsion is included in a geometric fashion as associated to dislocations. In the context of electron systems like graphene the dynamics of the lattice defects occurs at a much higher energy than the electronic processes (in graphene, lattice processes are related to the sigma bonds and have typical energies of the order of tens of eV while the continuum relevant low energy processes are of few tens of meV). It then makes perfect sense to consider the motion of electrons in a frozen geometry and hence the fact that torsion has no dynamics is not a problem. Moreover the suppression of the coupling of the torsion and electronics degrees of freedom will be of the order of magnitude of the quotient between electronic and elastic degrees of freedom that, although small, will be much bigger than the gravitational scale.

The purpose of these notes is twofold: First we describe a continuum description that allows to explore the effects of dislocations on the electronic properties of graphene and other materials, and, second, we propose an “analog” model that amplifies the effect of torsion in a General Relativity context.

## Ii Dislocations and torsion

Topological defects in crystal lattices can be classified in two kinds: Dislocations and disclinations. A dislocation can be thought of being formed by performing a cut in a bulk material and introducing extra rows of atoms (Volterra process). These types of point defects can be easily identified by going around a closed circuit enclosing the end of the line cut with discrete lattice steps. Without defects the path will end up at the same point. In the presence of a dislocation it ends up somewhere else and an extra vector has to be added to close the path. This is called the Burgers vector. Dislocations can be further classified by the direction of the Burgers vector with respect to the plane defined by the circuit: In an edge dislocation the Burgers vector lies in this plane, while in screw dislocation it is perpendicular to it. Fig. 1 shows an edge a and screw dislocation in the square lattice characterized by a Burgers vector parallel (perpendicular) to the displacement. A disclination is similarly formed by adding a wedge of atoms, and it can be identified by computing the total angle subtended by a closed path enclosing the defect (Frank angle). A disclination dipole can also be seen as a dislocation.

¿From the point of view of geometry, dislocations represent a source of translational mismatch “anholonomy”,(see Appendix A.2 for a rigorous description of torsion), and as we will see, are very naturally described in the continuum by a connection with torsion. In differential geometry, a connection defines the notion of parallel transport. When a vector is parallel transported from to , it experiences an infinitesimal variation given by

(1) |

The equivalent of the closed circuit in the definition of the Burgers vector can be defined by taking two infinitesimal vectors and , and trying to build a parallelogram with them. For this purpose, one parallel transports the vector along :

(2) |

and the vector along

(3) |

As depicted in Fig. 2, if there is a dislocation in the region enclosed by the path, the parallelogram obtained does not close, and the part that is missing is proportional to the antisymmetric part of the connection, defined as the torsion:

(4) |

Notice that the two paths would have closed in a curved space with a standard affine connection.

In the continuum limit, one can describe a density of dislocations by the density of Burgers vector, and we see that torsion plays precisely this role. Hence the natural identification of torsion with dislocations.

The connection of torsion with the continuum theory of crystal dislocations goes back to Kondo Kondo (1952) and has been formalized in Katanaev and Volovich (1992); Kleinert (1989). Various aspects of the problem have been explored in Furtado and Moraes (1999). A nice review on the relation of gravity with topological defects in solids is Moraes (2000). The geometric approach to defects in solids Kleinert (1989); Seung and Nelson (1988) relates the metric of the curved -crystalline- surface with the deformation tensor and establishes that disclinations are associated to the curvature tensor and a finite density of dislocations generates a torsion term.

In flat surfaces screw dislocations do not exist. Dislocations in graphene are made of pentagon–heptagon pairs and they have been widely studied in connection with the properties of carbon nanotubes Saito et al. (1998) and, more recently, in the flat graphene sheets or ribbons Anokhin et al. (1994); Carpio et al. (2008); López–Sancho et al. (2009). Observations of topological defects in graphene have been reported in An et al. (2001); Hashimoto et al. (2004); Duplock et al. (2004). The stability of dislocations in graphene has been explored recently in Carpio et al. (2008) and two types of stable dislocations were found. The most common are the so-called glide dislocations made by a pentagon-heptagon pair shown in the upper part of Fig. 3. In addition we can have shuffle dislocations (lower part of Fig. 3 made by a octagon with a dangling bond. Both are edge dislocations and add one or several lines of atoms to the lattice. The difference is that shuffle dislocations contain an extra atom with respect to the glide dislocations. In an infrared continuum formulation both are seen as point defects that do not disturb the lattice out of the core. The difference with conical defects as pentagons or heptagons is that when an electron circles such a defect it will found itself displaced from the original position a distance proportional to the Burgers vector as shown in Fig. 1.

## Iii Coupling Dirac fermions to a curved space with dislocations

The model that we propose to study the low energy properties of graphite and graphene in the presence of lattice defects (disclinations and dislocations) is based on the assumption that the low–energy description of the quasiparticles in terms of massless Dirac fermions remains valid in the presence of a diluted distribution of defects. The stability of the Fermi points in multilayer arrangements of the honeycomb lattice has been demonstrated in Mañes et al. (2007). As we have described in the previous section the presence of dislocations in a crystal can be described in the continuum limit by adding torsion to the space spanned by the lattice. We will here follow a covariant approach close to the general relativity formulation what allows to consistently couple matter to geometry. This description began with the early paper of ref. González et al. (1992) where the Dirac equation on a sphere was used to solve the electronic spectrum of fullerenes. The pentagonal defects (disclinations) were described by a fictitious gauge field. The formalism was later generalized to study the influence of disclinations Cortijo and Vozmediano (2007c, d) and of smooth curvature de Juan et al. (2007); Vozmediano et al. (2008) on the density of states of graphene. These previous works showed that the covariant formulation allows to obtain sensible predictions that are in agreement with the physical observations when the infrared properties are considered. In particular it was shown that the curvature induces a fictitious vector gauge field. We will here analyze the nature of the effective gauge fields created by the presence of dislocations within the covariant framework.

Before introducing the changes induced by having torsion we summarize the situation in a curved space.

The dynamics of a massless Dirac spinor in a curved spacetime is governed by the modified Dirac equation (see Appendix A):

(5) |

where are the curved space matrices that depend on the point of the space and can be computed from the generalized anticommutation relations

and the covariant derivative operator is defined as

where is the spin connection of the spinor field whose physical implications have been discussed in de Juan et al. (2007); Vozmediano et al. (2008).

Coupling a Dirac field to a curved space with torsion gives the so-called Einstein-Cartan-Dirac theory. The details of the derivation of the Dirac equation in a space with torsion are given in Appendix A.4. Here we will simplify the description putting it in a different, more physical way. The minimal coupling of any geometrical or physical fields to the Dirac spinors adopts always the form of a covariant derivative:

(6) |

where the given vector can be an electromagnetic potential induced by a real electromagnetic field or any other real or fictitious gauge field associated to deformations or to geometrical factors Vozmediano et al. (2009). In the case of having a density of dislocations in the graphene sheet modelled by torsion we can construct two potential candidates to gauge fields. In four dimensions, from the rank three torsion tensor (26) the following vectors can be built:

(7) |

The field is a real vector and can be associated to the density of edge dislocations while is an axial (pseudo) vector associated to the density of screw dislocations. These fields couple to the vector and axial current density respectively so the full Lagrangian is:

(8) |

where and are coupling constants related to the density of edge and screw dislocations respectively. The physical consequences of the vector torsion are then similar to these coming from curvature or elasticity. The coupling of the spinor to the axial field can have important consequences as it breaks time reversal symmetry.

### iii.1 Two dimensional systems

In (2+1) dimensions the completely antisymmetric part of the torsion tensor is proportional to the Levi-Civita tensor:

(9) |

where is a pseudoscalar field.

In 2+1 dimensions there are no screw dislocations, simply because there is no third spatial dimension, and the antisymmetric part of the torsion is rather related to time displacements. Since the time components of the torsion generated by a dislocation are zero, this term vanishes always. The remaining coupling to the trace part acts a new vector gauge field with special symmetry properties associated to the dislocations structure of the sample. This field will add up to these coming from curvature or to real electromagnetic fields.

A remaining interesting question is that of how this extra gauge field couples to the electronic excitations around different Fermi points (valleys). It is known that the fictitious gauge fields arising from elastic deformations as well as these coming from curvature, give rise to effective magnetic fields pointing in opposite directions in different valleys. The reason is that the sign of the coupling of the effective vector fields depends uniquely on the definition of the spin connection (18)

(10) |

and is determined by the product of gamma matrices in it. As the effective Hamiltonians around each of the Fermi points differ in the sign of only one of the gamma matrices (they are similar to parity conjugated) the spin connection has opposite signs in the two valleys. This is irrespective of whether the connection has or not an antisymmetric part and hence applies also to the coupling induced by torsion.

### iii.2 Three dimensional systems supporting Dirac fermions

We have seen in the previous sections that the effects of dislocations in two dimensional crystals like graphene are – at most – describable in terms of a fictitious gauge field similar to the one associated to curvature. The formalism outlined in this work acquires full power when applied to three dimensional crystals supporting screw dislocations and having Dirac fermions in the spectrum. Graphite is one of the most obvious candidates where screw dislocations are very common Henning (1965) and elementary excitations behaving as Dirac fermions were predicted near the K points of the band structure in the early calculations Wallace (1947); Slonczewski and Weiss (1958). These Dirac fermions have been experimentally confirmed in Zhou et al. (2006); Orlita et al. (2008).

The other perfect candidates are the the topological insulators
Moore (2009); Xia et al. (2009) specially the type that are protected
by time reversal symmetry. The axial part of the torsion will
couple to the three dimensional Dirac fermions, break the TRS and
affect the topological stability. The density of states associated
to the screw dislocations in Zhang et al. (2009) can be reproduced with
the techniques outlined in this work and will be reported
elsewhere ^{1}

## Iv Dislocations in conventional electronic systems

Another interesting question is how dislocations affect the electronic properties of crystal systems whose quasiparticles are described by the conventional Schrodinger equation. Within the present formalism the question can be posed as: Do scalar fields see torsion?

The question is interesting and relevant since the majority of electronic systems (Fermi liquids) obey a conventional Schrodinger equation and can be modelled with scalar fields. The point to notice here is that, although scalar fields transform in quite a trivial manner under rotations and translations, building a proper Hamiltonian means building derivatives and there is where curvature and perhaps torsion plays a role. The Lagrangian for a scalar field in a curved space is

(11) |

As we see it involves only the metric and hence it does not see the torsion. In the spirit of building phenomenological couplings of the electronic density to the dislocations the most natural scalar couplings that we can form from the torsion vectors (7) are

(12) |

where is the vector associated to the trace of the torsion, i. e., to the density of edge dislocations. The coupling being quadratic in the torsion will be suppressed with respect to by an extra power of the density of dislocations. In the two dimensional case we can form the same interactions with the scalar field of (9).

Landau levels in the presence of a screw dislocation have been explored with a similar formalism in Furtado and Moraes (1999).

## V Conclusions and discussion

The proposal of modelling dislocations in a crystal with a geometrical approach including torsion within the elasticity theory is very old. What is new in our approach is the suggestion to study the influence of dislocations on the electronic properties of graphene and other electronic systems having Dirac fermions in the spectrum, by coupling the Dirac equation to a generalized curved space with torsion. In what two dimensional graphene sheets are concerned the main results of this work are somehow disappointing: spinors in two spatial dimensions do not couple to the axial part of the torsion (since there are no screw dislocations). Such a coupling would have had very important consequences. The coupling to the trace part is in the form of a vector gauge field similar to the one produced by the curvature or by strain in a more conventional approach. The spacial distribution of this new gauge field will nevertheless depend on the density of dislocations.

The most interesting applications of the formalism outlined here will be in the area of topological insulators that often support massless chiral fermions in three dimensions. In Alimohammadia and Shariati (1999) it was shown that the coupling of spinors to the axial vector induced by torsion generates at higher order a four Fermi interaction on the fermions what had an influence on neutrino oscillations. In the context of condensed matter such a coupling would renormalize the Hubbard interaction or generate it if it was absent. The presence of the torsion axial vector in the Dirac equation will have also implications on various fundamental aspects of the physics of condensed matter system described by Dirac fermions. As it breaks time reversal symmetry, it can prevent weak localization. As for the conductivity, the disorder induced by a random distribution of dislocations may change the universality class of previous models based on vacancies or other T-preserved disorder. The breakdown of time reversal symmetry will also lift the quantum protectorate associated to the Fermi points Mañes et al. (2007) and allow the opening of a gap although the torsion itself will not explain the observations reported on gaps in graphene Zhou et al. (2007); Li et al. (2009). The physical implications of the effective gauge fields induced by dislocations on graphene have been discussed in González et al. (2001); Mañes (2007); Guinea et al. (2008); Fogler et al. (2008).

## Appendix A The Dirac equation in a space with torsion and curvature

### a.1 Dirac equation in a curved space without torsion

The behavior of spinors in curved spaces is more complicated than that of scalar or vector fields because their Lorentz transformation rules do not generalize easily to arbitrary coordinate systems Weinberg (1996). Instead of the usual metric we must introduce at each point described in arbitrary coordinates, a set of locally inertial coordinates and the vielbein fields , a set of orthonormal vectors labelled by that fixes the transformation between the local and the general coordinates:

(13) |

The curved space gamma matrices satisfying the commutation relations

(14) |

are related with the constant, flat space matrices by

(15) |

The spin connection is defined from the fielbein by

(16) |

where is the covariant derivative acting on the vectors as

(17) |

The spin connection for 1/2 spinors can be written as

(18) |

where are the spin connection coefficients

(19) |

In the flat case is the usual affine connection (Christoffel symbols) related to the metric tensor by

(20) |

Finally, the determinant of the metric needed to define a scalar density lagrangian is given by

(21) |

These formulas apply to a curved space without torsion. The Dirac equation is then written as

(22) |

where is the covariant derivative acting on the spinors in the curved space and the curved gamma matrices.

### a.2 Non-singular torsion in differential geometry

In differential geometry we consider a curved manifold, described by a set of coordinates . To this manifold we associate a metric which defines the distance locally. A notion of parallel transport is also needed to compare vectors in different points, in order to construct meaningful derivatives. The connection and the metric are independent geometric objects that are usually related by the so-called metricity condition ensuring that the relative angles and magnitudes of vectors are preserved under parallel transport. The metricity condition reads:

(23) |

this is, the covariant derivative of the metric vanishes. The metricity condition fixes the symmetric part of the connection and leaves the antisymmetric part undetermined. The most general connection, compatible with the Riemannian metric, is

(24) |

where: are the Christoffel symbols and is called the contortion tensor related with the torsion tensor by

(25) |

The torsion is the antisymmetric part of the connection:

(26) |

This means that when choosing a connection for a metric space, we are always free to choose the antisymmetric part of the connection (the torsion) at will. Only the symmetric part is fixed by the metricity condition. The usual choice is to take the torsion equal to zero, so that the metric alone determines all geometric properties. This is, for example, the case of General Relativity. The only restriction on K is that

(27) |

hence in a d-dimensional space it has components. The differential geometry of a general curved space is determined by the two tensor fields and (or, equivalently ).

In a space with torsion the shortest path (geodesics) and the covariantly paralell path do not coincide and, as was mentioned in Section II, parallelograms do not close.

### a.3 Torsion in a 2D flat space

Aiming to a description of dislocations alone, we study here the special case of a space with torsion but no curvature, in 2+0 dimensions. This will be the spatial part of the space-time used to describe graphene with dislocations as the Dirac equation in a space with spatial torsion.

As we will see, life is very simple in this space. First of all, since there is no curvature, there is a global coordinate system in which the metric tensor is the identity. Or seen from the tetrad perspective, there is a coordinate system in which the tetrads are trivial, and thus the Lorentz connection and the Einstein connection are equal. (Since the tetrads are trivial, there is no difference between latin and greek indices now, and moreover, since we are in 2+0, upper and lower indices are equal.) In this particular case, with the tetrads chosen in this way, we can say that the torsion in encoded in the connection.

Moreover, the antisymmetry of the torsion tensor in 2D immediately implies that it can be written as:

(28) |

where is an arbitrary vector field, basically the density of Burgers vector. It could be thought that this vector has some further constraints due to the Bianchi identities Kleinert (1989), but it can be checked that in 2D if these identities are trivial.

### a.4 Dirac equation in a curved space with torsion

As have seen in appendix A.2, inclusion of torsion in the geometrical description simply changes the connection and hence the covariant derivative.

It would be tempting then to propose that the Dirac equation in the space with torsion and curvature is given by eq. (5) with the appropriate connection but there are some important subtleties that we will specify Shapiro (2002); Hammond (2002).

Consider the Dirac field in flat space. The massless, manifestly hermitian Dirac Langrangian may be written as:

(29) |

Noting that:

(30) |

plus the fact that a total derivative doesn’t affect the equations of motion (with suitable boundary conditions), the Lagrangian (29) can be traded for:

(31) |

Indeed, the Dirac equation can be derived by variation with respect to of any of those Lagrangians. The simpler form (31) seems to be not hermitian, but this is not a problem since we can always introduce the boundary term back.

Now consider the same problem in a manifold with curvature and torsion. Consider the hermitian Lagrangian:

(32) |

with . Now the equivalent of the chain rule is:

(33) |

Note the following important point: since the covariant derivative contains gamma matrices, we can not just commute it to make it act directly on the spinor field. The “chain rule” that we have to use is really:

(34) |

where we have used:

(35) |

The commutation with the gamma matrices has introduced a new term, the trace of the connection. This means that the equivalence of the simpler form given in eq. (31) is not just obtained by promoting the derivative to a covariant one, but it also requires the introduction of this trace Hammond (2002).

In four dimensions, the following identity:

(36) |

and the definition of the connection (18) allows to rewrite this as:

(37) |

or:

(38) |

Which reveals the well known result that, in General Relativity, fermions only couple to the antisymmetric part of the torsion Hammond (2002). In three dimensions using the identity:

(39) |

we can write the action in the form

(40) |

from where the Dirac equation can be extracted directly to read

(41) |

Note added: After completing this work we became aware of ref. Mesaros et al. (2009) where similar issues are discussed. We think that, although there is some obvious overlapping, the two works entertain different points of view and can be seen as complementary.

###### Acknowledgements.

MAHV thanks the hospitality of the Centro de Fisica do Porto where this work began and the interesting conversations held there with the Cosmology group, specially with C. Herdeiro. AC acknowledges financial support from the European Commision, Marie Curie Excellence Grant MEXT-CT-2005-023778. This research was supported by the Spanish MECD grant FIS2005-05478-C02-01 and FIS2008-00124.### Footnotes

- F. de Juan, A. Cortijo, and M. A. H. Vozmediano, work in progress.

### References

- A. H. C. Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, Reviews of Modern Physics 81, 109 (2009).
- M. H. Gass, U. Bangert, A. L. Bleloch, P. Wang, R. R. Nair, and A. K. Geim, Nature Nanotechnology 3, 676 (2008).
- M. I. Katsnelson and A. K. Geim, Philosophical Transactions R. Soc. A 366, 195 (2008).
- C. Lee, X. Wei, J. W. Kysar, and J. Hone, Science 321, 5887 (2009).
- A. Deshpande, W. Bao, F. Miao, C. N. Lau, and B. J. LeRoy, Phys. Rev. B 79, 205411 (2009).
- W. Bao, F. Miao, Z. Chen, H. Zhang, W. Jang, C. Dames, and C. N. Lau (2009), eprint arXiv:0903.0414.
- M. Teague, A. Lai, J. Velasco, C. Hughes, A. Beyer, M. Bockrath, C. Lau, and N.-C. Yeh, Nano Letters 9, 2542 (2009).
- G. F. Sun, J. F. Jia, Q. K. Xue, and L. Li, Nanotechnology 20, 355701 (2009).
- J. C. Meyer, A. K. Geim, M. I. Katsnelson, K. S. Novoselov, T. J. Booth, and S. Roth, Nature 446, 60 (2007a).
- J. C. Meyer, A. K. Geim, M. I. Katsnelson, K. S. Novoselov, D. Obergfell, S. Roth, C. Girit, and A. Zettl, Solid State Commun. 143, 101 (2007b).
- E. Stolyarova, K. T. Rim, S. Ryu, J. Maultzsch, P. Kim, L. E. Brus, T. F. Heinz, M. S. Hybertsen, and G. W. Flynn, PNAS 104, 9211 (2007).
- M. Ishigami, J. H. Chen, W. G. Cullen, M. S. Fuhrer, and E. D. Williams, Nano Letters 7, 6 (2007).
- A. Fasolino, J. H. Los, and M. I. Katsnelson, Nat. Mater. 6, 858 (2007).
- A. Castro-Neto and E. Kim, Europhysics Letters 84, 57007 (2008).
- G.-D. Lee, C. Z. Wang, E. Yoon, N.-M. Hwang, D.-Y. Kim, and K. M. Ho, Phys. Rev. Lett. 95, 205501 (2005).
- A. Hashimoto, K. Suenaga, A. Gloter, K. Urita, and S. Iijima, Nature 430, 870 (2004).
- C. P. Ewels, M. Heggie, and P. R. Briddon, Chemical Physics Letters 351, 178 (2002).
- J. Coraux, A. T. N‘Diaye, C. Busse, and T. Michely, Nano Letters 2, 565 (2008).
- A. Cortijo and M. A. H. Vozmediano, Eur. Phys. Lett. 77, 47002 (2007a).
- A. Cortijo and M. A. H. Vozmediano, Nucl. Phys. B 763, 293 (2007b).
- F. de Juan, A. Cortijo, and M. A. H. Vozmediano, Phys. Rev. B 76, 165409 (2007).
- M. A. H. Vozmediano, F. de Juan, and A. Cortijo, Journal of Physics: Conference Series 129 (2008).
- A. Cortijo and M. A. H. Vozmediano, Phys. Rev. B 79, 184205 (2009).
- F. W. Hehl, P. von der Heyde, G. D. Kerlick, and J. M. Nester, Rev. Mod. Phys. 48, 393 (1976).
- I. L. Shapiro, Phys. Rep. 357, 113 (2002).
- K. Kondo, Proceedings of the 2nd Japan national congress for applied mechanics p. 41 (1952).
- M. O. Katanaev and I. V. Volovich, Annals of Physics 216, 1 (1992).
- H. Kleinert, Gauge fields in condensed matter, vol. 2 (World Scientific, Singapore, 1989).
- C. Furtado and F. Moraes, Europhys. Lett. 45, 279 (1999).
- F. Moraes, Braz. J. Phys. 30 (2000).
- H. S. Seung and D. R. Nelson, Phys. Rev. A 38, 1005 (1988).
- R. Saito, G. Dresselhaus, and M. S. Dresselhaus, Physical Properties of Carbon Nanotubes (World Scientific, 1998).
- A. Anokhin, M. N. Gal’perin, Y. N. Gornostyrev, and M. I. Katsnelson, JETP Lett. 59, 369 (1994).
- A. Carpio, L. L. Bonilla, F. de Juan, and M. A. H. Vozmediano, New J. Phys. 10, 053021 (2008).
- M. P. López–Sancho, F. de Juan, and M. A. H. Vozmediano, Phys. Rev. B 79, 075413 (2009).
- B. An, S. Fukuyama, K. Yokogawaa, M. Yoshimura, M. Egashira, Y. Korai, and I. Mochida, Appl. Phys. Lett. 23 (2001).
- E. J. Duplock, M. Scheffler, and P. J. D. Lindan, Phys. Rev. Lett. 92 (2004).
- J. L. Mañes, F. Guinea, and M. Vozmediano, Phys. Rev. B 75, 155424 (2007).
- J. González, F. Guinea, and M. A. H. Vozmediano, Phys. Rev. Lett. 69, 172 (1992).
- A. Cortijo and M. A. H. Vozmediano, Eur. Phys. J. 148, 83 (2007c).
- A. Cortijo and M. A. H. Vozmediano, Progress in Industrial Mathematics at ECMI 2006), Springer Verlag 12 (2007d).
- M. A. H. Vozmediano, M. I. Katsnelson, and F. Guinea, Physics Repports (2009), in preparation.
- G. R. Henning, Science 147, 3659 (1965).
- P. R. Wallace, Phys. Rev. 71, 622 (1947).
- J. C. Slonczewski and P. R. Weiss, Phys. Rev. 109, 272 (1958).
- S. Y. Zhou, G.-H. Gweon, J. Graf, A. V. Fedorov, C. D. Spataru, R. D. Diehl, Y. Kopelevich, D.-H. Lee, S. G. Louie, and A. Lanzara, Nature Phys. 2, 595 (2006).
- M. Orlita, C. Faugeras, G. Martinez, D. K. Maude, M. L. Sadowski, J. M. Schneider, and M. Potemski, J. Phys.: Condens. Matter 20, 454223 (2008).
- J. Moore, Nature Physics 5, 378 (2009).
- Y. Xia, D. Qian, D. Hsieh, L. Wray, A. Pal, H. Lin, A. Bansil, D. Grauer, Y. S. Hor, R. J. Cava, et al., Nature Physics 5, 398 (2009).
- Y. Zhang, Y. Ran, and A. Vishwanath, Nature Physics 5, 289 (2009).
- M. Alimohammadia and A. Shariati, Mod.Phys.Lett. A 14, 267 (1999).
- S. Y. Zhou, G. H. Gweon, A. V. Fedorov, P. N. First, W. A. De Heer, D. H. Lee, F. Guinea, A. H. Castro Neto, and A. Lanzara, Nature Mat. 6, 770 (2007).
- G. Li, A. Luican, and E. Andrei, Phys. Rev. Lett. 102, 176804 (2009).
- J. González, F. Guinea, and M. A. H. Vozmediano, Phys. Rev. B 63, 134421 (2001).
- J. L. Mañes, Phys. Rev. B 76, 045430 (2007).
- F. Guinea, B. Horovitz, and P. L. Doussal, Phys. Rev. B 77, 205421 (2008).
- M. M. Fogler, F. Guinea, and M. I. Katsnelson, Phys. Rev. Lett. 101, 226804 (2008).
- S. Weinberg, The quantum theory of fields (Cambridge University Press, 1996).
- R. T. Hammond, Rep. Prog. Phys. 65, 599 (2002).
- A. Mesaros, D. Sadri, and J. Zaanen (2009), eprint arXiv:0909.2703.