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
β, which we determine numerically.
We study the universal nonstationary evolution of wave turbulence (WT) in Bose-Einstein condensates (BECs). Their temporal evolution can exhibit different kinds of self-similar behavior corresponding to a large-time asymptotic of the system or to a finite-time blowup. We identify self-similar regimes in BECs by numerically simulating the forced and unforced Gross-Pitaevskii equation (GPE) and the associated wave kinetic equation (WKE) for the direct and inverse cascades, respectively. In both the GPE and the WKE simulations for the direct cascade, we observe the first-kind self-similarity that is fully determined by energy conservation. For the inverse cascade evolution, we verify the existence of a self-similar evolution of the second kind describing self-accelerating dynamics of the spectrum leading to blowup at the zero mode (condensate) at a finite time. We believe that the universal self-similar spectra found in the present paper are as important and relevant for understanding the BEC turbulence in past and future experiments as the commonly studied stationary Kolmogorov-Zakharov (KZ) spectra.
We investigate turbulent spin-down dynamics in self-gravitating Bose-Einstein condensates, comparing purely bosonic and axionic (higher-order interacting) systems. Through simulations of the Gross-Pitaevskii-Poisson system, we study condensates pinned to a crust potential undergoing rapid rotation slowdown. We find that axionic condensates exhibit more uniform density profiles and smaller sizes compared to their bosonic counterparts for similar interaction strengths, which facilitates earlier vortex entry. The sudden spin-down triggers vortex depinning and a turbulent cascade. For comparable sizes, both systems exhibit a short-lived Kolmogorov energy cascade (
k−5/3 scaling) followed by a transition to Vinen turbulence (
k−1 scaling). Crucially, their responses diverge with increasing interaction strength (and thus condensate size): the axionic system increasingly deviates from Kolmogorov scaling because of enhanced vortex retention, a trend quantitatively confirmed by analyzing the vortex fraction and its dependence on the final rotation frequency. Spectral analysis reveals that the growth of incompressible energy is primarily driven by quantum pressure during vortex detachment, rather than by compressible flows. The compressible spectrum shows thermalization (
k scaling). Our results demonstrate how distinct nonlinearities govern vortex dynamics and turbulent dissipation in self-gravitating quantum fluids.