1. Introduction
Layered quasi-two-dimensional (Q2D) metals (see
Figure 1 for illustration) represent a wide class of materials, intermediate between the various two-dimensional (2D) electron systems and usual three-dimensional (3D) compounds. The Q2D materials attract enormous research interest due to the variety of new electronic phenomena and diverse potential applications. They include high-temperature superconductors, organic metals, van der Waals crystals, artificial heterostructures, etc. An experimental study of their electronic structure, in addition to ab initio calculations, is crucial for the understanding and utilizing their properties. Although the angle-resolved photoemission spectroscopy (ARPES) provides visual data on electronic bands and Fermi surface (FS), more traditional tools, such as magnetic quantum oscillations (MQO) [
1], usually have higher accuracy and availability.
The MQO originate from the Landau-level (LL) quantization of electron spectrum in a magnetic field
and appear in all thermodynamic and transport quantities as their periodic dependence on
. The standard MQO theory, based on the Lifshitz-Kosevich (LK) formula [
1], is often applied even in Q2D metals and gives the FS extremal cross-section areas
from the fundamental MQO frequencies
. Such data, collected at various directions of the magnetic field
, helped to determine the FS of most known metals. Less accurate, the LK theory also describes the damping of MQO amplitude by temperature and disorder as a function of magnetic field strength
B, which allows the estimate of electron effective mass and mean free time for each band [
1].
The strong anisotropy of layered Q2D metals introduces special features to MQO. The Q2D electron dispersion, obtained in the tight-binding approximation for the interlayer
z direction, is given by
where
is the intralayer electron dispersion,
is the interlayer transfer integral of electrons, and
d is the lattice constant along the
z-axis. Often in Q2D metals the Fermi energy
, and the FS is a warped cylinder with two close cross-section areas corresponding to the ”neck” and ”belly” of the FS (see
Figure 1b). According to the LK theory, at
the MQO have two close fundamental frequencies
and
. As the MQO amplitudes with these two frequencies are also close, the resulting MQO at frequency
F have amplitude modulations with frequency
, called beats. The beat frequency
is proportional to the interlayer transfer integral:
. It non-monotonously depends on the direction of magnetic field, which follows the angular magnetoresistance oscillations (AMRO) [
2,
3,
4,
5] and the can be used to determine the in-plane Fermi momentum [
4,
5] and even high harmonics of FS warping [
6,
7].
Even when the MQO are weak, the interlayer electron transport in Q2D metals exhibits several unique properties, such as the aforementioned AMRO [
2,
3,
4,
5] and longitudinal interlayer magnetoresistance [
8,
9,
10]. There are even more new features of MQO in Q2D metals which are not described within the standard L-K theory. The interplay between the angular and quantum magnetoresistance oscillations is nontrivial [
11,
12] and leads to the angular modulations of MQO amplitude [
12]. The so-called difference or slow oscillations (SO) of magnetoresistance [
13,
14,
15,
16] at frequency
appear and survive at much higher temperature than the usual MQO. The SO help to determine the type of disorder and give other useful information, such as the interlayer transfer integral. Another feature beyond the L-K theory is the phase shift of beats between the MQO of thermodynamic and transport quantities [
14,
17]. This phase shift increases with magnetic field strength
B, as observed in various Q2D metals [
17,
18,
19] and substantiated theoretically [
14,
17]. The experimentally observed phase shift of beats
, usually, considerably exceeds the theoretical prediction
(e.g., see
Figure 3 of Ref. [
14] or Figure 6 of Ref. [
18]). In the present paper we explain this inconsistency and develop a more accurate theory of this effect.
The beatings of MQO amplitude due to the interlayer electron hopping were observed in many other layered materials [
4,
5,
20,
21,
22,
23,
24]. Note that the phase shift of beats is observed not only between the SdH and dHvA effect, but also between the oscillations of other quantities, e.g., the interlayer magnetoresistance and thermopower [
21], or magnetoresistance and specific heat [
25]. The MQO of magnetization and of interlayer conductivity were also measured simultaneously in the Dirac semimetals BaGa
2 [
26]. There are many other combined SdH and dHvA measurements in layered metals where the MQO beats are observed [
27,
28,
29,
30]. Similar beating effects may also appear in artificial heterostructures and multilayer topological insulators [
31,
32], as well as in layered Weyl semimetals [
25].
Another interesting feature of the 3D - 2D crossover, driven by the parameter
, is the phase inversion of MOQ of conductivity as compared to those of the electronic density of states (DoS) at the Fermi level and, hence, of all thermodynamic quantities. At low magnetic field the MQO of conductivity and of DoS have opposite phases, because the former is proportional to the electron mean free time which is inversely proportional to the DoS [
33,
34], according to the golden Fermi rule or the Born approximation. However, in a high magnetic field
in Q2D metals, when the Landau levels (LLs) become separated and the DoS between the LLs is zero, the MQO phase of conductivity and DoS coincide. The latter follows both from the direct observations [
9,
35,
36] and calculations [
8,
9,
10,
16,
37], being supported by the simple qualitative argument that if the DoS at the Fermi level is zero the conductivity must also be zero at low temperature. However, in spite of a general understanding, the quantitative description of this phase inversion is absent. Below we show how this phase inversion during the 3D - 2D crossover is described by analytical formulas. The understanding of this MQO phase inversion is helpful for experimental studies of the Berry phase using MQO measurements [
38,
39,
40,
41].
2. Available Experimental Observations and their Description
The MQO have been extensively investigated in various Q2D metals, including the cuprate [
42,
43,
44,
45] and iron-based [
46,
47] high-Tc superconductors, organic metals [
4,
5], van-der-Walls layered crystals, artificial heterostructures, and many other materials. The two most common MQO experiments are the measurements of magnetization and conductivity oscillations, i.e., the de Haas - van Alphen (dHvA) and the Shubnikov - de Haas (SdH) effect. However, the observation of the phase shift between their beats is less common because it requires the measurements of both these effects on the same sample.
In Ref. [
18] the Fermi surface of a layered organic metal (BEDT-TTF)
4[Ni(dto)
2] was studied by measuring the quantum oscillations of both magnetization and interlayer conductivity. The authors have shown that the standard Lifshitz-Kosevich theory (LK) [
1] well describes the results on magnetization oscillations. The detected positions of the beat nodes of the de Haas - van Alphen oscillation amplitude are well fitted by the formula
where
is the beat frequency giving the FS warping. However the observed position of the beat nodes of the Shubnikov oscillation amplitude in Ref. [
18] are strongly shifted and well fitted by another formula
A similar difference in the positions of beat nodes of de Haas van Alphen oscillations and Shubnikov oscillations was also observed in
[
19]. The LK theory does not explain this clearly observed phase shift of beat nodes about
between the oscillations of magnetization and interlayer conductivity.
The 3D → 2D crossover in MQO behavior happens [
9] when
, i.e., when the LL separation
becomes comparable to the interlayer bandwidth
, where
is the cyclotron frequency,
e – the elementary charge,
c – the speed of light,
– the effective electron mass,
ℏ – the Planck constant. The phase shift between the MQO beats of dHvA and SdH effects is considerable when
, i.e., during this 3D → 2D crossover.
The expressions for interlayer conductivity derived [
13,
17] in the tau approximation strongly underestimate the observed phase shift
between the beats of magnetization and interlayer conductivity, giving
This formula always predicts
, which does not explain the experimental data from [
18]. Note that the phase shift of beats observed in Ref. [
17] is also
times larger than the prediction of this formula. Indeed, substituting the known parameters in the quasi-two-dimensional organic metal
-(BEDT-TTF)
2IBr
2, namely, the cyclotron mass
, as obtained from the temperature dependence of MQO amplitude, and the interlayer bandwidth
meV, as extracted from the ratio between the beating and fundamental MQO frequencies, to Equation (
4) one obtains the
times smaller slope of the linear dependence of
on
than the experimental data in
Figure 4 of Ref. [
17].
The calculations of interlayer conductivity using the Kubo formula predict a larger phase shift of beats, given by [
14]
where
is the so-called Dingle temperature,
is the electron mean free time, and
is the dimensionless parameter. These quantities enter the Dingle factor
, describing the damping of MQO by crystal disorder (see Equations (
6)-(
8) below). However, even this enhanced phase shift of beats, given by Equation (
5), is still insufficient to describe quantitatively the effect observed in Refs. [
17,
18].
Note that the electron mean free time
measured from the MQO is, usually, shorter than the transport mean free time
and than the mean-free time
arising from short-range disorder only and measured from the damping of slow magnetoresistance oscillations [
13]. This difference,
, appears because the MQO are damped not only by short-range disorder but also by long-range sample inhomogeneities which smear the Fermi energy similar to the temperature effect [
9,
13,
14,
15,
16]. The opposite case when the transport mean free time
is also possible in heterogeneous conductors, where there are rare but strong inhomogeneities as, e.g., domain walls or linear crystal defects. Then inside each domain of the size larger than the Larmor radius the material is a clean metal with a large
, but the electronic transport across the sample is difficult because of its heterogeneity, which corresponds to a short
. As follows from the calculations in Refs. [
14,
16], it is
rather than
enters the Equation (
5).
4. Analysis and Results
Let us first consider the limit of
, corresponding to a rather large electron interlayer transfer integral
or weak magnetic field. Then one can use the large-argument asymptotic expansions of the Bessel functions in Equations (
7) and (
8). As a result, the beating factors of the MQO amplitudes in Equations (
7) and (
8) simplify to
At
the second term in Equation () is negligibly small, and we get the same beating factors of the MQO of magnetization and interlayer conductivity,
Hence, in this limit the phase shift of beats is absent, .
In the opposite limit of strong MQO damping by disorder,
, the first term in Equation () can be omitted, and one obtains
Comparing Equations (
9) and (
12) we see that at
the MQO beats of interlayer conductivity are shifted from those of magnetization by the phase
. As follows from Equation (
9), the beat nodes of magnetization oscillations
in Equation (
8) are located at
, and from Equation (
12) it follows that the beat nodes of the interlayer conductivity oscillations
given by Equation (
7) are shifted by
and located at
, where
is an integer number. These
give the same position of the beat nodes as the expressions (
2) and (
3) for
and
, corresponding to the fit of experimental data in Ref. [
18]. Hence, Equations (
9) and (
12) describe well the experimental observations in Ref. [
18]. However, Equation (
12) assumes
, while
in the experiment in Ref. [
18]. The Dingle temperature extracted from the experimental data in Ref. [
18] is
K. The beat frequency in the organic metal (BEDT-TTF)
4[Ni(dto)
2] studied in Ref. [
18] is
T, which corresponds to the interlayer transfer integral
meV
K. Hence, the ratio
in Ref. [
18]. Hence, to improve the quantitative description of the observed phase shift of MQO beats we need to consider the case
more accurately.
Figure 2.
The beating factors and , describing the MQO amplitudes of interlayer conductivity and magnetization (brown solid curve), as a function of the inverse magnetic field . Here is the beat frequency, and the ratios correspond to the blue thick long-dashed, red short-dashed, green long-dashed and black dotted curves for .
Figure 2.
The beating factors and , describing the MQO amplitudes of interlayer conductivity and magnetization (brown solid curve), as a function of the inverse magnetic field . Here is the beat frequency, and the ratios correspond to the blue thick long-dashed, red short-dashed, green long-dashed and black dotted curves for .
The phase shift
of MQO beats becomes even stronger than that given by Equations (
9) and (
12) if one applies Equations (
7) and (
8) without the asymptotic expansions of the Bessel functions at large argument, i.e., without the assumption
, which is more relevant to the experiment in Ref. [
18]. In
Figure 2, using Equations (
7) and (
8) at
, we plot the beating factors
and
, describing the MQO amplitudes of magnetization and conductivity without the Dingle and temperature damping factors. The theoretical curves in
Figure 2 correspond to four values
to describe the experimental data from Refs. [
17] and [
18] where
and
correspondingly. In
Figure 2 we take a special scale of the abscissa axis so that the zeros of both curves, called the beat nodes, fall into integer values. From
Figure 2 we see that the beat nodes for interlayer conductivity oscillations
are shifted relative to the magnetization
beat nodes by a quarter of the period
.
1 The beat nodes of magnetization oscillations
given by Equation (
8) in
Figure 2 are located at
. This coincides with the prediction of simplified limiting-case formula (
9) and agrees perfectly with the experimental data from Ref. [
18], fitted by Equation (
2). At
, the positions of beat nodes of the amplitude of conductivity oscillations are close to the position of the beat nodes of the amplitude of magnetization oscillations. The beat nodes of interlayer conductivity oscillations
given by Equation (
7) and plotted in
Figure 2 at
are located at
. For
they coincide with Equation (
3), which fits the experimental data from Ref. [
18].
Equations (
7) and (
8) also perfectly agree with the experimental data from Ref. [
17]. In
Figure 3 we plot the measured positions of magnetization and conductivity beat nodes from the experiments in Refs. [
17,
18] together with the theoretical predictions from formulas (
7) and (
8). The measured positions of the beat nodes in Ref. [
18] at
are very well described by Equations (
2) and (
3) and, therefore, denoted by ”dHvA fit” and ”SdH fit” in
Figure 3. The points ”dHvA exp” and ”SdH exp” in
Figure 3 are taken from the experimental data in
Figure 2 or 3 of [
17]. The prediction of Equations (
7) and (
8) are denoted by ”SdH thr” and ”dHvA thr” in
Figure 3. From
Figure 3 we see a perfect agreement between the above theory and the experimental data from Refs. [
17,
18] for the beat-node positions with number
. These beat-node positions for the MQO of magnetization and interlayer conductivity are fitted well by two straight lines parallel to each other. Hence, we conclude that Equations (
8) and (
7) describe very well the beating nodes with
.
Figure 3.
The dependence of the beat-node number
of the MQO of
and
on the inverse magnetic field
at
. Points <<dHvA fit>> and <<SdH fit>> refer to the fit expressions (
2) and (
3) from Ref. [
18]. Points <<dHvA exp>> and <<SdH exp>> refer to the experimental data from
Figure 2 and 3 of Ref. [
17]. Points <<dHvA thr>> are found from the zeros of expression (
9), while points <<SdH thr>> are found numerically from the solution of equation
.
Figure 3.
The dependence of the beat-node number
of the MQO of
and
on the inverse magnetic field
at
. Points <<dHvA fit>> and <<SdH fit>> refer to the fit expressions (
2) and (
3) from Ref. [
18]. Points <<dHvA exp>> and <<SdH exp>> refer to the experimental data from
Figure 2 and 3 of Ref. [
17]. Points <<dHvA thr>> are found from the zeros of expression (
9), while points <<SdH thr>> are found numerically from the solution of equation
.
The position
of the zeroth beating node given by Equation (
7) holds at any ratio
, as illustrated in
Figure 2 and can be shown analytically. Indeed, in a very high magnetic field
, the ratio
, and the beating nodes are given by the equation
is the root of this equation at any ratio
. The second root
of Equation (
13) is physically irrelevant as it gives a too large value
much beyond the limit
where the expansion (
13) holds.
The position of the very first beat node at
differs in Equations (
3) and (
7). According to Equation (
3), the beating node
of
oscillations should be at
, which does not coincide with the position
or
of the zeroth beating node given by Equation (
7) and shown in
Figure 2. As we just have shown, the origin of this difference is not a ”wrong” choice of the parameter
, entering Equation (
7).
In fact, the position
of the zeroth beating node
, predicted by the ”fit” Equation (
3), has not been observed experimentally. The beat frequency in the organic metal (BEDT-TTF)
4[Ni(dto)
2] studied in Ref. [
18] is
T. According to Equation (
3) the zeroth beating node is at
T. In Ref. [
18] the experimental data on the SdH oscillations in (BEDT-TTF)
4[Ni(dto)
2] are shown up to
T only in Figure 6 of Ref. [
18], while the data up to
T are given only for dHvA oscillations in
Figure 3 of Ref. [
18]. In Figure 6.1 of Ref. [
50] the SdH data up to
T are shown for the same organic metal (BEDT-TTF)
4[Ni(dto)
2], and no beating nodes are observed in the magnetic-field interval
T. A similar difference in the beat-node positions of de Haas - van Alphen and Shubnikov oscillations in
[
19] was also detected by the non-observation of the SdH beat node in the magnetic-field interval
T. Hence, Equation (
7) describes very well the available experimental data on the MQO beats of interlayer conductivity.
In Ref. [
17] the Dingle temperature extracted from the experimental data is
K, but the Dingle temperature from the scattering by short-range disorder is
K (see
Figure 4 of Ref. [
13]). The beat frequency
T in the organic metal
-(BEDT-TTF)
2IBr
2 studied in Ref. [
17] is much larger than in (BEDT-TTF)
4[Ni(dto)
2] studied in Ref. [
18]. It corresponds to the interlayer transfer integral
meV
K and
. For
the expression for
in Equation () is approximately equal to that in Equation (
11). Since formulas (
9) and (
11) are similar, this means that the beat nodes for magnetization
and conductivity oscillations
in this limit are close to each other. In
Figure 4 we show that the beat nodes found at
fit well with the experimental data on the position of beat nodes from
Figure 3 of Ref. [
17]. The higher is the number
of the MQO beat nodes of magnetization and interlayer conductivity, the shorter is the distance between them along the abscissa axis in the units of
.
Figure 4.
The dependence of the number
of the beat node of magnetic quantum oscillations of
and
on the inverse magnetic field
at
. The points <<dHvA exp>> and <<SdH exp>> refer to the experimental points taken from the
Figure 3 of [
17]. The points <<dHvA thr>> are found from the zeros of expression (
9), while the points <<SdH thr>> are found numerically from the solution of equation
. Here
.
Figure 4.
The dependence of the number
of the beat node of magnetic quantum oscillations of
and
on the inverse magnetic field
at
. The points <<dHvA exp>> and <<SdH exp>> refer to the experimental points taken from the
Figure 3 of [
17]. The points <<dHvA thr>> are found from the zeros of expression (
9), while the points <<SdH thr>> are found numerically from the solution of equation
. Here
.
5. Discussion
The above analytical formulas describe the interlayer conductivity
. They can also be applied to describe the first MQO harmonic of interlayer magnetoresistance
, if the MQO are not very strong so that the non-linear effects in MQO amplitude are not important [
51]. If only the first MQO harmonic is kept and the higher-order terms in the MQO amplitude are negligible, one can safely apply the relation
to describe the interlayer magnetoresistance [
51]. However, for strong MQO their averaging over temperature and long-range disorder differs [
9,
51] for
and
.
Above we have shown that the developed theory and Equation (
7) describe very well the available experimental data on the MQO beats of interlayer conductivity
. They also predict some new interesting and observable features. In particular, Equation (
7) predicts a strong decrease of the MQO beating factor
of interlayer conductivity with the increase of magnetic field, corresponding to
. This decrease of
is partially compensated by the Dingle factor
and by the temperature damping factor
. However, the product
, describing the MQO amplitude of interlayer conductivity, may still be non-monotonic in high field, in contrast to the amplitude of magnetization MQO described by Equation (
8), which always monotonically increases in a very high field at
. This nonmonotonic magnetic-field dependence of the MQO amplitude of interlayer conductivity in a high field is rather robust to the variation of material parameters, e.g, to variations of the ratio
. In
Figure 5 we show this non-monotonic dependence of the
MQO amplitude
on
for two different ratios
and temperatures
T. At the small ratio
the increase of temperature up to the Dingle temperature
has almost no effect on the graph. The predicted decrease of MQO amplitude of the interlayer conductivity in a high field is somewhat counterintuitive, being opposite to what is observed in magnetization, but it can be easily tested experimentally in various Q2D layered compounds, such as organic metals, van der Waals crystals, artificial heterostructures, etc. We suppose that the proposed decrease of MQO amplitude in a high field was actually observed in many experiments but mistaken for the non-existing beat node at a high magnetic field beyond the available range. For example, as we mentioned in the previous section, in Refs. [
18,
19,
50] the beating node
of the SdH effect was not observed but assumed from some decrease of MQO amplitude in a very high field. This decrease may be considered as the experimental indication of the predicted high-field non-monotonic behavior of the MQO amplitude of interlayer conductivity, but further experimental studies of this interesting effect are required to compare with our theoretical predictions and to determine the range of parameters where it is observed.
Another amazing prediction of our theoretical analysis is that all the curves in
Figure 2, describing the MQO amplitudes, periodically cross at the same points. It can also be checked experimentally. However, this experimental test is less obvious and convenient than the above prediction of the monotonic magnetic-field dependence of
MQO amplitude, because the corresponding experimental data for
and
magnetic oscillations must be properly normalized for this comparison.
Now we consider the MQO phase of interlayer conductivity. For the Q2D electron spectrum (
1) in a magnetic field one can easily calculate [
48,
49] the one-particle DoS
. If we neglect higher harmonics, the DoS oscillations are given by [
48,
49]
where
is its nonoscillating part per one spin component. The phase of DoS oscillations is strictly tied to that of magnetization given by Equation (
8), because the beating factors
. In a weak magnetic field, when
corresponding to the 3D limit, the second term in the beating factor
of interlayer conductivity given by Equation
7 is much smaller than its first term. Then we have
, and the MQO phases of the interlayer conductivity
and of the DoS
are strictly opposite, as predicted by the 3D theory of SdH effect [
33,
34]. In the opposite case
, corresponding to the 3D →2D crossover or nearly 2D limit, the SdH beating factor
is negative, as follows from Equation
7 and illustrated in
Figure 5. Since
is also negative at
, at the 3D →2D crossover
or in the almost 2D limit
the quantum oscillations of the electronic DoS and of the interlayer conductivity have the same phase. We see that our theory describes the phase inversion of the SdH oscillations during the 3D →2D crossover.
The above analysis is generic to layered Q2D metals and can be applied to a very large class of layered compounds, which are actively studied by MQO till nowadays [
24,
26,
27,
28,
29,
52,
53,
54,
55].