The role of the strain induced population imbalance in Valley polarization of graphene: Berry curvature perspective

Real magnetic and lattice deformation gauge fields have been investigated in honeycomb lattice of graphene. The coexistence of these two gauges will induce a gap difference between two valley points (K and K′) of system. This gap difference allows us to study the possible topological valley Hall current and valley polarization in the graphene sheet. In the absence of magnetic field, the strain alone could not generate a valley polarization when the Fermi energy coincides exactly with the Dirac points. Since in this case there is not any imbalance between the population of the valley points. In other words each of these gauges alone could not induce any topological valley-polarized current in the system at zero Fermi energy. Meanwhile at non-zero Fermi energies population imbalance can be generated as a result of the external strain even at zero magnetic field. In the context of Berry curvature within the linear response regime the valley polarization (both magnetic free polarization, Π0, and field dependent response function, χ α) in different values of gauge fields of lattice deformation has been obtained.

Real magnetic and lattice deformation gauge fields have been investigated in honeycomb lattice of graphene. The coexistence of these two gauges will induce a gap difference between two valley points (K and K′) of system. This gap difference allows us to study the possible topological valley Hall current and valley polarization in the graphene sheet. In the absence of magnetic field, the strain alone could not generate a valley polarization when the Fermi energy coincides exactly with the Dirac points. Since in this case there is not any imbalance between the population of the valley points. In other words each of these gauges alone could not induce any topological valley-polarized current in the system at zero Fermi energy. Meanwhile at non-zero Fermi energies population imbalance can be generated as a result of the external strain even at zero magnetic field. In the context of Berry curvature within the linear response regime the valley polarization (both magnetic free polarization, Π 0 , and field dependent response function, χ α ) in different values of gauge fields of lattice deformation has been obtained.
Valleytronics, the valley version of spintronics is based on quantum valley number which carries the information by valley degree of freedom 1 . The family of two dimensional materials with hexagonal lattice structures are potentially a good candidate to manipulate valley polarization [1][2][3][4][5][6][7] . Recently the first observation of valley Hall effect in graphene family with inversion symmetry breaking has been reported 8,9 . Various proposals have been presented for generation of the valley Hall conductivity 10 . Massive Dirac electrons in graphene has been suggested for valleytronics applications 11,12 . In the literatures, most of the proposals for valley-filtering realization are combinations of strain induced gauge fields and electromagnetic fields [13][14][15][16][17] . A valley filter has been proposed regarding the fact that the effective magnetic field (that results from combination of the strain and external magnetic field) is large in one of the Dirac cones and can be zero in the other 17 . In the presence of Rashba spin orbit coupling and magnetic barrier in strained graphene, it is possible to produce valley and spin polarized currents 18 . Jiang and et al. proposed a scheme for generation of valley current in strained graphene in which the strain could be described by cyclic adiabatic deformations and there is a chemical potential in suspended region 19 . Wang and et al. have used the method of adiabatic quantum pumping in a three barrier structure with strained graphene and a ferromagnetic layer to generate pure valley current 20 .
The heart of the field of straintronic in two dimensional materials is to manipulate the valley polarization by applying various deformations. The possibility of strain induced pseudomagnetic field in graphene family up to 300T is a remarkable result and getting attention to strain engineering in graphene-like systems [21][22][23] . Further challenges for generating pseudomagnetic fields in graphene and other two dimensional materials not restricted to locally strained graphene nanobubbles. There is recently another proposal by Zhu and coworkers in which a pseudomagnetic field around 200T has been obtained by an uniaxial stretched graphene 23 . Arias and coworkers considered the possibility of gauge field generation from an elastic deformation in graphene by a quantum field theory approach. They found a relation between the pseudomagnetic field and Riemann curvature 24 . Vaezi and coworkers by means of time dependent strain and consequently a time dependent gauge field showed that the charge current could be generated due to the time dependent elastic deformations in graphene 25 . By specific configurations of external potential 26 or strain 27 , it is possible to access valley polarized current in bilayer graphene. In addition, the valley dependent Hall transport in the presence of electric field was proposed before 7 . In other graphene-like structures such as silicene highly efficient spin-valley filtering has been reported in the presence of magnetic exchange interaction and external strain 28 . Meanwhile pure valley-polarized current in monolayer graphene has been obtained using uniaxial strain within a region that separated by magnetic barriers 29 . It has also been reported that spin-dependent magnetic gap could be generated when the magnetic exchange field is assumed to be different at A-and B-sublattices of silicene. Where in this case the spin-valley polarization can be controlled by an external electric field 30 . In this letter we propose a device based on magnetic field and typical nonuniform strain which could be applied as depicted in Fig. 1(a-c). We study the valley polarization of the sample in the context of Berry curvature. By considering the gap difference between two valley points (K and K′) the possible topological valley Hall current and valley polarization are studied.

