Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
The Impact of Herding Tendency in Modular Networks on the Diffusion of Internet Investment Products
Next Article in Special Issue
Non-Hermitian Chiral Magnetic Effect in Equilibrium
Previous Article in Journal
A Shift-Dependent Measure of Extended Cumulative Entropy and Its Applications in Blind Image Quality Assessment
Previous Article in Special Issue
Quantum Hall Conductivity in the Presence of Interactions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Elastic Deformations and Wigner–Weyl Formalism in Graphene

by
I.V. Fialkovsky
1,2,† and
M.A. Zubkov
1,3,*,†
1
Physics Department, Ariel University, Ariel 40700, Israel
2
On Leave of Absence from CMCC-Universidade Federal do ABC, 09210-170 Santo Andre, Brazil
3
On Leave of Absence from NRC “Kurchatov Institute”-ITEP, B. Cheremushkinskaya 25, Moscow 117259, Russia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Symmetry 2020, 12(2), 317; https://doi.org/10.3390/sym12020317
Submission received: 21 January 2020 / Revised: 12 February 2020 / Accepted: 17 February 2020 / Published: 23 February 2020

Abstract

:
We discuss the tight-binding models of solid state physics with the Z 2 sublattice symmetry in the presence of elastic deformations in an important particular case—the tight binding model of graphene. In order to describe the dynamics of electronic quasiparticles, the Wigner–Weyl formalism is explored. It allows the calculation of the two-point Green’s function in the presence of two slowly varying external electromagnetic fields and the inhomogeneous modification of the hopping parameters that result from elastic deformations. The developed formalism allows us to consider the influence of elastic deformations and the variations of magnetic field on the quantum Hall effect.

1. Introduction

Recently there has been a revival of interest in Wigner–Weyl formalism in both condensed matter and high energy physics. It was proposed long ago by H. Groenewold [1] and J. Moyal [2], primarily in the context of one-particle quantum mechanics. The main notions of this formalism are those of the Weyl symbol of the operator and the Wigner transformation of the functions. Correspondingly, the formalism accumulated the ideas of H. Weyl [3] and E. Wigner [4]. In quantum mechanics, the Wigner–Weyl formalism substitutes the notion of the wave function by the so called Wigner distribution, that is, the function of both coordinates and momenta. The operators of physical quantities are described by their Weyl symbols (that are also the complex-valued functions of momenta and coordinates). The product of the operators on the language of Wigner–Weyl formalism becomes the Moyal product of their Weyl symbols [5,6]. The Wigner–Weyl formalism has been applied to several problems in quantum mechanics [7,8]. Notice that certain modifications of this formalism were proposed [9,10,11,12,13,14,15,16], where the main notions were changed somehow.
Let us recall the basic notions of Wigner–Weyl formalism in quantum mechanics with the example of the one dimensional model. The Wigner distribution W ( x , p ) is a function of coordinate x and momentum p. It gives the probability that the coordinate x belongs to the interval [ a , b ] in the following way:
P [ a x b ] = 1 2 π a b W ( x , p ) d p d x .
Let A ^ be operator of a certain physical observable. Its Weyl symbol A W ( x , p ) is defined as the function in phase space, which gives the expectation value of the given quantity with respect to the Wigner distribution W ( x , p ) as follows [2,17]:
A ^ = 1 2 π A W ( x , p ) W ( x , p ) d p d x .
For the pure quantum state, the Wigner function is given by:
W ( x , p ) = d y e i p y ψ * ( x + y / 2 ) ψ ( x y / 2 ) ,
where ψ ( x ) is the wave function of the state. The formalism is readily generalized to multidimensional x and p .
The Schrodinger equation in the language of the Wigner–Weyl formalism acquires the form:
i t W ( x , p , t ) = H W ( x , p ) W ( x , p , t ) W ( x , p , t ) H W ( x , p ) .
The Moyal product of two functions f and g is given here by:
f g = f exp i 2 ( x p p x ) g .
The left arrow above the derivative shows that the derivative acts on f while the right arrow assumes the action of the derivative on g.
The definition of the Weyl symbol of an arbitrary operator A ^ is:
A W ( x , p ) = d y e i p y x + y / 2 | A ^ | x y / 2 .
We denote by H W ( x , p ) the Weyl symbol of Hamiltonian H ^ .
The Wigner–Weyl formalism was also modified somehow in order to be applied to quantum field theory. The analogue of the Wigner distribution was introduced in QCD [18,19]. It has been used in field-theoretic kinetic theory [20,21] and noncommutative field theories [22,23]. Certain applications of the Wigner–Weyl formalism were proposed to several fields of theoretical physics including cosmology [24,25,26]. Using this formalism, quantum field theory has been discussed in this aspect in relation to quantum information theory (see [27,28,29] and references therein).
Wigner–Weyl formalism has been widely applied to the study of nondissipative transport phenomena [30,31,32,33,34,35]. Using this formalism it has been shown that the response of nondissipative currents to the external field strength is expressed through the topological invariants that are robust to the smooth deformation of the system. This allows the calculation of nondissipative currents for certain complicated systems within the more simple ones connected to the original systems by a smooth deformation. Using this method, the absence of the equilibrium chiral magnetic effect [36] was demonstrated within lattice regularized field theory [34]. The anomalous quantum Hall effect was studied using Wigner–Weyl formalism for the Weyl semimetals and topological insulators [35]. In addition, Wigner–Weyl formalism allows one to derive the chiral separation effect [37] within lattice models [30,32]. The same method was also applied to the investigation of the hypothetical color-flavor locking phase in QCD [33], where the fermion zero modes on vortices were discussed. The scale magnetic effect [38] has also been investigated using the same technique [31].
Historically the momentum space topological invariants were treated mainly in the context of condensed matter physics theory [39,40,41,42,43]. They protect gapless fermions on the edges of topological insulators [44,45] and the gapless fermions in the bulk of Weyl semi-metals [43,46]. The fermion zero modes of various topological defects in 3 He are also governed by momentum space topology [47]. In high energy physics, the topological invariants in momentum space were considered in [46,48,49,50,51,52,53,54,55]. The Wigner–Weyl formalism in [30,32,34,35] was developed in the context of lattice field theory. The one-particle fermion Green’s function G ( p , q ) was considered in momentum space M ( p , q M ). It has been shown that the introduction of an Abelian external gauge field A ( x ) resulting in the Peierls substitution leads to the following equation:
Q ^ ( p A ( i p ) ) ) G ( p , q ) = δ ( p q ) .
Here Q ^ ( p ) is the lattice Dirac operator. Notice that in lattice field theory, the imaginary time is discretized on the same ground as space coordinates.
The Wigner transformation of Green’s function is defined as follows:
G W ( x , p ) d q e i x q G ( p + q / 2 , p q / 2 )
It was shown [30] that for slowly varying external fields it obeys the Groenewold equation:
G W ( x , p ) Q W ( x , p ) = 1
with the above defined Moyal product ∗ extended to the D-dimensional vectors of coordinates x and momentum p . Here Q W is the Weyl symbol of operator Q ^ ( p A ( i p ) ) .
In [56,57,58], the approach of [30,32,34,35] was further developed. In [56], the lattice model with Wilson fermions was investigated in detail. The precise expression for the Weyl symbol of the Wilson Dirac operator was derived in the presence of an arbitrarily varying external gauge field. In addition, the complete iterative solution of the Groenewold equation, Equation (2), was given. As a result the fermion propagator in the background of the arbitrary external electromagnetic field was calculated. We refer to [56] for the technical details of the Wigner–Weyl formalism in lattice models, which will be used in the present paper as well. In [57] it was shown that in the lattice models (i.e., in the tight-binding models) of solid state physics with essential inhomogeneity (caused by the varying external magnetic field), the Hall conductance integrated over the whole space is given by a topological invariant in phase space. This quantity is expressed through the Wigner transformation of the one-fermion Green’s functions. The expression for the phase space topological invariant repeats the form of the momentum space topological invariants of [30,44,45,59,60]. The difference is that Green’s functions entering this expression depend on both momenta and coordinates, and the ordinary product is substituted by the Moyal product, while the extra integration over the whole space is added (The topological invariants of [44,45,46,59,60] repeat the structure of the degree of mapping of the three-dimensional manyfold to a group of matrices). It has been shown that the value of the topological invariant in phase space responsible for the Hall conductance is robust to the introduction of disorder. Certain indications were found that it is also robust to weak Coulomb interactions.
The topological description of the Quantum Hall Effect (QHE) started from the discovery of the Thouless - Kohmoto - Nightingale - den Nijs (TKNN) invariant [61] defined in two-dimensional systems. Three dimensional topological invariants for the QHE were considered in [62]. This formalism allows one to deal with the intrinsic Anomalous Quantum Hall Effect (AQHE) and with the QHE in the presence of a constant magnetic field [63]. Unfortunately, the formalism that is based on the notion of Berry curvature does not admit the direct generalization to the QHE in the presence of a varying external magnetic field and elastic deformations when the system becomes essentially inhomogeneous. It is widely believed that the total QHE conductance is robust to the introduction of disorder and weak interactions. Expression for QHE conductivity through one-particle Green’s functions has been invented in [59,60]. In the presence of interactions the full two-point Green’s function should be substituted to the corresponding expression. It has been checked in [58], that the leading order contributions due to the Coulomb interactions do not change expression for AQHE conductivity in topological insulators. This expression has the form of an integral in momentum space over the certain composition of the interacting two-point Green’s function. It is worth mentioning, however, that there is still no proof in a general case to all orders in the perturbation theory that higher order full Green’s functions do not give contributions to the QHE. For a discussion of this issue see also [44,45]. The AQHE conductivity discussed in [58] may be applied, in particular, to Weyl semimetals [64,65,66,67,68,69].
It is difficult to overestimate the role of disorder in the Quantum Hall Effect [63,70,71,72]. One of its effects is the elimination of the Hall current from the bulk, and its concentration along the boundary. The formalism developed in [57] allows one to give an alternative to prove that the total conductance remains robust to the introduction of disorder in the majority of systems although the total current remains only along the boundary of the sample. However, for graphene there are certain complications to be discussed in the Conclusions section. Namely, when the Hall current remains only along the boundary, the QHE is absent at the half filling (neutrality point). According to the common lore, the Hall conductance is assumed to be robust to the introduction of weak interactions, at least in the presence of the sufficient amount of disorder. However, Coulomb interactions are able to give rise to the fractional QHE [63,72,73] for clean systems at very small temperatures. This, however, is out of the scope of the present paper.
Graphene [74,75,76,77,78,79,80,81,82] represents the two-dimensional Weyl semimetal. The low energy physics of its electronic quasiparticles is described by the massless Dirac equation. Therefore, it allows one to simulate in the contect of a laboratory certain features of high energy physics that cannot be observed directly. Examples of such effects are the Schwinger pair production and gravitational effects in the quantum-mechanical motion of particles. Gravity appears in graphene in the presence of elastic deformations [79]. One more exceptional feature of graphene is that (unlike the later discovered three-dimensional Weyl semimetals) it is described with a very good accuracy by the simple tight-binding model defined on the honeycomb lattice. The investigation of various features of this model (including the QHE) based on the Wigner–Weyl formalism constitutes the subject of the present paper. It is worth noting here that many phenomena in graphene can be adequately treated even with low energy continuum approximation, within appropriate (pseudo-relativistic) field-theoretical methods [83,84].
In graphene, the elastic deformations lead both to the appearance of the emergent gauge field and emergent gravitational field (see, for example, [53,79,85,86,87,88,89]). The emergent gauge field appears as the variation of the Fermi point position in momentum space while the emergent gravitational field comes as the slope’s variation of the dependence of energy on momentum (i.e., the anisotropic Fermi velocity).
Although our main aim is to investigate the tight-binding model of monolayer graphene, the paper is organized in such a way, that many of the obtained expressions may be applied to some other lattice models of solid state physics (though, only the model of graphene from this class describes quantitatively the really existing physical system). The paper is organized as follows. We start from the description of the almost arbitrary non-homogeneous lattice model in Section 2. We represent the formulation of such models in momentum space. Next, we reduce the considered class to the tight-binding models with the jumps of electrons between only the adjacent lattice sites. This section ends with the consideration of non-homogeneous tight-binding models with Z 2 sublattice symmetry. The tight-binding model of graphene belongs to this class. However it is much wider, in particular the tight-binding models defined on rectangular lattices in 2 D and 3 D remain in this class.
In Section 3, we introduce Wigner–Weyl formalism in nonhomogeneous lattice models with Z 2 sublattice symmetry. We explore the definition of the Weyl symbol of the lattice Dirac operator (entering the fermion action), which is defined through the integral in momentum space. We calculate the Weyl symbol of the Dirac operator for the considered models both in the presence of inhomogeneous hopping parameters and in the presence of the varying external electromagnetic field. Both the electromagnetic field and the hopping parameters are assumed to vary slowly, i.e., we neglect their variations on the distance of the lattice spacing. Next, we turn directly to the physics of graphene. We recall the relation between elastic deformations and the non-homogeneous hopping parameters. After that, we express the Weyl symbol of the lattice Dirac operator in graphene in the presence of elastic deformations and the electromagnetic field through the electromagnetic potential and tensor of elastic deformations.
In Section 4, we consider the relation between the Weyl symbol of the lattice Dirac operator in the considered systems and the Wigner transformation of Green’s functions. Additionally, we express the electric current through the quantities of the Wigner–Weyl formalism.
In Section 5, we again consider the lattice models of the general type. The corresponding calculations are of course applied to the case of graphene directly. Namely, we extend the results of [56] for the calculation of the Wigner transformation of the fermion Green’s function (obtained for lattice Wilson fermions on rectangular lattice) to the case of the non-homogeneous tight-binding models of the arbitrary form. It is also explained how to reconstruct the Green’s function both in momentum and coordinate representations from its known Weyl symbol.
In Section 6 we extend the consideration of [57] to the case when elastic deformations are present. Namely, we prove that for the noninteracting 2 D condensed matter model with slowly varying electromagnetic fields and elastic deformations, the Hall conductivity integrated over the whole area of the sample is given by the topological invariant in phase space composed of the Wigner transformed one-particle Green’s function. It is the same topological invariant proposed in [57]. It remains robust to the smooth modification of the model (if the modification remains local and bounded to the smooth modification of the Hamiltonian in the limited region of the sample that remains far from its boundary).
In Section 7, we apply the results of the previous sections to the discussion of Hall conductivity in graphene in the presence of both elastic deformations and inhomogeneous magnetic field. First, we recall the standard derivation of Hall conductance in the noninteracting 2 D models with constant magnetic field and constant hopping parameters. Next, this standard derivation is extended to the case of the weakly varying elastic deformations that cause varying hopping parameters that remain isotropic (i.e., their values are equal for all directions in the given point though they vary from point to point). We obtain the formula for the Hall conductance that allows one to express it through the total number of electrons in the occupied energy levels and external magnetic field. Next, we apply the topological invariant in phase space defined in Section 3 to the consideration of the QHE in graphene. The very existence of such a representation for the QHE conductance allows one to prove that it remains robust to the weak elastic deformations of arbitrary form and weak modification of magnetic field unless the topological phase transition is encountered. Both are assumed to be localized in the region that remains far from the boundaries of the sample. Finally, in this section we notice that the elastic deformations in graphene that do not cause the emergent magnetic field give rise to the isotropic hopping parameters. The corresponding displacement appears to be the analytical function of the atom coordinates of the unperturbed honeycomb lattice. For the constant external magnetic field, this allows one to derive the simple relation between the number of electrons in the occupied branches of spectrum and the value of magnetic field.
In Section 8, we end with the conclusions, discuss the obtained results, and directions for future research.
Throughout the paper the following notational conventions are used. Latin letters in subscript a , b , c numerate the spatial components of vectors. The Latin letters in superscript i , j enumerate the elementary translations. All momenta vectors are bold italic l , k , p , q from the middle of the alphabet, coordinate vectors are from its end, x , y , u , v . Operators are denoted by the Latin letters with hat Q ^ , G ^ and their matrix elements of operators by functions of two variables Q ( p , q ) . The Weyl symbols of operators are denoted by the sub-index W: ( Q ^ ) W Q W .

