Tortuosity measurement and the effects of finite pulse widths on xenon gas diffusion NMR studies of porous media.

We have extended the utility of NMR as a technique to probe porous media structure over length scales of approximately 100-2000 microm by using the spin 1/2 noble gas 129Xe imbibed into the system's pore space. Such length scales are much greater than can be probed with NMR diffusion studies of water-saturated porous media. We utilized Pulsed Gradient Spin Echo NMR measurements of the time-dependent diffusion coefficient, D(t), of the xenon gas filling the pore space to study further the measurements of both the pore surface-area-to-volume ratio, S/V(p), and the tortuosity (pore connectivity) of the medium. In uniform-size glass bead packs, we observed D(t) decreasing with increasing t, reaching an observed asymptote of approximately 0.62-0.65D(0), that could be measured over diffusion distances extending over multiple bead diameters. Measurements of D(t)/D(0) at differing gas pressures showed this tortuosity limit was not affected by changing the characteristic diffusion length of the spins during the diffusion encoding gradient pulse. This was not the case at the short time limit, where D(t)/D(0) was noticeably affected by the gas pressure in the sample. Increasing the gas pressure, and hence reducing D(0) and the diffusion during the gradient pulse served to reduce the previously observed deviation of D(t)/D(0) from the S/V(p) relation. The Pade approximation is used to interpolate between the long and short time limits in D(t). While the short time D(t) points lay above the interpolation line in the case of small beads, due to diffusion during the gradient pulse on the order of the pore size, it was also noted that the experimental D(t) data fell below the Pade line in the case of large beads, most likely due to finite size effects.     

Measurements of restricted diffusion using an oscillating gradient spin-echo sequence

An oscillating gradient spin-echo (OGSE) pulse sequence was used to measure the apparent diffusion coefficient (D(app)) of water in the short diffusion time regime in the presence of restrictions. The diffusion coefficients of water in a simple water sample and a water and oil mixture were measured to be the same for different periods of the gradient oscillation, as expected when there are no restriction effects. The D(app) of water in the spaces between closely packed beads was also measured as a function of the gradient oscillation periods in the range 11 to 80 ms. The D(app) of water in restricted systems varies with the period of the gradient oscillation and the dispersion depends on the scale of the restriction. For a sample of packed beads of diameter 9.1 +/- 0.7 microm, the pore surface-to-volume ratio was estimated experimentally by this method to be 1.3 +/- 0.1 microm(-1), corresponding to a mean pore diameter of 6.4 +/- 0.7 microm. A Monte Carlo computer simulation of the NMR OGSE signal from the spins diffusing in a system of compartments was also implemented and the D(app) demonstrated similar behavior with gradient oscillation periods.     

Spatial and temporal coarse graining for dispersion in randomly packed spheres

The propagator for molecular displacements P(zeta, t) and its first three cumulants were measured for Stokes flow in monodisperse bead packs with different sphere sizes d and molecular diffusion coefficients D(m). We systematically varied the normalized mean displacement /d and diffusion length L(D)=sqrt[2D(m)t]/d. The experimental results map onto each other with this scaling. For L(D)/d<0.2 the propagator remains non-Gaussian, and thus an advection diffusion equation is not obeyed, for mean displacements measured up to 10d. A Gaussian shape is approached for large mean displacements when L(D)>0.3d.     

MAG-PGSTE: a new STE-based PGSE NMR sequence for the determination of diffusion in magnetically inhomogeneous samples

A new stimulated echo based pulsed gradient spin-echo sequence, MAG-PGSTE, has been developed for the determination of self-diffusion in magnetically inhomogeneous samples. The sequence was tested on two glass bead samples (i.e., 212-300 and <106 microm glass bead packs). The MAG-PGSTE sequence was compared to the MAGSTE (or MPFG) (P.Z. Sun, J.G. Seland, D. Cory, Background gradient suppression in pulsed gradient stimulated echo measurements, J. Magn. Reson. 161 (2003) 168-173; P.Z. Sun, S.A. Smith, J. Zhou, Analysis of the magic asymmetric gradient stimulated echo sequence with shaped gradients, J. Magn. Reson. 171 (2004) 324-329; P.Z. Sun, Improved diffusion measurement in heterogeneous systems using the magic asymmetric gradient stimulated echo (MAGSTE) technique, J. Magn. Reson. 187 (2007) 177-183; P. Galvosas, F. Stallmach, J. K?rger, Background gradient suppression in stimulated echo NMR diffusion studies using magic pulsed field gradient ratios, J. Magn. Reson. 166 (2004) 164-173, P. Galvosas, PFG NMR-Diffusionsuntersuchungen mit ultra-hohen gepulsten magnetischen Feldgradienten an mikropor?sen Materialien, Ph.D. Thesis, Universit?t Leipzig, 2003, P.Z. Sun, Nuclear Magnetic Resonance Microscopy and Diffusion, Ph.D. Thesis, Massachusetts Institute of Technology, 2003] sequence and Cotts 13-interval [R.M. Cotts, M.J.R. Hoch, T. Sun, J.T. Marker, Pulsed field gradient stimulated echo methods for improved NMR diffusion measurements in heterogeneous systems, J. Magn. Reson. 83 (1989) 252-266] sequence using both glass bead samples. The MAG-PGSTE and MAGSTE (or MPFG) sequences outperformed the Cotts 13-interval sequence in the measurement of diffusion coefficients; more interestingly, for the sample with higher background gradients (i.e., the <106 microm glass bead sample), the MAG-PGSTE sequence provided higher signal-to-noise ratios and thus better diffusion measurements than the MAGSTE and Cotts 13-interval sequences. In addition, the MAG-PGSTE sequence provided good characterization of the surface-to-volume ratio for the glass bead samples.     

