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.     

³He lung morphometry technique: accuracy analysis and pulse sequence optimization

The (3)He lung morphometry technique (Yablonskiy et al., JAP, 2009), based on MRI measurements of hyperpolarized gas diffusion in lung airspaces, provides unique information on the lung microstructure at the alveolar level. 3D tomographic images of standard morphological parameters (mean airspace chord length, lung parenchyma surface-to-volume ratio, and the number of alveoli per unit lung volume) can be created from a rather short (several seconds) MRI scan. These parameters are most commonly used to characterize lung morphometry but were not previously available from in vivo studies. A background of the (3)He lung morphometry technique is based on a previously proposed model of lung acinar airways, treated as cylindrical passages of external radius R covered by alveolar sleeves of depth h, and on a theory of gas diffusion in these airways. The initial works approximated the acinar airways as very long cylinders, all with the same R and h. The present work aims at analyzing effects of realistic acinar airway structures, incorporating airway branching, physiological airway lengths, a physiological ratio of airway ducts and sacs, and distributions of R and h. By means of Monte-Carlo computer simulations, we demonstrate that our technique allows rather accurate measurements of geometrical and morphological parameters of acinar airways. In particular, the accuracy of determining one of the most important physiological parameter of lung parenchyma - surface-to-volume ratio - does not exceed several percent. Second, we analyze the effect of the susceptibility induced inhomogeneous magnetic field on the parameter estimate and demonstrate that this effect is rather negligible at B(0) ≤ 3T and becomes substantial only at higher B(0) Third, we theoretically derive an optimal choice of MR pulse sequence parameters, which should be used to acquire a series of diffusion-attenuated MR signals, allowing a substantial decrease in the acquisition time and improvement in accuracy of the results. It is demonstrated that the optimal choice represents three not equidistant b-values: b(1)=0, b(2)～2 s/cm(2), b(3)～8 s/cm(2).     

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.     

NMR Imaging of Thermally Polarized Helium-3 Gas

It is shown that thermally polarized 3He gas can be used to measure important physical parameters and to design, test, and tune imaging sequences. The bulk values of T1, T2, and the diffusion coefficient were measured in a glass cell containing a mixture of helium-3 (0.8 bar) and oxygen (0.2 bar). They were found to be T1 = 7 s, T2 = 2.4 s, and D = 1.6 cm2 s(-1). The relaxation times T2* and T1 and the apparent diffusion coefficient of thermally polarized helium-3 gas were measured in the rat lung, and these parameters were used to design a helium-3 optimized multi-spin-echo sequence which was shown to increase the signal-to-noise ratio sufficiently to obtain the first NMR-images of thermally polarized helium-3 in the rat lung.     

T2-shortening of 3He gas by magnetic microspheres

In a gas-filled material like the lung parenchyma, the transverse relaxation time (T2) for 3He is shortened by the deposition of magnetic microspheres and rapid molecular diffusion through induced field distortions. Here, this unique relaxation process is described theoretically and predicted T2-shortening is validated using pressurized 3He gas in a foam model of alveolar airways. Results demonstrate that: (1) significant T2-shortening is induced by microsphere deposition, (2) shortened 3He T2s are accurately predicted, and (3) measured relaxation times are exploitable for quantifying local deposition patterns. Based on these findings the feasibility of imaging inhaled particulates in vivo with hyperpolarized 3He is examined and performance projections are formulated.     

Effective interactions in dilute mixtures of3He in4He

Nonlocal pseudopotentials which describe the effective interaction between³He quasiparticles, and between these quasiparticles and the background⁴He liquid, are obtained as a function of concentration and pressure by generalizing the Aldrich-Pines pseudopotentials for pure³He and⁴He to dilute mixtures. The hierarchy of physical effects which determine these pseudopotentials is established. Interaction-induced short-range correlations are the dominant physical feature; next in order of importance is the greater zero point motion associated with the replacement of a⁴He atom by a³He atom, while spin-induced Pauli principle correlations play a significantly smaller, albeit still important role. We find a consistent trend in the change of the effective direct quasiparticle interactions with increasing concentration, and show how the Aldrich-Pines pseudopotentials for pure³He quasiparticles represent a natural extension of our results for dilute mixtures. Our calculated nonlocal pseudopotential for³He quasiparticles is qualitatively similar to that proposed by Bardeen, Baym, and Pines; it changes sign at somewhat lower momentum transfers than the BBP result, varies little with concentration, and provides a physical basis for understanding the BBP result. The effective interaction between quasiparticles of parallel spin, here determined for the first time, is essentially repulsive in the very dilute limit; as the concentration increases, it becomes increasingly attractive at low momentum transfers, and resembles closely that between antiparallel spin quasiparticles at 5% concentration. The concentration-dependent transport properties calculated from these pseudopotentials (which involve only one phenomenological parameter) are in good agreement with experiment at saturated vapor pressure (SVP), 10 atm, and 20 atm. Maxima in the thermal conductivity and spin diffusion are predicted to occur at concetrations somewhat less than 4%. Because the effective quasiparticle interactions are somewhat more repulsive than those previously proposed, we find the transition of the³He quasiparticles to the superfluid state takes place at significantly lower temperatures than many previous estimates; our predicted maximum superfluid transition temperature is 2×10^–8 K (for a 0.6% mixture at 20 atm).     

