Next Article in Journal
Detection of Tomato Leaf Miner Using Deep Neural Network
Previous Article in Journal
Application of Biosensors, Sensors, and Tags in Intelligent Packaging Used for Food Products—A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Non-Steady State NMR Effect and Application on Time-Varying Magnetic Field Measurement

1
School of Electronic Information and Communications, Huazhong University of Science and Technology, Wuhan 430074, China
2
College of Life Science and Technology, Huazhong University of Science and Technology, Wuhan 430074, China
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(24), 9960; https://doi.org/10.3390/s22249960
Submission received: 12 October 2022 / Revised: 11 December 2022 / Accepted: 15 December 2022 / Published: 17 December 2022
(This article belongs to the Section Physical Sensors)

Abstract

:
The measurement of a time-varying magnetic field is different from a constant magnetic field, due to its field intensity variation with time. Usually, the time-varying magnetic field measurement converts the solution of the magnetic induction intensity into the calculation of the induced electromotive force (EMF); then, the magnetic induction intensity is obtained by the time integration of the EMF, but the process is vulnerable to external interference. In this paper, a non-steady state nuclear magnetic resonance (NSS-NMR) scheme for the measurement of a time-varying magnetic field is proposed. In a time-varying magnetic field environment, an RF excitation signal with a certain frequency bandwidth is applied to excite the nuclear spin system. The NSS-NMR signal, which varies with time in the frequency range corresponding to the frequency bandwidth of the RF excitation, could finally be obtained after a series of processing of the probe output signal. During the NSS-NMR experiment, an orthogonal dual-coil probe is adopted to synchronously generate the RF excitation and induce the probe output signal. Moreover, a directional coupler that utilized in the experiment outputs a reference signal from the coupling port for the subsequent signal processing. The experimental results show that the weak NSS-NMR signal is indeed observed. The longitudinal time-varying magnetic field ranges from 0.576 T to 0.582 T, which is inverted by the Larmor precession relationship, have been successfully detected based on the so-called NSS-NMR effect.

1. Introduction

Magnetic field measurement provides reliable field-strength information for scientific research in a magnetic field environment. The common magnetic resonance methods, including the electron paramagnetic resonance (EPR) method, optical pump magnetic resonance (OPMR) method and nuclear magnetic resonance (NMR) method [1,2,3], could be utilized to measure a constant magnetic field [4,5]. The EPR method applies the microwave radiation resonant absorption effect of the paramagnetic system but it is susceptible to external interference, so the measurement accuracy mainly depends on the experimental environment [6,7]. The OPMR method is characterized by poor stability and is sensitive to ambient temperature fluctuation and mechanical vibration. In addition, the operation of the OPMR experimental system is rather complicated, resulting in a lower measurement accuracy [8]. The NMR method utilizes the inductive absorption effect for the incident electromagnetic wave in a nuclear spin system, for which the Larmor precession frequency is proportional to the magnetic induction intensity of the longitudinal magnetic field. Theoretically, since the accuracy of the gyromagnetic ratios of the common nucleus (usually 1H) and frequency measurement is extremely high, the accuracy of the magnetic field measurement by the NMR method will also be very high; for example, the relative accuracy is better than 10 ppm [9]. Therefore, the NMR measurement result is often used as the magnetic field measurement benchmark to calibrate other types of magnetometers [10,11].
Conventionally, the NMR method is applied to measure the magnetic induction intensity of a permanent magnet [12]. The radio frequency (RF) pulse sequence is utilized to excite the nuclear spin system in a conventional NMR experiment, the pulse width is inversely proportional to the measurement range of the time-varying magnetic field. It suggests that the wider the magnetic field measurement range is, the shorter the pulse width of the exciting pulse sequence will be. However, the shorter pulse width means the RF pulse power should be greater than that of the wider RF pulse for the effective excitation of the nuclear spin system. Wenjun Chen et al. had suggested that a continuous band-pass RF signal could be used to excite the nuclear spin system in an NMR experiment and simulated the so-called non-steady state nuclear magnetic resonance (NSS-NMR) effect in a pulsed high magnetic field [13,14,15], but they have not verified it in practice. In this paper, the NSS-NMR based time-varying magnetic field measurement scheme is proposed and verified in practice. After the data processing of the output signal of the orthogonal dual-coil probe, the NSS-NMR signal is successfully observed, and the longitudinal time-varying magnetic field is finally inverted by the Larmor precession relationship.
In an NSS-NMR experiment, the RF excitation signal is time continuous with a certain frequency bandwidth, so a pair of wideband dual-coil probes is designed and fabricated to excite the nuclear spin system and simultaneously induce the NSS-NMR signal. In this paper, the equivalent circuit of the probe and its extraction parameters, and two kinds of dual-port remote impedance matching methods are also described in detail. The experimental results show that the NSS-NMR effect indeed exists and could be used to measure the longitudinal time-varying magnetic field.

2. The NSS-NMR Effect and Its Numerical Simulation Result

2.1. The NSS-NMR Effect

A conventional NMR experiment is performed in a constant magnetic field while the nuclear spin system is simultaneously excited by an orthogonal-polarization RF magnetic field. The exciting condition of the NMR effect is that the angular frequency of the RF magnetic excitation field ω is equal to the Larmor precession frequency ω 0 [16], as depicted by the following Equation (1):
ω = ω 0 = γ B 0
where B 0 and γ are the longitudinal constant magnetic induction intensity and gyromagnetic ratio of the nucleus spin system, respectively. For example, the gyromagnetic ratio of 1H nucleus is 2.6752 × 10 8   r a d s 1 T 1 .
The Bloch equations systematically describe the NMR effect in the constant magnetic field from the perspective of classical mechanics [17]. Accordingly, if the magnetic induction intensity of the longitudinal magnetic field varies from B 01 to B 02 and the angular frequency ω of the exciting RF signal synchronously covers ω 1 (= γ B 01 ) to ω 2 (= γ B 02 ),the NMR effect may occur throughout the duration of the RF excitation signal. Since the nuclei are in the excitation and relaxation state alternately, the Bloch equations should be appropriately modified to describe the variation of the macroscopic magnetization and the proposed NSS-NMR effect. The NSS-NMR effect means the unstable and non-steady precession and nutation in the nuclear spin system. Therefore, the meaning of the NSS-NMR effect includes two aspects; since the macroscopic magnetization varies along with the longitudinal time-varying magnetic field and the orthogonal-polarization RF magnetic field, there is no stable equilibrium state or termination state for the precession and nutation.