Probing porous media with gas diffusion NMR.

We show that gas diffusion nuclear magnetic resonance (GD-NMR) provides a powerful technique for probing the structure of porous media. In random packs of glass beads, using both laser-polarized and thermally polarized xenon gas, we find that GD-NMR can accurately measure the pore space surface-area-to-volume ratio, S/V rho, and the tortuosity, alpha (the latter quantity being directly related to the system's transport properties). We also show that GD-NMR provides a good measure of the tortuosity of sandstone and complex carbonate rocks.     

Measurement of diffusion coefficients in polystyrene using xenon-129 NMR

A sample of polystyrene beads, 18 μm in diameter, has been sealed in an NMR tube under 10 atm of xenon gas. Two dimensional,¹²⁹Xe NMR spectra show cross peaks between the resonances corresponding to xenon in the free gas and the sorbed state, indicating that appreciable exchange occurs during the mixing time of the NMR experiment. Selective saturation of the free gas resonance attenuates the integrated intensity of the sorbed xenon resonance as a function of saturation time, thus allowing the accurate measurement of the exchange rates between the gas and the sorbed states. A model has been developed using a slightly modified form of Crank’s treatment of diffusion in a sphere which allows for the accurate determination of the diffusion coefficient for xenon in the sorbed state. The diffusion coefficient for xenon in polystyrene at 25°C is determined to be 2.9·10^?9 cm²/s.     

Time-dependent velocities in porous media dispersive flow.

Pulsed Gradient Spin Echo (PGSE) NMR methods may be used to measure the asymptotic dispersion coefficient as well as the velocity autocorrelation function (VACF) in porous media flow. The VACF can be measured in the frequency domain using repetitive gradient pulse trains, and in the time domain using double PGSE encoding. The one dimensional double PGSE method, and the two dimensional velocity exchange experiment (VEXSY) are briefly outlined and their application to flow in monodisperse 0.5 mm diameter beads packs described, both axial and transverse VACFs being examined. The measured correlation times are shown to agree well with calculated values. The asymptotic dispersion coefficients agree with literature values in the case of transverse flow while in axial flow it is shown that asymptotic conditions are not achieved, even for observation times longer than the correlation time for flow around a bead.     

Internal magnetic gradient fields in glass bead packs from numerical simulations and constant time diffusion spin echo measurements

Internal magnetic field gradients in water saturated glass bead packs were studied by numerical simulations and a constant time spin echo (CTSE) experiment. The CTSE is comprised of two spin echo refocusing periods where each of the two evolution periods, tau1 and tau2, is varied so that the total evolution, 2(tau1 + tau2), is held constant. The experiment is similar to that introduced by Norwood and Quilter and allows the effects of dephasing due to diffusion in a magnetic field gradient to be separated from other relaxation mechanisms. In our experiments, the magnetic susceptibility difference between the pore fluid and glass beads creates the internal field gradient. CTSE measurements were performed at 7 T (300 MHz 1H) for water saturated in 50 microm diameter glass bead pack. We find that the internal gradients in the center of the pore bodies, where free diffusion applies, is in the range of 10 to 100 G/cm. This fluid volume accounts for < or =10% of the total pore volume. From direct numerical simulations of the internal magnetic field based on a first principles calculation, we find that the major fraction, >90%, of the pore volume has internal gradients of order 500 to 5,000 G/cm. Signals from water in these large gradients is not observed in our CTSE measurements.     

Pulsed-Field-Gradient Measurements of Time-Dependent Gas Diffusion

Pulsed-field-gradient NMR techniques are demonstrated for measurements of time-dependent gas diffusion. The standard PGSE technique and variants, applied to a free gas mixture of thermally polarized xenon and O₂, are found to provide a reproducible measure of the xenon diffusion coefficient (5.71 × 10⁻⁶m²s⁻¹for 1 atm of pure xenon), in excellent agreement with previous, non-NMR measurements. The utility of pulsed-field-gradient NMR techniques is demonstrated by the first measurement of time-dependent (i.e., restricted) gas diffusion inside a porous medium (a random pack of glass beads), with results that agree well with theory. Two modified NMR pulse sequences derived from the PGSE technique (named the Pulsed Gradient Echo, or PGE, and the Pulsed Gradient Multiple Spin Echo, or PGMSE) are also applied to measurements of time dependent diffusion of laser polarized xenon gas, with results in good agreement with previous measurements on thermally polarized gas. The PGMSE technique is found to be superior to the PGE method, and to standard PGSE techniques and variants, for efficiently measuring laser polarized noble gas diffusion over a wide range of diffusion times.     

From random walk to single-file diffusion

We report an experimental study of diffusion in a quasi-one-dimensional (q1D) colloid suspension which behaves like a Tonks gas. The mean squared displacement as a function of time is described well with an ansatz encompassing a time regime that is both shorter and longer than the mean time between collisions. The ansatz asserts that the inverse mean squared displacement is the sum of the inverse mean squared displacement for short time normal diffusion (random walk) and the inverse mean squared displacement for asymptotic single-file diffusion (SFD). The dependence of the 1D mobility in the SFD on the concentration of the colloids agrees quantitatively with that derived for a hard rod model, which confirms for the first time the validity of the hard rod SFD theory. We also show that a recent SFD theory by Kollmann leads to the hard rod SFD theory for a Tonks gas.     

Measuring pore connectivity by pulsed field gradient diffusion editing with hydrocarbon gases

Measurements of time-dependent diffusion are performed on a rock sample saturated first with water, then methane and finally ethane. The gases were selected because their increased diffusivities and relaxation times allow probing greater length scales than water and because of their practical relevance. The nuclear magnetic resonance measurements employed pulse field gradient diffusion editing pulse sequences, allowing analysis of D(t) as a function of relaxation time. Very different D(t) behaviors are observed for different relaxation times, including indications of connected pore networks at moderate relaxation times.     

Development of the liquid level meters for the PandaX dark matter detector

The two-phase xenon detector is at the frontier of dark matter direct search. This kind of detector uses liquid xenon as the sensitive target and is operated in two-phase (liquid/gas) mode, where the liquid level needs to be monitored and controlled in sub-millimeter precision. In this paper, we present a detailed design and study of two kinds of level meters for the PandaX dark matter detector. The long level meter is used to monitor the overall liquid level while short level meters are used to monitor the inclination of the detector. These level meters are cylindrical capacitors that are custom-made from two concentric metal tubes. Their capacitance values are read out by a universal transducer interface chip and are recorded by the PandaX slow control system. We present the developments that lead to level meters with long-term stability and sub-millimeter precision. Fluctuations (standard deviations) of less than 0.02 mm for the short level meters and less than 0.2 mm for the long level meter were achieved during a few days of test operation.     

Time domain analysis: an alternative way to interpret PGSE experiment.

We report a theoretical development which aims at interpreting Pulsed Gradient Spin Echo data for diffusing fluids in saturated porous media. It consists in analyzing the time dependence of PGSE amplitudes for each single gradient strength, and introduces a new diffusion coefficient D(q) as being continuously dependent on the length scale in the material. Both experimental measurements on water saturated bead packings and simulated experiments in 2D and 3D model systems are interpreted with this approach. D(q) is shown to give a new insight in the micro-macro transition. Its evolution vs. the length scale is non trivial, and can be sensitive to local slow kinetic effects.     

Impact of weak localization in the time domain

We find a renormalized "time-dependent diffusion coefficient," D(t), for pulsed excitation of a nominally diffusive sample by solving the Bethe-Salpeter equation with recurrent scattering. We observe a crossover in dynamics in the transformation from a quasi-1D to a slab geometry implemented by varying the ratio of the radius, R, to the length, L, of the cylindrical sample with reflecting side walls and open ends. Immediately after the peak of the transmitted pulse, D(t) falls linearly with a nonuniversal slope that approaches an asymptotic value for R/L>1. The value of D(t) extrapolated to t=0 depends only upon the dimensionless conductance g for R/L<1 and only upon kl(0) for R/L>1, where k is the wave vector and l(0) is the bare mean free path.     

Anomalous viscosity exponent in a discretized model for a chain diffusion in one dimension

We present a one-dimensional Monte Carlo simulation for the diffusion motion of a chain of N beads. We found that the scaling exponent for the viscosity can be smaller or greater than 3. This anomalous behavior cannot be attributed to the diffusivity scaling or the length fluctuations but is due to the chain dynamics details during diffusion in which the end beads play the key role. The viscosity exponent 3 and its expected relation with the diffusivity exponent are recovered in the asymptotic regime (N ↦∞). Received 24 September 2001 and Received in final form 28 January 2002     

基于片层光照明的新型单分子横向磁镊

磁镊是一种高精度的单分子技术,它用磁场对连有生物大分子的超顺磁球产生磁力,通过追踪磁球的位置来测量生物大分子的长度信息.磁镊包括横向磁镊和纵向磁镊.纵向磁镊空间精度高,但昂贵;横向磁镊简单便宜,但由于受其成像原理的限制,一般情况下只能连接较长的DNA等生物大分子,且其空间精度较差,进而限制了其应用范围.为了解决这个问题,本文改进了横向磁镊,用片层光照明的方法使光线主要被磁球散射,从而能够直接观察到吸附在样品槽侧壁上的磁球,这使得测量短连接的底物成为可能.对于实际应用的检测,首先测试了包含270 bp发卡结构的0.5μm双链DNA,用其中发卡结构的"折叠-去折叠"跳变过程证明了改进后的横向磁镊的确可以追踪短DNA等生物大分子.然后,进一步用16μm的λ-DNA检验了实验系统.最后,将新型横向磁镊与普通横向磁镊及纵向磁镊在小力和大力条件下拉伸不同长度DNA的噪声进行了比较,发现改进后的横向磁镊在空间精度上明显优于普通横向磁镊,与纵向磁镊相比也无明显差异.以上结果证明了改进后的横向磁镊的精度优势,并扩展了横向磁镊的应用范围.     

Series solutions of stagnation slip flow and heat transfer by the homotopy analysis method

An analytical approximation for the similarity solutions of the two- and three-dimensional stagnation slip flow and heat transfer is obtained by using the homotopy analysis method. This method is a series expansion method, but it is different from the perturbation technique, because it is independent of small physical parameters at all. Instead, it is based on a continuous mapping in topology so that it is applicable for not only weakly but also strongly nonlinear flow phenomena. Convergent [m,m] homotopy Padé approximants are obtained and compared with the numerical results and the asymptotic approximations. It is found that the homotopy Padé approximants agree well with the numerical results. The effects of the slip length ℓ and the thermal slip constant β on the heat transfer characteristics are investigated and discussed. Supported by the National Natural Science Foundation of China (Grant No. 10872129)     

The effect of rotational angle and experimental parameters on the diffraction patterns and micro-structural information obtained from q-space diffusion NMR: implication for diffusion in white matter fibers

Diffusion NMR may provide, under certain experimental conditions, micro-structural information about confined compartments totally non-invasively. The influence of the rotational angle, the pulse gradient length and the diffusion time on the diffusion diffraction patterns and q-space displacement distribution profiles was evaluated for ensembles of long cylinders having a diameter of 9 and 20 microm. It was found that the diffraction patterns are sensitive to the rotational angle (alpha) and are observed only when diffusion is measured nearly perpendicular to the long axis of the cylinders i.e., when alpha= 90 degrees +/- 5 degrees under our experimental conditions. More importantly, we also found that the structural information extracted from the displacement distribution profiles and from the diffraction patterns are very similar and in good agreement with the experimental values for cylinders of 20 microm or even 9 microm, when data is acquired with parameters that satisfy the short gradient pulse (SGP) approximation (i.e., delta -->0) and the long diffusion time limit. Since these experimental conditions are hardly met in in vitro diffusion MRI of excised organs, and cannot be met in clinical MRI scanners, we evaluated the effect of the pulse gradient duration and the diffusion time on the structural information extracted from q-space diffusion MR experiments. Indeed it was found that, as expected, accurate structural information, and diffraction patterns are observed when Delta is large enough so that the spins reach the cylinders' boundaries. In addition, it was found that large delta results in extraction of a compartment size, which is smaller than the real one. The relevance of these results to q-space MRI of neuronal tissues and fiber tracking is discussed.     

A cumulant analysis for non-Gaussian displacement distributions in Newtonian and non-Newtonian flows through porous media

We use displacement encoding pulsed field gradient (PFG) nuclear magnetic resonance to measure Fourier components S(q) of flow displacement distributions P(zeta) with mean displacement (zeta) for Newtonian and non-Newtonian flows through rocks and bead packs. Displacement distributions are non-Gaussian; hence, there are finite terms above second order in the cumulant expansion of ln(S(q)). We describe an algorithm for an optimal self-consistent cumulant analysis of data, which can be used to obtain the first three (central) moments of a non-Gaussian P(zeta), with error bars. The analysis is applied to Newtonian and non-Newtonian flows in rocks and beads. Flow with shear-thinning xanthan solution produces a 15.6+/-2.3% enhancement of the variance sigma(2) of displacement distributions when compared to flow experiments with water.