Single-Shot Diffusion Measurement in Laser-Polarized Gas

A single-shot pulsed gradient stimulated echo sequence is introduced to address the challenges of diffusion measurements of laser polarized 3He and 129Xe gas. Laser polarization enhances the NMR sensitivity of these noble gases by >10(3), but creates an unstable, nonthermal polarization that is not readily renewable. A new method is presented which permits parallel acquisition of the several measurements required to determine a diffusive attenuation curve. The NMR characterization of a sample's diffusion behavior can be accomplished in a single measurement, using only a single polarization step. As a demonstration, the diffusion coefficient of a sample of laser-polarized 129Xe gas is measured via this method.     

Zero-temperature equation of state of two-dimensional 3He

The equation of state of two-dimensional 3He at zero temperature has been calculated using the diffusion Monte Carlo method. By means of a combination of the fixed-node and released-node techniques, it is shown that backflow correlations provide a very accurate equation of state. The results prove unambiguously the non-self-bound character of two-dimensional 3He due to its Fermi statistics. We present solid evidence that the gas phase, predicted for the two-dimensional system, can be extrapolated to the case of 3He adsorbed on a strong substrate such as graphite.     

In vivo NMR of hyperpolarized 3He in the human lung at very low magnetic fields

We present NMR measurements of the diffusion of hyperpolarized 3He in the human lung performed at fields much lower than those of conventional MRI scanners. The measurements were made on standing subjects using homebuilt apparatus operating at 3mT. O(2)-limited transverse relaxation (T(2) up to 15-35s) could be measured in vivo. Accurate global diffusion measurements have been performed in vivo and in a plastic bag; the average apparent diffusion coefficient (ADC) in vivo was 14.2+/-0.6mm(2)/s, whereas the diffusion coefficient in the bag (3He diluted in N(2)) was 79.5+/-1mm(2)/s. 1D ADC mapping with high SNR ( approximately 200-300) demonstrates the real possibility of performing quality lung imaging at extremely low fields.     

Measurements and modeling of long range 3He diffusion in the lung using a "slice-washout" method

In healthy lung tissue, pulsed-gradient-spin-echo (PGSE) methods reveal apparent diffusion coefficients (ADC) of the order 0.20 cm2 s(-1); for diffusion times of approximately 2 ms. For these short diffusion times the ADC is only sensitive to structures approximately (2Dt)1/2 approximately 0.6mm in size. Recent work, using magnetic tagging of the longitudinal magnetization has revealed much smaller ADC values for longer length scales. In this work, the in vivo ADC from within the air-spaces, was measured using a new technique. The signal from a series of images was analyzed from a slice that was repeatedly imaged. Diffusion tends to "top-up" the non-renewable polarization within the slice, which leads to a non-exponential decay in image signal. Image data were compared to 1D finite-difference simulations of diffusion to calculate a long range ADC value. The results yield values of the order 0.034 cm2 s(-1), which are nearly an order of magnitude smaller than those reported by PGSE measurements at shorter diffusion times.     

Investigating 3He diffusion NMR in the lungs using finite difference simulations and in vivo PGSE experiments

Finite difference simulations have been used to model (3)He gas diffusion in simulated lung tissue. The technique has the advantage that a wide range of structural models and diffusion-sensitizing gradient waveforms can be investigated, for which analytical methods would otherwise be virtually impossible. Results from simulations and in vivo pulsed-gradient-spin-echo (PGSE) experiments show that the apparent diffusion coefficient (ADC) is a function of diffusion time and gradient strength, and suggests diffusion is locally anisotropic. The simulations have been compared to recent work on an analytical model that characterizes lung tissue as a series of independent cylinders. The results presented may have clinical implications for (3)He ADC measurements in assessing lung diseases such as chronic-obstructive-pulmonary-disease.     

