NMR-based diffusion pore imaging

Diffusion pore imaging in the presence of extraporal water

Diffusion-weighted imaging (DWI) is a powerful, non-invasive tool which is widely used in clinical routine. Mostly, apparent diffusion coefficient maps are acquired, which cannot be related directly to cellular structure. More recently it was shown that DWI is able to reconstruct pore shapes using a specialized magnetic field gradient scheme so that cell size distributions may be obtained. So far, artificial systems have been used for experimental demonstration without extraporal signal components and relatively low gradient amplitudes. The aim of this study was to investigate the feasibility of diffusion pore imaging in the presence of extraporal fluids and to develop correction methods for the effects arising from extraporal signal contributions. Monte Carlo simulations and validation experiments on a 14.1 T NMR spectrometer equipped with a dedicated diffusion probe head were performed. Both by using a filter gradient approach suppressing extraporal signal components as well as by using post-processing methods relying on the Gaussian phase approximation, it was possible to reconstruct pore space functions in the presence of extraporal fluids with little to no deviations from the expectations. These results may be a significant step towards application of diffusion pore imaging to biological samples.

Double diffusion encoding and applications for biomedical imaging

Diffusion Magnetic Resonance Imaging (dMRI) is one of the most important contemporary non-invasive modalities for probing tissue structure at the microscopic scale. The majority of dMRI techniques employ standard single diffusion encoding (SDE) measurements, covering different sequence parameter ranges depending on the complexity of the method. Although many signal representations and biophysical models have been proposed for SDE data, they are intrinsically limited by a lack of specificity. Advanced dMRI methods have been proposed to provide additional microstructural information beyond what can be inferred from SDE. These enhanced contrasts can play important roles in characterizing biological tissues, for instance upon diseases (e.g. neurodegenerative, cancer, stroke), aging, learning, and development. In this review we focus on double diffusion encoding (DDE), which stands out among other advanced acquisitions for its versatility, ability to probe more specific diffusion correlations, and feasibility for preclinical and clinical applications. Various DDE methodologies have been employed to probe compartment sizes (Section 3), decouple the effects of microscopic diffusion anisotropy from orientation dispersion (Section 4), probe displacement correlations, study exchange, or suppress fast diffusing compartments (Section 6). DDE measurements can also be used to improve the robustness of biophysical models (Section 5) and study intra-cellular diffusion via magnetic resonance spectroscopy of metabolites (Section 7). This review discusses all these topics as well as important practical aspects related to the implementation and contrast in preclinical and clinical settings (Section 9) and aims to provide the readers a guide for deciding on the right DDE acquisition for their specific application.

Stimulated echo double diffusion encoded imaging of closed pores: Influence and removal of unbalanced terms

Nuclear magnetic resonance (NMR) diffusion pore imaging has been proposed to study the shape of arbitrary closed pores filled with an NMR-detectable medium by use of nonclassical diffusion encoding schemes. Potential applications can be found in biomedical imaging and porous media research. When studying non-point-symmetric pores, NMR signals with nonvanishing imaginary parts arise containing the pore shape information, which is lost for classical diffusion encoding schemes. Key limitations are the required high magnetic field gradient amplitudes and T2 relaxation while approaching the diffusion long-time limit. To benefit from the slower T1 decay, we demonstrate the feasibility of diffusion pore imaging with stimulated echoes using Monte Carlo simulations and experiments with hyperpolarized xenon-129 gas in well-defined geometries and show that the necessary complex-valued signals can be acquired. Analytical derivation of the stimulated echo double diffusion encoded signal was performed to investigate the effect of the additionally arising undesired terms on the complex phase information. These terms correspond to signals arising for spin-echo sequences with unbalanced gradients. For most possible applications, the unbalanced terms can be neglected. If non-negligible, selection of the appropriate signal component using a phase cycling scheme was demonstrated experimentally. Using stimulated echoes may be a step towards application of diffusion pore imaging to larger pores with gradient amplitudes available today in preclinical systems.

Multimodal integration of diffusion MRI for better characterization of tissue biology

The contrast in diffusion-weighted MR images is due to variations of diffusion properties within the examined specimen. Certain microstructural information on the underlying tissues can be inferred through quantitative analyses of the diffusion-sensitized MR signals. In the first part of the paper, we review two types of approach for characterizing diffusion MRI signals: Bloch's equations with diffusion terms, and statistical descriptions. Specifically, we discuss expansions in terms of cumulants and orthogonal basis functions, the confinement tensor formalism and tensor distribution models. Further insights into the tissue properties may be obtained by integrating diffusion MRI with other techniques, which is the subject of the second part of the paper. We review examples involving magnetic susceptibility, structural tensors, internal field gradients, transverse relaxation and functional MRI. Integrating information provided by other imaging modalities (MR based or otherwise) could be a key to improve our understanding of how diffusion MRI relates to physiology and biology.

