## Numerical Study of Quantum Oscillations of the Quasiparticle Lifetime: Impurity Spectroscopy, Novel Electric Field and Strain Effects

Valentin Leeb 1, 2 and Johannes Knolle 1, 2, 3

1 Technical University of Munich, TUM School of Natural Sciences, Physics Department, TQM, 85748 Garching, Germany 2 Munich Center for Quantum Science and Technology (MCQST), Schellingstr. 4, 80799 M¨nchen, Germany u 3 Blackett Laboratory, Imperial College London, London SW7 2AZ, United Kingdom (Dated: August 2, 2024)

Quantum oscillation (QOs) measurements constitute one of the most powerful methods for determining the Fermi surface (FS) of metals, exploiting the famous Onsager relation between the FS area and the QO frequency. The recent observation of non-Onsager QOs with a frequency set by the difference of two FS orbits in a bulk three-dimensional metal can be understood as the QO of the quasiparticle lifetime (QPL) due to interorbital scattering [Huber, Leeb, et al. , Nature 621 (2023)]. QPL oscillations (QPLOs) generalize magneto-intersubband oscillations (MISOs) known from coupled two-dimensional metals. They may provide a novel tool for extracting otherwise hardto-measure intra- versus interband scattering times of quasiparticles. Here, we provide a numerical lattice study of QPLOs comparing transport and thermodynamic observables. We explore the effect of different imperfections like general impurities, Hall effect-induced electric fields, various forms of strain from bending, and magnetic field inhomogeneities. We confirm the basic phenomenology of QPLOs as predicted in analytical calculations and identify additional novel, non-perturbative features. Remarkably, we find that some imperfections can stabilize, or even enhance, non-Onsager QPLOs in contrast to standard QO frequencies. We discuss various avenues for identifying QPLOs in experiments and how to use their dependence on imperfections to extract material properties.

## I. INTRODUCTION

Quantum oscillation (QO) measurements are an exceptionally sensitive tool for measuring Fermi surface (FS) geometries as well as the strength of interaction effects via extracting the effective masses m i from the temperature dependence [1-3]. Often, the observation of a QO frequency serves as a unique and sensitive indicator for the existence of a FS. For example, in underdoped cuprates, QO studies famously confirmed the presence of a closed FS pocket in a magnetic field [4, 5]. Similarly, QO experiments in the spin density wave parent phase of iron-based superconducting compounds showed the emergence of small pockets, in contrast with angleresolved photoemission spectroscopy (ARPES) results at the time [6-8]. Agreement on the electronic structure suggested by QOs was confirmed only much later [9].

argued that QPLOs are generic for bulk metals with multiple extremal orbits or multiple FS pockets. This claim is supported by the identification of many further candidate materials [12], ranging from topological semi-metals to iron-based superconducting parent compounds [14]. A key observation of Ref. [10, 14] was that QPLOs could be erroneously interpreted as FS orbits, leading to incorrect identification of the electronic structure of a given material.

The interpretation of QOs, as measured in transport or thermodynamic observables, is based on the famous Onsager relation, which ascribes each QO frequency to a semi-classical FS orbit [2, 3]. Recently, unusual frequencies have been reported in the multi-fold semi-metal CoSi [10], which are forbidden in the semi-classical Onsager theory and generalizations thereof [11]. In addition, they are at odds with known mechanisms generating non-Onsager QO frequencies like magnetic breakdown, chemical potential oscillations, magnetic interaction, and (Stark) quantum interference [10, 12]. The emergence of these forbidden frequencies as combinations of two semiclassical QO frequencies, originating from underlying FSs, has been explained by a non-linear coupling of the FSs via interband impurity scattering [13], dubbed quasiparticle lifetime oscillations (QPLOs) here. Ref. [10]

The basic ingredient leading to QPLOs is a coupling between two QO orbits, such as interband impurity scattering. The inverse lifetime of the quasiparticles on a given orbit then is a sum of three terms -a constant, the so-called Dingle temperature T D 1 , an oscillatory contribution with frequency F 1 from intraorbit scattering, and, crucially, an oscillatory term with frequency F 2 determined by the other orbit 2 with an amplitude proportional to interorbit scattering. The constant Dingle temperature term damps the oscillations with a universal factor, known as the Dingle factor R D = exp( -2 π mT 2 D /eB ), with e the electron charge and B the magnetic field. The third term leads to cross terms of QOs in the conductivity with frequency F 1 and F 2 which appear as the sum and difference of the basis frequencies F 1 and F 2 in the conductivity

$$\mathfrak { a } _ { \mathfrak { c } _ { \mathfrak { h o } } } ^ { \mathfrak { a } _ { \mathfrak { a } } } \quad \sigma _ { \pm } \infty \cos \left ( 2 \pi \frac { F _ { 1 } \pm F _ { 2 } } { B } \right ) R _ { D 1 } R _ { D 2 } R _ { T } \left ( m _ { 1 } \pm m _ { 2 } \right ), \ \ ( 1 )$$

where R T ( m ) is the Lifshitz-Kosevich dependence as introduced below in Eq. (17).

Conceptually, QPLOs generalize so-called magnetointersubband oscillations (MISOs) known in coupled 2D electron gases [15] as well as quasi-2D materials [16]

and observed in several experiments [17-23]. MISOs are nowadays claimed to be understood in great detail [24] which is supported by several analytic works [24-30]. The new insight of QPLOs is the generalization to bulk 3D materials [13] and a new qualitative picture of their origin [10]. Hereinafter, we refer to the phenomenon of additional QO frequencies, which originate in oscillatory contributions of the self-energy due to interorbit scattering, as QPLOs.

So far, all theories of QPLOs are based on analytical, perturbative calculations, requiring approximations like the self-consistent Born approximation for deriving Eq. (1) [13] or even more ad-hoc assumptions. This calls for a controlled check of the analytic predictions, which numerical simulations can provide. The numerical evaluation of QOs in the presence of disorder is numerically not challenging, in the sense that the Hilbert space dimension grows only linearly with system size, which allows for an efficient numerical simulation at least for 2D systems, e.g., for transport [14, 31]. Besides supporting the phenomena itself, numerical simulations offer the possibility to test predictions of analytic theories in a more controlled manner than experiments can, provding additional insights to advance our understanding of QPLOs. In addition, one can capture effects beyond the limit of perturbation theory.

In this work, we aim to benchmark analytical theories of QPLOs by full numerical lattice simulations of the Shubnikov-de Haas (SdH) effect, i.e., the conductance as function of the magnetic field. In addition, we show how QPLOs respond to different types of impurities compared to standard Onsager QO frequencies, which paves the way for using QPLO measurements to understand otherwise hard to obtain impurity properties of materials. Beyond disorder, we also study the effects of i) electric fields from current-induced charge accumulation, ii) strain engineering from bending nanostructures, and iii) magnetic field inhomogeneities. The unique response of QPLOs to such imperfections may, on the one hand, serve as an additional criterion simplifying their identification and, on the other hand, allow to extract material properties.

The work is structured as follows. In sec. II we introduce the lattice model and the numerical method. In the main part, sec. III, we discuss the response of QPLOs to different types of impurities. Next, in sec. IV, we analyze the properties of QPLOs like the appearance of higher harmonics, the dependence of QPLO on the nonlinear coupling of the FSs, the phase relation to the basis frequencies and the behavior in thermodynamic observables. In the following sections, we discuss the effect of finite temperature, sec. V; electric fields induced by a Hall effect from a varying current strength, sec. VI); crystal strain induced by bending, sec. VII; and inhomogeneities in the magnetic field, sec. VIII.

FIG. 1. Reduction scheme to describe the extremal QO orbits of layered quasi-2D materials with a basic 2D two-band model. (a) Quasi-2D materials consist of layers which are strongly coupled within each layer by t, t ′ whereas the transfer integral in z -direction t ⊥ is smaller. (c) The FSs are typically warped cylinders where the extremal QO orbits for B ∝ ˆ z are located at k z = 0 , π . (d) An analogous description in 2D is just a simple two-band model where bands are split by 2 t ⊥ . (b) The corresponding real space picture is just a system with two layers where the bands correspond to the bonding and antibonding states of the layers.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00115v1-figure-6d589f15-cea1-4b6f-9665-27e7132c93c9.jpg)

## II. TWO-BAND MODEL AND SET-UP

QPLOs appear in two and three-dimensions [13]. Here, we focus on a basic 2D two-band model that also captures the relevant physics of quasi-2D materials. The latter are characterized by weakly coupled layers and a nonzero interlayer tunneling leads to a warping of the FS such that 2 extremal orbits at k z = 0 and k z = π appear. To include these orbits determining the QOs, it is sufficient to model only 2 layers, which would for periodic boundary conditions in z -direction (which are trivial in that case) correspond to two discrete momenta.

Hence, we consider the Hamiltonian

$$H _ { \text{quasi-2D} } & = - \, t \sum _ { \langle i j \rangle, z } c _ { i, z } ^ { \dagger } c _ { j, z } - t ^ { \prime } \sum _ { \langle \langle i j \rangle \rangle, z } c _ { i, z } ^ { \dagger } c _ { j, z } \\ & \quad - \frac { t _ { \perp } } { 2 } \sum _ { \langle z, z ^ { \prime } \rangle, i } c _ { i, z } ^ { \dagger } c _ { i, z ^ { \prime } }$$

on a square lattice, where t ( t ′ ) are intralayer (next-) nearest-hoppings and t ⊥ ≪ t is the interlayer hopping, see Fig. 1 (a) and (c). For simplicity we focus on z = 0 1 , which is essentially the 2D limit, shown in Fig. 1 (b) and (d), and introduce the orbital (layer) basis Ψ i =

B

B

B

B

B

B

B

B

B

B

B

B

B

B

B

B

B

B

B

B

B

B

B

B

B

B

B

B

B

B

B

B

B

B

B

B

L

L

(b)

140

120

100

80

60

40

20

0

12

10

8

6

4

2

0

0.00

Conductance for independent impurities with

=0.4

t

fit with 5th order polynomial

50

100

150

200

250

0/

B

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00115v1-figure-0c9728fd-4d71-4980-afe8-2047ebedbf9a.jpg)

0.25

0.50

0.75

1.00

1.25

1.50

1.75

2.00

QO freq.

F

[

]

FIG. 2. (a) shows a schematic sketch of the implemented model used for conductance simulations. The scattering region (blue vertices) of system size L × L implements the tight-binding model with impurities H = H 0 + H imp in the presence of a magnetic field B . The leads (red vertices) are translational invariant, hence do not have impurities, and are not exposed to the magnetic field. The conductance G from the left lead to the right lead is computed. (b) shows the exemplary result for the conductance G for independent impurities. We fit the conductance in the field region I B = (20 250) (marked by the black dotted lines) , with a 5 th -order polynomial, which is then subtracted from the numerical data to obtain the oscillating part of the signal, see (c) for a part of the oscillating signal. To obtain a spectrum, (d) the oscillating signal is zero-padded, scaled with a windowinf function, and then Fourier transformed. The absolute of this, in general, complex-valued data is the fast-Fourier transformation (FFT) amplitude. The identified peaks are consistent for various window types.

( c i, 0 , c i, 1 ) such that

## A. Peierls substitution

$$H _ { 0 } = - \, t \sum _ { \langle i j \rangle } \Psi _ { i } ^ { \dagger } \Psi _ { j } - t ^ { \prime } \sum _ { \langle \langle i j \rangle \rangle } \Psi _ { i } ^ { \dagger } \Psi _ { j } - t _ { \perp } \sum _ { j } \Psi _ { j } ^ { \dagger } \tau ^ { x } \Psi _ { j } \quad \text{We} \\ ( 3 ) \quad \text{hypi}$$

$$= - \sum _ { k } \Psi _ { k } ^ { \dagger } [ 2 t ( \cos k _ { x } + \cos k _ { y } ) \\ + 4 t ^ { \prime } \cos k _ { x } \cos k _ { y } + t _ { \perp } \tau ^ { x } ] \Psi _ { k }$$

$$+ \, 4 t ^ { \prime } \cos k _ { x } \cos k _ { y } + t _ { \perp } \tau ^ { x } ] \Psi _ { k }$$

$$( 4 )$$

$$= - \sum E _ { \pm } ( k ) \Phi _ { k } ^ { \dagger } \Phi _ { k }$$

$$= - \sum _ { k } E _ { \pm } ( k ) \Phi _ { k } ^ { \dagger } \Phi _ { k } & & ( 5 ) & & \text{Her} \\ & & \text{gau} _ { \dots }$$

where we introduce the bands E ± ( k ) = -2 t cos k x -2 t cos k y -4 t ′ cos k x cos k y ± t ⊥ and the fields in the band basis Ψ k = U Φ k which are just bonding and antibonding states

$$U = \frac { 1 } { \sqrt { 2 } } \begin{pmatrix} 1 & 1 \\ 1 & - 1 \end{pmatrix}.$$

Weapply an orbital magnetic field B in the out of plane z -direction, by Peierls substitution, which transforms the hopping t r r , ′ from site r to site r ′

$$t _ { r, r ^ { \prime } } \rightarrow t _ { r, r ^ { \prime } } \exp \left ( - \text{i} e \int _ { r } ^ { r ^ { \prime } } A ( l ) \text{d} l \right ). \quad \ \ ( 7 )$$

Here, A = ( -yB, 0 0) is the vector potential. , For this gauge choice, only the hoppings t x in x-direction and the diagonal next-nearest neighbor hoppings t xy acquire a phase.

Additionally, we consider potential impurities

$$H _ { i m p } = \sum _ { j } \Psi _ { j } ^ { \dagger } \Lambda _ { j } ^ { \Psi } \Psi _ { j }$$

where Λ i is, in general, a random hermitian matrix. In rectangular

Hamming

Blackman-Harris

]

2

[

G e

(d)

FFT amplitude

sec. III A we specify different types of the scattering matrix Λ Ψ j and study their effect on QOs and QPLOs.

## B. Numerical method

Numerically, it is more convenient to work with dimensionless quantities. We set the lattice constant and ℏ to 1. The dimensionless magnetic field B/ Φ 0 can be interpreted as the amount of flux through a single plaquette measured in units of the flux quantum Φ 0 = 2 π e .

Our goal is a full numerical lattice calculation of the SdH effect. This requires computing the conductance G B ( ) of the model given by H = H 0 + H imp as function of the magnetic field B and averaged over impurity samples.

We evaluate the conductance through the built-in Landauer-B¨ uttiker algorithm [32-35] of the python software package kwant [36]. For this, we construct an open lattice, dubbed scattering region, of L × L sites with 2 orbitals per site and implement H including the Peierls phases. Additionally, we attach two leads of width L on the right and left side of the scattering region, see Fig. 2 (a). The leads realize H 0 without magnetic field and are characterized by translational invariance in x -direction. The Landauer-B¨ uttiker algorithm evaluates the transmission probability T nm ( E ) from an eigenstate n with energy E of the left lead to an eigenstate m with energy E of the right lead via an S -matrix approach. The conductance is the sum over all channels at the Fermi energy µ [32-35]

$$G = \frac { e ^ { 2 } } { 2 \pi } \sum _ { n, m } T _ { n m } ( \mu ). \quad \quad ( 9 ) \quad \text{cai} \\ \text{tur} \\ \text{an}.$$

We evaluated the conductance G numerically for 2500 equidistant values of Φ 0 / B from 10 to 250 and observed SdH oscillations as a function of Φ 0 / B , see Fig. 2 (b). We then analyzed the Fourier transform in Φ 0 / B with standard QO techniques, which include subtraction of a 5th-order polynomial background to extract the oscillatory signal, shown in Fig. 2 (c), zero padding to increase the point density of the spectrum and windowing with a Blackman-Harris window to decrease spectral leakage. The representative steps are shown in Fig. 2. Fourier spectra, see Fig. 2 (d), are shown with frequencies in units of the BZ area. They show sharp peaks at the contributing QO frequencies.

We also implemented a numerical simulation of the de Haas-van Alphen (dHvA) effect. For this, we use the scattering region without additional leads. We then compute the density of states (DOS) ρ E ( ) via the kernel polynomial method [37], which is based on a Chebyshev expansion of the spectral density. For sampling, we used 30 randomly chosen vectors such that only the bulk of the system is sampled. Cutting off the edges suppresses finite-size boundary effects. We defined the bulk by the set of all lattice points at least 40 sites away from the edges and used 7000 Chebyshev moments.

The DOS can be used to determine thermodynamic observables. First, the zero temperature thermodynamic potential is computed

$$\Omega ( \mu ) = \int _ { - \infty } ^ { \mu } \mathrm d E \rho ( E ) ( E - \mu ). \quad \quad ( 1 0 )$$

The magnetization M can then be computed via the derivative M = -∂ Ω /∂B . Here, we experienced that a numerical evaluation of the derivative of the data obtained with the kernel polynomial method is unstable. A different thermodynamic observable, the particle number N , is also defined via a derivative N = -∂ Ω( µ /∂µ ) , however it can be carried out analytically to obtain

$$N ( \mu ) = \int _ { - \infty } ^ { \mu } \mathrm d E \rho ( E ). \text{ \quad \ \ } ( 1 1 )$$

The particle number can hence be evaluated directly from ρ without increasing the noise on the data by an additional numerical derivative.

The oscillating part of particle number and magnetization are related by N ∝ B m F M for a single frequency F with effective mass m . Hence, the main difference is that the QO amplitude of N reduces for decreasing magnetic field B . This does not pose a problem in numerical simulations since the magnetic fields are relatively large.

## III. IMPURITY SPECTROSCOPY WITH QPLOs

Within our basic model, different types of impurities can be classified according to their orbital/layer structure. We first discuss various types of impurity vertices and then show how they lead to qualitatively different effects.

## A. Impurity models

## 1. Identical impurities

In the simplest impurity model, the vertex at each point is taken to be proportional to the identity Λ Ψ i = λ τ i 0 and the prefactor is drawn from λ i ∈ -[ Λ 0 / 2 Λ , 0 / 2] randomly and uniformly. The physical picture is that each layer is an identical copy of the other layer, which may arise from systematic defects when growing the material or atoms sitting between the layers, hence having the same effect on both. A transformation to the band basis Λ Φ i = U † Λ Ψ i U shows that the scattering vertex remains diagonal in the band basis, therefore leading only to pure intraband scattering.

## 2. Opposite identical impurities

The impurities at each lattice point are again perfectly correlated but with opposite signs in each layer, Λ Ψ i =

λ τ i z . Their strength is drawn from λ i ∈ -[ Λ z / 2 Λ , z / 2] randomly and uniformly. The corresponding physical picture could be dipoles located between the layers. A transformation to the band basis Λ Φ = λ τ i x shows that it gives rise to a pure interband coupling term.

## 3. Independent impurities

Both impurity models above are somewhat artificial. Therefore, we also consider impurities randomly distributed in each layer, i.e.,

$$\Lambda _ { i } ^ { \Psi } = \begin{pmatrix} \lambda _ { i } & 0 \\ 0 & \eta _ { i } \end{pmatrix} \quad \quad \ \ ( 1 2 )$$

where λ , η i i ∈ -[ Λ 2 Λ 2] are distributed randomly, uni-/ , / formly, and independently. Crucially, the vertex is a superposition Λ i = ( λ i + η i ) / 2 τ 0 +( λ i -η i ) / 2 τ z of the identical and the opposite identical impurities. The probability distribution of ( λ i ± η i ) / 2 is the convolution of the 2 uniform (rectangular) probability distributions and hence a triangular probability distribution. Transforming this to the band basis shows that independent impurities in the layers lead to equally strong intra- and interband scattering, as expected from the form of Λ Ψ i .

## 4. Random interlayer coupling

Finally, there are interlayer impurities Λ Ψ i = λ τ i x or Λ Ψ i = x τ i y where λ i ∈ [ -Λ 2 Λ 2] randomly and uni-/ , / formly. Unlike the other disorder models, random interlayer coupling constitutes a form of bond disorder and can hence be complex. A varying distance between the layers will only influence the real part of the hoppings. A complex disordered part of the hopping may, for example, arise from locally trapped fluxes or random spin-orbit coupling (i.e., heavy impurities).

Real interlayer impurities, i.e. Λ Ψ i = λ τ i x , transform into the band basis as Λ Φ ∝ τ z , hence, only leading to intraband contributions. In contrast, complex interlayer impurities, i.e., Λ Ψ i = λ τ i y , contribute only to interband scattering because the scattering vertex transforms to Λ Φ ∝ τ y .

## B. SdH oscillations

Next, we show numerical results of the QO spectra for different impurity types. QPLOs only emerge for an effective coupling of the semiclassical orbits [13]. In sec. III A, we argue that in our model, only opposite identical and independent impurities lead to an effective interpocket coupling, whereas identical or interlayer impurities only contribute to intrapocket scattering processes. Therefore, we expect a strong dependence of the QO spectra on the types of impurities.

FIG. 3. QO spectra of the different impurity models for impurity strengths Λ ν = 0 4 and for comparison in the absence . of impurities (black). Spectra are normalized and shifted for clarity. The main frequencies F , F 1 2 are marked by gray, solid lines, and the combination frequencies F 2 ± F 1 by gray, dashed lines. QPLOs, visible by peaks at F 2 ± F 1 , only appear for opposite identical, independent, and complex interlayer impurities. Additional peaks are higher harmonics. ( t ⊥ = 0 4 , . t µ = -0 5 , . t I B = (20 250), Blackman-Harris window) ,

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00115v1-figure-274ada30-6681-4fcc-94b8-14b702c5bca5.jpg)

We computed the SdH effect numerically for the different impurity models, following sec. II B. Fig. 3 demonstrates that the main frequencies F , F 1 2 consistently appear for any model whereas the combination frequencies F 2 ± F 1 only appear for sufficiently strong interpocket coupling, i.e., for opposite identical, independent and complex interlayer impurities. Indeed, the presence of a difference frequency peak from QPLOs may serve as tool discriminating between different forms of impurities.

In the remainder of this work, we will focus on opposite identical and identical impurities, where the strength of interband scattering is Λ z and the strength of intraband scattering Λ 0 .

## IV. ANALYSIS OF QPLOs

## A. Higher harmonics

QOs are not perfectly harmonic, i.e., their shape is not a perfect cosine. In the spectrum, this becomes visible

FIG. 4. (a) FFT amplitude for various fillings µ . (b) We extracted the position of all peaks above a detection limit defined by the 1.25 times the height of the low-frequency peak typically appearing in the data. The line indicate the analytic predictions for the basis frequencies and their higher harmonics (gray), and the combination frequencies (red) for different orders. ( t ⊥ = 0 4 , . t Λ z = 0 4 , . t I B = (50 250), Blackman-, Harris window)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00115v1-figure-d68744d8-eed4-4eab-bac7-6552de6473e5.jpg)

by additional peaks at integer multiples kF, k ∈ N of the basis frequency F . The k th harmonic always appears together with k th -order of the Dingle damping factor R k D and is, therefore, typically decreasing in magnitude.

Research on QPLOs has so far focused on 2nd order frequencies, i.e., F 2 ± F 1 which are of second order in the Dingle factor, namely governed by R D 1 R D 2 . Higher harmonics also exist for QPLOs as pointed out in Ref. [13]. For a model with two orbits, like our model here, higher harmonics of order | k 1 | + | k 2 | can be found at k F 1 1 + k F , k 2 2 i ∈ Z .

We analyzed the spectra for various different Fermi energies µ . Fig. 4 (a) shows that a complicated QO spectrum with several peaks arises. We extract the peak positions of the QO spectrum and compare the results to analytic predictions, which we obtain from calculating the area of the semiclassical orbits and their integer and sum combinations. In panel (b), we demonstrate that nearly all peaks may be explained by higher harmonics of QPLOs up to 5th order following the analytic predictions [13].

We note that for isolated parameters, additional peaks appear in the QO spectrum, which higher harmonics cannot explain. However, these do not seem to exist consistently for different Fermi energies. These unexplained peaks could, for example, be artifacts of the Fourier transformation, e.g., side lobes. We do not have an explanation for the peaks at low frequencies and high Fermi energies. Finally, we observe a weak splitting of some combination frequencies towards low Fermi energies, which we discuss in sec. IV D.

## B. dHvA oscillations

The appearance of QOs in thermodynamic quantities, historically typically the magnetization/susceptibility, are known as the dHvA effect. dHvA oscillations are in many aspects similar to the ones appearing in the SdH effect and for standard Onsager QOs most aspects carry over.

However, it turns out that QPLO are substantially different in transport compared to thermodynamic quantities. In analytic calculations for parabolic bands, all difference frequency combinations cancel exactly in the dHvA effect, whereas sum combinations do not [13]. This result also holds qualitatively true for relativistic (linear) dispersions, i.e., the amplitudes of difference frequency QOs are strongly suppressed. From the analytic theory, it is, therefore, unclear whether the absence of difference frequency QPLOs in the dHvA effect is generic.

In Fig. 5, we compare quantities from the SdH effect (conductance) and from the dHvA effect (particle number and DOS) for increasing interband scattering rate. Fig. 5 confirms that difference frequency combinations of QPLOs are strongly suppressed in thermodynamic quantities for generic band structures, if not entirely absent.

Our high-quality spectrum of the DOS, Fig.5 (c), demonstrates that QPLOs cannot be understood on the basis of the DOS alone. Additional knowledge about the lifetime, hence the full spectral function, is required to explain them.

## C. Amplitude relation

Next, we study the dependence of the QO amplitudes on the interband scattering Λ z . In order to rule out the effect of a varying Dingle temperature, i.e., effects induced by the changing broadening of the LLs, we change the intra- Λ 0 and interband scattering Λ z such that their square sum (Λ ) 0 2 + (Λ ) , which is proportional to the z 2 Dingle temperature [13], remains constant.

Even though the Dingle temperature is kept constant, we find an increase of the main frequencies by roughly a factor of 8 when comparing pure intraband coupling to pure interband coupling. Hence, we analyze the relative amplitudes of the higher harmonics with respect to the algebraic mean of the peak height of the two basis frequencies, see Fig. 6. The second harmonics 2 F , 1 2 F 2 remain roughly constant as expected. The combination frequencies originating from QPLO vanish for decreasing interband coupling Λ . z We fitted the amplitudes of the frequencies F 2 ± F 1 with (Λ ) z α from 0 to 0.6 and found α ≈ 1. This result is at odds with analytic predictions. QPLOs arise in second order in the Born approximation and should, therefore, show a scaling behavior of α = 2. The stronger scaling points towards the fact that QPLOs are either a lower order phenomenon, which does not hold in analytic calculations, or that higher order scattering contributes more strongly than expected. The stronger

FIG. 5. Comparison of the spectra of transport/ conductance, i.e., SdH effect (a) and thermodynamic quantities, i.e., dHvA effect, the particle number N (b) and the DOS ρ (c) as a function of increasing interband scattering rate Λ /t . The main frequencies F , F 1 2 and their higher harmonics kF , kF 1 2 are marked by gray, solid lines. Sum (Difference) frequency combinations are marked by red (gray) lines, dashed lines for second order, and dash-dotted lines for third order. Difference frequency combinations are only present in the conductance, i.e., in the SdH effect. ( t ⊥ = 0 4 , . t µ = -0 5 , . t I B = (20 250), , Blackman-Harris window)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00115v1-figure-e4b74a6e-e5e2-4473-af6f-b50f766ac007.jpg)

scaling indicates that the scattering expansion is nonperturbative, rendering the second-order approximation only qualitatively correct.

A second important observation from Fig. 6 is the nonmonotonic behavior of difference frequency combinations. The amplitudes of F 2 -F 1 (2 F 1 -F 2 ) go to zero at around Λ z ≈ 1 5Λ . 0 (Λ z ≈ 2Λ ). We have checked that the zeros 0 do not depend on the absolute strength of Λ 2 z +Λ . It is 2 0 unclear what the origin of these zeros is.

## D. Peak splitting

A detailed analysis of the major peaks of the spectrum shows that, in some cases, they are actually split into two subpeaks, shown in detail in Fig. A.1 (a). The peak splitting is also visible in Fig. 4 (b) F 2 -F 1 at low µ , in Fig. 5 (a) F 2 -F 1 for Λ /t = 0 2, in Fig. 9 (a) . F 1 , F 2 at eR H I = 0 06 . t and F 1 + F 2 at eR H I = 0 03 . . t

We show in the appendix sec. B that peak splitting is

FIG. 6. Relative amplitudes of main frequencies for different inter- vs. intraband scattering ratios. Λ z = 0 4 . t sin φ , Λ 0 = 0 4 . t cos φ such that the Dingle temperature ∝ (Λ ) 0 2 +(Λ ) z 2 is constant. The solid lines indicate fits with ∝ (Λ ) z α , we find α F 2 + F 1 = 1 4 . ± 0 1 and . α F 2 -F 1 = 1 0 . ± 0 2. . The gray shaded area indicates where amplitudes would be below the detection limit.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00115v1-figure-20400e86-ccc3-49ec-a920-7eb5e1e1fbf6.jpg)

independent of the impurity type. The physical origin is that due to the oscillating quasiparticle lifetime, the QO amplitudes are polynomial functions of the magnetic field [13]. At the zeros of the amplitudes, the phase of the QOs jumps. Therefore, the phase itself becomes a function of the magnetic field. In the spectrum, a phase jump shows up as a weak peak splitting. However, peak splitting is of minor experimental relevance, because the zeros are located at values of T D /B which are mostly not tractable in experiment, see appendix sec. B.

## E. Phase relation

The phase of QOs is often of interest because of its relation to the Berry phase [38, 39]. However, an absolute determination of the phase is usually difficult in experiment because large fields close to the quantum limit, i.e., B ∼ F , of a QO frequency F are required. In contrast, the relative phase between different frequencies can already be estimated by comparing the raw oscillatory data.

Ref. [10, 13] established, in contradiction with Ref. [28], that there is a fixed phase relation between the phase ϕ ± of the QPLOs and the phases ϕ i of the basis frequencies

$$\phi _ { \pm } = \phi _ { 2 } \pm \phi _ { 1 } \text{ \ or \ } \phi _ { \pm } = \phi _ { 2 } \pm \phi _ { 1 } + \pi. \text{ \ \ (13)}$$

The additional term π accounts for the fact that the amplitudes may change sign, as also discussed in the appendix sec. B.

We extracted the phases ϕ , ϕ i ± as function of the magnetic field by filtering, see Fig. 7 (d,e). We also compared the relative phase of the, beating pattern, which is typical for close frequencies, with the filtered sum and difference

FIG. 7. (a) Oscillations of the sum of the filtered basis frequencies (blue) show a beating pattern. The filtered oscillations of the difference frequency (red) envelope the beating pattern, indicating a fixed phase relation between the phase of the difference frequency ϕ -and the difference of the phases of the main frequencies ϕ 2 -ϕ 1 . There is also a fixed phase relation for filtered oscillations of the sum frequency (orange) and the basis frequencies. At high (low) fields (b) ((c)), the minima (maxima) of the sum align with extrema of the beating pattern. (d) The filtered phase of the difference frequency (red) can be described by the difference of the filtered phases of the basis frequencies (blue). (e) The filtered phase of the sum frequency (orange) can be described by the sum of the filtered phases of the basis frequencies (blue). (Λ z = 0 4 , . t µ = -2 ) t

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00115v1-figure-5dc8164e-50e0-42ba-849b-b27ba157573d.jpg)

frequencies, see Fig. 7 (a-c). For a description of the filtering method, see appendix sec. A. We observe constant phases over large magnetic field ranges. ϕ -experiences an abrupt phase change of π at ϕ /B 0 = 100. At the same magnetic field, ϕ + shows a ≈ 0 16 . π phase change.

Our analysis shows that the phase of QPLOs is indeed the sum/difference of the basis frequencies as described by Eq. (13). Our numerical results indicate that the even stronger relation

$$\phi _ { - } & = \phi _ { 2 } - \phi _ { 1 } \quad \text{and} \quad \phi _ { + } = \phi _ { 2 } + \phi _ { 1 } + \pi \quad \text{or} \\ \phi _ { + } & = \phi _ { 2 } + \phi _ { 1 } \quad \text{and} \quad \phi _ { - } = \phi _ { 2 } - \phi _ { 1 } + \pi \quad \quad ( 1 4 ) \quad \text{For} \ \, \varepsilon$$

might hold, i.e., that the amplitudes of difference and sum frequency have always reversed signs. To our knowledge, this has not been discussed in the literature before. We conclude that a relative phase analysis, as done in Fig. 7, constitutes a strong indicator to confirm that a frequency in the spectrum can be assigned to QPLOs and not as a standard Onsager frequency.

## V. FINITE TEMPERATURE

The finite temperature conductance can be obtained within the Landauer-B¨ttiker formalism from u

$$G _ { T } = \frac { 2 \pi } { e ^ { 2 } } \int _ { - \infty } ^ { \infty } \mathrm d E ( - n _ { F } ^ { \prime } ( E - \mu ) ) \sum _ { n, m } T _ { n m } ( E ), \quad ( 1 5 )$$

a convolution with the derivative of the Fermi distribution function [40]

$$- n _ { F } ^ { \prime } ( \epsilon ) = \frac { e ^ { \epsilon / T } } { T ( 1 + e ^ { \epsilon / T } ) ^ { 2 } }. \quad \quad ( 1 6 )$$

For a numerical evaluation of this integral, we use that n ′ F ( E -µ ) is strongly peaked around µ and limit the integration boundaries to ≈ µ ± 6 T to capture more than 99.5% of the spectral weight.

We evaluated the zero temperature conductance G E ( ) for several energies E sampled logarithmically around the Fermi energy µ = -t . Applying Eq. (15) we evaluated G T at µ = -t for several temperatures and then extracted the temperature dependence of the main peaks of the spectrum of G T , see Fig. 8.

FIG. 8. QO amplitudes of several different frequencies for different temperatures. The amplitudes are extracted from the peak maxima of the spectra at different temperatures. The solid lines show the expected LK dependence with R T ( m ) (Eq. (17)) where m are combinations of m 1 = 0 47 . / Φ 0 and m 2 = 0 57 Φ . / 0 which are calculated from the band structure and for Φ 0 /B we used the mean window point. ( t ⊥ = 0 4 . t, µ = -t, Λ z = 0 4 , . t I B = (100 250), Blackman-Harris , window)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00115v1-figure-b620c22f-bde3-4fb9-b74f-3aa00650fa1d.jpg)

Fig. 8 shows a perfect Lifshitz-Kosevich dependence with

$$R _ { T } ( m ) = \frac { \pi m T \Phi _ { 0 } / B } { \sinh \left ( \pi m T \Phi _ { 0 } / B \right ) } \quad \quad ( 1 7 ) \quad \begin{pmatrix} \text{FIG} \\ \text{ind} \\ \text{for} \ \end{pmatrix}$$

of all QO frequencies. The effective mass m = d e F E / E ( ) d can be computed from the band structure by taking the derivative of the FS area with respect to the Fermi energy, obtaining m 1 = 0 47 Φ . / 0 and m 2 = 0 57 Φ . . / 0 The theory of QPLO predicts for combination frequencies k F 1 1 + k F 2 2 a LifshitzKosevich temperature dependence with an effective mass k m 1 1 + k m 2 2 . The numerical data in Fig. 8 confirm this result with high precision. Most importantly, it confirms the presence of the nearly temperature-independent difference frequency for m 1 ≈ m 2 .

We note that small deviations between the numerical results and the expected theory curve appear for increasing temperature. The deviations increase with decreasing mean of I B , i.e, with increasing magnetic field. We suggest that this is due to our method of extracting the temperature dependence from the Fourier transform, which is to some extent inaccurate because of the dependence of R T ( m ) on the magnetic field B . We verified that higher order corrections to the Lifshitz-Kosevich factor are much smaller and are therefore not responsible for discrepancies between numerics and analytics. 1

1 Lifshitz-Kosevich dependence assumes that the frequency can be linearly expanded around the chemical potential F µ ( + ) = ϵ F µ ( ) + ϵF ′ ( µ ) where F ′ ( µ ) determines the effective mass and that higher order corrections are sufficiently small.

FIG. 9. QO spectrum in the presence of finite electric field induced by a classical Hall effect. (a) Exemplary QO spectra for various currents eR H I . The spectra are offset for clarity. The main frequencies F , F 1 2 are marked by gray, solid lines and the combination frequencies F 2 ± F 1 by gray, dashed lines. (b) Peak maxima of various frequencies as function of the current/Hall coefficient eR H I . ( t ⊥ = Λ z = 0 4 . t, µ = -0 5 , . t I B = (80 250), Blackman-Harris window) ,

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00115v1-figure-40316663-ca70-40fc-8ce9-f7133020eaf0.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00115v1-figure-08e15465-cd99-4c46-9ce6-8b34162d81d7.jpg)

## VI. ELECTRIC FIELD EFFECTS

Next, we explore the effect of different types of additional perturbations on the QO spectrum.

In a typical magneto-transport set-up to measure the SdH effect, a current I is sent through the sample to measure the resistivity. Naturally, the set-up is susceptible to a classical Hall effect, where charge accumulates on the edges parallel to the current directions. The strength of the Hall effect can be characterized by a sample-dependent Hall coefficient R H . The generated perpendicular Hall voltage U H = R IB/ H Φ 0 leads to an electric field affecting the QO spectrum.

We model the effect of a Hall voltage-induced electric field in our system by a spatially varying Fermi energy µ → µ y ( ) = µ + eU H y/L = eR H IBy/L Φ . 0 The leads are kept at constant Fermi energy.

A qualitative picture of the effect of a spatially varying Fermi energy can be obtained by a phase smearing argument [3]. Note that formally, the application of a phase smearing argument requires that variations appear only

on a length scale longer than the size of a cyclotron orbit. Here, µ y ( ) varies continuously within the cyclotron orbit. According to phase smearing, the damping of the QO frequencies can be evaluated by averaging over all present Fermi energies

$$\Gamma \text{conv} \, \text{sounds} \, \text{sounds} \\ \frac { 1 } { L } \int _ { - L / 2 } ^ { L / 2 } \text{d} y \cos \left ( 2 \pi \frac { F ( \mu + e U _ { H } y / L ) } { B } \right ) ^ { ( \text{c} ) } \\ \approx & \int _ { - 1 / 2 } ^ { 1 / 2 } \text{d} u \cos \left ( 2 \pi \frac { F + m U _ { H } u } { B } \right ) ^ { ( \text{c} ) } \\ = & \cos \left ( 2 \pi \frac { F } { B } \right ) \text{sinc} \left ( \pi \frac { m R _ { H } I } { \Phi _ { 0 } } \right ) ^ { ( \text{sinc} ) } \\ \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \quad \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \ \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \, \text{sblock} \                                                                                                                                                                                                         \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock}" \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock}\quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \\ \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \ \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \, \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \	 \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \                                                                                                                                                                                                 \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \opacity \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock } \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock} \quad \text{sblock})$$

where sinc x = sin x/x .

The numerical results are shown in Fig. 9. We observe a stronger decay of the main frequencies with respect to eR H I compared to the analytical prediction in Eq. Eq. (18). We attribute this to the fact that the cyclotron orbits are not sufficiently large compared to the variation of the Fermi energy. Hence, the assumptions for the derivation of Eq. (18) are not met. Nevertheless, the derivation provides an intuitive picture of the qualitative effect. With increasing electric field (or Hall coefficient/current), the QO amplitude should monotonically decrease to zero with a scale which is determined by the effective mass of the QO frequency. If this semiclassical phase smearing picture can be applied similarly to finite temperature, it implies that the amplitude of the difference frequency decreases weakly with m 2 -m 1 .

The main finding of this section is the anomalous behavior of the difference frequency combinations F 2 -F 1 and 2 F 1 -F 2 in Fig. 9. Both amplitudes show a significant non-monotonic increase of more than a factor 6 before dropping again. This behavior is at odds with the present understanding of QPLOs or QOs in general.

## VII. STRAIN EFFECTS FROM BENDING

Intense material research in recent years has led to new developments, which allow strain tuning of materials in a controlled manner [41]. Specifically, there exist new experimental set-ups for continuously bending materials [42, 43].

Bending may have different effects depending on the microscopic details of the set-up and material. However, two effects appear quite generically and are straightforward to incorporate in our numerical implementation: (i) Bending leads effectively to a spatially dependent magnetic field by deforming the area of every single plaquette of the crystal differently, which changes the flux piercing through a plaquette. This effect becomes dominant when bending around an axis perpendicular to the magnetic field, e.g., the y -axis, see Fig 10 (a) and (c). (ii) The distance between the atoms in the crystal and, hence, the overlap of the atomic orbitals changes. This leads to spatially dependent hoppings, see Fig. 10 (b) and (d).

FIG. 10. Different forms of bending: (a) Bending around the y -axis, here κ y = 0 7, leads to a decrease of the flux piercing . through a single plaquette, see panel (c). (b) Bending around the z -axis, here κ z = 0 6 varies the longitudinal distance . a between the atomic orbitals. We assume t ∝ 1 /a and obtain a spatial dependence of the longitudinal hoppings ( t x ) and the next-nearest neighbor hoppings ( t xy ), see panel (d).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00115v1-figure-c93fcedc-0d27-46f7-93b7-8b96d1f57622.jpg)

This effect becomes dominant when bending in the plane perpendicular to the magnetic field.

Note that bending inside a material cannot only arise from external application of forces. Crystal imperfections may lead to similar phenomena, which become apparent in domains of slightly different material parameters and grain boundaries.

Here, we consider the effects of bending types (i), (ii) separately in order to understand their effects on QOs in a controlled manner. For simplicity, we focus on lattices with uniform curvature. We quantify the bending strength by the curvature κ j around the axisj in the middle of the lattice. The leads are kept at zero curvature.

## A. Bending-induced effective magnetic field

A non-zero value of κ y leads to a spatially dependent magnetic field B x ( ) = B (0) cos( κ x/L y x ) through each plaquette, where B (0) = B is the flux per plaquette of the relaxed model, see Fig. 10 (a) and (c). Note that the spatial dependence of the magnetic field is weak on the scale of the lattice, hence the vector potential A inherits the spatial dependence of the magnetic field.

The naive expectation is that the electrons only feel the average magnetic field acting on the system ¯ B = ∫ 1 0 B x ( )d( x/L ) = B sinc κ y as it has already been discussed in the standard literature [3]. Better estimates can be made by employing a phase smearing argument as Eq. (18). However, even for κ y ≪ 1 only small analytic progress can be made favoring a qualitative picture. The varying magnetic field damps the oscillations. Higher frequencies are damped stronger than lower frequencies.

Our numerical observations contradict this naive ex-

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00115v1-figure-27adc04e-f092-46ab-a96c-522658100f0a.jpg)

(b)

10 1

10 0

0.0

0.1

0.2

0.3

0.4

2

2

F

F

F

F

1

2

1

2

F

0.5

2

0.6

y

FIG. 11. QO spectrum in the presence of finite bending in y -direction, i.e., spatially dependent magnetic field. (a) Exemplary QO spectra for various bendings κ y . The spectra are offset for clarity. The main frequencies F , F 1 2 are marked by gray, solid lines and the combination frequencies F 2 ± F 1 by gray, dashed lines. (b) Peak maxima of various frequencies as function of the curvature κ y in y -direction. ( t ⊥ = Λ z = 0 4 . t, µ = -0 5 , . t I B = (80 250), Blackman-Harris , window)

pectations, see Fig. 11. At first, the amplitudes of semiclassical frequencies and sum combinations of QPLO remain relatively stable with increasing curvature κ y before eventually decreasing. Most remarkably is the behavior of the amplitudes of the difference frequency combinations F 2 -F 1 and 2 F 1 -F 2 . The amplitudes first decrease until κ y = 0 1 but then increase by roughly a factor of 3. .

## B. Spatially varying hoppings

Next, we study spatially dependent hoppings, see Fig. 10 (b) and (d). Similar to above, we expect that bending leads to a monotonic decrease of all QO amplitudes.

We observe a strong, almost monotonic decrease of the QO amplitude of F 1 , whereas the amplitude of F 2 remains stable with the tendency to even increase (factor ≈ 1 3), see Fig. 12. . Similar to above, the difference frequency increases in amplitude by a maximum factor of

F

F

1

1

F

Ai

+

FIG. 12. QO spectrum in the presence of finite bending in z -direction, i.e., spatially varying hoppings. (a) Exemplary QOspectra for various bendings κ z . The spectra are offset for clarity. The main frequencies F , F 1 2 are marked by gray, solid lines and the combination frequencies F 2 ± F 1 by gray, dashed lines. (b) Peak maxima of various frequencies as function of the curvature κ z in z -direction. ( t ⊥ = Λ z = 0 4 . t, µ = -0 5 , . t I B = (80 250), Blackman-Harris window) ,

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00115v1-figure-0d325c23-f5a2-4720-9cd7-ab781d7e15a4.jpg)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00115v1-figure-f65434ba-6d71-4d1e-bbf5-d8bf3d6b192f.jpg)

roughly 2, peaking around κ z ≈ 0 35. .

We have verified that the stability of F 2 is also present without interband coupling. The stability probably arises due to the quantum Hall limit where transport is mediated by edge states and is not generic for QOs.

## VIII. MAGNETIC FIELD INHOMOGENEITIES

Finally, we study the effect of inhomogeneous magnetic fields. Historically, unwanted field inhomogeneity posed a technical challenge in the early days of QOs, as already noted in Landau's seminal work with the original prediction of QOs [44]. We explore the possibility of magnetic fields which are not constant over the sample. A possible reason for this is a spatially dependent magnetic susceptibility χ x, y ( ), which changes the magnetic field B → B (1 + χ x, y ( )). For simplicity we model χ x, y ( ) = χ 0 (sin(2 πx/L x ) + sin(2 πy/L y )).

An increasing magnetic field inhomogeneity should decrease the QO amplitude, as already discussed in the

FIG. 13. QO spectrum in the presence of magnetic field inhomogeneities, i.e., B → B (1 + χ x, y ( )). (a) Exemplary QO spectra for various strengths of the field inhomogeneity χ 0 . The spectra are offset for clarity. The main frequencies F , F 1 2 are marked by gray, solid lines and the combination frequencies F 2 ± F 1 by gray, dashed lines. (b) Peak maxima of various frequencies as function of the magnetic field inhomogeneity χ 0 . ( t ⊥ = Λ z = 0 4 . t, µ = -0 5 , . t I B = (80 250), Blackman-, Harris window)

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00115v1-figure-0a50b54b-797a-42d5-ae9e-88e7b95146f0.jpg)

literature [3]. The limit at which the damping becomes relevant is χ 0 ≈ B/πF . From this condition, it is apparent that higher frequencies are affected more strongly or earlier than lower frequencies.

Fig. 13 shows the QO spectrum and the amplitudes for various χ 0 . As expected, the QO amplitudes decrease monotonically. However, yet again, we find an unusual non-monotonic behavior for the difference frequency F 2 -F 1 . After an initial decay, it peaks at roughly the same value as for zero inhomogeneity before decaying below the detection limit.

## IX. SUMMARY AND DISCUSSION

In this work, we studied SdH and dHvA QOs in the presence of different types of impurities, electric fields, bending-induced strain, and magnetic field inhomogeneities. The motivation for this quantitative study comes from the analytical prediction [13] and observation [10] of a new mechanism of non-Onsager difference frequency QO dubbed QPLOs, which we benchmarked here in a full numerical lattice implementation.

We focused on a minimal but generic model featuring two FSs leading to two QO frequencies F 1 and F 2 . Different types of impurities are introduced and by transforming the impurity vertex into the band basis, we evaluated the strength of intra- and interband coupling between the two semiclassical QO orbits F 1 , F 2 . Weconfirmed numerically that new QPLO frequencies at F 2 ± F 1 appear only in the presence of nonzero interband coupling. We also observed higher harmonics k F 1 1 + k F 2 2 of the QPLO up to 5th order. We confirmed that the temperature dependence of k F 1 1 + k F 2 2 is an exact Lifshitz-Kosevich with an effective k m 1 1 + k m 2 2 . In particular, we observed the temperature stability of F 2 -F 1 , making it the only observable frequency for T &gt; 0 025 , a tell-tale prediction . t for their experimental identification.

In accordance with the theory for parabolic and linear bands [13], all difference frequency contributions F 2 -F 1 , 2 F 2 -F 1 , ... vanish in the dHvA effect. For parabolic bands, the cancellation of difference frequency contributions in analytics appears somewhat fine-tuned. Therefore, the exact absence in numerics points towards a more general reason that transport and thermodynamics quantities behave qualitatively differently for QPLOs.

In sec. IV C, we showed that the amplitudes of QPLOs scale linearly with interband coupling and not as expected quadratically. This shows that the expansion of the self-energy in the self-consistent Born approximation is non-perturbative, rendering the approximation only qualitatively correct. Additionally, we observe that difference frequencies may cancel even for finite interband couplings for particular parameter values.

In the remainder of our work, we showed that the relative amplitude of QPLOs increases, i.e., it has a qualitatively different behavior from the main Onsager frequencies under several different types of imperfections. The imperfections include electric fields induced by currents through the sample, i.e., a finite Hall effect, bending of the sample such that hoppings or the magnetic field acquire a spatial dependence and fluctuations of the magnetic field. In the case of the Hall effect and spatially fluctuating magnetic field, we find that the absolute amplitude of F 2 -F 1 increases substantially. The amplitude dependencies are at odds with existing analytical theories, i.e., they cannot be explained by simple phase smearing arguments, as for finite temperature, or by a combination of the amplitudes of the basis frequencies. We verified that none of the above imperfections induces combination frequencies on its own. The unique response of the frequency F 2 -F 1 to these imperfections demonstrates the non-perturbative nature of this QO frequency arising from intriguing interference effects not captured by semi-classical arguments.

In the absence of interband scattering, we never observe combination frequencies in our exact numerics, even for non-trivial imperfections. Thus, interband scattering is key for inducing these new non-Onsager frequencies.

In turn, an observation thereof can be used to perform 'impurity spectroscopy,' e.g., extracting otherwise hardto-measure properties of intra- versus interband/orbital impurity contributions [14].

Our work has established rigorously the key aspects of difference frequency QPLOs and unearthed a whole range of new phenomena beyond the perturbative theory of Ref. [13] pointing to effects beyond the lowest order Born approximation of interband scattering, which will be very worthwhile to explore in the future. In addition, interband scattering may also arise naturally from interactions and collective bosonic excitations like phonons or (para-)magnons. While similar QPLOs are expected in these cases neither analytical nor numerical works exist and are an important, yet challenging, avenue for future research. Again, the hope is that QPLOs may serve as a novel tool for extracting otherwise hard-to-obtain material properties.

Data and code availability.Code and data related to this paper are available on Zenodo [45] from the authors upon reasonable request.

## ACKNOWLEDGMENTS

We acknowledge helpful discussions and related collaborations with N. Huber, P. Bieniek, M. A. Wilde, and C. Pfleiderer. We would like to thank P Moll for pointing us to the idea of studying QO in strained materials. V. L. acknowledges support from the Studienstiftung des deutschen Volkes. We acknowledge support from the Imperial-TUM flagship partnership. JK acknowledges support from the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany's Excellence Strategy-EXC-2111-390814868, DFG grants No. KN1254/1-2, KN1254/2- 1, and TRR 360 - 492547816; as well as the Munich Quantum Valley, which is supported by the Bavarian state government with funds from the Hightech Agenda Bayern Plus.

- [1] W. J. de Haas and P. M. van Alphen, The dependence of the susceptibility of diamagnetic metals upon the field, Proceedings of the Academy of Science of Amsterdam 33 , 1106 (1930).
- [2] L. Onsager, Interpretation of the de haas-van alphen effect, The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science 43 , 1006 (1952).
- [3] D. Shoenberg, Magnetic Oscillations in Metals (Cambridge University Press, 1984).
- [4] N. Doiron-Leyraud, C. Proust, D. LeBoeuf, J. Levallois, J.-B. Bonnemaison, R. Liang, D. Bonn, W. Hardy, and L. Taillefer, Quantum oscillations and the fermi surface in an underdoped high-t c superconductor, Nature 447 , 565 (2007).
- [5] S. E. Sebastian, N. Harrison, and G. G. Lonzarich, Towards resolution of the fermi surface in underdoped hightc superconductors, Reports on Progress in Physics 75 , 102501 (2012).
- [6] S. E. Sebastian, J. Gillett, N. Harrison, P. Lau, D. J. Singh, C. Mielke, and G. Lonzarich, Quantum oscillations in the parent magnetic phase of an iron arsenide high temperature superconductor, Journal of Physics: Condensed Matter 20 , 422203 (2008).
- [7] T. Terashima, N. Kurita, M. Tomita, K. Kihou, C.-H. Lee, Y. Tomioka, T. Ito, A. Iyo, H. Eisaki, T. Liang, M. Nakajima, S. Ishida, S.-i. Uchida, H. Harima, and S. Uji, Complete Fermi Surface in BaFe2As2 Observed via Shubnikov-de Haas Oscillation Measurements on Detwinned Single Crystals, Phys. Rev. Lett. 107 , 176402 (2011).
- [8] A. I. Coldea, D. Braithwaite, and A. Carrington, Ironbased superconductors in high magnetic fields, Comptes Rendus Physique 14 , 94 (2013).
- [9] M. D. Watson, P. Dudin, L. C. Rhodes, D. V. Evtushinsky, H. Iwasawa, S. Aswartham, S. Wurmehl, B. B¨chner, u M. Hoesch, and T. K. Kim, Probing the reconstructed Fermi surface of antiferromagnetic BaFe2As2 in one domain, npj Quantum Mater. 4 , 1 (2019).
- [10] N. Huber, V. Leeb, A. Bauer, G. Benka, J. Knolle, C. Pfleiderer, and M. A. Wilde, Quantum oscillations of the quasiparticle lifetime in a metal, Nature 621 , 276 (2023).
- [11] A. Alexandradinata and L. Glazman, Fermiology of topological metals, Annual Review of Condensed Matter Physics 14 , 261 (2023).
- [12] V. Leeb, N. Huber, C. Pfleiderer, J. Knolle, and M. A. Wilde, A field guide to non-onsager quantum oscillations in metals, in preparation (2024).
- [13] V. Leeb and J. Knolle, Theory of difference-frequency quantum oscillations, Phys. Rev. B 108 , 054202 (2023).
- [14] V. Leeb and J. Knolle, Interband scatteringand nematicity-induced quantum oscillation frequency in FeSe, Phys. Rev. B 109 , L081109 (2024).
- [15] V. Polyanovsky, Magnetointersubband oscillations of conductivity in a two-dimensional electronic system, Fiz. Tekh. Poluprovodn. 22 , 1408 (1988).
- [16] V. Polyanovsky, High-temperature quantum oscillations of the magnetoresistance in layered systems, Phys. Rev. B 47 , 1985 (1993).
- [17] D. R. Leadley, R. J. Nicholas, J. J. Harris, and C. T. Foxon, Influence of acoustic phonons on inter-subband
- scattering in GaAs-GaAlAs heterojunctions, Semicond. Sci. Technol. 4 , 885 (1989).
- [18] P. T. Coleridge, Inter-subband scattering in a 2D electron gas, Semicond. Sci. Technol. 5 , 961 (1990).
- [19] D. R. Leadley, R. Fletcher, R. J. Nicholas, F. Tao, C. T. Foxon, and J. J. Harris, Intersubband resonant scattering in GaAs-Ga 1 -x Al As heterojunctions, Phys. Rev. B x 46 , 12439 (1992).
- [20] A. V. Goran, A. A. Bykov, A. I. Toropov, and S. A. Vitkalov, Effect of electron-electron scattering on magnetointersubband resistance oscillations of two-dimensional electrons in GaAs quantum wells, Phys. Rev. B 80 , 193305 (2009).
- [21] W. Mayer, A. Ghazaryan, P. Ghaemi, S. Vitkalov, and A. A. Bykov, Positive quantum magnetoresistance in tilted magnetic field, Phys. Rev. B 94 , 195312 (2016).
- [22] S. Abedi, S. A. Vitkalov, N. N. Mikhailov, and Z. D. Kvon, Anomalous Decay of Quantum Resistance Oscillations of 2D Helical Electrons in Magnetic Field, Sci Rep 10 , 7875 (2020).
- [23] G. M. Minkov, O. E. Rut, A. A. Sherstobitov, S. A. Dvoretski, N. N. Mikhailov, V. A. Solov'ev, M. Yu. Chernov, S. V. Ivanov, and A. V. Germanenko, Magnetointersubband oscillations in two-dimensional systems with an energy spectrum split due to spin-orbit interaction, Phys. Rev. B 101 , 245303 (2020).
- [24] M. E. Raikh, High-temperature magneto-inter-chirality oscillations in 2D systems with strong spin-orbit coupling (2023), arXiv:2310.00774 [cond-mat].
- [25] M. E. Raikh and T. V. Shahbazyan, Magnetointersubband oscillations of conductivity in a two-dimensional electronic system, Phys. Rev. B 49 , 5531 (1994).
- [26] N. S. Averkiev, L. E. Golub, S. A. Tarasenko, and M. Willander, Theory of magneto-oscillation effects in quasi-two-dimensional semiconductor structures, J. Phys.: Condens. Matter 13 , 2517 (2001).
- [27] T. Champel and V. P. Mineev, Magnetic quantum oscillations of the longitudinal conductivity $ \ sigma { zz } $ in quasi-two-dimensional metals, Phys. Rev. B 66 , 195111 (2002).
- [28] P. D. Grigoriev, Theory of the Shubnikov-de Haas effect in quasi-two-dimensional metals, Phys. Rev. B 67 , 144401 (2003).
- [29] I. O. Thomas, V. V. Kabanov, and A. S. Alexandrov, Shubnikov-de Haas effect in multiband quasi-twodimensional metals, Phys. Rev. B 77 , 075434 (2008).
- [30] T. I. Mogilyuk and P. D. Grigoriev, Magnetic oscillations of in-plane conductivity in quasi-two-dimensional metals, Phys. Rev. B 98 , 045118 (2018).
- [31] Z. E. Krix, O. A. Tkachenko, V. A. Tkachenko, D. Q. Wang, O. Klochan, A. R. Hamilton, and O. P. Sushkov, Quantum magnetic oscillations in the absence of closed electron trajectories (2024), arXiv:2404.04592 [condmat].
- [32] R. Landauer, Spatial Variation of Currents and Fields Due to Localized Scatterers in Metallic Conduction, IBM Journal of Research and Development 1 , 223 (1957).
- [33] R. Landauer, Electrical resistance of disordered onedimensional lattices, The Philosophical Magazine: A Journal of Theoretical Experimental and Applied Physics 21 , 863 (1970).

- [34] M. B¨ uttiker, Four-Terminal Phase-Coherent Conductance, Phys. Rev. Lett. 57 , 1761 (1986).
- [35] M. Buttiker, Symmetry of electrical conduction, IBM J. Res. &amp; Dev. 32 , 317 (1988).
- [36] C. W. Groth, M. Wimmer, A. R. Akhmerov, and X. Waintal, Kwant: a software package for quantum transport, New Journal of Physics 16 , 063065 (2014).
- [37] A. Weiße, G. Wellein, A. Alvermann, and H. Fehske, The kernel polynomial method, Rev. Mod. Phys. 78 , 275 (2006).
- [38] J. B. Keller, Corrected bohr-sommerfeld quantum conditions for nonseparable systems, Annals of Physics 4 , 180 (1958).
- [39] M. Wilkinson, An example of phase holonomy in WKB theory, J. Phys. A: Math. Gen. 17 , 3459 (1984).
- [40] D. Ryndyk, Theory of Quantum Transport at Nanoscale: An Introduction , Springer Series in Solid-State Sciences, Vol. 184 (Springer International Publishing, Cham, 2016).
- [41] C. W. Hicks, M. E. Barber, S. D. Edkins, D. O. Brodsky, and A. P. Mackenzie, Piezoelectric-based apparatus for strain tuning, Review of Scientific Instruments 85 (2014).
- [42] J. Diaz, C. Putzke, X. Huang, A. Estry, J. G. Analytis, D. Sabsovich, A. G. Grushin, R. Ilan, and P. J. W. Moll, Bending strain in 3D topological semi-metals, J. Phys. D: Appl. Phys. 55 , 084001 (2021).
- [43] Philip J.W. Moll, personal communication (2022).
- [44] L. Landau, Diamagnetismus der metalle, Zeitschrift f¨r u Physik 64 , 629 (1930).
- [45] V. Leeb, , and J. Knolle, Code and data for 'numerical study of quantum oscillations of the quasiparticle lifetime: Impurity spectroscopy, novel electric field and strain effects', Zenodo online repository (2024).

## Appendix A: Filtering

For the filtering process, we first compute a complexvalued spectrum s F ( ) with a rectangular window, as described in Fig. 2. To filter a frequency f , we scale the spectrum s F ( ) with a narrow Gaussian window exp( -( F -f ) 2 / ∆ f 2 ) centered around f , see Fig. A.1 (a), and compute the inverse Fourier transform of exp( -( F -f ) 2 / ∆ f 2 ) × s F ( ). The real part of the inverse Fourier transform is the filtered oscillation, see Fig. A.1 (b). The absolute (complex phase) of the inverse Fourier transform of exp( -( F -f ) 2 / ∆ f 2 ) × s F ( -f ) are the filtered envelope (phase) of the filtered oscillation, see Fig. A.1 (b) (Fig. A.1 (c)).

## Appendix B: Peak splitting

We filter the frequencies as described in the appendix sec. A. In Fig. A.1 we show that in the case of t ⊥ = 0, i.e., the case of a single degenerate QO frequency, there is a remarkable agreement of the conductance for interband and intraband scattering. Especially, peak splitting is not a result of interband scattering but of scattering in general, as expected from analytics [13]. For inter- and intraband scattering, a phase change around Φ 0 /B = 150 is the origin of the peak splitting. The phase change is continuous and roughly 0 6 . π . It is unclear why the phase change differs from π .

The magnetic field at which the phase jump occurs increases as scattering, i.e., T D increases. This is expected from the zeros of the polynomial amplitudes. E.g. for the first harmonic, the zeros are at T D ≈ eB/ πm 2 . Inserting the zero in R D ≈ 0 04 shows that the Dingle factor . strongly damps the amplitude of QOs around the zero.

FIG. A.1. Phase dependence of the filtered main frequencies and filtering method. (a) Zoom-in of the SdH QO spectra to the first harmonics F 1 = F 2 (gray solid line) for t ⊥ = 0 ( µ = -2 , t I B = (20 250), rectangular window). , The peak is split into two subpeaks. The spectrum is scaled with the Gauss window (gray area) before performing an inverse Fourier transformation. (b) The absolute of the inverse Fourier transformation signal is the filtered signal. Computing the inverse Fourier transformation of the shifted, windowed spectrum yields the filtered envelope, the shaded area in (b), and the filtered phase, shown in (c).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00115v1-figure-a069eaa8-cab5-434e-932e-c2e4daecd0a5.jpg)

FIG. B.2. The filtered phase as function of magnetic field for degenerate frequencies, i.e., t ⊥ = 0 ( µ = -2 ), t and various values of interband scattering. Increasing scattering, i.e., increasing T D moves the phase jump to higher magnetic fields.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00115v1-figure-2e412e77-952c-47fc-a136-8a79d5114214.jpg)
Valentin Leeb Johannes Knolle
Numerical Study of Quantum Oscillations of the Quasiparticle Lifetime: Impurity Spectroscopy, Novel Electric Field and Strain Effects
Quantum oscillation (QOs) measurements constitute one of the most powerful
methods for determining the Fermi surface (FS) of metals, exploiting the famous
Onsager relation between the FS area and the QO frequency. The recent
observation of non-Onsager QOs with a frequency set by the difference of two FS
orbits in a bulk three-dimensional metal can be understood as the QO of the
quasiparticle lifetime (QPL) due to interorbital scattering [Huber, Leeb, {\it
et al.}, Nature 621 (2023)]. QPL oscillations (QPLOs) generalize
magneto-intersubband oscillations (MISOs) known from coupled two-dimensional
metals. They may provide a novel tool for extracting otherwise hard-to-measure
intra-versus interband scattering times of quasiparticles. Here, we provide a
numerical lattice study of QPLOs comparing transport and thermodynamic
observables. We explore the effect of different imperfections like general
impurities, Hall effect-induced electric fields, various forms of strain from
bending, and magnetic field inhomogeneities. We confirm the basic phenomenology
of QPLOs as predicted in analytical calculations and identify additional novel,
non-perturbative features. Remarkably, we find that some imperfections can
stabilize, or even enhance, non-Onsager QPLOs in contrast to standard QO
frequencies. We discuss various avenues for identifying QPLOs in experiments
and how to use their dependence on imperfections to extract material
properties.