V. N.
Davydov
M. V. Lomonosov Moscow State University, Leninsky pr. 71, app. 121, 117296 Moscow, Russia. E-mail: vladimirdavydov@yahoo.com
First published on 22nd July 2020
The origin of the Lifshitz topological transition (LTT) and the 2D nature of the LTT in graphene has been established. The peculiarities of the Lifshitz topological transitions in graphene are described at the Brillouin zone centre Γ, the zone corners K in the vicinity of the Dirac points, and at the saddle point M. A general formulation of the thermodynamics at the LTT in graphene is given. The thermodynamic characteristics of graphene are investigated at the Lifshitz topological transitions. Anomalies are found in the electron specific heat Ce, the electron thermal coefficient of pressure, and the coefficients of electron compressibility and thermal expansion in graphene at the LTT. All the thermodynamic parameters possess the strongest singularities in graphene at the LTT near the saddle points. This opens exciting opportunities for inducing and exploring the Lifshitz topological transitions in graphene.
δΩ = −α|z|5/2. | (1.1) |
It can be seen that the third derivatives of the thermodynamic potential become infinite at the point z = μ − εc = 0 (the chemical potential μ is equal to the Fermi energy εF at T = 0); therefore, this change in the topology of the Fermi surface is also called the 2½-order phase transition (FT2½) according to the Ehrenfest terminology.6,7
The anomalies arising in the LTT manifest themselves at T ≪ εF not only in the thermodynamic characteristics of metals but also in other characteristics (e.g., magnetic-field dependence of the electrical resistance,8 sound absorption9–15). It has been shown in ref. 16 that the thermoelectric power has the square root divergency at the LTT. The discovery of graphene17 gave new inspiration to investigations of the Lifshitz transition properties. To date, numerous papers have been devoted to the Lifshitz transitions in bilayer (BLG) and multilayer (MLG) graphene.18–30 The results of these studies indicate that the effects of the LTT are appreciable and can be used to observe the 2½-order transitions as well as to investigate the degree of smearing of the Fermi surface in metals. In ref. 31, experimental evidence was obtained of the Lifshitz transition in the thermoelectric power of ultrahigh mobility bilayer graphene. Resolving low-energy features in the density of states (DOS) holds the key to understanding a wide variety of rich and novel phenomena in graphene based on 2D heterostructures. The Lifshitz transition in bilayer graphene (BLG) arising from trigonal warping has been established in31 theoretically and experimentally.
The 21st century has brought many new results related to graphene thermodynamics.32–45 Apparently, the thermodynamics of graphene has been characterized from many points of view; however, the role of the Lifshitz topological transitions in the thermodynamic properties of graphene has not been studied. This paper aims to close this gap. Because graphene is the first real two-dimensional solid, a general formulation of the thermodynamics at the LTT in graphene is given. The unusual thermodynamic properties of graphene stem from its 2D nature, forming a rich playground for new discoveries of heat flow physics and potentially leading to novel thermal management applications.
This paper is arranged as follows. In section 2, the origin of the Lifshitz topological transition and the 2D nature of the LTT in graphene are considered. Peculiarities of the Lifshitz topological transitions in graphene are then investigated in section 3. The thermodynamics of graphene at the Lifshitz topological transitions is proposed in section 4. Finally, conclusions are drawn in section 5.
(2.1) |
(2.2) |
At the neck rupture, the boundary isoenergetic surface ε(p) = εc contains the peculiar point of another type, named the cone point, at p = pc. In this case, the isoenergetic surface in the vicinity of the cone point p = pc is described by the hyperboloid of two sheets at energy ε < εc and the hyperboloid of one sheet at energy ε > εc (Fig. 1a):
(2.3) |
(2.4) |
Fig. 1 The cone point at p = pc (a), the saddle point (c), and the van Hove singularities of the state density δD(ε) (b). S, S1, and S2 are the saddle points. |
At energies close to εc, one can express the electron density of states (DOS) as
D(ε) = D0(ε) + δD(ε) | (2.5) |
(2.6) |
One can join eqn (2.1)–(2.4) into a single equation:
(2.7) |
Four types of singularities exist of the density of the electron states in three-dimensional space, and the singularity type depends on the signs of the coefficients aj.
The point M0 (min) corresponds to the case where all three coefficients aj = +1. This point corresponds to the minimum in the energy spectrum.
The point M1 (saddle) corresponds to the case where two coefficients aj are positive and the third one is negative. This is the saddle point (Fig. 1c).
The point M2 (saddle) is the case where two coefficients aj are negative and a third one is positive. This is the second saddle point (Fig. 1b).
The point M3 (max) corresponds to the case where all three coefficients aj = −1. This point corresponds to the maximum in the energy spectrum.
There are three types of singularities of the density of the electron states in two-dimensional space, and the singularity type depends on the signs of the coefficients aj (j = x, y).
The point P0 (min) corresponds to the case where both coefficients aj = +1. This point corresponds to the minimum in the energy spectrum.
The point P1 (saddle) corresponds to the case where one coefficient aj is positive and another one is negative. This is the saddle point.
The point P2 (max) corresponds to the case where both coefficients aj = −1. This point corresponds to the maximum in the energy spectrum.
There are two types of singularities of the density of the electron states in one-dimensional space:
• The point Q0 (min) corresponds to the case where coefficient a = +1. This point corresponds to the minimum in the energy spectrum.
• The point Q1 (max) corresponds to the case where coefficient a = −1. This point corresponds to the maximum in the energy spectrum.
The analytical behavior of the density of the electron states at LTT is given in Table 1. The results are computed from formulae (2.1)–(2.7); see also ref. 1, 2, and 4.
We can conclude based on the square root peculiarities of the density of the electron states in three dimensions (Table 1) at εc that nature of the Lifshitz topological transition stems from the van Hove singularities (VHS) of the state density δD(ε).48–52 In this connection, the van Hove topological theorem48 states that the spectrum must contain at least one of the saddle-points S1 and S2 (Fig. 1b), and the slope of D(ε) must tend to −∞ on the upper end.
The two-dimensional nature of graphene should exhibit special types of van Hove singularities and LTTs. This statement is illustrated by the general argument that in two dimensions, the saddle-points produce logarithmic singularities (Table 1), and the spectrum extrema produce finite discontinuities of the electron density of states [also see ref. 49, 53, and 54]. This is valid for the elementary excitations of the quasiparticles with values mi ≠ 0 of the effective mass tensor; however, it is not applied to the massless Dirac fermions in graphene. If the logarithmic singularity is a general property of 2D electronic systems in the saddle points, this statement should be valid for 2D graphene. However, the latter is in contradiction with because the value of C6 must be zero for the massless fermions in graphene. Here, the discrepancy arises because the logarithmic singularity vanishes. We will resolve this contradiction in section 3.
Therefore, the Lifshitz topological transitions in graphene require special investigation. The latter was also confirmed by recent work.55 It has been demonstrated in55 theoretically that the characteristic feature of a 2D system undergoing N consequent Lifshitz topological transitions is the occurrence of spikes of entropy per particle s of a magnitude ± ln2/(J − 1/2) with 2 ≤ J ≤ N at low temperatures.
Fig. 2 (Reproduced and adapted with permission from ref. 56). The crystal structure of the single graphite layer (a). The primitive translation vectors e1 and e2 form the rhombic unit cell, and the basis consists of two C atoms, shown as A and B. The Bravais lattice (consider, e.g., the lattice formed by the A atoms) is triangular with a Bravais lattice spacing of where a is the spacing between neighboring atoms. The graphene reciprocal lattice and the first Brillouin zone (b). |
From Fig. 2a, we see that two basis vectors of the Bravais lattice are
(3.1) |
The reciprocal Bravais lattices, b1 and b2, are defined such that
bi·ej = 2πδij | (3.2) |
The result is
(3.3) |
(3.4) |
These basis vectors are of equal length and at 120°; therefore, the reciprocal lattice is a hexagonal lattice (Fig. 2b). The first Brillouin zone is shown in Fig. 2b. The first Brillouin zone is a regular hexagon, whose most characteristic points are its centre Γ, the inequivalent corners K and K′, and the centres of the lateral edges M and M'. The distance ΓM is The distance ΓΚ is
Consider a state with amplitude ΨR for the electron to be at the site labeled by R. We look for the wave functions with amplitudes which vary, such as eik·r. There will be different amplitudes ΨA and ΨB on sublattices A and B, so ΨR = ΨAeik·r (R ∈ A), and ΨR = ΨBeik·r (R ∈ B). An electron at site R can hop to any of the neighboring sites. An atom on sublattice AAA has neighboring atoms (Fig. 2a), all on sublattice BBB, at displacements (0, a), and An atom on sublattice BBB has three neighbors on sublattice AAA at displacements (0, −a), and In the nearest-neighbour approximation, there are no hopping processes within the sublattices AAA and BBB; hopping only occurs between them. Hence, the eigenvalue ε and the amplitudes ΨA and ΨB are determined for each wavevector k from two equations:
(3.5) |
(3.6) |
We write the Hamiltonian in the tight-binding approximation for the wave vector k = p/ħ (where p is the electron quasimomentum), taking into account eqn (3.5) and (3.6), as
(3.7) |
(3.8) |
The electron energy eigenvalues of the Hamiltonian (3.7) are given by
(3.9) |
In expression (3.9), the plus and minus signs correspond to the conduction and valence bands, respectively.
To explore the possible Lifshitz topological transitions, consider the structure of the isoenergetic lines in the (kx, ky) plane for different positions of the wave vector k.
(i) The maximum and minimum energies are εc = ±3γ0, and they occur at k = 0, i.e. in the centre Γ (Fig. 2b) of the first Brillouin zone. The value of εc = 3γ0 corresponds to the maximum of energy in the conduction band (the point P2 (max) in Table 1). The value of εc = −3γ0 corresponds to the minimum of energy in the valence band (the point Po (min) in Table 1). Taking into account |kya|, |kxa| ≪ 1, from (3.9), we have accuracy of order (kya)2, (kxa)2:
(3.10) |
Thus, the isoenergetic lines are circles in the vicinity of the extrema of energy at k = 0. It follows from (3.10) that these circles are described by the following equation in p-space:
(3.11) |
(3.12) |
Consequently, there exist LTT at critical values of energy εc = ±3γ0, where the cavity (3.11) of the isoenergetic lines disappears (Fig. 3a) at εc = +3γ0 in the conduction band or appears at εc = −3γ0 in the valence band.
The number of the electron states inside the two-dimensional cavity (per spin direction and per unit area) is equal to
δN(ε) = (1/2πħ)2Δ(ε) | (3.13) |
We obtain from (3.11) and (3.13)
(3.14) |
One can represent the change of the density of states δD(ε) due to LTT as
(3.15) |
The discontinuity at the end points ε = ±3γ0 comes from a maximum or minimum of the dispersion relation in two dimensions. The analytic expression without derivation for the total DOS in graphene for model (3.7) was provided by Hobson and Nierenberg in 1953.53 In a recent paper,54 the expression for this model was derived for the total DOS (per unit hexagonal cell of area (Fig. 2a) and one spin component), valid for region 0 < |ε| < 3γ0:
(3.16) |
(3.17) |
(3.18) |
Expanding of (3.16) near the point Γ in the small region |ε| ≲ 3γ0 in the vicinity of the LTT points |εc| = 3γ0, one obtains (per spin direction and per unit area):
(3.19) |
The result (3.15) for δD (ε) coincides with the total D(ε) given by (3.19) in the region |ε| ≲ 3γ0 in the vicinity of the LTT points |εc| = 3γ0, where at |ε| ≲ 3γ0.
As follows from (3.19), for graphene in the vicinity of the point Γ, the coefficients C5 and C7 in Table 1 should be changed as
(3.20) |
Hence, in 2D graphene, the mass-factors of coefficients C5 and C7 in Table 1 acquire the substitutions of γ0/VF2 instead of respectively, and this substitution can be treated as the fermion effective mass in the vicinity of the point Γ:
(3.21) |
k = K + κ | (3.22) |
To the lowest order in κ we have, from (3.8) and (3.22),
(3.23) |
Substitution of (3.23) into Hamiltonian (3.7) gives the following equation for the isoenergetic lines in vicinity of the point K (the corner of the first Brillouin zone):
(3.24) |
Consequently, the isoenergetic lines are circles in the vicinity of the minimum of energy εc = 0 (for the conduction band) or in the vicinity of the maximum of energy εc = 0 (for the valence band) at the corner of the zone. It follows from (3.24) that these circles are described by the following equation in p-space (Fig. 3b):
(3.25) |
The value of εc = 0 corresponds to the minimum of energy in the conduction band (the point P0 (min) in Table 1). The value of εc = 0 corresponds to the maximum of energy in the valence band (the point P2 (max) in Table 1). Eqn (3.25) describes the linear dispersion of the massless Dirac fermions in the vicinity of the corners K(K′) of the first Brillouin zone:
ε = ±VFp | (3.26) |
Thus, the LTT exists at the critical value of energy εc = 0, where the new cavity (3.25) appears (Fig. 3b) in the conduction band or disappears in the valence band (Fig. 3b). Hence, there are six pockets of low energy excitations (Fig. 3b), one for each of the two inequivalent points K and K′ on the Brillouin zone boundary.
The area of the cavity inside the circle (3.25) is equal to.
(3.27) |
Using (3.13) and (3.27) to calculate the number of the electron states inside the two-dimensional cavity (3.25), one obtains (per spin direction and per unit area):
δN(ε) =πε2/[VF2(2πħ)2] | (3.28) |
(3.29) |
Expanding the expression (3.16) at ε = εc = 0, we obtain the same result for the total DOS (per spin direction and per unit area) in the vicinity of the zone corners K(K′):
D(ε) =|ε|/(2πħ2VF2) | (3.30) |
Therefore, in the points K(K′), we must make the following substitutions for the coefficients C6 and C7 for graphene in Table 1:
(3.31) |
(3.32) |
k = K + q | (3.33) |
(3.34) |
Taking into account |qya|, |qxa| ≪ 1, we have from (3.9), (3.33) and (3.34) to accuracy of order (qya), (qxa)
(3.35) |
Therefore, the isoenergetic lines are hyperboles in the vicinity of point M (the middle of the edge of the first Brillouin zone). One can say that middles of the edges of the first Brillouin zone are saddle points. We can consider the corresponding points similar to M (the middles of the edges of the first Brillouin zone) (Fig. 2b) as the saddle points or the “cone points” in two dimensions (points P1 (saddle) in Table 1). It follows from (3.35) that the isoenergetic lines are described by the following equations in p-space in the vicinity of these saddle points (Fig. 4):
(3.36) |
(3.37) |
(3.38) |
The Lifshitz topological transition is realized in the saddle point M by the variation of energy from ε < εc to ε > εc, and εc = γ0. In this transition, the isoenergetic lines in the vicinity of point M change from hyperbole (3.36) to hyperbole (3.37). Calculating the area Δ(ε) enclosed by hyperbole (3.36) or (3.37) and the corresponding boundary of the first Brillouin zone py = ħ/a, or px = ħ/a, one obtains:
(3.39) |
Now, we are able to resolve the contradiction with the coefficient used in Table 1. As follows from (3.39) and (3.13), the DOS change due to the LTT in the saddle point M (at energies close to εc) is equal to
(3.40) |
Also, the expansion of (3.16) yields the same VHS for the total DOS (per spin direction and per unit area) in the vicinity of the saddle point M, if |ε| → εc = γ0:
(3.41) |
Therefore, in the saddle point M, the coefficient C6 for graphene in Table 1 should be changed to
(3.42) |
Thereby, in 2D graphene, the mass-factor C6 in Table 1 acquires the substitution of γ0/VF2 instead of, and this substitution can be treated as the graphene effective mass in the vicinity of the saddle point M:
(3.43) |
It follows from (3.21) and (3.43) that the fermions are slowing down in the vicinities of the Γ and M points, becoming not massless but massive. A similar phenomenon was found in observations of Dirac node formation and mass acquisition in a topological crystalline insulator.57 Estimation of the values of mΓeffand mMeffyields (at γ0 = 4 eV, and VF= 108 cm s−1 (ref. 56)):
mΓeff ≈ 0.1me | (3.44) |
mMeff ≈ 0.1me | (3.45) |
Ω(μ, T) = Ω0 + δΩ | (4.1) |
(4.2) |
Introducing the Lifshitz parameter z = μ − εc in the case of the hyperbolic changes of the Fermi lines at the saddle point M, we obtain from (3.39) and (4.2), if T → 0:
(4.3) |
Transition from region I into region II corresponds to the appearance of a new cavity of the isoenergetic line ε(p) = μ, or to a decrease of the line connectivity.
Both formulae (4.3) are valid at T ≪ |z|. At absolute zero temperature, one has
(4.4) |
From comparison of (4.4) and (1.1), one can see the essential difference of δΩ given by (4.4) for the LTT in graphene and the LTT in three dimensions (1.1). The second derivative of δΩ diverges at the point z = 0 as ln(|z|/εc)in graphene, while only the third derivative of δΩ diverges at the point z = 0 as z−1/2 in three dimensions. As a result, completely different anomalies will be obtained for the thermodynamic parameters at the Lifshitz topological transitions in graphene.
When the electronic cavity disappears in the Brillouin zone centre Γ, one computes from (3.14) and (4.2) at T → 0:
(4.5) |
Considering the appearance of new cavities in the Brillouin zone corners K, we obtain from (3.28) and (4.2), if T → 0:
(4.6) |
Because the number of electrons in the conduction band is permanent (at least in vicinity of the point μ = εc), it is convenient to use the Helmholtz free energy potential F (T, S, N)2 instead of the grand thermodynamic potential Ω (T, S, μ). We assume that area S is the parameter connected with applied pressure, tensile stress, or other mechanical action. The critical energy is the function of the area, i.e. εc = εc(S), and the chemical potential μ is also a function of S because of constancy of the carriers:
N(μ, S) = N | (4.7) |
If Sc is the area at which the topology of the Fermi lines changes:
μ(Sc) = εc(Sc) | (4.8) |
According to (4.7) and (4.8), the value of |z| = |μ − εc| can be expressed via |S − Sc|:
|z| = η|S − Sc| | (4.9) |
The Helmholtz free energy potential F can be written in the form
F = F0 + δF | (4.10) |
One can readily be convinced that δF is quantitatively equal to the irregular contribution δΩ expressed in the variables S and T according to the Landau theorem5 about the small variations of the thermodynamic potentials due to small changes of some parameters of the solid state:
(δΩ)T,S,μ = (δF)T,S,N | (4.11) |
One can see that variations of δΩ given by the relations (4.3), (4.5), and (4.6) are small at T ≪ z in the vicinity of the Lifshitz electronic transition, where z → 0.
Thus, δF is given by eqn (4.3), (4.5), and (4.7), where one must set |z| = η|S − Sc|.
In the case of the hyperbolic changes of the Fermi lines at the saddle point M, we obtain from (4.3) and (4.9), if T → 0:
(4.12) |
When the electronic cavity disappears in the Brillouin zone centre, one obtains from (4.5) and (4.9) at T → 0:
(4.13) |
Considering the appearance of new cavities in the Brillouin zone corners, one obtains from (4.6) and (4.9), if T → 0:
(4.14) |
(4.15) |
We can explore the electron-specific heat peculiarities of monolayer graphene at the LTT based on expressions of the change of the Helmholtz free energy by eqn (4.12), (4.13), and (4.14).
In the case of the hyperbolic changes of the Fermi lines at the saddle point M, one obtains from (4.12), if T → 0:
(4.16) |
When the electronic cavity disappears in the Brillouin zone centre Γ, one computes from (4.13) at T → 0:
(4.17) |
Considering the appearance of new cavities in the Brillouin zone corners K, we obtain from (4.14), if T → 0:
(4.18) |
Thus, the formulae (4.16), (4.17) and (4.18) describe peculiarities of the specific heat in graphene at the LTT at the points M, Γ, and K of the first Brillouin zone, correspondingly (Fig. 5a).
The results for the specific heat (4.17) and (4.18) coincide with result (4.15) of the electron specific heat CMGeof the monolayer graphene (MG) obtained in paper.61 Indeed, substituting , and expressing T in Kelvins, one obtains from (4.17) and (4.18):
Ce ≈ SkB2μT/(VF2ħ2), kBT ≪ |z| ≪ μ | (4.19) |
Note that experimental observation of singularities (4.16) and (4.17) is possible under special conditions. The experimental value of the chemical potential μ of graphene does not exceed 1 eV,63 and the estimated value of the overlap integral γ0 is between 2 and 4 eV.64–66 Tuning of the chemical activity of graphene in a wide range could be achieved by the formation of a Moiré superstructure between the graphene and transition metal substrate.67 Thus, to observe the LTT at the graphene saddle point (where |z| = |μ − γ0| → 0), we must increase the chemical potential μ or decrease the overlap integral γ0 due to strain. Increasing μ can be realized by electron doping due to renormalization of the electron spectrum.68 It has been demonstrated in ref. 68 that the Coulomb interaction produces noticeable effects in increasing the graphene chemical potential μ at low temperature T ≪ μ, especially for high carrier concentrations of n0 > 1011 cm−2. The band structure of graphene has been determined under strain using density functional calculations.69 The ab initio band structure was then used to extract the best fit to the tight-binding hopping parameters used in a microscopic model of strained graphene.70 It was found that the hopping parameters may increase or decrease with increasing strain depending on the orientation of the applied stress. The fitted values were compared with an available parameterization for the dependence of the orbital overlap on the distance separating the two carbon atoms.
In case of the hyperbolic changes of the Fermi lines at the saddle point M, we obtain from (4.12) for the electron compressibility, if T → 0:
(4.20) |
When the electronic cavity disappears in the Brillouin zone centre, one obtains from (4.13) for the electron compressibility at T → 0:
(4.21) |
Considering the appearance of new cavities in the Brillouin zone corners, one obtains from (4.14) for the electron compressibility, if T → 0:
(4.22) |
As a result, the formulae (4.20), (4.21) and (4.22) describe peculiarities of the compressibility of graphene at the LTT at the points M, Γ, and K of the first Brillouin zone, correspondingly (Fig. 5b).
(4.23) |
When the electronic cavity disappears in the Brillouin zone centre, one obtains from (4.13) for the electron thermal coefficient of pressure at T → 0:
(4.24) |
Considering the appearance of new cavities in the Brillouin zone corners, one obtains from (4.14) for the electron thermal coefficient of pressure, if T → 0:
(4.25) |
The formulae (4.23), (4.24) and (4.25) describe peculiarities of the electron thermal coefficient of pressure in graphene at the LTT at the points M, Γ, and K of the first Brillouin zone, correspondingly (Fig. 5c).
In the case of the hyperbolic changes of the Fermi lines at the saddle point M, we obtain the change of the electron thermal expansion coefficient from (4.23) and (4.20), if T → 0:
(4.26) |
When the electronic cavity disappears in the Brillouin zone centre, one obtains from (4.24) and (4.21) for the electron thermal expansion coefficient, at T → 0:
(4.27) |
Considering the appearance of new cavities in the Brillouin zone corners, one obtains from (4.25) and (4.22) for the electron thermal expansion coefficient, if T → 0:
(4.28) |
It is interesting that the electron thermal expansion coefficient has a negative sign at T ≪ |z|, like the lattice thermal expansion coefficient of graphene (see ref. 72, 73, and 75).
The formulae (4.26), (4.27) and (4.28) describe the peculiarities of the electron thermal expansion coefficient in graphene at LTT at the points M, Γ, and K of the first Brillouin zone, correspondingly (Fig. 5d).
(i) Connection of the Lifshitz topological transition has been established with the van Hove singularities of the electron state density in graphene. There are three types of singularities of the density of the electron states in two dimensions. The point P0 (min) corresponds to the minimum in the energy spectrum. The point P1 (saddle) corresponds to the case where the Lifshitz topological transitions are realized by variation of the energy from ε < εc to ε >εc. The point P2 (max) corresponds to the maximum in the energy spectrum.
(ii) Peculiarities of the Lifshitz topological transitions in graphene are described at the Brillouin zone centre Γ, at the zone corners K, in the vicinity of the Dirac points, and at the saddle point M. It is found that LTT can be realized in the centre Γ at the critical energy value εc = 3γ0, where the cavity of the isoenergetic lines disappears. The existence of the LTT is shown at the critical energy value of εc = 0, where six pockets of low energy excitations appear, one for each of the two inequivalent Dirac points K and K′.
The Lifshitz topological transition is realized in the saddle point M by variation of the energy from ε < εc to ε > εc and at the critical value of energy εc = γ0. It is shown that the Dirac fermion slows down in the vicinities of the points Γ and M, becoming not massless but massive, and the values of the fermion effective mass in the vicinity of these points are mΓeff ≈ 0.1me, and mMeff ≈ 0.1me (where me is the free electron mass).
(iii) The thermodynamic characteristics of graphene were investigated at the Lifshitz topological transitions. A general formulation of the thermodynamics at the LTT in graphene is given. The anomalies are found at the LTT of the electron specific heat Ce, the electron compressibility , the electron thermal coefficient of pressure , and the electron thermal expansion coefficient . The anomalies are described in terms of the Lifshitz parameter z = μ − εc. The electron specific heat Ce diverges at the saddle point M as ln|z| and is proportional to |z| at points K and Γ of the first Brillouin zone. The electron compressibility diverges at the saddle point M as ln|z|, is proportional to |z| at points K, and takes a constant value at point Γ of the first Brillouin zone. The electron thermal coefficient of pressure diverges at the saddle point M as 1/|z|, and it becomes negative in the vicinities of points K and Γ of the first Brillouin zone. The electron thermal expansion coefficient diverges at the saddle point M as 1/(|z|ln|z|) becoming negative in the vicinity of point Γ, and it diverges as 1/|z| at point K of the first Brillouin zone.
One can conclude that all the thermodynamic parameters possess the strongest singularities in graphene at the LTT near the saddle points M. In 2D graphene, a saddle point M in the electronic band structure leads to a divergence in the density of states of the logarithmic-type van Hove singularities (VHS). This implies the possibility of experimental observation of the LTT by bringing the chemical potential μ and the VHS together. However, one cannot change the position of the VHS in the band structure. It is pointed out in section 4.2 that the accessible experimental value of the graphene chemical potential μ does not exceed 1 eV.63 Therefore, it is essential to tune μ through the VHS by chemical doping76,77 or by gating.17,57,78–81 In recent work,82 a simple technique of doping graphene by manipulating adsorbed impurities was reported, and a change in the electron mobility of 650% was observed. Also, it is worth paying experimental attention to the tuning of μ through the VHS by the following two methods. The first method is connected with deformation of the graphene monolayer to mimic twisted graphene. Rotation between two stacked graphene monolayers in twisted graphene83 can generate van Hove singularities, which can be brought arbitrarily close to the chemical potential μ by varying the angle of rotation.84 The second method consists of investigating the LTT in graphene under 3D high pressure.84 This opens exciting opportunities for inducing and exploring the Lifshitz topological transitions in graphene.
This journal is © The Royal Society of Chemistry 2020 |