---
title: High-Order van Hove Singularities
url: https://www.emergentmind.com/topics/high-order-van-hove-singularity
type: topic
---

# High-Order van Hove Singularities

A high-order van Hove singularity is a momentum-space critical point of a band dispersion at which the usual quadratic saddle-point description fails because the Hessian becomes degenerate or vanishes, so that the leading nonzero terms in the local expansion occur at cubic, quartic, or still higher order. In two dimensions this promotes the familiar logarithmic van Hove divergence of the density of states to a power law, thereby strongly enhancing the phase space for interaction effects. Recent work places high-order van Hove singularities at the intersection of several themes: quasi-flat and exactly flat bands, Lifshitz transitions, moiré minibands, kagome and ruthenate physics, topological Chern bands, and engineered synthetic platforms [2405.20226].

## 1. Definition and local classification

For a conventional two-dimensional van Hove singularity, one expands the band energy around a critical momentum \(k_0\) with \(\nabla \varepsilon(k_0)=0\) and nonvanishing Hessian determinant. In appropriate local coordinates,
\[
\varepsilon(k)-\varepsilon(k_0)\approx \alpha k_x^2-\beta k_y^2,
\]
which yields the standard logarithmic DOS singularity. A high-order van Hove singularity instead occurs when the quadratic description collapses: in the formulation emphasized in the 2024 review, all second derivatives vanish at \(k_0\), \(H=0\), and the first nonzero terms in the Taylor expansion appear at higher orders \(n,m>2\) [2405.20226].

A convenient classification labels the singularity by the pair \((n,m)\) in
\[
\varepsilon(k)=\varepsilon_{\mathrm{VHS}}+\alpha_x k_x^n+\alpha_y k_y^m+\text{(higher-order cross terms)}.
\]
The ordinary saddle is recovered as \((n,m)=(2,2)\), while larger \(n\) and \(m\) encode progressively flatter local dispersion. This \((n,m)\) language is especially useful for relating HOVHS to quasi-flat bands and for comparing different materials platforms [2405.20226].

A complementary classification distinguishes two mechanisms. In the type-I case, both quadratic coefficients vanish and the expansion begins at cubic order, so that three or more Fermi sheets meet at a multicritical saddle point. In the type-II case, exactly one Hessian eigenvalue vanishes while the other remains finite; then the low-energy expansion takes the form
\[
E(p)-E_v=-\alpha p_x^2+\beta p_y^2+\gamma p_x p_y^2+\kappa p_y^4+\cdots,
\]
and the high-order point is reached at \(\beta=0\), where two Fermi sheets touch tangentially rather than crossing at a finite angle [1901.05432]. This distinction is operationally important because type-II HOVHS arise naturally by tuning a single control parameter in moiré graphene and related systems [1901.05432].

## 2. Density of states and singular exponents

The DOS near a HOVHS follows from
\[
\rho(E)=\int d^2k\,\delta[E-\varepsilon(k)].
\]
Using the local form \(\varepsilon-\varepsilon_{\mathrm{VHS}}\sim \alpha_x k_x^n+\alpha_y k_y^m\), one obtains the scaling
\[
\rho(E)\propto |E-E_{\mathrm{VHS}}|^{(1/n+1/m)-1}
\equiv C|E-E_{\mathrm{VHS}}|^{-\mu},
\qquad
\mu=1-\left(\frac1n+\frac1m\right).
\]
The conventional saddle \((2,2)\) sits at the marginal case \(\mu=0\), corresponding to the logarithmic divergence, whereas \((p,2)\) or \((2,p)\) with \(p>2\) already produce algebraic singularities, and \(\mu\to 1\) as \(n,m\to\infty\) [2405.20226].

Several canonical exponents recur across the literature. A type-II HOVHS with \(\beta=0\) in the expansion above yields
\[
\rho(E)\propto |E-E_v|^{-1/4},
\]
with an intrinsically particle-hole-asymmetric peak in twisted bilayer graphene [1901.05432]. A cubic or “monkey-saddle” dispersion, such as \(\varepsilon\sim q_x^3-3q_x q_y^2\), gives \(\rho(E)\sim |E|^{-1/3}\), which appears in multiple settings including honeycomb models, twisted trilayer graphene, Haldane-type models, and kagome-derived band structures [2401.12384]. Quartic HOVHS with \(C_4\)-symmetric dispersion can produce \(\rho(E)\sim |E|^{-1/2}\), while anisotropic \(2+4\) saddles give \(\rho(E)\sim |E|^{-1/4}\) [2205.08828]. In mirror-symmetric twisted trilayer graphene, an additional topological Lifshitz transition leads to an anomalous \(|\omega|^{-2/5}\) singularity when the ratio \(2\gamma/\xi\) reaches unity [2106.14911].