2. Hamiltonian for the Nonlocal Tight–Binding Model

2.1. General Case

We start our discussion with a general case of the non-local tight-binding model in the presence of the external electromagnetic field A. The discussion of the present section is applicable, in principle, not only to the tight-binding model of graphene, but also to other 2 + 1 D and 3 + 1 D tight-binding models of solid state physics.
The Hamiltonian under consideration has the form:
H x , y Ψ ¯ ( y ) f ( y , x ) e i x y d v A ( v ) Ψ ( x ) = 1 | M | 2 d p d q d l d k Ψ ¯ ( p ) f ( q , l ) Ψ ( k ) x , y e i x y d v A ( v ) e i y ( q p ) + i x ( l + k )
and here the sum is over the lattice sites x , y , while f ( x , y ) is the matrix of hopping parameters. The lattice is assumed to be infinite, which means that we neglect the finite volume effects as well as the finite temperature corrections. Thus, the integrals in the second line are over momentum space M , which is the first Brillouin zone specific for the given lattice model.
In the following, we may absorb the electromagnetic field to the definition of f. Therefore, we omit it temporarily and will restore it in appropriate expressions.
Now let us consider a less general situation of the tight-binding model with the jumps of electrons between the adjacent sites only—the nearest neighbor approximation. We discuss the case of the inhomogeneous hopping parameters, which will allow us to discuss elastic deformations. Now,
f ( y , x ) = j = 1 M δ ( y ( x + b ( j ) ) ) f ( j ) ( y )
where b ( j ) are the vectors connecting each atom to its nearest M neighbors, j = 1 , , M , and f ( j ) ( y ) is the non-uniform varying hopping parameter.
Then,
f ( q , l ) = 1 | M | j = 1 M x , y e i q x + i l y δ ( y ( x + b ( j ) ) ) f ( j ) ( y ) = 1 | M | j = 1 M y e i ( q l ) y + i q b ( j ) f ( j ) ( y ) = j = 1 M f ( j ) ( q l ) e i q b ( j )
and
H = 1 | M | j = 1 M M d p d q Ψ ¯ ( p ) f ( j ) ( p q ) e i q b ( j ) Ψ ( q ) .

2.2. The Z 2 Sublattice Symmetry

Our next simplification is the consideration of a particular case, when crystal lattice exhibits Z 2 sublattice symmetry, i.e., there are two sublattices O 1 , 2 that constitute the crystal, and there is the one to one correspondence between them generated by shift x x + b ( j ) for any j = 1 , , M and x O 1 or x x b ( j ) for x O 2 .
The points of those two sublattices are to be considered independently, which gives the sublattice index α = 1 , 2 to Ψ . We identify Ψ ( t , x ) , where x O 1 , with Ψ 1 ( t , x ) , and Ψ ( t , x ) for x O 2 is identified with Ψ 2 ( t , x ) . We set Ψ 1 ( t , x ) = 0 for x O 2 and Ψ 2 ( t , x ) = 0 for x O 1 . Then,
Ψ α ( p ) = 1 | M | 1 / 2 x O α Ψ α ( x ) e i p x , α = 1 , 2 .
The inverse transform is similar:
Ψ α ( x ) = 1 | M | 1 / 2 M d p Ψ α ( p ) e i p x , α = 1 , 2 .
Note, that Brillouin zone M is the same for both sublattices since both of them are built over the same basis vectors, and thus have the same periodicity. M is formed as space of vectors p defined by modulo transformations:
p p + g ( k )
where g ( k ) are vectors of inverse lattice that solve the system of equations:
e i g ( k ) m ( j 1 j 2 ) = 1 , j 1 , j 2 = 1 , , M
while m ( j 1 , j 2 ) = b ( j 1 ) b ( j 2 ) form each of the two sublattices O 1 , 2 . Then the hoping parameters f ( j ) become 2 × 2 matrices. Besides, we assume that the spatial hoping parameters are coordinate dependent,
f 21 ( j ) ( x + b ( j ) ) = t ( j ) ( x + b ( j ) ) ,
i.e., the values of t ( j ) may vary independently but not with time. The diagonal ones are vanishing, f 11 ( j ) = f 22 ( j ) = 0 .
The Hamiltonian then receives the form:
H = j = 1 M x O 1 y = x + b ( j ) t ( j ) ( y ) Ψ ¯ 2 ( t , y ) Ψ 1 ( t , x ) + t ( j ) ( y ) Ψ ¯ 1 ( t , x ) Ψ 2 ( t , y ) .
In what follows we will omit the temporal argument of the wave function whenever no confusion is provoked.

3. Weyl Symbol for the Lattice Dirac Operator

3.1. Lattice Dirac Operator

Let us rewrite Equation (9) in the following way:
H = x O 1 y O 2 Ψ ¯ 1 ( x ) , Ψ ¯ 2 ( y ) H ( x , y ) Ψ 1 ( x ) , Ψ 2 ( y ) T x O 1 y O 2 Ψ ¯ 2 ( y ) H 21 ( y , x ) Ψ 1 ( x ) + Ψ ¯ 1 ( x ) H 12 ( x , y ) Ψ 2 ( y )
where
H 21 ( y 2 , y 1 ) = j = 1 M δ y 2 ( y 1 + b ( j ) ) t ( j ) y 1 + y 2 2 , y 1 O 1 y 2 O 2 H 12 ( y 1 , y 2 ) = H 21 ( y 2 , y 1 ) .
Note that we define the hoping parameter by its values in the middle of the lattice links, t ( j ) y 1 + y 2 2 for better readability of consequent formulas.
We will refer to the 2 by 2 matrix operator H as to the lattice Dirac Hamiltonian, although its geometrical symmetries will only be defined after specifying b ( j ) . Along with Hamiltonian, we also introduce the Dirac operator, which enters the action and consequently will be useful for an analysis of the system’s partition function:
Q i ω H = i ω H 12 H 21 i ω .
Let us consider the off-diagonal term 21 in the Hamiltonian. It can be written in terms of the Fourier transformation as:
H 21 = 1 | M | M d p d q Ψ ¯ 2 ( p ) H 21 ( p , q ) Ψ 1 ( q )
modifying Equation (5) for two sublattices, we have:
H 12 ( q , p ) = 1 | M | j = 1 M y 1 O 1 y 2 O 2 e i q y 1 + i p y 2 δ y 2 ( y 1 + b ( j ) ) t ( j ) y 1 + y 2 2 = 1 | M | j = 1 M y 1 O 1 e i ( q p ) y 1 + i q b ( j ) t ( j ) y 1 + b ( j ) / 2 .
We can identify now the points y 1 + b ( j ) / 2 , y 1 O 1 with those situated in the middle of the lattice links along the j-th direction. We call these sets of points by O 1 / 2 ( j ) . Now we can write,
H 21 ( p , q ) = j = 1 M t ( j ) ( p q ) e i ( p + q ) b ( j ) / 2
where
t ( j ) ( p ) = 1 | M | x O 1 / 2 ( j ) t ( j ) ( x ) e i x p .
Note, that the above expression is simply a Fourier transformation of a function shifted by b ( j ) / 2 since t ( j ) are only defined in the middle of the links. In this particular case, when t ( j ) ( x ) = t ( j ) , i.e., if it does not depend on x , we obtain:
t ( j ) ( p ) = t ( j ) δ ( p mod g ( j ) ) .

