Introduction

Non-centrosymmetric crystal structures, such as the B20 phase of FeGe1,2, can host chiral spin textures like magnetic skyrmions3,4, which have been proposed as new structures for memory storage applications5 at room temperature6. Chiral spin structures in general are of significant interest in current magnetic research7,8. Dzyaloshinsky-Morya-interaction (DMI)9,10 causes a chiral symmetry break of the magnetic interaction and influences the dynamic properties of the magnetic system. For example, the spin wave dispersion becomes non-reciprocal11,12, as experimentally confirmed by Brillouin spectroscopy13 and an additional phase shift between neighboring spins of a spin wave affects its resonance intensity14,15. The space group P213 of the FeGe B20 phase has an inherent broken inversion symmetry, but does not impose chirality. The chirality, in this case, results from the specific atomic sites occupied by Fe and Ge inside the unit cell16. The magnetic properties of FeGe were previously studied using the Mössbauer effect17, vibrating sample magnetometry18 and magnetic susceptibility measurements19 making FeGe a magnetically well characterized material.

In the Heisenberg model of direct nearest neighbour interactions, spin waves (magnons) have a dispersion relation proportional to the square of the wave vector k . An antisymmetric contribution to spin-spin interaction results in an additional term in the dispersion relation linearly proportional to k11,12 and therefore a shift with regard to the gamma point. Then spin waves propagating in opposite directions at the same frequency have different wavelengths leading to complex standing waves with a moving phase front. This allows to detect modes, which would otherwise cancel and not be detectable in FMR.

We measured the magnetodynamic properties of a millimeter-sized disk shaped sample and a micron-sized wedge shaped sample of B20 FeGe using ferromagnetic resonance (FMR)20,21. Previous FMR measurements on this material22,23,24,25 were performed with millimeter sized single crystalline samples. Solving the Landau-Lifshitz-Gilbert equation (LLG)26,27 for an FMR like excitation28, we determined magnetic material parameters in the usual way29.

Sample Preparation

Stoichiometric FeGe was melted, using induction heating and, to guarantee homogeneity, re-melted twice and annealed for 130 h at 1000 K. Cylinders were formed and a high pressure high temperature synthesis inside a Kawai-type30 multianvil apparatus with Walker-type31 module was applied. This resulted in 95% polycrystalline B20 FeGe, confirmed by X-ray diffraction. A maximum of 5% of the sample material could consist of secondary phase Iron Germanium. Energy-dispersive X-ray spectroscopy (EDX) measurements also reveal local composition variations with accumulation of iron (Fe:Ge 55:45). Further investigations with Lorentz microscopy show the formation of helices and skyrmions (Fig. 1(a)) in accordance to32. Micron sized samples (Fig. 1(b)) with wedge shaped geometries were cut using a Focused Ion Beam (FIB – FEI Helios nanolab 600) and placed inside an R-Type microresonator33,34,35 using standard lift-off FIB (Omniprobe manipulator with Pt gas insertion system) technique. During the lift-off process a carbon coating with up to 100 nm thickness and up to 15% platinum contamination36 could not be avoided. Furthermore, the lift-off process used Gallium as cutting ions and resulted in a localized deposition of a maximum of 2.6% of Ga (as verified by EDX).

Figure 1
figure 1

(a) Lorentz microscopy image at 190 K and 150 mT of an FeGe slice cut from the original sample using standard lift-off FIB. The magnetic field points perpendicular to the sample. The black and white dots represent an ordered skyrmion lattice. (b) Scanning electron micrograph of the specimen inside an R-Type microresonator. The inset shows a schematic representation of the geometry and the directions of the magnetic field B during the experiment. (Dimension of the sample: a = 11.3 ± 0.1 μm, b = 10.9 ± 0.1 μm, c = 5.9 ± 0.1 μm, d = 5.0 ± 0.1 μm, e = 0.9 ± 0.1 μm, f = 1.6 ± 0.1 μm).

Experimental

FMR spectra of a bulk polycrystalline, nearly disc shaped piece of FeGe with a diameter of 3.78 mm and a thickness of 0.78 mm (2(b) inset) was acquired in a range of 800 mT to 0 mT at a frequency of 9.517 GHz ± 0.006 GHz. The field was applied at angles of -8° to 172° (corresponding to the directions “up” and “down” in Fig. 2(a)) in steps of 0.5° from out-of-plane to in-plane and to the opposite out-of-plane orientation. The measurement of the uniform FMR mode can be seen in Fig. 2(a) shown as an amplitude contour plot. The temperature is 276 K ± 1 K, which is below the Curie Temperature of TC = 280 K19, where the sample is ferromagnetic32. The angular precision of our experimental setup is better than 0.05° and the precision of the magnetic field is better than 0.5 mT with a relative precision of 0.005 mT.

Figure 2
figure 2

(a) Angular dependent out-of-plane (differentiated) FMR spectra shown as an amplitude contour plot at 276 K ± 1 K and fMicrowave = 9.517 GHz ± 0.006 GHz. The yellow line marks the angular dependent resonance field position. The dotted white lines mark the position of the hard direction at 0° and 180°. They have been extended to the middle of the figure for better comparison of the 30 mT field difference. (b) shows the spectra of the same FMR measurement at the positions “up” and “down” marked in (a) . The resonance spectrum for the field applied along the “down” direction consists of two resonance lines. This is due to edge resonances inside the sample24,42 superimposing with the uniform FMR mode23. Additionally a schematic representation of the sample can be seen with the most important field positions marked.

Resonance lines in the FMR spectra are identified by a successive local maximum and minimum amplitude. We observe a distribution of resonances, which is in agreement to previous FMR investigations22 of single crystalline FeGe. Each crystallite in the sample is contributing to this resonance distribution. They are all influenced by the applied external field and the demagnetization field in the sample, due to its general shape. However, their resonance fields vary with respect to the applied magnetic field due to the different symmetry axis of the cubic anisotropy in each crystallite. We simulated the resonance distribution using the known magnetocrystalline anisotropy of FeGe22 and a random orientation of crystallites and compared it to the measurement. This can be found in the Supplementary Sec. S1. Figure 2(a) shows the differentiated angular dependent FMR spectra as a grey scale contour plot. The out-of-plane orientations are depicted in detail in Fig. 2(b). The measured resonance line exhibits a unidirectional anisotropy, indicated by a difference in the positon of maximum microwave absorption comparing opposite magnetic field directions. A similar anisotropy is observed in systems with exchange-bias37. Hence we performed additional magnetometry measurements, to exclude the presence of exchange bias in our system (Fig. 3). No such anisotropic behaviour is observed in static magnetometry using vibrating sample magnetometry (VSM). We therefore conclude, that this anisotropy must be dynamically induced under resonant excitation. Note that it cannot be equated with the linear contribution to the spin wave dispersion, as this changes directionality in accordance with the magnetic field direction. Due to the skin depth of approximately 10−3 mm22 one must, to fully reproduce the FMR lineshape, solve the non-uniform LLG24 taking the shape of the sample into account. However, we show exemplary in the Supplementary Sec. S1 that a Dysonian lineshape38,39 and the known magnetocrystalline anisotropy of FeGe18 are able to reproduce the measured FMR lineshape satisfyingly, which is sufficient for our needs. The position of resonance was obtained by subtracting the background and locating the zero crossing of the resonance line. We analyzed the angular dependent spectra using Eq. 1 as a model for the free energy density F. To account for the observed unidirectional symmetry in the angular dependent resonance field position, an additional unidiretional field contribution needs to be introduced. In this model an additional anisotropy field BU = KUD/M is used. This unidirectional contribution is merely a descriptive model to account for the observed phenomenon. It cannot be seen as an additional magnetocristalline anisotropy but rather as an emergent symmetry contribution which arises under dynamic excitation. In the Supplementary Sec. S2 the shape of such a unidirectional free energy density is shown. Additionally, a demagnetization and Zeeman term are considered.

$$F=-\,{K}_{UD}\cdot \,\cos (\theta )+\frac{{\mu }_{0}}{2}\cdot \overrightarrow{M}\cdot N\cdot \overrightarrow{M}-\overrightarrow{M}\cdot \overrightarrow{B}$$
(1)
Figure 3
figure 3

(a) Temperature dependent measurement of the magnetization at 310 mT external field. A part of the original sample was used for the measurement. (b) Hysteresis loop measured by vibrating sample magnetometry at 276 K. The sample is the same as in (a). The magnetisation M is plotted against the magnetic field B. The hysteresis shows no asymmetry or exchange bias.

