OPM magnetorelaxometry in the presence of a DC bias field

Spatial quantitative information about magnetic nanoparticle (MNP) distributions is a prerequisite for biomedical applications like magnetic hyperthermia and magnetic drug targeting. This information can be gathered by means of magnetorelaxometry (MRX) imaging, where the relaxation of previously aligned MNP’s magnetic moments is measured by sensitive magnetometers and an inverse problem is solved. To remove or minimize the magnetic shielding in which MRX imaging is carried out today, the knowledge of the influence of background magnetic fields on the MNP’s relaxation is a prerequisite. We show MRX measurements using an intensity-modulated optically pumped magnetometer (OPM) in background magnetic fields of up to 100μT\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$100~\upmu \mbox{T}$\end{document}. We show that the relaxation parameters alter or may be intentionally altered significantly by applying static fields parallel or antiparallel to the MNP’s alignment direction. Further, not only the relaxation process of the MNP’s magnetic moments could be measured with OPM, but also their alignment due to the MRX excitation field.


Introduction
Magnetic nanoparticles (MNP) offer a variety of applications in biomedicine, e.g. magnetic hyperthermia [1] and magnetic drug targeting [2]. All these applications require knowledge about the quantitative spatial distribution of MNP and magnetorelaxometry (MRX) is a promising technique for their quantitative imaging. In MRX, the magnetic moments of the MNP are aligned by an external magnetic field, forming a net magnetic moment. After switching-off the external field, the MNP's net magnetic moment decays, which can be measured by highly sensitive quantum magnetometers, e.g. superconducting quantum interference devices (SQUID) [3] and optically pumped magnetometers (OPM) [4][5][6], as well as fluxgates [7]. Whereas the relaxation amplitude contains quantitative MNP information, temporal properties of the relaxation curves can be exploited to estimate, e.g., the binding state of MNP [8] or their temperature [9]. MRX with spatially different excitation fields allows, after solving an appropriate inverse problem, to obtain quantitative spatial MNP information, e.g. with SQUID [10] or commercially available OPM [11]. Page 2 of 14 Currently, MRX imaging is performed within magnetic shieldings and thus at low background magnetic fields (<μT), although MNP quantification with MRX has been shown in noisy environments using SQUID gradiometers [12] and fluxgate gradiometers [7]. In Dolgovskiy et al. [13], OPM-MRX is performed within a weak magnetic shield. Novel OPM employing miniaturized vapor cells [14,15] are beneficial for magnetorelaxometry imaging measurements due to their small sensing volume and the possibility of flexible positioning. Moreover, these sensors have been shown to potentially reach sensitivities down to tens of fT/ √ Hz, even when operated within fields of geomagnetic strength (tens of μT) [16][17][18]. As magnetic shieldings are a crucial factor in terms of costs, availability and flexibility, it is desirable to avoid or minimize shielding in clinical applications. Unshielded magnetic particle spectroscopy (MPS), AC susceptometry (ACS) and magnetic particle imaging (MPI) with OPM were already shown [19,20]. Besides solving arising problems with environmental noise, a detailed understanding of how background magnetic fields influence the MNP relaxation processes is a prerequisite. It has been shown experimentally using SQUID, that Néel relaxation can be accelerated by applying external magnetic fields [21]. Transverse MRX of MNP wearing a viscoelastic corona, with an applied background magnetic field perpendicular to the excitation field has been modeled and investigated mathematically. It was shown in simulations that in transverse MRX, background magnetic fields alter the shape of the relaxation and decrease the relaxation time constant [22]. The dynamics of MNP's magnetic moments can be described by two concurrent processes: rotation of the whole particle is named Brownian motion, whereas the flipping of the single MNP's magnetic moments along their easy axis is described by so-called Néel relaxation [23]. As modeled mathematically by Fokker-Plank equations, it has been shown that the Brownian and Néel relaxation time constants decrease monotonically with increasing magnetic field strength, while for large fields (>mT) Néel relaxation is much more sensitive to magnetic fields [24]. Further, the dependence of the Néel relaxation time τ N on an applied field is well known [25,26]: where τ 0 is a constant and h = H H k , where H denotes the applied (background) field and H k the anisotropy field. Therefore, the Néel relaxation time decreases with an increase of the applied magnetic field.
In this paper, we present an OPM-MRX measurement setup for measuring within different defined background magnetic fields, ranging from 5 μT to 100 μT, and analyze the obtained data. We show, that the background magnetic field strongly influences MNP relaxation parameters. The relaxation parameters not only depend on the background magnetic field amplitude, but also on its direction. Finally, we show proof of principle measurements, where we detect not only the relaxation of MNP's magnetic moments, but also their alignment by the excitation field.