3.2. The Definition of the Weyl Symbol in Momentum Space

We propose the following definition of the Weyl symbol of operator A ^ :
( A ^ ) W ( x , p ) = M d q A ( p + q / 2 , p q / 2 ) e i q x .
Here, the integral is over momentum space M , in which the two vectors are equivalent if they differ by g ( j ) . In particular this means that p ± q / 2 do not span the whole Brillouin zone M .
For off-diagonal components of H from the above it gives:
H 21 , W ( x , p ) = M d q e i q x j = 1 M t ( j ) ( q ) e i p b ( j ) = j = 1 M e i p b ( j ) M d q t ( j ) ( q ) e i q x .
If the hopping parameters are homogeneous, then:
H 21 , W ( x , p ) = e i p b ( j ) t ( j ) .
On the other hand, when the hopping parameters vary, we use Equation (16):
H 21 , W ( x , p ) = 1 | M | j = 1 M e i p b ( j ) y O 1 / 2 ( j ) t ( j ) ( y ) M d q e i q ( x y ) = j = 1 M e i p b ( j ) y O 1 / 2 ( j ) t ( j ) ( y ) F ( x y ) ,
where
F ( x ) = 1 | M | M d q e i q x .
Notice that for x , y O 1 / 2 ( j ) , we have y x O and thus, the function F ( y x ) vanishes for all x O 1 / 2 ( j ) except for x = y . However, it remains nonzero and oscillates for all other values of x , including continuous ones, and gives unity if summed over O 1 / 2 ( j ) for any x ,
y O 1 / 2 ( j ) F ( x y ) = 1 .
Each term of the j-sum in Equation (19) receives a particular form if x O 1 / 2 ( j ) (with the same value of j):
H 21 , W ( j ) ( x , p ) | x O 1 / 2 ( j ) = e i p b ( j ) t ( j ) ( x )
However, Equation (19) also defines H W for continuous values of x .
The presence of an external electromagnetic field with a vector potential A may be introduced to the model via the modification of a hopping parameter in the term H 21 :
t ( j ) ( x ) t ( j ) ( x ) e i x b ( j ) / 2 x + b ( j ) / 2 A ( y ) d y .
In H 12 there should be a complex conjugate substitution:
t ( j ) ( x ) t ( j ) ( x ) e i x + b ( j ) / 2 x b ( j ) / 2 A ( y ) d y
From Equation (19) we see that it is simply the Pieirls substitution in the language of Weyl symbols.
Combining this substitution with Equation (12) we get:
Q W = j = 1 M Q W ( j )
where
Q W ( j ) ( x , p ) | x O 1 / 2 ( j ) = i ω / M t ( j ) ( x ) e i ( p b ( j ) A ( j ) ( x ) ) t ( j ) ( x ) e i ( p b ( j ) A ( j ) ( x ) ) i ω / M ,
and M is the number of the nearest neighbours. Here:
A ( j ) ( x ) = x b ( j ) / 2 x + b ( j ) / 2 A ( y ) d y .
For both t ( j ) and A that almost do not vary at the distances of order of lattice spacing we may use Equation (23) for arbitrary values of x , and get:
Q W ( x , p ) = j = 1 M i ω / N t ( j ) ( x ) e i ( p b ( j ) A ( j ) ( x ) ) t ( j ) ( x ) e i ( p b ( j ) A ( j ) ( x ) ) i ω / N .
This approximation corresponds to the situation, when the typical wavelength of electromagnetic field is much larger than the lattice spacing.

3.3. Elastic Deformation and Modification of Hoping Parameters

Now we are in the position to consider elastic deformations and Wigner–Weyl formalism in graphene. In this section we discuss the graphene monolayer in the presence of elastic deformations. The sheet of graphene is parametrized by coordinates x k , k = 1 , 2 . The displacements of each point have three components u a ( x ) , where a = 1 , 2 , 3 . The resulting coordinates of the graphene sheet embedded into three-dimensional space y a are given by:
y k ( x ) = x k + u k ( x ) , k = 1 , 2 y 3 ( x ) = u 3 ( x ) .
In the absence of the displacements, when u a = 0 , the graphene is flat. The metric of elasticity theory is given by:
g i k = δ i k + 2 u i k , u i k = 1 2 i u k + k u i + i u a k u a , a = 1 , 2 , 3 , i , k = 1 , 2 .
When the distance between the given atom and its neighbor is changed, and the change of the distance is in δ r , then the corresponding hopping parameter is modified as t t ( 1 β δ r ) , where β is the so-called Gruneisen parameter. As a result the elastic deformations change the spatial hopping parameters, which enter Equation (11), as follows:
t ( j ) ( x ) = t 1 β u i k ( x ) b i ( j ) b k ( j ) .
Here,
{ b ( j ) } j = 1 3 = { ( 1 , 0 ) ; ( 1 / 2 , 3 / 2 ) ; ( 1 / 2 , 3 / 2 ) } .
We imply that β | u i j | 1 .
The standard expression for the emergent electromagnetic potential has the form:
A 1 = β a u 12 A 2 = β 2 a ( u 22 u 11 ) .
For the arbitrarily varying field u we obtain the following expression for Q W :
Q W ( ω , p ; τ , x ) = i ω t j = 1 3 1 β u i k ( x ) b i ( j ) b k ( j ) 0 e i ( p b ( j ) A ( j ) ( x ) ) e i ( p b ( j ) A ( j ) ( x ) ) 0 .

4. Green’s Function and the Groenewold Equation

4.1. Appearance of the Moyal Product

Our definition of the Weyl symbol of operator A ^ in Equation (17) can also be written as:
A W ( x , p ) = M d P e i x P p + P 2 A ^ p P 2 .
The integral over P is over the Brillouin zone M , i.e., in M we identify the points that differ by a vector of reciprocal lattice g ( j ) .
Now let us consider the Weyl symbol ( A B ) W ( x , p ) of the product of two operators A ^ and B ^ such that their matrix elements p + q 2 A ^ p q 2 and p + q 2 B ^ p q 2 are nonzero only when q remains in the small vicinity of zero. Then:
( A B ) W ( x , p ) = M d P M d R e i x P p + P 2 A ^ R R B ^ p P 2 = 1 2 D M d P d K e i x P p + P 2 A ^ p K 2 p K 2 B ^ p P 2 = 2 D 2 D M d q d k e i x ( q + k ) p + q 2 + k 2 A ^ p q 2 + k 2 p q 2 + k 2 B ^ p q 2 k 2 = M d q d k e i x q p + q 2 A ^ p q 2 e k 2 p q 2 p e i x k p + k 2 B ^ p k 2 = M d q e i x q p + q 2 A ^ p q 2 e i 2 ( p x + x p ) M d k e i x k p + k 2 B ^ p k 2
Here the bra- and ket- vectors in momentum space are defined modulo vectors of reciprocal lattice g ( j ) , as it is inflicted by the periodicity of the lattice. In the second line we change variables:
P = q + k , K = q k
q = P + K 2 , k = P K 2
with the Jacobian:
J = 1 1 1 1 = 2 D .
This results in the factor 2 D in the third line. Here D is the dimension of space. In the present paper it may be either 2 or 3.
Hence, the Moyal product may be defined similarly to the case of continuous space:
( A B ) W ( x , p ) = A W ( x , p ) e i 2 ( x p p x ) B W ( x , p )
Notice, that for the chosen form of the Wigner transformation on a lattice above the equality is approximate and works only if the operators A ^ , B ^ are close to diagonal.

4.2. Lattice Groenewold Equation

Let us define the Fourier components of field Ψ ( τ , x ) that depends on both space coordinates x and imaginary time τ as:
Ψ α ( τ , x ) = 1 2 π | M | 1 / 2 R M d p d ω Ψ α ( ω , p ) e i p x , α = 1 , 2 .
The partition function of the considered models has the form:
Z = D Ψ ¯ D Ψ e S [ Ψ , Ψ ¯ ] S [ Ψ , Ψ ¯ ] = R M d ω d D p 2 π | M | Ψ ¯ T ( ω , p ) Q ^ Ψ ( ω , p ) .
As usual, we relate operators Q ^ and G ^ = Q ^ 1 defined in Hilbert space H of functions (on R M ) with their matrix elements Q ( p , q ) and G ( p , q ) , where the D + 1 dimensional vectors consist of the spatial parts p , q and frequencies p D + 1 , q D + 1 :
Q ( p , q ) = p | Q ^ | q , G ( p , q ) = p | Q ^ 1 | q .
It is implied that the basis of H is normalized as p | q = δ ( p D + 1 q D + 1 ) δ ( D ) ( p q ) . The mentioned operators satisfy:
Q ^ G ^ = 1
or, equivalently,
p | Q ^ G ^ | q = δ ( D + 1 ) ( p q ) .
Equation (34) may be rewritten as follows;
S [ Ψ , Ψ ¯ ] = R M d D + 1 p 1 2 π | M | R M d D + 1 p 2 2 π | M | Ψ ¯ T ( p 1 ) Q ( p 1 , p 2 ) Ψ ( p 2 )
while the Green’s function is:
G a b ( k 2 , k 1 ) = 1 Z D Ψ ¯ D Ψ W a b ( k 2 , k 1 ) e S [ Ψ , Ψ ¯ ]
where we introduced the Grassmann-valued Wigner function:
W a b ( p , q ) = Ψ b ( p ) 2 π | M | Ψ ¯ a ( q ) 2 π | M | .
Formally we may also define operator W ^ a b W ^ a b [ Ψ , Ψ ¯ ] , whose matrix elements are equal to the Wigner function, W a b ( p , q ) = p | W ^ a b [ Ψ , Ψ ¯ ] | q . Indices a , b enumerate the components of the fermionic fields and we will omit them for brevity.
We may consider the D + 1 dimensional version of the Wigner transformation of G ^ in a similar way to that of Equation (17):
G W ( x , p ) G W ( τ , x ; p D + 1 , p ) = R d q D + 1 M d q e i ( τ q D + 1 + x q ) G ( p + q / 2 , p q / 2 ) .
Its inverse then is:
G ( p + q / 2 , p q / 2 ) = 1 2 π | M | R d τ x O 1 e i ( τ ω + x q ) G W ( τ , x ; ω , p ) .
In the same way the D + 1 dimensional Weyl symbol of Q ^ may be defined. For Q ^ = τ H we obtain:
Q W ( x , p ) = i ω H W ( x , p )
where H W is the D-dimensional Weyl symbol of the Hamiltonian that was defined above.
For the slowly varying external electromagnetic field and/or in the presence of weak elastic deformations the function Q W ( x , p ) varies slowly as a function of x on the distances of order of lattice spacing. As a result, matrix elements p + q 2 Q ^ p q 2 and p + q 2 G ^ p q 2 are both nonzero in the small vicinity of q = 0 . This imposes the bounds on the value of external magnetic field B: It should be much smaller than 1 / a 2 (where a is the typical lattice spacing). In practice this means B 1000 Tesla. Then we are able to use Equations (33) and (35) becomes a lattice version of the Groenewold equation:
G W ( x , p ) Q W ( x , p ) = 1
that is,
1 = G W ( x , p ) e i 2 x p p x Q W ( x , p ) .
The Weyl symbol Q W of operator Q ^ has been calculated above and is given by Equation (30). For external fields that vary slowly on the distances of the order of lattice spacing, we are able to represent it as a function of t ( j ) ( x ) and combination p A ( x ) :
Q W ( x , p ) = Q W ( t ( j ) ( x ) , p A ( x ) ) .

