Improved diffusion-weighting efficiency of pulsed gradient stimulated echo MR measurements with background gradient cross-term suppression

Accurate diffusion measurements with pulsed gradient NMR are hampered by cross-terms of the diffusion-weighting and background gradients. For experiments based on a stimulated echo pulse sequence, that is preferred for samples with a T2 short compared to the diffusion time, a diffusion-weighting scheme has been presented that avoids these cross-terms in each of the en- and decoding periods separately. However, this approach suffers from a reduced diffusion-weighting efficiency because the two gradients applied in each of the periods have effectively opposite polarities leading to a partial cancellation. An extension of this scheme is presented that involves an additional gradient pulse in each period and delivers an improved diffusion-weighting efficiency without sacrificing the cross-term compensation. Analytical expressions for the gradient pulse lengths and amplitudes are given for arbitrary timing parameters. MR measurements with artificial (switched) background gradients were performed to test the cross-term compensation capability of the proposed extension. The results show that considerably higher q and b values can be achieved with the extension without changing the timing parameters. The MR measurements yielded identical diffusion coefficients without, with the same, and with different background gradients in the en- and decoding periods demonstrating the cross-term compensation of the presented approach.     

Cross-term-compensated pulsed-gradient stimulated echo MR with asymmetric gradient pulse lengths

The magic asymmetric gradient stimulated echo (MAGSTE) sequence developed to compensate background-gradient cross-terms in the preparation and readout interval independently, assumes identical lengths for the two gradient pulses applied in each interval. However, this approach is rather inefficient if some extra delay time is present in one half of an interval, e.g. as required for special RF excitations or spatial encoding prior to the stimulated echo in MR imaging. Therefore, a generalized version of the sequence is presented that considers different gradient pulse lengths within an interval. It is shown theoretically that (i) for any pulse lengths a "magic" amplitude ratio exists which ensures the desired cross-term compensation in each interval and that (ii) prolonging one of the gradients can deliver a considerably higher diffusion weighting efficiency. These results are confirmed in MR imaging experiments on phantoms and in vivo in the human brain at 3T using an echo-planar trajectory. In the examples shown, typically 10 times higher b values can be achieved or an echo time reduction with a 40% signal gain in brain white matter. Thus, in case of asymmetric timing requirements, the generalized MAGSTE sequence with different gradient pulse lengths may help to overcome signal-to-noise limitations in diffusion weighted MR.     

Crossing fibers, diffractions and nonhomogeneous magnetic field: correction of artifacts by bipolar gradient pulses

In recent years, diffusion tensor imaging (DTI) and its variants have been used to describe fiber orientations and q-space diffusion MR was proposed as a means to obtain structural information on a micron scale. Therefore, there is an increasing need for complex phantoms with predictable microcharacteristics to challenge different indices extracted from the different diffusion MR techniques used. The present study examines the effect of diffusion pulse sequence on the signal decay and diffraction patterns observed in q-space diffusion MR performed on micron-scale phantoms of different geometries and homogeneities. We evaluated the effect of the pulse gradient stimulated-echo, the longitudinal eddy current delay (LED) and the bipolar LED (BPLED) pulse sequences. Interestingly, in the less homogeneous samples, the expected diffraction patterns were observed only when diffusion was measured with the BPLED sequence. We demonstrated the correction ability of bipolar diffusion gradients and showed that more accurate physical parameters are obtained when such a diffusion gradient scheme is used. These results suggest that bipolar gradient pulses may result in more accurate data if incorporated into conventional diffusion-weighted imaging and DTI.     

Spin echo formation in the presence of stochastic dynamics

Spin echo formation in magnetic field gradients in the presence of fast stochastic motion is studied for hyperpolarized 3He gas at different diffusivities. The fast translational motion leads to frequency shifts already during echo formation, which can be described analytically for a linear gradient. Despite complete signal loss at the position of the spin echo itself, considerable intensity can be preserved at an earlier time (sqrt[2]tau rather than 2tau, where tau is the pulse delay). Hence, the phenomenon is designated as a pseudo spin echo.     

Steady-state free precession with hyperpolarized 3He: experiments and theory