Setup overview
Our setup consists of a single channel OPM, an excitation coil for MNP alignment and a 3D Helmholtz coil system for the generation of defined background magnetic fields. The setup is positioned within a three-layer magnetic shielding barrel [27]. A MNP sample is placed on top of the OPM vapor cell, surrounded by the excitation coil (Fig. 1).

MNP excitation circuit
The MNP's magnetic moments are aligned by a magnetic field of about 1 mT (unless otherwise noted), generated by a 65-turn, 48 mm diameter coil (Fig. 1). After 200 ms (unless otherwise noted), this excitation field is shut-off within a few μs, which is achieved by a fast de-energizing of the coil. This is ensured by a low coil inductance and a fast-switching high-voltage MOSFET, combined with a high-voltage rated transient-voltage-suppression (TVS) diode (Fig. 2). The TVS diode clamps the back EMF to a (maximum) fixed level, protecting the MOSFET and decreasing the magnetic field linearly over time.

Intensity-modulated OPM
A schematic overview of the OPM is shown in Fig. 3. The OPM employs a microfabricated Cs vapor cell with nitrogen buffer gas (38 mbar) and a sensing volume of approx. 50 mm 3 . The cell consists of a 4 mm thick silicon body with a cylindrical cavity, closed by anodically bonded glass faces [28]. The vapor cell is operated at 70 • C by electrical thin-film heaters glued to the side faces of the Si body and driven by a 10 kHz ac current. Circularly polarized light from a distributed Bragg reflector (DBR) laser diode at the Cs D1 transition (895 nm) is intensity-modulated at the chopping frequency f c resonantly tuned to the Larmor frequency f L of the Cs atoms. It should be noted that f c is not feedback-controlled to avoid an increased settling time and transient signals after switching-off the MNP's excitation field. In this all-optical Bell-Bloom configuration, the magnetometer sensitivity is optimized when operated in static magnetic fields oriented perpendicular to the laser beam propagation direction [29,30]. The laser light power transmitted through the vapor cell is detected by a Si photodiode.

Data acquisition and preprocessing
The OPM photodiode signal is passed through a current voltage converter (I/U) and amplifier. Then, the signal is mixed with the OPM laser chopping frequency f c by the lock-in amplifier (LIA), and the resulting in-phase component (LIA-X) and quadrature component (LIA-Y) are directly digitized by the LIA after 4th order low-pass filtering with a -3 dB bandwidth of 10 kHz. A sample rate of 107.1 kHz is used. Additionally, the MNP excitation coil current is recorded by the same device, which serves as trigger for the data processing. In order to remove spurious frequency components emerging from the electrical ac heating, the OPM data is preprocessed in software by low-pass filtering with a cutoff-frequency of 1 kHz for liquid MNP and 100 Hz for MNP immobilized in gypsum. The digital filter is realized as bessel filter to preserve sharp edges in the data.

MNP and MRX model
Magnetic nanoparticles with a hydrodynamic diameter of 100 nm (plain BNF-Dextran) from Micromod (Micromod Partikeltechnologie GmbH, Rostock, Germany) were used in this experiment. The MNP are dextran based, water suspended and their iron-oxide core with a diameter of 45 nm consists of multiple 15-20 nm iron-oxide crystallites [31]. Two samples of 140 μl were prepared for the experiment. One sample is composed of 100 μl liquid (factory supplied) MNP and diluted in 40 μl distilled water, whereas the second sample consists of 100 μl MNP, embedded and thus immobilized in gypsum, resulting in a total sample volume of 140 μl. The undiluted iron concentration is about 15 mg/ml, resulting in an iron amount of 1.5 mg in each sample.
The MNP relaxation signal B(t) can be described in several ways. For an MNP system with uniform particle size and non-interacting particles, a single exponential decay of the associated magnetic flux density at the OPM's location is expected [32]: with the relaxation amplitude B R , the relaxation time constant τ R and the time t after the start of the relaxation. However, real MNP systems usually show a non-uniform particle size distribution. Assuming a magnetometer bandwidth of 500 Hz and a typical MNP anisotropy of 10 4 J/m 3 , particles with a core diameter in the interval [17.5 nm, 21 nm], while having a hydrodynamic diameter ≥100 nm, contribute to the MRX signal and therefore can be detected [2]. For MNP ensembles with an equally distributed diameter in the detectable size range, the relaxation can be described by [3,26]: with the amplitude B R and a time constant τ which depends on the excitation parameters and the MNP's anisotropy. The relaxation signal can also be described phenomenologically by a stretched exponential [32]: with the relaxation amplitude B R , the relaxation time constant τ R and the stretching parameter β. The offset B offset is introduced to compensate for magnetometer heading error, background magnetic fields and static MNP magnetization (e.g. due to the applied background magnetic field).
In immobilized MNP, Brownian motion is suppressed, while in water suspended samples both, Brownian motion and Néel motion may contribute to the MRX signal. Therefore, it reasonable to fit for both fractions [32]: with the contribution of bound (immobilized) MNP B b and unbound MNP B ub , where for each, one or a linear combination of the previously presented approaches may be well suited.
When fitting the measured data to such a model, the high number of parameters leads to an ill-conditioned problem. Even for simple models, the parameter variance and mutual interdependence is often very high, not only due to the contribution of environmental noise. To overcome this uncertainty, first a model is fitted to the experimental data and the absolute amplitude difference at two fixed time points is used as measure for the relaxation amplitude.
A robust parameter for the relaxation time is the integral relaxation time [33], often called correlation time [34]. It is denoted as the area under the amplitude-normalized relaxation curve. For the estimation of this parameter no curve-fitting was involved.

