---
title: Majorana dHvA in Strained Kitaev Liquids
url: https://www.emergentmind.com/topics/majorana-de-haas-van-alphen-effect
type: topic
---

# Majorana dHvA in Strained Kitaev Liquids

Searching arXiv for the cited paper and directly related references.
arXiv Search Query: 2512.12570
The Majorana de Haas–van Alphen effect denotes a strain-driven quantum-oscillation phenomenon for charge-neutral Majorana quasiparticles in a Kitaev spin liquid with a Majorana Fermi surface. In the realization studied in "Strain-induced quantum oscillation in Kitaev spin liquid with Majorana-Fermi surface" [2512.12570], triaxial lattice strain generates an effective vector potential and a uniform pseudo-magnetic field that Landau-quantize itinerant Majorana fermions. The resulting pseudo-Landau levels produce oscillations of the density of states and the specific heat at very low temperatures, in close analogy to the conventional de Haas–van Alphen effect in metals, but with a neutral, valley-odd, strain-induced gauge coupling rather than minimal coupling to an electromagnetic field.

## 1. Microscopic setting in the Kitaev model

The starting point is the isotropic spin-\(1/2\) Kitaev model on the honeycomb lattice,
\[
H_K=-J\sum_{\langle ij\rangle_\gamma} S_i^\gamma S_j^\gamma,
\]
where \(\gamma\in\{x,y,z\}\) labels bond types and \(\langle ij\rangle_\gamma\) are nearest neighbors on \(\gamma\)-bonds. In the Majorana representation used in the paper, the spin operator is written as
\[
S_i^\alpha = i\, b_i^\alpha c_i,
\]
with one itinerant Majorana \(c_i\) and three gauge Majoranas \(b_i^\alpha\) per site, together with a static \(\mathbb{Z}_2\) link variable \(u_{ij}=ib_i^\alpha b_j^\alpha=\pm 1\). In a fixed flux sector, and in particular in the flux-free ground state, the model reduces to free Majorana hopping,
\[
H_K=i\sum_{\langle ij\rangle} J_{ij}u_{ij}c_i c_j,
\]
which yields Dirac cones at the honeycomb \(K\) and \(K'\) points for \(u_{ij}=+1\) [2512.12570].

The paper considers perturbations that generate Majorana Fermi surfaces near these Dirac points. Two microscopic routes are identified. A staggered Zeeman field \(h_{\mathrm{st}}\), associated with proximity to weak antiferromagnetism, induces an effective three-spin term that becomes next-nearest-neighbor Majorana hopping at low energies. A uniform electromagnetic field with a cross term \(E\cdot B\) likewise generates a three-spin interaction that maps to next-nearest-neighbor Majorana hopping. In both cases, time-reversal and inversion symmetries are broken.

In momentum space, these perturbations shift the Dirac cones at \(K\) and \(K'\) in opposite directions in energy: the \(K\) cone moves up and the \(K'\) cone moves down by a valley-dependent offset \(\Delta\varepsilon\propto \pm h_{\mathrm{st}}^3\). This produces small electron-like and hole-like Majorana pockets around \(K\) and \(K'\), respectively. The corresponding low-energy linearized Hamiltonian near valley \(\tau=\pm\) is
\[
H_\tau = v_F(p_x\sigma_x+p_y\sigma_y)+\tau\,\Delta,
\]
where \(v_F\) is the Dirac velocity, \(\sigma_{x,y}\) act in sublattice space, and \(\Delta\propto h_{\mathrm{st}}^3/J^2\). The simultaneous presence of \(+\Delta\) and \(-\Delta\) at zero energy is the immediate origin of the Majorana Fermi pockets.

## 2. Strain as an axial gauge field

The quantizing field in the Majorana de Haas–van Alphen effect is not an external magnetic field but a pseudogauge field generated by lattice strain. In the model, strain modulates the Kitaev couplings on each nearest-neighbor bond according to
\[
J_{ij}=J\left[1-b\left(\frac{|\boldsymbol{\delta}_{ij}|}{a_0}-1\right)\right],
\]
with magnetoelastic constant \(b\), lattice constant \(a_0=1\), and bond vector
\[
\boldsymbol{\delta}_{ij}=\boldsymbol{R}_i+\boldsymbol{U}_i-\boldsymbol{R}_j-\boldsymbol{U}_j.
\]
In the Dirac continuum limit, the nonuniform \(J_{ij}\) appears as an effective valley-odd gauge field \(\boldsymbol{A}_{\mathrm{eff}}(\boldsymbol{r})\) that minimally couples to the Majorana Dirac Hamiltonian, in direct analogy to strained graphene [2512.12570].

The effective vector potential has the form
\[
\boldsymbol{A}_{\mathrm{eff}} \propto (\epsilon_{xx}-\epsilon_{yy},-2\epsilon_{xy}),
\]
where
\[
\epsilon_{ij}=\tfrac{1}{2}(\partial_i U_j+\partial_j U_i)
\]
is the strain tensor and \(\boldsymbol{U}(\boldsymbol{r})\) is the displacement field. The specific strain pattern used is the triaxial profile
\[
\boldsymbol{U}(x,y)=C\,(2xy,x^2-y^2),
\]
with strain amplitude \(C\). This yields
\[
\boldsymbol{A}_{\mathrm{eff}}(x,y)=C\,(4y,-4x),
\]
and therefore a uniform pseudo-magnetic field perpendicular to the plane,
\[
B_{\mathrm{eff}}=(\nabla\times \boldsymbol{A}_{\mathrm{eff}})_z=-8C,
\]
so that \(|B_{\mathrm{eff}}|=8C\) in lattice units.

A defining feature is the axial character of the coupling: quasiparticles near \(K\) and \(K'\) experience opposite signs of \(B_{\mathrm{eff}}\). The strained system therefore preserves overall time-reversal invariance while still exhibiting Landau quantization within each valley. This valley-odd structure is central to the Majorana version of de Haas–van Alphen physics, because it produces mirror-related electron-like and hole-like sectors rather than a single charged carrier species.

## 3. Landau quantization of Majorana quasiparticles

Under a uniform axial field, the low-energy Majorana Dirac fermions form discrete pseudo-Landau levels. For each valley \(\tau=\pm\), the spectrum is
\[
E_{n,\tau}=\tau\,\Delta \pm v_F\sqrt{2|n|\,\hbar\,|B_{\mathrm{eff}}|},
\qquad n=0,1,2,\dots .
\]
The zeroth Landau level is pinned at
\[
E_{0,\tau}=\tau\,\Delta.
\]
When \(\Delta=0\), the \(n=0\) level sits at zero energy; with finite \(\Delta\), it shifts valley-dependently but persists under triaxial strain. The paper’s local density of states shows that the zero-mode peak survives and moves slightly away from zero when \(\Delta\neq 0\) [2512.12570].

For sufficiently large \(|\Delta|\), the dispersion near the small pockets becomes approximately parabolic, and a nonrelativistic description is appropriate:
\[
E_{n,\tau}^{\rm NR}\simeq \hbar\omega_{c,\tau}\left(n+\tfrac{1}{2}\right)-\mu_\tau,
\]
with cyclotron frequency \(\omega_{c,\tau}=|B_{\mathrm{eff}}|/m_\tau^*\) and valley-dependent chemical potential \(\mu_\tau\). The paper also gives the areal degeneracy per Landau level as
\[
\mathcal{N}_{\mathrm{LL}}=\frac{A}{2\pi \ell_B^2},\qquad \ell_B^{-2}\equiv |B_{\mathrm{eff}}|,
\]
for each valley, where \(A\) is the sample area.

Because Majorana fermions are real rather than complex, the spectrum is symmetric in \(\pm E\), and thermodynamic counting involves both signs. In the flux-free sector the Landau levels are sharp, whereas in the two-flux sector bound states appear between Landau levels. Semiclassically, the oscillations are governed by Onsager quantization,
\[
S(E)=2\pi (n+\gamma)\,\frac{\hbar}{\ell_B^2}
      =2\pi (n+\gamma)\,\hbar |B_{\mathrm{eff}}|,
\]
where \(S(E)\) is the extremal orbit area in momentum space. The phase offset is \(\gamma=0\) for Dirac pockets because of the Berry phase, and \(\gamma=\tfrac{1}{2}\) for nonrelativistic pockets.

## 4. Oscillatory thermodynamics and the de Haas–van Alphen analogy

The defining observable consequence of the pseudo-Landau quantization is an oscillatory contribution to the density of states and to the specific heat as functions of \(1/|B_{\mathrm{eff}}|\), or equivalently \(1/C\). The paper describes this using a Lifshitz–Kosevich-type form adapted to the pseudo-field,
\[
\delta \rho(E_F,B_{\mathrm{eff}})\propto
\sum_{p=1}^{\infty} R_T(p)\,R_D(p)\,
\cos\!\left(2\pi p\,\frac{F}{|B_{\mathrm{eff}}|}+\phi\right),
\]
with harmonic index \(p\), phase \(\phi\), temperature damping
\[
R_T(p)=\frac{X_p}{\sinh X_p},
\qquad
X_p\equiv \frac{2\pi^2 k_B T\,p}{\hbar\,\omega_c},
\]
and Dingle factor
\[
R_D(p)=\exp\!\left(-\frac{2\pi^2 k_B T_D\,p}{\hbar\,\omega_c}\right).
\]
The oscillatory part of the specific heat follows from the grand potential through
\[
\delta C=-T\,\partial^2 \delta\Omega/\partial T^2,
\]
and inherits the same periodicity and damping factors [2512.12570].

The fundamental frequency is set by the extremal area of the Majorana Fermi surface at zero energy,
\[
F\equiv \frac{S(E_F)}{2\pi},
\qquad E_F=0.
\]
This differs from the conventional metallic expression \(F=(\hbar/2\pi e)S(E_F)\): for neutral Majoranas the effective coupling is to the strain-generated pseudofield, so \(e\) does not appear, and quantization is controlled by \(\ell_B^{-2}=|B_{\mathrm{eff}}|\). In the numerical results, the density of states at \(\varepsilon=0\) shows clear oscillations versus strain \(C\) when \(h_{\mathrm{st}}>0\), and the very-low-\(T\) specific heat \(C_v(T)\) likewise oscillates with \(C\), with increasing amplitude as \(h_{\mathrm{st}}\) grows.

A common misconception is that de Haas–van Alphen-type oscillations require charged carriers. In this setting, the oscillatory response is instead driven by a pseudogauge field derived from strain; no Lorentz force on electric charge is involved. The analogy to conventional dHvA is therefore thermodynamic and semiclassical rather than electromagnetic.

## 5. Distinctive features relative to conventional dHvA

Several features distinguish the Majorana effect from ordinary electron dHvA while preserving the central logic of Landau quantization of closed Fermi-surface orbits [2512.12570].

| Aspect | Majorana de Haas–van Alphen effect | Conventional de Haas–van Alphen effect |
|---|---|---|
| Source of quantizing field | Pseudo-magnetic field \(B_{\mathrm{eff}}\) from triaxial strain; axial coupling with opposite signs in valleys | Real magnetic field \(B\); minimal coupling with electric charge \(e\) |
| Carriers | Charge-neutral Majorana quasiparticles; spectrum symmetric in \(\pm E\); Dirac zero-mode persists | Charged electrons; no zero-mode in parabolic bands |
| Frequency | \(F=S(E_F)/(2\pi)\) with \(E_F=0\); periodic in \(1/|B_{\mathrm{eff}}|\) | \(F=(\hbar/2\pi e)S(E_F)\); periodic in \(1/B\) |
| Phase | Berry phase \(\gamma=0\) for Dirac pockets; valley-dependent contributions and possible bound-state features in flux sectors | \(\gamma=\tfrac{1}{2}\) for parabolic bands; \(\gamma=0\) for Dirac electrons |
| Tuning parameter | Strain amplitude \(C\) | Magnetic field \(B\) |

The Dirac zero-mode is especially distinctive. In the strained Majorana system a robust zeroth Landau level remains present, shifted by \(\Delta\) when the perturbation-induced pocket offset is finite. The axial coupling is likewise nonstandard: the two valleys contribute with opposite pseudo-field signs, so the oscillatory response combines electron-like and hole-like pockets that are mirror images in spectrum but identical in level spacing.

The paper also notes that valley mixing and flux-sector bound states can modulate the oscillation phase and amplitude. This places the Majorana effect closer to a topological Dirac-Landau problem than to the textbook case of a parabolic electron band.

## 6. Quantitative scale, measurement protocols, and limitations

For the triaxial profile
\[
\boldsymbol{U}=(2Cxy,C(x^2-y^2)),
\]
the pseudo-field magnitude is \(|B_{\mathrm{eff}}|=8C\), and the magnetic length is
\[
\ell_B=1/\sqrt{|B_{\mathrm{eff}}|}.
\]
A representative numerical choice is flake radius \(r=20\) and \(C=0.01\), which gives \(|B_{\mathrm{eff}}|\approx 0.08\) and \(\ell_B\approx 3.5a_0\). The Dirac Landau-level spacing near \(n=1\) scales as
\[
\Delta E_1 \approx v_F\sqrt{2\hbar |B_{\mathrm{eff}}|}\propto v_F\sqrt{C}.
\]
Taking \(v_F\sim O(J)\) and \(J\sim 5\)–10 meV in candidate materials as an order-of-magnitude estimate yields \(\Delta E_1\sim 2\)–6 meV for \(C\sim 0.01\), corresponding to 20–70 K. In the calculations, the oscillations in \(D(0)\) and \(C_v\) are resolved at very low temperature, with \(T=10^{-3}J\ll \Delta E_1\) [2512.12570].

The oscillation period in inverse pseudo-field is
\[
\Delta\!\left(\frac{1}{|B_{\mathrm{eff}}|}\right)=\frac{2\pi}{S(E_F)}.
\]
Because \(S(E_F)\) is controlled by the pocket size near \(K/K'\), and the pocket size is controlled by \(\Delta\propto h_{\mathrm{st}}^3/J^2\), increasing \(h_{\mathrm{st}}\) enlarges the pockets and increases the oscillation amplitude. Tiny pockets imply a small \(S(E_F)\) and therefore slow oscillations in \(1/C\); larger pockets imply faster oscillations.

The proposed measurements are specific heat versus strain amplitude \(C\), seeking periodic oscillations in \(1/C\) at very low temperature, and STM measurements of the local density of states in strained monolayer flakes, seeking pseudo-Landau levels near the flake center. Nanobubble-like triaxial strain profiles are proposed as a way to produce a nearly uniform \(B_{\mathrm{eff}}\) in a central region. Smooth phonon contributions to specific heat provide a background that varies smoothly in \(C\) and \(T\), without \(1/C\) periodicity; magnetic backgrounds such as magnons are not expected to show axial-field-induced quantization at zero applied magnetic field. Valley-dependent electron-like and hole-like pockets and the shifted Dirac zero-mode are therefore identified as distinctive signatures.

The analysis is subject to explicit assumptions and limitations. The flux-free sector is assumed to dominate at low \(T\), with flux excitations neglected except when explicitly studied. The model is isotropic, and Heisenberg and \(\Gamma\) interactions and other perturbations are ignored. The calculations use open-boundary hexagonal flakes and impose \(C<0.3/r\) to avoid sign changes of \(J_{ij}\) at the edges. The mapping to \(\boldsymbol{A}_{\mathrm{eff}}\) is a linearized Dirac description valid for small strains and smooth displacement fields, and the precise mapping from magnetoelastic coefficients to \(B_{\mathrm{eff}}\) is material dependent. Within those assumptions, the Majorana de Haas–van Alphen effect provides a route to probing a charge-neutral Majorana Fermi surface through strain-controlled Landau quantization and the resulting oscillatory density of states and thermodynamics.

Source: https://www.emergentmind.com/topics/majorana-de-haas-van-alphen-effect