Increasing the pump-up rate to polarize He gas using spin-exchange optical pumping method

In recent years, polarized ³He gas has increasingly been used as neutron polarizers and polarization analyzers. Two of the leading methods to polarize the ³He gas are the spin-exchange optical pumping (SEOP) method and the meta-stable exchange optical pumping (MEOP) method. At present, the SEOP setup is comparatively compact due to the fact that it does not require the sophisticated compressor system used in the MEOP method. The temperature and the laser power available determine the speed, at which the SEOP method polarizes the ³He gas. For the quantity of gas typically used in neutron scattering work, this speed is independent of the quantity of the gas required, whereas the polarizing time using the MEOP method is proportional to the quantity of gas required. Currently, using the SEOP method to polarize several bar-liters of ³He to 70% polarization would require 20−40 h. This is an order of magnitude longer than the MEOP method for the same quantity of gas and polarization. It would therefore be advantageous to speed up the SEOP process. In this article, we analyze the requirements for temperature, laser power, and the type of alkali used in order to shorten the time required to polarize ³He gas using the SEOP method.     

A wide angle neutron spin filter system using polarized He

Polarized ³He neutron spin filters can operate over a wide neutron energy range and provide a large angular acceptance. A compact ³He neutron spin filter system has been developed for the Multi-Axis Crystal Spectrometer at the National Institute of Standards and Technology (NIST) Center for Neutron Research. Sealed ³He cells, polarized by spin-exchange optical pumping, are used as polarizer and analyzer. The polarization of the neutrons incident on the sample is inverted by flipping the polarization of the ³He gas in the polarizer, with only a small effect on the analyzer cells. The cell fabrication process, ³He spin flipper, and the holding magnetic field are discussed and we present the results of a first on-linetest.     

Phase-aligned multiple spin-echo averaging: a simple way to improve signal-to-noise ratio of in vivo mouse spinal cord diffusion tensor image

Purpose

To improve signal-noise-ratio of in vivo mouse spinal cord diffusion tensor imaging using-phase aligned multiple spin-echo technique.

Material and methods

In vivo mouse spinal cord diffusion tensor imaging maps generated by multiple spin-echo and conventional spin-echo diffusion weighting were examined to demonstrate the efficacy of multiple spin-echo diffusion sequence to improve image quality and throughput. Effects of signal averaging using complex, magnitude and phased images from multiple spin-echo diffusion weighting were also assessed. Bayesian probability theory was used to generate phased images by moving the coherent signals to the real channel to eliminate the effect of phase variation between echoes while preserving the Gaussian noise distribution. Signal averaging of phased multiple spin-echo images potentially solves both the phase incoherence problem and the bias of the elevated Rician noise distribution in magnitude image. The proposed signal averaging with Bayesian phase-aligned multiple spin-echo images approach was compared to the conventional spin-echo data acquired with doubling the scan time. The diffusion tensor imaging parameters were compared in the mouse contusion spinal cord injury. Significance level (p-value) and effect size (Cohen’s d) were reported between the control and contused spinal cord to inspect the sensitivity of each approach in detecting white matter pathology.

Results

Compared to the spin-echo image, the signal-noise-ratio increased to 1.84-fold using the phased image averaging and to 1.30-fold using magnitude image averaging in the spinal cord white matter. Multiple spin-echo phased image averaging showed improved image quality of the mouse spinal cord among the tested methods. Diffusion tensor imaging metrics obtained from multiple spin-echo phased images using three echoes and two averages closely agreed with those derived by spin-echo magnitude data with four averages (two times more in acquisition time). The phased image averaging correctly reflected pathological features in contusion spinal cord injury.

Conclusion

Our in vivo imaging results indicate that averaging the phased multiple spin-echo images yields an 84% signal-noise-ratio increase over the spin-echo images and a 41% gain over the magnitude averaged multiple spin-echo images with equal acquisition time. Current results from the animal model of spinal cord injury suggest that the phased multiple spin-echo images could be used to improve signal-noise-ratio.     

Applications of He neutron spin filters at the NCNR

At the NIST Center for Neutron Research (NCNR), we have applied ³He neutron spin filters (NSFs) to the instruments where ³He NSFs are advantageous, such as thermal triple-axis spectrometry, small-angle neutron scattering, and diffuse reflectometry. We present the status of our development and application of this method, including polarized gas production by spin-exchange optical pumping, magnetostatic cavities for storage of the polarized gas on the beam line, and nuclear magnetic resonance (NMR)-based, on-line monitoring and reversal of the ³He polarization. We present the status of developing user-friendly interfaces incorporated into the instrument software to handle these ³He neutron spin filters while taking data and performing data analysis. Finally we discuss the status of development of a polarization capability on the multi-axis crystal spectrometer, which requires polarization analysis over a 220° angular range.     