4.3. Expression for the Electric Current

Let us consider the variation of the partition functionin Equation (34), corresponding to the variation of the external field A.
We note first that the action can be written as an operator trace,
S [ Ψ , Ψ ¯ ] = Tr W ^ [ Ψ , Ψ ¯ ] Q ^ ,
where W ^ [ Ψ , Ψ ¯ ] is the Wigner operator corresponding to Equation (38). Vacuum expectation value, defined in the usual way,
O ^ = 1 Z D Ψ ¯ D Ψ O ^ e S [ Ψ , Ψ ¯ ] ,
gives then for the variation of the action:
δ S = R M d p 2 π | M | tr G W ( p , x ) p k Q W ( t ( j ) ( x ) , p A ( x ) ) δ A k ( x )
where we used Equation (37) for the expectation value of W ^ and expressed trace of (almost diagonal) operators through a trace of their Weyl symbols:
Tr A ^ B ^ = Tr ( A W B W ) = x d p ( 2 π ) D + 1 tr ( A W B W ) .
Now we obtain:
δ log Z = R D + 1 d x R M d p 2 π | M | tr G W ( p , x ) p k Q W ( t ( j ) ( x ) , p A ( x ) ) δ A k ( x ) .
We used that for the slow varying fields:
x O 1 R D d x | V |
where | V | is the volume of the lattice cell. Additionally, we used the following relation between | V | and | M | :
| V | | M | = ( 2 π ) D .
The total current, i.e., the current density integrated over the whole volume of the system, appears as the response to the variation of A ( x ) integrated over the whole space–time:
J k = T R D + 1 d x R M d p ( 2 π ) D + 1 tr G W ( p , x ) p k Q W ( p , x )
Here T is temperature that is assumed to be small. We rewrite the last equation in the following way:
J k = T ( 2 π ) D + 1 Tr G W ( p , x ) p k Q W ( p A ( x ) ) .
This expression for the total current is a topological invariant, i.e., it is not changed when the system is modified continuously and provided that there are periodic boundary conditions in space. Here Tr of a Weyl symbol of an operator stands for integration over whole phase space and summation over spinor indices, if any:
Tr A W ( x , p ) R D + 1 d x R M d p tr A W ( x , p ) .

5. Calculation of the Green’s Function in the Inhomogeneous Lattice Models

5.1. Calculation of the Wigner Transformation of the Green’s Function

In this section, we propose method of calculation of electronic Green’s function in lattice models. This method is based on solving the Groenewold equation of Equation (42):
Q W G W = 1
for the Wigner transformation G W as defined in Equation (39). In the following we use Equation (30) as the definition of the Weyl transform of Q ^ . Let us also introduce the following notation:
Δ = i 2 x p p x .
The solution may be written as follows:
G W = Q W 1 + n = 1 M = 1 j k j = n k i 0 n ( 1 ) M Π i = 1 M k i ! Q W 1 Δ k 1 Q W Q W 1 Δ k 2 Q W Q W 1 Δ k M Q W ] Q W 1 = M = 0 Q W 1 ( 1 e Δ ) Q W Q W 1 ( 1 e Δ ) Q W ( 1 e Δ ) Q W ] Q W 1 M b r a c k e t s = M = 0 Q W 1 ( 1 ) Q W Q W 1 ( 1 ) Q W Q W 1 ( 1 ) Q W ] Q W 1 M b r a c k e t s
In the first row, the sum may be extended to the values M = n = 0 , then the first term will be equal to Q W 1 . Let us introduce the product operator •, which works as follows being combined with the star product introduced above:
A B C = ( A B ) C , A B C = ( A B ) C .
In the first equation, * acts both on A B and on C while in the second equation it acts only on A and B. These rules allow us to write the above equation in a compact way:
G W ( x , p ) = M = 0 Q W 1 ( 1 ) Q W Q W 1 ( 1 ) Q W Q W 1 ( 1 ) Q W Q W 1 M p r o d u c t s = M = 0 Q W 1 ( 1 ) Q W M Q W 1
We may write symbolically:
G W ( x , p ) = 1 Q W 1 ( 1 ) Q W 1 Q W 1 = Q W 1 Q W 1 Q W 1 .
In order to show that Equation (52) is indeed the solution of the Groenewold equation, let us substitute Equation (55) to the star product G W Q W and obtain:
G W Q W = M = 0 Q W 1 ( 1 ) Q W M Q W 1 Q W = M = 0 Q W 1 ( 1 ) Q W M Q W 1 ( 1 ) Q W + M = 0 Q W 1 ( 1 ) Q W M = M = 0 Q W 1 ( 1 ) Q W M + M = 0 Q W 1 ( 1 ) Q W M = Q W 1 ( 1 ) Q W 0 = 1 .

5.2. Reconstruction of Fermion Propagator from Its Wigner Transformation

Using the definitions of the Wigner transform in Equation (39) and its inverse in Equation (40), we find the Green’s function in discrete coordinate space:
G ( x 1 , x 2 ) = 1 2 π | M | R M d p 1 d p 2 G ( p 1 , p 2 ) e i p 1 x 1 i p 2 x 2 = 1 4 π 2 | M | 2 d ω 1 d ω 2 M d p 1 d p 2 d τ x O 1 e i ( p 1 p 2 ) x + i p 1 x 1 i p 2 x 2 G W x , p 1 + p 2 2 = 1 2 π | M | d ω M d p x O 1 D ( x ( x 1 + x 2 ) / 2 | p ) G W ( x , p ) e i p ( x 1 x 2 ) .
It is assumed that p i = ( ω i , p i ) and x = ( τ , x ) . Here,
D ( y | p ) = 1 | M | M d p 1 d p 2 R D d q e i q y δ ( p ( p 1 + p 2 ) / 2 ) δ ( q ( p 1 p 2 ) ) = 1 | M | M d p 1 d p 2 e i ( p 1 p 2 ) y δ ( p ( p 1 + p 2 ) / 2 ) .
Notice, that function D ( y | p ) is not equal to the lattice delta function.
In this particular case, when both hopping parameters and the external electromagnetic field vary slowly, we may substitute the sum over x by an integral, and D ( y | p ) by δ ( y ) . This gives:
G ( x 1 , x 2 ) 1 2 π | M | R M d p G W ( ( x 1 + x 2 ) / 2 , p ) e i p ( x 1 x 2 ) .

6. Total Hall Conductance as the Topological Invariant in Phase Space

6.1. Derivation in the Framework of Wigner–Weyl Formalism

We discuss here the case when D = 2 and slightly modify the derivation presented in [57]. Let us start from Equation (49) for the electric current. We represent the electromagnetic potential as a sum of the two contributions:
A = A ( M ) + A ( E )
where A ( E ) is responsible for the electric field while A ( M ) is responsible for the magnetic field. The former is assumed to be weak, and we will keep in Equation (49) the term linear in A ( E ) .
The Groenewold equation for G W may be solved iteratively. We will keep in this solution the terms linear in A ( E ) and in its first derivative. The zeroth order term (that does not contain A ( E ) at all) is denoted G W ( 0 ) . Then,
G W G W ( 0 ) + G W ( 0 ) ( p m Q W A m ) G W ( 0 ) .
Next, we further expand the second term in derivatives of A and write symbolically:
G W G W ( 0 ) + G W ( 1 ) + G W ( 2 )
where,
G W ( 1 ) = ( G W ( 0 ) p m Q W ( 0 ) G W ( 0 ) ) A m ( E ) G W ( 2 ) = i 2 ( G W ( 0 ) p m Q W p l G W ( 0 ) ) x l A m ( E ) i 2 ( p l G W ( 0 ) p m Q W ( 0 ) G W ( 0 ) ) x l A m ( E ) = i 2 ( G W ( 0 ) p l Q W ( 0 ) G W ( 0 ) p m Q W ( 0 ) G W ( 0 ) ) F l m ( E )
where we used that p l G W ( 0 ) = G W ( 0 ) p l Q W ( 0 ) G W ( 0 ) . The Q W does not depend on the derivatives of A, therefore, it is given by:
Q W = Q W ( 0 ) + p m Q W ( 0 ) A m ( E )
Upon the substitution of Equations (60) and (62) in Equation (49) the terms proportional to A ( E ) (i.e., with no derivatives) cancel each other. The remaining term proportional to the field strength F ( E ) is:
Tr G W ( 2 ) k Q W ( 0 ) = i 2 Tr ( G W ( 0 ) p l Q W ( 0 ) G W ( 0 ) p m Q W ( 0 ) G W ( 0 ) ) F l m ( E ) k Q W ( 0 ) = i 2 Tr ( G W ( 0 ) p l Q W ( 0 ) G W ( 0 ) p m Q W ( 0 ) G W ( 0 ) k Q W ( 0 ) ) F l m ( E ) .
Here we assume that the periodic boundary conditions in space are imposed in the absence of the electric field. Then the star may be inserted to the product inside the integral: G W ( 2 ) k Q W ( 0 ) G W ( 2 ) k Q W ( 0 ) . We come to the following representation of the average Hall current (i.e., the Hall current integrated over the whole area of the sample divided by this area A ) in the presence of electric field along the x 2 axis:
J 1 = N 2 π E 2
Here,
N = T ϵ i j k A 3 ! 4 π 2 Tr G W ( x , p ) Q W ( x , p ) p i G W ( x , p ) p j Q W ( x , p ) p k A ( E ) = 0
with Tr defined in Equation (50). This expression for N is a topological invariant in phase space, i.e., it is not changed if the system is modified smoothly within a finite region distant from the boundary of the sample or from infinity if the sample is infinite. This may be checked via the direct consideration of a variation of Equation (64) with respect to the variation of Q W .
Note that Equation (64) is distinct and has a wider application than the classical TKNN invariant [61] since it is also applicable to non homogeneous systems.

6.2. From Topological Invariant in Phase Space Expressed Through G W , Q W to the Standard Expression for Hall Conductance