Methodology
Both strain and magnetic field could be realized as experimental tuning parameters [31][32][33][34] . Although the magnetic field breaks the time reversal symmetry, however, it cannot lift the valley degeneracy. Simultaneous presence of magnetic field and strain in gapped graphene leads to a gap difference between two valleys. Despite the various proposals where have been based on combinations of strain, magnetic and electric fields [13][14][15][16]18 in the current study, the proposed sample is based on a magnetic field and a nonuniform strain. The valley-resolved Hamiltonian in k-space reads where v f is Fermi velocity, the valley index η is +(−) for K(K′) and the matrices of τ act on sub-lattice indices (A and B).The strain and real magnetic field contributions in k-space given by , m = 5 meV, eV 0 0206 ′ the gap difference is ΔEVGD = 0.041 eV. For more information about two gauges of strain A S and magnetic field A M see the supplementary information.
where A S is the gauge field due to the strain and A M is the real magnetic gauge fields (see supplementary information). we have defined ( ) stands for the magnetic field induced momentum transfer of q. Finally the k-space mass term is expressed as From theoretical point of view a gap opening in the presence of non-uniform strain is 35 , y xx y y x xy z . The energy gap as function of pseudomagnetic field is the Zeeman coupling of pseudospin to the associated pseudomagnetic field 35 , E Zeeman = (3)/(8)V′a 2 B ps where a is the lattice constant, and V′ = 6(eV/Å) 36 .
It should be noted that the nature of the magnetic and strain dependent gauges are quite different. The magnetic gauge field depends on electronic positions, however, the strain induced gauge is a function of the atomic positions (u α ) that cannot act on electronic states. Accordingly results of the present approach can simply be generalized to uniform strains. Unlike the other contributing terms in the Hamiltonian since the magnetic vector potential depends on electrons position in real spac ™ e, therefore, the Hamiltonian of magnetic field, as indicated in Eq. (3), contains momentum transfer contributions and cannot be represented in block diagonal form of independent k subspaces. However, within the first order perturbation approach, in which the momentum transferring terms cannot contribute, it can be shown that the previous results could be considered reliable at the qualitative level 18 . Unlike the first order correction, momentum transferring terms contribute in the higher orders of the perturbation and generally used block diagonalization approach cannot give rise to an exact answer. Then the non-perturbative part of the strained graphene Hamiltonian is The matrix representation of H (0) in the following k-space basis , an be written as . First order perturbation accounts for the zero momentum transfer contribution of the magnetic field. Therefore it can be described as It should be noted that as mentioned before H (1) corresponds to the first order corrections and the higher orders of perturbative corrections can be achieved by taking into account the non-zero momentum transfer contributions. The gap difference between two valleys, K and K′, can be obtained within the perturbation theory (see supplementary information). It can be shown that in strained graphene this valley gap difference depends on both magnetic field and applied strain. Within the first order perturbation valley gap difference, which has been defined as E E E The external magnetic field is denoted by B and L stands for the size of the system. The α and β are the components of the unperturbed Dirac Hamiltonian eigenvectors, (0) ψ η and (see supplementary information for more details), As illustrated in Fig. 1d and e the valley gap difference has been appeared as a result of the magnetic field and strain. It is very important to note that Eq. (12) indicates that for a given k ≠ 0 state we have K k ≠ ± ′ ±   even at zero magnetic field. This means that there is a strain induced population imbalance between different valleys at B = 0. It has also been realized that the strain moves the Dirac points slightly. Meanwhile, this movement does not change the direct band gap of graphene. A valley Hall current can be generated by the population imbalance as a consequence of the valley gap difference. We have proposed a practical method for measuring the valley Hall current in Fig. 1c that could be available using the present experimental techniques. By means of this setup, we demonstrated that the valley polarized current can be achieved in strained graphene in the presence of magnetic field which can be expressed in terms of Berry curvature. The topological response to an external gauge field could be obtained by integrating of Berry curvature of filled bands over the momentum space 37 The Berry curvature k k Dependence of Berry curvature to magnetic and pseudo magnetic fields is shown in Fig. 2a,b and c.

