Abstract
A macroscopic collective motion of a Bose–Einstein condensate (BEC) is commonly associated with phenomena such as superconductivity and superfluidity, often generalised by the term supercurrent. Another type of motion of a quantum condensate is second sound—a wave of condensate’s parameters. Recently, we reported on the decay of a BEC of magnons caused by a supercurrent outflow of the BEC from the locally heated area of a room temperature magnetic film. Here, we present the observation of a macroscopic BEC transport mechanism related to the excitation of second sound. The condensed magnons, being propelled out of the heated area, form compact humps of BEC density, which propagate many hundreds of micrometers in the form of distinct second sound—Bogoliubov waves. This discovery advances the physics of quasiparticles and allows for the application of related transport phenomena for low-loss data transfer in magnon spintronics devices.
Similar content being viewed by others
Introduction
Supercurrent is a macroscopic quantum phenomenon, which appears, when many bosons (real- or quasiparticles) being self-assembled in one quantum state with minimum energy and zero velocity—a Bose–Einstein condensate (BEC)1,2,3,4,5,6,7,8,9—move as a whole due to a phase gradient imposed on their joint wave function. This phenomenon, although mostly associated with resistant-free electric currents of Cooper pairs10 in superconductors and superfluidity of liquid Helium11,12 is, however, much more widespread13,14,15. It is experimentally confirmed in the quantum condensates of diluted ultracold gases16,17, of nuclear magnons in liquid 3He18,19,20, of polaritons in semiconductor microcavities21 and, recently, of electron magnons in room-temperature ferrimagnetic films22. Supercurrents being topologically confined often manifest themselves in a form of quantum vortices16,23,24.
Another form of motion of the quantum condensates is a second sound11,25,26,27,28,29—wave motions of the superfluid parameters, which can propagate in continuous media with an almost linear dispersion law in the long-wavelength limit. From a general viewpoint, this is classified as a specific example of Nambu–Goldstone (NG) modes that appear necessarily in media exhibiting spontaneous breakdown of continuous (for example, translational and/or rotational) symmetries30.
The term second sound stems from an analogy with the ordinary sound waves or the first sound—the wave motions of media density and mechanical momentum. The most well-known example is the second sound in the superfluid 4He (and in a BEC of cold trapped atoms)—the anti-phase oscillations of the BEC parameters (describing the superfluid component) and of the parameters of excited Bose-particles (describing the normal-fluid component). During second sound propagation, the total density of the supporting media (e.g. superfluid 4He) remains constant, ρHe = const11.
Some solid dielectrics at low temperatures support another type of second sound—the waves of the phonon density31,32,33,34, with phonons being quanta of the first sound. Here, the phonon density is an analogue of the normal-fluid components in the superfluids. Similarly to the case of 4He, during the propagation of this type of second sound, the crystal density ρcr remains constant at large scales.
In the BECs of quasiparticles (e.g. magnons in ferromagnetic dielectrics) one more type of sound may exist—wave motions of the BEC parameters, similar to the fluctuations of the 4He superfluid component parameters during the second sound propagation. In the simplest case of the complex order parameter ψ(r, t), these motions are the Bogoliubov waves of the amplitude and the phase of ψ(r, t). Unlike the densities of superfluid and normal components of 4He, changes in the density of quasiparticles are unrelated to the density of the supporting media. As a result, during propagation of the magnon Bogoliubov waves, the density of the ferromagnetic crystal remains constant, ρfm = const. Therefore, the Bogoliubov waves in the magnon BEC cannot be associated with the first sound and should be understood as a new type of second sound.
In this paper, by exploring the spatio-temporal dynamics of a magnon BEC prepared by microwave parametric pumping in a room-temperature single-crystal film of a magneto-dielectric material35,36,37, we present the experimental discovery and analytical description of this new type of second sound. Furthermore, we demonstrate a transition from the supercurrent-type to the second-sound-type motion of the magnon BEC. In our experiment, a magnon supercurrent, flowing out from a thermally induced magnetic inhomogeneity22, creates perturbations in the BEC density, which travel hundreds of micrometers through the thermally homogeneous film areas almost without changing their form. A detailed analysis of the propagation features of these perturbations in the framework of the Gross-Pitaevskii equation allows us to describe the observed phenomenon as a solitary NG-wave, consisting of the magnon Bogoliubov waves, spontaneously breaking translational symmetry of the magnon BEC. In the long-wavelength limit, realized in our experiment, the Bogoliubov waves have a linear dispersion law and, thus, can be considered as a magnon second sound potentially featuring viscosity-free propagation through the magnon condensate. We consider the discovery of the magnon sound that supports the BEC propagation outside of the area of a strong thermal gradient as the main achievement of this paper.
Results
Experiment
The experiment is carried out in a sample (see Methods) cut out from a single-crystal ferrimagnetic film of Yttrium Iron Garnet (YIG)38. This magneto-dielectric material possesses the lowest known magnetic damping in nature and is one of the favourite magnetic media for fundamental and applied studies in modern magnonics and spintronics39. The detection of magnon dynamics is performed by frequency- and time-resolved Brillouin Light Scattering (BLS) spectroscopy with wavevector sensitivity (see Methods).
Recently, using this method, we succeeded in revealing the presence of a magnon supercurrent by the analysis of an enhanced decay of the magnon BEC in a heated spot created by a probing laser beam22. In ref. 22 we used a single laser beam both for the heating of the film and for the BLS probing of the magnon BEC. By varying the duration and the power of the probing light we measured and compared time evolutions of the magnon BEC density both in cold and in hot focal laser spots. It was concluded that the only reasonable explanation for the difference between the observed BEC evolutions is based on the assumption of existence of the thermally driven magnon supercurrent flowing from the hot spot to the cold film. This experiment can be considered as indirect observation of the magnon supercurrent, without explicit analysis of its parameters.
In our current work, we significantly improved the experimental setup, used in ref. 22. Now it enables the detection of temporal changes in the spatial distribution of the magnon density caused by the magnon supercurrent. The main idea underlying the direct observation of the motion of the magnon BEC is to separate the area of supercurrent formation from the area in which it was observed. Two independent light sources are used, one for the heating and the other for the BLS probing of the magnon gas as is shown in Fig. 1a. Light from a powerful blue laser, which is completely absorbed by a 5.6-μm-thick YIG film, is used for the film heating. The green light of low power, which is able to penetrate the YIG film with moderate attenuation, is used for the probing. Both the blue laser source and the YIG sample are mounted together on a movable stage to ensure that the predefined position of the heated area on the YIG-film surface remains unchanged. The stage motion relative to the focal green spot allows for the probing of the magnon gas density in different points of the sample.
In order to achieve BEC40,41, magnons are injected into the spin system of the YIG film via parallel parametric pumping35,37, which is currently considered the most efficient technique for magnon excitation over a large wavevector range. The process can be described by the splitting of a photon of a pumping electromagnetic wave with a wavevector of nearly zero and a pumping frequency ωp into two magnons with opposite wavevectors ±q and a frequency of ωp/2 (see Fig. 2a). The strength of the bias magnetic field H of 1690 Oe is chosen to allow for the injection of the magnon pairs slightly above the ferromagnetic resonance frequency ωFMR, where the parallel pumping process is most efficient42,43. The injected quasi-particles thermalize by way of four-magnon scattering processes conserving both their number and the total energy44,45,46,47,48. Finally, when the total number of magnons reaches a threshold value, a magnon BEC forms40,41. Due to the spatially confined microwave excitation, the BEC is formed in an oblong YIG region located just above the pumping resonator. This quasi-one-dimensional region extends for 50 μm along and for 1 mm across the direction of the bias magnetic field H (see Fig. 1b). Due to the large anisotropy of the magnon spectrum (see Fig. 2b), the effective mass of the magnons, which is inversely proportional to the magnon dispersion coefficient D(qx, qy) = d2ω(q)/(2dq2), strongly depends on the orientation of the magnon wavevector q. The mass of the magnons with wavevectors placed perpendicular to the bias magnetic field is ~21 times smaller than for magnons having wavevectors along the field22. The magnon supercurrent is therefore expected to be ~21 times stronger across the biased magnetic field and thus, along a large extent of the condensate. Taking into consideration these facts, we currently restrict ourselves only to one-dimensional scanning, along the pumping resonator (Fig. 1b).
The results of the spatially-resolved measurements of the magnon density evolution at the bottom of the spin-wave spectrum are shown in Fig. 3 for different heating conditions. To understand the dynamics of the magnon BEC, we first establish how it behaves in a spatially uniform room-temperature profile. The reference measurement, performed without heating, shows that the spatial distribution of the magnon condensate along the microstrip resonator is uniform (see BEC time profiles shown by the blue curves in Fig. 3a–e and a time-space diagram in Fig. 3f). Similar to the results of ref. 41, the BLS signal, which is proportional to the magnon density, rises sharply after the microwave pumping pulse is switched off due to the intensification of the BEC formation process caused by the so-called supercooling mechanism of the magnon gas41. Due to the intrinsic magnon relaxation to the phonon bath, the density of the freely-evolving magnon BEC exponentially decreases and we observe a spatially uniform decay of the BLS signal.
We now focus on the temperature-gradient-dependent behaviour of the magnon condensate. The outcome of the experiment is strongly changed when an additional local heating by the blue laser light with power of 15, 30, 44, 81, or 116 mW is applied. The lowest heating power of 15 mW leads to the result expected from the BEC’s decay dynamics22: the initially spatially uniform BEC density distribution, created after the termination of the microwave pumping, becomes progressively spatially inhomogeneous due to the enhanced decay of the BEC density in the centre of the heated area (see the magenta curve in Fig. 3a and time-space diagram in Fig. 3g). According to the model proposed in ref. 22, this decay is caused by the thermally driven magnon supercurrent flowing to the colder film areas. The experimental data presented in Fig. 3g unambiguously corroborate this model: being pushed out from the centre of the hot spot in the directions indicated by two empty arrows, the BEC phase accumulates around the hot spot. This results in two spatially localized magnon humps. The position of these humps does not change over time. Thus, we can conclude that in this case the magnon supercurrents vanish in the close vicinity of the hot spot.
At higher heating power, the BEC behaviour changes drastically: after some time elapses, two magnon density humps are identified travelling outwards from the hot spot. The red curves in Fig. 3a–e and time-space diagrams in Fig. 3h–k exemplify the BEC dynamics measured for higher laser powers. Note that the heating is applied continuously, thus unlike in ref. 22 the condensate forms in a predefined stable and strong non-uniform temperature profile, which does not change in the process of the BEC formation and the development of the supercurrent magnon flow. As a result, all effects become more pronounced (cf. the relaxation rates of the blue, magenta and red curves in Fig. 3a). In particular, we observe a two and a half orders of magnitude larger decay of the magnon density in the BEC regime (for the maximal heating power of 116 mW, see the red curve in Fig. 3a). To compare, in refs. 22,48 the BEC has decayed over only about two decades during the same time interval. Under these heating conditions, the magnon density distribution becomes spatially inhomogeneous already during the condensate formation and the deep magnon density hole surrounded by the dense magnon humps (see Fig. 3h–k) rapidly appear due to the strong magnon outflow. Notably, these humps do not remain at the border of the hot spot but continue to move in the sample areas with rather uniform temperature landscapes. The directions of this motion is marked in Fig. 3h–k by the filled arrows. As is shown by the simulated temperature distributions (see Methods) represented by the red shadings around x = 0 in Fig. 3g–k, even the strongest heating by the blue laser is rather local and does not extend beyond ±50 μm from the centre of the focal spot. At the same time, the humps, visible in Fig. 3j, k, propagate to more than 400 μm. The maximum propagation distance observed in our experiment is 600 μm, being restricted by the right edge of the YIG sample. Moreover, one of the density humps is efficiently reflected from the nearby YIG-film edge on the left side of the sample (see Fig. 1b) and propagates hundreds of microns back towards the hot spot (see the corresponding pulses on green curves in Fig. 3b–d). Such a dynamics is common for all heating powers from 30 mW to 116 mW shown in Fig. 3h–k. Stronger laser power leads to a faster formation of more intense BEC humps but does not changes the initial velocities of these humps (see Fig. 4), which is ~400 ms−1.
The observed distant propagation of magnon density humps through the uniform BEC cannot be understood using the model assuming the accumulation of a BEC phase in a thermally created magnetic potential and a subsequent supercurrent22. Indeed, (i) these humps propagate in the film areas which are subject to rather weak temperature gradients; (ii) the reflected magnon humps move towards the hot spot and thus cannot be driven by the temperature gradient.
Bogoliubov waves and time evolution of magnon BEC
To discuss the observed time- and space-evolution of the total number of magnons near the lower end of the frequency spectrum Ntot = Nb + Nc, which comprises the gaseous bottom magnons Nb and condensed magnons Nc, we need to clarify the initial conditions for Nb and Nc at the moment of time t = +0, immediately after the termination of pumping power. The previous wavevector- and time-resolved BLS studies in ref. 41 revealed that strong electromagnetic pumping overheats the magnon gas and can thus prevent the condensation process, which develops rapidly after the pumping is switched off41,44,45,46. In addition, the experimentally observed magnon density is not affected by the local heating during pumping (see the dashed line at t = +0 in Fig. 3g-k). It indicates that in the entire pumping area no BEC exists under given experimental conditions before the pumping is switched off. Thus, we can assume that for t ≤ 0 Ntot ≈ Nb and \(N_{\mathrm{c}} \ll N_{\mathrm{b}}\) (or equal to zero).
Our understanding of what happens for \(t \,\,\gtrsim \,\,0\) in the presence of the hot spot, where the locally reduced saturation magnetization creates a frequency well, is at an infancy state. Essential progress in this direction requires purposeful experimental studies, which will then serve as a basis for a microscopic analytical theory. For the time being, the experiments can be reasonably well interpreted under a conjecture that fast thermalization of Nb for t > 0 creates a high-intensity condensate mainly (or only) in the hot spot. To find possible reasons for such localization we can speculate that the almost-space-homogeneous chemical potential touches the bottom of the frequency spectrum initially at the hot spot. If so, the magnon BECs with ±qmin (see Fig. 2a) will try to escape from this area due to their nonlinear repulsion49. As is evident in our analysis of the propagation data (see below), the repulsive frequency shift ΩNL is about one order of magnitude larger than the depth of the frequency well δω(T), induced by the heating (see Methods). Therefore, the magnon BEC can successfully escape the well, creating second sound pulses that propagate away from the hot spot.
These pulses can be formally represented by a Fourier series of monochromatic waves, which propagate with almost the same velocity due to the linearity of the dispersion law (3) in the small wavevector limit. This explains why the observed density humps propagate practically without changing their shape. Moreover, the travelling monochromatic waves should reflect simultaneously from the sample edge, creating a hump of the same form, propagating in the opposite direction. This is exactly what we observe as a reflected pulse on the green curves in Fig. 3e.
We would like to note that the massive transition of the gaseous bottom magnons to an efflux magnon condensate should significantly reduce their concentration Nb in the centre of the hot spot, and some diffusion of these quasiparticles towards the hot spot to replace the escaping BEC magnons is expected. One can guess that the effective diffusion coefficient is given by the dispersion coefficient \(D_x = \partial ^2\omega ({\mathbf{q}})/(2\partial q_x^2)\) at the minimum frequency, which under the conditions of our experiment is about 7.4 cm2 s−1. If so, the characteristic speed of the diffusive broadening of the magnon density hole of width L (about 100 μm or broader, see Fig. 3) will be approximately Dx/L ≈ 7.4 ms−1 or even smaller for the maximum laser power P = 116 mW. This is much smaller than the experimentally determined velocity ≈ 400 ms−1 of the propagating BEC humps. Therefore, the diffusion of the gaseous bottom magnons is negligible and cannot compensate for the hump run-away motion, leading to the appearance of the hole in the total magnon density Ntot, as is indeed observed in our experiment (see Fig. 3g).
The problem of propagation of the magnon density humps outside of the hot spot can be solved by analysis of the dynamic behaviour of the magnon BEC confined in a long but narrow area along the microstrip resonator (Fig. 1b) using a one-dimensional version of the Gross-Pitaevskii equation (GPE)50:
Here C(x, t) is the BEC wave function and W is the effective amplitude of the four-wave interaction, responsible for the nonlinear frequency shift. The first term in the right-hand side of Eq. (1) is the density of kinetic energy, and the second one is the density of potential energy in mean field approximation. In our case W is positive, corresponding to the repulsive interaction in the system of two interacting BECs49.
The stationary solution of GPE (1) has the form
where |C0|2 = Nc is the density of magnons in the BEC state. A small perturbation c(x, t) ∝ exp i[kx − Ω(k)t] of the stationary solution (2) is an NG-mode. Its dispersion law found by Bogoliubov has the form51:
where Ω and k are the angular frequency and the wavenumber of the perturbation, respectively. In the long-wavelength limit \(D_xk^2 \ll 2{\mathrm{\Omega }}_{{\mathrm{NL}}}\) this perturbation obeys a linear dispersion law
and is understood as a second sound propagating through the BEC with the velocity cs. At the same time, for large k, the BEC contribution may be neglected and the standard quadratic dispersion law applies near the minimum of the spectra: Ω(k) = Dxk2.
To check, which of these two limiting cases is closer to our experimental situation, we measured the width Δ of the magnon density humps at half of their height for different heating powers. These measurements enabled us to estimate the characteristic value of the second sound wavenumbers \(k \simeq {\mathrm{\pi }}/{\mathrm{\Delta }}\) (the first spatial Fourier harmonic), which give the main contributions to the propagating pulse, see Table 1. Additionally, for each laser power used for heating, we measured the position x(t) of the hump maximum as a function of the propagation time τ = t − t0, where t0 is some initial moment of time, for which the density hump is already well shaped and such analysis is therefore possible. In order to reveal a possible dependence of the sound velocity on the propagation time we fitted the experimental data by a second order polynomial function x(t) = csτ − δcsτ2/2τmax (see Fig. 4). Here τmax is the maximum propagation time, limited by the length of our sample, cs and δcs are fitting parameters. Using these fits we can estimate the initial velocity of the propagating hump cs (at t = t0, τ = 0) and its final velocity cs − δcs at τ = τmax. The results are given in Table 1.
Our main finding in Table 1 is that the initial velocity cs(t0) of the hump only weakly depends on its wavenumber k: it varies by about 8% as k changes twice. A natural explanation of this fact is based on the assumption that the propagating hump consists of Bogoliubov waves over the background of the magnon BEC with a linear dispersion law (4). If so, the velocity cs(t0) should indeed be k-independent. To check whether the long-wavelength limit assumption is really valid in our case, we estimate the value of the nonlinear frequency shift ΩNL (2) and the product Dxk2/2 by taking cs = 410 m s−1 and k ≈ 714 cm−1 from Table 1 for PL = 116 mW. The resulting estimates are \({\mathrm{\Omega }}_{{\mathrm{NL}}} \approx 2{\mathrm{\pi }} \cdot 18.1\,{\mathrm{MHz}}\) and Dxk2/2 ≈ 2π ⋅ 0.3 MHz. It is evident, that the long-wavelength limit is well satisfied, supporting the suggested second sound scenario.
The second argument in favor of the suggested scenario for the propagating humps is the time dependence of their velocity. According to Eqs. (2) and (4), the time dependence of the velocity of magnon second sound is \(c_{\mathrm{s}}(t) \propto \sqrt {{\mathrm{\Omega }}_{NL}(t)} \propto \sqrt {N_{\mathrm{c}}(t)}\). During hump propagation, the amplitude of the background condensate number Nc decays and consequently leads to a decrease in the sound velocity, as indeed observed in the experiment (see Table 1 and Fig. 4). Note, that a similar dynamics of a propagating localized pulse of the Bogoliubov wave was observed in a BEC of ultra-cold trapped atoms26. Interestingly, the sound velocity decreases much less than expected from the observed decay rate of the cold condensate. This may be explained by peculiarities of the propagation of an intensive hump, which gathers the condensed magnons in front of it and leaves the deep hole behind. Such an accumulation of condensed magnons in the hump may compensate its natural decay to a large extent. We are not yet in the position to discuss this complicated phenomenon in detail. These studies are on our agenda.
Using the estimated value of ΩNL, we can also estimate the coherence length ξ that determines the characteristic wavenumber k* = 1/ξ, at which the linear dispersion law (4) is changed to the quadratic one for k > k*. This length \(\xi = \sqrt {D_x/\Omega _{{\mathrm{NL}}}} \simeq 2.5\,{\rm{\mu}} {\mathrm{m}}\) is much smaller than the minimal hump width \({\mathrm{\Delta }} \simeq 44\,{\rm{\mu}} {\mathrm{m}}\), which gives the maximal second sound wavelength in our experiment (see Table 1). Therefore, the experimental conditions fully correspond to the linear part (4) of the dispersion law (3). This conclusion is consistent with our observations that the second sound pulse propagates as a whole, without dispersive spreading proportional to d2Ω(k)/dk2. On the other hand, the coherent length determines the sizes of possible topological singularities in the magnon BEC. For example, the diameter a0 of vortex cores in the magnon condensate measured under similar experimental conditions in ref. 24 is about 1 μm, which is comparable with our estimation. It is also instructive to compare ξ and a0 with the total propagation length about 600 μm. The order-of-magnitude difference means that we are really observing the distant transport of the magnon BEC. Note that in the experiments with a BEC of diluted cold atoms16 the ratio of the system size to the core radius is much smaller.
Discussion
The results presented in this article, addressing the spatially-resolved probing of the dynamics of a magnon BEC, provide direct evidence of supercurrent-related motion of the condensate outwards from the heated spot. The observed occurrence of the magnon BEC propagation outside of the temperature gradient is associated with the excitation of a new type of second sound: Bogoliubov waves in the magnon BEC condensate. The newly discovered magnon second sound differs from the second sound in dielectrics31,33,34, in which the phonons can be described in terms of their occupation numbers only, not taking into account their phases. It also differs from the second sound in superfluid 4He and in the BEC in diluted atomic systems, where the wave function describes the distribution of real atoms and not of the quasiparticles, as in our case.
The terms supercurrent and second sound reflect different manifestations of the very same phenomenon—the superfluidity, and both result from the solutions of the same GP equation. Note also that being the solution of GPE, the second sound wave package transports the densities of its integrals of motion: energy, momentum, number of magnetic quasiparticles and, associated with them, spin and magnetic moment. In a space-homogeneous coherent state of BEC in an ideal Bose gas, the entropy is equal to zero. In real physical systems the situation is more complicated and the role of entropy requires a detailed study in each particular case52, including the BEC of an interacting magnon gas. Our preliminary educated guess is that the entropy transport by the considered second sound (long Bogoliubov waves of magnon BEC) is negligibly small provided its wavelength is much larger than the coherence length. In this respect it differs from other types of second sounds, such as temperature waves in superfluid 4He and solids11,29,31,32,33, that do transport significant amount of entropy. It is important to note that according to Landau’s criterion51 Ω(k) > csk (see Eqs. 3 and 4), a BEC with repulsive interaction, as realised in our case, is an inviscid superfluid.
Clearly, the magnon second sound requires further detailed experimental and theoretical investigations. For example, one needs to account for the interactions of the bottom gaseous magnons with the BEC magnons, by formulating a model, analogous to the two-fluid model of superfluid helium11. In this model, another NG-type magnon second sound—the waves of the bottom gaseous magnons density—can be found53,54. Having said that, we should emphasize that from a practical point of view, the three observed phenomena: (i) the transition from the supercurrent-type to the second-sound-type propagation regime, (ii) the excitation of the second-sound-pulses, and (iii) the possibility of a long-distance spin-transport in the magnon BEC, has already paved a way for the application of magnon macroscopic quantum states for low-loss data transfer and information processing in perspective magnon spintronic devices55,56,57.
Methods
Sample
The Yttrium Iron Garnet (YIG, Y3Fe5O12)38 sample is 5 mm long and 1 mm wide. The single-crystal YIG film of 5.6 μm thickness has been grown in the (111) crystallographic plane on a Gadolinium Gallium Garnet (GGG, Gd3Ga5O12) substrate by liquid-phase epitaxy58 at a Department of Crystal Physics and Technology of the Scientific Research Company Carat (Lviv, Ukraine).
Experimental setup
A sketch of the experimental setup is shown in Fig. 1a. It consists of microwave and optical parts. The microwave part includes a microwave generator, which is used as a source for pumping pulses (pulse duration 2 μs, repetition time 1 ms, carrier frequency 13.6 GHz) followed by a power amplifier, which drives a microstrip resonator circuit with a peak power of 12.6 W. The 50 μm wide half-wavelength microstrip resonator, fabricated on top of an alumina substrate, is used to further increase the amplitude of the pumping microwave magnetic field and its spatial localization. The YIG sample is positioned on top of the resonator, in the area of the maximum microwave magnetic field.
The optical part (Fig. 1a) is used both for the probing of the magnon BEC by means of Brillouin light scattering spectroscopy as well as for the controlled local heating of the YIG sample. Its main parts are the probing green laser (single-mode, 532 nm wavelength), an acousto-optic modulator (AOM), blue laser for heating (multi-mode, 405 nm wavelength), and a tandem Fabry–Pérot interferometer (TFPI)59,60,61.
An excessive heating of the sample by the probing laser can lead to the formation of a magnon supercurrent22. In order to minimize the influence of the probing beam on the magnon dynamics, we utilize an acousto-optic modulator, which is used for chopping the probing beam into pulses, in order to reduce a parasitic heating of the sample. The pulsed probing beam (pulse duration 6 μs, peak power 9 mW) is then focused onto the sample surface into a focal spot of 20 μm in diameter. The scattered light is directed to the multipass TFPI for frequency selection with resolution of 100 MHz. A single photon counting avalanche diode detector is placed at the output of the interferometer. The output of the detector is then connected to a counter module synchronized with a sequence of microwave pulses62. Every time the detector registers a photon, this event is recorded to a database which collects the number of photons ensuring a time resolution of 250 ps. The frequency of the interferometer transmission is also recorded, thus providing frequency information for each detected photon.
Frequency- and wavevector-resolved BLS spectroscopy
Brillouin light scattering (BLS) can be understood as the diffraction of the probing light from a moving Bragg grating created by a magnon mode. Some portion of the scattered light, which is proportional to the number of magnons in this mode, is shifted in frequency by an amount equal to the frequency of the mode. In addition, the diffraction from the grating leads to a transfer of momentum during this process. The in-plane component of the wavevector qL of the incident light is inverted by a magnon mode, if the magnon wavenumber q satisfies the Bragg condition q = −2qL sin (θ), where θ is the angle of incidence. Simultaneously, the out-plane wavevector of the probing light is inverted due to its reflection from the metal microstrip underlying a semi-transparent YIG sample. By setting the angle θ equal to 9.7°, the selection of in-plane magnons with wavenumbers qy ≈ 4 ⋅ 104 rad cm−1 around one of the minima of magnon spectra (Fig. 2) is implemented in our setup63. The wavevector resolution was ±2 ⋅ 103 rad cm−1 allowing for a rather selective observation of the magnon dynamics around the energy minimum at qy = qmin. It enables us to avoid a possible spurious contribution of travelling magnon-phonon hybrid quasiparticles47 in our transport measurements.
Optical heating of the sample
In order to locally heat the sample in a spatial point separated from the probing spot an additional continuous 405 nm wavelength laser is used. The reasoning behind the choice of this laser wavelength is twofold. Firstly, it allows an efficient heating of the thin YIG film, since the absorption of the light in the YIG layer is inversely proportional to the wavelength of the light. Secondly, the chosen wavelength is well separated from the wavelength of the probing laser, and therefore can be rather easily filtered out to exclude any possible influence on the detection system. The heating beam focused into a focal spot of 80 μm in diameter provides a local increase of the sample temperature by 40 K at the maximum laser power of 116 mW. Spatial-resolved probing of the magnon dynamics is performed by the controlled displacement of the sample using a precise linear positioning stage. The entire stage is placed directly between the poles of the electromagnet, ensuring high field uniformity and stability. Since the heating laser is located on the same stage, the position of the heating spot is fixed relative to the sample and to the excitation circuit (and, thus, relative to the created magnon BEC). As a result, the described setup allows for space-resolved measurements of the magnon dynamics across the heated and cold areas of the sample as it is shown in Fig. 1. The spatial resolution (scanning step) is set to 10 μm, which corresponds to half the probing focal spot size.
Temperature and BEC frequency shift in the hot spot
The temperature profiles of the heated sample were determined by solving a 3D heat-transfer model of the experimental setup using the COMSOL Multiphysics software (https://www.comsol.com/). Hereby, the conventional heat conduction differential equation is solved under consideration of the boundary conditions applied to the model, the material parameters of the used materials and the applied heat source22. Heat is deposited exponentially along the film thickness and has a Gaussian distribution in the film plane reflecting the shape of the laser focal spot. The temperature distributions along the long side of the YIG film sample (see Fig. 1b) are shown by red shadings at the bottoms of Fig. 3g–k. The calculated temperature difference ΔT between the centre of the laser focal spot and the cold film for the maximal laser power PL of 116 mW is about 40 K. The corresponding frequency difference δω(T) caused by the temperature induced decrease in the saturation magnetization of an YIG film22,38,64 is about 2π ⋅ 4.6 MHz. This value is much smaller than the nonlinear frequency shift ΩNL in our experiment.
Data availability
The data that support the plots within this paper and other findings of this study are available from the corresponding author upon reasonable request.
References
Borovik-Romanov, A. S., Bun’kov, Yu. M., Dmitriev, V. V. & Mukharskii, Yu. M. Long-lived induction signal in superfluid 3He-B. JETP Lett. 40, 1033–1037 (1984).
Anderson, M. H., Ensher, J. R., Matthews, M. R., Wieman, C. E. & Cirnell, E. A. Observation of Bose-Einstein condensation in a dilute atomic vapor. Science 269, 198–201 (1995).
Davis, K. B. et al. Bose-Einstein condensation in a gas of sodium atoms. Phys. Rev. Lett. 75, 3969–3973 (1995).
Butov, L. V. et al. Stimulated scattering of indirect excitons in coupled quantum wells: signature of a degenerate Bose-gas of excitons. Phys. Rev. Lett. 86, 5608–5611 (2001).
Kasprzak, J. et al. Bose-Einstein condensation of exciton polaritons. Nature 443, 409–414 (2006).
Kalafati, Yu. D. & Safonov, V. L. Possibility of Bose condensation of magnons excited by incoherent pump. JETP Lett. 50, 149 (1989).
Demokritov, S. O. et al. Bose-Einstein condensation of quasi-equilibrium magnons at room temperature under pumping. Nature 443, 430–433 (2006).
Balili, R., Hartwell, V., Snoke, D., Pfeiffer, L. & West, K. Bose-Einstein condensation of microcavity polaritons in a trap. Science 316, 1007–1010 (2007).
Klaers, J., Schmitt, J., Vewinger, F. & Weitz, M. Bose-Einstein condensation of photons in an optical microcavity. Nature 468, 545–548 (2010).
Bardeen, J., Cooper, L. N. & Schrieffer, J. R. Microscopic theory of superconductivity. Phys. Rev. 106, 162–164 (1957).
Landau, L. D. Theory of the superfluidity of Helium II. Phys. Rev. 60, 356–358 (1941).
Leggett, A. J. Nobel lecture: superfluid 3He: the early days as seen by the theorist. Rev. Mod. Phys. 76, 999 (2004).
Page, D., Lattimer, J. M., Prakash, M. & Steiner, A. W. in Novel Superfluids, Vol. 2 (eds Bennemann, K.H. & Ketterson, J.B.) 505–579 (Oxford University Press, Oxford, 2014).
Bunkov, Yu. M. & Safonov, V. L. Magnon condensation and spin superfluidity. J. Magn. Magn. Mater. 452, 30 (2018).
Sonin, E. B. Spin currents and spin superfluidity. Adv. Phys. 59, 181–255 (2010).
Matthews, M. R. et al. Vortices in a Bose-Einstein condensate. Phys. Rev. Lett. 83, 2498–2501 (1999).
Raman, C. et al. Evidence for a critical velocity in a Bose-Einstein condensed gas. Phys. Rev. Lett. 83, 2502–2505 (1999).
Borovik-Romanov, A. S. et al. Observation of a spin-current analog of the Josephson effect. JETP Lett. 47, 478–482 (1988).
Bunkov, Yu. M. Spin superfluidity and coherent spin precession (Fritz London Memorial Prize Lecture). J. Phys. Condens. Matter 21, 164201 (2009).
Bunkov, Yu. M. et al. High-T c spin superfluidity in antiferromagnets. Phys. Rev. Lett. 108, 177002 (2012).
Amo, A. et al. Superfluidity of polaritons in semiconductor microcavities. Nat. Phys. 5, 805–810 (2009).
Bozhko, D. A. et al. Supercurrent in a room-temperature Bose-Einstein magnon condensate. Nat. Phys. 12, 1027 (2016).
Bunkov, Yu. M. & Volovik, G. E. Spin vortex in magnon BEC of superfluid 3He-B. Phys. C. 468, 609–612 (2008).
Nowik-Boltyk, P., Dzyapko, O., Demidov, V. E., Berloff, N. G. & Demokritov, S. O. Spatially non-uniform ground state and quantized vortices in a two-component Bose-Einstein condensate of magnons. Sci. Rep. 2, 1 (2012).
Pitaevskii, L. P. Second sound in solids. Usp. Fiz. Nauk 95, 139–144 (1968).
Andrews, M. R. et al. Propagation of sound in a Bose-Einstein condensate. Phys. Rev. Lett. 79, 553–556 (1997).
Steinhauer, J., Ozeri, R., Katz, N. & Davidson, N. Excitation spectrum of a Bose-Einstein condensate. Phys. Rev. Lett. 88, 120407 (2002).
Utsunomiya, S. et al. Observation of Bogoliubov excitations in exciton-polariton condensates. Nat. Phys. 4, 700–705 (2008).
Safonov, V. L. Nonequilibrium Magnons: Theory, Experiment, and Applications. (Wiley-VCH, New York, 2013).
Goldstone, J. Field theories with superconductor solutions. Nuovo Cim. 19, 154 (1961).
Ackerman, C. C., Bertman, B., Fairbank, H. A. & Guyer, R. A. Second sound in solid helium. Phys. Rev. Lett. 16, 789–791 (1966).
Nielsen, H. & Shklovskii, B. I. Heat transfer and second sound in dielectrics at large drift velocities. J. Exp. Theor. Phys. 56, 709–717 (1969).
Gurevich, V. L. Transport in Phonon Systems. (North-Holland Publ., Amsterdam, 1986).
Lee, S., Broido, D., Esfarjani, K. & Chen, G. Hydrodynamic phonon transport in suspended graphene. Nat. Comm. 6, 6290 (2015).
Rezende, S. M. & de Agnuiar, F. M. Spin-wave instabilities, auto-oscillations, and chaos in yttrium-iron-garnet. Proc. IEEE 78, 893–908 (1990).
L’vov, V. S. Wave Turbulence Under Parametric Excitation. (Springer-Verlag, Berlin, 1994).
Gurevich, A. G. & Melkov, G. A. Magnetization Oscillations and Waves. (CRC Press, New York, 1996).
Cherepanov, V., Kolokolov, I. & L’vov, V. The saga of YIG: spectra, thermodynamics, interaction and relaxation of magnons in a complex magnet. Phys. Rep. 229, 81–144 (1993).
Serga, A. A., Chumak, A. V. & Hillebrands, B. YIG magnonics. J. Phys. D. Appl. Phys. 43, 264002 (2010).
Rezende, S. M. Theory of coherence in Bose-Einstein condensation phenomena in a microwave-driven interacting magnon gas. Phys. Rev. B 79, 174411 (2009).
Serga, A. A. et al. Bose-Einstein condensation in an ultra-hot gas of pumped magnons. Nat. Commun. 5, 3452 (2014).
Neumann, T., Serga, A. A., Vasyuchka, V. I. & Hillebrands, B. Field-induced transition from parallel to perpendicular parametric pumping for a microstrip transducer. Appl. Phys. Lett. 94, 192502 (2009).
Serga, A. A. et al. Brillouin light scattering spectroscopy of parametrically excited dipole-exchange magnons. Phys. Rev. B 86, 134403 (2012).
Demidov, V. E. et al. Magnon kinetics and Bose-Einstein condensation studied in phase space. Phys. Rev. Lett. 101, 257201 (2008).
Clausen, P. et al. Stimulated thermalization of a parametrically driven magnon gas as a prerequisite for Bose-Einstein magnon condensation. Phys. Rev. B 91, 220402(R) (2015).
Bozhko, D. A. et al. Formation of Bose-Einstein magnon condensate via dipolar and exchange thermalization channels. Low. Temp. Phys. 41, 1024–1029 (2015).
Bozhko, D. A. et al. Bottleneck accumulation of hybrid magnetoelastic bosons. Phys. Rev. Lett. 118, 237201 (2017).
Kreil, A. J. E. et al. From kinetic instability to Bose-Einstein condensation and magnon supercurrents. Phys. Rev. Lett. 121, 077203 (2018).
Dzyapko, O. et al. Magnon-magnon interactions in a room-temperature magnonic Bose-Einstein condensate. Phys. Rev. B 96, 064438 (2017).
Rezende, S. M. Wave function of a microwave-driven Bose-Einstein magnon condensate. Phys. Rev. B 81, 020414(R) (2010).
Pitaevskii, L. P. & Stringari, S. Bose-Einstein Condensation. (Oxford University Press, New York, 2003).
Kim, M. B., Svidzinsky, A., Agarwal, G. S. & Scully, M. O. Entropy of the Bose-Einstein-condensate ground state: Correlation versus ground-state entropy. Phys. Rev. A. 97, 013605 (2018).
Nowik-Boltyk, P., Demidov, V. E., Demokritov, S. O., Tyberkevych V. & Slavin, A. Sound-like excitations in a dense gas of thermalized parametrically pumped magnons. In Proc. IEEE International Magnetics Conference. (Beijing, 2015).
Slavin, A. N. Magnonic second sound. In Proc. Sol-SkyMag 2018 conference. (San Sebastian, 2018).
Nakata, K., van Hoogdalem, K. A., Simon, P. & Loss, D. Josephson and persistent spin currents in Bose-Einstein condensates of magnons. Phys. Rev. B 90, 144419 (2014).
Chumak, A. V., Vasyuchka, V. I., Serga, A. A. & Hillebrands, B. Magnon spintronics. Nat. Phys. 11, 453–461 (2015).
Schneider, M. et al. Bose-Einstein condensation of quasi-particles by rapid cooling. Preprint at https://arxiv.org/abs/1612.07305v2 (2018).
Glass, H. L. Ferrite films for microwave and millimeter-wave devices. Proc. IEEE 76, 151–158 (1988).
Lindsay, S. M., Anderson, M. W. & Sandercock, J. R. Construction and alignment of a high performance multipass vernier tandem Fabry-Pérot interferometer. Rev. Sci. Instrum. 52, 1478–1486 (1981).
Mock, R., Hillebrands, B. & Sandercock, J. R. Construction and performance of a Brillouin scattering set-up using a triple-pass tandem Fabry-Pérot interferometer. J. Phys. E: Sci. Instrum. 20, 656–659 (1987).
Hillebrands, B. Progress in multipass tandem Fabry-Perot interferometry: I. A fully automated, easy to use, self-aligning spectrometer with increased stability and flexibility. Rev. Sci. Instrum. 70, 1589–1598 (1999).
Büttner, O. et al. Linear and nonlinear diffraction of dipolar spin waves in yttrium iron garnet films observed by space- and time-resolved Brillouin light scattering. Phys. Rev. B 61, 11576–11587 (2000).
Sandweg, C. W. et al. Wide-range wavevector selectivity of magnon gases in Brillouin light scattering spectroscopy. Rev. Sci. Instrum. 81, 073902 (2010).
Mihalceanu, L. et al. Temperature dependent relaxation of dipolar-exchange magnons in yttrium-iron-garnet films. Phys. Rev. B 97, 214405 (2018).
Acknowledgements
Financial support of the European Research Council within the Advanced Grant 694709 SuperMagnonics—Supercurrents of Magnon Condensates for Advanced Magnonics as well as financial support of the Deutsche Forschungsgemeinschaft (DFG) through the Collaborative Research Center SFB/TRR-49 Condensed Matter Systems with Variable Many-Body Interactions (project INST 161/544-3) and within the DFG project INST 248/178-1 is gratefully acknowledged.
Author information
Authors and Affiliations
Contributions
A.A.S. and B.H. contributed to the experimental idea, planned and supervised the project. D.A.B. and A.J.E.K. contributed to the experimental setup and carried out the experiments. A.P. and V.S.L. developed the theoretical model. H.Y.M-S. carried out the numerical analysis. A.J.E.K. performed temperature simulations. All authors analysed the experimental data and discussed the results.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Journal peer review information: Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work.
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Bozhko, D.A., Kreil, A.J.E., Musiienko-Shmarova, H.Y. et al. Bogoliubov waves and distant transport of magnon condensate at room temperature. Nat Commun 10, 2460 (2019). https://doi.org/10.1038/s41467-019-10118-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-019-10118-y
This article is cited by
-
Classical analog of qubit logic based on a magnon Bose–Einstein condensate
Communications Physics (2022)
-
Advances in coherent magnonics
Nature Reviews Materials (2021)
-
AC Josephson effect between two superfluid time crystals
Nature Materials (2021)
-
Magnon Quantization in the Magnetic Field Qradient
Applied Magnetic Resonance (2021)
-
Reflection-less width-modulated magnonic crystal
Communications Physics (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.