MNP relaxation data processing
Since the relaxation signal is of an exponential form, a high OPM bandwidth is required to capture early parts of the relaxation. Latter parts of the relaxation, however, do not require high bandwidths, but would benefit from a high magnetometer sensitivity. To satisfy both requirements, adaptive filtering or resampling can be implemented [32]. With a similar effect, in this work, the data is weighted exponentially during curve fitting. The fits of the data to the relaxation model, which was selected as the sum (5) of two stretched exponentials (4), were performed using the trust-region-reflective least squares algorithm [35] provided by Matlab®. For the extraction of the relaxation amplitude and the integral relaxation  Figure 4 shows the noise level of the OPM at different background magnetic fields, which was varied from 5 μT to 100 μT. The noise floor was around 500 fT/ √ Hz at 500 Hz for all configurations, limited by excess laser intensity noise, while the shot-noise limited OPM sensitivity level was 51 fT/ √ Hz. In previous detailed studies of a similar setup we were able to suppress excess intensity noise by subtraction of a reference photo diode signal by a factor of 3 [30]. Common visible spikes in the noise spectrum arise at the mains frequency and its harmonics, whereas the spikes occuring at a single background magnetic field only, e.g., the spike at 390 Hz at B 0 = 5 μT, arise from LIA mixing (leakage) with the ac current driving the OPM heater.

OPM characterization
The OPM frequency response was estimated using a spectrum analyzer (HP 35460A) and a random-noise-fed pancake coil, placed near the OPM (Fig. 5). The OPM bandwidth was estimated to be about 500 Hz. It is noted, that the bandwidth does not significantly change at different background magnetic fields ≤100 μT.
Next, the region of linear OPM response was estimated. At a fixed background magnetic field B 0 , and thus at a fixed Larmor frequency f L , the OPM laser chopping frequency was swept over a range of f L ± 3.5 kHz corresponding to B 0 ± 1 μT. The raw LIA-Y and LIA-magnitude (LIA-R) data are depicted in Fig. 6. Each of the measured LIA-Y transfer functions was used to convert the LIA-output to an equivalent magnetic field. While the region of linear response depends on the background magnetic field, only sensor data within a region of ±200 Hz or ±57 nT, respectively, was used for the MNP relaxation fits, while still keeping in mind the sensor bandwidth.

Estimation of relaxation parameters: dependence of liquid MNP's relaxation behavior on background magnetic fields
First, the sample of liquid 100 nm MNP was positioned in the system. Then, the excitation coil was pulsed (200 ms, 1 mT) and the relaxation of the MNP was measured. This was repeated at different background magnetic field strengths. As it can be seen in Fig. 7(a), the relaxation curve is changing in dependence of the B 0 field modulus. Regarding the excitation pulse direction with respect to the B 0 field, a parallel and an antiparallel configuration was investigated. It should be noted that the direction of the B 0 field was kept constant during the whole experiment, but instead the excitation pulse direction was flipped. When investigating the raw OPM signals, care must be taken of the nonlinear OPM response and the limited OPM bandwidth. Therefore, for this configuration the first 10 ms of data were clipped. The estimated relaxation parameters are shown in Fig. 8. It can be seen that for both, parallel and antiparallel excitation, the relaxation amplitude decreases with an increase of the background magnetic field. In contrast, the integral relaxation time decreases for antiparallel background magnetic fields, whereas it increases for parallel background

Estimation of relaxation parameters: dependence of immobilized MNP's relaxation behavior on background magnetic fields
The immobilized MNP sample was positioned in the system and the MNP relaxation was measured after pulsing the excitation coil. This was repeated at different background magnetic field strengths. Figure 9 shows data averaged from three measurement cycles and the estimated relaxation parameters are depicted in Fig. 10. It can be seen that the relaxation amplitude is decreased for parallel excitations, when increasing the background magnetic field. For antiparallel excitations, the amplitude varies only by a few percent. The integral relaxation time is increased and decreased for antiparallel and parallel excitations, respectively, when increasing the background field.

Proof of principle inverse MRX
While MRX at different background magnetic fields was investigated, also the inverse process, namely the alignment of the MNP's magnetic moments with the excitation magnetic field may be of interest for various applications. Especially for MRX imaging, it is desired to measure both, the alignment and relaxation of the MNP's magnetic moments, to reduce the data acquisition time or to improve the ill-posed inverse problem.
Here, we turn on and off an external 65 μT magnetic field periodically (at 0.1 Hz), while a constant background field of 5 μT is applied continuously to the system. The static field and the excitation field are applied antiparallel to each other. Thus, to track both processes (alignment and decay) simultaneously, it is desirable to have a magnetometer which is sensitive at both 5 μT and 60 μT, without being limited by a possible dead time or small sensor slew rate during the stepwise magnetic field change. Our OPM usually has a dynamic range of approx. 60 nT and so it cannot cover the whole range. To avoid this limitation, one possibility is to locally compensate the excitation field, e.g., using self-compensating solenoids [20]. However, the practical realization of field compensations might be tricky. A well known way of extending the OPM's dynamic range is to employ a feedback loop [36]. As stated above, we did not want to follow this approach as we wanted to avoid transient signals after switching-off the excitation field.
A novel approach, the "bi-chromatic intensity modulation", is presented: The center of the sensitive range of the sensor is set by the laser chopping (and thus lock-in amplifier) frequency f c , which corresponds to a Larmor frequency at a specific magnetic field amplitude. To obtain a sensitivity at both magnetic field levels, we apply a bichromatic signal to the intensity modulator, chopping the light at two frequencies corresponding to the Cs vapor Larmor frequency during alignment and decay of MNP's magnetic moments, respectively. We demodulate the photodiode signal with two time synchronized LIA at each individual reference frequency.
In Fig. 11 it can be seen that one of the LIA detects the relaxation, while the other one detects the alignment of the MNP's magnetic moments. Both signals are analyzed separately in the time intervals ranging from 2 ms to 18 ms. The relaxation amplitude and the integral relaxation time were estimated as B R = 3.7 nT and T R = 9.75 ms, respectively, while the alignment amplitude and the integral alignment time were estimated as B A = 3.1 nT and T A = 8.56 ms, respectively. Thus, the integration time constant and the amplitude decrease with increased magnetic field.

