pith. sign in

arxiv: 2407.08738 · v3 · submitted 2024-07-11 · ❄️ cond-mat.quant-gas · cond-mat.stat-mech· physics.atom-ph· physics.flu-dyn

An Equation of State for Turbulence in the Gross-Pitaevskii model

Gevorg Martirosyan , Kazuya Fujimoto , Nir Navon This is my paper

Pith reviewed 2026-05-23 23:18 UTC · model grok-4.3

classification ❄️ cond-mat.quant-gas cond-mat.stat-mechphysics.atom-phphysics.flu-dyn
keywords Gross-Pitaevskii modelquantum turbulenceequation of stateenergy cascademomentum distributionBose-Einstein condensatefar-from-equilibrium
0
0 comments X

The pith

The Gross-Pitaevskii model encodes a universal equation of state for mixed turbulence that links momentum distribution amplitude to energy flux as n0 proportional to epsilon to the 0.67 power.

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

The paper shows that simulations of the Gross-Pitaevskii model produce a far-from-equilibrium equation of state for turbulence. In the regime where vortices and waves both contribute, the amplitude of the momentum distribution follows a specific power-law relation to the energy flux that existing turbulence theories do not predict. The same simulations confirm a direct energy cascade through explicit flux calculations and indicate that quasi-static thermodynamic processes can connect different far-from-equilibrium steady states.

Core claim

In the Gross-Pitaevskii model the turbulent cascade in the mixed regime obeys an equation of state in which the momentum distribution amplitude n0 scales with the energy flux ε according to n0 ∝ ε^{0.67(2)} in the steady state. This scaling is obtained after verifying that the state exhibits the signatures of a direct energy cascade and that the momentum distribution is consistent with wave-turbulent kinetic theory in the relevant limits.

What carries the argument

The scaling relation n0 ∝ ε^{0.67(2)} between the amplitude of the turbulent momentum distribution and the energy flux, extracted from steady-state simulations that also compute the energy and particle fluxes.

If this is right

  • The turbulent state possesses the hallmarks of a direct energy cascade as shown by the signs of the computed energy and particle fluxes.
  • The observed scaling matches a recent experimental measurement performed at large energy fluxes.
  • Quasi-static thermodynamic processes between equilibrium states extend to far-from-equilibrium steady states.

Where Pith is reading between the lines

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

  • The same scaling may appear in other models that support both vortex and wave dynamics at comparable amplitudes.
  • The exponent 0.67 could be recovered analytically by combining wave-kinetic and vortex contributions in a single closure for the Gross-Pitaevskii equation.
  • Control of the energy flux in an experiment would allow direct tuning of the turbulent amplitude according to the reported power law.

Load-bearing premise

The calculated fluxes correctly identify a direct energy cascade and the simulated state is free of significant numerical artifacts or misidentification of the mixed-turbulence regime.

What would settle it

A higher-resolution simulation or controlled experiment that measures a statistically different exponent in the n0-versus-ε relation would falsify the claimed scaling.

Figures

Figures reproduced from arXiv: 2407.08738 by Gevorg Martirosyan, Kazuya Fujimoto, Nir Navon.

