Universal self-similar evolution of two-dimensional acoustic turbulence in Bose-Einstein condensates

Giorgio Krstulovic; Guillaume Costa; Sergey Nazarenko

Universal constant β≈8.9 governs BEC turbulence cascade

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · glm-5.2

2026-07-08 17:36 UTC pith:P5UWQAT4

load-bearing objection 2D acoustic BEC turbulence exhibits hybrid first/second-kind self-similarity governed by a universal constant β≈8.9, but the eigenvalue extraction methodology has a gap that needs addressing. the 2 major comments →

arxiv 2607.06062 v2 pith:P5UWQAT4 submitted 2026-07-07 cond-mat.quant-gas nlin.CD

Universal self-similar evolution of two-dimensional acoustic turbulence in Bose-Einstein condensates

Guillaume Costa , Sergey Nazarenko , Giorgio Krstulovic This is my paper

classification cond-mat.quant-gas nlin.CD PACS 67.85.Hj47.27.-i05.45.-a

keywords Bose-Einstein condensatewave turbulenceself-similarityGross-Pitaevskii equationwave kinetic equationacoustic wavesKolmogorov-Zakharov spectrumnon-equilibrium dynamics

verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper claims that when a two-dimensional Bose-Einstein condensate is driven out of equilibrium in its acoustic regime, the resulting turbulent energy cascade evolves in a self-similar manner governed by a single universal dimensionless constant β≈8.9. The authors derive a self-similarity equation from the wave kinetic equation (WKE), which describes the statistical evolution of the energy spectrum of weakly nonlinear acoustic waves. The constant β emerges as a nonlinear eigenvalue of this equation, meaning it cannot be determined by dimensional analysis alone — it must be solved for. The front of the energy cascade, which propagates from large to small spatial scales, moves as k*(t) = k_0 exp[√(2βt/τ)] in the conservative case, where τ is a characteristic time set by the system's physical parameters and initial energy. The self-similar spectral shape Φ is shown to be independent of dissipation: whether the system is a lossless atomic BEC or a damped polariton BEC, the same universal profile governs the cascade. The authors verify these predictions through numerical simulations of both the full Gross-Pitaevskii equation (GPE) and the WKE, finding quantitative agreement. The system is identified as a hybrid case: it has an infinite-capacity Kolmogorov-Zakharov spectrum (which normally implies first-kind self-similarity determined by dimensional analysis), yet the scaling exponent is dynamically selected (a hallmark of second-kind self-similarity). The authors attribute this to the borderline nature of 2D acoustic turbulence, where the energy integral diverges only logarithmically.

Core claim

The self-similar evolution of the 2D acoustic BEC energy spectrum is fully characterized by a universal dimensionless constant β≈8.9 and a universal spectral shape Φ. The constant β arises as a nonlinear eigenvalue from the self-similarity equation −β[Φ(η)+ηΦ′(η)] = S_t[Φ], where η = k/k*(t) is the rescaled wavevector and S_t is the wave kinetic collision integral. The front propagates as k*(t) = k_0 exp[√(2βt/τ)] without dissipation, and as k*(t) = k_0 exp[√(2βτ_D/τ (1−e^{−t/τ_D}))] with dissipation, where τ_D is the dissipation timescale. The self-similar form Φ is independent of the damping rate. The spectrum behind the front follows the Kolmogorov-Zakharov scaling e_k ∝ k^{−1}, and ahead

What carries the argument

The wave kinetic equation (WKE) for 2D acoustic waves, the self-similarity ansatz e_k(t) = E_0 f(t) Φ(k/k*(t)), and the nonlinear eigenvalue problem for β

Load-bearing premise

The wave kinetic equation is a quantitatively accurate description of the 2D Gross-Pitaevskii equation dynamics in the acoustic regime, even though the authors note that the weak wave turbulence theory has known mathematical issues for acoustic waves in 2D, and the damping term is added in an ad-hoc manner.

What would settle it

If direct numerical simulations of the GPE with substantially larger domains or higher resolution yield a value of β that differs from 8.9, or if the self-similar collapse of the spectrum fails when tested over a wider range of initial conditions, the universality claim would be undermined.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

