REVIEW 3 major objections 5 minor 39 references
Universal self-similar evolution of two-dimensional acoustic turbulence in Bose-Einstein condensates
T0 review · 3 major / 5 minor · reviewed 2026-08-02 · deepseek-v4-flash
Pith's one-line read A single universal constant β≈8.9 governs how 2D acoustic turbulence in Bose-Einstein condensates evolves toward small scales, with the energy spectrum collapsing onto a master curve independent of initial conditions and damping.
desk verdict New front law for 2D acoustic turbulence is well supported in the undamped case; the dissipation-independence claim is under-tested. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the 2D acoustic wave kinetic equation (WKE) for Bogoliubov waves, whose collision integral is taken from weak wave turbulence theory in the small-healing-length limit, with an added linear damping term. Substituting the self-similar ansatz ek(t)=E0 f(t) Φ(k/k*(t)) into the WKE and requiring time dependence to cancel yields a nonlinear eigenvalue equation for Φ and a differential equation for the front k*(t); the dimensionless constant β appears as the eigenvalue. A differential approximation to the collision integral provides the exponential large-η tail Φ~exp(-βη/D), while the small-η asymptotic Φ~η^{-1} follows from the KZ spectrum.
What would settle it
Measure k*(t) in a 2D BEC experiment or in a high-resolution GPE simulation starting from a different initial spectral shape (e.g., a sharp ring at intermediate wavenumbers instead of large-scale equipartition). If ln(k*/k0) plotted against sqrt(t) does not become a straight line with slope sqrt(2β/τ), or if the extracted β differs from 8.9 by more than the numerical uncertainty, the universality claim fails. A direct GPE simulation with nonzero linear damping that shows a k*(t) deviating from the predicted factor sqrt(τ_D/τ)(1-e^{-t/τ_D}) would also falsify the damping-independence claim.
Extended reading notes
Core claim
In the nonstationary acoustic regime, the energy spectrum obeys ek(t)=E0 f(t) Φ(k/k*(t)) with Φ~η^{-1} in the wake of the front and an exponential tail at large η. The front itself obeys (1/k*) dk*/dt ln(k*/k0) = (β/τ) exp(-t/τ_D), so that k*(t) grows as exp(sqrt(2βt/τ)) when undamped and as exp(sqrt(2βτ_D/τ)(1-e^{-t/τ_D})) with damping. The dimensionless constant β is a universal eigenvalue emerging from the self-similar equation, taking the value β≈8.9, identical for GPE and WKE simulations, with and without dissipation. This places 2D acoustic turbulence at the crossover between first- and second-kind self-similarity, where the energy integral for the Kolmogorov-Zakharov spectrum is only
Load-bearing premise
The universal constant β≈8.9 is derived from the 2D acoustic wave kinetic equation, which is rigorously justified only in the weakly nonlinear, no-forcing, no-dissipation limit; if that collision integral is not quantitatively faithful to the Gross-Pitaevskii dynamics in the acoustic regime, the universality of β may not transfer to real condensates.
Editorial extensions
If this is right
- The front propagation law k*(t)=k0 exp(sqrt(2βt/τ)) provides a quantitative, parameter-free prediction for how rapidly small-scale density structures develop in a 2D atomic BEC, testable in current experiments.
- Because damping does not change the self-similar shape Φ and only modifies the front law through the factor sqrt(τ_D/τ)(1-e^{-t/τ_D}), polariton condensates with their finite lifetime should still exhibit the same universal spectrum.
- The wake spectrum is the Kolmogorov-Zakharov spectrum k^{-1}, now derived as the unique asymptotic solution of the similarity equation, resolving the infrared divergence through the logarithmic cutoff k0.
- The system exemplifies a borderline infinite-capacity case where the energy integral is only logarithmically divergent, explaining why dimensional analysis alone fails and β must be computed as an eigenvalue.
- Measuring k*(t) in an experiment or high-resolution simulation directly yields β, allowing the universality claim to be checked against the predicted master curve.
Reading between the lines
- If the universality extends beyond the weakly-nonlinear regime tested here, analogous self-similar front laws with a logarithmic-capacity KZ spectrum may appear in other systems, such as 3D acoustic turbulence or weak gravitational-wave turbulence, though the value of β would need to be computed separately for each.
- The predicted exponential tail Φ~exp(-βη/D) at large k is independently measurable; its slope in a log-linear plot provides a cross-check of β against the value extracted from front propagation.
- The damping-independence claim is tested only at the WKE level; a direct GPE simulation with linear damping would strengthen the case that polariton condensates exhibit the same universal constant β≈8.9.
- A natural testable extension is to perturb the condensate with initial conditions of different spectral shape (not just large-scale equipartition): universality predicts the same Φ and β should emerge after a transient, offering a clean experimental signature.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the nonstationary evolution of two-dimensional acoustic turbulence in Bose-Einstein condensates, using numerical simulations of the Gross-Pitaevskii equation (GPE) and the corresponding wave kinetic equation (WKE). It proposes that the energy spectrum exhibits self-similar evolution of a mixed first/second-kind type, with a universal front law k*(t) = k0 exp[sqrt(2βt/τ)] in the undamped case and k*(t) = k0 exp[sqrt(2βτ_D/τ (1 - e^{-t/τ_D}))] in the damped case, where β is a dimensionless universal constant measured as β ≈ 8.9. The paper further claims that this universality is independent of dissipation and represents a new type of non-thermal fixed point.
Significance. If correct, the identification of a new, genuinely two-dimensional acoustic-turbulence self-similar regime with a universal constant β would be a valuable contribution to wave turbulence and non-thermal fixed points. The main strengths are: (i) the direct comparison between GPE and WKE simulations, which supports the undamped self-similar collapse; (ii) the nontrivial functional form of the front law, exp(sqrt(t)), which is much more specific than a generic power law; and (iii) the transparent presentation of the self-similar ansatz and the differential-approximation analysis giving the exponential tail. However, the quantitative universality of β and the dissipation-independence claim rest on weaker evidence than the undamped core result.
major comments (3)
- [Eq. (10) and Fig. 3b] β ≈ 8.9 is obtained by fitting Eq. (10) to the same k*(t) time series that is used to construct the self-similar collapse in Fig. 3a. This is a fit, not a parameter-free prediction: the eigenvalue problem Eq. (8) is not solved independently for β. The paper should either solve Eq. (8) numerically to provide a genuine prediction, or at least report the fitting procedure, the statistical uncertainty of β, and a test of robustness (e.g., fitting on different time windows, different initial conditions, and checking that the fit residual is consistent). Without this, the 'universal constant' claim is not quantitatively grounded.
- [Eq. (3) and damped front law] The dissipation independence of β is established only at the WKE level, using an ad hoc linear damping term −γ e_k. The paper itself states that the WKE derivation is rigorous only in the absence of forcing and dissipation, and that 2D acoustic WWTT has serious mathematical issues. Since the headline claim includes 'regardless of dissipative effect', a damped GPE simulation is needed to verify the damped front law and β. Moreover, the damping in the GPE, −iγΨ in Eq. (1), gives an energy decay e^{-2γt}, not the e^{-γt} used in Eq. (4); if the damped front law were applied to the GPE, τ_D would be (2γ)^{-1}. This factor-of-two discrepancy is not discussed and is directly relevant to the polariton context.
- [SM Eq. (19) and Eq. (7)] The derivation of the front law Eq. (7) neglects the 1/k_* term in Eq. (19) of the SM, relying on k_* >> k0. In the GPE simulation with a 1024ξ box, the dynamic range of k_* is limited, and the collapse in Fig. 3a mixes runs with very different k_*/k0 ranges. The paper should verify that the fitted β is insensitive to the lower cut-off k0 and that the neglected term is indeed negligible over the fitted time window. This is important because the untrapped GPE has the shortest inertial range and yet is used to support β≈8.9.
minor comments (5)
- [Introduction] Typo: 'Perhaps, one the most simple' should read 'one of the most simple'.
- [Conclusion] Typo: 'abscence' should be 'absence'.
- [Eq. (2)] The definition of ek uses the notation ∫_{|k|=k} k^2 |ψ̂_k|^2 dk, which is a line integral over a circle; this is clear in context but could be stated explicitly for the non-specialist.
- [Fig. 3b inset] The caption refers to a 'double exponential regime' but the horizontal axis is sqrt(τ_D/τ)(1 − e^{−t/τ_D}). This is a stretched-exponential-like variable, not a double exponential; the phrasing should be adjusted.
- [Code availability] The paper names the codes FROST and WavKinS but does not state a data/code availability statement. Since the central result rests on numerics, a brief statement on availability would be helpful.
Circularity Check
No significant circularity: the self-similar front law is derived from the WKE with β honestly presented as a numerically fitted eigenvalue; the main fragility is an unvalidated damped-GPE extrapolation, not a circular reduction.
full rationale
The derivation chain is largely self-contained: the self-similar ansatz (6) is substituted into the WKE (3), yielding the front equation (7) and the eigenvalue equation (8) by algebra; the functional forms E(t)=E0 e^{-t/τ_D}, f(t)=e^{-t/τ_D}/[k_* ln(k_*/k_0)], and the exponential tail (9) all follow from the paper's own equations. The collapse in Fig. 3a and the consistency of Φ(η) across GPE and WKE runs provide an independent test of the self-similar functional form that does not use β. The constant β≈8.9 is not derived from first principles; the paper states explicitly: 'The numerical evolution of k*(t) is then fitted to Eq. (10), leading to the universal value β≈8.9.' Thus the close agreement in Fig. 3b for that constant is in-sample fitting, but it is transparently reported as a numerical determination rather than a parameter-free prediction, so it does not reduce the central claim by construction. The WKE itself is taken from the authors' prior work (refs. [29-31]), and the codes FROST and WavKinS are also author tools; however, the undamped predictions are checked directly against GPE simulations, so these self-citations are corroborated rather than purely load-bearing. The paper's own caveats—'the derivation of the WKE is rigorous only in the absence of forcing and dissipation' and 'WWTT presents serious mathematical issues for acoustic waves, notably in 2D'—identify a real fragility: the damping independence of β is tested only at the WKE level, not in damped GPE runs, and the GPE damping term −iγΨ gives energy decay e^{-2γt} while the WKE uses e^{-γt}. These are validation gaps/correctness risks, not evidence that an output is equivalent to an input by definition. No circular step meeting the quoted-evidence standard was found.
Assumptions & free parameters
free parameters (1)
- β (universal front constant) =
≈8.9
assumptions (6)
- domain assumption The 2D acoustic WKE (Eq. 3) with small-ξ prefactor describes GPE dynamics in the acoustic regime.
- domain assumption Self-similar ansatz e_k(t) = E0 f(t) Φ(η), η = k/k*(t).
- domain assumption Asymptotic Φ(η) ~ η^{-1} as η→0.
- domain assumption Limit k* >> k0 for dropping the 1/k* term in the self-similar reduction of the WKE.
- ad hoc to paper Damping enters the WKE as -γ e_k.
- domain assumption At large η the collision integral reduces to the differential diffusion operator D k^3 ∂_k(k^{-2} ∂_k e_k).
Cite this review
Pith. "Pith review of Universal self-similar evolution of two-dimensional acoustic turbulence in Bose-Einstein condensates." pith.science (2026). https://pith.science/paper/P5UWQAT4
@misc{pith2026260706062,
author = {Pith},
title = {Pith review of: Universal self-similar evolution of two-dimensional acoustic turbulence in Bose-Einstein condensates},
year = {2026},
howpublished = {\url{https://pith.science/paper/P5UWQAT4}},
note = {Machine review of arXiv:2607.06062}
}
abstract
When driven out of equilibrium, a Bose-Einstein condensate develops nonlinearly interacting density waves that trigger a turbulent cascade, transferring energy toward small scales. In this Letter, we investigate the nonstationary evolution of solutions to the two-dimensional Gross-Pitaevskii equation (GPE). Through numerical simulations of both the GPE and the corresponding Wave Kinetic Equation (WKE), we identify self-similar solutions relevant to turbulence in atomic and polariton Bose-Einstein Condensates. These solutions correspond to a new type of non-thermal fixed point and exhibit characteristics of both first and second kind self-similarity. In particular, we show that the dynamics of the propagating front is universal, governed by a dimensionless universal constant $\beta$, which we determine numerically.
Figures
Reference graph
Works this paper leans on
-
[1]
G. I. Taylor, The formation of a blast wave by a very intense explosion i. theoretical discussion, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences201, 159 (1950)
1950
-
[2]
G. I. Taylor, The formation of a blast wave by a very intense explosion.-ii. the atomic explosion of 1945, Pro- ceedings of the Royal Society of London. Series A. Math- ematical and Physical Sciences201, 175 (1950)
1945
-
[3]
Y. B. Zel’Dovich and Y. P. Raizer,Physics of shock waves and high-temperature hydrodynamic phenomena(Courier Corporation, 2002)
2002
-
[4]
Schmied, A
C.-M. Schmied, A. N. Mikheev, and T. Gasenzer, Non- thermal fixed points: Universal dynamics far from equi- librium, International Journal of Modern Physics A34, 1941006 (2019)
2019
-
[5]
Connaughton and S
C. Connaughton and S. Nazarenko, Warm cascades and anomalous scaling in a diffusion model of turbulence, Physical review letters92, 044501 (2004)
2004
-
[6]
Grebenev, S
V. Grebenev, S. Nazarenko, S. Medvedev, I. Schwab, and Y. A. Chirkunov, Self-similar solution in the leith model of turbulence: anomalous power law and asymp- totic analysis, Journal of Physics A: Mathematical and Theoretical47, 025501 (2014)
2014
-
[7]
Thalabard, S
S. Thalabard, S. Nazarenko, S. Galtier, and S. Medvedev, Anomalous spectral laws in differential models of turbu- lence, Journal of Physics A: Mathematical and Theoret- ical48, 285501 (2015)
2015
-
[8]
W. J. Bos, C. Connaughton, and F. Godeferd, Developing homogeneous isotropic turbulence, Physica D: Nonlinear Phenomena241, 232 (2012)
2012
Show all 39 references
-
[9]
Costa, A
G. Costa, A. Barral, A. Lopez, Q. Pikeroen, and B. Dubrulle, Behind the mirror: The hidden dissipative singular solutions of ideal reversible fluids on log-lattices, Physical Review Fluids10, 114603 (2025)
2025
-
[10]
C. S. Campolina and A. A. Mailybaev, Chaotic blowup in the 3d incompressible euler equations on a logarithmic lattice, Physical review letters121, 064501 (2018)
2018
-
[11]
Pikeroen, A
Q. Pikeroen, A. Barral, G. Costa, C. Campolina, A. Mai- lybaev, and B. Dubrulle, Tracking complex singularities of fluids on log-lattices, Nonlinearity37, 115003 (2024). 6
2024
-
[12]
D. V. Semikoz and I. I. Tkachev, Kinetics of bose con- densation, Physical Review Letters74, 3093 (1995)
1995
-
[13]
Y. Zhu, B. Semisalov, G. Krstulovic, and S. Nazarenko, Self-similar evolution of wave turbulence in gross- pitaevskii system, Physical Review E108, 064207 (2023)
2023
-
[14]
V. E. Zakharov, Weak turbulence in media with a decay spectrum, Journal of Applied Mechanics and Technical Physics6, 22 (1965)
1965
-
[15]
Nazarenko,Wave turbulence, Vol
S. Nazarenko,Wave turbulence, Vol. 825 (Springer Sci- ence & Business Media, 2011)
2011
-
[16]
Galtier and S
S. Galtier and S. V. Nazarenko, Turbulence of weak grav- itational waves in the early universe, Physical review let- ters119, 221101 (2017)
2017
-
[17]
D¨ uring, C
G. D¨ uring, C. Josserand, and S. Rica, Weak turbulence for a vibrating plate: Can one hear a kolmogorov spec- trum?, Physical review letters97, 025503 (2006)
2006
-
[18]
Galtier, S
S. Galtier, S. Nazarenko, A. C. Newell, and A. Pouquet, A weak turbulence theory for incompressible magnetohy- drodynamics, Journal of plasma physics63, 447 (2000)
2000
-
[19]
Caillol and V
P. Caillol and V. Zeitlin, Kinetic equations and station- ary energy spectra of weakly nonlinear internal grav- ity waves, Dynamics of atmospheres and oceans32, 81 (2000)
2000
-
[20]
Galtier, Weak inertial-wave turbulence theory, Physi- cal Review E68, 015301 (2003)
S. Galtier, Weak inertial-wave turbulence theory, Physi- cal Review E68, 015301 (2003)
2003
-
[21]
Y. Zhu, B. Semisalov, G. Krstulovic, and S. Nazarenko, Direct and inverse cascades in turbulent bose-einstein condensates, Physical Review Letters130, 133001 (2023)
2023
-
[22]
Ga lka, P
M. Ga lka, P. Christodoulou, M. Gazo, A. Karailiev, N. Dogra, J. Schmitt, and Z. Hadzibabic, Emergence of isotropy and dynamic scaling in 2d wave turbulence in a homogeneous bose gas, Physical Review Letters129, 190402 (2022)
2022
-
[23]
Navon, A
N. Navon, A. L. Gaunt, R. P. Smith, and Z. Hadzibabic, Emergence of a turbulent cascade in a quantum gas, Na- ture539, 72 (2016)
2016
-
[24]
K. G. Lagoudakis, M. Wouters, M. Richard, A. Baas, I. Carusotto, R. Andr´ e, L. S. Dang, and B. Deveaud- Pl´ edran, Quantized vortices in an exciton–polariton con- densate, Nature physics4, 706 (2008)
2008
-
[25]
Panico, P
R. Panico, P. Comaron, M. Matuszewski, A. Lanotte, D. Trypogeorgos, G. Gigli, M. D. Giorgi, V. Ardizzone, D. Sanvitto, and D. Ballarini, Onset of vortex clustering and inverse energy cascade in dissipative quantum fluids, Nature Photonics17, 451 (2023)
2023
-
[26]
Byrnes, N
T. Byrnes, N. Y. Kim, and Y. Yamamoto, Exciton– polariton condensates, Nature Physics10, 803 (2014)
2014
-
[27]
Carusotto and C
I. Carusotto and C. Ciuti, Quantum fluids of light, Re- views of Modern Physics85, 299 (2013)
2013
-
[28]
Krstulovic,A theoretical description of vortex dy- namics in superfluids
G. Krstulovic,A theoretical description of vortex dy- namics in superfluids. Kelvin waves, reconnections and particle-vortex interaction, Habilitation ` a diriger des recherches, Universite Cˆ ote d’Azur (2020)
2020
-
[29]
V. E. Zakharov and R. Z. Sagdeev, Spectrum of acoustic turbulence, inSoviet Physics Doklady, Vol. 15 (1970) p. 439
1970
-
[30]
Costa, G
G. Costa, G. Krstulovic, and S. Nazarenko, Stability of stationary solutions in acoustic wave turbulence, arXiv preprint arXiv:2508.09799 (2025)
2025
-
[31]
Griffin, G
A. Griffin, G. Krstulovic, V. S. L’vov, and S. Nazarenko, Energy spectrum of two-dimensional acoustic turbulence, Physical review letters128, 224501 (2022)
2022
-
[32]
V. E. Zakharov, V. S. L’vov, and G. Falkovich, Kolmogorov spectra of turbulence I: Wave turbulence (Springer Science & Business Media, 2012)
2012
-
[33]
Y. Zhu, G. Krstulovic, and S. Nazarenko, Turbulence and far-from-equilibrium equation of state of bogoli- ubov waves in bose-einstein condensates, arXiv preprint arXiv:2408.15163 (2024)
2024 arXiv
-
[34]
Krstulovic and V
G. Krstulovic and V. Labarre, Wavkins. jl: an efficient and modular julia software for solving wave kinetic equa- tions, arXiv preprint arXiv:2504.00252 (2025)
2025 arXiv
-
[35]
We shall also notice that this KZ solution also diverges in the infrared (IR), which is inconsistent with the finite energyE 0 of the ini- tial condition
for known scaling relations). We shall also notice that this KZ solution also diverges in the infrared (IR), which is inconsistent with the finite energyE 0 of the ini- tial condition. Nevertheless, the evolution proceeds only towards the UV, and it is then natural to introduc...
-
[36]
Chantesana, A
I. Chantesana, A. P. Orioli, and T. Gasenzer, Kinetic theory of nonthermal fixed points in a bose gas, Physical Review A99, 043620 (2019)
2019
-
[37]
Hasselmann, K
S. Hasselmann, K. Hasselmann, J. H. Allender, and T. P. Barnett, Computations and parameterizations of the nonlinear energy transfer in a gravity-wave spectrum. part ii: Parameterizations of the nonlinear energy trans- fer for application in wave models, Journal of Physical Oc...
1985
-
[38]
Nazarenko, Sandpile behaviour in discrete water-wave turbulence, Journal of Statistical Mechanics: Theory and Experiment2006, L02002 (2006)
S. Nazarenko, Sandpile behaviour in discrete water-wave turbulence, Journal of Statistical Mechanics: Theory and Experiment2006, L02002 (2006)
2006
-
[39]
Nazarenko, Differential approximation for kelvin wave turbulence, Journal of Experimental and Theoretical Physics Letters83, 198 (2006)
S. Nazarenko, Differential approximation for kelvin wave turbulence, Journal of Experimental and Theoretical Physics Letters83, 198 (2006). 7 Supplemental material : Universal self-similar evolution of two-dimensional acoustic turbulence in Bose-Einstein condensates. NUMERICAL...
2006
Reviewed August 2, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.