Precise Identification and Analysis of Micro/Nano-Sized Pore Structure in Shale with Fe3O4/Au Hybrid Nanocomposite

Analysis and characterization of micro/nano-sized pore structure are critical issues in shale geology and engineering. Scanning electron microscopy (SEM) imaging is one of the most widespread methods for the analysis of the micro/nano-sized pores in shale, but precise identification of the ultrafine pore structure in shale is still a big challenge because shale is so complex that some components may have overlap with pores based on the simple discrimination of gray scale under SEM microscopy. Here, Fe₃O₄/Au nanocomposite with magnetic properties is synthesized, characterized, and introduced as a novel pore-marker to improve SEM identification and quantitation of micro/nano-sized pores in shale. Due to the superparamagnetic property, the nanomarker is conveniently controlled by an external magnetic field to fill into pores and offers a sharp contrast imaging between matrix of shale (various gray) and pores (bright), which makes the identification of micro/nano-sized pores in shale much more straightforward and reliable. Furthermore, because gold, as a noble metal, is particularly rare in shale, energy-dispersive X-ray spectroscopy mapping of Au is delicately used to precisely calculate area porosity in shale. Combining with the aforementioned merits of the nanomarker, a precise and practical technique is proposed to promote characterization of micro/nano-sized pores in shale.

Experimental determination of pore shapes using phase retrieval from q-space NMR diffraction

This paper presents an approach to solving the phase problem in nuclear magnetic resonance (NMR) diffusion pore imaging, a method that allows imaging the shape of arbitrary closed pores filled with an NMR-detectable medium for investigation of the microstructure of biological tissue and porous materials. Classical q-space imaging composed of two short diffusion-encoding gradient pulses yields, analogously to diffraction experiments, the modulus squared of the Fourier transform of the pore image which entails an inversion problem: An unambiguous reconstruction of the pore image requires both magnitude and phase. Here the phase information is recovered from the Fourier modulus by applying a phase retrieval algorithm. This allows omitting experimentally challenging phase measurements using specialized temporal gradient profiles. A combination of the hybrid input-output algorithm and the error reduction algorithm was used with dynamically adapting support (shrinkwrap extension). No a priori knowledge on the pore shape was fed to the algorithm except for a finite pore extent. The phase retrieval approach proved successful for simulated data with and without noise and was validated in phantom experiments with well-defined pores using hyperpolarized xenon gas.

Dynamics of local magnetization in the eigenbasis of the Bloch-Torrey operator

We consider diffusion within pores with general shapes in the presence of spatially linear magnetic field profiles. The evolution of local magnetization of the spin bearing particles can be described by the Bloch-Torrey equation. We study the diffusive process in the eigenbasis of the non-Hermitian Bloch-Torrey operator. It is possible to find expressions for some special temporal gradient waveforms employed to sensitize the nuclear magnetic resonance (NMR) signal to diffusion. For more general gradient waveforms, we derive an efficient numerical solution by introducing a novel matrix formalism. Compared to previous methods, this new approach requires a fewer number of eigenfunctions to achieve the same accuracy. This shows that these basis functions are better suited to the problem studied. The new framework could provide new important insights into the fundamentals of diffusion sensitization, which could further the development of the field of NMR.

Imaging the Pore Structure of Geological Materials with Bifunctional Nanoparticles

Analysis of complex pore structure of geomaterials is a fundamental issue in geoscience. Here bifunctional nanoparticles with magnetic and fluorescent properties are introduced as novel markers for optical imaging of pore structure in geomaterials. Using the paramagnetic property, powder of the nanoparticle is driven into pores under an external magnetic field, avoiding a tedious sample preparation and eliminating artificial damage of sample preparation in conventional methods. Meanwhile, the fluorescent nanoparticle marker offers a sharp contrast imaging between the rock matrix (black) and pores (bright) under microscopy. Furthermore, fluorescent nanoparticles with different sizes and colors are designed to demonstrate the potential of the method for describing pore throat sizes. Combining the merits of the paramagnetic and fluorescent properties of nanoparticles, a convenient and practical sample preparation is proposed to promote optical imaging analysis of the pore structure in geomaterials.

