Strong and weak wave turbulence regimes in Bose–Einstein condensates

Decaying turbulence in rotating Magneto-hydrodynamic systems is studied theoretically and numerically. In the linear limit, when the velocity and magnetic perturbations are small, the system supports two types of waves. When the rotation effects are stronger than the ones of the external magnetic field, one of these waves contains most of the kinetic energy (inertial wave) and the other—most of the magnetic energy (magnetostrophic wave). The weak wave turbulence (WWT) theory for decoupled inertial and magnetospheric wave systems was previously derived in Galtier (2014 J. Fluid Mech. 757 114–54). In the present paper, we suggest theory of strong turbulence for such waves based on the critical balance (CB) approach conjecturing that the linear and nonlinear timescales are of similar magnitudes in a wide range of turbulent scales. Regimes of weak and strong wave turbulence are simulated numerically. The results appear to be in good agreement with the WWT and CB predictions, particularly for the exponents of the kinetic and magnetic energy spectra.

A heuristic model is given for anisotropic magnetohydrodynamics turbulence in the presence of a uniform external magnetic field B0ê‖. The model is valid for both moderate and strong B0 and is able to describe both the strong and weak wave turbulence regimes as well as the transition between them. The main ingredient of the model is the assumption of constant ratio at all scales between the linear wave period and the nonlinear turnover time scale. Contrary to the model of critical balance introduced by Goldreich and Sridhar [Astrophys. J. 438, 763 (1995)], it is not assumed, in addition, that this ratio be equal to unity at all scales. This allows us to make use of the Iroshnikov–Kraichnan phenomenology; it is then possible to recover the widely observed anisotropic scaling law k‖∝k⊥2∕3 between parallel and perpendicular wave numbers (with reference to B0ê‖) and to obtain for the total-energy spectrum E(k⊥,k‖)∼k⊥−αk‖−β the universal prediction, 3α+2β=7. In particular, with such a prediction, the weak Alfvén wave turbulence constant-flux solution is recovered and, for the first time, a possible explanation to its precursor found numerically by Galtier et al. [J. Plasma Phys. 63, 447 (2000)] is given.

Universal dynamics of a dilute Bose gas is studied within the kinetic regime where the time evolution of the mode occupation number is governed by a wave-Boltzmann equation. The universality manifests itself in the form of the dynamical evolution which can be a cascade in the context of Kolmogorov-Zakharov wave turbulence or a self-similar shift in time and space. Which is the case depends on the relevant global conservation laws and on the particular range of scales in which the respective dynamics takes place. A nonperturbative kinetic equation is derived by applying the Schwinger-Keldysh closed-time functional integral and the 2-particle-irreducible formalism to an action of complex scalar Bose fields with quartic interaction.The resulting dynamic equation for two-point correlations are reduced to a wave-Boltzmann-type kinetic equation with an effective $T$-matrix which depends explicitly on the mode occupation numbers. This nonperturbative dependence on the solution occurs in the infrared regime where wave numbers are large and collective scattering of many particles prevails. Thus, the scaling analysis of wave turbulence theory can be applied. We explicitly calculate the effective $T$-matrix analytically taking into account an infrared cutoff for dealing with the infrared divergences. Our results show that the scaling behaviour of the $T$-matrix departs from the one predicted by naive dimensional counting due to the presence of the infrared cutoff. The Kolmogorov-Zakharov and the self-similar exponents are evaluated by power counting using our $T$-matrix and the results are confirmed by numerical integration of the wave-Boltzmann scattering integral. The scaling exponents governing the time evolution are determined by means of the global conservation laws and the kinetic equation. Depending on the scaling properties of the quasiparticle spectrum, the momentum cutoff scale in the infrared evolves critically slowed in time. The respective scaling exponent is universal in the nonperturbative regime regardless of whether the process is an inverse cascade or a self-similar shift towards the infrared. Thus, our results provide a general framework for classifying nonthermal fixed points in dilute ultracold Bose gases. They pave the way to a straightforward generalisation and application to trapped systems in different dimensionalities and to systems with more than one internal degree of freedom. Our results furthermore provide a possible interpretation of recent experimental results on wave turbulence in an ultracold Bose gas. %We studied the universal dynamics of diluted Bose gas within the kinetic regime where time evolution of occupation number is governed by the Boltzmann's scattering integral. %There are three possible scenarios, depending on the scaling behaviour of occupation number. They are the particle and energy cascade which are the steady state of turbulence transport and self-similar evolution where local occupation evolves in time. %We evaluated the scaling exponents in both momentum and time scales in each case under the constraints of global number and energy conservation. Furthermore, we can go beyond the perturbative regime by employing the closed-time-path contour and the 2-particle-irreducible formalism to the Gross-Pitaevskii action. The kinetic equation is modified by an effective T-matrix as a result of the resummation scheme. However, the structure of kinetic equation remains the same as in the perturbative case. Thus, all exponents in the nonperturbative regime depart from those in the perturbative regime because of the scaling of the effective T-matrix.

The performance of free space optical (FSO) links in a clear atmosphere is affected by the non-ideal characteristics of the communication channel. Atmospheric turbulence causes fluctuations in the received signal level, which increase the bit errors in a digital communication link. In order to quantify performance limitations, a better understanding of the effect of the intensity fluctuations on the received signal at all turbulence levels is needed. Theory reliably describes the behavior in the weak turbulence regime, but theoretical descriptions in the intermediate and strong turbulence regimes are less well developed. We have developed a flexible empirical approach for characterizing link performance in strong turbulence conditions through image analysis of intensity scintillation patterns coupled with frame aperture averaging on an FSO communication link. These measurements are complemented with direct measurements of temporal and spatial correlation functions. A He-Ne laser beam propagates 106 meters in free-space over flat terrain about a meter above the ground to provide strong atmospheric turbulence conditions. A high performance digital camera with a frame-grabbing computer interface is used to capture received laser intensity distributions at rates up to 30 frames per second and various short shutter speeds, down to 1/16,000s per frame. A scintillometer is used for accurate measurements of the turbulence parameter Cn 2 . Laboratory measurements use a local strong turbulence generator, which mimics a strong phase screen. Spatial correlation functions are measured using laterally separated point detectors placed in the receiver plane. Correlations and captured image frames are analyzed in Labview to evaluate correlation functions, Cn 2 , and the aperture averaging factor. Our correlation measurement technique represents a probe into the length scales of different flows and make for a “fingerprint” of the air conditions at the time. The aperture averaging results demonstrate the expected reduction in intensity fluctuations with increasing aperture diameter, and show quantitatively the differences in behavior between various strengths of turbulence. This paper will present accurate empirical data in the strong turbulence regime. Such results can help build upon existing empirical data and lead to the development of new theories.

When a Bose-Einstein condensate (BEC) is driven out of equilibrium, density waves interact nonlinearly and trigger turbulent cascades. In a turbulent BEC, energy is transferred toward small scales by a direct cascade, whereas the number of particles displays an inverse cascade toward large scales. In this work, we study analytically and numerically the direct and inverse cascades in wave-turbulent BECs. We analytically derive the Kolmogorov-Zakharov spectra, including the log correction to the direct cascade scaling and the universal prefactor constants for both cascades. We test and corroborate our predictions using high-resolution numerical simulations of the forced-dissipated Gross-Pitaevskii model in a periodic box and the corresponding wave-kinetic equation. Theoretical predictions and data are in excellent agreement, without adjustable parameters. Moreover, in order to connect with experiments, we test and validate our theoretical predictions using the Gross-Pitaevskii model with a confining cubic trap. Our results explain previous experimental observations and suggest new settings for future studies.