In the previous section we showed that Hall conductance (i.e., the conductivity integrated over the whole area of the sample) is given by σ x y = N / 2 π where N is the topological invariant in phase space. Our derivation is applicable to the general case of the inhomogeneous one-particle Hamiltonian including the case when elastic deformations are present. Our next purpose is to bring Equation (64) to the conventional expression for the Hall conductance in the case, when the non-interacting charged fermions have Hamiltonian H .
First of all, one may show that Equation (64) is equivalent to the following representation for N in terms of the Green’s ’s function written in momentum representation:
N = T ( 2 π ) 3 A 3 ! 4 π 2 ϵ i j k l = 1 4 d 3 p ( l ) tr G ( p ( 1 ) , p ( 2 ) ) [ p i ( 2 ) + p i ( 3 ) ] Q ( p ( 2 ) , p ( 3 ) ) × [ p j ( 3 ) + p j ( 4 ) ] G ( p ( 3 ) , p ( 4 ) ) [ p k ( 4 ) + p k ( 1 ) ] Q ( p ( 4 ) , p ( 1 ) ) A = 0 .
This may be proved noticing that the functional trace of a product of two operators is expressed through their Weyl symbols as follows:
Tr A ^ B ^ = Tr ( A W B W ) = d 3 x d 3 p ( 2 π ) 3 tr ( A W B W ) .
Again, we need the matrix elements p + q 2 A ^ p q 2 and p + q 2 B ^ p q 2 to be nonzero only when q remains in the small vicinity of zero. Applying this formula several times to Equation (64) we come to Equation (65).
Secondly, for non-interacting fermions described by H with energy eigenstates | n : H | n = E n | n , function Q ( p ( 1 ) , p ( 2 ) ) in momentum space has the following form:
Q ( p ( 1 ) , p ( 2 ) ) p ( 1 ) | Q ^ | p ( 2 ) = δ ( 2 ) ( p ( 1 ) p ( 2 ) ) i ω ( 1 ) p ( 1 ) | H | p ( 2 ) δ ( ω ( 1 ) ω ( 2 ) )
where p = ( p 1 , p 2 , p 3 ) = ( p , ω ) . At the same time:
G ( p ( 1 ) , p ( 2 ) ) = n 1 i ω ( 1 ) E n p ( 1 ) | n n | p ( 2 ) δ ( ω ( 1 ) ω ( 2 ) )
This way we obtain:
N = i ( 2 π ) 2 8 π 2 A n , k R d ω l = 1 4 d 2 p ( l ) ϵ i j tr [ 1 ( i ω E n ) 2 p ( 1 ) | n n | p ( 2 ) [ p i ( 2 ) + p i ( 3 ) ] p ( 2 ) | H | p ( 3 ) 1 ( i ω E k ) p ( 3 ) | k k | p ( 4 ) [ p j ( 4 ) + p j ( 1 ) ] p ( 4 ) | H | p ( 1 ) ] A = 0 .
One may represent:
[ p j ( 4 ) + p j ( 1 ) ] p ( 4 ) | H | p ( 1 ) = i p ( 4 ) | H x ^ j x ^ j H | p ( 1 ) = i p ( 4 ) | [ H , x ^ j ] | p ( 1 ) .
By operator x ^ we understand i p acting on the wavefunction written in momentum representation:
x ^ j Ψ ( p ) = p | x ^ j | Ψ = i p j p | Ψ = i p j Ψ ( p ) .
Then, for example,
x ^ j δ ( 2 ) ( q p ) = p | x ^ j | q = i p j p | q = i p j δ ( 2 ) ( p q ) = i p j q | p .
Therefore, we can write
x ^ j | p = i p j | p .
Notice, that the sign minus here is counter-intuitive because the operator x ^ is typically associated with + i p . We should remember, however, that with this latter representation the derivative acts on p in the bra-vector p | rather than on p in | p . Above we have shown that the sign is changed when the derivative is transmitted to p of | p .
Thus we have:
N = i ( 2 π ) 2 8 π 2 A n , k R d ω ϵ i j n | [ H , x ^ i ] | k k | [ H , x ^ j ] | n ( i ω E n ) 2 ( i ω E k ) = 2 i ( 2 π ) 3 8 π 2 A n , k ϵ i j θ ( E n ) θ ( E k ) ( E k E n ) 2 n | [ H , x ^ i ] | k k | [ H , x ^ j ] | n .
The last expression is just the conventional expression for the Hall conductance (multiplied by 2 π ) for the given system. Notice, that it is valid for the slowly varying electromagnetic potential only (the potential almost does not vary at the distance of order of lattice spacing). Then operator x ^ = i p has the meaning of the coordinate operator.

7. Integer Quantum Hall Effect in the Presence of Varying Magnetic Field and Elastic Deformations

7.1. Constant Magnetic Field and Constant Hopping Parameters

In this subsection we repeat the standard derivation of the Hall conductance in the noninteracting 2 D models with a constant magnetic field perpendicular to the surface, and constant hopping parameters. It is assumed here that the magnetic field B is sufficiently weak, so that | B | a 2 1 , where a is the lattice spacing. Then the Hall conductivity may be represented as N / ( 2 π ) , where N is given above in Equation (67).
In order to calculate the value of N we decompose the coordinates x 1 , x 2 in relative coordinates ξ i (with bounded values) and center coordinates X i (the unbounded part):
x ^ 1 = ξ ^ 1 + X ^ 1 , x ^ 2 = ξ ^ 2 + X ^ 2
where,
ξ ^ 1 = p ^ 2 B x 1 B , X ^ 1 = p ^ 2 B
ξ ^ 2 = p ^ 1 B , X ^ 2 = p ^ 1 B x 2 B .
Then the commutation relations follow:
[ ξ ^ 1 , ξ ^ 2 ] = i B , [ X ^ 1 , X ^ 2 ] = i B , [ ξ ^ i , X ^ j ] = 0 i , j .
Since the Hamiltonian is a function of ξ i only (in Landau gauge)
H H ( ξ 1 , ξ 2 ) ,
its commutator with X j vanishes:
[ H , X ^ 1 ] = [ H , X ^ 2 ] = 0 .
We use these relations to obtain:
N = 2 i ( 2 π ) 3 8 π 2 A n , k 1 ( E k E n ) 2 n | [ H , ξ ^ i ] | k k | [ H , ξ ^ j ] | n A = 0 ϵ i j θ ( E n ) θ ( E k ) = 2 i π A n , k ϵ i j n | ξ ^ i | k k | ξ ^ j | n A = 0 θ ( E n ) θ ( E k ) = 2 i π A n n | [ ξ ^ 1 , ξ ^ 2 ] | n A = 0 θ ( E n ) = 2 i π A B n n | n θ ( E n ) .
Momentum p 2 is a good quantum number, and it enumerates the eigenstates of the Hamiltonian:
H | n = H ( p ^ 1 , p ^ 2 B x ^ 1 ) | p 2 , m = E m ( p 2 ) | p 2 , m , m Z .
We assume that the size of the system is L × L . This gives,
N = ( 2 π ) A m d p 2 L 2 π 1 B θ ( E m ( p 2 ) ) .
Average value x = p y / B plays the role of the center of orbit, and this center should belong to the interval ( L / 2 , L / 2 ) . This gives,
N = N sign ( B ) .
Here A = L 2 is the area of the system while N is the number of occupied branches of the spectrum. In conventional systems they are counted from the neutrality point, and this way we come to the standard expression for the Hall conductance of the fermionic system in the presence of constant magnetic and electric fields. On the other hand, in the exact treatment of the honeycomb lattices (e.g., for graphene) the levels are to be counted from the edge of the band. This introduces into consideration levels with large negative Chern numbers [90]. As a result, close to half filling the conductivity is given by the ‘Dirac-Landau level index’ [91] counted from zero energy.
However, our approximation in Equation (47) is probably valid only up to | E F | t , i.e., in the region between the innermost van Hove singularities [91]. Thus we cannot correctly take into account the above mentioned deep lying levels. On the other hand, for the constant magnetic field in the absence of elastic deformations their contribution is known as it cancels precisely that of N / ( 2 π ) at the half filling. We denote this term by σ x y ( 0 ) = N ( 0 ) / ( 2 π ) , and the final expression for the Hall conductivity becomes:
σ x y = N 2 π σ x y ( 0 ) ,
which we can also rewrite as:
σ x y = N 2 π sign ( B )
where N is counted from the half filling (the LLL being occupied contributes with the factor 1 / 2 ).

7.2. Constant Magnetic Field and Weakly Varying Hopping Parameters

Let us consider the case when:
t ( j ) ( x ) = t 0 ( j ) e x f ,
(with some constant spatial vector f ) for specific Hamiltonian, which we define by its Weyl symbol:
H W ( x , p ) = j = 1 M μ / N t ( j ) ( x ) e i ( p b ( j ) A ( j ) ( x ) ) t ( j ) ( x ) e i ( p b ( j ) A ( j ) ( x ) ) μ / N
and here μ is the chemical potential. As above, we decompose the coordinates x 1 , x 2 into relative coordinates ξ i and the center coordinates X i using Equation (68). We still have the commutation relations in Equation (71). However, since t ( j ) ( x ) now depend on coordinates, the Hamiltonian does not commute anymore with X 1 , X 2 . Instead we have:
H W X W , j X W , j H W = i 2 B ϵ j i f i ( H W μ ) .
Here,
X ^ 1 = p ^ 2 B , X ^ 2 = p ^ 1 B x ^ 2 B
and X W , i is their Weyl symbols. This gives,
H X ^ j X ^ j H = i 2 B ϵ j i f i ( H μ ) .
Then for n k we obtain:
n | H X ^ j X ^ j H | k = 0
and we come again to Equation (74)
N = 2 π A B n n | n θ ( E n ) = 2 π B ρ
where now ρ is the average density of occupied states.
If we require, in addition, that t ( j ) does not depend on y and depends on x only, then momentum p y is still a good quantum number, and Equation (75) is applicable. As above, we will obtain:
N = N sign ( B ) .
Recall that A = L 2 is the area of the system while N is the number of the occupied branches of spectrum. One can see, that in the presence of constant magnetic field and the hopping parameters that depend on x 1 and do not depend on x 2 (i.e., t ( j ) ( x ) = t 0 ( j ) e x 1 f 1 ) the Hall conductivity is given by the same standard value as for the constant hopping parameters. Here we also assume that weak elastic deformations are not able to modify the contribution σ x y ( 0 ) of the deep Landau Levels, so that Equation (77) remains valid, and we come finally to the conductivity (averaged over the area of the sample) as:
σ x y = N 2 π sign ( B ) ,
where the number of occupied branches of spectrum N is counted from the half-filling.

7.3. Weak Variations of Magnetic Field and Hopping Parameters

