Imaging non-collinear antiferromagnetic textures via single spin relaxometry

Antiferromagnetic materials are promising platforms for next-generation spintronics owing to their fast dynamics and high robustness against parasitic magnetic fields. However, nanoscale imaging of the magnetic order in such materials with zero net magnetization remains a major experimental challenge. Here we show that non-collinear antiferromagnetic spin textures can be imaged by probing the magnetic noise they locally produce via thermal populations of magnons. To this end, we perform nanoscale, all-optical relaxometry with a scanning quantum sensor based on a single nitrogen-vacancy (NV) defect in diamond. Magnetic noise is detected through an increase of the spin relaxation rate of the NV defect, which results in an overall reduction of its photoluminescence signal under continuous laser illumination. As a proof-of-concept, the efficiency of the method is demonstrated by imaging various spin textures in synthetic antiferromagnets, including domain walls, spin spirals and antiferromagnetic skyrmions. This imaging procedure could be extended to a large class of intrinsic antiferromagnets and opens up new opportunities for studying the physics of localized spin wave modes for magnonics.

Despite decades of continuous progress in magnetic microscopy, nanoscale imaging of the magnetic order in antiferromagnets still remains a notoriously difficult task 3 .A first strategy consists in using probe particles, like photons or electrons, whose properties are modified by their interaction with the sample magnetization.This approach can be pursued with several techniques including X-ray photoemission electron microscopy 8,9 and spinpolarized scanning tunneling microscopy 10,11 , the latter providing direct imaging of the Néel order with atomic scale resolution.Yet, these highly complex techniques require either synchrotron-based experimental facilities or conductive samples with nearly perfect surfaces in ultrahigh vacuum.A more versatile approach to magnetic microscopy relies on the detection of static magnetic fields generated outside the sample.In antiferromagnets, weak magnetic fields are produced by the rotation of tiny uncompensated moments whose orientations are linked to the Néel order.Mapping such fields requires the use of scanning magnetometers combining high field sensitivity with nanoscale spatial resolution.Besides magnetic force microscopy (MFM) pushed at its limits 7,12 , these performances are offered by a new generation of table-top magnetic microscopes employing a single Nitrogen-Vacancy (NV) defect in diamond as a non-invasive, atomic-size quantum sensor [13][14][15] .This technique was recently used to image the stray field distribution produced by spin textures in magnetoelectric antiferromagnets under ambient conditions [16][17][18] .However, it remains limited to the study of a special class of antiferromagnets producing strong enough static magnetic fields.
Following a recent proposal 4 , we demonstrate here that antiferromagnetic spin textures can be imaged by sensing the magnetic noise they locally produce, rather than their static magnetic stray fields.To this end, we employ a single NV defect in diamond for nanoscale noise sensing.Magnetic noise with a spectral component resonant with the electron spin transition of the NV defect is usually detected by recording variations of its longitudinal spin relaxation time T 1 19 .This method, commonly referred to as relaxometry, has been used for various purposes in the past years, including the study of Johnson noise in met- als 20,21 , current fluctuations in graphene devices 22 , and spin waves in ferromagnets [23][24][25] .In contrast to these previous works, we introduce a relaxometry-based imaging mode, which relies on the simple measurements of the photoluminescence (PL) signal of the NV defect under continuous laser illumination.
The imaging mechanism is illustrated by Fig. 1.The spin-dependent PL response of the NV defect is modeled by considering a closed two-level system corresponding to the ground states with spin projection m s = 0 (bright state) and m s = ±1 (dark state) 26 (Fig. 1a).Under continuous laser illumination, the PL intensity is linked to steady-state spin populations, which only result from the competition between the spin relaxation rate Γ 1 = 1 /2T1 and optically-induced spin polarisation in the m s = 0 sublevel with a rate Γ p that depends on the optical excitation power.Within this simple model, the evolution of the PL intensity with the spin relaxation time T 1 is shown in Fig. 1b for two different optical powers (see Supplementary Information).Under standard experimental conditions corresponding to Γ 1 Γ p , the NV defect is efficiently polarized in m s = 0 by optical pumping, leading to a high PL signal.However, if magnetic noise increases the spin relaxation rate such that Γ 1 ∼ Γ p , NV spin polarization is degraded and the overall PL intensity is reduced.Magnetic noise is therefore directly imprinted into the PL signal of the NV defect.In addition, depending on the optical excitation power, the PL drop occurs at different spin relaxation times (Fig. 1b), so that the laser power can be adjusted to optimize the imaging contrast.Here we show that this procedure can be employed for nanoscale, all-optical imaging of non-collinear antiferromagnetic spin textures by probing the magnetic noise they locally produce via thermally-activated spin waves 4 .
To demonstrate the efficiency of the method, we image the magnetic order in synthetic antiferromagnets (SAFs).These magnetic systems are made of thin ferromagnetic layers coupled antiferromagnetically through Ruderman-Kittel-Kasuya-Yoshida (RKKY) interlayer exchange coupling.The high tunability of magnetic parameters in SAFs offers many opportunities for antiferromagnetic spintronics 6 , from efficient spin-torque induced domain wall motion in racetrack memory devices 27 to the recent stabilization of antiferromagnetic skyrmions 7 .The SAF structure used in this work consists of a sputtered [Pt/Co/Ru] ×2 multilayer stack with broken inversion symmetry, in which the non-magnetic Ru-Pt spacer layer ensures RKKY-based antiferromagnetic coupling between the two identical Co layers (see Fig. 2a).In such samples, the Pt/Co interface provides a sizable interfacial Dzyaloshinskii-Moriya interaction (DMI) combined with an effective perpendicular magnetic anisotropy, that can be finely tuned by varying the Co thickness t Co 7 .Magnetic imaging is performed with a scanning-NV magnetometer operating under ambient conditions.A commercial diamond tip hosting a single NV defect at its apex (Qnami, Quantilever MX) is integrated into an atomic force microscope and scanned above the SAF surface (see Fig. 2a).At each point of the scan, a confocal optical microscope is used to record the PL signal of the NV defect under continuous optical illumination with a green laser.For the experiments described below, the flying distance of the NV spin sensor is d NV = 80 nm, as measured through an independent calibration procedure 28 .
We first focus on a SAF structure with a significant effective perpendicular magnetic anisotropy, and thus large antiferromagnetic domains.A typical PL map recorded with scanning-NV magnetometry is shown in Fig. 2b.It reveals a network of sharp PL quenching lines, which correspond to domain walls in the SAF.We note that very similar PL maps can be obtained above ultrathin ferromagnetic samples producing sufficiently strong magnetic fields (> 5 mT) to induce an efficient mixing of the NV defect's spin sublevels [29][30][31] .In a SAF, however, the static field generated above a domain wall is too weak to produce such a spin mixing, since this field mainly results from the slight vertical shift between the two compensated magnetic layers.This is confirmed by quantitative magnetic field imaging above a domain wall (Fig. 2c,d).For this measurement, a microwave excitation is combined with optical illumination to record the Zeemanshift of the electron spin resonance (ESR) frequency of the NV defect at each point of the scan 15 .The maximum stray field above the domain wall is around 500 µT, in agreement with the value expected from micromagnetic simulations (see Supplementary Information).Such a static field is not strong enough to induce a detectable PL quenching 29 .
In our experiments, PL quenching is rather induced by magnetic noise localized at SAF domain walls.This is verified through measurements of the longitudinal spin relaxation time T 1 of the NV defect using the experimental sequence shown in Fig. 3a.The NV defect is first initialized into the m s = 0 spin sublevel with a laser pulse.After relaxation in the dark during a time τ , a second laser pulse is used to read-out the final population in m s = 0 by recording the spin-dependent PL signal 32 .As illustrated in Fig. 3b-e, the measurement is performed for three different positions of the NV spin sensor.When the tip is retracted far away from the sample surface, we obtain T 1 = 860 ± 300 µs, a value smaller than the one commonly inferred for single NV defects in bulk diamond samples.We attribute this observation to magnetic noise generated by paramagnetic impurities lying on the diamond tip surface 32 .By engaging the tip above an antiferromagnetic domain, the spin relaxation time drops to T 1 = 120 ± 10 µs, indicating the presence of an additional source of magnetic noise with a spectral component at the ESR frequency of the NV defect.We checked that this noise comes from the magnetization of the sample by performing control experiments with the tip engaged on non-magnetic surfaces (see Supplementary Information).Finally, when the tip is placed above a SAF domain wall, T 1 is further shortened down to 22 ± 2 µs, revealing that magnetic noise is stronger at domain walls.As shown in Fig. 1b, a large PL quenching contrast is indeed expected at low optical excitation power when the T 1 time drops from 120 µs to 20 µs.In addition, the quenching contrast should almost vanish when the excitation power reaches the saturation of the NV defect optical transition.Mag- netic images recorded at different optical power confirm this effect [Fig.3f-h], illustrating that the imaging contrast can be optimized by adjusting the excitation power.
The observed variations in NV spin relaxation time are explained by magnetic noise generated by thermal populations of magnons.This is supported by micromagnetic simulations of the dispersion relation of spin waves and their power spectral density.In Fig. 4a we show the simulation geometry, which consists of an equilibrium magnetic configuration comprising Néel-type domain walls with an antiparallel alignment between the layers.In the uniform SAF domain and in the absence of applied magnetic fields, the dipolar coupling between the two ferromagnetic layers lifts the degeneracy between the acoustic and optic spin wave modes, as indicated by the solid lines in Fig. 4b.Propagation of spin waves is characterized by a quadratic dispersion in the wave vector k with a frequency gap determined by the perpendicular magnetic anisotropy and the RKKY interlayer exchange coupling.By including the measured value of the Gilbert damping (α ∼ 0.1), however, these two branches are smeared out as shown in the colour map of the simulated power spectral density in Fig. 4b.The bottom of the spin wave band, which is determined by the frequency of the uniform mode (k = 0), remains above the ESR frequency of the NV defect at zero field, f 0 = 2.87 GHz [Fig.4b].However, given the line broadening induced by damping, the NV frequency falls into the tail of the thermal power spectrum (see Supplementary Information), which results in a reduction in T 1 time when the tip is engaged above a uniform SAF domain.For a domain wall, on the other hand, we find as in the ferromagnetic case 33,34 a gapless mode that propagates parallel to the domain wall, with a nonreciprocity that arises from a combination of the Dzyaloshinskii-Moriya and dipolar interactions 35,36 (Fig. 4c).These modes have been recently observed experimentally with time-resolved scanning transmission X-ray microscopy (STXM) 37 .To illustrate how the existence of such gapless modes provide the spatial contrast in NV-based relaxometry, we compute a spatial map of the magnetic noise intensity at f 0 = 2.87 GHz in the presence of a domain wall in a realistic, disordered SAF system.To achieve this, we evoke the fluctuation-dissipation theorem; rather than computing the noise power spectrum directly from the magnetic response to fluctuating thermal fields, which represents a Langevin dynamics problem that is computationally intensive, we calculate instead the harmonic response of the magnetization dynamics to spatially random fields at the frequency f 0 (see Supplementary Information).Fig. 4d represents the resulting ensemble-averaged map of the magnetic noise intensity at the flying distance d NV of the NV defect.As observed in the experimental data, the strongest noise intensity is obtained at the domain wall, due to the thermally-excited gapless modes.We note that the NV defect is mainly sensitive to magnetic excitations with a wave vector on the order of the inverse of the flying distance.More precisely, this filter in k-space has a shape given by 23 k 2 e −2kdNV .For the experiments shown in Fig. 3, d NV ∼ 80 nm and the NV center is therefore mostly sensitive to thermal magnons with wave vectors lying between 5 and 30 µm −1 .From the dispersion relation in Fig. 4c, it appears that we are mainly detecting excitations which are slightly off the maximum power density, but nevertheless sufficient to obtain a significant imaging contrast.
So far, we have shown that all-optical NV-based relaxometry provides nanoscale imaging of domain walls in a SAF with perpendicular magnetic anisotropy.We now go a step further by demonstrating that this simple method can also be used to image more complex textures, such as antiferromagnetic spin spirals and skyrmions 7 .For this purpose, we exploit the high tunability of magnetic parameters in SAF systems.Keeping the very same stack structure, the Co thickness is first carefully adjusted in order to obtain a SAF with vanishing effective magnetic anisotropy [Fig.5a].In this case, the magnetic ground state corresponds to an antiferromagnetic spin spiral, whose period is given by λ = 4π A /D, with A the Heisenberg exchange parameter and D the effective DMI constant.A typical PL quenching image recorded with the scanning-NV sensor above a SAF with vanishing magnetic anisotropy is shown in Fig. 5b.A disordered spin spiral structure is clearly visible through alternat- ing bright and dark PL stripes.As for the domain wall case, micromagnetic simulations confirm that the magnetic noise intensity at f 0 is stronger in the part of the spin spiral in which the magnetization lies in the film plane, which provides the imaging contrast in NV-based relaxometry (see Supplementary Information).Analysis of line profiles of the PL quenching map indicates a spin spiral period λ ∼ 250 nm [Fig.5c].Considering an exchange parameter A = 20 pJ m −1 , this period corresponds to an effective DMI constant D = 1 mJ m −2 , in good agreement with earlier studies 7 .
Starting from a spin spiral state, antiferromagnetic skyrmions can then be stabilized by using a bias interaction 7 .This is realized by adding a uniformly magnetized bias layer coupled ferromagnetically to the bottom layer of the SAF with vanishing anisotropy.This coupling creates an effective out-of-plane bias field µ 0 H bias [Fig.5d], which plays the same role as an external magnetic field for the stabilization of magnetic skyrmions in ferromagnetic materials.A PL quenching image recorded above a biased-SAF is shown in Fig. 5e.It features isolated PL quenching spots which correspond to antiferromagnetic skyrmions.A statistical analysis of line profiles over a set of 52 isolated skyrmions indicate an average width at half maximum of about 76 nm [Fig.5f ], which is likely limited by the flying distance of the NV spin sensor.These results are in line with MFM measurements performed under vacuum on similar samples 7 .Despite the fact that the internal breathing mode frequencies of the skyrmion are larger than f 0 , the calculated noise map at this frequency indicates a stronger intensity above skyrmions (see Supplementary Information).We attribute this to the scattering of propagating spin waves by the skyrmions, which results in small displacements and deformations of their equilibrium profiles.
The present work introduces a new way to image non-collinear antiferromagnetic spin textures with nanoscale spatial resolution, which relies on the detection of magnetic noise locally produced by thermal populations of magnons.This is achieved by adding a relaxometry-based imaging mode to the scanning-NV magnetometry toolbox, which is based on measurements of variations in the NV defect PL signal induced by magnetic noise.
Beyond ordered antiferromagnetic structures like domain walls, spirals, and skyrmions, this imaging procedure could be extended to study magnetic order and disorder in other low-moment materials, such as domain structures in two-dimensional van der Waals systems with low Curie temperature in which spin fluctuations would become dominant under ambient conditions.The sensitivity to magnetic noise with nanoscale spatial resolution also opens up possibility of mapping out magnon-driven spin currents in different magnetic textures, which would provide valuable insights into magnon transport that are inaccessible through other experimental means.

