Spatial separation of degenerate components of magnon Bose–Einstein condensate by using a local acceleration potential

Bose–Einstein condensation (BEC) of magnons is one of the few macroscopic quantum phenomena observable at room temperature. Due to the competition of the exchange and the magnetic dipole interactions, the minimum-energy magnon state is doubly degenerate and corresponds to two antiparallel non-zero wavevectors. Correspondingly, the room-temperature magnon BEC differs essentially from other condensates, since it takes place simultaneously at ± kmin. The degeneracy of BEC and interaction between its two components have significant impact on condensate properties. Phase locking of the two condensates causes formation of a standing wave of the condensate density and quantized vortices. Additionally, interaction between the two components is believed to be important for stabilization of the condensate with respect to a real-space collapse. Thus, the possibility to create a non-degenerate, single-component condensate is decisive for understanding of underlying physics of magnon BEC. Here, we experimentally demonstrate an approach, which allows one to accomplish this challenging task. We show that this can be achieved by using a separation of the two components of the degenerate condensate in the real space by applying a local pulsed magnetic field, which causes their motion in the opposite directions. Thus, after a certain delay, the two clouds corresponding to different components become well separated in the real space. We find that motion of the clouds can be described well based on the peculiarities of magnon dispersion characteristics. Additionally, we show that, during the motion, the condensate cloud harvests non-condensed magnons, which results in a partial compensation of condensate depletion.

www.nature.com/scientificreports/ with respect to a collapse in the real space 18,19 , the possibility of the superfluid behavior 16,20 , as well as the existence of the second sound 21 and Bogolubov waves 22 have been addressed. Another unique feature of the room-temperature magnon BEC in YIG, which makes its study particularly interesting and challenging, is its double degeneracy. Since magnon dispersion in YIG films is governed by the competition of two interactions-the dipole interaction dominating at small wavevectors and the exchange interaction dominating at large wavevectors,-the spectrum of magnons in magnetic films possess two energy minima corresponding to non-zero wavevectors ± k min 23 . Consequently, the magnon BEC is doubly degenerate: two condensate components are simultaneously formed in the two spectral minima at k BEC = ± k min . Below we refer to these components described by the wave functions proportional to exp(+ ik min z) and exp(− ik min z), as ( +) and (-) components, respectively. We emphasize that, although the phase velocities of magnons at ± k min can be as large, as several kilometers per second, the group velocity at these spectral points is zero. Therefore, the created condensate clouds do not move and occupy the same area in the real space.
The degeneracy has a significant impact on properties of magnon BEC and on interpretation of experimental results. For example, in 24 , it was experimentally shown that the two components of the condensate can be phase-locked to each other, resulting in formation of a standing wave of the total condensate density in the real space. From the theoretical point of view, on one hand, it is believed that the interaction of the two components is important for spatial stability of the entire condensate 18 . On the other hand, the same interaction is predicted to cause symmetry breaking resulting in different densities of the two condensate components 15 . In order to address these phenomena in detail, it is necessary to elaborate an approach, which would allow one to create single-component magnon BEC experimentally. At first glance, the most straightforward approach could be based on utilization of magnon injection mechanisms breaking the symmetry between +k and −k. However, all magnon pumping mechanisms used up to now are believed to satisfy intrinsically this symmetry, making the direct creation of single-component magnon BEC very challenging.
Here, we propose and experimentally test a different approach to the solution of this long-standing problem. Instead of the direct creation of single-component BEC, we implement a spatial separation of the two components of the condensate corresponding to ± k min . We first create a confined degenerate condensate cloud and then apply a pulse of localized magnetic field, which causes the two components of the condensate to move in the opposite directions. Thus, after a certain delay, the two clouds corresponding to different condensate components became well separated in the real space. We show that the observed velocities of the condensate clouds and their dependencies on the amplitude of the accelerating pulse are in good agreement with theoretical values obtained from calculations of magnon dispersion. Moreover, we find experimentally that the observed separation process conserves the total number of magnons in the split clouds.