| Local form | DOS behavior | Representative setting |
|---|---|---|
| \((2,2)\) quadratic saddle | \(\rho\sim \ln|E-E_{\mathrm{VHS}}|\) | ordinary VHS |
| \((2,4)\) or \((4,2)\) anisotropic HOVHS | \(\rho\sim |E-E_{\mathrm{VHS}}|^{-1/4}\) | extended saddles, distorted kagome |
| cubic “monkey saddle” | \(\rho\sim |E|^{-1/3}\) | moiré graphene, Haldane, kagome |
| quartic isotropic HOVHS | \(\rho\sim |E|^{-1/2}\) | \(C_4\)-symmetric quartic saddles |
| special Lifshitz critical point | \(\rho\sim |\omega|^{-2/5}\) | mirror-symmetric twisted trilayer graphene |

This hierarchy of exponents formalizes the notion that HOVHS are not a single universality class but a family of critical points whose low-energy singularity is set by the first nonzero homogeneous term in the dispersion [2405.20226].

## 3. Relation to flat bands, Lifshitz structure, and topology

The 2024 review makes the flat-band connection explicit: an exactly flat band corresponds formally to the limit \(n=m\to\infty\), where all derivatives vanish to infinite order and \(\varepsilon(k)\equiv \text{const}\) in a neighborhood of \(k_0\). In practice, finite but large \(n\) and \(m\) produce a quasi-flat local pocket. Quantitatively, the local bandwidth shrinks as
\[
W\sim \Lambda^{\min(n,m)},
\]
while the DOS peak height grows as
\[
\rho_{\max}\sim \Lambda^{-2+1/n+1/m}.
\]
This places HOVHS as an interpolation between ordinary saddles and true flat bands rather than as a disconnected phenomenon [2405.20226].

Because a van Hove singularity marks a Lifshitz transition, HOVHS typically encode the merging, annihilation, or higher-order tangency of Fermi-surface sheets. In twisted bilayer graphene, the type-II scenario describes the coalescence of ordinary VHS into tangential touching points as a single tuning parameter is varied [1901.05432]. In Bernal bilayer and rhombohedral trilayer graphene, a one-parameter interpolation from three separate ordinary VHS to a single cubic HOVHS captures the Lifshitz transition from three pockets to one larger pocket, with distinct parquet-RG fixed-point structures on the two sides of the crossover [2401.12384].

Topological band structures supply further variants. Mirror-symmetric twisted trilayer graphene hosts a zero-energy HOVHS with exponent \(-1/3\), protected by \(C_3\) rotation and a combined mirror-particle-hole symmetry, and tunable by twist angle and perpendicular electric field [2106.14911]. Topological moiré surface states on three-dimensional topological insulators realize \(C=\pm 1\) Chern bands whose valley saddles merge into a cubic HOVHS, and the low-temperature intrinsic anomalous Hall response satisfies
\[
\frac{d\sigma_{xy}^{\mathrm{int}}}{d\mu}\propto \rho(\mu)\langle \Omega\rangle_{\mathrm{FS}},
\]
so the Hall anomaly inherits the HOVHS power law [2402.16772]. Kagome topological bands with complex next-nearest-neighbor hopping similarly exhibit HOVHS with exponents \(\nu=1/2,1/3,1/4\), alongside Chern numbers ranging from \(C=\pm 1\) to \(\pm 4\) [2410.07000]. This suggests that HOVHS often mediate between band flattening and topological response rather than belonging exclusively to either category.

## 4. Interaction effects and ordered phases