Figure 1
Figure 1. Figure 1: FIG. 1. The direct turbulent cascade in the GP model. (a) Cartoon of the simulation geometry and the driving protocol. We use a cylindrical [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Energetics of the direct turbulent cascade. (a) Energy input and dissipation rates. We show [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. (a) A universal equation of state of the GP model. Expressing the state variables [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
read the original abstract

We report the numerical observation of a far-from-equilibrium equation of state (EOS) in the Gross-Pitaevskii model. We first show that the momentum distribution of the turbulent cascade is well described by wave-turbulent kinetic theory in the appropriate limits. Calculating the energy and particle fluxes $\Pi_\varepsilon(k)$ and $\Pi_N(k)$, we show that the turbulent state possesses the hallmarks of a direct energy cascade. Building on this, we show that the GP model encodes a universal EOS in the form of a relationship between the turbulent cascade's momentum distribution amplitude $n_0$ and the energy flux $\epsilon$ in the steady state. We find that in our regime of `mixed' turbulence - where both vortices and waves play a significant role - $n_0\propto \epsilon^{0.67(2)}$, a result that is not captured by any existing theory of turbulence but that agrees with a recent experimental measurement for large energy fluxes. Finally, we find that the concept of quasi-static thermodynamic processes between equilibrium states extends to far-from-equilibrium steady states.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper numerically observes a far-from-equilibrium equation of state in the Gross-Pitaevskii model. It first verifies that the momentum distribution of the turbulent cascade matches wave-turbulent kinetic theory in appropriate limits, then computes energy and particle fluxes Π_ε(k) and Π_N(k) to establish the hallmarks of a direct energy cascade. Building on this, it reports a universal relation n0 ∝ ε^{0.67(2)} between the momentum distribution amplitude n0 and energy flux ε in the steady state for the 'mixed' turbulence regime (vortices + waves), a result not captured by existing theory but consistent with a recent experiment at large fluxes. It further claims that quasi-static thermodynamic processes extend to far-from-equilibrium steady states.

Significance. If the central numerical relation holds under robust checks, the work supplies an empirical EOS specific to mixed turbulence in the GP model, directly linking numerics to experiment and extending thermodynamic concepts beyond equilibrium. The agreement with experiment for large fluxes and the use of flux diagnostics are positive features; the result is empirical rather than derived from first principles.

major comments (2)
  1. [Section on flux calculations and cascade identification (discussion of Π_ε(k) and Π_N(k))] The extraction of n0 ∝ ε^{0.67(2)} rests on identifying an inertial interval where Π_ε(k) is constant and Π_N(k) ≈ 0. The manuscript does not quantify the robustness of the chosen k-interval: no tests are shown for how modest shifts in the lower/upper cutoffs affect the fitted exponent relative to the quoted (2) uncertainty. This is load-bearing for the universality claim.
  2. [Section characterizing the turbulent state as mixed] The regime is labeled 'mixed' turbulence on the basis that both vortices and waves play a significant role, yet the diagnostics used to establish this classification (e.g., relative vortex vs. wave contributions, or momentum-distribution diagnostics) and their sensitivity to simulation parameters are not specified. Without this, it is unclear whether the reported exponent is tied to a well-defined regime or could shift under reclassification.
minor comments (2)
  1. [Abstract] The abstract states the exponent with uncertainty but omits simulation parameters, grid sizes, or data-exclusion criteria; adding these would improve verifiability.
  2. [Notation and figure captions] Clarify the precise definition of the amplitude n0 (e.g., prefactor of the k^{-x} fit to the momentum distribution) and ensure consistent notation between text, figures, and equations.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments on our manuscript. We address each major comment below and will make revisions to improve the robustness and clarity of our presentation.

read point-by-point responses

  1. Referee: The extraction of n0 ∝ ε^{0.67(2)} rests on identifying an inertial interval where Π_ε(k) is constant and Π_N(k) ≈ 0. The manuscript does not quantify the robustness of the chosen k-interval: no tests are shown for how modest shifts in the lower/upper cutoffs affect the fitted exponent relative to the quoted (2) uncertainty. This is load-bearing for the universality claim.

    Authors: We agree that explicit quantification of the sensitivity to the inertial-interval cutoffs would strengthen the result. In the revised manuscript we will add a supplementary analysis (or main-text figure) that varies the lower and upper k cutoffs over a modest range and shows the resulting variation in the fitted exponent. Our internal checks indicate that the exponent remains within the quoted 0.67(2) uncertainty for reasonable choices, but we will document this explicitly. revision: yes

  2. Referee: The regime is labeled 'mixed' turbulence on the basis that both vortices and waves play a significant role, yet the diagnostics used to establish this classification (e.g., relative vortex vs. wave contributions, or momentum-distribution diagnostics) and their sensitivity to simulation parameters are not specified. Without this, it is unclear whether the reported exponent is tied to a well-defined regime or could shift under reclassification.

    Authors: We accept that a more precise operational definition of the 'mixed' regime is required. We will expand the relevant section to specify the quantitative diagnostics (e.g., decomposition of the momentum distribution into vortex and wave components, real-space vortex identification, and their dependence on forcing amplitude and dissipation) and will discuss the parameter range over which the classification remains stable. This will clarify the regime in which the reported EOS applies. revision: yes

Circularity Check

0 steps flagged

No circularity: empirical numerical observation of power-law relation

full rationale

The paper reports a numerical observation of the relation n0 ∝ ε^{0.67(2)} extracted from GP simulations after verifying flux constancy for a direct cascade. This is framed as an empirical finding with explicit uncertainty that agrees with an external experiment, not as a derivation or prediction that reduces to fitted inputs, self-citations, or definitions by construction. The provided text contains no load-bearing self-citations, no ansatz smuggled via prior work, and no renaming of known results as new derivations. The central claim remains an independent numerical result.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The claim rests on the standard Gross-Pitaevskii model plus empirical fitting of simulation outputs; no new entities are introduced.

free parameters (1)
  • EOS exponent = 0.67(2)
    Fitted numerically from the mixed-turbulence data to obtain 0.67(2)
axioms (2)
  • domain assumption Gross-Pitaevskii equation governs the condensate dynamics
    Invoked as the underlying model for the simulations
  • domain assumption Wave-turbulent kinetic theory describes the momentum distribution in appropriate limits
    Used to interpret the cascade spectrum

pith-pipeline@v0.9.0 · 5738 in / 1360 out tokens · 29097 ms · 2026-05-23T23:18:23.667177+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Turbulence and far-from-equilibrium equation of state of Bogoliubov waves in Bose-Einstein Condensates

    cond-mat.quant-gas 2024-08 unverdicted novelty 6.0

    Derives new stationary spectra for short and long Bogoliubov waves in 3D BEC turbulence and uses them to interpret the far-from-equilibrium equation of state from Dogra et al. 2023.

Reference graph

Works this paper leans on

45 extracted references · 45 canonical work pages · cited by 1 Pith paper

  1. [1]

    P. C. Hohenberg and B. I. Halperin, Theory of dynamic critical phenomena, Rev. Mod. Phys. 49, 435 (1977)

    work page 1977

  2. [2]

    A. J. Bray, Theory of phase-ordering kinetics, Advances in Physics 43, 357 (1994)

    work page 1994

  3. [3]

    Forster, Hydrodynamic fluctuations, broken symmetry, and correlation functions (CRC Press, 2018)

    D. Forster, Hydrodynamic fluctuations, broken symmetry, and correlation functions (CRC Press, 2018)

    work page 2018

  4. [4]

    Berges, A

    J. Berges, A. Rothkopf, and J. Schmidt, Nonthermal Fixed Points: Effective Weak Coupling for Strongly Correlated Sys- tems Far from Equilibrium, Phys. Rev. Lett. 101, 041603 (2008)

    work page 2008

  5. [5]

    S. Erne, R. Bücker, T. Gasenzer, J. Berges, and J. Schmied- mayer, Universal dynamics in an isolated one-dimensional Bose gas far from equilibrium, Nature 563, 225 (2018)

    work page 2018

  6. [6]

    Prüfer, P

    M. Prüfer, P. Kunkel, H. Strobel, S. Lannig, D. Linnemann, C.-M. Schmied, J. Berges, T. Gasenzer, and M. K. Oberthaler, Observation of universal dynamics in a spinor Bose gas far from equilibrium, Nature 563, 217 (2018)

    work page 2018

  7. [7]

    J. A. P. Glidden, C. Eigen, L. H. Dogra, T. A. Hilker, R. P. Smith, and Z. Hadzibabic, Bidirectional dynamic scaling in an isolated Bose gas far from equilibrium, Nat. Phys. 17, 457 (2021)

    work page 2021

  8. [8]

    Gałka, P

    M. Gałka, P. Christodoulou, M. Gazo, A. Karailiev, N. Dogra, J. Schmitt, and Z. Hadzibabic, Emergence of Isotropy and Dy- namic Scaling in 2D Wave Turbulence in a Homogeneous Bose Gas, Phys. Rev. Lett. 129, 190402 (2022)

    work page 2022

  9. [9]

    M. Gazo, A. Karailiev, T. Satoor, C. Eigen, M. Gałka, and Z. Hadzibabic, Universal Coarsening in a Homogeneous Two- Dimensional Bose Gas, arXiv:2312.09248 (2023)

    work page arXiv 2023

  10. [10]

    S. Huh, K. Mukherjee, K. Kwon, J. Seo, J. Hur, S. I. Mistakidis, H. R. Sadeghpour, and J.-y. Choi, Universality class of a spinor Bose-Einstein condensate far from equilibrium, Nat. Phys. 20, 402 (2024)

    work page 2024

  11. [11]

    S. C. Takatori and J. F. Brady, Towards a thermodynamics of active matter, Phys. Rev. E91, 032117 (2015)

    work page 2015

  12. [12]

    Ginot, I

    F. Ginot, I. Theurkauff, D. Levis, C. Ybert, L. Bocquet, L. Berthier, and C. Cottin-Bizonne, Nonequilibrium Equation of State in Suspensions of Active Colloids, Phys. Rev. X 5, 011004 (2015)

    work page 2015

  13. [13]

    L. H. Dogra, G. Martirosyan, T. A. Hilker, J. A. P. Glidden, J. Etrych, A. Cao, C. Eigen, R. P. Smith, and Z. Hadzibabic, Universal equation of state for wave turbulence in a quantum gas, Nature 620, 521 (2023)

    work page 2023

  14. [14]

    Navon, A

    N. Navon, A. L. Gaunt, R. P. Smith, and Z. Hadzibabic, Emer- gence of a turbulent cascade in a quantum gas, Nature 539, 72 (2016)

    work page 2016

  15. [15]

    Navon, C

    N. Navon, C. Eigen, J. Zhang, R. Lopes, A. L. Gaunt, K. Fu- jimoto, M. Tsubota, R. P. Smith, and Z. Hadzibabic, Synthetic dissipation and cascade fluxes in a turbulent quantum gas, Sci- ence 366, 382 (2019)

    work page 2019

  16. [16]

    Fibich, The nonlinear Schrödinger equation , V ol

    G. Fibich, The nonlinear Schrödinger equation , V ol. 192 (Springer, 2015)

    work page 2015

  17. [17]

    Pitaevskii and S

    L. Pitaevskii and S. Stringari, Bose–Einstein condensation and superfluidity (Oxford University Press, Oxford, 2016)

    work page 2016

  18. [18]

    [19] and references therein)

    There also exists a vast body of work on the topic of vortex tur- bulence in the GP model (see e.g. [19] and references therein)

    work page

  19. [19]

    Kobayashi, P

    M. Kobayashi, P. Parnaudeau, F. Luddens, C. Lothodé, L. Danaila, M. Brachet, and I. Danaila, Quantum turbu- lence simulations using the gross–pitaevskii equation: high- performance computing and new numerical benchmarks, Com- puter Physics Communications 258, 107579 (2021)

    work page 2021

  20. [20]

    Dyachenko, A

    S. Dyachenko, A. Newell, A. Pushkarev, and V . Zakharov, Op- tical turbulence: weak turbulence, condensates and collapsing filaments in the nonlinear schrödinger equation, Physica D 57, 96 (1992)

    work page 1992

  21. [21]

    Falkovich and I

    G. Falkovich and I. Ryzhenkova, Kolmogorov spectra of lang- muir and optical turbulence, Physics of Fluids B: Plasma Physics 4, 594 (1992)

    work page 1992

  22. [22]

    Y . Zhu, B. Semisalov, G. Krstulovic, and S. Nazarenko, Direct and inverse cascades in turbulent bose-einstein condensates, Phys. Rev. Lett. 130, 133001 (2023)

    work page 2023

  23. [23]

    V . E. Zakharov, V . S. L’V ov, and G. Falkovich, Kolmogorov spectra of turbulence (Springer Berlin, 1992)

    work page 1992

  24. [24]

    G. P. Agrawal, Nonlinear fiber optics, in Nonlinear Science at the Dawn of the 21st Century (Springer, 2000) pp. 195–211

    work page 2000

  25. [25]

    P. G. Kevrekidis, D. J. Frantzeskakis, and R. Carretero- González, Emergent nonlinear phenomena in Bose-Einstein condensates: theory and experiment, V ol. 45 (Springer, 2008)

    work page 2008

  26. [26]

    Craig, C

    W. Craig, C. Sulem, and P.-L. Sulem, Nonlinear modulation of gravity waves: a rigorous approach, Nonlinearity5, 497 (1992)

    work page 1992

  27. [27]

    Supplemental Material

    work page

  28. [28]

    Navon, R

    N. Navon, R. P. Smith, and Z. Hadzibabic, Quantum gases in optical boxes, Nat. Phys. 17, 1334 (2021)

    work page 2021

  29. [29]

    Note that even though the trap shape affects large-scale prop- erties, the far-from-equilibrium state at smaller scales is largely insensitive to the details of the confining potential

    work page

  30. [30]

    S. J. Garratt, C. Eigen, J. Zhang, P. Turzák, R. Lopes, R. P. Smith, Z. Hadzibabic, and N. Navon, From single-particle ex- citations to sound waves in a box-trapped atomic bose-einstein condensate, Phys. Rev. A 99, 021601 (2019)

    work page 2019

  31. [31]

    In [27], we derive an analytic expression for tD in terms of the cascade properties

    work page

  32. [32]

    Nazarenko, Wave turbulence (Springer, 2011)

    S. Nazarenko, Wave turbulence (Springer, 2011)

    work page 2011

  33. [33]

    In our simulations, a description in terms of an equilibrium-like condensate is meaningless since for our strong drives, the Bogoliubov approximation is invalid

    In the presence of a condensate, the proper quasiparticles to consider for kξ ≲ 1 are the Bogoliubov phonons, for which the appropriate WWT description is a wave-kinetic equation with 3-wave interactions. In our simulations, a description in terms of an equilibrium-like condensate is meaningless since for our strong drives, the Bogoliubov approximation is invalid

    work page

  34. [34]

    Interestingly, a recent experiment has observed that the mo- mentum distribution of 2D shaken Bose gases also becomes isotropic at k ≈ kξ [8]

    work page

  35. [35]

    Experimentally, kD is bounded, so the separation between the injection scale k0 ∝ kξ and kD – and hence the region where the 4-wave WWT theory is applicable – is reduced for stronger interactions

    work page

  36. [36]

    Even though α ̸= 1 means that the energy-density flux is not ˙NUD/V , the fact that α has only a very weak dependence on a (∝ a0.08(1)) and is independent ofUs suggests that the determi- nation of the energy flux via ˙NUD/V , as done in [13, 15], will only introduce a small, presumably system-specific, systematic error

    work page

  37. [37]

    Note that the slow decrease of the plateau of Πε is due to a pro- gressive depletion of the condensate, which results in a (slow) decrease of ϵin [corresponding horizontal dashed lines in the upper panel of Fig. 2(c)]

    work page

  38. [38]

    Note that it is only in the limit kD → ∞ that one recovers the expected limits Πε → ϵ with ΠN → 0 for the direct energy cas- cade [20]

    work page

  39. [39]

    3, but is used in the calculation of n0 = Nk k−3(kξt )−0.5

    Note that ξt is not directly present in the axes of Fig. 3, but is used in the calculation of n0 = Nk k−3(kξt )−0.5

    work page

  40. [40]

    The independence of our EoS with kD cannot rigorously hold as for large enough kD (presumably beyond what we explored), WWT results should be recovered

    work page

  41. [41]

    The factor of ξ0.5 on the y-axis is a finite-kD effect, but even in 6 its absence our EoS would depend on n; within the GP model, only the WWT prediction of a power-law EoS with an exponent of 1/3 is independent of n

    work page

  42. [42]

    Frisch, Turbulence: the legacy of A.N

    U. Frisch, Turbulence: the legacy of A.N. Kolmogorov (Cam- bridge University Press, 1995)

    work page 1995

  43. [43]

    As shown in [27], n0 and ϵ/n are roughly constant over time, resulting in the increase of n0/n and τt ϵ/(nζt )

    work page

  44. [44]

    H. B. Callen, Thermodynamics and an Introduction to Thermo- statistics (John Wiley & Sons, 1991)

    work page 1991

  45. [45]

    SUPPLEMENTAL MA TERIAL I

    We benchmarked our simulations for different choices of A and w and have verified their choice within an appropriate range does not affect our results. SUPPLEMENTAL MA TERIAL I. S IMULATION DETAILS Our numerical simulations solve the time-dependent GP equation [Eq. (1)] using a pseudo-spectral method with the fourth- order Runge-Kutta time evolution. We a...

    work page