The magnetization response of hyperpolarized 3He gas to a steady-state free precession (SSFP) sequence was simulated using matrix product operators. The simulations included the effects of flip angle (alpha), sequence timings, resonant frequency, gas diffusion coefficient, imaging gradients, T1 and T2. Experiments performed at 1.5 T, on gas phantoms and with healthy human subjects, confirm the predicted theory, and indicate increased SNR with SSFP through use of higher flip angles when compared to optimized spoiled gradient echo (SPGR). Simulations and experiments show some compromise to the SNR and some point spread function broadening at high alpha due to the incomplete refocusing of transverse magnetization, caused by diffusion dephasing from the readout gradient. Mixing of gas polarization levels by diffusion between slices is also identified as a source of signal loss in SSFP at higher alpha through incomplete refocusing. Nevertheless, in the sample experiments, a SSFP sequence with an optimized flip angle of alpha=20 degrees, and 128 sequential phase encoding views, showed a higher SNR when compared to SPGR (alpha=7.2 degrees) with the same bandwidth. Some of the gas sample experiments demonstrated a transient signal response that deviates from theory in the initial phase. This was identified as being caused by radiation damping interactions between the large initial transverse magnetization and the high quality factor (Q=250) birdcage resonator. In 3He NMR experiments, performed without imaging gradients, diffusion dephasing can be mitigated, and the effective T2 is relatively long (1 s). Under these circumstances the SSFP sequence behaves like a CPMG sequence with sinalpha/2 weighting of SNR. Experiments and simulations were also performed to characterize the off-resonance behaviour of the SSFP HP 3He signal. Characteristic banding artifacts due to off-resonance harmonic beating were observed in some of the in vivo SSFP images, for instance in axial slices close to the diaphragm where B0 inhomogeneity is highest. Despite these artifacts, a higher SNR was observed with SSFP in vivo when compared to the SPGR sequence. The trends predicted by theory of increasing SSFP SNR with increasing flip angle were observed in the range alpha=10-20 degrees without compromise to image quality through blurring caused by excessive k-space filtering.     

Generalized MAGSTE with bipolar diffusion-weighting gradient pulses

The generalized magic asymmetric gradient stimulated echo (generalized MAGSTE) sequence compensates background gradient cross-terms and can be adjusted to asymmetric timing boundary conditions which for instance are present in echo-planar MR imaging. However, its efficiency is not optimal because one of the two diffusion-weighting gradients applied in each interval usually must have a reduced amplitude to ensure the desired cross-term compensation. In this work, a modification of generalized MAGSTE is investigated where this gradient pulse is replaced by two gradient pulses with full amplitude but opposite polarities. It is shown that with these bipolar gradients (i) the sequence retains the cross-term compensation capability for an appropriate choice of the gradient pulse durations and (ii) the diffusion-weighting efficiency is improved, i.e. higher k and b values can be achieved without prolonging the echo time. These results are confirmed in MR imaging experiments on phantoms and in vivo in the human brain at 3 T using spin-echo and echo-planar MR imaging. In the examples shown, the b value could be increased between about 30% and 200% when using the bipolar gradient pulses. Thus, bipolar gradients may help to improve the applicability of the generalized MAGSTE sequence.     

Radiofrequency magnetic field gradient echoes have reduced sensitivity to susceptibility gradients

The amplitudes of gradient-echoes produced using static field gradients are sensitive to diffusion of tissue water during the echo evolution time. Gradient-echoes have been used to produce MR images in which image intensity is proportional to the self-diffusion coefficient of water. However, such measurements are subject to error due to the presence of background magnetic field gradients caused by variations in local magnetic susceptibility. These local gradients add to the applied gradients. The use of radiofrequency (RF) gradients to produce gradient-echoes may avoid this problem. The RF magnetic field is orthogonal to the offset field produced by local magnetic susceptibility gradients. Thus, the effect of the local gradients on RF gradient-echo amplitude is small if the RF field is strong enough to minimize resonance offset effects. The effects of susceptibility gradients can be further reduced by storing magnetization longitudinally during the echo evolution period. A water phantom was used to evaluate the effects of background gradients on the amplitudes of RF gradient-echoes. A surface coil was used to produce an RF gradient of between 1.3 and 1.6 gauss/cm. Gradient-echoes were detected with and without a 0.16 gauss/cm static magnetic field gradient applied along the same direction as the RF gradient. The background static field gradient had no significant effect on the decay of RF gradient-echo amplitude as a function of echo evolution time. In contrast, the effect of the background gradient on echoes produced using a 1.6 gauss/cm static field gradient is calculated to be significant. This analysis suggests that RF gradient-echoes can produce MR images in which signal intensity is a function of the self-diffusion coefficient of water, but is not significantly affected by background gradients.     

MR imaging of apparent 3He gas transport in narrow pipes and rodent airways

High sensitivity makes hyperpolarized ³He an attractive signal source for visualizing gas flow with magnetic resonance (MR) imaging. Its rapid Brownian motion, however, can blur observed flow lamina and alter measured diffusion rates when excited nuclei traverse shear-induced velocity gradients during data acquisition. Here, both effects are described analytically, and predicted values for measured transport during laminar flow through a straight, 3.2-mm diameter pipe are validated using two-dimensional (2D) constant-time images of different binary gas mixtures. Results show explicitly how measured transport in narrow conduits is characterized by apparent values that depend on underlying gas dynamics and imaging time. In ventilated rats, this is found to obscure acquired airflow images. Nevertheless, flow splitting at airway branches is still evident and use of 3D vector flow mapping is shown to reveal surprising detail that highlights the correlation between gas dynamics and lung structure.     

3D MRI of non-Gaussian (3)He gas diffusion in the rat lung

In (3)He magnetic resonance images of pulmonary air spaces, the confining architecture of the parenchymal tissue results in a non-Gaussian distribution of signal phase that non-exponentially attenuates image intensity as diffusion weighting is increased. Here, two approaches previously used for the analysis of non-Gaussian effects in the lung are compared and related using diffusion-weighted (3)He MR images of mechanically ventilated rats. One approach is model-based and was presented by Yablonskiy et al., while the other approach utilizes the second order decay contribution that is predicted from the cumulant expansion theorem. Total lung coverage is achieved using a hybrid 3D pulse sequence that combines conventional phase encoding with sparse radial sampling for efficient gas usage. This enables the acquisition of nine 3D images using a total of only approximately 1 L of hyperpolarized (3)He gas. Diffusion weighting ranges from 0 s/cm(2) to 40 s/cm(2). Results show that the non-Gaussian effects of (3)He gas diffusion in healthy rat lungs are directly attributed to the anisotropic geometry of lung microstructure as predicted by the Yablonskiy model, and that quantitative analysis over the entire lung can be reliably repeated in time-course studies of the same animal.     

Relaxation and diffusion of perfluorocarbon gas mixtures with oxygen for lung MRI

We report measurements of free diffusivity D(0) and relaxation times T(1) and T(2) for pure C(2)F(6) and C(3)F(8) and their mixtures with oxygen. A simplified relaxation theory is presented and used to fit the data. The results enable spatially localized relaxation time measurements to determine the local gas concentration in lung MR images, so the free diffusivity D(0) is then known. Comparison of the measured diffusion to D(0) will express the extent of diffusion restriction and allow the local surface-to-volume ratio to be found.     

MR detection of mechanical vibrations using a radiofrequency field gradient

A new method for NMR characterization of mechanical waves, based upon radiofrequency field gradient for motion encoding, is proposed. A binomial B1 gradient excitation scheme was used to visualize the mobile spins undergoing a periodic transverse mechanical excitation. A simple model was designed to simulate the NMR signal as a function of the wave frequency excitation and the periodicity of the NMR pulse sequence. The preliminary results were obtained on a gel phantom at low vibration frequencies (0-200 Hz) by using a ladder-shaped coil generating a nearly constant RF field gradient along a specific known direction. For very small displacements and/or B1 gradients, the NMR signal measured on a gel phantom was proportional to the vibration amplitude and the pulse sequence was shown to be selective with respect to the vibration frequency. A good estimation of the direction of vibrations was obtained by varying the angle between the motion direction and the B1 gradient. The method and its use in parallel to more conventional MR elastography techniques are discussed. The presented approach might be of interest for noninvasive investigation of elastic properties of soft tissues and other materials.     

Assessment of in vitro vs. in vivo lung structure using hyperpolarized helium-3 diffusion magnetic resonance imaging

The purpose of this study was to assess the properties of a model system for hyperpolarized He-3 (HHe) diffusion MR imaging created from the lungs of New Zealand white rabbits by drying the lungs while inflated at constant pressure. The dried lungs were prepared by sacrificing the animal, harvesting the lungs en bloc and dehydrating the lungs for several days using dry compressed air. In four rabbits, the apparent diffusion coefficient (ADC) of HHe gas was measured in vivo and, within 1 week, in vitro in the dried lungs. To assess long-term repeatability, in vitro ADC values were measured again 3 months later. Dried lungs from four additional rabbits were imaged twice on the same day to assess the short-term repeatability of ADC measurements, and tissue samples from these lungs were then removed for histology. In vivo and in vitro ADC maps showed similar features and similar distributions of ADC values; mean in vivo and in vitro ADC values differed by less than 12%. The in vitro mean ADC values were highly reproducible, with no more than 5% difference between measurements for the short-term repeatability and less than 17% difference between measurements for the long-term repeatability. Histological samples from the dried lungs demonstrated that the lung structure remained intact. These results suggest that the dried lungs are a useful and inexpensive alternative to human or in vivo animal studies for HHe diffusion MR sequence development, testing and optimization.     

一种改进的核磁共振测井测量脉冲序列

提出了一种可以提高核磁共振测井测量速度的新型测量脉冲序列,它基于对原有CPMG测量脉冲序列的改进.在原有的测量序列中,在一个CPMG脉冲序列结束之后,直接进行等待(T₁延时),以使磁化强度矢量能够完全恢复到平衡位置,然后再进行下一个CPMG脉冲序列的测量.由于测井仪器的连续运动,在测量含有长T₁成分的地层时,等待时间需要取得较长,使得测井速度很低.但在实际地层中,恰恰又是油气的T₁较长,这无疑给核磁共振测井带来了麻烦.作者发现,如果在一个CPMG脉冲序列结束之后在-x方向先施加一个90°脉冲(或在x方向施加一个270°脉冲),再进行等待,就能更好地使长T₁成分得到极化,从而缩短等待时间、提高测量速度.文中对这种改进后的脉冲序列进行了计算和分析,证明它确实具有一定的优越性.     

Anomalous diffusion expressed through fractional order differential operators in the Bloch-Torrey equation

Diffusion weighted MRI is used clinically to detect and characterize neurodegenerative, malignant and ischemic diseases. The correlation between developing pathology and localized diffusion relies on diffusion-weighted pulse sequences to probe biophysical models of molecular diffusion-typically exp[-(bD)]-where D is the apparent diffusion coefficient (mm(2)/s) and b depends on the specific gradient pulse sequence parameters. Several recent studies have investigated the so-called anomalous diffusion stretched exponential model-exp[-(bD)(alpha)], where alpha is a measure of tissue complexity that can be derived from fractal models of tissue structure. In this paper we propose an alternative derivation for the stretched exponential model using fractional order space and time derivatives. First, we consider the case where the spatial Laplacian in the Bloch-Torrey equation is generalized to incorporate a fractional order Brownian model of diffusivity. Second, we consider the case where the time derivative in the Bloch-Torrey equation is replaced by a Riemann-Liouville fractional order time derivative expressed in the Caputo form. Both cases revert to the classical results for integer order operations. Fractional order dynamics derived for the first case were observed to fit the signal attenuation in diffusion-weighted images obtained from Sephadex gels, human articular cartilage and human brain. Future developments of this approach may be useful for classifying anomalous diffusion in tissues with developing pathology.     

Diffusion and relaxation effects in general stray field NMR experiments

We analyze the evolution of magnetization following any series of radiofrequency pulses in strongly inhomogeneous fields, with particular attention to diffusion and relaxation effects. When the inhomogeneity of the static magnetic field approaches or exceeds the strength of the RF field, the magnetization has contributions from different coherence pathways. The diffusion or relaxation induced decay of the signal amplitude is in general nonexponential, even if the sample has single relaxation times T(1), T(2) and a single diffusion coefficient D. In addition, the shape of the echo depends on diffusion and relaxation. It is possible to separate contributions from different coherence pathways by phase cycling of the RF pulses. The general analysis is tested on stray field measurements using two different pulse sequences. We find excellent agreement between measurements and calculations. The inversion recovery sequence is used to study the relaxation effects. We demonstrate two different approaches of data analysis to extract the relaxation time T(1). Finite pulse width effects on the timing of the echo formation are also studied. Diffusion effects are analyzed using the Carr--Purcell--Meiboom--Gill sequence. In a stray field of a constant gradient g, we find that unrestricted diffusion leads to nonexponential signal decay versus echo number N, but within experimental error the diffusion attenuation is still only a function of g(2)Dt(3)(E)N, where t(E) is the echo spacing.     

Rapid surface-to-volume ratio and tortuosity measurement using Difftrain

Analysis of diffusion measurements as a function of observation time (Delta), to calculate surface-to-volume ratios (S/V) and tortuosities (kappa), is a useful tool in the characterisation of porous media using NMR. However, using conventional pulsed field gradient (PFG) measurements, this requires long total experiment times (typically hours). Here, we show how the rapid diffusion measurement pulse sequence, Difftrain, can be used to provide the required experimental data much more rapidly (typically within minutes) with a consequential reduction in total experiment time of typically over an order of magnitude. Several novel modifications to the Difftrain pulse sequence are also presented to tailor it to this particular application; these include a variable delay between echoes (to ensure optimal echo position with respect to Delta) and a variable tip angle for the refocusing pulse (to ensure optimal use of available signal). Difftrain is applied to measure both S/V and kappa for a model glass bead pack; excellent agreement is found with both a conventional PFG measurement and with a bulk gravimetric measurement of S/V.     

Dynamic NMR q-space studies of microstructure with the multigrade CPMG sequence

A new approach to q-space studies of microstructure is proposed, which exploits the combined information contained in the water proton transverse relaxation time distribution and the frequency dependence of the apparent water diffusivity in heterogeneous systems. Using an automated two-dimensional multigrade CPMG sequence, both the pulse spacing and the amplitude of the applied field gradient are varied systematically and used to measure the frequency and wave vector dependence of the multiple exponential echo decay constants and amplitudes. Undesirable crossterms in the applied and background field gradients are eliminated by a simple procedure involving a sign reversal in the applied gradient. Nonlinear, local susceptibility-induced field gradients are shown to lead to enhanced, frequency-dependent apparent water diffusivities that are sensitive to the local microstructure.     

A Pulsed Field Gradient Spin-Echo Method for Diffusion Measurements in the Presence of Internal Gradients

Over the past decade several pulsed field gradient stimulated-echo methods have been presented for diffusion measurements in heterogeneous media. These methods have reduced or eliminated the coupling between the applied magnetic field gradient and a constant internal magnetic field gradient caused by susceptibility changes throughout the sample. For many research purposes thez-storage delay between the second and third π/2 RF pulse has been included in order to increase the decay of the echo attenuation to an appropriate level and to increase the signal-to-noise ratio by avoidingT₂relaxation of the magnetization in parts of the pulse sequence. For these reasons a stimulated-echo method has been applied instead of a spin-echo method. When studying systems where it is necessary to keep the duration of the pulse sequence at a minimum, and one is not dependent on usingz-storage time to increase the echo attenuation or to study diffusion as a function of observation time, a spin-echo method should be chosen. Here we propose a bipolar pulsed field gradient spin-echo method which is well suited to this purpose, and preliminary diffusion measurements are presented as illustration.     

MR diffusion - "diffraction" phenomenon in multi-pulse-field-gradient experiments

Using pulsed-field-gradient (PFG) experiments, the sizes of the pores in ordered porous media can be estimated from the "diffraction" pattern that the signal attenuation curves exhibit. A different diffraction pattern is observed when the experiment is extended to a larger number (N) of diffusion gradient pulse pairs. Simulations to calculate signal values from arbitrary gradient waveforms are performed for diffusion in restricted geometries using a matrix operator formalism. The simulations suggest that the differences in the characteristics of the attenuation curves are expected to make it possible to measure smaller pore sizes, to improve the accuracy of pore size measurements and potentially to distinguish different pore shapes using the N-PFG technique. Moreover, when an even number of PFG pairs is used, it is possible to observe the diffraction pattern at shorter diffusion times and measure an approximation to the average pore size even when the sample contains pores with a broad distribution of sizes.