Nonlinear Landau levels in the almost-bosonic anyon gas

The gist: The authors derive and analyze an effective Chern–Simons–Schrödinger (CSS) model for interacting many-anyon gases in the bosonic limit, revealing a spectrum of energies well described by exact nonlinear Landau levels (NLLs).

Model Derivation and Functional Theory

The study considers a many-particle gas of interacting abelian anyons confined in a trapping potential, modeling them as bosons with magnetic flux attachment where the total magnetic flux is small compared to the number of particles. By means of a Hartree–Jastrow ansatz, they derive a two-parameter Chern–Simons–Schrödinger energy functional which extends the well-known Gross–Pitaevskii / nonlinear Schr¨odinger density functional theory to the magnetic (anyonic) self-interaction. One parameter determines the total number of self-generated magnetic flux units in the system, and the other the effective strength of spin-orbit self-interaction. This latter interaction can be either attractive/focusing or repulsive/defocusing, and depends both on the intrinsic spin-orbit interaction and the relative length scale of the flux profile of the anyons.

Microscopic Hamiltonian and Effective Functional

The microscopic description is given by an Rextended, N-anyon Hamiltonian. If N ≫ 1 while the total flux αN ≃ β is of order 1, then on intermediate length scales R << l << ω − 1/2 the interacting anyon gas is well described by the following Chern–Simons–Schr¨odinger (CSS) energy functional. This functional involves a one-body collective state ϕ ∈ L2(R2) and a self-consistently generated magnetic field Bϱ = curl Aϱ = 2πϱ where the density distribution is ϱ = ϕ2. The ground-state energy of this functional is denoted by E(β, γ, V).

Stability and Critical Coupling

Stability of the system holds if and only if E(β, γ, V) > −∞. NLLs exist at β ∈ 2Z ≠ 0 with the critical coupling γ = −2πβ. The critical attractive coupling for stability is denoted by γ∗(β), which is exactly the smallest (by infimum) ratio between the magnetic-kinetic self-energy Eβ,0,0[ϕ] and the scalar/electrostatic self-interaction Rϕ4 among all normalized and regular ϕ. For any (β, γ) s.t. γ ≥ −γ∗(β), we have E(β, γ, V) > −∞ (stability), while if γ < −γ∗(β) then E(β, γ, V) = −∞ (instability).

Exact Nonlinear Landau Levels (NLLs)

The analysis reveals a spectrum of energies/states that are well describable using a family of exact NLLs. The general formula for NLL states is parametrized by the space of solutions (P, Q) of the inverse Wronskian problem W(P, Q):= P′Q − P Q′ = f. For a fixed n = β/2, the 4n-dimensional manifold of NLL states foliates over the degree d = deg f, with n − 1 ≤ d ≤ 2(n − 1). Specific solutions like the “vortex ring” state are described by P(z) = (z - w)n and Q(z) = 1.

Analytical Energies of Trial States

The paper calculates the exact energy E[ϕ] for certain trial CSS states in the harmonic trap V(z) = Vh(z) = z2. For a radially symmetric, Gaussian state, the total energy at the optimal scaling parameter λ∗ is given by Eβ,γ,Vh[ϕλ∗] = 2r1 + β2 log 4/3 + γ2/π. The analysis of excited harmonic oscillator eigenstates yields minimal energies Eβ,γ,Vh[ϕ±λ±] = 2r4 + γ2/π + 7/27 β2 ± β. Finally, for a symmetric NLL “vortex ring” state of degree n − 1 > 0 centered at w = 0, the minimal total energy becomes Eβ,γ,Vh[ϕλ∗] = r23 s(β − 2n)2 / (2 - 1/n) + (2πβ + γ)n + 1/π × r1 − 1/n π/n sin(π/n).

Conclusion and Outlook

The authors conclude that important keys to understanding the complex interpolation between the almostbosonic and almost-fermionic anyon spectra will come via the NLLs and their possible extension to more general functions. Exploring the two-parameter energy landscape may uncover fruitful and highly nontrivial relations between the CSS, GP, TF, and other models. The study of exact NLLs is mainly limited by the increasing complexity of the inverse Wronskian problem, however they may extend relatively easily to some symmetryreduced cases for β ≥ 10. Ongoing work also includes external magnetic fields and the Abrikosov problem for self-generated vortex lattices. The authors believe that in order to obtain even lower energies, to eventually match the ideal Fermi energy at α → 1, we must relax these assumptions. The study of systems of CSS functions and Liouville–Toda-type equations will be addressed in future work. Some new insight into the fermionic side of the problem, via fermionic Hartree and magnetic TF theory, has been recently discussed in [78] (see also [79]).

