
==== Front
J Low Temp Phys
J Low Temp Phys
Journal of Low Temperature Physics
0022-2291
1573-7357
Springer US New York

3051
10.1007/s10909-024-03051-y
Article
Spectroscopy of Oscillation Modes in Homogeneously Precessing Domain of Superfluid 3He-B
http://orcid.org/0000-0001-9960-9323
Zavjalov V. V. v.zavjalov@lancaster.ac.uk

123
http://orcid.org/0000-0001-9894-1180
Savin A. 12
http://orcid.org/0000-0002-8668-396X
Sergeicheva E. 12
http://orcid.org/0000-0002-8247-4108
Hakonen P. J. 12
1 https://ror.org/020hwjq30 grid.5373.2 0000 0001 0838 9418 Low Temperature Laboratory, Aalto University School of Science, Espoo, Finland
2 https://ror.org/020hwjq30 grid.5373.2 0000 0001 0838 9418 Department of Applied Physics, Aalto University School of Science, Espoo, Finland
3 https://ror.org/04f2nsd36 grid.9835.7 0000 0000 8190 6402 Present Address: Department of Physics, Lancaster University, Lancaster, UK
17 2 2024
17 2 2024
2024
214 5-6 385398
20 10 2023
14 1 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ 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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/.
We study the homogeneously precessing domain (HPD) in superfluid 3He-B in a regular continuous-wave nuclear magnetic resonance (CW NMR) experiment. Using Fourier analysis of CW NMR time traces, we identify several oscillation modes with frequency monotonically increasing with the frequency shift of the HPD. Some of these modes are localized near the cell walls, while others are localized in bulk liquid and can be interpreted as oscillations of ϑ-solitons. We also observe chaotic motion of the HPD in a certain range of temperatures and frequency shifts.

http://dx.doi.org/10.13039/501100002341 Academy of Finland 341913 341913 341913 341913 Zavjalov V. V. Savin A. Sergeicheva E. Hakonen P. J. http://dx.doi.org/10.13039/501100000781 European Research Council 670743 670743 670743 670743 Zavjalov V. V. Savin A. Sergeicheva E. Hakonen P. J. http://dx.doi.org/10.13039/100010661 Horizon 2020 Framework Programme 824109 824109 824109 824109 Zavjalov V. V. Savin A. Sergeicheva E. Hakonen P. J. issue-copyright-statement© Springer Science+Business Media, LLC, part of Springer Nature 2024
==== Body
pmcIntroduction

Nuclear magnetic resonance (NMR) in superfluid 3He has played an important role in studies of its properties since the discovery of superfluid phases in 1972 [1]. The condensate in superfluid 3He is formed by p-wave pairing of fermions (pairs having spin 1 and orbital angular momentum 1), which allows the order parameter of the system to be written as 3×3 complex matrix. A few superfluid phases with distinct broken symmetries have been observed [2]; the most studied phases are the A phase and the B phase [3]. In the B phase, the order parameter structure involves an arbitrary 3D rotation matrix which sets the mutual orientation of spin and orbital spaces. The degrees of freedom of this matrix lead to three possible spin-wave modes, strongly affected by spin-orbit interaction. As a result, many unique phenomena can be observed in NMR experiments, including signatures of non-uniform spatial distribution and topological defects in the order parameter field [4–6], longitudinal NMR [7], and homogeneously precessing domain (HPD) [10]. Fundamental connections to physics in other areas have also been established [8, 9].

The HPD was discovered in 1985 [10, 11]. Since then, it has been used as a convenient probe in experiments with spin supercurrents [12] and vortices [13]. Various oscillation modes of HPD have been studied, including uniform rotational oscillations around magnetic field [14, 15], corresponding non-uniform oscillations [16], and oscillations of the HPD surface [17–19].

One of possible 2D topological defects in 3He-B is the so-called ϑ-soliton. It appears between two regions where the order parameter rotation matrix stays in different configurations at the minimum of spin-orbit interaction energy. In the regime of linear NMR, the analysis of the ϑ-soliton is quite straightforward [5]. The solitons can also be present in the HPD state [20], but their structure is more involved. ϑ-solitons in HPD have been observed experimentally in conjunction with spin-mass vortices in Ref. [21]. Recently, dynamics of a ϑ-soliton in the HPD has been considered in numerical simulations [22]. In this work, we experimentally observe spatially localized HPD oscillation modes that can be compared with simulations and attributed to oscillations of the ϑ-soliton, which so far have remained unidentified in experiments.

Theoretical Background

In the B phase of superfluid 3He, broken symmetries of the order parameter can be described by a 3×3 matrix1 Aaj=Raj(n^,ϑ)eiφ,

where Raj is a rotation matrix with axis n^ and angle ϑ, and eiφ is a complex phase factor. Various spatial distributions of φ and Raj are possible, as well as a few types of topological defects. Spin dynamics of 3He-B is described by nonlinear Leggett equations [3, 23] for spin S and matrix Raj:2 S˙a=[S×γB]a+415ΩB2χγ2sinϑ(4cosϑ+1)na-∇kJak,R˙aj=eabcRcj(γ2χS-γB)b.

