Mobility Crossover in Two-Dimensional Berry Crystals

A Berry crystal is a random superposition of N plane waves of equal amplitude and fixed wavevector magnitude, propagating in different directions. Using numerical simulations of wavepacket dynamics, spectral analysis based on autocorrelation functions, and scaling of Inverse Participation Ratio, the nature of eigenstates across the energy spectrum of a two-dimensional Berry crystal is characterized. It exhibits Anderson localization and critical (extended but non-ergodic) states, reminiscent of quasicrystals, which sit in the middle ground between periodic and disordered systems and can host critical states. However, in contrast to quasicrystals that display sharp mobility edges separating extended and localized phases, the Berry crystal exhibits an extended regimes of critical states. We name this a “mobility crossover”. At weak potential strength, low-energy states are extended while higher-energy states near the backscattering momentum are critical. As the potential strength increases to become comparable with the recoil energy, these critical states evolve into localized states, yielding a transition from extended non-ergodic to localized behavior near the backscattering momentum. An estimate for the boundaries of ergodic extended regimes is given by the mapping onto an effective Anderson model on a Bethe Lattice. The results shed light on the relation between backscattering and Anderson localization in continuous two-dimensional aperiodic systems.

💡 Research Summary

**
In this work the authors introduce and thoroughly investigate a novel two‑dimensional aperiodic potential they call a “Berry crystal.” A Berry crystal is constructed by superposing N plane waves of identical amplitude and fixed wave‑vector magnitude q, each propagating in a random direction θj and carrying a random phase φj. The resulting potential U(r)=β/√(2N)∑jcos(qj·r+φj) possesses both short‑range randomness and long‑range power‑law correlations (in the limit N→∞ the autocorrelation decays as a power of distance). Because the constituent plane waves have a fixed magnitude, the potential supports Bragg‑type scattering while lacking translational symmetry.

The authors study single‑particle dynamics in this potential using three complementary numerical approaches:

1. Wave‑packet propagation – Gaussian minimum‑uncertainty packets (σx=σy=2) are evolved on a 2‑D grid (L=320) with an absorbing boundary. A third‑order split‑operator scheme (e−iHτ=e−iUτ/2 e−iKτ e−iUτ/2) is used for a total dimensionless time T=1.1×10⁴. The autocorrelation ⟨ψ(0)|ψ(t)⟩ is Fourier‑transformed to obtain a weighted density of states ρ(ε). By varying the integration window T′ the authors show that low‑energy peaks are insensitive to T′ (continuous spectrum), whereas high‑energy peaks sharpen but grow sub‑linearly, indicating a singular‑continuous component.

2. Exact diagonalization + Lanczos – To access eigenstates directly, the Hamiltonian is diagonalized with Dirichlet boundaries. Edge‑localized states are removed using the “binary‑map” filter of Zhu et al. The inverse participation ratio (IPR = ∫|ψ(r)|⁴ dr) is computed for system sizes L=200–600 and fitted to IPR∝L⁻ᵞ. The scaling exponent γ serves as an order parameter: γ≈2 for fully ergodic extended states, γ≈0 for fully localized states, and 0<γ<2 for multifractal critical (non‑ergodic extended) states.

3. Momentum‑space mapping – In the plane‑wave basis the Hamiltonian reads ⟨k|H|k′⟩ = δ_{k,k′}ε_k + (β/2)∑j(e^{iφj}δ_{k′,k+qj}+e^{-iφj}δ_{k′,k−qj}). The potential therefore induces hopping in momentum space by vectors ±qj. For a given energy shell |k|≈k₀ the set of resonant momenta forms a graph where each node is connected to four nearest neighbours (two forward and two backward scattering processes). Because the random set {qj} prevents closed loops, the graph is a Cayley tree with branching number K≈3, i.e. a Bethe lattice. The onsite energies are set by the Bragg angle θ₀ and small angular detuning δθ, while the hopping amplitude is t=β/√(2N) at first order (or t_eff≈β²/(8N) at second order when k₀<q/2). Thus the low‑energy sector of the Berry crystal maps onto an Anderson model on a Bethe lattice, a system whose localization properties are analytically tractable.

Key numerical findings

• For weak potential strength (β≤0.7) the low‑energy part of the spectrum is ergodic (γ≈2). As energy approaches the back‑scattering momentum k=½q, γ drops to ≈1.5, signalling a transition to critical, non‑ergodic extended states. At still higher energies γ rises again, indicating a re‑entrant extended regime.
• When β≈0.8 (comparable to the recoil energy E₀=ℏ²q²/8m) the critical region widens and eventually a localized regime (γ≈0) appears around the back‑scattering energy. Above this localized band, higher‑energy states become critical again.
• The scaling exponent γ is essentially independent of the number of plane‑wave components N (tested for N=20–400), confirming that the observed phases are intrinsic to the Berry crystal rather than finite‑N artifacts.
• By scanning the (E,β) plane and extracting the fractal dimension D₂ from IPR scaling, the authors construct a phase diagram: an extended (D₂≈2) region at low β, a broad critical band (1<D₂<2) that expands with β, and a localized island (D₂≈0) that emerges near the back‑scattering energy for β≳0.8.

Conceptual significance
The Berry crystal exhibits a “mobility crossover” rather than a sharp mobility edge. In conventional disordered 2‑D systems all states localize, while quasicrystals can display well‑defined mobility edges separating extended and localized phases. The Berry crystal, however, hosts an extended regime of critical (non‑ergodic) states whose extent is controlled by the disorder strength β. The crossover is tied to back‑scattering: states with momentum close to k=½q experience near‑resonant scattering that generates the Cayley‑tree connectivity in momentum space, leading to multifractal wavefunctions. As β increases, the effective disorder on the Bethe lattice pushes the system across the known Anderson transition on the tree, turning critical states into truly localized ones.

Experimental relevance – The potential can be realized in a single‑mode laser cavity with chaotic or rough walls, where the intracavity field naturally forms a random superposition of plane waves. Similar aperiodic potentials could also be engineered with ultracold atoms in optical speckle patterns or with photonic structures that combine Bragg scattering and disorder.

Overall contribution – By combining time‑dependent wave‑packet simulations, finite‑size scaling of IPR, and an analytical mapping to an Anderson model on a Bethe lattice, the paper provides a comprehensive picture of how back‑scattering and long‑range correlations give rise to a novel mobility crossover in a continuous 2‑D aperiodic system. It bridges the gap between Anderson localization in random media and critical phenomena in quasicrystals, opening new avenues for exploring non‑ergodic extended phases in both matter‑wave and photonic platforms.