The central many-body consequence of a HOVHS is the enhancement of any susceptibility that is weighted by the DOS. In the review formulation, even a simple mean-field criterion such as \(U\rho(E_F)>1\) becomes easier to satisfy when \(\rho(E)\) diverges with \(\mu>0\). Candidate orders include ferro- and antiferromagnetism, charge- and spin-density waves, Pomeranchuk distortions, and unconventional superconductivity; for \((n,m)=(3,3)\), the pairing susceptibility diverges as \(T^{-1/3}\) in the normal state, which strongly boosts weak-coupling pairing tendencies [2405.20226].

More detailed treatments corroborate that general picture. In magic-angle twisted bilayer graphene, the intervalley density-wave susceptibility diverges faster than \(\ln(1/T)\), and weak strain or intervalley interactions can split the VHS peak into two, consistent with scanning-tunneling observations of peak splitting under doping [1901.05432]. In a honeycomb-lattice extended Hubbard model, determinant and constrained-path quantum Monte Carlo locate the HOVHS at the critical hopping relation
\[
t''_c=\frac{t-2t'}{4},
\]
where the quadratic mass tensor collapses in one direction and the DOS becomes \(\rho(E)\sim |E|^{-1/3}\). Near the corresponding filling, \(\langle n\rangle_{\mathrm{HOVH}}\simeq 0.881\) for \(t=1.0\), \(t'=0.10\), \(t''=0.20\), the system exhibits a crossover between ferromagnetic and antiferromagnetic fluctuations, and long-distance \(f_n\)-wave correlations dominate the pairing sector; nearest-neighbor Coulomb interactions suppress that pairing approximately in proportion to \(|V|\) [2505.10358].

Renormalization-group studies reach similar conclusions from a different direction. For a model with nested Fermi-surface patches plus a quartic HOVHS patch, the additional power-law bubbles can raise the spin-density-wave critical temperature from \(T_c\approx 7\times 10^{-7}\) to \(T_c\approx 4\times 10^{-5}\) for the quoted bare couplings, i.e. by roughly two orders of magnitude, and can also favor charge-density-wave order for other couplings [2205.08828]. On the distorted kagome surface of Co\(_3\)Sn\(_2\)S\(_2\), a quartic HOvHS pinned to the Fermi energy is reported to precipitate a Pomeranchuk instability and a cascade of nematic states over an energy shell of about \(100\) meV, without generating additional translational symmetry breaking [2410.01994]. More specialized analyses further argue that an \(X_9\) quartic HOVHS with \(\rho(\varepsilon)\propto |\varepsilon|^{-1/2}\) can support triplet superconductivity with \(T_c(U)\propto U^2\), with an upper estimate of about \(40\) mK for Sr\(_3\)Ru\(_2\)O\(_7\) in the quoted parameter set [2603.22645].

## 5. Materials platforms and experimental signatures

The materials survey in the 2024 review identifies twisted bilayer and multilayer graphene near magic angles, ABC-stacked graphene trilayers under displacement field, moiré transition-metal dichalcogenides, kagome metals, and photonic or cold-atom simulators as established or promising HOVHS platforms [2405.20226]. In magic-angle twisted bilayer graphene, scanning tunneling spectroscopy observes an asymmetric DOS peak whose log-log plot shows two parallel lines of slope \(-1/4\), with peak width of tens of meV; tuning by twist angle or hydrostatic pressure moves the system through the critical regime where the type-II HOVHS emerges [1901.05432]. The same work gives a critical coupling \(g_c\approx 1.995\) for \(g'/g=1.2\), corresponding to \(\theta_c\approx 0.95^\circ\)–\(1.0^\circ\) for the quoted continuum-model parameters [1901.05432].

ABC-stacked trilayer graphene aligned with hBN provides a field-tuned type-I realization: at \(V_c\approx 0.06\) eV the quadratic terms vanish and three Fermi pockets meet at \(K'\) [1901.05432]. Mirror-symmetric twisted trilayer graphene hosts a zero-energy HOVHS tunable by twist angle and perpendicular field; for realistic corrugation \(r\approx 0.8\) and \(\theta\approx 1.59^\circ\), one finds \(U_c/\varepsilon_0\approx 0.37\) V/nm, and scanning tunneling spectroscopy is predicted to distinguish the \(|\omega|^{-1/3}\) and \(|\omega|^{-2/5}\) regimes [2106.14911].

Kagome-derived systems show several distinct mechanisms. The review notes that purely two-dimensional kagome dispersion can support \((n,m)=(2,4)\) or higher saddles at high-symmetry points [2405.20226]. On the distorted Co\(_3\)Sn\(_2\)S\(_2\) surface, scanning tunneling spectroscopy reveals a sharp zero-bias peak well fitted by \(|E|^{-1/4}\) after subtracting a constant background, and quasiparticle interference at the Fermi energy retains only one \(\Gamma\)–\(M\) scattering branch, consistent with a nematic reconstruction driven by the HOvHS [2410.01994]. Bilayer kagome borophene has been proposed to host both a conventional and a monkey-saddle VHS in the same band: the conventional one lies \(0.065\) eV below the Fermi level, the high-order one \(0.385\) eV below, and the adjacent Dirac-like cone reaches a Fermi velocity of \(1.34\times 10^6\) m/s [2307.07137].

Synthetic settings broaden the scope beyond solid-state systems. Photonic and cold-atom lattices can realize tailored \((n,m)\) dispersions with tunable exponents [2405.20226]. In a checkerboard optical-lattice insulator designed for cavity polaritons, the interband gap edge realizes an effective HOVHS in the joint DOS, with the divergence engineered without sub-gap absorption by maintaining a full direct gap and using spin-polarized noninteracting fermions [2509.15849]. Such examples indicate that the relevant singular object need not be only the electronic DOS; analogous constructions can target the joint DOS or cavity self-energy.

## 6. Detection, engineering, and open directions

A general method for diagnosing HOVHS in multiband Hamiltonians is provided by the extended Feynman-Hellmann framework. Rather than diagonalizing a model analytically near every candidate momentum, one computes derivatives of the band eigenvalue along a path \(k(\lambda)=k_0+\lambda\Delta k\) using recursive formulas for \(\partial^M \varepsilon_n\). This yields the Taylor expansion directly and supports a concrete workflow: identify critical points from \(\nabla \varepsilon_n(k_0)=0\), inspect the Hessian, test whether the third-order tensor has support on null directions, and continue to quartic or higher order when necessary [2207.06099]. The same formalism incorporates tuning parameters by promoting the path to \((k(\lambda),p(\lambda))=(k_0+\lambda\Delta k,p_0+\lambda\delta p)\), so that the coefficient whose zero defines the HOVHS can be solved as a function of strain, bias, staggered potential, or other control parameters [2207.06099].

The engineering program already has concrete case studies. In the Haldane model, tuning the staggered potential to
\[
M_0=\frac{t_1^2-18t_2^2}{2\sqrt{3}\,t_2}
\]
produces a monkey saddle at \(K_-\) with
\[
\varepsilon_+(K_-+q)\simeq \varepsilon_+(K_-)-\frac{3\sqrt{3}t_2}{2}\,[q_y^3-3q_x^2 q_y]+O(q^4),
\]
and small detuning restores quadratic curvature and splits the multicritical point into ordinary extrema and saddles [2207.06099]. In the surface layer of Sr\(_2\)RuO\(_4\), tight-binding models constrained by ARPES and QPI show how octahedral rotation and weak nematicity move an ordinary VHS toward an \(A_3\)-type HOVHS; the bare DFT model yields \(\theta_c\approx 9.7^\circ\), while the experimentally refined surface lies close to the critical regime [2310.15331]. Non-Hermitian Floquet interfaces offer a different route: when an exceptional ring grazes the Fermi line, paired ordinary VHS can coalesce into a single higher-order DOS peak with power-law exponents \(-1/2\) and \(-2/3\), a mechanism explicitly distinguished from Hermitian saddle merging [2309.11565].

The present literature therefore frames HOVHS as a unifying concept rather than a niche singularity class. It links local band flattening to experimentally tunable Lifshitz transitions, provides a systematic route from ordinary saddles to quasi-flat bands, and organizes a broad range of interaction-driven phenomena across moiré graphene, kagome systems, ruthenates, Chern bands, and synthetic lattices [2405.20226]. A plausible implication is that future progress will depend less on identifying isolated examples and more on building controlled classification, detection, and tuning schemes across multiband and topological settings.

Source: https://www.emergentmind.com/topics/high-order-van-hove-singularity