Results and Discussion
In the absence of the strain, there is not any gap difference between two valleys ( = ′ E E K gap K gap ). Using the relations which were presented in Eqs 13 and 14 it can be inferred that the charge conductivity of inequivalent valleys satisfy xy KC xy K C σ σ = − ′ . Where σ η xy C is the contribution of the η-valley in the charge Hall conductivity. Meanwhile, the valley resolved current with a given valley index, η, can be defined by where v is the velocity operator.
where j C is the charge current. The valley polarization in the presence of the external magnetic field is then given by M As long as the population balance between the two valley points is maintained, valley current cannot be induced in the sample. Valley population imbalance can be achieved when the chemical potential of inequivalent Dirac points are not the same where for the strained graphene and in the presence of magnetic field this is really the case. It should be mentioned that valley population imbalance has been employed for generation of valley current in the presence of magnetic field in strained graphene as described in ref. 18 . Furthermore, realization of the valley Hall current strongly depends on population imbalance at different valleys 7 . Accordingly, the population imbalance that comes from the coexistence of real and pseudo-magnetic fields 18 automatically guaranties the generation of valley polarization in strained valley Hall systems. In the present work it was shown that when the Fermi energy has not exactly been located at Dirac points the strain itself could results in population imbalance and magnetic free valley polarization. Since the Berry curvatures of different valleys have opposite sign then the Berry curvature dependent quantities such as Hall conductivity vanishes identically once the population balance is established. This is the key point for generation of topological valley current in the system. Using the Eq. 14 it can be shown that, in the absence of the strain, Berry curvature of different Dirac points in the presence of an external magnetic field are opposite. In the other words the total charge conductivity vanishes exactly i.e. This gap difference leads population imbalance in two inequivalent Dirac points. Meanwhile unlike the previous case, the Berry curvature of each valleys doesn't contribute oppositely and therefore σ σ In the absence of the magnetic field the total conductivity vanishes just when the Fermi energy is exactly located at the Dirac points (which has been indicated by k x = 0 and k y = 0 in Eq. 14 which indicates the contribution of the Dirac points). In Fig. 3 the valley polarization of strained graphene has been depicted in terms of magnetic field. By increasing the value of magnetic fields, functionality of valley polarization for each value of B ps changes in different way. It should be noted that the zero magnetic polarization which has been shown in this figure comes from the strain induced population imbalance in nonzero Fermi circle as discussed before. We have chosen E 0 1 F  . eV that corresponds to n-type doping with conduction electron density about Π is the magnetic free polarization and  χ α is linear response function of the strained system. After expansion, x y k k s 0 Within a numerical calculation it can be shown that, depending on the value of the Fermi energy, the magnetic free valley polarization Π  0 cannot go all the way to zero. As discussed before this can be generalized to the uniform strains and calculations show that this is really the case. Figure 4 (inset) shows the magnetic free polarization in terms of pseudo magnetic field. The linear response function of the magnetic field induced polarization, χ α , in terms of the B ps has also been illustrated in Fig. 4. It should be mentioned that because of the symmetric form of magnetic gauge field within the Dirac cone approximation, in which the band anisotropy of the sample has been ignored, the behavior of χ y is expected to be identical with χ x . Strong magnetic field modifies the electronic spectrum and forms landau levels. In this case, the magnetic field itself induces valley density polarization in the presence of non-vanishing intrinsic orbital magnetic moment 38 . Since the landau levels of two-dimensional  structures are localized states, therefore this polarized density could not be extracted form the system. However the external magnetic field here in the present case is weak enough that cannot change the spectrum significantly. In addition we have shown that the population imbalance can be induced by deformation gauge even at zero magnetic field. Berry curvature dependent quantities can arise provided that the inversion symmetry is broken. Meanwhile since the contribution of each valley has been canceled out by its counterpart there is no valley polarization in the graphene based Hall current system. Meanwhile as discussed in the current study strain could induce population imbalance between different Dirac points when the radius of the Fermi circle is not zero. Since the valley degeneracy has been lifted by the strain for k ≠ 0 even at B = 0. This is due to the fact that in the presence of the strain we have ≠ , (0)   . This means that when the energy of the Fermi level does not exactly meet the Dirac points there is a population imbalance between two nonequivalent Dirac cones. Therefore in this case valley polarized Berry curvature dependent quantities, such as valley polarized Hall conductivity, could be realized as a result of this population imbalance.