2.2. Numerical Simulation Result of the NSS-NMR Effect

In order to acquire the variety of macroscopic magnetization in the NSS-NMR, the numerical integration method is performed to solve the modified Bloch equations under the longitudinal time-varying magnetic field and RF field co-excitation. Assuming that the amplitude of the longitudinal time-varying magnetic field B ( t ) is expressed as:
B ( t ) = k B 0 ( e m 1 t e m 2 t )
where B 0 is the peak value of the time-varying magnetic induction intensity, the parameter k is an adjustment coefficient, m 1 and m 2 are the tuning parameters corresponding to the peak value appearing time. In the simulation process, the polarization direction of B ( t ) is along the +z axis, the duration time is 0.01 s and the peak value is 0.6 T (at the peak value, the corresponding Larmor frequency of the 1H nucleus is 25.55 MHz). The experimental sample is distilled water and its volume is 1.57 × 10 9 m 3 . The polarization direction of a modulated broadband RF excitation signal B ( t ) is along the +x axis. The magnetic flux density passing through the intersecting surface of the receiving coil varies with time and, thus, generates an induced electromotive force (EMF) in the solenoid coil. Since the longitudinal magnetic field, the RF excitation and the macroscopic magnetization M are all time-varying parameters, and there is no magnetic field along the y axis, the classic Bloch equations could be modified as Equation (3):
{ d M x d t = γ M z ( t ) B z ( t ) M x ( t ) T 2 d M y d t = γ ( M z ( t ) B x ( t ) M x ( t ) B z ( t ) ) M y ( t ) T 2 d M z d t = γ M y ( t ) B x ( t ) M z ( t ) M z 0 ( t ) T 1
where T 1 and T 2 are the longitudinal relaxation time and transverse relaxation time, respectively, and the parameter M z 0 ( t ) represents the magnetization component of the thermal equilibrium state before the RF excitation signal is applied. Additionally, the total magnetic induction intensity B t o t a l ( t ) = B ( t ) + B 1 ( t ) = B x ( t ) a x + B y ( t ) a y + B z ( t ) a z is satisfied in the above equations.
In order to acquire the instantaneous components M x ( t ) , M y ( t ) and M z ( t ) of the macroscopic magnetization, the four-order Runge–Kutta (RK-4) method is adopted to solve the first-order linear, nonhomogeneous and variable coefficient differential equations. Before solving the equations, three conditions should be satisfied in advance. First of all, since the magnetic induction intensity of the longitudinal magnetic field is greater than that of the RF excitation field, the influence of the RF excitation on the longitudinal magnetic field could be ignored. Secondly, the spin relaxation time T 1 is independent from the magnetic induction intensity of the time-varying magnetic field B z ( t ) due to the narrow variation range of the time-varying magnetic field (~0.6T). Finally, the nuclear spin system should not be magnetized before the time-varying magnetic field is applied, which means that the macroscopic magnetization M z 0 ( t ) and M ( 0 , 0 , 0 ) are zeroes at the initial time t = 0 .
In the process of the numerical simulation, the pulsed time-varying magnetic field and sinusoidal time-varying magnetic field are employed for performing the NSS-NMR experiment to verify the NSS-NMR effect. The former is an aperiodic magnetic field with dramatic changes and a short duration (about 0.01 s). The latter is a periodic time-varying magnetic field, whose amplitude and period (T = 0.02 s) are consistent with that of the mixed magnets in the actual experiments. In the numerical simulation, the excitation signal is a band-limited white noise signal, which is a Gaussian white noise signal that processed by an equal ripple band-pass filter. Suppose the Larmor precession frequency corresponding to the peak value B 0 of the pulsed time-varying magnetic field is f 0 (= γ B 0 ), the parameters of the ripple band-pass filter utilized are as follows: the frequency range of the pass-band covers 0.2 f 0 to 1.0 f 0 with the pass-band fluctuation as 0.01 dB, the transition band ranges from 0.1 f 0 to 0.2 f 0 and 1.0 f 0 to 1.1 f 0 respectively, and the stop-band attenuation is 80 dB. Figure 1 shows the numerical simulation results of the pulsed time-varying magnetic field; the black curve in Figure 1a is the waveform of the longitudinal pulsed time-varying magnetic field. Figure 1b is the waveform of the broadband RF excitation signal whose amplitude and frequency band are about 0.006 T and 20 MHz, respectively. Figure 1c describes the NSS-NMR time-domain response of the macroscopic magnetization component M y ( t ) , whose waveform is completely different from that of the conventional NMR response. It can be concluded from Figure 1c that the parameter M y ( t ) is time continuing and lasts from 0.1 ms to 4.7 ms, which is consistent with the expected time duration of the NSS-NMR effect. Since the macroscopic magnetization increases as the magnetic induction intensity increases, the amplitude of M y is correspondingly very weak at the beginning and end stages. Figure 1d displays the time–frequency diagram of the M y obtained by short-time Fourier Transform (STFT), which is exactly in accordance with the time-varying magnetic field waveform shown in Figure 1a.
The numerical simulation results of the NSS-NMR effect in the sinusoidal time-varying magnetic field are shown in Figure 2. Figure 2a is the waveform of the sinusoidal time-varying magnetic field with a period of 0.02 s, amplitude range of 0.576 T to 0.582 T and the data length selected for the simulation as0.04 s. The red dotted line represents the magnetic field range corresponding to the bandwidth of the RF excitation signal. It should be emphasized that the broadband RF excitation applied in the simulation is identical to that in the pulsed time-varying magnetic field numerical simulation. Figure 2c displays the NSS-NMR time-domain response of the macroscopic magnetization component M y ( t ) of the sinusoidal time-varying magnetic field after digital down converters (DDC) due to its low frequency resolution. Figure 2c shows that the M y ( t ) is time-varying and its duration is 0.04 s. Figure 2d describes the time–frequency diagram versus the offset frequency from the corresponding median value B m e d i a n (=0.579 T) of the upper and lower ranges of the sinusoidal time-varying magnetic field. It can be seen from Figure 2d that a sinusoidal waveform is observed, which is consistent with the sinusoidal time-varying magnetic field waveform displayed in Figure 2a. The simulation of the above aperiodic and periodic time-varying magnetic field NMR experiment verifies the NSS-NMR effect and shows the feasibility of the magnetic field measurement by the NSS-NMR theory.

3. The NSS-NMR Experimental Scheme

3.1. The Time-Varying Magnetic Field and the Experimental System

The time-varying magnetic field environment is generated by a mixed magnet, which includes a permanent magnet and an energized multi-turn Hermholtz coil connecting to a tunable AC voltage source. The AC voltage can be adjusted from 0 V to 250 V by a knob installed on top of the energized coil. Therefore, the intensity of the magnetic field of the mixed magnet could be controlled by the AC voltage source. The magnetic induction intensity of the mixed magnet could be expressed in Equation (4) as:
B ( t ) = B 0 + B max cos ( 2 π f 0 t + φ 0 )
where B 0 corresponds to the magnetic induction intensity of the permanent magnet, and B max , f 0 and φ 0 represent the amplitude, frequency and initial phase of the magnetic induction intensity produced by the AC energized coil, respectively. The waveform of the time-varying magnetic field B ( t ) is shown in Figure 3, its DC component is about 0.579 T and the AC amplitude is about 0.003 T while the AC controlling voltage is set to 250 V. Corresponding to the magnetic induction intensity of the mixed magnet B ( t ) , the Larmor precession frequency should cover from 24.526 MHz to 24.782 MHz for the 1H nuclei; therefore, the frequency bandwidth of the NMR probe should exceed at least 260 kHz.
An NSS-NMR experimental system mainly includes a mixed magnet, an RF broadband signal source, an orthogonal dual-coil probe, a directional coupler, an attenuator and a set of data acquisition subsystems, as displayed in Figure 4. The RF excitation signal, which is a band-limited white noise signal with a frequency bandwidth of 300 kHz that is generated by the signal source, is sent to the input of the directional coupler after amplified by the power amplifier (PA). The two-way output signals from the directional coupler feed into the transmitting coil of the orthogonal dual-coil probe and the reference receiving port “1” of the data acquisition subsystem, respectively. The NSS-NMR signal, the leakage part of the transmitting signal and the noise outputting from the induction coil of the probe feed into port “0” of the data acquisition subsystem. The port “1” signal can serve as a reference signal for the following subsequent signal processing. The sampling rate is 100 MHz and the total acquisition time is set to 50 ms to cover multiple periods of the longitudinal time-varying magnetic field.

3.2. The Probe Design and Signal Processing Scheme

3.2.1. The Coil Structure Optimization

The signal-to-noise ratio (SNR) of an NMR signal is related to the quality of the NMR spectroscopy. In a conventional NMR experiment, the nuclear spin system is excited by a 90° pulse sequence [18,19], and the SNR of the free induction attenuation (FID) signal is illustrated in Equation (5).
S N R = P e a k   S i g n a l R M S   N o i s e = k 0 ω 0 2 ( B 1 / i ) V s N γ I ( I + 1 ) / 3 2 k T s V n = k 0 ω 0 2 ( B 1 / i ) M 0 V s 2 V n
In the above equation, the SNR is relevant to the distribution homogeneity k 0 of the RF excitation field, the resonance frequency ω 0 , the excitation efficiency B 1 / i , the macroscopic magnetization M 0 , the sample volume V s and the noise RMS voltage V n . Since the NMR signal is a macroscopic manifestation produced by the precession and nutation of all the spin nuclei with a non-zero magnetic moment in the sample space, the poor distribution homogeneity of the RF excitation field would eventually result in the gradual phase dispersion of the magnetization in the nutation plane [20,21]. The homogeneity parameter k 0 is related to many factors, such as the coil type, the length and diameter and the turn gap of the exciting/inducing coil, the gap between the coil and the wall of the magnet.
A common coil, such as a solenoid, Helmholtz, saddle and birdcage coil, could be used to fabricate an NMR probe [22,23,24,25]. The solenoid coil possesses a higher excitation efficiency and better distribution homogeneity of the magnetic field [26]; therefore, it was selected to design the transmitting coil and the receiving coil of the orthogonal dual-coil installed in the NSS-NMR experimental system. The orthogonal placement of the transmitting coil and receiving coil could obviously increase the isolation between the two coils. As revealed in Figure 5, the 1H nuclear sample that is encapsulated inside a short tube is positioned in the internal area of the vertical polyformaldehyde (POM) stick. The solenoid coil wrapping around the vertical POM stick is defined as the receiving coil. The six-turn transmitting solenoid coil is divided into two parts and wrapped on both sides of a thicker horizontal POM stick. The diameter of the enameled wire is 0.5mm, and the net length of the transmitting solenoid is 4.5 mm.
Since the turn gap and the length-to-diameter ratio are two significant parameters that affect the RF excitation field homogeneity [27,28,29,30], the magnetic field homogeneity parameter δ B 1 is defined to quantify the magnetic field distribution homogeneity in the center region of the transmitting coil. Assuming that the total number of discrete grids on each coordinate axis is N, then the average value of the magnetic induction intensity in the center region of the RF excitation field can be expressed in Equation (6) as:
B 1 m e a n = 1 N ( i N x B 1 x , i + j N y B 1 y , j + k N z B 1 z , k )
where B 1 x , i , B 1 y , j and B 1 z , k are the magnetic induction intensity at the i-th grid along the x-axis, y-axis and z-axis, respectively, N x , N y and N z are the numbers of the grids along each axis. Under the condition of i + j + k = N , the homogeneity of the RF excitation field δ B 1 is defined as Equation (7) as:
δ B 1 = 1 B 1 m e a n ( 1 N x i N x ( B 1 x , i 1 N x i N x B 1 x , i ) ) 2 + ( 1 N y j N y ( B 1 y , j 1 N y j N y B 1 y , j ) ) 2 + ( 1 N z k N z ( B 1 z , k 1 N z k N z B 1 z , k ) ) 2
The above definition is the same as the standard deviation of the RF excitation field distribution. It means thatthe closer δ B 1 is to zero, the better homogeneity of the RF excitation field will be. In the case that the diameter, length and position of the transmitting coil is fixed, a different turn gap and length-to-diameter ratio should be observed for the influencing degree of the parameter δ B 1 while the observed area is set as a sphere with a diameter of 1 mm. The front view of the observed area for the calculation of the RF homogeneity is interpreted as the red dotted region in Figure 6.
The homogeneity δ B 1 of the RF excitation field corresponding to the turn gap is displayed in Table 1. It can be concluded that δ B 1 reaches a minimum value when the equation d g a p = d c o i l is established.
The homogeneity δ B 1 of the RF excitation field corresponding to the length-to-diameter ratio is shown in Table 2. The results show that δ B 1 reaches a minimum value while the length-to-diameter ratio l / d is 0.5 rather than the maximum ratio 0.8. The reason is that the two-part structure of the transmitting solenoid is more similar to a multi-turn Helmholtz coil.

3.2.2. The Equivalent Circuit Parameters Extraction for the Orthogonal Dual-Coil Probe

The orthogonal dual-coil probe proposed in this paper consists of two mutually orthogonal solenoids; therefore, the equivalent lumped-parameter circuit could be depicted as Figure 7, where the mutual inductance M and mutual capacitance C 1 _ 2 between the two coils, the self-inductances L 1 and L 2 , the turn capacitances C 1 and C 2 and the ohm resistance R 1 and R 2 are included in the equivalent circuit.
The above lumped-circuit parameters can be calculated by many analytical methods, such as the method of moments (MoM) [31], the boundary element method [32], the simulated charge method [33] and the finite element method (FEM) [34]. The MoM could be used to calculate the distribution inductance between the wire loops, but it depletes the amount of memory and time to solve the multiple loops inductance. The boundary element method could effectively reduce the number of variables, but it is not convenient to solve the electromagnetic field problems. The simulated charge method replaces the free charge that is continuously distributed on the electrode surface with the simulated charge, but the setting of the optimal position of the charge is normally determined by experience, resulting in a low calculation accuracy. The FEM can analyze the distribution inductance among any complex wire loops, and its calculation process is relatively simple. Therefore, the FEM combined with the electric field energy and the magnetic field energy is adopted to extract the equivalent lumped-circuit inductances and capacitances of the orthogonal dual-coil probe by applying the correct boundary condition to the simulation region.
Assume that the coil number of a multi-loop system is N , the loop current is I i ( i = 1 , 2 , 3 N ), the total magnetic field energy of the multi-loop system W m , as shown in Equation (8) as:
W m = i = 1 N 1 2 ψ i I i = = τ 1 2 μ H 2 d τ
where H is the magnetic field intensity distributed in the whole region τ excited by the currents I i ( i = 1 , 2 , 3 N ), ψ i is the total magnetic flux crossing over the i-th coil, as described in Equation (9) as:
ψ i = M i 1 I 1 + M i 2 I 2 + + M i N I N = j = 1 N M i j I j
where M i j is the mutual inductance between the j-th ( j = 1 , 2 , 3 N ) coil and the i-th coil, and M i j = M j i , M i i = L i is established.
Therefore, the magnetic field energy W m can be expressed in Equation (10) as:
W m = τ 1 2 μ H 2 d τ = i = 1 N j = 1 N ( 1 2 I i M i j I j ) = ( 1 2 I 1 2 , 1 2 I 2 2 1 2 I N 2 , I 1 I 2 , I 1 I 3 , , I 1 I N , I 2 I 3 , I 2 I 4 , , I 2 I N , , I N 1 I N ) · ( L 1 , L 2 , , L N , M 12 , M 13 , , M 1 N , M 23 , M 24 , , M i N , , M N 1 , N ) T
In the above equation, the total number of inductances to be solved is N u m = N + C N 2 = N ( N + 1 ) 2 . Therefore, N ( N + 1 ) 2 groups of loop currents I ( i ) , i = 1 , 2 , , N ( N + 1 ) 2 , should be set to excite the multi-loop system for calculating the magnetic field distribution H ( i ) and the corresponding energy W M ( i ) in the whole space region. Where every exciting current group is I ( i ) = ( I 1 ( i ) , I 2 ( i ) , I 3 ( i ) I N ( i ) ) , i = 1 , 2 , , N ( N + 1 ) 2 , the magnetic field energy W m ( 1 ) , W m ( 2 ) , , W m ( i ) , , W m ( N ( N + 1 ) 2 ) , is as illustrated in the matrix Equation (11) below.
( W m ( 1 ) W m ( N ) W m ( N + 1 ) W m ( N ( N + 1 ) 2 ) ) = ( τ 1 2 μ ( H ( 1 ) ) 2 d τ τ 1 2 μ ( H ( N ) ) 2 d τ τ 1 2 μ ( H ( N + 1 ) ) 2 d τ τ 1 2 μ ( H ( N ( N + 1 ) 2 ) ) 2 d τ ) = ( 1 2 ( I 1 ( 1 ) ) 2 1 2 ( I 2 ( 1 ) ) 2 I n ( 1 ) I n + 1 ( 1 ) I N 1 ( 1 ) I N ( 1 ) 1 2 ( I 1 ( N ) ) 2 1 2 ( I 2 ( N ) ) 2 I n ( N ) I n + 1 ( N ) I N 1 ( N ) I N ( N ) 1 2 ( I 1 ( N + 1 ) ) 2 1 2 ( I 2 ( N + 1 ) ) 2 I n ( N + 1 ) I n + 1 ( N + 1 ) I N 1 ( N + 1 ) I N ( N + 1 ) 1 2 ( I 1 ( N ( N + 1 ) 2 ) ) 2 1 2 ( I 2 ( N ( N + 1 ) 2 ) ) 2 I n ( N ( N + 1 ) 2 ) I n + 1 ( N ( N + 1 ) 2 ) I N 1 ( N ( N + 1 ) 2 ) I N ( N ( N + 1 ) 2 ) ) ( L 1 L N M 1 , 2 M N 1 ,   N )
The above matrix equation could be abbreviated as W = I M , so all the inductances could be solved from the equation M = I 1 W . Here, every magnetic energy W m ( i ) is numerically integrated by using the magnetic field distribution H ( i ) calculated by the 3D FEM. It is very important to properly set up the value and direction of every loop current so as to guarantee the non-singularity of the matrix I .
The turn capacitances and the mutual capacitances of the dual-coil probe could also be solved by a similar method in the electric field environment. As a result, you only need to replace the current excitations with the multiple group charge excitations and calculate the electric field distribution by FEM and the energy density to calculate the electrostatic field energy under the different charge excitation groups. The extracted equivalent circuit parameters of the orthogonal dual-coil probe are displayed in Table 3. It can be concluded that the self-inductance, turn capacitance and ohm resistance values of the transmitting coil are all greater than those of the receiving coil, which accord with the facts that the turn number and diameter of the transmitting coil are all greater than those of the receiving coil. Moreover, the mutual electromagnetic coupling between the transmitting coil and the receiving coil is very weak due to the orthogonal placement of the coils reflected from the small mutual inductance and mutual capacitance, which indicates that the orthogonal dual-coil probe possesses a naturally high isolation between the input and output ports.

3.2.3. Dual-Port Matching for the Orthogonal Dual-Coil Probe

Since the RF excitation signal is broadband and time continuous in the NSS-NMR experiment, it is necessary to cascade the two impendence matching circuits on the two ports of the orthogonal dual-coil structure. The basic requirement of the matching circuits is to broadly match the impendence of the T/R coils to the impendence of the external source or load on the premise of ensuring enough isolation between the two ports. Due to the limitation of the available space (L × W × H = 54 × 11 × 60 mm3) for placing the matching circuits inside of the mixed magnet, a remote matching strategy is utilized practically, i.e., two long coaxial lines that link the orthogonal dual coils and matching circuits are placed outside of the magnet.
There are two kinds of remote matching methods, i.e., remote resistance matching and remote reactance multiple Π-stage matching for the dual-coil probe. The former possesses a wide frequency band, high isolation and simple circuit but a weak current excitation efficiency of the RF magnetic field. The latter possesses a strong current excitation efficiency but a narrow frequency band, poor isolation, complicated circuit and painful debugging process. If the continuous-wave average output power of the RF exciting source is enough or affluent, the remote resistance matching circuits are recommended for achieving a wider frequency band and higher isolation. Otherwise, the remote reactance matching circuits should be designed to achieve a high excitation efficiency under the conditions of a lower average power output of the excitation source.
Following the simulation of the aforementioned equivalent lumped-parameter circuit of the orthogonal dual-coil probe, only a 50 Ohm resistor and a capacitor are required in the remote resistance matching circuit, as depicted in Figure 8. The capacitors C t 1 and C t 2 are for tuning the imaginary part of the input impedance of the T/R coils, respectively, and the two 50 Ohm resistors are employed to match the very small real part of the input impedance of the T/R coils to the source impedance or load impedance.
The practical performance of the remote resistance matching method is revealed in Figure 9. Figure 9a,d show that the 20 dB and 10 dB impedance matching bandwidths are 5 MHz and 12.77 MHz, respectively. Figure 9b,c show that the 60 dB isolation bandwidth between the T/R ports exceeds 15 MHz, at a center frequency of 24.66 MHz of the NSS-NMR experiment on the condition of the 1H nuclear sample and the average value 0.579 T of the mixed magnet.
Figure 10 describes the remote reactance three-stage Π-type matching circuit, which is also designed and tuned practically for the actual probe prototype. Since the one-stage capacitance Π-type matching circuit only possesses a narrow bandwidth, the three-stage Π-type circuit is used in each remote port of the probe [35].
The prototype of the proposed wideband orthogonal dual-coil probe with the remote reactance matching circuits is devised and fabricated, as shown in Figure 11. The probe consists of a coil skeleton, nuclear sample, transmitting coil, receiving coil and matching circuits. Meanwhile, a thickness of up to 5 mm copper baffle (not shown in the figure) is designed to reduce mutual coupling between the outside T/R matching circuits.
Due to the non-magnetism of the resistors and capacitors, the latter matching circuits could also be installed in neighboring positions of the dual-coil structure when the inner space of the magnet is roomy enough. This near matching circuit may achieve a better performance and wider bandwidth than the remote one.
The S-parameter measurement results of the prototype probe are presented in Figure 12. The transmitting coil is connected to port 1 and the receiving coil is connected to port 2 of the vector network analyzer (VNA). A 10-dB impedance matching bandwidth is about 300 kHz, which coversfrom 24.50 MHz to 24.80 MHz. The isolation between the T/R ports is greater than 45.09 dB, which indicates that the performance of the very weak leakage from the RF exciting port to the receiving port will help to easily extract the NSS-NMR signal from the output compound signal in the subsequent signal processing process.

3.2.4. Signal Processing Scheme

Due to the finite isolation between the T/R coils, there is a small part of the RF excitation signal leaking from the input port to the output port. Although the leakage is very small, it is still much greater than the power of the NSS-NMR signal. Therefore, it is necessary to cancel the RF excitation leakage from the compound of the output signal to finally extract the NSS-NMR signal. As a result, a Wiener adaptive filtering interference cancellation algorithm is used in the data processing scheme. The reference signal comes from the coupling port of the directional coupler. The overall signal processing scheme includes band-pass filtering (BPF), delay correction, quadrature demodulation, low-pass filtering (LPF), sample extraction, adaptive interference cancellation (Wiener filtering) and time–frequency analysis, as depicted in Figure 13. The delay correction processing is for aligning the group-delay-time (GDT) of the probe output signal to that of the reference signal. The time-difference value corresponding to the peak of the cross-correlation function between the reference signal and output signal of the probe might be used to align the two signals. The quadrature demodulation, LPF and sample extraction are for outputting the digital baseband I/Q signals to conduct the subsequent adaptive filtering and time–frequency analysis.

4. Processing Results

Here, a segment of 30 ms time-length acquisition data with original sampling rate 100 MHz is utilized for verifying the NSS-NMR effect and the signal processing scheme. After quadrature demodulation, LPF and sample extraction, the sampling rate has been decreased to 2 MHz yet is covering the efficient frequency range of the NSS-NMR signal. Since the Wiener filter’s order is 800, the front 800 sampling points should be removed to avoid time–frequency analysis distortion. Figure 14 displays the waveform and corresponding power spectrum of the probe output signal before and after the interference cancellation. According to the time- and frequency-domain results, the interference cancellation degree of the Wiener filtering algorithm is about 46 dB. The residuary part in Figure 14d shows the suspected NSS-NMR signal.
After the interference cancellation, the processed data are divided into 228 sub-segments for carrying out STFT. Each sub-segment consists of 1024 sampling points with 75% overlapped points. The time–frequency diagram of the data of these segments versus the offset frequency from the NMR frequency corresponding to parameter B 0 is displayed in Figure 15a. A sinusoidal waveform could be clearly observed, which demonstrates the existence of the NSS-NMR effect, and the NSS-NMR signal has been successfully extracted. In addition, when the excitation RF power is low, the NSS-NMR signal is also weak in certain time segments. The inverted magnetic flux density of the mixed magnet by the Larmor precession relationship is shown in Figure 15b, which indicates that the range of time-varying magnetic field range is from 0.576 T to 0.582 T and the time-varying period is 20 ms. The experimental results are identical to the actual parameters of the mixed magnet. There is only a slight difference in the display clarity between the simulation and the actual experiment, which is mainly due to the spectral peak of the weak signals in certain segments submerged by the background noise.

5. Conclusions

A novel NSS-NMR effect is forecasted, and an NSS-NMR experimental scheme is proposed to verify the principle. In the NSS-NMR experiment, the nuclear spin system is excited by the continuous broadband RF signal; therefore, the excited NSS-NMR signal is composed of numerous signal components with different frequencies and amplitudes which vary with time. The NSS-NMR effect could be used to invert the magnetic induction intensity of the longitudinal time-varying magnetic field by the Larmor precession relationship. The NSS-NMR signal and the magnetic field inversion are successfully validated in a mixed sinusoidal time-varying magnet environment.
The next validation experiment would be carried out in a pulsed high magnet environment, where a more compact orthogonal dual-coil probe and more efficient matching circuits are needed to accommodate the new requirements. The NSS-NMR effect is even expected to provide a new observation method for studying the field-induced phase transition phenomenon and material structure with broad spectrum characteristics in a pulsed high magnetic field environment. In addition, the proposed NSS-NMR effect may be applied in material analysis (such as 27Al and 29Si nuclei in clay minerals) and liver disease detection.

Author Contributions

Conceptualization, X.Z., H.M., and J.M.; methodology, X.Z. and H.M.; validation, X.Z. and H.M.; formal analysis, X.Z. and H.M.; investigation, H.M., J.J. and H.Z.; resources, H.M., J.J. and H.Z.; Data curation, X.Z. and H.M.; writing—original draft preparation, X.Z. and H.M.; writing—review and editing, X.Z., H.M. and J.M.; visualization, H.M.; supervision, H.M.; project administration, H.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, Grant No. 10975056, 11475067; the Innovative Research Program of Huazhong University of Science and Technology, Grant No. 2015ZDTD017 and Grant No. 2017KFXKJC012; and the Experimental Apparatus Research Project of Wuhan Pulsed High Magnetic Field Center, Grant No. 2015KF17.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Zou, S.; Zhang, H.; Chen, X.; Fang, J. In-situ triaxial residual magnetic field measurement based on optically-detected electron paramagnetic resonance of spin-polarized potassium. Measurement 2022, 187, 110388. [Google Scholar] [CrossRef]
  2. Oida, T.; Ito, Y.; Kamada, K. Detecting rotating magnetic fields using optically pumped atomic magnetometers for measuring ultra-low-field magnetic resonance signals. J. Magn. Reson. 2012, 217, 6–9. [Google Scholar] [CrossRef]
  3. Kim, C.G.; Ryu, K.S.; Woo, B.C.; Kim, C.S. Low magnetic field measurement by NMR using polarized flowing water. IEEE. Trans. Magn. 1993, 29, 3198–3200. [Google Scholar] [CrossRef]
  4. Wang, X.; Zhu, M.; Xiao, K.; Guo, J.; Wang, L. Static weak magnetic field measurements based on low-field nuclear magnetic resonance. J. Magn. Reson. 2019, 307, 106580. [Google Scholar] [CrossRef]
  5. Frait, Z.; Fraitová, D. Measurement of magnetic field intensity by means of NMR. I. Accuracy of measurement. Czechoslov. J. Phys. B 1977, 27, 1292–1305. [Google Scholar] [CrossRef]
  6. Switala, L.E.; Black, B.E.; Mercovich, C.A.; Seshadri, A.; Hornak, J. An electron paramagnetic resonance mobile universal surface explorer. J. Magn. Reson. 2017, 285, 18–25. [Google Scholar] [CrossRef]
  7. Kittell, A.W.; Camenisch, T.G.; Ratke, J.J.; Sidabras, J.W.; Hyde, J.S. Detection of undistorted continuous wave (CW) electron paramagnetic resonance (EPR) spectra with non-adiabatic rapid sweep (NARS) of the magnetic field. J. Magn. Reson. 2001, 211, 228–233. [Google Scholar] [CrossRef] [Green Version]
  8. West, M.E.; Sesti, E.L.; Willmering, M.M.; Wheeler, D.D.; Ma, Z.L.; Hayes, S.E. Describing angular momentum conventions in circularly polarized optically pumped NMR in GaAs and CdTe. J. Magn. Reson. 2021, 327, 106980. [Google Scholar] [CrossRef]
  9. Zanche, N.D.; Barmet, C.; Nordmeyer-Massner, J.A.; Pruessmann, K.P. NMR probes for measuring magnetic fields and field dynamics in MR systems. Magn. Reson. Med. Sci. 2008, 60, 176–186. [Google Scholar] [CrossRef]
  10. Blanchard, J.W.; Wu, T.; Eills, J.; Hu, Y.; Budker, D. Zero-to ultralow-field nuclear magnetic resonance J-spectroscopy with commercial atomic magnetometers. J. Magn. Reson. 2020, 314, 106723. [Google Scholar] [CrossRef]
  11. Bloch, F.; Hansen, W.W.; Packard, M. Nuclear induction. Phys. Rev. 1946, 70, 460–475. [Google Scholar] [CrossRef] [Green Version]
  12. Stork, H.; Bontemps, P.; Rikken, G.L. NMR in pulsed high-field magnets and application to high-Tc superconductors. J. Magn. Reson. 2014, 234, 30–34. [Google Scholar] [CrossRef]
  13. Chen, W.; Ma, H.; Yu, D. Broadband Continuous Nuclear Magnetic Resonance Signal in a Pulsed Magnetic Field: Numerical Solutions of Bloch Equations under Radio Frequency Irradiation. Appl. Magn. Reson. 2016, 47, 41–52. [Google Scholar] [CrossRef]
  14. Chen, W.; Ma, H.; Yu, D.; Zhang, H. SVD-Based Technique for Interference Cancellation and Noise Reduction in NMR Measurement of Time-Dependent Magnetic Fields. Sensors 2016, 16, 323. [Google Scholar] [CrossRef] [Green Version]
  15. Chen, W.; Ma, H.; Yu, D.; Zeng, X. Numerical solution to the Bloch equations: Paramagnetic solutions under wideband continuous radio frequency irradiation in a pulsed magnetic field. Chin. Phys. C 2016, 40, 088201. [Google Scholar] [CrossRef] [Green Version]
  16. Rabi, I.I.; Millman, S.; Kusch, P.; Zacharias, J.R. The molecular beam resonance method for measuring nuclear magnetic moments. Phys. Rev. 1938, 53, 318. [Google Scholar] [CrossRef]
  17. Dougherty, W.M.; Bruland, K.J.; Chao, S.H.; Garbini, J.L.; Jensen, S.E.; Sidlesb, J.A. The Bloch Equations in High-Gradient Magnetic Resonance Force Microscopy: Theory and Experiment. J. Magn. Reson. 2000, 143, 106–119. [Google Scholar] [CrossRef]
  18. Hoult, D.I.; Richards, R.E. The signal-to-noise ratio of the nuclear magnetic resonance experiment. J. Magn. Reson. 1976, 24, 71–85. [Google Scholar] [CrossRef]
  19. Chen, C.N.; Sank, V.J.; Cohen, S.M.; Hoult, D.I. The field dependence of NMR imaging. I. Laboratory assessment of signal-to-noise ratio and power deposition. Magn. Reson. Med. 1986, 3, 722–729. [Google Scholar] [CrossRef]
  20. Podkorytov, I. Multipulse NMR part I. The simplest case: Quantum mechanical treatment of a single-spin system. Concepts Magn. Reson. 1996, 8, 17–32. [Google Scholar] [CrossRef]
  21. Hoult, D.I. The principle of reciprocity in signal strength calculations—A mathematical guide, Concepts. Magn. Reson. 2000, 12, 173–187. [Google Scholar]
  22. Roméo, F.; Hoult, D.I. Magnet field profiling: Analysis and correcting coil design. Magn. Reson. Med. 1984, 1, 44–65. [Google Scholar] [CrossRef]
  23. Strattan, R.D.; Young, F.J. Fields in square helmholtz coils. Appl. Sci. Res. Sec. B 1961, 9, 117–124. [Google Scholar] [CrossRef]
  24. Carlson, J.W. Currents and fields of thin conductors in rf saddle coils. Magn. Reson. Med. 1986, 3, 778–790. [Google Scholar] [CrossRef]
  25. Temnikov, A.N. An NMR Probe Coil. Instrum. Exp. Tech+. 2005, 48, 636–637. [Google Scholar] [CrossRef]
  26. Mispelter, J.; Lupu, M.; Briguet, A. NMR Probeheads for Biophysical and Biomedical Experiments; Imperial College Press: London, UK, 2006; pp. 382–394. [Google Scholar]
  27. Peck, T.L.; Magin, R.L.; Lauterbur, P.C. Design and analysis of microcoils for NMR microscopy. J. Magn. Reson. 1995, 108, 114–124. [Google Scholar] [CrossRef]
  28. Idziak, S.; Haeberlen, U. Design and construction of a high homogeneity rf coil for solid-state multiple-pulse NMR. J. Magn. Reson. 1982, 50, 281–288. [Google Scholar] [CrossRef]
  29. Privalov, A.F.; Dvinskikh, S.V.; Vieth, H.M. Coil design for large-volume high B1 homogeneity for solid-state NMR applications. J. Magn. Reson. 1996, 123, 157–160. [Google Scholar] [CrossRef]
  30. Mispelter, J.; Lupu, M. Homogeneous resonators for magnetic resonance: A review. Comptes Rendus. Chim. 2008, 11, 340–355. [Google Scholar] [CrossRef]
  31. Zheng, J.; Li, Z.F.; Qian, X.N. An efficient solver for the three-dimensional capacitance of the interconnects in high speed digital circuit by the multiresolution method of moments. IEEE Trans. Adv. Packag. 1999, 22, 9–15. [Google Scholar] [CrossRef]
  32. Shi, W.P.; Yu, F.Q. A divide-and-conquer algorithm for 3-D capacitance extraction. Comput. Des. Integr. Circuits Syst. 2004, 23, 1157–1163. [Google Scholar] [CrossRef]
  33. Harimurugan, D.; Punekar, G.S.; Krishna, K.N. Design of an HV capacitor using the inherent advantage of charge simulation method and experimentations. IET Sci. Meas. Technol. 2018, 12, 126–131. [Google Scholar]
  34. Liu, X.; Wang, Y.; Zhu, J. Calculation of Capacitance in High-Frequency Transformer Windings. IEEE Trans. Magn. 2016, 52, 2003204. [Google Scholar]
  35. Pozar, D.M. Microwave Engineering, 3rd ed.; John Wiley & Sons: New York, NY, USA, 2004; pp. 209–216. [Google Scholar]
Figure 1. The numerical simulation results: (a) The waveform of the pulsed time−varying magnetic field; (b) The waveform of the RF excitation signal applied; (c) The waveform of the macroscopic magnetization component M y ( t ) ; (d) The STFT result of the macroscopic magnetization component M y ( t ) .
Figure 1. The numerical simulation results: (a) The waveform of the pulsed time−varying magnetic field; (b) The waveform of the RF excitation signal applied; (c) The waveform of the macroscopic magnetization component M y ( t ) ; (d) The STFT result of the macroscopic magnetization component M y ( t ) .
Sensors 22 09960 g001
Figure 2. The numerical simulation results: (a) The waveform of the sinusoidal time−varying magnetic field; (b) The waveform of the RF excitation signal applied; (c) The waveform of the macroscopic magnetization component M y ( t ) ; (d) The STFT result of the macroscopic magnetization component M y ( t ) .
Figure 2. The numerical simulation results: (a) The waveform of the sinusoidal time−varying magnetic field; (b) The waveform of the RF excitation signal applied; (c) The waveform of the macroscopic magnetization component M y ( t ) ; (d) The STFT result of the macroscopic magnetization component M y ( t ) .
Sensors 22 09960 g002
Figure 3. The waveform of the mixed magnet.
Figure 3. The waveform of the mixed magnet.
Sensors 22 09960 g003
Figure 4. An NSS−NMR experimental system based on a mixed magnet, an orthogonal dual-coil probe and data acquisition subsystem.
Figure 4. An NSS−NMR experimental system based on a mixed magnet, an orthogonal dual-coil probe and data acquisition subsystem.
Sensors 22 09960 g004
Figure 5. The section view of the orthogonal dual-coil probe.
Figure 5. The section view of the orthogonal dual-coil probe.
Sensors 22 09960 g005
Figure 6. The front view of the observed area for the calculation of the RF excitation field homogeneity.
Figure 6. The front view of the observed area for the calculation of the RF excitation field homogeneity.
Sensors 22 09960 g006
Figure 7. The equivalent lumped-parameter circuit of the orthogonal dual-coil probe.
Figure 7. The equivalent lumped-parameter circuit of the orthogonal dual-coil probe.
Sensors 22 09960 g007
Figure 8. The remote resistance matching for the orthogonal dual-coil probe.
Figure 8. The remote resistance matching for the orthogonal dual-coil probe.
Sensors 22 09960 g008
Figure 9. The S−parameters of the remote resistance matching for the orthogonal dual-coil probe (a) S11; (b) S12; (c) S21; (d) S22.
Figure 9. The S−parameters of the remote resistance matching for the orthogonal dual-coil probe (a) S11; (b) S12; (c) S21; (d) S22.
Sensors 22 09960 g009
Figure 10. The remote reactance three-stage Π-type matching for the orthogonal dual-coil probe.
Figure 10. The remote reactance three-stage Π-type matching for the orthogonal dual-coil probe.
Sensors 22 09960 g010
Figure 11. The prototype of the orthogonal dual-coil probe with remote reactance matching circuits.
Figure 11. The prototype of the orthogonal dual-coil probe with remote reactance matching circuits.
Sensors 22 09960 g011
Figure 12. The measurement performances of the prototype probe: (a) S11; (b) S12; (c) S21; (d) S22.
Figure 12. The measurement performances of the prototype probe: (a) S11; (b) S12; (c) S21; (d) S22.
Sensors 22 09960 g012
Figure 13. The signal processing scheme for achieving the NMR signal of the NSS-NMR effect.
Figure 13. The signal processing scheme for achieving the NMR signal of the NSS-NMR effect.
Sensors 22 09960 g013
Figure 14. Results of: (a) The waveform of the original probe output signal. (b) The power spectrum of the original probe output signal. (c) The waveform of the probe output signal after interference cancellation by adaptive filtering. (d) The power spectrum of the probe output signal after interference cancellation.
Figure 14. Results of: (a) The waveform of the original probe output signal. (b) The power spectrum of the original probe output signal. (c) The waveform of the probe output signal after interference cancellation by adaptive filtering. (d) The power spectrum of the probe output signal after interference cancellation.
Sensors 22 09960 g014
Figure 15. The time−frequency analysis results: (a) The 3D STFT diagram. (b) The inverted magnetic flux density of the mixed magnet vs. time variation.
Figure 15. The time−frequency analysis results: (a) The 3D STFT diagram. (b) The inverted magnetic flux density of the mixed magnet vs. time variation.
Sensors 22 09960 g015
Table 1. Homogeneity of the RF excitation field corresponding to the turn gap.
Table 1. Homogeneity of the RF excitation field corresponding to the turn gap.
dgap0.5 rcoilrcoil1.5 rcoil
δ B 1 0.00950.00560.0102
Table 2. Homogeneity of the RF excitation field corresponding to the length-to-diameter ratio.
Table 2. Homogeneity of the RF excitation field corresponding to the length-to-diameter ratio.
l / d 0.20.40.50.60.8
δ B 1 0.0110.01540.01030.01550.0318
Table 3. The extracted equivalent circuit parameters of the orthogonal dual-coil probe.
Table 3. The extracted equivalent circuit parameters of the orthogonal dual-coil probe.
ParametersSelf-
Inductance
(Unit: μH)
Turn
Capacitance
(Unit: pF)
Ohm Resistance
of Coil
(Unit: mΩ)
Mutual
Inductance
(Unit: nH)
Mutual
Capacitance
(Unit: pF)
Transmitting Coil0.26251.14518.70.290.449
Receiving Coil0.13730.4918.6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zeng, X.; Ma, H.; Jin, J.; Zhang, H.; Ma, J. Non-Steady State NMR Effect and Application on Time-Varying Magnetic Field Measurement. Sensors 2022, 22, 9960. https://doi.org/10.3390/s22249960

AMA Style

Zeng X, Ma H, Jin J, Zhang H, Ma J. Non-Steady State NMR Effect and Application on Time-Varying Magnetic Field Measurement. Sensors. 2022; 22(24):9960. https://doi.org/10.3390/s22249960

Chicago/Turabian Style

Zeng, Xiaohu, Hong Ma, Jiang Jin, Hua Zhang, and Jingwen Ma. 2022. "Non-Steady State NMR Effect and Application on Time-Varying Magnetic Field Measurement" Sensors 22, no. 24: 9960. https://doi.org/10.3390/s22249960

APA Style

Zeng, X., Ma, H., Jin, J., Zhang, H., & Ma, J. (2022). Non-Steady State NMR Effect and Application on Time-Varying Magnetic Field Measurement. Sensors, 22(24), 9960. https://doi.org/10.3390/s22249960

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop