Cover Page The handle http://hdl.handle.net/1887/50492 holds various files of this Leiden University dissertation Author: Wagenaar, Jelmer J.T. Title: Magnetic resonance force microscopy for condensed matter Issue Date: 2017-07-05
5 Probing the nuclear spin-lattice relaxation time at the nanoscale Nuclear spin-lattice relaxation times are measured on copper using magnetic resonance force microscopy performed at temperatures down to 42 mk. The low temperature is verified by comparison with the Korringa relation. Measuring spin-lattice relaxation times locally at very low temperatures opens up the possibility to measure the magnetic properties of inhomogeneous electron systems realized in oxide interfaces, topological insulators and other strongly correlated electron systems such as high-t c superconductors. This chapter has been published as J. J. T. Wagenaar et al., Probing the nuclear spin-lattice relaxation time at the nanoscale, Phys. Rev. Appl., vol. 6, p. 014007, Jul. 2016
68 magnetic resonance force microscopy for condensed matter 5.1 Introduction 1 Korringa 1950 2 Alloul 2014 3 Glover and Mansfield 2002 4 Degen et al. 2009 5 Alexson et al. 2012 Among the most informative probes of electron systems in solids is the nuclear spin-lattice relaxation rate 1/T 1. This value quantifies the damping of the nuclear spin precession due to the coupling to the electron spins. This in turn can be related to the momentum-averaged imaginary part of the dynamical spin susceptibility of the electron system, measured at the very low Larmor frequency of the nuclear spins. This is a most useful quantity revealing basic physics of the electron system. For instance, in Fermi liquids like copper, the Korringa relation 1, 1/(T 1 T) = constant, universally holds, while non-korringa behaviors play a pivotal role in establishing the unconventional nature of various strongly interacting electron systems 2. This mainstay of traditional NMR methods is hampered by the fact that it is a very weak signal that usually can be detected only in bulk samples 3. We demonstrate here an important step forward towards turning this into a nanoscopic probe, by delivering a proof of principle that 1/T 1 can be measured at subkelvin temperatures using magnetic-resonance force microscopy (MRFM). State-ofthe-art MRFM demonstrates an imaging resolution of < 10 nm, a 100-million-fold improvement in volume resolution over bulk NMR, by detecting the proton spins in a virus particle 4. In this experiment the statistical polarization of the protons is measured, since the Boltzmann polarization is too small to detect. Measurements of the recovery of the Boltzmann polarization after radio-frequent pulses relate directly to spinlattice relaxation times. The force-gradient detection of T 1 is measured before 5 in GaAs at temperatures down to 4.8 K. In this chapter we demonstrate T 1 measurements (see Fig. 5.3a) at temperatures which are lower by 2 orders of magnitude and a volume sensitivity of 3 orders of magnitude larger. Furthermore, we obtain the temperature dependence of the relaxation time T 1 satisfying quantitatively the Korringa relation of copper; see Fig. 5.4. This is a major step forward towards the ultimate goal of a device that can measure variations of T 1 at the nanoscale to study the properties of the inhomogeneous electron sys-
probing the nuclear spin-lattice relaxation time at the nanoscale 69 tems which are at the forefront of modern condensedmatter physics, such as the surface states of topological insulators 6, oxide interfaces 7 and other strongly correlated electron systems 8 including the high-t c superconductors 9. MRFM uses a magnetic particle attached to an ultrasoft cantilever (Fig. 5.1a) to detect the force gradients between sample spins and the tip. In most MRFM apparatuses a laser is used to read out the motion of the cantilever. The spins in the sample are manipulated using radio-frequent (rf) currents that are applied using a copper or gold nanowire 10. The laser and large rf currents in the copper or gold nanowires prevent experiments at millikelvin temperatures. To our knowledge, sample temperatures below 4.2 K have not been unequivocally demonstrated, only that the working temperature of the cryostat can remain at 300 mk 11. In our setup, a superconducting quantum-interference device (SQUID) 12 is used combined with a superconducting rf wire (Fig. 5.1b). This enables us to perform nuclear MRFM experiments down to temperatures of 42 mk, and we are able to verify that the sample indeed is at this temperature. 6 Chen et al. 2009 7 Mannhart and Schlom 2010; Richter et al. 2013; Kalisky et al. 2013; Scopigno et al. 2016 8 Dagotto 2005 9 Keimer et al. 2015 10 Poggio et al. 2007; Isaac et al. 2016 11 Degen et al. 2009 12 Usenko et al. 2011; Vinante et al. 2011a, 2012 5.2 Methods 5.2.1 Experimental setup The experimental setup closely resembles the one that was recently used to measure the dangling bonds on a silicon substrate 13. The NdFeB magnetic particle is 3.43 ± 0.02 µm in diameter, and has a saturation magnetization of 1.15 T. The cantilever s natural resonance frequency is f 0 = 3.0 khz with a mechanical quality factor of Q 0 = 2.8 10 4, when far away from the surface. An in-housedeveloped cryopositioning system is used to move the cantilever relative to the sample with a range of 1 mm in all dimensions. The absolute position of the cantilever is measured using three capacitive sensors. The sample holder is placed on a fine stage which is used to change the height of the sample in a range of 2.3 µm. For the fabrication of the detection chip consisting of 13 Den Haan et al. 2015
70 magnetic resonance force microscopy for condensed matter Figure 5.1: a) Scanning-electronmicroscope image (false color) of the magnetic particle after it is glued to the cantilever using electron-beam-induced deposition from a platinum-containing precursor gas. b) Optical-microscope image (false colors) of the detection chip. The blue horizontal wire at the top is the NbTiN rf wire. The gold-capped copper sample is connected to the mixing chamber for thermalization. The small black dot is where we perform the MRFM experiments. Next to the sample and below the rf wire is the pickup coil for detection of the cantilever s motion. c) Schematic drawing of a side view of the local NMR experiment. When the rf frequency meets the resonance γ condition f r f = 2π B0, copper nuclei at magnetic-field strength B0 are in resonance. (b) (a) y RF wire Cu x pick up 3.43 μm SQUID 30 μm (c) z y m 3.43 μm 625 khz RF-wire 55 mt height copper silicon a superconducting rf wire, a pickup coil, and the copper sample, we start with a NbTiN film on a silicon wafer with its natural oxide. The film has an average thickness of 378 nm, a resistivity of 92 µω cm and a critical temperature of 15.3 K. Using standard electron-beam-lithography techniques, with a negative resist, we pattern the rf wire structure and pickup coil. We use reactive ion etching in a SF6 /O2 plasma to etch the film. The copper is sputtered with a thickness of 300 nm and a roughness of 10 nm. In order to prevent oxidation of the copper, a gold capping layer of 15 nm is sputtered while the sample remains in high vacuum. The copper sample is thermalized using a patterned copper wire leading to a large copper area, which is connected to the sample holder via gold-wire bonds. The sample holder is connected to the mixing chamber of the dilution refrigerator via a welded silver wire. We position the cantilever above the copper, as indicated by a small black dot in Fig. 5.1b, close (7 ± 1 µm) to
probing the nuclear spin-lattice relaxation time at the nanoscale 71 the center of the rf wire and close (5 ± 1 µm) to the pickup coil for sufficient signal strength. While approaching the copper, we measure an enormous drop in the quality factor of the cantilever towards Q = 300 which strongly depends on the distance to the sample. This is caused by eddy currents in the copper induced by the magnetic fields of the magnetic particle attached to the moving cantilever. The drop in the quality factor limits the minimal height above the surface, reducing the signal strength. Natural copper exists in two isotopes. Both isotopes have a spin S = 3/2, 63 Cu has a natural abundance of 69% and a gyromagnetic ratio 63 γ/(2π) = 11.3 MHz/T and 65 Cu has a natural abundance of 31% and a gyromagnetic ratio 65 γ/(2π) = 12.1 MHz/T. In fcc copper, there is no quadrupolar coupling, since electric-field gradients are zero given the cubic symmetry. The magnetic particle on the cantilever couples to all nuclear spins in the sample via its field gradient. Every spin effectively shifts the resonance frequency of the cantilever by an amount f s. The alignment of the nuclear spins is perturbed when a radio-frequent pulse is applied to the sample. The spins that meet the resonance condition f r f = γ 2π B 0 are said to be in the resonant-slice, where f r f is the rf frequency and B 0 the local magnetic field. The total frequency shift due to all the spins within a resonant-slice can be calculated by the summation of the single-spin contributions f 0 = p s f s, where p is the net Boltzmann polarization of the nuclear spins. Those spins losing their net Boltzmann polarization (i.e., p 0) will cause a total frequency shift of f 0. This effect can be obtained at low radio-frequent magnetic fields in a saturation experiment with low enough currents to prevent heating of the sample. A rf current of 0.2 ma is applied in all measurements. Taking the approximation of a circular wire, we estimate the typical field strength in the rotating frame for a current of 0.2 ma at a distance r = 7 µm to be B 1 = µ 0 I 4πr 2.9 µt. The rf-field is mostly perpendicular to B 0 (see Fig. 5.1c). The saturation parameter s γ 2 B1 2T 1T 14 2 is used to calculate the ratio of total saturation which is given by s 1+s. For temperatures 14 Abragam 1961
72 magnetic resonance force microscopy for condensed matter 15 Pobell 2007 T < 250 mk, the expected values for the relaxation times are T 1 > 5 s. Given a T 2 = 0.15 ms 15, the saturation parameter becomes s > 30, indicating that the spin saturation is at least 95%, for those spins that satisfy the resonance condition. Since we work with field gradients, there will always be spins which are not fully saturated. For the sake of simplicity, we assume that we have a resonant-slice thickness d, within which the spins are fully saturated. The Rabi frequency is γ 2π B 1 30 Hz, so we assume that for pulse lengths of 1 s all levels are fully saturated. When all levels saturate, the magnetization after a pulse restores according to a single exponential, with a decay time equal to the spin-lattice relaxation time T 1. The frequency shift of the cantilever is proportional to the magnetization, and the time-dependent resonance frequency f (t) becomes: f (t) = f 0 + f 0 e (t t 0)/T 1 (5.1) 16 Poggio et al. 2007 We use the phase-locked loop (PLL) of a Zurich Instruments lock-in amplifier to measure the shifts in resonance frequency δ f = f (t) f 0 of the cantilever at a bandwidth of 40 Hz. The measurement scheme is as follows: First we measure the natural resonance frequency f 0 using a PLL. The PLL is turned off, and the rf current is turned on for 1 s. At t 0, the rf current is turned off. Shortly thereafter, the PLL is turned on, and we measure the frequency shift relative to f 0. The PLL is switched off during and shortly after the rf pulse in order to avoid cross talk. At frequencies larger than 1 MHz and currents higher than 1 ma, we measure about a few millikelvins increase in the temperature of the sample. The dissipation of the rf wire will be subject of further study, since for nanoscale imaging based on (proton) density, the magnetic rf fields need to be large (> 3 mt) enough for adiabatic rapid passages 16. In this experiment, we avoid the large currents in order to prevent heating. 5.2.2 Frequency noise We determine the shifts in resonance frequency δ f using a phase-ocked loop. The noise spectrum of the frequency shifts S δ f ( f ), which has units Hz 2 /Hz, consists of the
probing the nuclear spin-lattice relaxation time at the nanoscale 73 S δf (f) (Hz 2 /Hz) 10-0 Figure 5.2: The frequency noise Q=1767 at h=2.46 μm Q= 317 at h=1.31 μm S δ f ( f ) for two different cantilever 10-1 heights above the copper sample. The solid line is the expected frequency noise according to Eq. 5.2 10-2 with an added 1/f noise. The Q dependence of the 1/f noise suggests 10-3 that the additional noise is created by the oscillating currents in the 10-4 copper, which also cause the drop in the quality factor due to eddy currents. 10-5 10 10-6 -1 10 0 10 1 10 2 Frequency (Hz) thermal noise of the cantilever and detector noise S xdet ( f ). Albrecht et al. (1991); Kobayashi et al. (2009); Isaac et al. (2016) derive the frequency noise for frequency-modulation detection. We obtain the following expression, without the commonly used approximation f f 0 2Q, and with A the rms amplitude for the cantilever displacement: S δ f ( f ) = 4k BTQ f 2 2πA 2 ( ) k 0 f 0 1 + 2Q f 2 + f 0 S xdet ( f ) A 2 ( f 2 + f 2 0 4Q 2 ) (5.2) The detector noise S xdet ( f ) is determined by measuring the voltage noise S Vdet ( f ) of the SQUID. The voltage noise corresponds to a certain position noise of the cantilever. By measuring the thermal spectrum of the cantilever motion at different temperatures and by using the equipartition theorem x 2 = k BT k, we calibrated the transfer between the position 0 and the SQUID s output voltage. With this, we obtain S xdet ( f ) = 82 pm/ Hz. The thermodynamic temperature of the cantilever is determined to be 139 mk and the cantilever motion is kept constant at a rms amplitude of A = 110 nm. Representative measurements for two different heights are shown in Fig. 5.2. The solid line in Fig. 5.2 is Eq. (5.2) with no fitting parameters but with an added noise source equal to S add ( f ) = 1Hz 3. The numerical factor of 1 is chosen just for convenience, but the dependence on Q indicates that the Q 2 f 2 ad-
74 magnetic resonance force microscopy for condensed matter 17 Isaac et al. 2016; Kuehn et al. 2006 ditional noise is created by the oscillating currents in the copper, which also cause the drop in the quality factor due to eddy currents. In other MRFM apparatuses, the 1/f noise is shown to be caused by dielectric fluctuations 17, but in those setups the tip-sample distance is much smaller than in our case. We hope that nuclear MRFM on other systems than copper result in a smaller drop in the quality factor and therefore in a much better signal-to-noise ratio. Table 5.1: Overview of the setup parameters and NMR constants and parameters for copper (Lounasmaa, 1974; Huiku et al., 1986; Oja and Lounasmaa, 1997; Pobell, 2007). Parameter Variable Value Stiffness cantilever k 0 7.0 10 5 Nm 1 Bare resonance frequency f 0 3.0 khz Intrinsic quality factor Q 0 2.8 10 4 Radius magnet R 0 3.43/2 µm Magnetic dipole M 1.9 10 11 Am 2 Copper-layer thickness d Cu 300 nm Copper spin density ρ 85 nm 3 Spin copper nuclei S 3/2 Gyromagnetic ratio 63 Cu 63 γ/(2π) 11.3 MHz/T Gyromagnetic ratio 65 Cu 65 γ/(2π) 12.1 MHz/T Natural abundance 63 Cu 69 % Natural abundance 65 Cu 31 % Korringa constant 63 Cu κ 63 1.27 sk Korringa constant 65 Cu κ 65 1.09 sk Spin-spin relaxation time T 2 0.15 ms 5.2.3 Numerical calculation of frequency shifts 18 Isaac et al. 2016 To calculate the frequency shifts f 0 due to the spins in resonance [Eq. 5.1], the contributions of all spins meeting the resonance condition are integrated. Traditionally, the measured frequency shift of a mechanical resonator in MRFM is simulated using k = m z B z, the stiffness shift due to the coupling to a single spin 18. However, a recent theoretical analysis by De Voogd et al. (2017), supported by experiments performed by Den Haan et al. (2015), suggests that this is only an approximation. Following the new analysis, using f r f T 1 2 6 khz, we
probing the nuclear spin-lattice relaxation time at the nanoscale 75 obtain for the single-spin contribution f s f s = 1 2 f 0 k 0 Re((F 1 + F 2 + F 34 )) (5.3) F 1 = m B B 0 (5.4) F 2 = m B 0 B B 0 2 1 1 + iωt 1 (5.5) F 34 = m B B B 0 2 0 (5.6) The primes and double primes refer to the first and second derivative, respectively, with respect to the fundamental direction of motion of the cantilever x, which is in our experimental setup the x direction. B B 0 is the component along B 0. B B 0 is the perpendicular component. The term F 2 vanishes because ωt 1 1. m is the mean (Boltzmann) polarization in units J/T. For arbitrary spin S we find for m ( ) ( m = γ h S + 1 ) coth γ hb 0 S + 1 2 2 k B T 1 ( ) ] 2 coth γ hb0. (5.7) 2k B T In order to find the total direct frequency shift f 0 of a saturated resonant-slice with thickness d, we need to integrate all terms over the surface of the resonant-slice: f 0 = ρd f s ds (5.8) slice with ρ the spin density. To simplify the integration, we switch to spherical coordinates, with r the radius from the center of the magnet to the position in the sample, φ the angle running from 0 to 2π within the plane and θ the angle with the z axis, which is perpendicular to the copper surface. When the magnet is centered above a uniform sample that extends sufficiently far in the plane, we can assume that φ can be integrated from 0 to 2π. The distance to the resonant-slice r res can be expressed as a function of B 0, θ and φ. The upper boundary of θ equals π for an infinitely thick sample. The lower boundary θ min can be expressed analytically using a geometric relation
76 magnetic resonance force microscopy for condensed matter and the expression for the magnetic field of a magnetic dipole M: h = r res (θ min, φ) cos θ min (5.9) B 0 = µ 0M 4πr 3 3 cos 2 θ sin 2 φ + 1 (5.10) Combining the above equations we obtain an expression for θ min : ( ) 0 = 3 cos 2 θ min sin 2 φ + 1 cos 6 θ min ( 4πB0 h 3 µ 0 M ) 2 (5.11) Solving this equation gives us the final unknown integration boundary θ min. π 2π f 0 = ρd dθ f s (θ, φ, B 0 )dφ (5.12) θ min 0 The signal coming from a sample of thickness d Cu at height h is equal to an infinite sample at height h minus an infinite sample at distance h + d Cu. All parameters relevant for the experiment on copper are listed in Table 5.1. The only free parameters are the height above the surface h, the absolute value of which we are not able to determine accurately because of the eddy currents close at small heights, and the thickness of the resonant-slice. 5.3 Results and discussion 5.3.1 Direct frequency shifts In Fig. 5.3a we show a representative plot of a saturation recovery measurement performed at a temperature of 42 mk and a rf frequency of 542 khz. The cantilever is positioned at a height of approximately 1.5 µm above the surface of the copper. At f r f = 542 khz, the measured frequency shifts are larger than those observed at other frequencies. This is caused by the presence of a higher resonance mode of the cantilever which can be driven by the magnetic field of the rf wire. The resulting increased motion of the magnetic particle at the end of the cantilever gives an additional rf-field that will increase the number
probing the nuclear spin-lattice relaxation time at the nanoscale 77 (a) δf (mhz) (c) f 0 (mhz) 10 0 5 0-5 -10-15 -20-25 -30 5 0-5 -10-15 -20-25 -30-35 -40 f 0 single measurement 50 average exponential fit T -40 1 10 5 12.5 20 20 30 40 Time (s) h=1.31 μm and d=31 nm h=1.54 μm and d=38 nm h=1.99 μm and d=90 nm (b) f 0 (mhz) -10-20 -30 h=1.54 μm and 542 khz linear fit T -1 (K -1 ) 0 100 200 300 400 500 600 700 800 RF frequency f RF (khz) of spins that are saturated. The exact mechanism 19 for this can be found in Ch. 6. We measure the frequency shift up to t = 50 s and average 50 times. Subsequently, the data are fitted to Eq. (5.1) to extract f 0, T 1 and a residual frequency offset δ f 0. The direct frequency shift f 0 is plotted versus the temperature in Fig. 5.3b and fitted with a straight line according to the expected Curie s law f 0 m T 1, which implicitly assumes that the resonant-slice thickness is temperature independent. f 0 shows a saturation at the lowest temperatures. This indicates that the (electron) spins are difficult to cool, although the saturation could also be caused by the approach of the minimum temperature of the refrigerator. For every temperature, we collect at least three sets of 50 curves. The error bars are the standard deviations calculated from the separate fits. In Fig. 5.3c the direct frequency shift is plotted versus rf-frequency for three different cantilever heights. Every data point resembles a data set of 50 averaged single mea- Figure 5.3: a) The saturating rf pulse results in a direct frequency shift f 0 of the cantilever s resonance frequency. The shift relaxes with the nuclear spin s relaxation time T 1. The dashed orange line shows a single measurement of the frequency shift (with a moving average of 1 s); the solid blue line shows 50 averages. The red solid line is an exponential fit according to equation 5.1. From the fit we extract the direct shift f 0, the relaxation time T 1, and a residual frequency offset δ f 0. b) The direct frequency shift f 0 versus the temperature. The linear fit is expected according to Curie s law m T 1, which implicitly assumes that the resonant-slice thickness is temperature independent. c) f 0 as function of rf-frequency f r f for three different cantilever heights at T = 42 mk. The solid line is the calculated signal with the cantilever height h and resonant-slice width d as parameters. The two resonant-slices of the two copper isotopes largely overlap each other. The dashed vertical line at 542 khz indicates the location of a higher resonance mode, at which the cantilever can be driven using the rffrequency resulting in an amplification of the B 1 field and, therefore, the spin signal (see the main text).
78 magnetic resonance force microscopy for condensed matter surements. The solid line is a numerical calculation, with the only fitting parameters the resonant-slice width and height of the cantilever above the surface. The height of the cantilever is left as a fitting parameter, because the absolute height is not known with sufficient accuracy, due to the nonlinear behavior of the piezostack in the fine stage. Several mechanisms will contribute to the resonantslice thickness; for example, the width of the resonantslice is determined by the NMR linewidth, i.e., the saturation parameter, and may be further broadened by the nuclear spin diffusion and by the resonant displacement of the cantilever. These various contribution are discussed in more detail below. 5.3.2 The Korringa relation 20 Pobell 2007 For two different rf-frequencies, 542 and 625 khz, at two different heights, 1.54 and 1.31 µm, we measure the relaxation time T 1 versus the temperature. The results are shown in Fig. 5.4. We observe a linear dependence of T 1 on the temperature, following the Korringa law T 1 T = κ. T is the electron temperature, and κ is the Korringa constant. From the linear fits, we extract κ = 1.0 ± 0.1 sk and κ = 0.9 ± 0.2 sk, which is close to the expected value of the combined 63 Cu and 65 Cu Korringa constants 20, κ = 1.2 sk, which is measured in bulk copper. 5.3.3 The resonant-slice thickness The thickness of the resonant-slice is important for the resolution of possible future imaging or spectroscopic experiments or when similar measurements are to be performed on much thinner films. We also need the resonantslice width when modeling the signals we measured in Fig. 5.3c, because they are determined by the number of spins taking part. This results in a fit parameter of the resonant-slice thicknesses of 90, 38 and 31 nm (for heights of 1.31, 1.54, and 1.99 µm, respectively). In the case of our experiment on copper, we believe that the width of the resonant-slice is determined by the NMR linewidth, i.e., the saturation parameter, and may be further broad-
probing the nuclear spin-lattice relaxation time at the nanoscale 79 Relaxation time T 1 (s) 30 25 20 15 10 5 at 625 khz and h=1.31 um at 542 khz and h=1.54 um linear fit with κ=0.9±0.2 sk linear fit with κ=1.0±0.1 sk 50 100 150 200 250 Temperature T (mk) Figure 5.4: The inverse of the relaxation time of the copper nuclei measured by applying saturation pulses versus the temperature. We find the same linear dependence with temperature as is found in bulk copper samples, the so-called Korringa relation. In blue are data measured at a height of 1.54 µm at 542 khz, where a larger signal is achieved due to amplification of the rf-field by a higher cantilever mode. In green are data recorded at 1.31 µm at a frequency of 625 khz. We find κ = 1.0 ± 0.1 sk and κ = 0.9 ± 0.2 sk, respectively, with a linear fit with the error indicating the 95% confidence intervals. Every data point is an average of at least three sets of averaged data. The error bars give the standard deviation of the found relaxation times for the averaged data sets. ened by the nuclear spin diffusion and by the resonant displacement of the cantilever. Below, we discuss these various contributions. First, consider the steady-state solution for the z component of the magnetization m z of the spins in a conventional magnetic-resonance saturation experiment 21 : 21 Abragam 1961 m z = 1 + (2πT 2 f ) 2 1 + (2πT 2 f ) 2 + γ 2 B 2 1 T 1T 2 m, (5.13) where f is the detuning of the rf frequency from the Larmor frequency and s = γ 2 B1 2T 1T 2 is the saturation parameter. m is the magnetization in equilibrium. The value f for which m z = 1 2 m is a measure for the resonance thickness by d = 4π f γ r B 0 with r B 0 the field gradient in the radial direction. We find a thickness of 98, 58 and 45 nm for, respectively, the fitted thicknesses of 90, 38, and 31 nm. One additional possibility to broaden the slice in the radial direction is the effect of nuclear diffusion. For polycrystalline copper, we find a transition rate 22 of W = 30T 1 2. Using T 2 = 0.15 ms, we find W = 2.2 10 2 s 1. Our field gradients are too small to quench the spin-diffusion 23. New analysis after publication: After publication, we did a more detailed analysis on the first concept here discussed, as can be found in Ch. 4. With Eq. 4.26 we find a thickness of 30, 19 and 14 nm for the three found thicknesses. When we include an additional 15 nm for spin diffusion, and we neglect the possible influence of the motion of the cantilever, we come close to the found values of 90, 38, and 31 nm. Note that especially the 90 nm value is extracted with large uncertainty. 22 Abragam 1961; Bloembergen 1949 23 Isaac et al. 2016
80 magnetic resonance force microscopy for condensed matter The nearest-neighbor distance a is 0.256 nm. This gives us a spin-diffusion constant of D = Wa 2 = 15 nm 2 s 1. This yields a diffusion length l D = 2 Dt. With a pulse length of 1 s, we obtain a diffusion length of 7.6 nm. Since this diffusion happens at both the inside and the outside of the resonant-slice, the actual expected broadening is approximately 15 nm. From the quality factor at each height combined with the piezovoltage we used to drive the cantilever, we estimate the drive amplitudes to be 60, 37 and 100 nm, respectively, for the fitted thicknesses of 90, 38 and 31 nm. Note that the cantilever s motion is in the x direction, while in the calculation the resonant-slice thickness is in the radial direction. When the resonant-slice is broadened only in the x direction, the spin signal is 3 times lower than when it is broadened in the radial direction. All the above effects cannot be easily summed, which makes it difficult to precisely determine which of the above processes is the limiting factor in the thickness of the resonantslice. Altogether, we conclude that the NMR linewidth, and the motion of the cantilever together with spin diffusion qualitatively explain the experimentally found slice thicknesses. In most solids, T 2 will be larger than in copper. Therefore it is reasonable to assume that the resonantslice width in these future experiments will be smaller. For the calculations of this experiment, we take the resonant-slice thickness to be uniform. 5.4 Summary and outlook In conclusion, by performing nuclear magnetic-resonance force microscopy experiments down to 42 mk, we demonstrate that the nuclear spin-lattice relaxation time T 1 can be detected by measuring the spin polarization in a volume of about 3000 3000 30 nm 3 ; see Fig. 5.1c. A much smaller volume gives the same signal-to-noise ratio when a different sample than copper is measured and a smaller magnetic particle is used. First, copper is such a good conductor that our experiment is hampered by the eddy currents generated in the copper due to the motion of the
probing the nuclear spin-lattice relaxation time at the nanoscale 81 magnetic-force sensor. The subsequent drop in quality factor increases the frequency noise in the measurements by an order of magnitude. This prevents us from performing experiments at smaller tip-sample heights. When the magnet is closer to the sample, the magnetization of the nuclei is also larger. Second, by improving the operation temperature of the fridge, future experiments can be performed at 10 mk. By taking the ratio between the frequency noise within a bandwidth of 1 Hz at a Q > 1000 and the frequency shift of a single spin when the height between the magnet and sample is minimal, we find a volume sensitivity of (30 nm) 3. Further improvements are possible when using smaller magnets on the force sensor 24, which provide much higher field gradients. Previous measurements of the relaxation time with MRFM 25 have a volume sensitivity 3 orders of magnitude smaller 26. Our volume sensitivity is at least 4 orders of magnitude larger than what can theoretically be achieved in bulk NMR and even 10 orders larger than is usually reported in solid-state NMR 27. To substantiate our claim that using our technique new systems can be investigated, we discuss two examples. The first example is the measurements on nuclei that are in the vicinity of the surface state of a topological insulator. Bulk NMR on powdered samples show that the nuclei close to the surface state exhibit a Korringa-like relation 28. Our technique could be used to directly measure inhomogeneities in the surface state, even if it dives below the surface. As a second example, we turn to superconducting interfaces between two oxides, such as the measurement of the two-dimensional electron gas between the interface of lanthanum aluminate (LaAlO 3 ) and strontium titanate (SrTiO 3 ) that becomes superconducting below 300 mk 29. The pairing symmetry of this superconducting state is still an open question, and by measuring T 1 of the nuclei in the vicinity of the interface it could be possible to obtain information about the pairing symmetry, as was done in 1989 with bulk NMR to show that the high Tc superconductor YBa 2 Cu 3 O 7 is an unconventional superconductor 30. We believe that our technique to probe the electronic state 24 Overweg et al. 2015 25 Alexson et al. 2012 26 Garner et al. 2004 27 Glover and Mansfield 2002 28 Koumoulis et al. 2013 29 Reyren et al. 2007 30 Hammel et al. 1989
82 magnetic resonance force microscopy for condensed matter through T 1 measurements with nanometer resolution will lead to new physics in the field of condensed matter at the ground state of strongly correlated electrons. 5.5 Acknowledgments We thank F. Schenkel, J. P. Koning, G. Koning and D. J. van der Zalm for technical support. This work is supported by the Dutch Foundation for Fundamental Research on Matter (FOM), by the Netherlands Organization for Scientific Research (NWO) through a VICI fellowship to T. H. O., and through the Nanofront program. T. M. K. acknowledges financial support from the Ministry of Science and Education of Russia under Contract No. 14.B25.31.0007 and from the European Research Council Advanced Grant No. 339306 (METIQUM). J. J. T. W. and A. M. J. d. H. contributed equally.