Here γ is the gyromagnetic ratio of 3He, B denotes the external magnetic field, χ is the magnetic susceptibility, eabc is a permutation tensor, Jak specifies the spin current which carries component a of spin in the direction k, and ΩB is the Leggett frequency, a measure of spin-orbit interaction strength in superfluid 3He. All relaxation terms are neglected for simplicity.

Homogeneously precessing domain is a solution of these nonlinear equations [11]. Consider a coordinate system having axis z^ along the magnetic field, and rotating around this direction with some frequency ω≥γB. Then, the HPD state is given by3 n=y^,cosϑ=-14-1516ω(ω-γB)ΩB2,S=ωχγ2x^sinϑ+z^cosϑ-ω-γBω.

Note that this solution can have an arbitrary orientation around the z^ axis because the system is symmetric, but in a regular continuous-wave nuclear magnetic resonance (CW NMR) experiment, the presence of a radio-frequency field along x^ axis stabilizes HPD in n‖y^ orientation in the rotating frame. The important parameter here is the frequency shift δω=ω-γB. Usually, the shift is small, i.e., angle ϑ is slightly bigger than the so-called Leggett angle, ϑL=cos-1(-1/4)≈104∘, and the spin is tilted by approximately the same angle in the direction of the axis x^. The fact that the deflection of the spin is connected with the precession frequency ω makes an HPD state stable in a non-uniform magnetic field. If some spatial gradient appears in the precession frequency, spin currents transfer magnetization and compensate the difference. As a result, homogeneous precession takes place in the whole volume where the HPD exists. In fact, the HPD behaves like a liquid which fills low-field parts of the experimental cell up to the level ω=γB. In the case of free precession, this level is determined by the total system energy which decreases because of relaxation. In a continuous-wave NMR experiment with a sufficient feed of energy, the level can be controlled by the pumping frequency or magnetic field.

The topology of the HPD state is more complicated than that of an equilibrium state of 3He-B because, in addition to the orbital order parameter distribution, one can have a non-trivial distribution of the spin. One possible 2D topological defect in this state is the ϑ-soliton [20]. The minimum of spin-orbit interaction energy in 3He-B is achieved at the Leggett angle ϑl=cos-1(-1/4), the soliton appears between two different energy minima, ϑL and 2π-ϑL. There are two characteristic length scales in the ϑ-soliton in the HPD. First, a small core region of size ξD≈10 μm, across which the angle ϑ changes; ξD is referred to as “the dipolar length” because its scale is governed by the dipolar energy in comparison to gradient terms. Second, outside the core region, the magnetization varies on the length scale ξω that is inversely proportional to the square root of the frequency shift [20]:4 ξω=c‖ω(ω-γB)

where c‖ is the spin-wave velocity. The typical range of frequency shifts in our experiments is 10–100 Hz. This corresponds to ξω=0.15–0.5 mm.

An analytical calculation of ϑ-solitons in HPD is quite an involved task. In Ref. [22], we have performed numerical simulations in one-dimensional geometry, to obtain the structure and dynamics of the ϑ-soliton in HPD. A few low-frequency modes were identified in the simulation, and their dependencies on temperature, frequency shift, and radio-frequency field were determined. We use the simulated signatures and their characteristics to identify the measured oscillation modes in this work.

Experiment

Our experiments were performed on a pulse-tube-based nuclear demagnetization cryostat with the minimum temperature of 0.2 mK [24]. The experimental chamber was made using epoxy Stycast-1266 (see Fig. 1). It has a cylindrical shape with an inner diameter of 7.8 mm and a length of 9 mm. The experimental volume is connected to the heat exchanger of the nuclear stage through a channel having a diameter of 1 mm in its most narrow section.Fig. 1 Experimental cell with NMR magnet (left), close-up of the NMR coils and the cylindrical volume (center), and schematics of the NMR spectrometer (right). Detailed description of the setup is given in the text

The magnet system for a uniform static field includes a solenoid with a bore diameter of 36 mm, a gradient coil and a quadratic field coil. The magnet system, thermally anchored to the mixing chamber of the cryostat, is surrounded by a superconductive niobium shield. By contact to the nuclear stage, this shield provides mechanical stabilization, while keeping thermal isolation between the magnet system and the demagnetization cooling stage. The field homogeneity in the cell volume was measured using NMR linewidth in normal 3He. By adjusting the current in the gradient and quadratic field coils to their optimal values, it was possible to achieve a homogeneity of 1.5×10-4 over the sample volume. These normal-state NMR measurements also gave us a calibration of the gradient coil: 0.157 (T/m)/A.

Our NMR spectrometer consists of two pairs of RF coils, aligned perpendicular to each other, located symmetrically at a distance of 7 mm from the central axis of the cell. These RF coils with a diameter of 12 mm are made of 50 μm copper wire. The inductance of each coil pair is 55.5 μH, and the calculated value of the RF field in the center is 1.66 mT/A. One coil pair is used for applying NMR excitation from a signal generator, while the other one is embedded into an LC tank circuit with a resonant frequency of 1.124 MHz and Q=200. The signal from the receiving coils is amplified by a home-made HEMT amplifier [25]. Differential input of the amplifier is used to compensate for the background signal received away from the NMR resonance. By using the excitation and compensation voltages and the parameters of the electric circuit, one can calculate currents in both coil pairs and the total amplitude of RF-field. At an RF excitation of 1 Vpp (on the generator output), we need compensation voltage of 5.73 Vpp which corresponds to the RF-field rms amplitude BRF=297.6 nT. The NMR signal is recorded by a lock-in amplifier tuned at the excitation frequency and, in parallel, by a digital oscilloscope.

Our experiments were done at a pressure of 25.7 bar in a magnetic field of 34.67 mT (the NMR frequency equals to the resonant frequency of the tank circuit). Cooling capacity of the cryostat allowed us to measure a few hours in superfluid 3He in each demagnetization cycle. Measurements were done both during demagnetization to the lowest temperature and while warming up.

Temperature was measured using a SQUID-based noise thermometer attached to the nuclear stage. Since it is possible to observe the superfluid transition temperature Tc by means of NMR, we measured temperature difference between liquid helium and the nuclear stage at temperature Tc in the sample volume. The result can be described using a simple model with a thermalization time τ0: (THe-KTNS)/τ0=-KTNS˙, where THe and TNS are temperatures of the helium volume and the nuclear stage, respectively, and K≃1 is an adjustment factor for the noise thermometer calibration. Normally, the noise thermometer is calibrated at higher temperature against a thermometer at the mixing chamber; the factor K fixes the inaccuracy of this calibration. Measurements at different rates of cooling give us a time constant τ0=1540 s. A rough estimation of the thermalization time using the Kapitza resistance RK=900/T [K2m2/W] [26] and the normal 3He heat capacity CN=34.5T [J/K/mol] [27] for our amount of helium (approximately 1 mole) and sinter area (20 m2) give a quite similar value τ0=1550 s.

Below the superfluid transition temperature, the Kapitza resistance is expected to follow a 1/T law [28], on the basis of which the thermalization time can be estimated as τ=CB/CNτ0 where CN and CB denote the heat capacities of normal phase and B phase of 3He, respectively. By integrating the heat transfer model, we may estimate the temperature of 3He, THe, during the whole experiment as a function of time. The use of THe significantly reduces the observed hysteresis in temperature-dependent frequency shifts measured during cool-down and warm-up, which indicates that our simple model works well.

There is also another thermometric uncertainty: during HPD measurements with strong NMR excitation we have a significant (about 0.2 Tc) overheating of the 3He sample. It arises because of the narrow channel between the HPD volume and the heat exchanger volume. A correction of this overheating is implemented in a simple fashion (see Fig. 2A): we determine temperature-dependent frequencies of HPD oscillations (details are given later) and assume that they should extrapolate to zero at Tc. We do extrapolation using a smooth rational function and obtain a temperature shift, which is naively assumed to be temperature-independent. We do this separately for each sequence of measurements, because different measurement parameters (RF-field amplitude, gradient, range of magnetic field sweep) result in different overheating. The correction has been applied to all given values of T except those in Fig. 2A where the uncorrected value THe is used. The accuracy of such temperature corrections and our thermometry in general is quite poor, of the order of 0.1Tc. Nevertheless, this does not affect the main results of the work.

The HPD in our experiments was created in the usual continuous-wave method by sweeping the magnetic field in the presence of large enough RF-pumping at the resonance frequency of the LC circuit. The sweeping rate was always 0.807 μT/s (26.2 Hz/s in frequency units). The appearance of the HPD could easily be observed in the NMR signal recorded by the lock-in amplifier. In parallel, we recorded the same signal by an oscilloscope to observe the response in a wide range of frequencies. Oscillations of HPD can be seen as modulation of the main signal, with side bands separated from the main signal by the frequency of the oscillations. The oscillations were visible without any special excitation. We did measurements as a function of temperature, field gradient ∇B, and RF-field amplitude BRF. Because of overheating, we were able to create the HPD only in a small range of RF fields and only at small field gradients. Note that the field homogeneity of the NMR magnet 1.5×10-4 corresponds to field variations of about 5μT across the cell, while measurements were done at the absolute value of applied gradient from 0 to 20μT/cm. This means that the actual field profile and the process of HPD growth in our experiment are not well-defined. There could be a few field minima and maxima in the cell which lead to a complicated topology of the HPD surface, and this could be a possible source of ϑ-solitons. However, the main results in this work are obtained on a well-defined coherent HPD that is obtained in the regime where the HPD fills the whole cell.

We were unable to observe uniform oscillations of HPD [14, 15], presumably because we could not apply a rapid step change in B to excite them. It would have been useful to have a separate longitudinal coil around the cell to excite uniform oscillations and use their frequency for better temperature and BRF calibration.

An important experimental parameter is the frequency shift δω=ω-γB, the difference between pumping frequency ω which is always constant, and the Larmor frequency proportional to the magnetic field B, which changes during the experiments. The frequency shift varies across the sample because of the field inhomogeneity and the applied field gradient ∇B. The inhomogeneity limits the accuracy of calibration of B and the absolute value of the frequency shift δω. However, relative changes in the frequency shift are well-defined and accurate. Normally, we calibrate the magnetic field using NMR in liquid helium above Tc, and calculate the frequency shift using this calibration. We can only say that, at zero frequency shift condition, ω=γB is valid in some region inside the sample chamber. Such frequency shift values are presented in Figs. 2, 3, and 4. For high-temperature modes (displayed in Fig. 5), we use a more specific calibration of the magnetic field: we assume that the field at which the modes appear corresponds to zero frequency shift at the exact position of the soliton (details will be given in the Results section).Fig. 2 A Overheating correction for HPD measurements. We measure high-frequency mode gaps (will be explain later in the Results section) as a function of THe temperature, assume that they are proportional to each other, and extrapolate them to zero in order to find the actual Tc in the sample volume. Four different data sets are displayed, taken at different excitation levels. B Typical range of corrected temperatures and frequency shifts over which we detect HPD using BRF=357 nT (1.2 Vpp) and a field gradient of ∇B=-0.94μT/cm. The dashed red line is a theoretical curve for the maximum frequency shift using Eq. (6), with the value scaled down by a factor of 0.58 to match the measured data

Results

Figure 2B displays a typical range of temperatures and frequency shifts, over which we can observe the HPD state. During each measurement, the magnetic field is swept continuously upwards, the HPD appears when the frequency shift crosses zero in some part of the cell and exists until a maximum frequency shift is reached. We found that this maximum value can be estimated in the following way. Using Leggett’s equations with radio-frequency field BRF and including the Leggett–Takagi relaxation term [29], one can find the steady state that corresponds to the HPD. In the presence of Leggett–Takagi relaxation, the vector n is deflected from the direction given by Eq. (3), in such a way that5 nx=154δω2ΩB21γBRF1τ

where τ is the Leggett–Takagi relaxation time which is of the order of the mean Bogolyubov quasiparticle relaxation time (we use values calculated in Ref. [30]). We do not know the exact conditions at which the HPD becomes unstable, but we can say that it should happen before the vector n is rotated by 90 degrees, i.e., nx≈1. Due to these approximations, we cannot have an exact expression, but only an estimation for the maximum frequency shift:6 δωmax≈ΩBγBRFτ

In Fig. 2B, we display this estimation which matches the experimental results quite well when scaled by a factor of 0.58. Furthermore, the good agreement in the temperature dependence indicates that our thermometry corrections described in the Experiment section are realistic.Fig. 3 Typical recorded signals. Upper pictures display spectrograms obtained by Fourier transform on time-domain magnetization signals captured by a digital oscilloscope. The vertical axis is the frequency (difference from constant NMR frequency f0), horizontal axis is time (bottom axis) or frequency shift (upper axis), which is changing because of constant ramping of magnetic field, below them are signals measured with lock-in amplifier at the excitation frequency. A T=0.54Tc, HPD with chaotic region and “low-temperature modes” (side bands marked with a blue frame). B T=0.62Tc, “high-temperature modes” (marked with a blue frame). Data in A and B were taken at pressure 25.7 bar, BRF=357 nT, and ∇B=-0.94 μT/cm

We find an HPD state at temperatures above approximately 0.5Tc. Below this temperature the HPD is unstable because of catastrophic relaxation effects ([31, 32]). We distinguish two temperature ranges, below and above T≈0.59Tc, with completely different HPD oscillation modes. We call them “low-temperature” and “high-temperature” modes. At low temperatures and high frequency shifts we also observe chaotic motion of the HPD. This effect was studied numerically by Y. Bunkov in Ref. [33]. All these features can be seen in the spectrograms displayed in Fig. 3. The frequency spectrum is plotted as a function of time (bottom scale) and the frequency shift corresponding to the change of magnetic field in frequency units (upper scale). The left plot of Fig. 3A illustrates data recorded at 0.54Tc on the low-temperature modes and chaotic motion of the HPD; the right plot Fig. 3B presents data of the high-temperature modes recorded at 0.62Tc. All modes are seen as symmetric side bands of the NMR frequency ω. From the measured signals, we extract the mode frequencies Ω as a function of the frequency shift δω. We can see a clear dependence of the mode frequencies on the frequency shift in the presence of field inhomogeneity and different values of the applied field gradient. This indicates that the oscillation modes are localized.

Low-temperature modes are seen at T<0.59Tc. We observe two different modes and their harmonics at 2×, 3× the fundamental frequency (see the blue frame in Fig. 3A). The mode frequencies are proportional to the square root of the frequency shift and can be written as7 Ω2=SLT(δω-δω0),

where the slope S and offset δω0 are mode-dependent parameters as seen for the data in Fig. 4A. We can assume that the low-temperature modes do not have any intrinsic frequency gap and the offset δω0 is determined only by the location inside the cell in the presence of a field gradient. We found that δω0 does not depend on temperature nor on the amplitude of the RF field, but it does depend on the field gradient in a way as if the modes are localized at the upper and lower ends of the experimental cell (Fig. 4B). The temperature and RF-field dependence of the slope parameter SLT of mode frequencies is illustrated in Fig. 4C.Fig. 4 Characteristics of low-temperature modes. A Example of mode frequencies as a function of the frequency shift, extracted from two signals measured using two different field gradients: (a) ∇B=-13 μT/cm, (b) ∇B=-0.6 μT/cm. Both signals contain two modes which are indicated in green (mode 1) and magenta (mode 2). One can see a linear dependence of Ω2 on δω. Both signals have been measured with BRF=297 nT and temperature T=0.58Tc. B Mode offset position δω0 as a function of the field gradient. Solid lines are frequency shifts at cell ends calculated using the known field profile of the gradient coil and the cell size. Signals were measured with BRF=298 and 358 nT in the whole temperature range over which the low-temperature modes were visible. C Temperature dependence of mode slopes SLT. Different symbols are used for signals measured with different RF-field amplitudes: (∙), (∘), (∗) correspond to BRF=238,297,357 nT, respectively. ∇B=-0.6 μT/cm. The two colors are the same as in the other panels for the two different modes

High-temperature modes At higher temperature we observe a different set of oscillation modes (see the blue frame in Fig. 3B). In some rare cases they coexist with low-temperature modes, but usually there is a clear transition between these two regimes. There are many high-temperature modes, some are more common, while some appear only in a few measurements. Examples of two measurements are shown in Fig. 5, and different modes are marked by numbers 1...8.Fig. 5 Examples of high-temperature modes. BRF=297 nT, ∇B=-0.9μT/cm. The frequency shift is counted from the point where mode 2 appears, as explained in the text. A T=0.66Tc, B: T=0.61Tc. Different modes are marked by numbers 1...8. Black lines is a quadratic fit of modes 1 and 2 with a single frequency shift dependence and different gaps. Red dashed lines “num 1” and “num 2” are results of numerical simulations in Ref. [22]

Modes 1 and 2 are the most common ones, they appear in almost all measurements in the high-temperature regime. If we try to extrapolate their frequencies to zero, like we did with low-temperature modes, we find a big negative frequency shift which cannot be explained by field inhomogeneity. It means that the modes have some intrinsic frequency gap. We found that at the point where mode 2 becomes visible, a tiny but clear feature on the HPD signal can be observed ("oscillations start" line in Fig. 3B). We assume that this corresponds to the localization point of the modes, and the frequency shift should be measured from this value. We find that this point is always located in a random place inside the cell. (We can convert the magnetic field value to a position using the field gradient.) This allows us to determine the mode gaps Δ1 and Δ2. We see that the modes 1 and 2 have exactly the same frequency shift dependence, but there is a constant separation specified by the difference of Δ1 and Δ2 (see Fig. 5).

This information leads us to formulate our central conjecture: We reason that the high-temperature modes are oscillations of a ϑ-soliton localized somewhere in the cell volume. Low-frequency oscillations of the ϑ-soliton are calculated using 1D geometry as discussed in Ref. [22]. In our experiment, the soliton is comprised of a 3D circular membrane, in which case the radial part of oscillations should give an additional frequency gap8 Δn=anC/r,

where C is the wave velocity along the membrane, r is the radius, and an denote the zeros of the Bessel function of the first kind. For the two first radial modes a1=2.405, a2=5.520. Note that Ref. [22] discusses another type of oscillation of the circular membrane, namely one in which the whole soliton is moving. These modes should have much larger frequency, but we do not find them in the experiment. In Fig. 5, the dashed lines marked as “num 1” and “num 2” display results of the 1D numerical simulation, given by the approximative formulas Eqs. (9) and (10) of Ref. [22], shifted by Δ1. Apart from the gaps Δ1 and Δ2, we see that the measured modes 1 and 2 behave closely to the second calculated mode for the ϑ-soliton. The behavior observed for modes 7 or 8 could correspond to the first mode, but it’s hard to conclude with certainty owing to the small visibility range of the mode in our experiment.

Further basic information on the high-temperature modes 1 and 2 is presented in Fig. 6. The mode gaps Δ1 and Δ2 are plotted as a function of temperature in Fig. 6A. These data are the same as in Fig. 2A, but with the temperature correction applied. We do not observe any noticeable dependence of the gaps on BRF and ∇B. In Fig. 6B, the ratio Δ2/Δ1 is plotted for measurements where both modes exist. The observed ratio is pretty close to the theoretical ratio given by a2/a1=2.295. The values deduced for the wave velocity C=Δir/an are depicted in Fig. 6C. In Fig. 6D we compare the linear term in the mode frequency SHT=dΩdδω(δω=0) with numerical simulation. We find good agreement with the calculated second mode of soliton oscillations, with only a weak dependence on temperature and BRF. The fact that we can have a separate temperature-dependent gap and an almost temperature-independent soliton mode supports our choice of δω=0 for the frequency shift scale.Fig. 6 High-temperature modes 1 and 2 vs. T. Colors refer to different data sets measured in different cooling cycles using BRF=238,297,357 nT and ∇B=-0.9μT/cm. A Gaps Δ1 and Δ2 vs. temperature. The data are the same as in Fig. 2A but now with temperature correction. B Ratio Δ2/Δ1 for measurements where both modes exist. C Wave velocity C evaluated from Δ1 and Δ2 using Eq. (8). D Slope SHT of mode frequency Ω vs. frequency shift δω in comparison with numerical simulations (lines “num1” and “num2”). Dashed lines are calculated using BRF=238 nT and the solid ones using BRF=357 nT

We leave the detailed analysis of the other high-temperature modes outside the scope of this paper. Modes 3 and 4 (see Fig. 5A) were quite common, but they existed only in a small frequency shift range. They look like the next two modes of ϑ-soliton oscillations obtained in the numerical simulation, but at the moment we cannot confirm it. Modes 5 and 6 (see Fig. 5B) appeared in pairs with a frequency ratio of 2. They do not extrapolate to the same energy gap as modes 1 and 2, and probably they do not have a gap. Modes 7 and 8 were visible only in a few signals, which is unfortunate as they seem to bear resemblance with the first mode in our simulations. However, due to scarce information we could not analyze and classify them accurately.

Conclusions

In this work, we measured the response of the HPD state in an about 10 kHz band around the NMR frequency and could make interesting observations: new oscillation modes localized near the cell walls, additional modes in the bulk volume of the HPD that we attribute to oscillations of ϑ-solitons, and chaotic motion at high frequency shifts. We also estimated the maximum frequency shift at which the HPD state is stable and found that our experimental result matches our estimations based on Leggett–Takagi relaxation.

The richness of the observed soliton oscillation modes makes the most interesting finding of our experiments. Our work shows that, overall, the observed mode frequencies match with numerical simulations, but there are still more questions than answers: how these modes are excited, how the solitons are created, how to classify them and explain all the observed modes. Further investigation is certainly needed, in experiment, theory, and numerical simulations. It would be interesting to study how solitons are created at different magnetic field ramping rates, to isolate a single soliton by tuning experimental conditions without destroying the HPD, to populate certain modes by using additional RF-field excitation pulses, etc. It would be important to study how the creation of a soliton is affected by the magnetic field profile and sweeping rate. Altogether, the soliton dynamics would be a sophisticated and quite fundamental problem for further investigations.

Acknowledgements

This work was supported by the Academy of Finland projects 341913 (EFT) and 312295, 352926 (CoE, Quantum Technology Finland) as well as by ERC (Grant No. 670743). The research leading to these results has received funding from the European Union’s Horizon 2020 Research and Innovation Programme, under Grant Agreement no 824109. The experimental work benefited from the Aalto University OtaNano/LTL infrastructure. We would like to thank Vladimir Dmitriev, Yury Bunkov, and Grigory Volovik for useful discussions.

Data Availability

