Electronic structure and transport in graphene: quasi-relativistic Dirac – Hartry – Fock self-consistent field approximation
H. V. Grushevskaya and G. G. Krylov
Physics Department, Belarusian State University, 4 Nezalezhnasti Ave., 220030 Minsk, BELARUS
Graphene and graphene-like materials are considered today as a prominent candidates to be used in new devices with functionality based on quantum effects and(or) spin-dependent phenomena in low-dimensional systems. Technically, the main obstacle to such devices implementation is the lack of methods that provide minimizatioon of distortion of these material unique properties in bulk nanoheterostructures. Theoretical approaches and computer simulation play an important role in systematic search of this kind of nanostructures for nanoelectronic applications.
However, significant part of all theoretical consideration of graphene-like materials as well as modern model representations of charge transport in these systems on pseudo-Dirac massless fermion model, originally based on tight binding approximation and applied to the description of graphite [1, 2, 3, 4], which is a bulk material.
According to this approach , -electrons in graphene are massless fermion type quasiparticle excitations moving with the Fermi velocity. The approach has been seriously developed and successfully applied to a number of experimental situation, some related review papers on the topics and further references are in [17, 18]. There are few known key points where (at our knowledge) one could expect the necessity of some generalized consideration. The first one is the cyclotron mass dependence upon the carriers concentration. Due to weakness of the signal, modern experimental techniques can register cyclotron mass of charge carriers which is just a little smaller than 0.02 of a free electron mass [5, 6, 7, 8]. Assessments are absent whether this mechanism of conductivity prevails in the region of very small values of charge carriers concentration.
The another point is the experimentally observable carrier asymmetry in graphene. According to modern theoretical concepts the bands for pseudo-relativistic electrons and holes in graphene must be symmetrical. With this in mind in the paper  in the generalized gradient approximation there were modelled the hexagonal Si and Ge, with the same structure as in graphene. But the electrons and holes bands near the Dirac points in the Brillouin zone turned out to be strongly asymmetric ones for both cases. Firstly, a Dirac cone deformation takes place far away from a circular shape, as for the second, the Dirac velocities for the valent and conduction bands are different . Therefore, one can assume the existence of some asymmetry in the behavior of pseudo-relativistic electrons and holes of the graphene as well. Since the value of the asymmetry seems to be very small, for its experimental observation one should use a highly sensitive method, such e.g., as based on the measurement of noises . For graphene, such a method is based on large amplitudes of non-universal fluctuations of charge carriers current in the form of nonmonotonic noise in the crossover region of the scattering . At high charge densities, the contribution to the resistance of clean graphene basically gives the scattering of charge carriers on long-range impurities at ordinary (symplectic) diffusion. Regime of pseudo-diffusion with a charge carriers scattering on short-range impurities is realized in the vicinity of the Dirac points of the Brillouin zone. In the papers [12, 13] measurements of quantum interference noise in a crossover between a pseudodiffusive and symplectic regime and magnetoresistance measurements in graphene p-n junctions have been performed, which established asymmetric behaviour of pseudo-relativistic electrons and holes based on asymmetric form of non-monotonic dependence of noise and magnetoresistance. And the third known key point needed theoretical explanation is the replicas existence. A weakly interacting epitaxial graphene on the surface of Ir(111) has a non-distorted hexagonal symmetry due to the weakness of the interaction with the substrate in temperature range up to a room temperature. Therefore in ARPES (angle-resolved photo-electron spectroscopy) spectra, the perturbation of the band structure is manifested in the form of replica of the inverted Dirac cone and mini-gaps in places of quasi-crossings of replicas and the Dirac cone . Authors of paper  propose replicas existence explanation such that replicas are produced only in the areas of convergence of C and Ir atoms, that explains the weak intensity of photoelectron emission replicas and brightness of the main cone. However, besides different intensity, the asymmetry of photoelectron emission spectra also manifests itself in the fact that maxima from the replicas are below then zero maximum () of Dirac cone in the ARPES spectrum at the same angle of incidence of photons and with are equal to maxima of ARPES spectrum of replicas oppositely arranged on the hexagon. And conversely, if the corners of oppositely disposed replicas are in the neighborhood , the Dirac cone in the ARPES spectrum is located lower. The above described is possible if axes of the Dirac cone and its replicas are not parallel. It means that the top of the replica does not correspond to the corners of the hexagonal mini Brillouin zone, centered at the Dirac cone corner for the epitaxial graphene on the surface of Ir(111). It has been also demonstrated that epitaxial graphene on SiC(0001) holds a hexagonal mini Brillouin zone near the Dirac points [15, 16].
And at last, a bit more philosophical but also important comment. The majority of modern software for ab initio band structure simulations uses models being some variant of the Dirac equation or at least take into account the known relativistic corrections to the Schrödinger equation when attacking the problems. The quasi-Dirac massless fermion approach based purely on tight binding non-relativistic Hamiltonian seems to be oversimplified and hardly extendable.
With the goal to investigate the balance of exchange and correlation interactions the ab initio band structure simulations of loosely-packed solids(it was used the developed generalization of the LMTO method) have been performed and demonstrated that the strong exchange leads to appearance of an energy gap in the spectrum whereas strong correlation interactions leads to tightening of this gap . The spin-unpolarized ab initio simulations of partial electron densities of two-dimensional graphite have shown that the material is a semiconductor. Interlayer correlations tighten the energy gap that results in semi-metal behaviour of three-dimensional graphite . This means that in the absence of correlation holes, the correlation interaction in a monoatomic carbon layer (monolayer) is weak in comparison with the exchange. This theoretical prediction for spin-nonpolarized graphene were confirmed experimentally in [15, 16], where it was demonstrated a bandgap in bilayer graphene on SiC(0001) and its diminishing up to vanishing in multilayer graphene.
Experimentally manufactured quasi-two-dimensional systems, such as graphene, carbon armchair nanotubes and ribbons as well as some types of zigzag carbon nanotubes manifest metallic properties (see, for example, [17, 19]). In this regard, there are discrepancies between theoretical predictions and experimental data.
Therefore, based on the results of ab initio spin-unpolarized simulations of two-dimensional and three-dimensional graphite , we can make the following assumption. Carbon low-dimensional systems having the properties similar to graphite-like materials should possess spin-polarized electronic states with the correlation holes (a magnetic ordering). Enhancing of correlation interaction due to correlation-hole contribution leads to tightening of the energy gap and, as a consequence, the emergence of semi-metal conductivity.
The goal of this chapter is to represent a Dirac – Hartree – Fock self-consistent field quasirelativistic approximation for quasi-two-dimensional systems and to describe the origin of asymmetry of electron – correlation hole carriers in graphene-like materials.
2 Graphene model bands with correlation holes
Single atomic layer of carbon atoms (two-dimensional graphite) is called monolayer graphene. Its hexagonal structure can be represented by two triangular sublattices A, B . The primitive unit cell of the graphene contains two carbon atoms C and C belonging to the sublattices A and B respectively.
The carbon atom has four valent electrons s, p, p, p. Electrons s, p, p are hybridized in the plane of the monolayer, p-electron orbitals form a half filled band of -electron orbitals on a hexagonal lattice.
Let electrons with spin ”down” (”up”) are placed on the sublattice , and the electrons with spin ”up” (”down”) – on the sublattice , as shown in Fig. 1.
With such a symmetry of the problem, all the relevant bands of sublattices are half-filled and are formed due to correlation holes. In the representation of secondary quantization and Hartree – Fock self-consistent field approximation, when not accounting for the electron density fluctuation correlation, the hole energy is simply added to the electron energy :
because the sum of projection operators in in parentheses equals to the identity operator : .
We denote spinor wave functions of the valent electrons of graphene as and . From Fig. 1 it follows that the spinor quantum fields and are transformed into each other under the mirror reflection . Therefore, a quasi-particle excitation in the proposed model of graphene is a pair of an electron and a correlation hole. As an electron-hole pairs at the same time represent themselves their proper antiparticle, the wave functions belong to the space of Majorana bispinors , and upper and lower spin components , are transformed via different representations of the Lorenz group
one gets the following expression for the bispinor wave function of an electron in graphene:
where is a vacuum vector which consists of uncorrelated vacuum states with spin “down” and “up” : .
The density matrix is expressed through the components of bispinor (4) as
3 Equation for the density matrix
In description we will consider only valent electrons. We denote by the number of atoms in two sublattices. For valent electrons,the Dirac hamiltonian has the following form:
Here is a set of Dirac matrices, is the set of Pauli matrices, indices and enumerate s-, p-, p and p electron orbitals, indices and enumerate sublattices and atoms within them respectively, is the electron radius-vector, is the radius-vector of -th carbon atom without valent electrons (atomic core), is the charge of the atomic core, is the electron charge, is the free electron mass, is the speed of light.
where is the kinetic energy operator for a single-particle state, is the self-consistent potential, is the exchange interaction, () is the eigenvalue of a non-excited single-particle state (energy of an electron orbital for an isolated atom), , is the number of valent electrons. At the equation (8) can be considered as the equation for the Green function of the quasiparticle excitations [29, 30]:
Further, the ”electron” will be uses in the sense of a quasiparticle.
A Dirac – Hartree – Fock Hamiltonian for quasiparticle excitations in graphene can be obtained by the procedure of secondary quantization of the Dirac Hamiltonian (LABEL:Dirac-Hamiltonian):
where are relativistic analogs of operators .
The equation (10) can be rewritten for the quasiparticle field as
Here is the energy of the quasiparticle excitation.
Now, one can write the relativistic equation (12) in an explicit form:
and exchange interaction term 
4 Quasirelativistic corrections
In quasirelativistic limit it is possible to neglect lower components of bispinor on respect to upper ones, as the components of the bispinor have an order of . So, it is sufficient to find upper components to describe the behavior of the system.
Expanding the factor in curly brackets in a power series on a small parameter
we obtain the quasirelativistic Dirac – Hartree – Fock approximation for graphene:
Let us find the non-relativistic limit. With this goal in eq. (26) we write down
and leave only first order terms on :
After some elementary algebra, we transform the equation (28) to the form
Since at replacements and , first two terms do not change the form of equation, then they give the non-relativistic contributions. Quadratic summand
is a quasirelativistic correction, because its form is sensitive to the above mentioned change.
Presenting as a difference of -th energy eigenvalue for one-electron non-excited state and the energy eigenvalue for the hole : and taking into account the chain of equalities
one can rewrite the equation (32) as
As mentioned above, in non-relativistic limit the indices can be omitted and eq. (34) can be written in a final form