1D and 2D diffusion pore imaging on a preclinical MR system using adaptive rephasing: Feasibility and pulse sequence comparison

Diffusion pore imaging (DPI) has recently been proposed as a means to acquire images of the average pore shape in an image voxel or region of interest. The highly asymmetric gradient scheme of its sequence makes it substantially demanding in terms of the hardware of the NMR system. The aim of this work is to show the feasibility of DPI on a preclinical 9.4T animal scanner. Using water-filled capillaries with an inner radius of 10μm, four different variants of the DPI sequence were compared in 1D and 2D measurements. The pulse sequences applied cover the basic implementation using one long and one temporally narrow gradient pulse, a CPMG-like variant with multiple refocusing RF pulses as well as two variants splitting up the long gradient and distributing it on either side of the refocusing pulse. Substantial differences between the methods were found in terms of signal-to-noise ratio, contrast, blurring, deviations from the expected results and sensitivity to gradient imperfections. Each of the tested sequences was found to produce characteristic gradient mismatches dependent on the absolute value, direction and sign of the applied q-value. Read gradients were applied to compensate these mismatches translating them into time shifts, which enabled 1D DPI yielding capillary radius estimations within the tolerances specified by the manufacturer. For a successful DPI application in 2D, a novel gradient amplitude adaption scheme was implemented to correct for the occurring time shifts. Using this adaption, higher conformity to the expected pore shape, reduced blurring and enhanced contrast were achieved. Images of the phantom's pore shape could be acquired with a nominal resolution of 2.2μm.

Nuclear magnetic resonance diffusion pore imaging: Experimental phase detection by double diffusion encoding

Diffusion pore imaging is an extension of diffusion-weighted nuclear magnetic resonance imaging enabling the direct measurement of the shape of arbitrarily formed, closed pores by probing diffusion restrictions using the motion of spin-bearing particles. Examples of such pores comprise cells in biological tissue or oil containing cavities in porous rocks. All pores contained in the measurement volume contribute to one reconstructed image, which reduces the problem of vanishing signal at increasing resolution present in conventional magnetic resonance imaging. It has been previously experimentally demonstrated that pore imaging using a combination of a long and a narrow magnetic field gradient pulse is feasible. In this work, an experimental verification is presented showing that pores can be imaged using short gradient pulses only. Experiments were carried out using hyperpolarized xenon gas in well-defined pores. The phase required for pore image reconstruction was retrieved from double diffusion encoded (DDE) measurements, while the magnitude could either be obtained from DDE signals or classical diffusion measurements with single encoding. The occurring image artifacts caused by restrictions of the gradient system, insufficient diffusion time, and by the phase reconstruction approach were investigated. Employing short gradient pulses only is advantageous compared to the initial long-narrow approach due to a more flexible sequence design when omitting the long gradient and due to faster convergence to the diffusion long-time limit, which may enable application to larger pores.

Phase incremented echo train acquisition applied to magnetic resonance pore imaging

Efficient phase cycling schemes remain a challenge for NMR techniques if the pulse sequences involve a large number of rf-pulses. Especially complex is the Carr Purcell Meiboom Gill (CPMG) pulse sequence where the number of rf-pulses can range from hundreds to several thousands. Our recent implementation of Magnetic Resonance Pore Imaging (MRPI) is based on a CPMG rf-pulse sequence in order to refocus the effect of internal gradients inherent in porous media. While the spin dynamics for spin-1/2 systems in CPMG like experiments are well understood it is still not straight forward to separate the desired pathway from the spectrum of unwanted coherence pathways. In this contribution we apply Phase Incremented Echo Train Acquisition (PIETA) to MRPI. We show how PIETA offers a convenient way to implement a working phase cycling scheme and how it allows one to gain deeper insights into the amplitudes of undesired pathways.

Effective potential for magnetic resonance measurements of restricted diffusion

The signature of diffusive motion on the NMR signal has been exploited to characterize the mesoscopic structure of specimens in numerous applications. For compartmentalized specimens comprising isolated subdomains, a representation of individual pores is necessary for describing restricted diffusion within them. When gradient waveforms with long pulse durations are employed, a quadratic potential profile is identified as an effective energy landscape for restricted diffusion. The dependence of the stochastic effective force on the center-of-mass position is indeed found to be approximately linear (Hookean) for restricted diffusion even when the walls are sticky. We outline the theoretical basis and practical advantages of our picture involving effective potentials.