FIG. 1 .
FIG. 1. Principle of the measurement.(a), Simplified two-level system used to model the spin-dependent photoluminescence (PL) of the NV defect.(b), Normalized PL intensity of the NV defect as a function of the spin relaxation time T1 for two different regimes of optical excitation power P .The parameter Psat denotes the saturation power of the NV defect optical transition.

FIG. 2 .
FIG. 2. All-optical imaging of domain walls in a synthetic antiferromagnet.(a), Simplified sketch of the scanning-NV magnetometer and the SAF structure.The two magnetic Co layers are antiferromagnetically coupled through the Ru/Pt spacer layer by RKKY exchange interaction.The exact composition of the stack is Ta(10)/Pt(8)/Co(tCo)/Ru(0.75)/Pt(0.6)/Co(tCo)/Ru(0.75)/Pt(3), with the thicknesses in nm.The effective perpendicular magnetic anisotropy of the SAF is controlled by varying the Co thickness tCo.(b), PL quenching image recorded with the scanning-NV magnetometer above a SAF sample with perpendicular magnetic anisotropy (tCo = 1.41 nm), showing an overview of the magnetic domain structure.(c), Closer view of a domain wall.(d), Static magnetic field distribution recorded above the area marked with the dashed rectangle in c.Only the component of the stray field along the NV quantization axis, BNV, is measured.