Without elastic deformations the constant external magnetic field results in the appearance of Landau Levels. Each Landau level is highly degenerate. Weak elastic deformations result in the appearance of the emergent gauge field and anisotropic Fermi velocity. Both emergent gauge field and anisotropic Fermi velocity vary in space. The latter may be interpreted, in addition, as the emergent gravitational field (see [80,81]). Both emergent fields result in a certain modification of the band structure. It is discussed, for example, in [82], where particular cases were considered. In these particular cases the elastic deformations modify the Fermi velocity and degeneracy of Landau Levels. In a general case it is expected that the degeneracy of the Landau Levels may be lifted by both spatial variations of magnetic field and elastic deformations, so that the energy levels become thick. The overall number of the one-electron states in each Landau Level may be affected as well. When the elastic deformations and variations of magnetic field become more strong, we expect that the band structure may be modified more crucially. In particular, certain Landau Levels may merge and form the common strip of the states. The gaps between different strips may disappear and, possibly, appear again in the different ranges of energies if the variations of fields become stronger.
In Section 7.1 we considered a constant magnetic field and hopping parameters. It was noticed there that although the deep LL of the honeycomb lattices (those near the edge of the band) could not be taken correctly in our approximation, their contribution was known and could be used, given by the value of Equation (64) at half filling. Since those levels are deep, they cannot be affected by weak elastic deformations and weak variations of magnetic field, which modify only the contributions of the energy levels sufficiently close to half filling.
Let us suppose, that the magnetic field is slowly varying in the limited region of the sample, such that it approaches constant value B close to the boundary of the sample. Then we come to the following result valid in the presence of varying elastic deformations as well. The total average Hall conductivity (i.e., the Hall conductivity integrated over the area A of the sample and divided by this area) has the form of:
σ x y = 1 2 π N N ( 0 )
where N is the topological invariant in phase space:
N = T A 3 ! 4 π 2 ϵ i j k R D + 1 d x R M d p tr G W ( p , x ) Q W ( p , x ) p i G W ( p , x ) p j Q W ( p , x ) p k .
In graphene, N ( 0 ) is the value of Equation (87) at half filling and unlike N it is to be calculated with a constant magnetic field and without elastic deformations.
The value of N may be affected by variations of magnetic field and by elastic deformations. The resulting Hall conductivity becomes vanishing at half filling when those variations of magnetic field are sufficiently weak as well as the elastic deformations. At the same time, we do not have any reasons to suppose that the sufficiently strong elastic deformations/variations of magnetic field cannot lead to the deviation of Hall conductivity from zero at half filling.

7.4. Analytical Elastic Deformations In Graphene

In graphene, the relation between t, A, and u is given by [79],
t ( j ) ( x ) = t [ 1 β u k l ( x ) b k ( j ) b l ( j ) ] , A 1 = β a u 12 A 2 = β 2 a ( u 22 u 11 ) .
With elementary translations given by Equation (28), the nontrivial part of t ( j ) is:
u k l ( x ) b k ( j ) b l ( j ) = a 2 4 4 u 11 u 11 + 2 3 u 12 + 3 u 22 u 11 2 3 u 12 + 3 u 22 .
Requiring that,
t ( 1 ) ( x ) = t ( 2 ) ( x ) = t ( 3 ) ( x ) ,
is consistent with Section 7.2, thus we come to the Cauchy–Riemann conditions:
u 1 x 1 = u 2 x 2 , u 2 x 1 = u 1 x 2
that is h ( z ) u 1 ( z ) + i u 2 ( z ) is analytic as a function of z = x 1 + i x 2 . (There is another solution of (90)
1 u 1 = 2 2 u 2 , 1 u 2 = 2 u 1 ,
which however brakes the smallness condition β | u i j | 1 .)
In this case we have a vanishing emergent gauge field, while up to the terms linear in the derivatives of the hopping parameters the results of Section 7.2 give:
N = 2 π B ρ
where ρ is the average density of occupied states. On the other hand, the results of Section 7.2 ensure that any weak variations of both hopping parameters and magnetic field give:
N = N sign B
where N is the number of occupied Landau Levels (now instead of the degenerate Landau level we may have the energy band parametrized by certain parameters). Comparing this result with Equation (91) we obtain:
ρ = | B | 2 π N
for the elastic deformations given by an analytical function of coordinates (i.e., when the emergent magnetic field is absent).

8. Conclusions and Discussions

In the present paper we proceeded with the development of Wigner–Weyl formalism for tight-binding models of solid state physics (or, equivalently, for the lattice regularized quantum field theory). We extended the previous works made in this direction [30,31,32,33,34,35,56,57,58]. The developed technique was applied to the class of inhomogeneous models that include, in particular, the tight-binding model of graphene in the presence of both inhomogeneous magnetic field and nontrivial elastic deformations. It is worth mentioning that the majority of our results may be applied to other models of solid state physics. Apart from the family of two-dimensional honeycomb lattice materials (graphene, germanene, silicene, etc), all rectangular lattice crystals, both in two and three dimensions could be treated with developed methods, if described within nearest-neighbor approximation. In these cases the electrons may jump only between the nearest neighbors, and there is the Z 2 sublattice symmetry. The lattice consists of the two sublattices O 1 and O 2 . For each x O 1 site x + b ( j ) O 2 with fixed vectors b ( j ) , where j = 1 , 2 , , M . For the honeycomb lattice M = 3 , for the 2 D rectangular lattice M = 4 , for the 3 D rectangular lattice M = 8 . Among the mentioned models only the two-dimensional model on the honeycomb lattice sufficiently and accurately describe the real system (graphene). Therefore, the emphasis is on the application to the physics of graphene. The particular interest in our study was the consideration of the arbitrarily varying external magnetic field and nonhomogeneous elastic deformations.
We obtained the following main results:
  • We calculated the Weyl symbol of the lattice Dirac operator (i.e., the operator Q ^ that enters the action x , y Ψ ¯ x Q x , y Ψ y ) in the presence of both elastic deformations and slowly varying external electromagnetic field:
    Q W = i ω t j 1 β u k l ( x ) b k ( j ) b l ( j ) 0 e i ( p b ( j ) A ( j ) ( r ) ) e i ( p b ( j ) A ( j ) ( r ) ) 0
    where u i j is the tensor of elastic deformations while,
    A ( j ) ( x ) = x b ( j ) / 2 x + b ( j ) / 2 A ( y ) d y .
    It was assumed that the variation of electromagnetic field A ( x ) at the distances of order of the lattice spacing may be neglected. In practice this corresponded to magnetic fields B that obeyed B a 2 1 . In practice this bound read B 1000 Tesla. Additionally, we required that the typical wavelenth of the external electromagnetic field was much larger than the lattice spacing. This did not allow the use of Equation (92) for matter interacting with the X-rays with the wavelengths of the order of several Angstroms and smaller;
  • The Wigner transformation of the electron propagator in the presence of the slowly varying magnetic field and arbitrary elastic deformations may be calculated using the following expression:
    G W ( x , p ) = M = 0 Q W 1 ( 1 e Δ ) Q W Q W 1 ( 1 e Δ ) Q W ( 1 e Δ ) Q W ] Q W 1 M b r a c k e t s
    where
    Δ = i 2 x p p x ;
  • The electron propagator in the presence of a slowly varying electromagnetic field and elastic deformations may be expressed through the Wigner transformed Green’s function as follows:
    G ( x 1 , x 2 ) 1 2 π | M | d p G W ( ( x 1 + x 2 ) / 2 , p ) e i p ( x 1 x 2 )
    where | M | is the area of the Brillouin zone.
  • The total average Hall conductivity (i.e., the Hall conductivity integrated over the area of the sample and divided by this area) in the presence of varied weak magnetic field B 1 / a 2 and elastic deformations had the form of:
    σ x y = N 2 π N ( 0 ) 2 π
    where N is the topological invariant in phase space, which is the generalization of the classical TKNN invariant [61] since it is also applicable to non homogeneous systems:
    N = T A 3 ! 4 π 2 ϵ i j k R D + 1 d x R M d p tr G W ( p , x ) Q W ( p , x ) p i G W ( p , x ) p j Q W ( p , x ) p k ,
    here A is the area of the sample. For the conventional systems, like those with the one particle Hamiltonian p 2 / 2 m (in the absence of magnetic field), when the spectrum was bounded from below, N ( 0 ) = 0 , but in graphene it was given by the value of N at half filling (with a constant magnetic field and without elastic deformations). Thus, in the presence of a constant magnetic field at half filling the Hall conductivity vanished. The value of N ( 0 ) originated from the contribution of the deep Landau Levels was not affected by variations of magnetic field and by elastic deformations. At the same time, the value of N may, in principle, be affected by them. The resulting expression worked for the magnetic field slowly varying in the limited region of the sample, such that it approached constant value B close to the boundary of the sample. The applied electric field resulted from the electric potential varying in space. The latter may be considered also as a variation in space of chemical potential or Fermi energy. Voltage U appeared as the difference of the Fermi energy between the edges of the sample. Then the Hall current is given by:
    I x y = N N ( 0 ) 2 π U
    (Recall that we used the relativistic system of units. To obtain expression for the Hall current in an ordinary system of units we have to multiply the above expression by the unity of conductance e 2 );
  • The above mentioned representation of the average Hall conductivity through the topological invariant in phase space allowed us to prove that in graphene it was robust to both sufficiently weak variations of magnetic field and sufficiently weak elastic deformations. It is worth mentioning that both mentioned variations of magnetic field and elastic deformations were to be concentrated within the finite region far from the boundary of the sample. Under these conditions Equation (96) was not changed for the smooth variations of lattice Hamiltonian (for the proof see Appendix D in [57]);
  • The special case of elastic deformations was considered, when the emergent gauge field in graphene was absent. It was shown that the corresponding deformations were given by the arbitrary analytical functions of coordinates. Namely, the condition of the absence of emergent gauge field was equivalent to the Riemann–Cauchy conditions for the displacement function u i , i = 1 , 2 . As a result the function u ( z ) = u 1 ( z ) + i u 2 ( z ) appeared to be an analytical function of z = x 1 + i x 2 , where x i were the coordinates of the carbon atoms in the unperturbed honeycomb lattice. Under these circumstances for the constant magnetic field B the Hall current was given by:
    I x y = N U 2 π sign B
    where U is voltage while N is the number of occupied Landau Levels (counted from the half filling). Now unlike the case of the unperturbed graphene the Landau Levels may already not be degenerated.