If β is truly universal, experiments with 2D atomic BECs or polariton condensates should observe the same front-propagation law regardless of the specific initial perturbation, provided the acoustic regime holds.
The independence of Φ from dissipation means polariton BECs, despite their short particle lifetimes, should exhibit the same self-similar spectral shape as lossless atomic BECs — a testable prediction for optical experiments.
The hybrid first/second-kind classification suggests that systems with logarithmically divergent capacity spectra may generically exhibit dynamically selected exponents, extending beyond BECs to other wave turbulence systems.
The exponential front propagation (double-exponential in k* with dissipation) is a distinctive signature that could be used to identify the acoustic turbulent regime experimentally, distinguishing it from other cascade mechanisms.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The borderline nature of 2D acoustic turbulence — where the KZ spectrum's energy integral diverges only logarithmically — may imply that perturbations to the dispersion relation or dimensionality could push the system toward either pure first-kind or pure second-kind self-similarity, making β a sensitive probe of the effective dimensionality and dispersion.
If the WKE quantitatively fails to capture GPE dynamics at late times (as the authors acknowledge mathematical issues exist), the measured β from GPE simulations could deviate from the WKE-predicted value, which would reveal where the weak turbulence approximation breaks down for acoustic waves in 2D.
The structure of the self-similarity equation suggests that β encodes information about the efficiency of nonlinear energy transfer across scales, and its numerical value could potentially be related to other universal constants in turbulent cascade problems.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Desk Editor's Note

2D acoustic BEC turbulence exhibits hybrid first/second-kind self-similarity governed by a universal constant β≈8.9, but the eigenvalue extraction methodology has a gap that needs addressing.

read the letter

The main thing to know: this paper identifies a new self-similarity class for 2D acoustic wave turbulence in BECs, sitting at the borderline between first and second kind. The front propagation is governed by a universal constant β≈8.9, extracted from both GPE and WKE simulations. The conceptual framing — that the logarithmic divergence of the KZ spectrum in 2D acoustics places the system at a crossover between dimensional-analysis-determined and eigenvalue-determined self-similarity — is genuinely new and well-motivated. The derivation of the self-similar ansatz (6) and the eigenvalue equation (8) from the WKE is clean. The spectral collapse in Fig. 3a across GPE and WKE, with and without damping, is convincing evidence for a universal self-similar form Φ. The differential approximation yielding the exponential UV decay (9) is a nice analytical touch, and the dissipation-independence of Φ is a nontrivial and elegant result. Credit is due for running both GPE and WKE simulations and showing quantitative agreement on the front dynamics — that's the right comparison to make. The WavKinS code is publicly available, which matters for reproducibility. Now the soft spots. The stress-test concern about circularity in β extraction is real but somewhat overstated. The standard approach for second-kind self-similarity is exactly this: measure the front dynamics, fit to the predicted form, extract the eigenvalue. What's missing is not the fitting method but a consistency check — the paper never verifies that the measured Φ(η) and the fitted β jointly satisfy the eigenvalue equation (8). That's a one-figure addition that would close the loop. Without it, a reader can't confirm the system is in the asymptotic self-similar regime rather than a transient. The WKE validity concern is acknowledged honestly by the authors (page 3) and is inherent to 2D acoustic wave turbulence — not a flaw of this paper specifically. The GPE domain (1024ξ) is reasonable but the asymptotic regime may not be fully reached. No uncertainty on β is reported, which is minor but should be fixed. The vortex-free assumption is stated but not explicitly verified in the GPE runs — worth a sentence. This paper is for researchers in wave turbulence, non-equilibrium BEC dynamics, and self-similar solutions. It deserves a serious referee who can check the SI derivation and push for the Φ–β consistency verification. The central result is likely correct; the gaps are addressable without changing the conclusions.

Referee Report

2 major / 6 minor

Summary. This manuscript studies the self-similar evolution of the energy spectrum of a 2D Bose-Einstein condensate in the acoustic (weakly nonlinear) regime. Starting from the Gross-Pitaevskii equation (GPE), the authors invoke the wave kinetic equation (WKE) to derive a self-similar ansatz for the energy spectrum. They show that the propagating front $k^*(t)$ obeys a universal law governed by a dimensionless constant $beta$, which arises as a nonlinear eigenvalue from the self-similarity equation. The front dynamics are of second-kind self-similarity (exponential in time for the conservative case), while the system sits at a borderline of infinite capacity due to the logarithmic divergence of the KZ spectrum. Numerical simulations of both the GPE and the WKE are presented, showing collapse of the compensated spectra onto a universal function $Phi(eta)$ and quantitative agreement with the predicted front propagation law, yielding $beta approx 8.9$. The self-similar form is shown to be independent of dissipation, making the result applicable to both atomic and polariton BECs.

Significance. The paper addresses a well-posed and timely question in wave turbulence and non-equilibrium BEC dynamics. The identification of a universal constant $beta$ characterizing the self-similar front propagation in 2D acoustic turbulence is a concrete, falsifiable prediction. The demonstration that the self-similar form $Phi$ is independent of dissipation is a non-trivial result with direct experimental relevance for polariton systems. The use of two independent numerical frameworks (GPE via FROST and WKE via WavKinS) to cross-validate the theory is a strength. The analytical derivation of the UV asymptotic behavior via a differential approximation provides a clean, parameter-free prediction for $Phi(eta to infty)$ that is confirmed numerically. The hybrid first/second-kind self-similarity classification is an interesting conceptual contribution to the wave turbulence literature.

major comments (2)

The central claim is that $beta approx 8.9$ is a universal constant arising as a nonlinear eigenvalue from Eq. (8): $-beta[Phi(eta)+eta Phi'(eta)] = mathrm{St}[Phi]$. However, $beta$ is never determined by independently solving this eigenvalue problem. Instead, $beta$ is extracted by fitting the measured $k^*(t)$ to Eq. (10), which is itself derived from the self-similar ansatz combined with Eqs. (7)-(8). This introduces a methodological gap: the fit assumes self-similarity holds, and the resulting $beta$ is only meaningful as an eigenvalue if the measured $Phi(eta)$ and $beta$ jointly satisfy Eq. (8). The paper presents the spectral collapse (Fig. 3a) as confirming $Phi$ and the $k^*(t)$ fit (Fig. 3b) as determining $beta$, but never verifies that these two independently extracted quantities are consistent with the eigenvalue equation that links them. A direct check — substituting theme
The GPE simulation domain (1024$xi$) and resolution may not be sufficient to fully resolve the asymptotic self-similar regime, particularly the exponential UV tail of $Phi(eta)$ predicted by Eq. (9). The paper acknowledges (page 3) that the WKE has 'serious mathematical issues' for 2D acoustic waves, and the agreement between GPE and WKE in Fig. 3 is the primary evidence for the universality claim. However, the GPE data in Fig. 3a appears to span a more limited range of $eta$ than the WKE data. The authors should clarify the range of $eta$ over which the GPE and WKE data overlap and quantitatively agree, and discuss whether finite-size or finite-resolution effects in the GPE simulation could affect the extracted value of $beta$. A quantitative comparison of $beta$ extracted from GPE alone versus WKE alone, with error bars or uncertainty estimates, would strengthen the universality claim.

minor comments (6)

Page 3: The statement that 'the WWTT presents serious mathematical issues for acoustic waves, notably in 2D' is important context but is briefly mentioned and then set aside. A more explicit discussion of what these issues are and why they do not affect the quantitative predictions would help the reader.
Eq. (2): The choice of powers ($k^3$ and $k$) in the definition of $k^*$ is stated to be 'the lowest ones such that the average does not contain the logarithmic term.' This reasoning is somewhat opaque on first reading; a brief expansion in the SI would improve readability.
Fig. 3a: The axis labels and legend are difficult to read in the current rendering. The inset showing the exponential asymptotic behavior would benefit from clearer annotation of the predicted slope $-beta/D$.
The value $beta approx 8.9$ is reported without an uncertainty estimate. Given that it is extracted from numerical fits, an estimate of the statistical or systematic error would be appropriate.
Reference [30] (Costa et al., arXiv:2508.09799) appears to be a companion or prior work by the same authors. The relationship between the results presented here and that work should be clarified.
SI, Eq. (13): The diffusion coefficient $D = 2int_0^epsilon eta Phi(eta) deta$ depends on the measured $Phi$, but the exponential prediction $Phi sim exp[-beta eta / D]$ is then used to confirm the numerical $Phi$. The interplay between the measured $D$ and the predicted form should be made explicit (i.e., is $D$ measured independently or self-consistently?). Please check.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for a careful reading and constructive comments. Both major points are well-taken. On the first, we agree that a direct verification of the eigenvalue equation (8) using the independently extracted Phi and beta is a necessary consistency check and will add it. On the second, we agree that a more quantitative comparison of beta extracted separately from GPE and WKE, with error estimates and a discussion of finite-size/resolution effects, is needed. We will revise the manuscript accordingly.

read point-by-point responses

Referee: The central claim is that beta ≈ 8.9 is a universal constant arising as a nonlinear eigenvalue from Eq. (8). However, beta is never determined by independently solving this eigenvalue problem. Instead, beta is extracted by fitting the measured k*(t) to Eq. (10), which is itself derived from the self-similar ansatz combined with Eqs. (7)-(8). This introduces a methodological gap: the fit assumes self-similarity holds, and the resulting beta is only meaningful as an eigenvalue if the measured Phi(eta) and beta jointly satisfy Eq. (8). The paper presents the spectral collapse (Fig. 3a) as confirming Phi and the k*(t) fit (Fig. 3b) as determining beta, but never verifies that these two independently extracted quantities are consistent with the eigenvalue equation that links them. A direct check — substituting the measured Phi and beta into Eq. (8) — is missing.

Authors: The referee is correct that the manuscript does not currently include a direct verification that the independently extracted Phi and beta jointly satisfy the eigenvalue equation (8). This is a genuine methodological gap: we determine Phi from the spectral collapse and beta from the front-propagation fit, but we do not explicitly close the loop by substituting both back into Eq. (8) to check consistency. We agree this check is needed and will add it in the revised manuscript. Specifically, we will substitute the numerically measured Phi(eta) and the fitted beta into the left-hand side, -beta[Phi + eta Phi'], and compare it to the collision integral St[Phi] evaluated on the same measured Phi. This will provide a direct, quantitative test that the two independently extracted quantities are mutually consistent. We expect agreement within numerical resolution, but if discrepancies arise (e.g., in the far UV where finite-resolution effects are strongest), we will report them honestly. We note that solving the eigenvalue problem (8) independently as a boundary-value problem to extract beta without reference to the front dynamics would be an even stronger test, and we will explore whether this is numerically feasible within the revision timeframe. However, we cannot guarantee at this stage that such a direct eigenvalue solve will converge, given the known mathematical subtleties of the 2D acoustic collision integral. At minimum, the substitution check will be performed and reported. revision: yes
Referee: The GPE simulation domain (1024 xi) and resolution may not be sufficient to fully resolve the asymptotic self-similar regime, particularly the exponential UV tail of Phi(eta) predicted by Eq. (9). The paper acknowledges (page 3) that the WKE has 'serious mathematical issues' for 2D acoustic waves, and the agreement between GPE and WKE in Fig. 3 is the primary evidence for the universality claim. However, the GPE data in Fig. 3a appears to span a more limited range of eta than the WKE data. The authors should clarify the range of eta over which the GPE and WKE data overlap and quantitatively agree, and discuss whether finite-size or finite-resolution effects in the GPE simulation could affect the extracted value of beta. A quantitative comparison of beta extracted from GPE alone versus WKE alone, with error bars or uncertainty estimates, would strengthen the universality claim.

Authors: This is a fair and important point. The GPE simulation is indeed more limited in dynamic range than the WKE simulation, which spans four decades in wavenumber. The GPE domain of 1024 xi constrains the accessible range of k, and finite resolution limits the UV tail that can be resolved. We will address this in the revision by: (1) explicitly stating the range of eta over which GPE and WKE data overlap and quantitatively agree in Fig. 3a; (2) discussing the finite-size and finite-resolution effects that limit the GPE's ability to resolve the exponential UV tail; (3) providing separate fits of beta from GPE-only and WKE-only data with uncertainty estimates. We expect the WKE value to be more reliable given its larger spectral range, and the GPE value to be consistent within error bars but with larger uncertainty. We will state this transparently. We acknowledge that we cannot fully rule out systematic effects in the GPE extraction due to the limited dynamic range, and the revised manuscript will make this limitation clear rather than overstating the GPE-WKE agreement. revision: yes

Circularity Check

0 steps flagged

Minor self-citations in the WKE derivation chain; the central self-similarity derivation is self-contained and β is transparently fitted, not circularly predicted.

full rationale

The paper's derivation chain proceeds as follows: (1) the 2D GPE is mapped to the WKE (Eq. 3), whose original derivation traces to Zakharov–Sagdeev 1970 [29, external]; (2) a self-similar ansatz (Eq. 6) is substituted into the WKE, yielding the front equation (Eq. 7) and the eigenvalue equation (Eq. 8) — this step is algebraic and self-contained; (3) Eq. (7) is integrated to give the front law (Eq. 10); (4) β is extracted by fitting numerical k*(t) data to Eq. (10). The skeptic's concern that β is fitted rather than independently solved from the eigenvalue problem (8) is a legitimate verification gap — the paper never checks that the measured Φ and fitted β jointly satisfy Eq. (8). However, this is not circularity in the strict sense: the fit of k*(t) to Eq. (10) tests a non-trivial functional prediction (exponential/double-exponential front propagation) that could have failed, and β is a parameter within that prediction, not a quantity defined in terms of itself. The paper is transparent that β is 'determined numerically' and 'fitted,' not predicted from first principles. The self-citations [30, 31] (Costa/Krstulovic/Nazarenko) appear in the WKE derivation chain, but the core WKE framework originates from external work [29], and the self-cited results [31, PRL] are published and externally reviewed extensions of that framework. These self-citations provide the specific form of the collision integral and KZ solution but are not unverified claims on which the argument's validity depends. No step in the derivation reduces to its own inputs by construction. The verification gap regarding Eq. (8) is a correctness/completeness risk, not a circularity pattern.

Axiom & Free-Parameter Ledger

3 free parameters · 4 axioms · 0 invented entities

The paper introduces no new physical entities, particles, forces, or dimensions. The universal constant β is a new numerical result, not a new entity. The WKE, KZ spectrum, and self-similar ansatz are all established constructs in wave turbulence theory. The axiom ledger captures the key assumptions: WKE validity (acknowledged as problematic in 2D), the self-similar ansatz (checked a posteriori), the differential approximation (standard but approximate), and the vortex-free condition (assumed but not explicitly verified).

free parameters (3)

β = 8.9
Dimensionless eigenvalue arising from the self-similarity equation (8). Not determined analytically; extracted by fitting numerical k*(t) data to Eq. (10). This is the central universal constant of the paper.
k_0
IR cutoff set by the initial condition, introduced in Eq. (5) to regularize the logarithmic divergence of the KZ spectrum. Its value depends on the initial spectrum but is not a fitting parameter for the self-similar form.
D
Diffusion coefficient in the differential approximation, defined as D = 2∫_0^ε ηΦ(η)dη. Depends on the self-similar function Φ and is not independently fitted; it appears in the UV asymptotic prediction (Eq. 9) but does not affect the front dynamics.

axioms (4)

domain assumption The 2D wave kinetic equation (3) quantitatively describes the energy spectrum evolution of the 2D GPE in the acoustic regime.
Invoked throughout the paper as the basis for the self-similar derivation. The authors acknowledge (page 3) that 'the WWTT presents serious mathematical issues for acoustic waves, notably in 2D' and that the WKE derivation is rigorous only without forcing/dissipation. The damping term is added ad-hoc.
domain assumption The self-similar ansatz (6) with Φ(η) ~ η^{-1} as η→0 is consistent with the WKE dynamics.
Stated as an assumption to be checked a posteriori (page 3-4). The a posteriori check is the numerical collapse in Fig. 3a and the asymptotic analysis in the SI showing that the collision integral vanishes faster than the LHS for power-law ansätze.
ad hoc to paper The differential approximation used to derive the UV asymptotic behavior (Eq. 9) is valid for large η.
The scale separation η_1 << η, Φ(η_1) >> Φ(η) is assumed in the SI to simplify the collision integral. This is a standard technique in wave turbulence [34-36] but is an approximation, not an exact result.
domain assumption The GPE simulations are vortex-free.
The mapping from GPE to acoustic wave turbulence requires the absence of vortices. The paper states (page 5) 'in the absence of vortices, the GPE can be mapped into an acoustic wave turbulence problem.' The simulations use weak random Bogoliubov wave perturbations, but vortex formation is not explicitly monitored or excluded.

pith-pipeline@v1.1.0-glm · 15246 in / 2764 out tokens · 547878 ms · 2026-07-08T17:36:39.349028+00:00 · methodology

0 comments

read the original 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.