Our measured data are available on 10.5281/zenodo.8431264. It includes oscilloscope records, spectrograms, and extracted oscillation modes for about six hundred measurements.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Osheroff DD Richardson RC Lee DM Evidence for a new phase of solid He3 Phys. Rev. Lett. 1972 28 14 885 10.1103/PhysRevLett.28.885
D.D. Osheroff, R.C. Richardson, D.M. Lee, Evidence for a New Phase of Solid He3. Phys. Rev. Lett. 28(14), 885 (1972). 10.1103/PhysRevLett.28.885
2. Volovik GE Exotic Properties of Superfluid Helium 3 1992 Singapore World Scientific
G.E. Volovik, Exotic Properties of Superfluid Helium 3 (World Scientific, Singapore, 1992). 10.1142/1439
3. D. Vollhardt, P. Wölfle, The Superfluid Phases of Helium 3. Dover Books on Physics (Dover Publications, 2013)
4. Hakonen PJ Krusius M Salomaa MM Salmelin RH Simola JT Gongadze AD Vachnadze GE Kharadze GA NMR and axial magnetic field textures in stationary and rotating superfluid 3He-B J. Low Temp. Phys. 1989 76 3 225 10.1007/BF00681586
P.J. Hakonen, M. Krusius, M.M. Salomaa, R.H. Salmelin, J.T. Simola, A.D. Gongadze, G.E. Vachnadze, G.A. Kharadze, NMR and axial magnetic field textures in stationary and rotating superfluid He-B. J. Low Temp. Phys. 76(3), 225 (1989). 10.1007/BF00681586
5. Maki K Kumar P Magnetic solitons in superfluid 3He Phys. Rev. B 1976 14 118 10.1103/PhysRevB.14.118
K. Maki, P. Kumar, Magnetic solitons in superfluid . Phys. Rev. B 14, 118 (1976). 10.1103/PhysRevB.14.118
6. Maki K Kumar P n→ solitons (planarlike n textures) in superfluid 3He-B Phys. Rev. B 1977 16 4805 10.1103/PhysRevB.16.4805
K. Maki, P. Kumar, solitons (planarlike textures) in superfluid . Phys. Rev. B 16, 4805 (1977). 10.1103/PhysRevB.16.4805
7. Osheroff DD Brinkman WF Longitudinal resonance and domain effects in the a and b phases of liquid helium three Phys. Rev. Lett. 1974 32 584 10.1103/PhysRevLett.32.584
D.D. Osheroff, W.F. Brinkman, Longitudinal resonance and domain effects in the and phases of liquid helium three. Phys. Rev. Lett. 32, 584 (1974). 10.1103/PhysRevLett.32.584
8. Zavjalov VV Autti S Eltsov VB Heikkinen PJ Volovik GE Light Higgs channel of the resonant decay of magnon condensate in superfluid 3He-B Abstract Nat. Commun. 2016 7 1 10294 10.1038/ncomms10294 26743951
V.V. Zavjalov, S. Autti, V.B. Eltsov, P.J. Heikkinen, G.E. Volovik, Light Higgs channel of the resonant decay of magnon condensate in superfluid 3He-B Abstract. Nat. Commun. 7(1), 10294 (2016). 10.1038/ncomms1029426743951
9. Autti S Heikkinen PJ Mäkinen JT Volovik GE Zavjalov VV Eltsov VB AC Josephson effect between two superfluid time crystals Nat. Mater. 2021 20 2 171 174 10.1038/s41563-020-0780-y 32807922
S. Autti, P.J. Heikkinen, J.T. Mäkinen, G.E. Volovik, V.V. Zavjalov, V.B. Eltsov, AC Josephson effect between two superfluid time crystals. Nat. Mater. 20(2), 171–174 (2021). 10.1038/s41563-020-0780-y32807922
10. Borovik-Romanov A Bun’kov Y Dmitriev V Mukharskii Y Flachbart K Experimental study of separation of magnetization precession in 3He-B into two magnetic domains JETP 1985 61 1199
A. Borovik-Romanov, Y. Bun’kov, V. Dmitriev, Y. Mukharskii, K. Flachbart, Experimental study of separation of magnetization precession in He-B into two magnetic domains. JETP 61, 1199 (1985)
11. Fomin I Separation of magnetization precession in 3He-B into two magnetic domains. Theory JETP 1985 61 1207
I. Fomin, Separation of magnetization precession in He-B into two magnetic domains. Theory. JETP 61, 1207 (1985)
12. A.S. Borovik-Romanov, Y.M. Bunkov, V.V. Dmitriev, Y.M. Mukharsky, in Proceedings of 18th International Conference on Low Temperature Physics, Kyoto, vol. 26 Suppl. 3 (1987), vol. 26 Suppl. 3, pp. 175–176. 10.7567/JJAPS.26S3.175
13. Kondo Y Korhonen JS Krusius M Dmitriev VV Mukharsky YM Sonin EB Volovik GE Direct observation of the nonaxisymmetric vortex in superfluid 3He-B Phys. Rev. Lett. 1991 67 1 81 10.1103/PhysRevLett.67.81 10044057
Y. Kondo, J.S. Korhonen, M. Krusius, V.V. Dmitriev, Y.M. Mukharsky, E.B. Sonin, G.E. Volovik, Direct observation of the nonaxisymmetric vortex in superfluid He-B. Phys. Rev. Lett. 67(1), 81 (1991). 10.1103/PhysRevLett.67.8110044057
14. Dmitriev V Zavjalov V Zmeev D Spatially homogeneous oscillations of homogeneously precessing domain in 3He-B J. Low Temp. Phys. 2005 138 765 10.1007/s10909-005-2300-5
V. Dmitriev, V. Zavjalov, D. Zmeev, Spatially homogeneous oscillations of homogeneously precessing domain in He-B. J. Low Temp. Phys. 138, 765 (2005). 10.1007/s10909-005-2300-5
15. Kupka M Skyba P Bec of magnons in superfluid 3he-b and symmetry breaking fields Phys. Rev. B 2012 85 184529 10.1103/PhysRevB.85.184529
M. Kupka, P. Skyba, Bec of magnons in superfluid he- and symmetry breaking fields. Phys. Rev. B 85, 184529 (2012). 10.1103/PhysRevB.85.184529
16. Človečko M Gažo E Kupka M Skyba P New non-goldstone collective mode of BEC of magnons in superfluid 3He-B Phys. Rev. Lett. 2008 100 15 155301 10.1103/PhysRevLett.100.155301 18518119
M. Človečko, E. Gažo, M. Kupka, P. Skyba, New non-goldstone collective mode of BEC of magnons in superfluid He-B. Phys. Rev. Lett. 100(15), 155301 (2008). 10.1103/PhysRevLett.100.15530118518119
17. Bunkov Y Dmitriev V Mukharskii Y Low frequency oscillations of the homogeneously precessing domain in 3He-B Phys. B Condens. Matter. 1992 178 1–4 196 10.1016/0921-4526(92)90198-2
Y. Bunkov, V. Dmitriev, Y. Mukharskii, Low frequency oscillations of the homogeneously precessing domain in He-B. Phys. B Condens. Matter. 178(1–4), 196 (1992). 10.1016/0921-4526(92)90198-2
18. Lokner Ľ Feher A Kupka M Harakály R Scheibel R Bunkov YM Skyba P Surface oscillations of homogeneously precessing domain with axial symmetry Europhys. Lett. 1997 40 5 539 10.1209/epl/i1997-00501-8
Ľ Lokner, A. Feher, M. Kupka, R. Harakály, R. Scheibel, Y.M. Bunkov, P. Skyba, Surface oscillations of homogeneously precessing domain with axial symmetry. Europhys. Lett. 40(5), 539 (1997). 10.1209/epl/i1997-00501-8
19. Gažo E Kupka M Medeová M Skyba P Spin precession waves in superfluid 3he-b Phys. Rev. Lett. 2003 91 055301 10.1103/PhysRevLett.91.055301 12906602
E. Gažo, M. Kupka, M. Medeová, P. Skyba, Spin precession waves in superfluid he-b. Phys. Rev. Lett. 91, 055301 (2003). 10.1103/PhysRevLett.91.05530112906602
20. Misirpashaev TS Volovik G Topology of coherent precession in superfluid 3He-B JETP 1992 75 650
T.S. Misirpashaev, G. Volovik, Topology of coherent precession in superfluid He-B. JETP 75, 650 (1992)
21. Kondo Y Korhonen JS Krusius M Dmitriev VV Thuneberg EV Volovik GE Combined spin-mass vortex with soliton tail in superfluid 3He-B Phys. Rev. Lett. 1992 68 3331 10.1103/PhysRevLett.68.3331 10045675
Y. Kondo, J.S. Korhonen, M. Krusius, V.V. Dmitriev, E.V. Thuneberg, G.E. Volovik, Combined spin-mass vortex with soliton tail in superfluid . Phys. Rev. Lett. 68, 3331 (1992). 10.1103/PhysRevLett.68.333110045675
22. V.V. Zavjalov, Dynamics of -solitons in the HPD state of superfluid He-B (2022). 10.48550/arXiv.2212.12232
23. Leggett AJ A theoretical description of the new phases of liquid 3He Rev. Mod. Phys. 1975 47 2 331 10.1103/RevModPhys.47.331
A.J. Leggett, A theoretical description of the new phases of liquid He. Rev. Mod. Phys. 47(2), 331 (1975). 10.1103/RevModPhys.47.331
24. Todoshchenko I Kaikkonen JP Blaauwgeers R Hakonen PJ Savin A Dry demagnetization cryostat for sub-millikelvin helium experiments: refrigeration and thermometry Rev. Sci. Instrum. 2014 85 8 085106 10.1063/1.4891619 25173311
I. Todoshchenko, J.P. Kaikkonen, R. Blaauwgeers, P.J. Hakonen, A. Savin, Dry demagnetization cryostat for sub-millikelvin helium experiments: refrigeration and thermometry. Rev. Sci. Instrum. 85(8), 085106 (2014). 10.1063/1.489161925173311
25. Zavjalov VV Savin AM Hakonen PJ Cryogenic differential amplifier for NMR applications J. Low Temp. Phys. 2019 195 72 10.1007/s10909-018-02130-1 31073248
V.V. Zavjalov, A.M. Savin, P.J. Hakonen, Cryogenic differential amplifier for NMR applications. J. Low Temp. Phys. 195, 72 (2019). 10.1007/s10909-018-02130-131073248
26. Franco H Bossy J Godfrin H Properties of sintered silver powders and their application in heat exchangers at millikelvin temperatures Cryogenics 1984 24 477 10.1016/0011-2275(84)90006-7
H. Franco, J. Bossy, H. Godfrin, Properties of sintered silver powders and their application in heat exchangers at millikelvin temperatures. Cryogenics 24, 477 (1984). 10.1016/0011-2275(84)90006-7
27. Greywall DS Specific heat of normal liquid 3He Phys. Rev. B 1983 27 2747 10.1103/PhysRevB.27.2747
D.S. Greywall, Specific heat of normal liquid He. Phys. Rev. B 27, 2747 (1983). 10.1103/PhysRevB.27.2747
28. Ahonen AI Lounasmaa OV Veuro MC Measurements of the thermal boundary resistance between 3He and silver from 0.4 to 10 mk J. Phys. Colloques 1978 39 265 10.1051/jphyscol:19786117
A.I. Ahonen, O.V. Lounasmaa, M.C. Veuro, Measurements of the thermal boundary resistance between He and silver from 0.4 to 10 mk. J. Phys. Colloques 39, 265 (1978). 10.1051/jphyscol:19786117
29. Leggett A Takagi S Orientational dynamics of superfluid 3He: A “two-fluid” model. I. Spin dynamics with relaxation Ann. Phys. 1977 106 1 79 10.1016/0003-4916(77)90005-7
A. Leggett, S. Takagi, Orientational dynamics of superfluid He: A “two-fluid’’ model. I. Spin dynamics with relaxation. Ann. Phys. 106(1), 79 (1977). 10.1016/0003-4916(77)90005-7
30. Einzel D Wölfle P Transport and relaxation properties of superfluid 3He. I. Kinetic equation and Bogoliubov quasiparticle relaxation rate J. Low Temp. Phys. 1978 10.1007/bf00116904
D. Einzel, P. Wölfle, Transport and relaxation properties of superfluid He. I. Kinetic equation and Bogoliubov quasiparticle relaxation rate. J. Low Temp. Phys. (1978). 10.1007/bf00116904
31. Bunkov YM Dmitriev VV Mukharskiy YM Nyeki J Sergatskov DA Catastrophic relaxation in 3He-B at 0.4 Tc Europhys. Lett. 1989 8 7 645 10.1209/0295-5075/8/7/011
Y.M. Bunkov, V.V. Dmitriev, Y.M. Mukharskiy, J. Nyeki, D.A. Sergatskov, Catastrophic relaxation in He-B at 0.4 . Europhys. Lett. 8(7), 645 (1989). 10.1209/0295-5075/8/7/011
32. Surovtsev EV Suhl instability as a mechanism of the catastrophic relaxation in the superfluid 3He-B J. Low Temp. Phys. 2007 10.1007/s10909-007-9416-8
E.V. Surovtsev, Suhl instability as a mechanism of the catastrophic relaxation in the superfluid He-B. J. Low Temp. Phys. (2007). 10.1007/s10909-007-9416-8
33. Bunkov YM Golo VL A chaotic regime of internal precession in 3He-B J. Low Temp. Phys. 1993 10.1007/bf00681998
Y.M. Bunkov, V.L. Golo, A chaotic regime of internal precession in He-B. J. Low Temp. Phys. (1993). 10.1007/bf00681998