NMR-based diffusion lattice imaging

Nuclear magnetic resonance (NMR) diffusion experiments are widely employed as they yield information about structures hindering the diffusion process, e.g., about cell membranes. While it has been shown in recent articles that these experiments can be used to determine the shape of closed pores averaged over a volume of interest, it is still an open question how much information can be gained in open well-connected systems. In this theoretical work, it is shown that the full structure information of connected periodic systems is accessible. To this end, the so-called "SEquential Rephasing by Pulsed field-gradient Encoding N Time intervals" (SERPENT) sequence is used, which employs several diffusion encoding gradient pulses with different amplitudes. Two two-dimensional solid matrices that are surrounded by an NMR-visible medium are considered: a hexagonal lattice of cylinders and a rectangular lattice of isosceles triangles.

Symmetry of the gradient profile as second experimental dimension in the short-time expansion of the apparent diffusion coefficient as measured with NMR diffusometry

The time-dependent apparent diffusion coefficient as measured by pulsed gradient NMR can be used to estimate parameters of porous structures including the surface-to-volume ratio and the mean curvature of pores. In this work, the short-time diffusion limit and in particular the influence of the temporal profile of diffusion gradients on the expansion as proposed by Mitra et al. (1993) is investigated. It is shown that flow-compensated waveforms, i.e. those whose first moment is zero, are blind to the term linear in observation time, which is the term that is proportional to mean curvature and surface permeability. A gradient waveform that smoothly interpolates between flow-compensated and bipolar waveform is proposed and the degree of flow-compensation is used as a second experimental dimension. This two-dimensional ansatz is shown to yield an improved precision when characterizing the confining domain. This technique is demonstrated with simulations and in experiments performed with cylindrical capillaries of 100 μm radius.

Effects of pore-size and shape distributions on diffusion pore imaging by nuclear magnetic resonance

In medical imaging and porous media research, NMR diffusion measurements are extensively used to investigate the structure of diffusion restrictions such as cell membranes. Recently, several methods have been proposed to unambiguously determine the shape of arbitrary closed pores or cells filled with an NMR-visible medium by diffusion experiments. The first approach uses a combination of a long and a short diffusion-weighting gradient pulse, while the other techniques employ short gradient pulses only. While the eventual aim of these methods is to determine pore-size and shape distributions, the focus has been so far on identical pores. Thus, the aim of this work is to investigate the ability of these different methods to resolve pore-size and orientation distributions. Simulations were performed comparing the various pore imaging techniques employing different distributions of pore size and orientation and varying timing parameters. The long-narrow gradient profile is most advantageous to investigate pore distributions, because average pore images can be directly obtained. The short-gradient methods suppress larger pores or induce a considerable blurring. Moreover, pore-shape-specific artifacts occur; for example, the central part of a distribution of cylinders may be largely underestimated. Depending on the actual pore distribution, short-gradient methods may nonetheless yield good approximations of the average pore shape. Furthermore, the application of short-gradient methods can be advantageous to differentiate whether pore-size distributions or intensity distributions, e.g., due to surface relaxation, are predominant.

Characterizing magnetic resonance signal decay due to Gaussian diffusion: the path integral approach and a convenient computational method

The influence of Gaussian diffusion on the magnetic resonance signal is determined by the apparent diffusion coefficient (ADC) and tensor (ADT) of the diffusing fluid as well as the gradient waveform applied to sensitize the signal to diffusion. Estimations of ADC and ADT from diffusion-weighted acquisitions necessitate computations of, respectively, the b-value and b-matrix associated with the employed pulse sequence. We establish the relationship between these quantities and the gradient waveform by expressing the problem as a path integral and explicitly evaluating it. Further, we show that these important quantities can be conveniently computed for any gradient waveform using a simple algorithm that requires a few lines of code. With this representation, our technique complements the multiple correlation function (MCF) method commonly used to compute the effects of restricted diffusion, and provides a consistent and convenient framework for studies that aim to infer the microstructural features of the specimen.

Magnetic-resonance pore imaging of nonsymmetric microscopic pore shapes

Imaging of the microstructure of porous media such as biological tissue or porous solids is of high interest in health science and technology, engineering and material science. Magnetic resonance pore imaging (MRPI) is a recent technique based on nuclear magnetic resonance (NMR) which allows us to acquire images of the average pore shape in a given sample. Here we provide details on the experimental design, challenges, and requirements of MRPI, including its calibration procedures. Utilizing a laser-machined phantom sample, we present images of microscopic pores with a hemiequilateral triangular shape even in the presence of NMR relaxation effects at the pore walls. We therefore show that MRPI is applicable to porous samples without a priori knowledge about their pore shape and symmetry. Furthermore, we introduce "MRPI mapping," which combines MRPI with conventional magnetic resonance imaging (MRI). This enables one to resolve microscopic pore sizes and shapes spatially, thus expanding the application of MRPI to samples with heterogeneous distributions of pores.

Two-dimensional diffusion time correlation experiment using a single direction gradient

The time dependence of the diffusion coefficient is a well known property of porous media and commonly obtained by pulsed field gradient (PFG) NMR. In practical materials, its analysis can be complicated by the presence of a broad pore size distribution and multiple fluid phases with different diffusion coefficients. We propose a two-dimensional Diffusion Time Correlation experiment (DTC), which utilizes the double-PFG with a single-direction gradient to yield a two-dimensional correlation function of the diffusion coefficient for two different diffusion times. This correlation map separates out restricted diffusion from the bulk diffusion process and we demonstrate this on a plant and bulk water sample. In its development, we show that the d-PFG should then be thought of as correlating two apparent diffusion coefficients measured by two overlapping gradient waveforms.

Parallel acquisition of q-space using second order magnetic fields for single-shot diffusion measurements

A proof of concept is presented for the parallel acquisition of q-space under diffusion using a second order magnetic field. The second order field produces a gradient strength which varies in space, allowing a range of gradients to be applied in a single pulse, and q-space encoded into real space. With the use of a read gradient, the spatial information is regained from the NMR signal, and real space mapped onto q-space for a thin slice excitation volume. As the diffusion encoded image for a thin slice can be mapped onto q-space, and the average propagator is the inverse Fourier transform of the q-space data, it follows that the acquisition of the echo is a direct measurement of the average propagator. In the absence of a thin slice selection, the real space to q-space mapping is lost, but the ability to measure the diffusion coefficient retained with an increase in signal to noise.

Measuring small compartment dimensions by probing diffusion dynamics via Non-uniform Oscillating-Gradient Spin-Echo (NOGSE) NMR

Noninvasive measurements of microstructure in materials, cells, and in biological tissues, constitute a unique capability of gradient-assisted NMR. Diffusion-diffraction MR approaches pioneered by Callaghan demonstrated this ability; Oscillating-Gradient Spin-Echo (OGSE) methodologies tackle the demanding gradient amplitudes required for observing diffraction patterns by utilizing constant-frequency oscillating gradient pairs that probe the diffusion spectrum, D(ω). Here we present a new class of diffusion MR experiments, termed Non-uniform Oscillating-Gradient Spin-Echo (NOGSE), which dynamically probe multiple frequencies of the diffusion spectral density at once, thus affording direct microstructural information on the compartment's dimension. The NOGSE methodology applies N constant-amplitude gradient oscillations; N-1 of these oscillations are spaced by a characteristic time x, followed by a single gradient oscillation characterized by a time y, such that the diffusion dynamics is probed while keeping (N-1)x+y≡TNOGSE constant. These constant-time, fixed-gradient-amplitude, multi-frequency attributes render NOGSE particularly useful for probing small compartment dimensions with relatively weak gradients - alleviating difficulties associated with probing D(ω) frequency-by-frequency or with varying relaxation weightings, as in other diffusion-monitoring experiments. Analytical descriptions of the NOGSE signal are given, and the sequence's ability to extract small compartment sizes with a sensitivity towards length to the sixth power, is demonstrated using a microstructural phantom. Excellent agreement between theory and experiments was evidenced even upon applying weak gradient amplitudes. An MR imaging version of NOGSE was also implemented in ex vivo pig spinal cords and mouse brains, affording maps based on compartment sizes. The effects of size distributions on NOGSE are also briefly analyzed.

Diffusion pore imaging by hyperpolarized xenon-129 nuclear magnetic resonance

While NMR diffusion measurements are widely used to derive parameters indirectly related to the microstructure of biological tissues and porous media, direct imaging of pore shapes would be of high interest. Here we demonstrate experimentally that complexly shaped closed pores can be imaged by diffusion acquisitions. Collecting the signal from the whole sample eliminates the problem of vanishing signal at increasing resolution of conventional NMR imaging. This approach may be used to noninvasively obtain structural information inaccessible so far such as pore or cell shapes, cell density, or axon integrity.