We report a numerical investigation of three dimensional, incompressible, Hall magnetohydrodynamic turbulence with a relatively strong mean magnetic field. Using helicity decomposition and cross-bicoherence analysis, we observe that the resonant three--wave coupling is substantial among ion cyclotron and whistler waves. A detailed study of the degree of non-linearity of these two populations shows that the ion cyclotron component experiences a transition from weak to strong wave turbulence going from large to small scales, while the whistler fluctuations display a weak wave turbulence character for all scales. This non-trivial coexistence of the two regimes with the two populations of waves gives rise to anomalous anisotropy and scaling properties. The weak and strong wave turbulence components can be distinguished rather efficiently using spatio-temporal Fourier transforms. The analysis shows that while resonant triadic interactions survive the highly non-linear bath of ion cyclotron fluctuations at large scales for which the degree of non-linearity is low for both populations of waves, whistler waves tend to be killed by the non-linear cross-coupling at smaller scales where the ion cyclotron component is in the strong wave turbulent regime. Such situation may have far-reaching implications for the physics of magnetized turbulence in many astrophysical and space plasmas where different waves coexist and compete to transfer non-linearly energy across scales.

It is widely accepted that the primordial universe experienced a brief period of accelerated expansion called inflation. This scenario provides a plausible solution to the horizon and flatness problems. However, the particle physics mechanism responsible for inflation remains speculative with, in particular, the assumption of a scalar field called inflaton. Furthermore, the comparison with the most recent data raises new questions that encourage the consideration of alternative hypotheses. Here, we propose a completely different scenario based on a mechanism whose origins lie in the nonlinearities of the Einstein field equations. We use the analytical results of weak gravitational wave turbulence to develop a phenomenological theory of strong gravitational wave turbulence where the inverse cascade of wave action plays a key role. In this scenario, the space-time metric excitation triggers an explosive inverse cascade followed by the formation of a condensate in Fourier space whose growth is interpreted as an expansion of the universe. Contrary to the idea that gravitation can only produce a decelerating expansion, our study reveals that strong gravitational wave turbulence could be a source of inflation. The fossil spectrum that emerges from this scenario is shown to be in agreement with the cosmic microwave background radiation measured by the Planck mission. Direct numerical simulations can be used to check our predictions and to investigate the question of non-Gaussianity through the measure of intermittency.

In a two-dimensional (2D) classical fluid, a large-scale flow structure emerges out of turbulence, which is known as the inverse energy cascade where energy flows from small to large length scales. An interesting question is whether this phenomenon can occur in a superfluid, which is inviscid and irrotational by nature. Atomic Bose-Einstein condensates (BECs) of highly oblate geometry provide an experimental venue for studying 2D superfluid turbulence, but their full investigation has been hindered due to a lack of the circulation sign information of individual quantum vortices in a turbulent sample. Here, we demonstrate a vortex sign detection method by using Bragg scattering, and we investigate decaying turbulence in a highly oblate BEC at low temperatures, with our lowest being ~0.5Tc, where Tc is the superfluid critical temperature. We observe that weak spatial pairing between vortices and antivortices develops in the turbulent BEC, which corresponds to the vortex-dipole gas regime predicted for high dissipation. Our results provide a direct quantitative marker for the survey of various 2D turbulence regimes in the BEC system.

Large-scale turbulence in fluid layers and other quasi-two-dimensional compressible systems consists of planar vortices and waves. Separately, wave turbulence usually produces a direct energy cascade, while solenoidal planar turbulence transports energy to large scales by an inverse cascade. Here, we consider turbulence at finite Mach numbers when the interaction between acoustic waves and vortices is substantial. We employ solenoidal pumping at intermediate scales and show how both direct and inverse energy cascades are formed starting from the pumping scale. We show that there is an inverse cascade of kinetic energy up to a scale $\ell$, where a typical velocity reaches the speed of sound; this creates shock waves, which provide for a compensating direct cascade. When the system size is less than $\ell$, the steady state contains a system-size pair of long-living condensate vortices connected by a system of shocks. Thus turbulence in fluid layers processes energy via a loop: Most energy first goes to large scales via vortices and is then transported by waves to small-scale dissipation.

In the present work, we investigate how single- and multi-vortex-ring states can emerge from a planar dark soliton in three-dimensional (3D) Bose-Einstein condensates (confined in isotropic or anisotropic traps) through bifurcations. We characterize such bifurcations quantitatively using a Galerkin-type approach, and find good qualitative and quantitative agreement with our Bogoliubov-de Gennes (BdG) analysis. We also systematically characterize the BdG spectrum of the dark solitons, using perturbation theory, and obtain a quantitative match with our 3D BdG numerical calculations. We then turn our attention to the emergence of single- and multi-vortex-ring states. We systematically capture these as stationary states of the system and quantify their BdG spectra numerically. We find that although the vortex ring may be unstable when bifurcating, its instabilities weaken and may even eventually disappear, for sufficiently large chemical potentials and suitable trap settings. For instance, we demonstrate the stability of the vortex ring for an isotropic trap in the large chemical potential regime.