FIG. 3 .
FIG. 3. Measurement of the NV spin relaxation time.(a), Experimental sequence used to measure T1.The spindependent PL signal is integrated at the beginning of the reading pulse (300 ns window), and normalized using a reference PL value obtained at the end of the pulse 32 .(b-d), Sketches of the position of the NV sensor with respect to the magnetic structure.The tip is either (b) retracted, (c) engaged above a domain, or (d) engaged above a domain wall.(e),NV spin relaxation curves measured for the three different tip positions.The T1 time is obtained from a fit to an exponential decay (dashed lines).(f-h), PL quenching images recorded above the same domain wall with increasing excitation laser power P .In these experiments the saturation power of the NV defect optical transition is Psat = 540 µW.

FIG. 4 .
FIG. 4. Origin of the magnetic noise.(a), Equilibrium domain wall configuration obtained by micromagnetics simulation.xw and x d indicate lines along which spin wave propagation is characterized.(b), Dispersion relation of the propagating spin wave modes in the SAF domain, taken at x = x d in (a).The solid lines indicate a theoretical prediction, while the colour map represents the power spectral density (PSD) computed from the damped transient response in micromagnetics simulations.(c), Dispersion relation of the channelled domain wall spin waves computed at x = xw in (a).The white dashed lines in (b,c) indicate the position of ESR frequency of the NV defect at zero field, f0 = 2.87 GHz.(d), Simulated map of the magnetic noise intensity at the frequency f0, computed at a flying distance dNV = 80 nm above a domain wall configuration in the presence of disorder.The quantity δB ⊥,i 2 denotes the noise intensity perpendicular to the NV quantization axis, which is obtained from an average over 500 different realizations of the random driving field.

FIG. 5 .
FIG. 5. Imaging antiferromagnetic spin spirals and skyrmions via single spin relaxometry.(a), Sketch of a SAF sample hosting a spin spiral.The composition of the stack is identical to the one shown in Fig. 2(a), with a Co thickness (tCo = 1.47 nm) leading to a SAF with vanishing magnetic anisotropy.(b), PL quenching image of a spin spiral state.(c), Line profile along the white dashed line in (b).(d), Sketch of a SAF sample including an additional bias layer which enables to switch the spin spiral into antiferromagnetic skyrmions.The composition of the bias layer is Ta(10)/Pt(8)/[Co(0.6)/Pt(0.45)]3.(e), PL quenching image showing antiferromagnetic skyrmions in the biased-SAF sample.(f ), Line profile across the skyrmion marked with the circle in (e).The dashed line is data fitting with a Gaussian function.(g), Statistical distribution of the full width at half maximum w of the PL quenching spots, inferred from measurements over a set of 52 isolated skyrmions.