An efficient diffusion approximation for 3D radiative transfer parameterization: application to cloudy atmospheres

The three-dimensional (3D) diffusion radiative transfer equation, which utilizes a four-term spherical harmonics expansion for the scattering phase function and intensity, has been efficiently solved by using the full multigrid numerical method. This approach can simulate the transfer of solar and thermal infrared radiation in inhomogeneous cloudy conditions with different boundary conditions and sharp boundary discontinuity. The correlated k-distribution method is used in this model for incorporation of the gaseous absorption in multiple-scattering atmospheres for the calculation of broadband fluxes and heating rates in the solar and infrared spectra. Comparison of the results computed from this approach with those computed from plane-parallel and 3D Monte Carlo models shows excellent agreement. This 3D radiative transfer approach is well suited for radiation parameterization involving 3D and inhomogeneous clouds in climate models.     

超极化Xenon 对慢阻肺的可视化加权成像

超极化气体3He 或者¹²⁹Xe 扩散加权成像已经被证明了能够有效检测慢性阻塞性肺部疾病(COPD)中肺部微结构的改变．相比于³He,¹²⁹Xe 更便宜而且更容易获得,但是¹²⁹Xe 成像中较低的信噪比致使¹²⁹Xe 的肺部表面扩散系数(ADC)的测量面临着许多困难．在该研究中,为了得到更高的图像信噪比,作者对气球模型,健康大鼠和COPD大鼠进行了单个b 值(14 cm2/s)的扩散加权超极化¹²⁹Xe 磁共振成像(MRI)．所有的COPD模型大鼠是通过烟熏和注射内毒素(LPS)进行诱导得到的．在7 T 磁共振成像仪上面获得了大鼠肺实质的超极化¹²⁹Xe ADC 值分布图．COPD 大鼠肺实质的¹²⁹Xe ADC 值是0.044 22±0.002 9 和0.042 34±0.002 3 cm2/s (Δ = 0.8/1.2 ms),远大于健康大鼠肺实质的¹²⁹Xe ADC 值0.037 7±0.002 3 和0.036 7±0.001 3 cm2/s．而且COPD 大鼠肺实质相关的¹²⁹Xe ADC 直方图也表现出了一定的展宽．这些结果说明了COPD 大鼠肺泡空腔的增大能够通过129Xe 在肺里面的ADC 增长和相关直方图的拓宽反应出来,从而证明了单个b 值的扩散加权MRI 方法可以有效地对COPD 大鼠进行检测．     

An MR imaging method for simultaneous measurement of gaseous diffusion constant and longitudinal relaxation time

A magnetic resonance imaging method for simultaneous and accurate determination of gaseous diffusion constant and longitudinal relaxation time is presented. The method is based on direct observation of diffusive motion. Initially, a slice-selective saturation of helium-3 (3He) spins was performed on a 3He/O2 phantom (9 atm/2 atm). A time-delay interval was introduced after saturation, allowing spins to diffuse in and out of the labeled slice. Following the delay interval a one-dimensional (1-D) projection image of the phantom was acquired. A series of 21 images was collected, each subsequent image having been acquired with an increased delay interval. Gradual spreading of the slice boundaries due to diffusion was thus observed. The projection profiles were fit to a solution of the Bloch equation corrected for diffusive motion. The fitting procedure yielded a value of D3He = 0.1562+/-0.0013 cm2/s, in good agreement with a measurement obtained with a modified version of the standard pulsed-field gradient technique. The method also enabled us to accurately measure the longitudinal relaxation of 3He spins by fitting the change of the total area under the projection profiles to an exponential. A value of T1 = 1.67 s (2 T field) was recorded, in excellent agreement with an inversion recovery measurement.     

旋风分离器固体颗粒浓度三维电容层析成像分布

本文采用12电极电容传感器测量得到旋风分离器锥体部分的电容信号,利用线性逆反推算法得到断面固体颗粒的二维浓度分布。基于二维断面浓度分布,采用电容三维图像重建算法,得到锥体部分的三维颗粒浓度分布。图像重建结果与旋风分离器实际运行状况和数值计算结果吻合,说明电容层析成像及三维图像重建方法的准确性和可行性。