We study the statistical properties of an ensemble of weak gravitational waves interacting nonlinearly in a flat space-time. We show that the resonant three-wave interactions are absent and develop a theory for four-wave interactions in the reduced case of a 2.5+1 diagonal metric tensor. In this limit, where only plus-polarized gravitational waves are present, we derive the interaction Hamiltonian and consider the asymptotic regime of weak gravitational wave turbulence. Both direct and inverse cascades are found for the energy and the wave action, respectively, and the corresponding wave spectra are derived. The inverse cascade is characterized by a finite-time propagation of the metric excitations-a process similar to an explosive nonequilibrium Bose-Einstein condensation, which provides an efficient mechanism to ironing out small-scale inhomogeneities. The direct cascade leads to an accumulation of the radiation energy in the system. These processes might be important for understanding the early Universe where a background of weak nonlinear gravitational waves is expected.

In this paper, we have investigated the performance of an underwater vertical wireless optical communication (UVWOC) link employing multiple input-multiple output (MIMO) operating in conjunction with equal gain combing (EGC) techniques perturbed by weak and strong turbulence in the presence of pointing errors and attenuation losses. Vertical underwater turbulence, which varies from layer to layer due to temperature and salinity variation connected to depth, is modeled using hyperbolic tangent log-normal (HTLN) distribution in the case of weak underwater turbulence and gamma-gamma (GG) distribution in the case of strong underwater turbulence. Novel closed-form expressions quantifying the average bit error rate (BER) have been derived for the UVWOC MIMO EGC system for weak and strong turbulence regimes. The expression for the average BER associated with the UVWOC link for different values of pointing error, differing vertical layer depth, modulation types, and differing numbers of sources and detectors have been determined. In addition, closed-form expressions for the outage probability (OP) and ergodic channel capacity (ECC) have been derived for the UVWOC MIMO EGC system. The accuracy of all closed-form expressions derived in the paper has been verified using Monte Carlo simulations.

A refined cascade model for kinetic turbulence in weakly collisional astrophysical plasmas is presented that includes both the transition between weak and strong turbulence and the effect of nonlocal interactions on the nonlinear transfer of energy. The model describes the transition between weak and strong MHD turbulence and the complementary transition from strong kinetic Alfvén wave (KAW) turbulence to weak dissipating KAW turbulence, a new regime of weak turbulence in which the effects of shearing by large scale motions and kinetic dissipation play an important role. The inclusion of the effect of nonlocal motions on the nonlinear energy cascade rate in the dissipation range, specifically the shearing by large-scale motions, is proposed to explain the nearly power-law energy spectra observed in the dissipation range of both kinetic numerical simulations and solar wind observations.

In turbulence phenomena, including the quantum turbulence in superfluids, an energy flux flows from large to small length scales, composing a cascade of energy. A universal characteristic of turbulent flows is the existence of a range of scales where the energy flux is scale-invariant: this interval of scales is often referred to as inertial region. This property is fundamental as, for instance, in turbulence of classical fluids it characterizes the behavior of statistical features such as spectra and structure functions. Here we show that also in decaying quantum turbulence generated in trapped Bose-Einstein condensates (BECs), intervals of momentum space where the energy flux is constant can be identified. Indeed, we present a procedure to measure the energy flux using both the energy spectrum and the continuity equation. A range of scales where the flux is constant is then determined employing two distinct protocols and in the same range, the momentum distribution measured is consistent with previous work. The successful identification of a region with constant flux in turbulent BECs is a manifestation of the universal character of turbulence in these quantum systems. These measurements pave the way for studies of energy conservation and dissipation in trapped atomic superfluids, and also analogies with the related processes that take place in ordinary fluids.

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.