How it works

The paper employs several mathematical tools to bridge the gap between different physical descriptions. First, it models the anyon gas using a microscopic Hamiltonian and then approximates the interacting gas using an effective CSS energy functional. This functional is derived from a Hartree–Jastrow ansatz for the microscopic Hamiltonian. The resulting energy functional is then minimized with respect to an optimal length scale λ∗, leading to minimal energies Eλ∗ = 2p K1U1.

Key Findings Enumerated

The paper highlights several key findings regarding the spectrum and stability of the system:

1. The existence of a two-parameter Chern–Simons–Schrödinger energy functional that extends GP theory to magnetic self-interaction.

2. The identification of NLLs at β ∈ 2Z ≠ 0 with the critical coupling γ = −2πβ.

3. The derivation of an effective CSS model with precise relations for the physical parameters, its stability and supersymmetric properties.

4. The analysis of the emergent spectrum, numerically and analytically studied with the help of Jackiw–Pi vortex solitons or “nonlinear Landau levels” (NLLs).

5. The identification of specific NLL states through solutions to the inverse Wronskian problem, which yield exact analytical energies for these states.

The paper demonstrates that the minimal energy E(β, γ, Vh) is achieved at an optimal scaling parameter λ∗ = (U1/K1)1/4. This optimization leads to a sequence of states whose energies align well with a sequence of exact nonlinear Landau levels describing Jackiw–Pi self-dual solitons. The analysis also shows that the negative mode (counterrotating vortex) gives the lowest energy for β > 0 and is shown as the dashed purple curve in the plots.

Numerical Validation

The numerical data, computed in Matlab and optimized using a steepest descent algorithm, provides approximate energies of (2). These numerical results are shown in various tables (Table I through Table XII) and plots (Fig. 4). The convergence to the desired precision is typically achieved after 50-200 iterations of the algorithm. The resulting energies are shown in the plots using solid markers (colors and shapes indicate numbers of visible vortices).

Summary of Key Results

The research provides a quantitative description of the low-energy collective state, showing that for increasing flux, counter-rotating vortices are formed, enhancing the stability of the gas against collapse. The study clarifies relations between various different anyon models and exemplifies a novel supersymmetry-breaking phenomenon. The findings suggest that understanding the complex interpolation between almostbosonic and almost-fermionic anyon spectra will come via the NLLs and their possible extension to more general functions. The two-parameter energy landscape may uncover fruitful and highly nontrivial relations between the CSS, GP, TF, and other models.

1. Improve quantum state modeling by implementing a two-parameter Chern–Simons–Schrödinger energy functional derived from a Hartree–Jastrow ansatz, allowing for explicit control over the total number of self-generated magnetic flux units in the system and the effective strength of spin-orbit self-interaction.

2. Enable analytical calculation of ground and excited states by utilizing the sequence of exact nonlinear Landau levels describing Jackiw–Pi self-dual solitons, which are shown to align with results from the effective CSS model derived for the interacting manyanyon gas.

3. Enhance simulation accuracy by employing numerical methods based on an approximate minimizer for E(β = 10, γ = 20π, Vh) and a steepest descent algorithm adopted from [62] to find the global minimum energy of the functional, specifically addressing issues of geometrically and topologically induced metastability.

4. Develop a predictive tool for system stability by defining a criterion: "stability holds if and only if E(β, γ, V) > −∞, allowing researchers to map out regions where the gas exhibits collapse or stability based on the coupling parameter γ relative to the critical value γ∗(β)."

5. Create a spectral analysis engine capable of enumerating NLL solutions by solving the inverse Wronskian problem for specific polynomial pairs (P, Q) derived from a given flux profile, thereby predicting the number of vortices in an NLL state using the constraint: the number of vortices for a given NLL state ϕP,Q (counting multiplicities) is precisely deg W(P, Q).

6. Implement an energy minimization routine that can optimize trial states by rescaling according to the optimal parameter λ∗:= (U1/K1) 1/4, ensuring the resulting energy Eλ∗ achieves its minimum, as shown in equation (35).

7. Integrate a library of analytical solutions for specific NLL manifolds, such as the vortex ring state described by ϕ(z) = r n π z¯n−1 z2n + 1, allowing the AI to rapidly retrieve or interpolate energy values corresponding to known exact states in the stable regime.