Bound states of quasiparticles with quartic dispersion in an external potential: WKB approach

This paper formulates a WKB approach for quasiparticles with quartic dispersion, demonstrating that higher-order Airy-type functions and their hyperasymptotic corrections are essential for matching wave functions at turning points, leading to a generalized Bohr-Sommerfeld quantization condition that includes non-perturbative corrections even in the absence of tunneling.

The Gist

Here is an explanation of the paper "Bound states of quasiparticles with quartic dispersion in an external potential: WKB approach," translated into simple, everyday language with creative analogies.

The Big Picture: A New Kind of "Bouncy Ball"

Imagine you are a physicist studying tiny particles (quasiparticles) trapped inside a valley created by a force field (a potential). In the standard world of quantum mechanics (like electrons in a normal metal), these particles behave like bouncy balls. If you throw a ball, its energy goes up as the square of its speed ($E \sim v^2$). This is the "quadratic" world we are used to.

But in this paper, the authors are studying a very strange, exotic world—like ABC-stacked graphene (a super-thin material made of carbon layers). In this world, the particles don't just bounce; they behave like heavy, sluggish trucks or marshmallows. Their energy doesn't go up with the square of the speed; it goes up with the fourth power ($E \sim v^4$).

This "quartic" behavior makes the math much harder. It's like trying to predict the path of a truck that gets heavier the faster it goes, compared to a normal ball.

The Problem: The "Foggy" Turning Points

To figure out where these particles can stay trapped (bound states), physicists use a tool called the WKB method. Think of this as a map-making technique.

• The Allowed Region: Where the particle has enough energy to roll around. On the map, this is clear, sunny terrain.
• The Forbidden Region: Where the particle doesn't have enough energy. This is a deep fog or a cliff edge.
• The Turning Point: The exact spot where the particle stops rolling and turns back.

In the normal "bouncy ball" world, the map is easy to draw. You just use a standard tool (the Airy function) to connect the sunny side to the foggy side. It's like using a standard bridge to cross a river.

But in this "quartic truck" world, the standard bridge doesn't work. The math gets messy near the turning point. The authors realized that to cross this river, they needed a super-bridge made of a more complex material called Fourth-Order Airy Functions.

The Secret Sauce: The "Ghost" Paths (Hyperasymptotics)

Here is the most surprising part of the paper.

In the normal world, when you are in the "sunny" allowed region, the particle just oscillates back and forth. It's like a pendulum swinging.
However, in this "quartic truck" world, the math says there are also "ghost" paths. Even in the sunny region, there are invisible, exponentially decaying waves (like a faint echo) and exponentially growing waves (like a whisper getting louder).

Usually, physicists ignore these ghosts because they seem too small to matter. But the authors discovered that in this specific "quartic" world, ignoring the ghosts leads to the wrong answer.

They used a mathematical technique called Steepest Descents (imagine finding the steepest path down a mountain to get to the bottom quickly) to analyze these "ghost" paths. They found that these tiny, hidden contributions are actually crucial. They call this "Hyperasymptotics."

The Analogy:
Imagine you are trying to tune a radio to a specific station.

• Standard Physics: You turn the dial until you hear the music clearly.
• This Paper: You turn the dial, but there is a tiny, almost silent static noise (the ghost path) underneath the music. If you ignore the static, your radio is slightly out of tune. If you account for the static (the hyperasymptotics), the music becomes crystal clear.

The Result: A New Rulebook

By properly connecting the "sunny" side and the "foggy" side using these complex "super-bridges" and accounting for the "ghost" paths, the authors derived a new rulebook (a quantization condition) for calculating the energy levels of these particles.

This new rulebook is a generalization of the famous Bohr-Sommerfeld rule (the old rulebook for bouncy balls).

• The Old Rule: "Count the waves, add a half-step, and you have your energy."
• The New Rule: "Count the waves, add a half-step, AND add a tiny, invisible correction term that comes from the ghost paths."

Why Does This Matter?

1. It's Accurate: When they tested this new rule on a "harmonic oscillator" (a simple valley shape), they found that the old rule was off by about 11% for the lowest energy state. The new rule, with the "ghost" correction, matched the exact answer almost perfectly.
2. It's Universal: This isn't just about graphene. It applies to any system where particles behave like these "quartic trucks."
3. It Changes How We Think: It proves that even in places where a particle should be moving freely (the allowed region), there are hidden, decaying waves that are essential for the particle to stay trapped. You can't just look at the main wave; you have to listen to the whispers.

Summary in One Sentence

The authors figured out how to accurately predict the energy of exotic, slow-moving particles in graphene by realizing that even in the "allowed" zones, tiny, invisible "ghost" waves are secretly holding the whole system together, and ignoring them leads to big mistakes.

Technical Summary

Here is a detailed technical summary of the paper "Bound states of quasiparticles with quartic dispersion in an external potential: WKB approach" by E.V. Gorbar and V.P. Gusynin.

1. Problem Statement

The paper addresses the challenge of determining bound state energies for quasiparticles with quartic energy-momentum dispersion ($E(p) \sim p^4$) moving in an external potential $V(x)$.

• Context: Such systems arise in condensed matter physics, specifically in ABC-stacked $N$-layer graphene (where $N=4$), where the low-energy dispersion is $E(p) \propto p^N$.
• Mathematical Challenge: Unlike the standard quadratic dispersion ($E \sim p^2$) which leads to a second-order Schrödinger equation, quartic dispersion leads to a fourth-order ordinary differential equation (ODE).
• Specific Difficulty: The standard Wentzel-Kramers-Brillouin (WKB) semiclassical method requires matching wave functions across "turning points" (where $E = V(x)$). For quadratic dispersion, this is done using Airy functions. For quartic dispersion, the matching requires higher-order Airy functions. Crucially, the presence of quartic dispersion introduces purely imaginary momenta even in classically allowed regions, leading to exponentially growing and decaying components in the wave function alongside the standard oscillatory terms. Standard WKB matching fails to account for these non-oscillatory components correctly without advanced asymptotic analysis.

2. Methodology

The authors employ a semiclassical WKB approach extended to fourth-order differential equations, focusing on the connection formulas at turning points.

• WKB Expansion: The wave function is assumed in the form $\psi(x) = \exp[iS(x)/\hbar]$. Expanding the action $S$ in powers of $\hbar$ yields recursive relations for the momentum $p(x)$ and correction terms $S_n$.
• Linearization at Turning Points: Near a turning point $x_0$, the potential is linearized ($V(x) - E \approx C(x-x_0)$). This reduces the fourth-order Schrödinger-like equation to a canonical form:
  $$\psi^{(4)}(z) + z\psi(z) = 0$$
  where $z$ is a scaled coordinate.
• Higher-Order Airy Functions: The solutions to the above equation are identified as fourth-order Airy functions ($Ai_4$ and $\tilde{Ai}_4$). These are defined via Laplace-type integrals over contours in the complex plane.
• Steepest Descent & Hyperasymptotics:
  • The authors analyze the asymptotic behavior of these functions as $z \to \pm \infty$ using the method of steepest descents.
  • A critical innovation is the calculation of hyperasymptotics: exponentially suppressed terms (e.g., $e^{-1/\hbar}$) that are usually neglected in standard asymptotic expansions.
  • These terms are essential because they arise from saddle points that are exponentially distant but contribute to the correct matching of the exponentially decaying/growing components in the classically allowed region.
• Matching Procedure: The wave functions in the classically allowed region (oscillatory + exponential) and forbidden region (exponential decay) are matched using the asymptotic forms of the fourth-order Airy functions. This yields a system of equations linking the coefficients of the wave function.

3. Key Contributions

1. Generalized Quantization Condition: The authors derive a generalized Bohr-Sommerfeld quantization condition for quartic dispersion:
   $$\frac{1}{\hbar a} \int_{x_a}^{x_b} (E - V(x))^{1/4} dx = 2(-1)^n \arctan\left( \frac{1}{2} e^{-I/\hbar} \right) + \pi\left(n + \frac{1}{2}\right)$$
   where $I$ is the phase integral.
2. Non-Perturbative Correction: The term involving the arctangent represents a non-perturbative correction in $\hbar$ (specifically, exponentially small in $1/\hbar$).
   • Unlike the standard case where such corrections arise from tunneling between complex turning points, this correction appears here even for real turning points (e.g., in a harmonic potential) due to the specific structure of the quartic dispersion.
   • This correction is directly linked to the hyperasymptotic terms of the fourth-order Airy functions.
3. Handling of Exponential Components: The paper demonstrates that for quartic dispersion, the wave function in the classically allowed region must include exponentially increasing and decreasing terms to satisfy boundary conditions, a feature absent in quadratic dispersion.
4. Analytic Regularization: The authors address the divergence of higher-order WKB corrections (specifically the second-order term $S_2$) at turning points by using analytic regularization and contour deformation, allowing for the inclusion of higher-order perturbative corrections.

4. Results

The derived quantization condition was applied to two specific systems:

• Harmonic Potential ($V \propto x^2$):
  • The system is equivalent via Fourier transform to a Schrödinger equation with a quartic potential ($V \propto k^4$).
  • Numerical results show that including the hyperasymptotic correction significantly improves the accuracy of the ground state energy ($n=0$), reducing the error from ~11% (leading order) to ~3% (compared to exact numerical values).
  • For higher energy levels ($n > 1$), the correction becomes negligible, and the leading-order WKB approximation is already highly accurate.
• Double Quartic System ($H = a^4 p^4 + b^4 x^4$):
  • The authors calculated bound state energies for this system.
  • Similar to the harmonic case, the non-perturbative correction is most significant for the ground state ($n=0$), bringing the WKB estimate ($\epsilon \approx 1.4241$) closer to direct numerical calculations ($\epsilon \approx 1.3967$).

Table I & II Analysis: The paper provides tables comparing:

• $\epsilon_{lead}$: Leading order WKB.
• $\epsilon_{hyper}$: Leading order + hyperasymptotics.
• $\epsilon_{2nd}$: Leading + 2nd order perturbative.
• $\epsilon_{2ndhyper}$: Leading + 2nd order + hyperasymptotics.
• $\epsilon_{exact}$: Exact numerical values.
  The results confirm that hyperasymptotics are crucial for low-lying states, while higher-order perturbative corrections dominate for excited states.

5. Significance

• Theoretical Advancement: This work extends the WKB formalism to higher-order differential equations, a field where global properties and connection formulas were previously underdeveloped. It establishes the necessity of hyperasymptotics for correct matching in systems with soft dispersion relations ($E \sim p^{2n}, n \ge 2$).
• Physical Insight: It reveals that "bound states in the continuum" (BIC) or specific tunneling-like effects can occur naturally in quartic dispersion systems without fine-tuning, due to the presence of decaying components in classically allowed regions.
• Material Science Application: The results are directly applicable to multilayer graphene (specifically ABC-stacked tetralayer), providing a tool to calculate bound state energies and transport properties in these materials where quartic dispersion is dominant.
• Generalizability: The methodology can be extended to higher-order dispersions ($p^6, p^8$, etc.) and higher spatial dimensions, provided the corresponding higher-order Airy functions and their hyperasymptotics can be determined.

In summary, the paper provides a rigorous semiclassical framework for quartic dispersion systems, highlighting that non-perturbative corrections arising from hyperasymptotics are essential for accurate ground-state energy predictions, even in the absence of complex turning points.