Discussion
We presented a novel experimental investigation of the influence of parallel and antiparallel background magnetic fields on the relaxation of MNP.
By increasing the antiparallel background magnetic field from 5 μT to 100 μT, the integral relaxation time decreased by 60% (from 15.2 ms to 6.1 ms) in our liquid MNP sample. In the parallel configuration, the integral relaxation time increased from 15.4 ms to 23.7 ms, which corresponds to an increase of 54%. This decrease of the integral relaxation time for antiparallel excitation and the increase of the integral relaxation time for parallel excitaton is in good agreement with the literature: As calculated by Deissler et al. [24], the magnetic energy density (for Néel and Brownian relaxation combined) of MNP for parallel fields decreases, while it increases for antiparallel fields, respectively. The magnetic energy density is proportional to the relaxation time [24].
Regarding the relaxation amplitude, three counteracting effects have to be considered. Firstly, the effectively applied excitation field is increased or decreased by the background magnetic field. It is well known that the relaxation amplitude increases linearly with an increase of the excitation field (for excitation fields up to several mT) [37]. Secondly, due to the extraction of the relaxation amplitude at two fixed time points, this parameter is coupled with the integral relaxation time. Thirdly and finally, the background magnetic field, if applied parallel to the excitation field, prevents some MNP from relaxing. In our liquid MNP sample, the absolute value of the relaxation amplitude decreased by 31% or 20% by increasing the background magnetic field for antiparallel or parallel excitation, respectively.
For the MNP response in our immobilized sample, the partially non-monotonical response of the extracted parameters to the background magnetic field may be explained by considering the non-homogeneous size distribution in our sample. As Coffey et al. [38] emphasize, the influence of the parameter h (cp. Eq. (1)) depends strongly on the energy barrier height of the Néel relaxation. Therefore, by increasing the background magnetic fields, the different size fractions in the MNP sample behave differently. According to the previously mentioned particle size effect study it is likely that hereby the MNP relaxation shape and as a consequence thereof the integral relaxation time changes.
The integral relaxation times of our liquid MNP sample are below 30 ms and the integral relaxation times of our immobilized MNP sample are in the region of several 100 ms. While the immobilized sample is governed by Néel relaxation only, in the liquid sample both Néel and Brownian relaxation occur. In view of the different time scales for the two types of samples, one might be inclined to neglect Néel relaxation in the liquid samples. However, this is only reasonable if Brownian and Néel relaxation can be considered to be uncoupled phenomena, but different groups have shown that in fact they are not [24,39]. As a consequence, modeling and experimental characterization practice has to take both effects into account in most cases. Improved understanding in mechanisms affecting Néel and Browian relaxation in different ways, like, e.g., orientation and strength of static magnetic fields, may help to disentangle intricate MNP dynamics. The importance of the interplay of the relaxation mechanisms is substantiated, e.g., by the study of Dolgovskiy et al. [13], where it has been shown that small variations in the MNP diameter distribution might have a big impact on the measurement results obtained by MRX.
The results from our proof of principle inverse MRX, where the integral time and amplitude decrease with increased magnetic field, are in accordance to the theory [24,40] and our other measurement results for antiparallel excitation.

Conclusion and outlook
In conclusion, we demonstrated MRX of MNP at background fields up to 100 μT using an intensity-modulated OPM. We have investigated experimentally, how the relaxation processes of liquid and immobilized MNP's magnetic moments alter or can be altered on purpose by applying static background magnetic fields in parallel and antiparallel to the excitation field. For biomedical applications it is particularly important to investigate and differentiate both liquid and immobilized MNP to obtain information about the environment of the particles [41]. In this respect, also intermediate states and environments will be of interest with the aim of a multicolor imaging [42]. The knowledge of the background magnetic field's effect on the relaxation is an important prerequisite for unshielded MRX and may lead to new spatial encoding schemes in MRX imaging, similar to [21]. In the future, background magnetic fields perpendicular to the excitation field will be investigated. Hand in hand, MNP model estimation studies with larger parameter variation in terms of background fields and excitation fields (amplitude, duration and direction) will be undertaken. The variation of excitation fields is especially important for unshielded MRX imaging or novel MRX imaging spatial encoding schemes. In unshielded MRX (imaging), the effect of the background magnetic field (e.g. Earth's magnetic field) has to be considered and modeled. The inverse problem in MRX imaging significantly benefits from inhomogeneous excitation fileds [43]. In a similar way, an additional spatial encoding scheme using (temporaly active) homogeneous or inhomogeneous background magnetic fields could potentially improve the imaging resolution and accuracy. When applying inhomogeneous background fields, the MNP sample's response will be a superposition of responses at various background magnetic fields, which have to be entangled in post-processing, e.g. using a multicolor approach [42].
In this paper we have shown that not only the relaxation, but also the alignment of MNP's magnetic moments can be measured with our intensity-modulated OPM. This additional information might decrease the critical data acquisition time, e.g., during a treatment, of MRX imaging by at least a factor of two. Furthermore, such an approach may open the way towards new MRX imaging modalities in the future, e.g. using stepwise excitation and relaxation schemes. For unshielded measurements, a feedback loop of the laser modulation frequency will be embedded to enhance the OPM bandwidth and measurement range [36,44]. The OPM sensitivity might be improved by adding a (magnetically shielded) reference cell [30] or by using multipass cells [45]. Obviously, a gradiometric OPM setup would facilitate background noise suppression in unshielded environments [46].
To sum up, we believe that given previous work achieving unshielded OPM measurements of biomagnetic signals (e.g. magnetoencephalography [46]) and the understanding of MNP's relaxation processes at different background magnetic fields, unshielded MRX (imaging) will be feasible. An unshielded setup will foster the impact and availability of MNP imaging, in particular with respect to monitoring of therapeutic applications such as magnetic hyperthermia and magnetic drug targeting.