The demagnetisation tensor N (\({N}_{zz}=0.676,{N}_{xx,yy}=0.162\)) was deduced, using the demagnetisation tensor of a cylinder as described in40. The known g-factor of FeGe (g = 2.07)23 was used. θ is the out-of-plane angle of the magnetisation M, and B is the external magnetic field. Additionally, the magnetisation M is considered as a fit parameter. The obtained parameters are \({K}_{UD}=960\,{{\rm{J}}/{\rm{m}}}^{3}\pm 10\,{{\rm{J}}/{\rm{m}}}^{3}\), and \(M=82580\,{{\rm{Am}}}^{-1}\pm 200\,{{\rm{Am}}}^{-1}\). The magnetization matches the magnetization measured by VSM at 281 K, 5 K above the temperature measured by a sensor below the sample. This offset is likely due to microwave heating.

Figure 4 shows the angular dependent FMR spectra (293 K ± 2 K, fMicrowave = 9.134 GHz ± 0.006 GHz) of the wedge shaped FeGe sample (Fig. 1(b)) measured inside a microresonator as a grey scale contour plot. Multiple resonances are visible in the spectra, which exhibit anisotropic behavior. The anisotropy is directed such that the resonance field increases when the static field is applied parallel to the long (dipolar-easy) axis of the sample. This suggests that these modes are spinwaves with energies below that of the gamma point (FMR mode), which may be induced by strong dipolar coupling41 or DMI. Around ±90°, the number of superimposing resonances and the complex mode intesity distribution15 make it difficult to separate individual lines. We assume that these resonances arise due to geometrical confinement of the modes in our specimen (Fig. 1). Consequently, the inclined surface of our wedge results in different geometrical boundary conditions at the same time. Bidirectional measurements along the ±81° direction as shown in Fig. 5, however, reveal a clear unidirectional shift of the resonances under field reversal. Figure 5(a) shows the reproducibility of resonances for field sweep up and field sweep down, whereas Fig. 5(b) illustrates that under field reversal, the resonance position of the spinwaves has shifted. Hence, we find a unidirectional anisotropy. Figure 5(c) shows the noise floor of our spectrometer in a field region where no resonances are observed.

Figure 4
figure 4

A grey scale contour plot of the (differentiated) FMR signal amplitude of the micron sized wedge shaped sample as a function of applied magnetic field for different orientations of the magnetic field between −93° and 99° (compare Fig. 1(b)) at 9.134 GHz ± 0.006 GHz. The scale bar is depicted on the left. The dotted white lines indicate spin wave modes.

Figure 5
figure 5

Analysis of the bidirectional (differentiated) FMR measurements (~293 K) at 81° and 3.587 GHz ± 0.006 GHz, of the specimen shown in Fig. 1(b). (a,ii) and (b,ii) comparison between two different (differentiated) FMR measurements, the former with the same magnetic field direction, the latter at opposite field directions. (a,i) and (b,i) show a schematic representation of the sample and the respective magnetic field directions of the compared measurements. (a,iii) and (b,iii) depict the highlighted areas in (a,ii) and (b,ii) in detail. (a,iv) and (b,iv) are the plotted differences between the compared measurements. (a,v) noise floor of the measurement in a magnetic field region without resonances.

Conclusion

From angular dependent ferromagnetic resonance spectroscopy we find an unexpected dynamic unidirectional anisotropy (Fig. 2) in the magnetic excitation of FeGe just below the Curie temperature. This anisotropy is of a dynamic character, since it is not detectable in static hysteresis measurements. The magnitude of the unidirectional anisotropy of the bulk resonance line is \({K}_{UD}=960\,{{\rm{J}}/{\rm{m}}}^{3}\pm 10\,{{\rm{J}}/{\rm{m}}}^{3}\). Spin waves, detected at 293 K ± 2 K for sample sizes with micrometer dimensions, also exhibit unidirectional anisotropy (Fig. 5).

Methods

A conventional Bruker X-band FMR spectrometer was used for FMR measurements on the millimetre sized FeGe sample (see Fig. 1) inside a cylindrical TE011 cavity. FMR measurements on the micron sized FeGe sample (Figs. 1(b), 4 and 5) were performed inside an R-Type microresonator33,34,35. The resonator was connected to a Varian E102 microwave bridge. The modulated microwave reflection was recovered using a SRS SR830DSP lock-in amplifier.