It is worth mentioning that our results were obtained in the absence of both disorder and Coulomb interactions. In the presence of disorder the Hall current density is pushed towards the boundary. However, its total value is generally believed to remain the same. This is proved, at least, in the absence of elastic deformations when the magnetic field is homogeneous It appears that the neutrality point (when chemical potential is in the middle of the Lowest Landau Level) corresponds to vanishing Hall conductivity. Again, the Landau Levels (LL) participate in the QHE being counted from the neutrality point. This occurs now because close to the boundary, the branches of energy spectrum above and below the neutrality point behave differently. Energies of those branches above the neutrality point are increased while energies of the branches situated below it are decreased. As a result there is no crossing of the energy levels with the Fermi level on the boundary at neutrality point [72], and, consequently, there are no gapless edge states that are to be the carriers of the Hall current. This shows that weak disorder does not cause a jump in the value of total conductance. The average conductivity (the conductivity integrated over the area of the sample divided by this area) is given by σ x y = N 2 π sign B , where N is the number of occupied electronic energy levels counted from the half filling. Therefore, for graphene, N may be both negative and positive. Moreover, for the chemical potential just above zero only half of the Lowest Landau Level contributed to the Hall conductance. Therefore, in this case N = 1 2 g s g v = 4 2 = 2 . We also supposed that weak variations of magnetic field and weak elastic deformations do not change this property, and the value of the average Hall conductivity should remain the same until the topological phase transition to the state with a different value of Hall conductivity is encountered. For sufficiently strong elastic deformations and/or variations of magnetic field the very notion of Landau Levels may lose its sense, but the values of Hall conductivity may still remain nonzero.
We expected that the Hall current would be robust to the weak Coulomb interactions (at least in the presence of a sufficient amount of disorder) although the detailed investigation of this issue is still to be performed (see [58] and references therein), especially in the presence of elastic deformations and variations of magnetic field. At the same time, the clean samples of graphene (very weak disorder) exhibit the Fractional Quantum Hall Effect (FQHE) due to the Coulomb interactions. The investigation of this issue remains out of the scope of the present paper although we expect that Equation (96) may still be related somehow to the description of the FQHE.
We suppose that the results obtained here may be used further in the investigation of various properties of graphene. In particular, Equations (92) and (93) determined electron propagator in the complete tight-binding model in the presence of both elastic deformations and slowly varying external electromagnetic field. This propagator may be used in those investigations of transport properties that require use of the complete tight-binding model, i.e., when the low energy effective continuum field theory of graphene is not sufficient for the solution of a particular problem. Since the form of the obtained expressions is rather complicated, and the result of Equation (93) is represented in the form of the infinite series, the practical applications of the obtained formulas are likely to require certain numerical techniques.
We also expect that the practical calculation of Hall conductivity using Equation (96) may require the application of certain numerical procedures. The possible problem to be solved using this expression is the calculation of Hall conductivity in the presence of varying magnetic field and/or varying elastic deformations. For the constant external magnetic field and without elastic deformations the result for the Hall conductance is well known. According to our results, weak elastic deformations and weak variations of magnetic field could not affect the value of the total Hall current. However, when the variations become stronger, the system may undergo a topological phase transition to the state with different value of Hall conductance. We may determine the critical values of magnetic field variation and/or deformation tensor variation using the direct evaluation of an integral in Equation (96). Both numerical and analytical methods of this evaluation await for their development.
It is worth mentioning, that the simplified version of Equation (96) (discussed in [30]) that appeared when G W did not depend on coordinates, represents the generator of the co-homology group H ( 3 ) ( M ) , where M is momentum space. Equation (96) also awaits for the interpretation using the language of algebraic topology. At the present moment we noticed only that this topological invariant certainly played a role in the classification of the homotopic classes of maps G : M R G L ( 2 , C ) , where R is the coordinate space with certain boundary conditions while M is momentum space.
We would like to notice again, that the theory presented here is valid for the slowly varying potentials, which is consistent with the requirement B a 2 1 . It would be interesting to extend the Wigner–Weyl formalism to the precise consideration of the tight-binding model of graphene in the presence of strong magnetic fields B 1 / a 2 . Another challenge is an extension of our results to the investigation of the fractional Hall effect.
The authors are grateful for the ideas and comments, as well as the collaboration in the adjacent fields, of M. Suleymanov, Xi Wu, and Chunxu Zhang. M.A.Z. is indebted for the valuable discussions to G.E.Volovik.

Author Contributions

Authors contribute equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

The present research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Groenewold, H.J. On the Principles of elementary quantum mechanics. Physica 1946, 12, 405–460. [Google Scholar] [CrossRef]
  2. Moyal, J.E. Quantum mechanics as a statistical theory. Proc. Philos. Soc. 1949, 45, 99–124. [Google Scholar] [CrossRef]
  3. Weyl, H. Quantenmechanik und Gruppentheorie. Z. Fur Phys. 1927, 46, 1–46. [Google Scholar] [CrossRef]
  4. Wigner, E.P. On the quantum correction for thermodynamic equilibrium. Phys. Rev. 1932, 40, 749–759. [Google Scholar] [CrossRef]
  5. Ali, S.T.; Englis, M. Quantization Methods: A Guide for Physicists and Analysts. Rev. Math. Phys. 2005, 17, 391–490. [Google Scholar] [CrossRef] [Green Version]
  6. Berezin, F.A.; Shubin, M.A. Colloquia Mathematica Societatis Janos Bolyai; North-Holland: Amsterdam, The Netherlands, 1972; p. 21. [Google Scholar]
  7. Curtright, T.L.; Zachos, C.K. Quantum Mechanics in Phase Space. Asia Pac. Phys. Newsl. 2012, 1, 37. [Google Scholar] [CrossRef] [Green Version]
  8. Zachos, C.; Fairlie, D.; Curtright, T. Quantum Mechanics in Phase Space; World Scientific: Singapore, 2015. [Google Scholar]
  9. Cohen, L. Generalized Phase-Space Distribution Functions. J. Math. Phys. 1966, 7, 781. [Google Scholar] [CrossRef]
  10. Agarwal, G.S.; Wolf, E. Calculus for Functions of Noncommuting Operators and General Phase-Space Methods in Quantum Mechanics. II. Quantum Mechanics in Phase Space. Phys. Rev. D 1970, 2, 2187–2205. [Google Scholar] [CrossRef]
  11. Sudarshan, E.C. Sudarshan Equivalence of Semiclassical and Quantum Mechanical Descriptions of Statistical Light Beams. Phys. Rev. Lett. 1963, 10, 277–279. [Google Scholar] [CrossRef]
  12. Glauber, R.J. Coherent and Incoherent States of the Radiation Field. Phys. Rev. 1963, 131, 2766–2788. [Google Scholar] [CrossRef]
  13. Husimi, K. Some Formal Properties of the Density Matrix. Proc. Phys. Math. Soc. Jpn. 1940, 22, 264–314. [Google Scholar]
  14. Agarwal, G.S.; Wolf, E. Calculus for Functions of Noncommuting Operators and General Phase-Space Methods in Quantum Mechanics. I. Mapping Theorems and Ordering of Functions of Noncommuting Operators. Phys. Rev. D 1970, 2, 2161–2186. [Google Scholar] [CrossRef]
  15. Cahill, K.E.; Glauber, R.J. Glauber Ordered Expansions in Boson Amplitude Operators. Phys. Rev. 1969, 177, 1882–1902. [Google Scholar] [CrossRef] [Green Version]
  16. Buot, F.A. Nonequilibrium Quantum Transport Physics in Nanosystems; World Scientific: Singapore, 2009. [Google Scholar]
  17. Lax, M. Quantum Noise. XI. Multitime correspondence between Quantum and Classical stochastic processes. Phys. Rev. 1968, 172, 350–361. [Google Scholar] [CrossRef]
  18. Lorce, C.; Pasquini, B. Quark Wigner distributions and orbital angular momentum. Phys. Rev. D 2011, 84, 014015. [Google Scholar] [CrossRef] [Green Version]
  19. Elze, H.T.; Gyulassy, M.; Vasak, D. Transport equations for the QCD Quark Wigner Operator. Nucl. Phys. B 1986, 706, 276. [Google Scholar] [CrossRef] [Green Version]
  20. Hebenstreit, F.; Alkofer, R.; Gies, H. Schwinger pair production in space and time-dependent electric fields: Relating the Wigner formalism to quantum kinetic theory. Phys. Rev. D 2010, 82, 105026. [Google Scholar] [CrossRef] [Green Version]
  21. Calzetta, E.; Habib, S.; Hu, B.L. Quantum Kinetic Field Theory in curved space-time: Covariant Wigner function and Liouville-Vlasov equation. Phys. Rev. D 1988, 37, 2901. [Google Scholar] [CrossRef]
  22. Bastos, C.; Bertolami, O.; Dias, N.C.; Prata, J.N. Weyl-Wigner formulation of noncommutative quantum mechanics. J. Math. Phys. 2008, 49, 072101. [Google Scholar] [CrossRef] [Green Version]
  23. Dayi, O.F.; Kelleyane, L.T. Wigner functions for the Landau problem in noncommutative spaces. Mod. Phys. Lett. A 2002, 17, 1937. [Google Scholar] [CrossRef] [Green Version]
  24. Habib, S.; Laflamme, R. Wigner function and decoherence in quantum cosmology. Phys. Rev. D 1990, 42, 4056. [Google Scholar] [CrossRef] [PubMed]
  25. Chapman, S.; Heinz, U.W. HBT correlators: Current formalism versus Wigner function formulation. Phys. Lett. B 1994, 340, 250. [Google Scholar] [CrossRef] [Green Version]
  26. Berry, M.V. Semi-classical mechanics in phase space: A study of Wigner’s function. Philos. Trans. R. Soc. Lond. A 1977, 287, 0145. [Google Scholar] [CrossRef]
  27. Bastos, C.; Bernardini, A.; Santos, J. Probing phase-space noncommutativity through quantum mechanics and thermodynamics of free particles and quantum rotors. Physica A 2015, 438, 340–354. [Google Scholar]
  28. Bernardini, A.E.; Bertolami, O. Non-classicality from the phase-space flow analysis of the Weyl-Wigner quantum mechanics. Eur. Phys. Lett. 2017, 120, 20002. [Google Scholar] [CrossRef] [Green Version]
  29. Bernardini, A.E. Testing nonclassicality with exact Wigner currents for an anharmonic quantum system. Phys. Rev. A 2018, 98, 052128. [Google Scholar] [CrossRef] [Green Version]
  30. Zubkov, M.A.; Khaidukov, Z.V. Topology of the momentum space, Wigner transformations, and a chiral anomaly in lattice models. JETP Lett. 2017, 106, 166. [Google Scholar] [CrossRef]
  31. Chernodub, M.N.; Zubkov, M.A. Scale magnetic effect in Quantum Electrodynamics and the Wigner-Weyl formalism. Phys. Rev. D 2017, 96, 056006. [Google Scholar] [CrossRef] [Green Version]
  32. Khaidukov, Z.V.; Zubkov, M.A. Chiral Separation Effect in lattice regularization. Phys. Rev. D 2017, 95, 074502. [Google Scholar] [CrossRef] [Green Version]
  33. Zubkov, M.A. Momentum space topology of QCD. Ann. Phys. 2018, 393, 264. [Google Scholar] [CrossRef] [Green Version]
  34. Zubkov, M.A. Absence of equilibrium chiral magnetic effect. Phys. Rev. D 2016, 93, 105036. [Google Scholar] [CrossRef] [Green Version]
  35. Zubkov, M.A. Wigner transformation, momentum space topology, and anomalous transport. Ann. Phys. 2016, 373, 298. [Google Scholar] [CrossRef] [Green Version]
  36. Kharzeev, D.E. The Chiral Magnetic Effect and Anomaly-Induced Transport. Prog. Part. Nucl. Phys. 2014, 75, 133. [Google Scholar] [CrossRef] [Green Version]
  37. Metlitski, M.A.; Zhitnitsky, A.R. Anomalous Axion Interactions and Topological Currents in Dense Matter. Phys. Rev. D 2005, 72, 045011. [Google Scholar] [CrossRef] [Green Version]
  38. Chernodub, M.N. Anomalous Transport Due to the Conformal Anomaly. Phys. Rev. Lett. 2016, 117, 141601. [Google Scholar] [CrossRef]
  39. Hasan, M.Z.; Kane, C.L. Topological Insulators. Rev. Mod. Phys. 2010, 82, 3045–3067. [Google Scholar] [CrossRef] [Green Version]
  40. Qi, X.L.; Zhang, S.C. Topological insulators and superconductors. Rev. Mod. Phys. 2011, 83, 1057–1110. [Google Scholar] [CrossRef] [Green Version]
  41. Volovik, G.E. Topology of Quantum Vacuum; Draft for Chapter in Proceedings of the Como Summer School on Analogue Gravity, Lecture Notes in Physics; Cornell University: Ithac, NY, USA, 2013. [Google Scholar]
  42. Volovik, G.E. Quantum phase transitions from topology in momentum space. In Quantum Analogues: From Phase Transitions to Black Holes and Cosmology; Unruh, W.G., Schutzhold, R., Eds.; Springer Lecture Notes in Physics; Springer: Berlin/Heidelberg, Germany, 2007; Volume 718, pp. 31–73. [Google Scholar]
  43. Volovik, G.E. Topological invariants for Standard Model: From semi-metal to topological insulator. JETP Lett. 2010, 91, 55. [Google Scholar] [CrossRef] [Green Version]
  44. Gurarie, V. Single-particle Green’s functions and interacting topological insulators. Phys. Rev. B 2011, 83, 085426. [Google Scholar] [CrossRef] [Green Version]
  45. Essin, A.M.; Gurarie, V. Bulk-boundary correspondence of topological insulators from their Green’s functions. Phys. Rev. B 2011, 84, 125132. [Google Scholar] [CrossRef] [Green Version]
  46. Volovik, G.E. The Universe in a Helium Droplet; Clarendon Press: Oxford, UK, 2003. [Google Scholar]
  47. Volovik, G.E. Topological Superfluids. arXiv 2016, arXiv:1602.02595. [Google Scholar]
  48. Nielsen, H.B.; Ninomiya, M. Absence of neutrinos on a lattice: I-Proof by homotopy theory. Nucl. Phys. B 1981, 185, 20. [Google Scholar] [CrossRef]
  49. So, H. Induced topological invariants by lattice fermions in odd dimensions. Prog. Theor. Phys. 1985, 74, 585–593. [Google Scholar] [CrossRef] [Green Version]
  50. Ishikawa, K.; Matsuyama, T. Magnetic field induced multi component QED in three-dimensions and quantum Hall effect. Z. Phys. C 1986, 33, 41–45. [Google Scholar] [CrossRef]
  51. Kaplan, D.B. Method for simulating chiral fermions on the lattice. Phys. Lett. B 1992, 288, 342–347. [Google Scholar] [CrossRef] [Green Version]
  52. Golterman, M.F.L.; Jansen, K.; Kaplan, D.B. Chern-Simons currents and chiral fermions on the lattice. Phys. Lett. B 1993, 301, 219–223. [Google Scholar] [CrossRef] [Green Version]
  53. Hořava, P. Stability of Fermi surfaces and K-theory. Phys. Rev. Lett. 2005, 95, 016405. [Google Scholar] [CrossRef] [Green Version]
  54. Creutz, M. Four-dimensional graphene and chiral fermions. JHEP 2008, 2008, 017. [Google Scholar] [CrossRef]
  55. Kaplan, D.B.; Sun, S. Spacetime as a topological insulator. arXiv 2011, arXiv:1112.0302. [Google Scholar]
  56. Suleymanov, M.; Zubkov, M.A. Wigner-Weyl formalism and the propagator of Wilson fermions in the presence of varying external electromagnetic field. Nucl. Phys. B 2019, 938, 171. [Google Scholar] [CrossRef]
  57. Zubkov, M.A.; Wu, X. Topological invariant in terms of the Green’s functions for the Quantum Hall Effect in the presence of varying magnetic field. arXiv 2019, arXiv:1901.06661. [Google Scholar]
  58. Zhang, C.X.; Zubkov, M.A. Influence of interactions on the anomalous quantum Hall effect. arXiv 2019, arXiv:1902.06545. [Google Scholar]
  59. Matsuyama, T. Quantization of Conductivity Induced by Topological Structure of Energy Momentum Space in Generalized QED in Three-dimensions. Prog. Theor. Phys. 1987, 77, 711. [Google Scholar] [CrossRef] [Green Version]
  60. Volovik, G.E. An analog of the quantum Hall effect in a superfluid 3 He film. Sov. Phys. JETP 1988, 67, 1804–1811. [Google Scholar]
  61. Thouless, D.J.; Kohmoto, M.; Nightingale, M.P.; den Nijs, M. Quantized Hall Conductance in a Two-Dimensional Periodic Potential. Phys. Rev. Lett. 1982, 49, 405–408. [Google Scholar] [CrossRef] [Green Version]
  62. Qi, X.L.; Hughes, T.L.; Zhang, S.C. Topological field theory of time-reversal invariant insulators. Phys. Rev. B 2008, 78, 195424. [Google Scholar] [CrossRef] [Green Version]
  63. Fradkin, E. Field Theories of Condensed Matter Physics, 2nd ed.; Cambridge University Press: Cambridge, UK, 2013. [Google Scholar]
  64. Vazifeh, M.; Franz, M. Electromagnetic response of Weyl semimetals. Phys. Rev. Lett. 2013, 111, 027201. [Google Scholar] [CrossRef] [Green Version]
  65. Chen, Y.; Wu, S.; Burkov, A. Axion response in Weyl semimetals. Phys. Rev. B 2013, 88, 125105. [Google Scholar] [CrossRef] [Green Version]
  66. Chen, Y.; Bergman, D.; Burkov, A. Weyl fermions and the anomalous Hall effect in metallic ferromagnets. Phys. Rev. B 2013, 88, 125110. [Google Scholar] [CrossRef] [Green Version]
  67. Ramamurthy, S.T.; Hughes, T.L. Patterns of electro-magnetic response in topological semi-metals. Phys. Rev. B 2015, 92, 085105. [Google Scholar] [CrossRef] [Green Version]
  68. Zyuzin, A.A.; Burkov, A.A. Topological response in Weyl semimetals and the chiral anomaly. Phys. Rev. B 2012, 86, 115133. [Google Scholar] [CrossRef] [Green Version]
  69. Goswami, P.; Tewari, S. Axionic field theory of (3+1)-dimensional Weyl semi-metals. Phys. Rev. B 2013, 88, 24510. [Google Scholar] [CrossRef] [Green Version]
  70. Niu, Q.; Thouless, D.J.; Wu, Y.S. Quantized Hall conductance as a topological invariant. Phys. Rev. B 1985, 31, 3372. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Shiba, H.; Kanada, K.; Hasegawa, H.; Fukuyama, H. Galvanomagnetic Effects in Impurity Band Conductions. J. Phys. Soc. Jpn. 1971, 30, 972–987. [Google Scholar] [CrossRef]
  72. Tong, D. Lectures on the Quantum Hall Effect. arXiv 2016, arXiv:1606.06687. [Google Scholar]
  73. Hatsugai, Y. Topological aspects of the quantum Hall effect. J. Phys. Condens. Matter 1997, 9, 2507–2549. [Google Scholar] [CrossRef]
  74. de Juan, F.; Manes, J.L.; Vozmediano, M.A.H. Gauge fields from strain in graphene. Phys. Rev. B 2013, 87, 165131. [Google Scholar] [CrossRef] [Green Version]
  75. Vozmediano, M.A.H.; Katsnelson, M.I.; Guinea, F. Gauge fields in graphene. Phys. Rep. 2010, 496, 109. [Google Scholar] [CrossRef] [Green Version]
  76. Cortijo, A.; Guinea, F.; Vozmediano, M.A.H. Geometrical and topological aspects of graphene and related materials. J. Phys. A: Math. Theor 2012, 45, 383001. [Google Scholar] [CrossRef] [Green Version]
  77. Manes, J.L.; de Juan, F.; Sturla, M.; Vozmediano, M.A.H. Generalized effective Hamiltonian for graphene under nonuniform strain. Phys. Rev. B 2013, 88, 155405. [Google Scholar] [CrossRef] [Green Version]
  78. Oliva-Leyva, M.; Naumis, G.G. Understanding electron behavior in strained graphene as a reciprocal space distortion. Phys. Rev. B 2013, 88, 085430. [Google Scholar] [CrossRef] [Green Version]
  79. Volovik, G.E.; Zubkov, M.A. Emergent gravity in graphene. In Talk Presented at the International Moscow Phenomenology Workshop; Cornell University: Ithac, NY, USA, 2015. [Google Scholar]
  80. Volovik, G.E.; Zubkov, M.A. Emergent Hořava gravity in graphene. Ann. Phys. 2014, 340, 352. [Google Scholar] [CrossRef] [Green Version]
  81. Volovik, G.E.; Zubkov, M.A. Emergent geometry experienced by fermions in graphene in the presence of dislocations. Ann. Phys. 2015, 356, 255. [Google Scholar] [CrossRef] [Green Version]
  82. Khaidukov, Z.V.; Zubkov, M.A. Landau Levels in graphene in the presence of emergent gravity. Eur. Phys. J. B 2016, 89, 213. [Google Scholar] [CrossRef] [Green Version]
  83. Fialkovsky, I.V.; Vassilevich, D.V. Graphene through the looking glass of QFT. Mod. Phys. Lett. A 2016, 31, 1630047. [Google Scholar] [CrossRef] [Green Version]
  84. Fialkovsky, I.V.; Vassilevich, D.V. Quantum Field Theory in Graphene. Int. J. Mod. Phys. A 2012, 27, 1260007. [Google Scholar] [CrossRef] [Green Version]
  85. Volovik, G.E.; Zubkov, M.A. Emergent Weyl spinors in multi-fermion systems. Nucl. Phys. B 2014, 881, 514. [Google Scholar] [CrossRef]
  86. Shapourian, H.; Hughes, T.L.; Ryu, S. The viscoelastic response of topological tight-binding models in 2d and 3d. Phys. Rev. B 2015, 92, 165131. [Google Scholar] [CrossRef] [Green Version]
  87. Cortijo, A.; Zubkov, M.A. Emergent gravity in the cubic tight-binding model of Weyl semimetal in the presence of elastic deformations. Ann. Phys. 2016, 366, 45. [Google Scholar] [CrossRef] [Green Version]
  88. Chernodub, M.N.; Zubkov, M. Chiral anomaly in Dirac semimetals due to dislocations. Phys. Rev. B 2017, 95, 115410. [Google Scholar] [CrossRef] [Green Version]
  89. Zubkov, M.A. Emergent gravity and chiral anomaly in Dirac semimetals in the presence of dislocations. Ann. Phys. 2015, 360, 655. [Google Scholar] [CrossRef] [Green Version]
  90. Sheng, D.N.; Sheng, L.; Weng, Z.Y. Quantum Hall effect in graphene: Disorder effect and phase diagram. Phys. Rev. B 2006, 73, 233406. [Google Scholar] [CrossRef] [Green Version]
  91. Hatsugai, Y.; Fukui, T.; Aoki, H. Topological analysis of the quantum Hall effect in graphene: Dirac-Fermi transition across van Hove singularities and edge versus bulk quantum numbers. Phys. Rev. B 2006, 74, 20541. [Google Scholar] [CrossRef] [Green Version]

Share and Cite

MDPI and ACS Style

Fialkovsky, I.V.; Zubkov, M.A. Elastic Deformations and Wigner–Weyl Formalism in Graphene. Symmetry 2020, 12, 317. https://doi.org/10.3390/sym12020317

AMA Style

Fialkovsky IV, Zubkov MA. Elastic Deformations and Wigner–Weyl Formalism in Graphene. Symmetry. 2020; 12(2):317. https://doi.org/10.3390/sym12020317

Chicago/Turabian Style

Fialkovsky, I.V., and M.A. Zubkov. 2020. "Elastic Deformations and Wigner–Weyl Formalism in Graphene" Symmetry 12, no. 2: 317. https://doi.org/10.3390/sym12020317

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop