Abstract
The observation of Bose–Einstein condensation, in which particle interactions lead to a thermodynamic transition into a single, macroscopically populated coherent state, is a triumph of modern physics1,2,3,4,5. It is commonly assumed that this transition is a quantum process, relying on quantum statistics, but recent studies in wave turbulence theory have suggested that classical waves with random phases can condense in a formally identical manner6,7,8,9. In complete analogy with gas kinetics, particle velocities map to wavepacket k-vectors, collisions are mimicked by four-wave mixing, and entropy principles drive the system towards an equipartition of energy. Here, we use classical light in a self-defocusing photorefractive crystal to give the first observation of classical wave condensation, including the growth of a coherent state, the spectral redistribution towards equilibrium, and the formal reversibility of the interactions. The results confirm fundamental predictions of kinetic wave theory and hold relevance for a variety of fields, ranging from Bose–Einstein condensation to information transfer and imaging.
Similar content being viewed by others
Main
Weakly interacting bosons may exhibit, under extremely low temperatures, a Bose–Einstein transition1,2,3,4,5. This transition is characterized by a macroscopic occupation of the ground state. The quantum nature of the bosons involved is crucial for this process, and indeed Bose–Einstein condensation (BEC) is a wonderful manifestation of quantum mechanics at a macroscopic scale. On the other hand, a growing body of theoretical work has predicted that completely classical waves may undergo an analogous condensation process6,7,8,9. The requirements are a random ensemble of waves, so that statistical arguments apply, and a means of interaction between modes. Examples include sea waves stirred by wind, vibrations on elastic plates, interacting oscillators, and diffracting light that propagates in a nonlinear medium6,10,11. As these systems evolve, turbulent wave mixing leads to a self-organized redistribution of energy: an inverse cascade increases the ‘number of waves’ in the lowest allowed mode while a normal cascade transfers energy towards higher momenta. The process is ruled by the natural thermalization of a conservative Hamiltonian wave system. As with the collisions of particles, each wave interaction is formally reversible, yet entropy principles mandate that the ensemble evolves towards an equilibrium state of maximum disorder. In this way, a large-scale coherent structure grows and becomes immersed in a sea of small-scale fluctuations (‘uncondensed particles’), which store the information necessary for reversible evolution of the waves. Here, we directly observe this process of wave condensation, as well as its formal reversibility, by imaging classical light dynamics in a photorefractive crystal.
As in ordinary kinetic theory, wave dynamics can often be expressed as a hierarchy of interaction processes. Known as wave turbulence theory10,11,12,13,14,15, this approach focuses on the spectral distribution of modes as the field ψ(x,y,z) evolves. We apply the theory here to the nonlinear propagation of an optical beam, where ψ describes the envelope of the electric field and the dynamics occur in the x–y plane as the field propagates along z. To consider the interaction of modes, we concentrate on the (ensemble-averaged) power spectrum n(k,z), defined through , where is the Fourier transform of ψ(x,y,z) and * denotes complex conjugation. Formally, following this quantity results in an infinite hierarchy of moment equations. However, the hierarchy can be closed if the interactions are sufficiently weak (the ‘gas’ is sufficiently dilute), so that low-order collisions dominate. For example, for systems in which four-wave mixing is the dominant interaction, the evolution of n(k,z) follows the Boltzmann-like kinetic equation10,11,13,14,15
where g is the strength of the nonlinear interaction, is shorthand for n(ki,z), and W0 = σ δ(k1+k2−k3−k4)δ(k12+k22−k32−k42) is the collision term. In this latter term, σ = σ(k1,k2,k3,k4) is a scattering matrix, which in general depends on the complete set of k vectors, and the Dirac delta functions ensure the conservation of momentum and energy. For many systems, including and especially the homogeneous optical system considered below, σ reduces to a simple constant factor. This scalar form is therefore foundational, in the same sense that hard sphere or s-wave models form the basis of more complex theories16.
In complete analogy with Boltzmann’s equation, the kinetic wave equation is not reversible with respect to the propagation variable z. This happens despite the reversibility of the individual interactions, due to the truncation of higher-order terms and the averaging of the overall dynamics14,15. Formally, this is expressed by the H-theorem of entropy growth dS/dz≥0, where the nonequilibrium entropy reads . The equilibrium spectrum neq(k) is determined from the condition of zero entropy growth (maximum entropy11). One obtains the celebrated Rayleigh–Jeans distribution for classical nonlinear waves:
where k2/2kL may be regarded as the (transverse) kinetic energy of waves with wavelength 2π/kL and the parameters 1/T and −μ/T are Lagrange multipliers introduced to ensure conservation of energy and power (‘particle number’) . By analogy with thermodynamics, we call T the temperature and μ the chemical potential of the system. Interestingly, the one-to-one correspondence between {N,E} and {μ,T} implies that the field evolution can be described in the framework of the microcanonical statistical ensemble8,9 (in contrast with the canonical treatment using a thermal bath16). In what follows, we will use both sets of variables interchangeably.
The equilibrium distribution (2) is a Lorentzian spectrum, with statistics characterized by the correlation length . The correlation length diverges as μ →0, indicating that in an infinite system the k = 0 (‘plane-wave’) mode become macroscopically populated. In two dimensions, however, the requirements of finite temperature, density and energy mean that condensation can occur only in bounded systems8,9,17; in this case, waves accumulate in the lowest available wavenumber k0≠0 as μ→μc = k02/2kL. The other ‘excited’ waves occupy an algebraic k−2 spectrum, given by the tails of the distribution (2), indicating that there is an equipartition of energy among the uncondensed modes. For the experiment, a relatively large beam size was chosen (3 mm diameter), so that beam spreading was negligible over the length of the crystal. The surface area S of the beam thus bounds the region of wave interactions, and the system evolves with statistically homogeneous fluctuations (consistent with a basic assumption of wave turbulence theory). For the circular beam, the fundamental mode is characterized by , where is the first zero of the Bessel function J0(r). More details are given in the Supplementary Information.
We realize the condensation of classical waves experimentally using the set-up shown in Fig. 1. A finite-size, 10 mW beam is projected onto two sequential spatial light modulators (SLMs). The SLMs create a random-phase field (taken from a Gaussian (Maxwell–Boltzmann) distribution) with a user-defined correlation length lc. Nonlinear wave action occurs in a 5×5×10 mm SBN:75 (Sr0.75Ba0.25Nb2O6) photorefractive crystal with a self-defocusing (repulsive) coupling strength g controlled by applying a voltage across the c axis18. Output from the crystal is then imaged directly in both position (x) and momentum (k) space.
The condensation of classical waves is shown in Figs 2, 3, 4. Figure 2 shows experimental pictures when the nonlinear coupling strength is fixed at −500 V and the energy per particle E/N is varied. As the initial kinetic energy decreases, waves accumulate in the mode with the largest available spatial scale. Considering that the initial input beam was centred around k = 0, all modes with spatial frequencies belong to the condensate mode. Using the base width of the spectral peaks in Fig. 2d,e as the characteristic wavenumber scale, the measured value of k0 = 2×10−3 μm−1 matches very well with the theoretical prediction. The approach to this value as condensation progresses is seen more clearly by considering the corresponding chemical potential, plotted in Fig. 3a.
The amount of light in the condensate can be determined by a mode-counting argument. To do this, we limit the wavenumber range to , where the ultraviolet cut-off kD corresponds to the Debye length in the photorefractive crystal (∼2 μm). Following the usual interpretation of BEC, we separate the ground state contribution from that of the ‘excited states’16. In the limit μ→μc = k02/2kL, nkeq exhibits a divergence, showing that the fundamental mode k0 becomes macroscopically populated. The fraction of condensed power follows directly from equation (2) (a detailed derivation is given in the Supplementary Information) and reads , where N*is the number of modes and E0 = N k02/2kL refers to the lowest energy, that is the energy of the optical field when all quasi-particles occupy the ground state. Condensate fractions for μ = μc and μ≠μc are plotted in Fig. 3b. We note that the dynamics follows the ideal limit over most of the observed energy range, despite the finite size of the beam.
According to equation (2), the equilibrium state should be a Rayleigh–Jeans spectrum. Numerical and experimental proof of this is given in Fig. 3c,d respectively. The figures show explicitly the conservative nature of the spectral redistribution, with modes at both low and high wavenumbers growing at the expense of those in the middle of the range. This bidirectional energy transfer is required by the formal reversibility of the wave condensation process: as waves condense into the fundamental mode, information about their progress is recorded in the dynamics of the uncondensed, high- k modes8,9,19,20.
Nonlinear interactions are necessary for condensation to occur. Figure 4 shows experimental pictures when the initial kinetic energy is fixed (lc = 200 μm) while the nonlinear wave coupling (applied voltage) is increased. For linear propagation, the output spectrum (Fig. 4b) is statistically similar to the input spectrum (Fig. 4a). As the nonlinearity is increased (Fig. 4c–g), the spectrum concentrates in the fundamental k0≈0 mode. This mode reached a steady-state height at approximately −400 V, verifying that changing the interaction strength only changes the rate, not the asymptotic state, of the condensation process.
We emphasize that the observed spectral redistribution and condensation of light arises spontaneously from the natural thermalization of the conservative wave system. This contrasts with the familiar case of blackbody radiation, in which interaction with cavity walls prevents photon condensation from taking place5,21. The thermalization also contrasts with similar wave dynamics observed in other systems, such as optical mode coupling/locking in lasers22,23 and self-focusing turbulence in liquid crystals24, which are inherently dissipative systems driven far from equilibrium. In a similar way, optical wave thermalization through supercontinuum generation is affected by the dissipative Raman effect25. We note, in particular, that although wave turbulence theory is independent of the sign of the interactions11, in practice modulation instability dominates in the self-focusing case; the dynamics are controlled by the potential (rather than kinetic) energy of the system, the (inverse) momentum cascade is unidirectional, self-similar spectra are transient, and ‘condensation’ appears in the form of coherent soliton structures24,26. Localized intensity structures (zeros in the field) are present here from the random-phase initial conditions, but these did not participate actively in the dynamics. We did not observe waveguide filtering from each zero27 or an increase in their number during propagation, for example from vortex unbinding in a Berezinskii–Kosterlitz–Thouless transition16.
For long propagation distances (analogous to long evolution times), dynamics should consist of fluctuations around the equilibrium state. Unfortunately, the finite length of the crystal inhibits a direct measurement of this. Indeed, all nonlinear propagation experiments are limited by the finite extent of their media. To get around this, we introduce a novel feedback mechanism which takes advantage of the digital nature of our input/output devices (Fig. 1): the output field is recorded in the CCD (charge-coupled device) cameras and then used as the input field of the SLMs. This field recycling creates an effectively infinite crystal, subject to discretization errors, band-pass filtering, and bit-depth limitations. After a second pass through the crystal, the peak of the condensate remains within 10% of the single-pass value (with a correlation coefficient of 0.95 between the two spectral profiles). These results, in combination with longer-propagation simulations (Fig. 3c), confirm that the optical field has reached a (quasi-) equilibrium state.
The feedback method enables a scenario not possible in classical gas dynamics: selective and controlled adjustment of each particle’s momentum. To show reversibility, we change the sign of all the particle velocities; in the wave case, this is achieved by phase conjugation. As shown in Fig. 5b, the SLM can act as an ideal ‘Maxwell demon,’28 reversing the flow of condensation and recovering the initial ‘thermal’ cloud (with a correlation coefficient of 0.88). This experimentally proves that although the ensemble dynamics described by equation (1) always evolves towards maximum entropy, the individual wave collisions are formally reversible (conservative); as wavepackets collapse, the uncondensed ones retain all the information necessary to reconstruct the dynamical history.
The dynamics has interesting consequences when the waves have signal content. For example, recent methods of imaging using spatial nonlinearity have shown that wave mixing can greatly increase signal coherence29 and resolution19 while the information content of the system remains constant. More generally, information capacity is defined by entropy, yet the ability to reach capacity and the corresponding transmission dynamics are typically treated as separate issues28. Here, we have demonstrated that nonlinear interactions link the parameters of propagating beams, with entropy both constraining and driving the thermodynamics of classical waves.
Methods
The initial condition consists of a random-phase beam with Gaussian statistics and a correlation length lc = 2π/Δk:
where is delta-correlated. This input approximates a Maxwell–Boltzmann distribution and, in particular, ensures that all initial odd moments are zero (necessary to maintain a zero mean field). Inside the nonlinear crystal, the slowly varying field ψ obeys the Schrödinger equation , where the photorefractive screening nonlinearity Δn = −(1/2)n03ri jE0(I/1+I) (ref. 18). In this expression, n0 = 2.3, ri jis the electro-optic coefficient, I = |ψ|2 is in units of background illumination, and E0 is generated by applying a bias across the crystal. At the maximum applied voltage of −500 V, the nonlinear index change is measured to be Δn/n0 = −6.4×10−4. These experimental parameters were used for the numerical results in Fig. 3c, obtained by direct simulation of the Schrödinger equation for a spatially confined speckle beam of radius 1.5 mm (using a numerical grid of 1,024×1,024 points and an average over 20 realizations).
The real and Fourier space intensity at the output face of the crystal are imaged using CCD cameras. Phase information is retrieved experimentally using a phase shifting algorithm30, where the reference beam is shifted using a piezo-actuated mirror in increments of λ0/4±5 nm. The range of the cameras is limited to 12 bits. In Fig. 3d, the output power spectrum was measured using a series of different intensity filters, re-scaled and summed together to produce the entire power spectrum over 4.5 decades of intensity (2.7 decades of wavenumber). However, the input power spectrum was Gaussian only within two decades of intensity, below which it deviates from a quadratic falloff. Despite this, the system dynamically self-adjusts to the asymptotic k−2 spectrum over the full recorded range.
References
Davis, K. B. et al. Bose–Einstein condensation in a gas of sodium atoms. Phys. Rev. Lett. 75, 3969–3973 (1995).
Bradley, C. C., Sackett, C. A., Tollett, J. J. & Hulet, R. G. Evidence of Bose–Einstein condensation in an atomic gas with attractive interactions. Phys. Rev. Lett. 75, 1687–1690 (1995).
Anderson, M. H., Ensher, J. R., Matthews, M. R., Wieman, C. E. & Cornell, E. A. Observation of Bose–Einstein condensation in dilute atomic vapor. Science 269, 198–201 (1995).
Jochim, S. et al. Bose–Einstein condensation of molecules. Science 302, 2101–2103 (2003).
Klaers, J., Schmitt, J., Vewinger, F. & Weitz, M. Bose–Einstein condensation of photons in an optical microcavity. Nature 468, 545–548 (2010).
Dyachenko, S., Newell, A. C., Pushkarev, A. & Zakharov, V. E. Optical turbulence: Weak turbulence, condensates and collapsing filaments in the nonlinear Schrödinger equation. Physica D 57, 96–160 (1992).
Pomeau, Y. Long time behavior of solutions of nonlinear classical field equations: The example of NLS defocusing. Physica D 61, 227–239 (1992).
Connaughton, C., Josserand, C., Picozzi, A., Pomeau, Y. & Rica, S. Condensation of classical nonlinear waves. Phys. Rev. Lett. 95, 263901 (2005).
During, G., Picozzi, A. & Rica, S. Breakdown of weak-turbulence and nonlinear wave condensation. Physica D 238, 1524–1549 (2009).
Newell, A. C., Nazarenko, S. & Biven, L. Wave turbulence and intermittency. Physica D 152–153, 520–550 (2001).
Zakharov, V., L’vov, V. & Falkovich, G. Kolmogorov Spectra of Turbulence I (Springer, 1992).
Hasselmann, K. On the non-linear energy transfer in a gravity-wave spectrum. J. Fluid Mech. 12, 481–500 (1962).
Zakharov, V. E. Weak turbulence in plasma without magnetic field. Sov. Phys. JETP 24, 455–459 (1967).
Benney, D. J. & Saffman, P. G. Nonlinear interactions of random waves in a dispersive medium. Proc. R. Soc. Lond. A 289, 301–320 (1966).
Benney, D. J. & Newell, A. C. Random wave closures. Stud. Appl. Math. 48, 29–53 (1969).
Pitaevskii, L. & Stringari, S. Bose–Einstein Condensation (Oxford Science Publications, 2003).
Ketterle, W. & van Druten, N. J. Bose–Einstein condensation of a finite number of particles trapped in one or three dimensions. Phys. Rev. A 54, 656–660 (1996).
Kukhtarev, N. et al. Holographic storage in electrooptic crystals. Ferroelectrics 22, 949–960 (1979).
Barsi, C., Wan, W. & Fleischer, J. W. Imaging through nonlinear media via digital holography. Nature Photon. 3, 211–215 (2009).
Bromberg, Y., Lahini, Y., Small, E. & Silberberg, Y. Hanbury Brown and Twiss interferometry with interacting photons. Nature Photon. 4, 721–726 (2010).
Klaers, J., Vewinger, F. & Weitz, M. Thermalization of a two-dimensional photonic gas in a ‘white wall’ photon box. Nature Phys. 6, 512–515 (2010).
Conti, C., Leonetti, M., Fratalocchi, A., Angelani, L. & Ruocco, G. Condensation in disordered lasers: theory, 3D+1 simulations and experiments. Phys. Rev. Lett. 101, 143901 (2008).
Weill, R., Fischer, B. & Gat, O. Light-mode condensation in actively-mode-locked lasers. Phys. Rev. Lett. 104, 173901 (2010).
Bortolozzo, U., Laurie, J., Nazarenko, S. & Residori, S. Optical wave turbulence and the condensation of light. J. Opt. Soc. Am. B 26, 2280–2284 (2009).
Barviau, B. et al. Experimental signature of optical wave thermalization through supercontinuum generation in photonic crystal fibers. Opt. Express 17, 7392–7406 (2009).
Rumpf, B. & Newell, A. C. Coherent structures and entropy in constrained, modulationally unstable, nonintegrable systems. Phys. Rev. Lett. 87, 054102 (2001).
Coskun, T. H., Grandpierre, A. G., Christodoulides, D. N. & Segev, M. Coherence enhancement of spatially incoherent light beams through soliton interactions. Opt. Lett. 25, 826–828 (2000).
Brillouin, L. Science and Information Theory (Academic, 1962).
Dylov, D. V. & Fleischer, J. W. Nonlinear self-filtering of noisy images via dynamical stochastic resonance. Nature Photon. 4, 323–328 (2010).
Yamaguchi, I. & Zhang, T. Phase-shifting digital holography. Opt. Lett. 22, 1268–1270 (1997).
Acknowledgements
This work was supported by the Air Force Office of Scientific Research (USA), the Department of Energy (USA), Fondecyt N 1100289 (Chile), and the Grant ANR SYSCOM COSTUME (France).
Author information
Authors and Affiliations
Contributions
C.S. designed and performed experiments, performed simulations, analysed data and wrote the paper. S.J. and C.B. designed experiments and analysed data. S.R. and A.P. developed analytical tools, analysed data and wrote the paper. J.W.F. designed experiments, analysed data and wrote the paper. A.P. thanks R. Kaiser and G. Millot for fruitful discussions and valuable comments.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Information
Supplementary Information (PDF 610 kb)
Rights and permissions
About this article
Cite this article
Sun, C., Jia, S., Barsi, C. et al. Observation of the kinetic condensation of classical waves. Nature Phys 8, 470–474 (2012). https://doi.org/10.1038/nphys2278
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/nphys2278
This article is cited by
-
Non-equilibrium Bose–Einstein condensation in photonic systems
Nature Reviews Physics (2022)
-
Real-time multispeckle spectral-temporal measurement unveils the complexity of spatiotemporal solitons
Nature Communications (2021)
-
Ultra-light dark matter
The Astronomy and Astrophysics Review (2021)
-
Spatial Beam Self-Cleaning in Second-Harmonic Generation
Scientific Reports (2020)
-
Disordering a superfluid of light
Nature Photonics (2020)