Results
The schematics of our experiment is illustrated in Fig. 1a. The magnon condensate is created at room temperature in the yttrium iron garnet film with the thickness of 5.1 μm and the lateral dimensions of 2 × 10 mm 2 by using microwave parametric pumping. The pumping field is produced by a dielectric resonator with the resonant frequency of f MW = 9 GHz. The pumping injects primary magnons at the frequency f p = f MW /2 = 4.5 GHz, which thermalize and create BEC in the lowest-energy spectral states at k BEC = ± k min . The frequency of BEC is mainly determined by the applied static magnetic field H 0 and is equal to 3.5 GHz at H 0 = 1.2 kOe. As discussed in detail below, in order to separate condensate clouds, they are put in motion by applying an additional localized field ΔH(z) produced by an electric current I flowing through a control line with the width of 10 µm and the thickness of 400 nm, adjacent to the surface of YIG (Fig. 1a). We note that a non-uniform static magnetic field was also demonstrated to enable propagation of magnon BEC in 3 He 25 . Due to the shielding of the microwave field of the resonator by the metallic control line, the maximum pumping efficiency is reached in the area above the line. The microwave power applied to the resonator is chosen in such a way that the BEC threshold is exceeded in the vicinity of the control line only. Thus, we create a BEC cloud with the width of about 50 μm along the z-axis. In the absence of the acceleration field, the cloud remains centered at the control line.
We study spatial and temporal evolution the condensate density using micro-focus Brillouin light scattering (BLS) technique 26,27 with selectivity with respect to the wavevector of the scattering magnons 12,28,29 . As depicted in Fig. 1a, due to the conservation of the linear momentum, photons of the incident probing light are scattered by the ( +) and (−) component of the condensate at different angles. Correspondingly, by selecting for the analysis one or the other part of the scattered beam (shown by different colors in Fig. 1a), one can address one or the other component. Figure 1b,c illustrate the idea of the experiment. In order to put the condensate clouds in motion, we apply a pulse of additional localized magnetic field ΔH(z) with the rise time 5-10 ns. As a result, the initial condensate cloud is brought to the top of a potential hill defined by the magnetic field profile H 0 + �H(z) . Figure 1c shows dispersion curves for magnons at H 0 and at H 0 + H . While the former curve corresponds to the magnon spectrum far from the control line, the latter illustrates the magnon spectrum close to its center. The frequency of the condensate above the control line is higher than that of the lowest-energy magnon state in the surrounding film. Experiencing influence of the energy gradient caused by �H(z) , the condensate cloud starts to move in the real space. As a result, according to the energy conservation law, the condensate evolves into states with the absolute value of the wavenumber, which is larger (exchange-dominated magnons) and/or smaller (dipole-dominated magnons) than k BEC , as indicated in Fig. 1c by the dashed lines. Thus, two fractions (dipole-and exchangedominated) are formed. Let us first consider the evolution of the ( +) component of the condensate. It is obvious from Fig. 1c that, since the group velocity is determined by the slope of the dispersion curve, v gr = 2πdf/dk, for a non-monotonic dispersion spectrum f(k), two fractions of the ( +) component move in the opposite directions. Moreover, since, even in the vicinity of k BEC , f(k) is strongly non-parabolic, the cloud of dipole-dominated Scientific RepoRtS | (2020) 10:14881 | https://doi.org/10.1038/s41598-020-71525-6 www.nature.com/scientificreports/ magnons, corresponding to k < k BEC , possesses larger velocity than that built by exchange-dominated magnons, k > k BEC . Thus, the ( +) component is split in the real space into two clouds: a fast cloud moving in the (-z)direction and a slow cloud moving into the (+ z)-direction. By analyzing the motion of the (-) component, one obtains a fully symmetrical picture: its fast fraction moves into the (+ z)-direction, whereas the slow one moves into (-z)-direction (see Fig. 1b). Using spatial and temporal resolution of the BLS technique, we examine the above concept by recording temporal and spatial profiles of the condensate density separately for the two components, as shown in Fig. 2. The data obtained for ( +) and (-) components are shown as color-coded space-time maps in Fig. 2a,c, respectively. Note here that, since the horizontal axis in Fig. 2a,c shows the propagation distance, steeper trajectories in the maps correspond to slower motion. Figure 2b,d illustrate spatial profiles of the density recorded at different time delays with respect to the start of the accelerating field pulse, as labelled. One can definitely see that, indeed, each component is split into two clouds moving in the opposite directions with different velocities. Moreover, the already mentioned symmetry connecting the direction of the wavevectors and that of the velocities is also clearly seen in Fig. 2. It is important to notice that the magnetic field induced by the current is strongly localized around the control line. Thus, the majority of the travel range shown in Fig. 2 is free from the field gradient, which allows us to consider the motion of condensate clouds as a free propagation. We also note that, at large delay times, magnons start to build a new low-intensity stationary condensate in the area of the control line. Since this process is slow and the newly build condensate is at rest, it has no influence on the moving clouds.

Discussion
To get deeper insight into the studied process, we plot in Fig. 3 Figure 3 clearly demonstrates that, within the experimental accuracy, the clouds move with constant velocities: slow clouds move along the corresponding wavevectors k BEC , whereas fast clouds move in the opposite direction, in agreement with the above discussion based on the magnon dispersion spectra shown in Fig. 1c. As also seen from Fig. 3, the trajectories of the clouds of ( +) and (−) component start at t = 0 at slightly different positions. The same can be seen, if one compares one-dimensional traces in Fig. 2b,d for t = 0: the center of the ( +) component at t = 0 is shifted in the positive direction, whereas the center of the (-) component is shifted in the negative one. We associate these small shifts with the motion of magnons in the real space during their thermalization and formation of the condensate. In fact, although the distribution of the pumping field is symmetric, non-zero group velocities  Figure 4 illustrates temporal evolution of the spatial width and of the number of magnons in the moving clouds. Figure 4a shows variation of the width for the clouds formed by fast (dipole-dominated) and slow (exchange-dominated) magnons. 0 = 50 μm indicates the initial width of the cloud before application of the splitting field pulse. As seen from Fig. 4a, the splitting causes constriction of both clouds: the width of the exchange-dominated cloud is smaller than 0 during the entire propagation interval, whereas that of the dipole-dominated cloud is smaller than 0 at t < 100 ns. Additionally, the moving clouds experience dispersive broadening. For dipole-dominated magnons, this broadening is very strong, whereas for exchange-dominated magnons it is much weaker. These behaviors are consistent with the fact that the dispersion coefficient D = 2πd 2f / dk 2 of the dipole-dominated part of the magnon spectrum is much larger than that for the exchange-dominated part: D dip = 3.9 cm 2 /(rad s), D ex = 0.24 cm 2 /(rad s) 18 . In agreement with this interpretation, the dashed straight lines representing the best linear fit of the experimental data with the ratio of two slopes equal to D dip /D ex = 16.2, reproduce the observed broadening very well.
More intriguing is the evolution of the total number of magnons in the clouds, N, characterized in Fig. 4b. The shown data are normalized integrals of recorded spatial profiles of the condensate density versus propagation www.nature.com/scientificreports/ time for clouds formed by dipole-dominated, N dip and exchange-dominated magnons, N ex , correspondingly. The data are normalized by the value of the integral for the initial condensate cloud existing before application of the splitting field pulse. Results of two types of experiments are presented in Fig. 4b. In the first case (open symbols), the microwave pumping was switched off at t = 0. In the second case (solid symbols), the pumping was applied continuously. Keeping in mind short thermalization time of magnons 10 , the former case corresponds to free propagation of condensate clouds, whereas the latter case corresponds to the propagation of condensate clouds interacting with thermalizing non-coherent magnons. Before discussing the temporal decay of magnon clouds, we would like to stress, that the data of Fig. 4 clearly indicate that the splitting process does not significantly reduce the total number of magnons in the condensate. Indeed, the sum N dip + N ex remains close to N tot for t < 50 ns. From the theoretical point of view, the main process changing the total number of magnons in the condensate cloud is the Cherenkov emission and absorption of condensate magnons by thermal magnons. However, this process is strongly suppressed under conditions of our experiment due to the small Mach angle of the process, which results in the small statistical weight of final states.  Fig. 4b), one observes a clear exponential decay of the number of magnons N ~ exp(−t/τ) (note the log-linear scale of the graph) for both exchange-dominated and dipole-dominated magnons, although with slightly different decay times: τ ex = 175 ns, τ dip = 230 ns. In the case of continuous pumping (solid symbols in Fig. 4b), the observed behaviors differ strongly. While the existence of non-condensed magnons practically does not influence the decay of the exchange-dominated condensate, N dip (t) changes dramatically. For t < 170 ns (propagation distance < 80 µm) one observes a rather slow decay with the decay time τ * = 290 ns. This indicates existence of a mechanism, which partially compensates the intrinsic decay of N dip . We attribute these behaviors to a harvesting of non-condensed magnons by the cloud. As mentioned above, the control line modifies pumping conditions around it, leading to a formation of a condensate cloud with the width of 0 =50 µm (situated in the area ± 25 µm around the center of the control line). Away from this area, the pumping strength is below the threshold of BEC. However, it is still large enough to create an essential number of primary, and, as a consequence, of intermediate non-coherent magnons. Non-zero group velocity of the thermalizing magnons further increases the spatial extent of this zone. Thus, for some time, the accelerated condensate cloud propagates through the area, where these magnons are present. To understand the consequences, one should keep in mind, that the probability for a boson to scatter into a given spectral state is proportional to the population of this state. Outside of the condensation area, the density of non-coherent magnons is not high enough to create a condensate. However, as soon as the condensate cloud enters this area, the non-coherent magnons can efficiently scatter into the condensate and, in this way, partially compensate the decay of N. After some propagation time, however, the decay of the cloud becomes even faster than in the case of free propagation. This indicates presence of additional interaction mechanisms, which result in the scattering of magnons out of the condensate. It is interesting to note, that, as follows from the experimental data, the exchange-and the dipole-dominated magnons are influenced differently by the interaction with the non-condensed magnons.
Finally, we discuss the effects of the amplitude of the accelerating magnetic field ΔH, on the velocities of the clouds. Performing an analysis of the condensate movement at different ΔH values in the same way, as demonstrated in Fig. 3, we obtain dependences of the velocities on ΔH (Fig. 5) for the dipole-dominated (circles) and the exchange-dominated (squares) cloud. The solid curves in Fig. 5 present results of calculations based on the magnon dispersion curves shown in Fig. 1c assuming that the velocities of the clouds are equal to the group velocities of magnons v gr = 2πdf/dk corresponding to the spectral states marked in Fig. 1c by circles. A good agreement between the results of calculations and the experimental data validates our interpretation of the observed motion of the condensate clouds as caused by the accelerating potential.

conclusion
We have experimentally demonstrated that doubly degenerate condensate of magnons can be split in the real space into two components (( +) and (-)), corresponding to the opposite wavevectors ± k BEC , by applying a spatially localized pulse of magnetic field, which puts the components of the initially degenerate condensate in motion in the opposite directions. We find that this process is additionally accompanied by a splitting of each condensate component into two fractions constituted by dipole-and exchange-dominated magnons moving with different velocities. Thus, the initially localized condensate splits into four moving clouds. During the propagation, the dipole-and exchange-dominated clouds evolve differently, as the broadening of the cloud and the decay of the total numbers of magnons in the cloud are concerned. The faster, dipole-dominated condensate was found www.nature.com/scientificreports/ to harvest non-condensed magnons, which partially compensates the decay of the total number of magnons in the cloud. Our findings open a route for the efficient manipulation of magnon Bose-Einstein condensates in the real space and contribute to the deep understanding of the physics of degenerate condensates in general.

Methods
test system. The magnon BEC was studied in a 5.1 µm-thick YIG film epitaxially grown on a Gadolinium Gallium Garnet (GGG) substrate. A strong microwave field necessary for magnon injection was created by using a dielectric microwave resonator with the resonant frequency fp = 9 GHz fed by a microwave signal with the power of 4.5 W, which is by a factor of about 10 4 larger than the power threshold of the parametric instability 0.2 mW. The resonator was made of microwave ceramics with the dielectric constant ε = 80 and had the height of 3.5 mm and the diameter of 5 mm. The microwave magnetic field created by the resonator was directed perpendicular to the resonator top surface. The device was placed into a uniform static magnetic field H 0 = 1.2 kOe. An additional spatially inhomogeneous magnetic field was created by using an Au control line (width 10 µm, thickness 400 nm, and length 8 mm) carrying dc electric current. The line was fabricated on top of a 0.5 mm thick sapphire substrate. The substrate was placed between the resonator and the YIG sample in such a way, that the conductor was directly attached to the surface of the YIG film (see Fig. 1a). An electrical current with the peak value of 100 mA was applied in the form of pulses with the duration of 1 μs, the repetition period of 10 μs, and the rise/fall times of 5-10 ns. Since the condensate clouds quickly leave the area of the localized field, the time of interaction between the field and the clouds was 10-50 ns, depending on the velocity of the cloud determined by the amplitude of the field pulse. Estimates show that, during this time, the Joule heating of the sample by the current in the control line does not exceed 0.5 K, which ensures negligible influence of the heating on the studied processes. The absence of the influence of the heating was additionally confirmed by varying the duration and the repetition period of the current pulses, which was proven to have no effect on the results.
Space, time-and k-resolved BLS measurements. All the measurements were performed at room temperature. We focused the single-frequency probing laser light with the wavelength of 532 nm onto the surface of the YIG film through the transparent substrate and analyzed the light inelastically scattered from magnons. The measured signal-the BLS intensity-is directly proportional to the magnon density. By varying the lateral position of the laser spot, we recorded spatial profiles of the magnon density. The integral of the BLS intensity over the cloud in the real space is proportional to the total number of magnons in the cloud. To analyze the spatio-temporal dynamics of the condensate, the BLS measurement was synchronized with the field pulse: the BLS signal was recorded as a function of the temporal delay with respect to the rising edge of the pulse. To achieve the wavevector selectivity, the scattered light beam was divided into two parts, as indicated by different colors in Fig. 1a. By obscuring one half of the scattered-light beam, we discriminated the BLS signal caused by one component of the condensate only.

Data availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.