Role of the Lifshitz topological transitions in the thermodynamic properties of graphene

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.


Introduction
Anomalies in the thermodynamic quantities of metals due to a change in the topology of the Fermi surface 3(p) ¼ 3 F are customarily called the Lifshitz topological transition (LTT) (I. M. Lifshitz; 1 see also ref. 2, 3 and 4). It is understood that a change in the topology of the Fermi boundary surface is due to the continuous variation of some parameter (e.g., pressure), as a result of which the difference z ¼ 3 F À 3 c (where 3 F is the Fermi energy and 3 c is the critical energy at which the topology of the constant energy surface changes) passes through zero continuously. This leads to a change in the connectivity of the Fermi surface (the appearance of a new cavity, the rupture of a connecting neck, etc.), and at absolute zero temperature T ¼ 0 K, the grand thermodynamic potential U (oen called the Landau free energy or Landau potential 5 ) acquires an irregular correction: (1.1) It can be seen that the third derivatives of the thermodynamic potential become innite at the point z ¼ m À 3 c ¼ 0 (the chemical potential m is equal to the Fermi energy 3 F at T ¼ 0); therefore, this change in the topology of the Fermi surface is also called the 2 1 2 -order phase transition (FT2 1 2 ) according to the Ehrenfest terminology. 6,7 The anomalies arising in the LTT manifest themselves at T ( 3 F not only in the thermodynamic characteristics of metals but also in other characteristics (e.g., magnetic-eld dependence of the electrical resistance, 8 sound absorption [9][10][11][12][13][14][15] ). It has been shown in ref. 16 that the thermoelectric power has the square root divergency at the LTT. The discovery of graphene 17 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][19][20][21][22][23][24][25][26][27][28][29][30] The results of these studies indicate that the effects of the LTT are appreciable and can be used to observe the 2 1 2order 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 lowenergy 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 in 31 theoretically and experimentally.
The 21st century has brought many new results related to graphene thermodynamics. [32][33][34][35][36][37][38][39][40][41][42][43][44][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 rst 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 ow 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 2 Origin of the Lifshitz topological transitions and connection of the twodimensional nature of graphene with the van Hove singularities It is known [see ref. 1, 2, and 3] that at the LTT, a new cavity of the electronic isoenergetic surface appears (or disappears) at the critical energy 3 c in the critical point of the momentum space p ¼ p c , where the electron energy as a function of quasimomentum 3 ¼ 3(p) has a minimum 3 min or maximum 3 max . In this case, the isoenergetic surfaces in the vicinity of p ¼ p c are well described by the ellipsoid equations: where m x , m y , and m z are the main values of the effective mass tensor in the vicinity of 3 min (see ref. 46 45 At the neck rupture, the boundary isoenergetic surface 3(p) ¼ 3 c contains the peculiar point of another type, named the cone point, at p ¼ p c . In this case, the isoenergetic surface in the vicinity of the cone point p ¼ p c is described by the hyperboloid of two sheets at energy 3 < 3 c and the hyperboloid of one sheet at energy 3 > 3 c (Fig. 1a): At energies close to 3 c , one can express the electron density of states (DOS) as where D 0 (3) is the regular smooth function of the energy, and dD(3) depends on the type of LTT. The latter is computed by the relation (per spin direction and per unit volume in 3D-space, or per unit area in 2D-space): where ħ is the reduced Planck-Dirac constant, and D(3) is volume of the ellipsoid (2.1) or (2.2). At the neck rupture ( Fig. 1a), D(3) is the change of the volume enclosed by the plane p 3 ¼ p 0 and the hyperboloid of one (2.3) or two (2.4) sheets ( Fig. 1a) in the disruption of the isoenergetic surface neck in the vicinity of the cone point p ¼ p c .
One can join eqn (2.1)-(2.4) into a single equation: where m i (i ¼ x, y, z) are the positive values of the effective mass tensor, m i s0, and integers a j (j ¼ x, y, z) are equal to AE1. 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 a j .
The point M 0 (min) corresponds to the case where all three coefficients a j ¼ +1. This point corresponds to the minimum in the energy spectrum.
The point M 1 (saddle) corresponds to the case where two coefficients a j are positive and the third one is negative. This is the saddle point (Fig. 1c).
The point M 2 (saddle) is the case where two coefficients a j are negative and a third one is positive. This is the second saddle point (Fig. 1b).
The point M 3 (max) corresponds to the case where all three coefficients a j ¼ À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 a j (j ¼ x, y).
The point P 0 (min) corresponds to the case where both coefficients a j ¼ +1. This point corresponds to the minimum in the energy spectrum.
The point P 1 (saddle) corresponds to the case where one coefficient a j is positive and another one is negative. This is the saddle point. The point P 2 (max) corresponds to the case where both coefficients a j ¼ À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 Q 0 (min) corresponds to the case where coefficient a ¼ +1. This point corresponds to the minimum in the energy spectrum.
The point Q 1 (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 3 c that nature of the Lifshitz topological transition stems from the van Hove singularities (VHS) of the state density dD(3). [48][49][50][51][52] In this connection, the van Hove topological theorem 48 states that the spectrum must contain at least one of the saddle-points S 1 and S 2 (Fig. 1b), and the slope of D(3) must tend to ÀN 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 nite 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 m i s 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 C 6 ¼ ffiffiffiffiffiffiffiffiffiffiffi ffi m s x m s y p =ð2pħÞ 2 because the value of C 6 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 conrmed by recent work. 55 It has been demonstrated in 55 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 AE ln 2/(J À 1/2) with 2 # J # N at low temperatures.

Peculiarities of the Lifshitz topological transitions in graphene
A single layer of graphene consists of carbon atoms in the form of a honeycomb lattice (Fig. 2). The primitive translation vectors e 1 and e 2 form the rhombic unit cell. The hexagonal lattice consists of two trigonal sublattices AAA and BBB. 56 There are four valence electrons (two 2s and two 2p electrons). Three of those participate in the chemical bonding and thus are in bands well below the Fermi energy. We therefore consider the bands formed by the one remaining electron. We assume a tightbinding model in which the electron hops between neighboring atoms. We denote the spacing between neighbouring atoms by a.
From Fig. 2a, we see that two basis vectors of the Bravais lattice are The reciprocal Bravais lattices, b 1 and b 2 , are dened such that The result is This journal is © The Royal Society of Chemistry 2020 RSC Adv., 2020, 10, 27387-27400 | 27389 These basis vectors are of equal length and at 120 ; therefore, the reciprocal lattice is a hexagonal lattice (Fig. 2b). The rst Brillouin zone is shown in Fig. 2b. The rst Brillouin zone is a regular hexagon, whose most characteristic points are its centre G, the inequivalent corners K and K 0 , and the centres of the lateral edges M and M'. The distance GM is The distance GΚ is GM=sin 60 Consider a state with amplitude J R for the electron to be at the site labeled by R. We look for the wave functions with amplitudes which vary, such as e ik$r . There will be different amplitudes J A and J B on sublattices A and B, so J R ¼ J A e ik$r (R˛A), and J R ¼ J B e ik$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), ða ffiffiffi 3 p =2; À a=2Þ; and ðÀa ffiffiffi 3 p =2; À a=2Þ: An atom on sublattice BBB has three neighbors on sublattice AAA at displacements (0, Àa), ð ffiffiffi 3 p a=2; a=2Þ; and ðÀ ffiffiffi 3 p a=2; a=2Þ: In the nearest-neighbour approximation, there are no hopping processes within the sublattices AAA and BBB; hopping only occurs between them. Hence, the eigenvalue 3 and the amplitudes J A and J B are determined for each wavevector k from two equations: Àg 0 e Àikya þ 2e ikya=2 cos where g 0 is the hopping parameter between nearest neighbors.
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 The electron energy eigenvalues of the Hamiltonian (3.7) are given by 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 (k x , k y ) plane for different positions of the wave vector k.
(i) The maximum and minimum energies are 3 c ¼ AE3g 0 , and they occur at k ¼ 0, i.e. in the centre G (Fig. 2b) of the rst Brillouin zone. The value of 3 c ¼ 3g 0 corresponds to the maximum of energy in the conduction band (the point P 2 (max) in Table 1). The value of 3 c ¼ À3g 0 corresponds to the minimum of energy in the valence band (the point P o (min) in Table 1). Taking into account |k y a|, |k x a| ( 1, from (3.9), we have accuracy of order (k y a) 2 , (k x a) 2 : where a is the spacing between neighboring atoms. The graphene reciprocal lattice and the first Brillouin zone (b).
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: Consequently, there exist LTT at critical values of energy 3 c ¼ AE3g 0 , where the cavity (3.11) of the isoenergetic lines disappears (Fig. 3a) at 3 c ¼ +3g 0 in the conduction band or appears at 3 c ¼ À3g 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 is the area of the cavity in the momentum space. We obtain from (3.11) and (3.13) (3.14) where 3 c ¼ 3g 0 . One can represent the change of the density of states dD(3) due to LTT as The discontinuity at the end points 3 ¼ AE3g 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| < 3g 0 : where K(x) is the complete elliptic integral of the rst kind, i.e.
and Re KðxÞ ¼ Expanding of (3.16) near the point G in the small region |3| ( 3g 0 in the vicinity of the LTT points |3 c | ¼ 3g 0 , one obtains (per spin direction and per unit area): The result (3.15) for dD (3) coincides with the total D(3) given by (3.19) in the region |3| ( 3g 0 in the vicinity of the LTT points |3 c | ¼ 3g 0 , where j3j z ffiffiffiffiffiffiffiffiffiffiffiffi 3g 0 j3j p at |3| ( 3g 0 . As follows from (3.19), for graphene in the vicinity of the point G, the coefficients C 5 and C 7 in Table 1 should be changed as Hence, in 2D graphene, the mass-factors of coefficients C 5 and C 7 in Table 1  ; respectively, and this substitution can be treated as the fermion effective mass in the vicinity of the point G: (ii) Now, we compute the dispersion relation in the vicinity of the zone corners K(K 0 ), where the energy tends to zero. We write where K is the wavevector at the zone corner, K ¼ ð4p=ð3 ffiffiffi 3 p aÞ; 0Þ for example, and we will assume that k is small. To the lowest order in k we have, from (3.8) and (3.22), 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 rst Brillouin zone): Consequently, the isoenergetic lines are circles in the vicinity of the minimum of energy 3 c ¼ 0 (for the conduction band) or in the vicinity of the maximum of energy 3 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): The value of 3 c ¼ 0 corresponds to the minimum of energy in the conduction band (the point P 0 (min) in Table 1). The value of 3 c ¼ 0 corresponds to the maximum of energy in the valence band (the point P 2 (max) in Table 1). Eqn (3.25) describes the linear dispersion of the massless Dirac fermions in the vicinity of the corners K(K 0 ) of the rst Brillouin zone: Thus, the LTT exists at the critical value of energy 3 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 0 on the Brillouin zone boundary.
The area of the cavity inside the circle (3.25) is equal to.
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): Computing DOS from (3.28), one has Expanding the expression (3.16) at 3 ¼ 3 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 0 ): Therefore, in the points K(K 0 ), we must make the following substitutions for the coefficients C 6 and C 7 for graphene in Table 1: (iii) Let us compute the dispersion relation in the vicinity of the zone edges, e.g. near the middle M of the boundary of the rst Brillouin zone. We write where K is the wavevector at point M of the zone boundary (Fig. 2b). The distance GM to the center of the edge of the zone is 2p 3a . Therefore: and we will assume that q is small. Taking into account |q y a|, |q x a| ( 1, we have from (3.9), (3.33) and (3.34) to accuracy of order (q y a), (q x a) Therefore, the isoenergetic lines are hyperboles in the vicinity of point M (the middle of the edge of the rst Brillouin zone). One can say that middles of the edges of the rst Brillouin zone are saddle points. We can consider the corresponding points similar to M (the middles of the edges of the rst Brillouin zone) (Fig. 2b) as the saddle points or the "cone points" in two dimensions (points P 1 (saddle) in Table 1). It follows from (3.35) that the isoenergetic lines are described by the following equations in pspace in the vicinity of these saddle points (Fig. 4): The Lifshitz topological transition is realized in the saddle point M by the variation of energy from 3 < 3 c to 3 > 3 c , and 3 c ¼ g 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 D(3) enclosed by hyperbole (3.36) or (3.37) and the corresponding boundary of the rst Brillouin zone p y ¼ ħ/a, or p x ¼ ħ/a, one obtains: where 3 c ¼ g 0 . Now, we are able to resolve the contradiction with the coef-cient C 6 ¼ ffiffiffiffiffiffiffiffiffiffiffi ffi m s x m s y p =ð2pħÞ 2 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 3 c ) is equal to 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 |3| / 3 c ¼ g 0 : Therefore, in the saddle point M, the coefficient C 6 for graphene in Table 1 should be changed to Thereby, in 2D graphene, the mass-factor C 6 in Table 1    From comparison of (4.4) and (1.1), one can see the essential difference of dU given by (4.4) for the LTT in graphene and the LTT in three dimensions (1.1). The second derivative of dU diverges at the point z ¼ 0 as ln(|z|/3 c) in graphene, while only the third derivative of dU 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 G, one computes from (3.14) and (4.2) at T / 0: Considering the appearance of new cavities in the Brillouin zone corners K, we obtain from (3.28) and (4.2), if T / 0: Because the number of electrons in the conduction band is permanent (at least in vicinity of the point m ¼ 3 c ), it is convenient to use the Helmholtz free energy potential F (T, S, N) 2 instead of the grand thermodynamic potential U (T, S, m). 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. 3 c ¼ 3 c (S), and the chemical potential m is also a function of S because of constancy of the carriers: If S c is the area at which the topology of the Fermi lines changes: According to (4.7) and (4.8), the value of |z| ¼ |m À 3 c | can be expressed via |S À S c |: where h is independent of S, and h ¼ h(m). The Helmholtz free energy potential F can be written in the form where F 0 is the smooth part of the free energy. One can readily be convinced that dF is quantitatively equal to the irregular contribution dU expressed in the variables S and T according to the Landau theorem 5 about the small variations of the thermodynamic potentials due to small changes of some parameters of the solid state: (dU) T,S,m ¼ (dF) T,S,N (4.11) One can see that variations of dU 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, dF is given by eqn (4.3), (4.5), and (4.7), where one must set |z| ¼ h|S À S c |.
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: When the electronic cavity disappears in the Brillouin zone centre, one obtains from (4.5) and (4.9) at T / 0: Considering the appearance of new cavities in the Brillouin zone corners, one obtains from (4.6) and (4.9), if T / 0:

The specic heat of the graphene monolayer at the Lifshitz topological transition
The specic heat of graphene is stored in the lattice vibrations (phonons) and the free conduction electrons of graphene, C ¼ C ph + C e . However, phonons dominate the specic heat of graphene 58 at all practical temperatures (>1 K), 59,60,and61,62 and the phonon specic heat increases with temperature. 62 The electron specic heat C e of monolayer graphene (MG) is given by the following relations: 62 where S is the sample area and k B is the Boltzmann constant. We can explore the electron-specic 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: When the electronic cavity disappears in the Brillouin zone centre G, one computes from (4.13) at T / 0: Considering the appearance of new cavities in the Brillouin zone corners K, we obtain from (4.14), if T / 0: Thus, the formulae (4.16), (4.17) and (4.18) describe peculiarities of the specic heat in graphene at the LTT at the points M, G, and K of the rst Brillouin zone, correspondingly (Fig. 5a).
The results for the specic heat (4.17) and (4.18) coincide with result (4.15) of the electron specic heat C MG e of the monolayer graphene (MG) obtained in paper. 61 Indeed, substituting a ¼ ffiffiffi 3 p ð1=2pħV F Þ 2 , and expressing T in Kelvins, one obtains from (4.17) and (4.18): Note that experimental observation of singularities (4.16) and (4.17) is possible under special conditions. The experimental value of the chemical potential m of graphene does not exceed 1 eV, 63 and the estimated value of the overlap integral g 0 is between 2 and 4 eV. [64][65][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| ¼ |m À g 0 | / 0), we must increase the chemical potential m or decrease the overlap integral g 0 due to strain. Increasing m 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 m at low temperature T ( m, especially for high carrier concentrations of n 0 > 10 11 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 t 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 tted values were compared with an available parameterization for the dependence of the orbital overlap on the distance separating the two carbon atoms.

The compressibility of the graphene monolayer at the Lifshitz topological transition
The compressibility of graphene is an important object of study in condensed matter physics because it gives information about the intrinsic nature of the electron structure of graphene and its interactions with external elds. 71,72 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: À2ah g 0 2 3 c ln jzj 3 c ðIIÞ (4.20) where z ¼ h g 0 |S À S c |. When the electronic cavity disappears in the Brillouin zone centre, one obtains from (4.13) for the electron compressibility at T / 0: Considering the appearance of new cavities in the Brillouin zone corners, one obtains from (4.14) for the electron compressibility, if T / 0: where |z| ¼ h 0 |S À S c |, h 0 ¼ h(m / 0). 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, G, and K of the rst Brillouin zone, correspondingly (Fig. 5b). In case of the hyperbolic changes of the Fermi lines at the saddle point A, we obtain the change of the electron thermal coefficient of pressure from (4.12), if T / 0: 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: 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: 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, G, and K of the rst Brillouin zone, correspondingly (Fig. 5c).

The electron thermal expansion coefficient of the graphene monolayer at the Lifshitz topological transition
The thermal properties of graphene have been investigated in recent years; in particular, its thermal expansion and heat conduction have been studied by various theoretical and experimental techniques in ref. 73 and 74. Some theoretical studies that have been carried out to study the thermodynamic properties of graphene (e.g., specic heat and thermal expansion) are based on density-functional theory (DFT) calculations combined with a quantum quasi-harmonic approximation (QHA) for the vibrational modes. 73,74 This is expected to yield reliable results at low temperature for the graphene lattice contribution to the thermal properties; however, it may be questioned with respect to the graphene electron part. The answer to the latter can give the Lifshitz topological transitions.
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 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: 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, G, and K of the rst Brillouin zone, correspondingly (Fig. 5d).

Conclusion and outlook
To summarize the results of the present paper, one can emphasize the following.
(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 P 0 (min) corresponds to the minimum in the energy spectrum. The point P 1 (saddle) corresponds to the case where the Lifshitz topological transitions are realized by variation of the energy from 3 < 3 c to 3 >3 c . The point P 2 (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 G, 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 G at the critical energy value 3 c ¼ 3g 0 , where the cavity of the isoenergetic lines disappears. The existence of the LTT is shown at the critical energy value of 3 c ¼ 0, where six pockets of low energy excitations appear, one for each of the two inequivalent Dirac points K and K 0 .
The Lifshitz topological transition is realized in the saddle point M by variation of the energy from 3 < 3 c to 3 > 3 c and at the critical value of energy 3 c ¼ g 0 . It is shown that the Dirac fermion slows down in the vicinities of the points G and M, becoming not massless but massive, and the values of the fermion effective mass in the vicinity of these points are m G eff z 0.1m e , and m M eff z 0.1m e (where m e is the free electron mass). 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 m 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 m does not exceed 1 eV. 63 Therefore, it is essential to tune m through the VHS by chemical doping 76,77 or by gating. 17,57,[78][79][80][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 m through the VHS by the following two methods. The rst method is connected with deformation of the graphene monolayer to mimic twisted graphene. Rotation between two stacked graphene monolayers in twisted graphene 83 can generate van Hove singularities, which can be brought arbitrarily close to the chemical potential m 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.

Data accessibility
This work includes theoretical investigations.

Authors' contributions
The paper has been written without any assistance, and I gave nal approval for publication.

Funding statement
The work has no source of funding.

Ethics statement
This work did not involve any ethics statements.

Conflicts of interest
I have no competing interests.