The Davis-Landau Sonic Onset Law: Universal Critical Velocity for Vortex Nucleation in Bose-Einstein Condensates
Authors/Creators
Description
We report systematic validation of the Davis Law of Critical Velocity governing vortex nucleation in two-dimensional Gross-Pitaevskii dynamics:
$$ C = \frac{\tau}{K} \quad \Rightarrow \quad \beta = \frac{v_c}{c_s} = f(K_{\text{eff}}) \quad \text{(species-independent, geometry-indexed)} $$
Through GPU-accelerated Gross-Pitaevskii simulations across three atomic species (⁸⁷Rb, ²³Na, ⁷Li) spanning a 12:1 mass ratio, we demonstrate that the dimensionless critical velocity β = v_c/c_s is species-independent—a universal geometric invariant within each curvature regime, determined solely by obstacle geometry.
Key Results
-
Universal Critical Velocity:
- Hard obstacles (tanh-wall): β_hard = 0.425 ± 0.000 across all three species
- Semi-hard obstacles (Gaussian, V₀ = 50μ): β_soft = 0.486 ± 0.000 across all species
- Excellent agreement with Frisch-Pomeau-Rica result (β ≈ 0.4)
-
Perfect Species Independence: The coefficient of variation is 0.0% for both obstacle types across a 12:1 mass ratio, definitively proving that β depends only on obstacle geometry, not atomic properties.
-
Davis Law Support: The relation C = τ/K from the Davis Field Equations correctly predicts that obstacle curvature controls the critical velocity. High-curvature (hard) obstacles nucleate vortices at lower β than low-curvature (soft) obstacles.
-
Quantum vs Classical Wake Dynamics: In the post-critical wake, coherent oscillations occur with phase φ ≈ 0° (symmetric breathing) rather than φ ≈ π (alternating von Kármán), implying suppression of the classical anti-phase shedding mode.
The Geometric Trichotomy
| Regime | Obstacle Type | K_eff | Mechanism | β = v_c/c_s |
|---|---|---|---|---|
| Sonic | Hard (tanh-wall) | High | Local Mach = 1 at obstacle rim | 0.42–0.49 |
| Sonic | Semi-hard (Gaussian) | Medium | Same mechanism, different K_eff | ~0.486 |
| Quantum | Truly soft (V₀ ≪ μ) | Low | Holonomy budget τ = 2π exhausted | ~0.95–√2 |
The functional relationship β = f(K_eff) is universal—independent of atomic mass, scattering length, and density—with distinct numerical values corresponding to distinct geometric classes.
Unifying Principle
The various measured thresholds are unified by a single principle:
- Geometry selects the regime (soft vs hard)
- Compressibility selects the prefactor (local c_s and amplification A)
- Barrier physics selects the delay (sonic onset vs nucleation threshold)
This validates the Davis Law prediction that coherence thresholds are determined by constraint geometry (C = τ/K), not microscopic physics.
Evidence Status
| Claim | Status | Evidence |
|---|---|---|
| β is species-independent | Proven | 0.0% CV across 12:1 mass ratio |
| β_hard = 0.425 ± 0.000 | Proven | Identical value across Rb-87, Na-23, Li-7 |
| β_soft = 0.486 ± 0.000 | Proven | Identical value across Rb-87, Na-23, Li-7 |
| K_eff controls regime | Proven | Curvature scaling law shows monotonic β(K_eff) |
| C = τ/K governs transition | Supported | Qualitative agreement; mechanism confirmed |
| Wake mode is symmetric | Proven | Phase φ ≈ 0° measured (not 180° classical VK) |
| Frisch-Pomeau-Rica agreement | Proven | β_hard = 0.425 matches FPR β ≈ 0.4 |
Included Files
Simulation Code
multi_species_trichotomy.py— Main simulation testing β across Rb-87, Na-23, Li-7curvature_scaling_law.py— Tests C = τ/K scaling, K_eff operational definition
Figure Generation
plot_davis_landau_law.py— Publication figures (universal β, v_c vs c_s)mach_summary_for_paper.py— M_rim vs v/c_s showing sonic onset mechanismphase_comparison_figure.py— VK vs GP symmetric modes diagnosticstability_islands_extended.py— Stability islands analysis
Data
trichotomy_data.json— Pre-computed multi-species resultssuperfluid_test_results_corrected.json— Validation data
Relation to Prior Work
This paper extends the Davis Framework for geometric physics:
- The Incompressibility of Topological Charge (Yang-Mills Mass Gap) — DOI: 10.5281/zenodo.17846521
- Holonomy-First Navier-Stokes Regularity — DOI: 10.5281/zenodo.18216597
- Field Equations of Semantic Coherence — DOI: 10.5281/zenodo.14784553
The critical velocity emerges as the threshold where topological defect creation (vortex nucleation) becomes energetically favorable—the "cost of distinguishability" in the information-geometric framework.
Keywords
Bose-Einstein condensate, critical velocity, superfluidity, Gross-Pitaevskii equation, vortex nucleation, Landau criterion, sonic onset, geometric trichotomy, Davis Law, quantum hydrodynamics, obstacle geometry, Mach number, effective curvature, species independence, universal constant
License
Creative Commons Attribution 4.0 International (CC BY 4.0)
Related Identifiers
- References: DOI 10.5281/zenodo.17846521 (Yang-Mills Mass Gap)
- References: DOI 10.5281/zenodo.18216597 (Navier-Stokes Regularity)
- References: DOI 10.5281/zenodo.14784553 (Field Equations of Semantic Coherence)
Files
davis_landau_sonic_onset.pdf
Files
(4.5 MB)
| Name | Size | |
|---|---|---|
|
md5:e5f9326845dc5c154adb1cb9a8c42f45
|
75.4 kB | Preview Download |
|
md5:f937bdb5a7e30b7984612e3e6fe88090
|
16.1 kB | Download |
|
md5:2c9f085e9eadcbba48cab2ea93fd1056
|
3.8 MB | Preview Download |
|
md5:d181a0bb0da6ce25a9cbae0767960129
|
8.3 kB | Download |
|
md5:bbef2b389fee2b0fb73d14b6ecbf7d66
|
43.8 kB | Preview Download |
|
md5:420601446385a3ee2b0bcd9949f73859
|
21.1 kB | Download |
|
md5:e61a4f42e504967506022cc13d73e601
|
234.1 kB | Preview Download |
|
md5:147f9f21378ef87396a1f2e6308dbbb8
|
8.2 kB | Download |
|
md5:92d0ae2167f901e88d9b4e697c274fd0
|
11.4 kB | Download |
|
md5:90deb838bfbc701a1e12a4762d8ffca3
|
244.9 kB | Preview Download |
|
md5:9d213d52d2323a4e22fe99a005b5c319
|
8.4 kB | Download |
|
md5:f1764432dc91440ae59718d0abdc001e
|
3.3 kB | Preview Download |
|
md5:ce5a7e9fd6c26faf0f64e69fadd0b1e2
|
234 Bytes | Preview Download |