Uniform-penalty inversion of multiexponential decay data. II. Data spacing, T(2) data, systemic data errors, and diagnostics

The basic method of UPEN (uniform penalty inversion of multiexponential decay data) is given in an earlier publication (Borgia et al., J. Magn. Reson. 132, 65-77 (1998)), which also discusses the effects of noise, constraints, and smoothing on the resolution or apparent resolution of features of a computed distribution of relaxation times. UPEN applies negative feedback to a regularization penalty, allowing stronger smoothing for a broad feature than for a sharp line. This avoids unnecessarily broadening the sharp line and/or breaking the wide peak or tail into several peaks that the relaxation data do not demand to be separate. The experimental and artificial data presented earlier were T(1) data, and all had fixed data spacings, uniform in log-time. However, for T(2) data, usually spaced uniformly in linear time, or for data spaced in any manner, we have found that the data spacing does not enter explicitly into the computation. The present work shows the extension of UPEN to T(2) data, including the averaging of data in windows and the use of the corresponding weighting factors in the computation. Measures are implemented to control portions of computed distributions extending beyond the data range. The input smoothing parameters in UPEN are normally fixed, rather than data dependent. A major problem arises, especially at high signal-to-noise ratios, when UPEN is applied to data sets with systematic errors due to instrumental nonidealities or adjustment problems. For instance, a relaxation curve for a wide line can be narrowed by an artificial downward bending of the relaxation curve. Diagnostic parameters are generated to help identify data problems, and the diagnostics are applied in several examples, with particular attention to the meaningful resolution of two closely spaced peaks in a distribution of relaxation times. Where feasible, processing with UPEN in nearly real time should help identify data problems while further instrument adjustments can still be made. The need for the nonnegative constraint is greatly reduced in UPEN, and preliminary processing without this constraint helps identify data sets for which application of the nonnegative constraint is too expensive in terms of error of fit for the data set to represent sums of decaying positive exponentials plus random noise.     

Examples of marginal resolution of NMR relaxation peaks using UPEN and diagnostics.

The multiexponential inversion program UPEN by the authors [J. Magn. Reson. 1998; 132: 65-77; Ibid. 2000;147:273-85] employs negative feedback to a regularization penalty to implement variable smoothing when both sharp and broad features appear on a single distribution of relaxation times. This allows a good fit to relaxation data that correspond to a sum of decaying exponentials plus random noise, but it usually does not give a good fit to data that are distorted by systematic errors from instrument problems, which can cause erroneous "resolution" or erroneous non-resolution of peaks. UPEN provides a series of diagnostic parameters to help identify such data problems that can lead to interpretation errors, and, in particular, to warn when a close call on the resolution or non-resolution of nearby peaks might be questionable. Examples are given from a series of T(2) data sets from desiccated bone samples, with examples where the presence of two peaks is required by good data, examples where the presence of two peaks is negated by good data, and examples where the resolution or non-resolution of peaks cannot be trusted because of instrumental distortions revealed by UPEN diagnostic parameters. It is suggested that processing relaxation data with UPEN in nearly real time could permit retaking data while a sample is still available if the diagnostic parameters show instrumental problems.     

基于差分进化算法的核磁共振T2谱多指数反演

提出了一种基于差分进化(DE)算法的核磁共振弛豫信号多指数反演新方法. 将核磁共振T₂谱反演问题转化为带非负约束的非线性优化问题,不需要预先给定横向弛豫时间T_{2分布,直接利用差分进化算法进行反演计算. 在测量信号低信噪比的情况下,计算机模拟和实验数据反演都表明了该方法在分析处理NMR弛豫信号中的有效性. 关键词： 核磁共振 多指数反演 差分进化 岩心分析}     

T1–T2 Correlation Spectra Obtained Using a Fast Two-Dimensional Laplace Inversion

Spin relaxation is a sensitive probe of molecular structure and dynamics. Correlation of relaxation time constants, such as T₁ and T₂, conceptually similar to the conventional multidimensional spectroscopy, have been difficult to determine primarily due to the absense of an efficient multidimensional Laplace inversion program. We demonstrate the use of a novel computer algorithm for fast two-dimensional inverse Laplace transformation to obtain T₁–T₂ correlation functions. The algorithm efficiently performs a least-squares fit on two-dimensional data with a nonnegativity constraint. We use a regularization method to find a balance between the residual fitting errors and the known noise amplitude, thus producing a result that is found to be stable in the presence of noise. This algorithm can be extended to include functional forms other than exponential kernels. We demonstrate the performance of the algorithm at different signal-to-noise ratios and with different T₁–T₂ spectral characteristics using several brine-saturated rock samples.     

A simulation based method to assess inversion algorithms for transverse relaxation data

NMR relaxometry is a very useful tool for understanding various chemical and physical phenomena in complex multiphase systems. A Carr-Purcell-Meiboom-Gill (CPMG) [P.T. Callaghan, Principles of Nuclear Magnetic Resonance Microscopy, Clarendon Press, Oxford, 1991] experiment is an easy and quick way to obtain transverse relaxation constant (T₂) in low field. Most of the samples usually have a distribution of T₂ values. Extraction of this distribution of T₂s from the noisy decay data is essentially an ill-posed inverse problem. Various inversion approaches have been used to solve this problem, to date. A major issue in using an inversion algorithm is determining how accurate the computed distribution is. A systematic analysis of an inversion algorithm, UPEN [G.C. Borgia, R.J.S. Brown, P. Fantazzini, Uniform-penalty inversion of multiexponential decay data, Journal of Magnetic Resonance 132 (1998) 65–77; G.C. Borgia, R.J.S. Brown, P. Fantazzini, Uniform-penalty inversion of multiexponential decay data II. Data spacing, T₂ data, systematic data errors, and diagnostics, Journal of Magnetic Resonance 147 (2000) 273–285] was performed by means of simulated CPMG data generation. Through our simulation technique and statistical analyses, the effects of various experimental parameters on the computed distribution were evaluated. We converged to the true distribution by matching up the inversion results from a series of true decay data and a noisy simulated data. In addition to simulation studies, the same approach was also applied on real experimental data to support the simulation results.     

Properties of (Bi1/9Na2/3)(Mn1/3Nb2/3)O3 analysed within dielectric permittivity,conductivity, electric modulus and derivative techniques approach

The (Bi_1/9Na_2/3)(Mn_1/3Nb_2/3)O₃ ceramics with perovskite structure were sintered. The XRD test proved that the samples are cubic (a?=?3.920?±?0.001?Å). Microstructure and atomic composition were determined with a SEM (JSM-5410) equipped with energy dispersion X-ray analyser (ISIS-300). The fluctuation in the chemical composition was found indicating on local disorder. Broadband dielectric spectroscopy in the range 10^?1–3?·?10⁷?Hz was applied within the range of 100–650?K. The real, ?′(f,?T), and imaginary, ?″(f,?T), parts of complex dielectric permitivity characteristics, both in the temperature and frequency domain, show relaxation processes partially covered by electric conductivity. At high temperatures the electric conductivity exhibits a thermally activated behaviour σ(f,?T)?∝?exp(?E _a/kT) but the variable range hopping (VRH) dependence σ?∝?exp[?(T ₀/T)^1/4] is manifested at low temperatures. The derivatives technique in the frequency (??log??/??log?ω) and temperature (??log??/?T) domain enabled various relaxation processes to be distinguished. The data converted to electric modulus representation, M*(f,?T)?=?1/?*, exhibited clearly resolved relaxation peaks. The relaxation times obtained from the peaks position show a slightly non-Arrhenius temperature behaviour with the activation energy varying in 0.4–0.6?eV range and characteristic time of the electric conductivity relaxation of the order of 10^?12?s. The relaxation times can be fitted at better accuracy with the VRH dependence where T ₀ are of the order of 10⁸?K. It is shown that the low frequency ac-conductivity converges to dc-conductivity and the relation σ(0)?～?ω_m?～?τ_m ^?1 typical for the disordered solids applies. The conduction current relaxation relationship behaves in accord with the VRH system: σ_dc?∝?(T/T ₀)^q (e ²/kT) ω_c, where ω_c?=?ν_ph exp[?(T ₀/T)^1/4] is valid for the locally disordered (Bi_1/9Na_2/3)(Mn_1/3Nb_2/3)O₃ compound.     

Relaxation times and line widths of isotopically-substituted nitroxides in aqueous solution at X-band

Optimization of nitroxides as probes for EPR imaging requires detailed understanding of spectral properties. Spin lattice relaxation times, spin packet line widths, nuclear hyperfine splitting, and overall lineshapes were characterized for six low molecular weight nitroxides in dilute deoxygenated aqueous solution at X-band. The nitroxides included 6-member, unsaturated 5-member, or saturated 5-member rings, most of which were isotopically labeled. The spectra are near the fast tumbling limit with T₁∼T₂ in the range of 0.50–1.1 μs at ambient temperature. Both spin–lattice relaxation T₁ and spin–spin relaxation T₂ are longer for ¹⁵N- than for ¹⁴N-nitroxides. The dominant contributions to T₁ are modulation of nitrogen hyperfine anisotropy and spin rotation. Dependence of T₁ on nitrogen nuclear spin state m_I was observed for both ¹⁴N and ¹⁵N. Unresolved hydrogen/deuterium hyperfine couplings dominate overall line widths. Lineshapes were simulated by including all nuclear hyperfine couplings and spin packet line widths that agreed with values obtained by electron spin echo. Line widths and relaxation times are predicted to be about the same at 250 MHz as at X-band.     

Linewidth-Resolved N HSQC, a Simple 3D Method to Measure N Relaxation Times from T1 and T2 Linewidths

A three-dimensional approach for measuring ¹⁵N relaxation times is described. Instead of selecting particular values for the relaxation period, in the proposed method the relaxation period is incremented periodically in order to create a 3D spectrum. This additional frequency domain of the transformed spectrum contains the relaxation time information in the T₁ and T₂ linewidths, and thus the longitudinal and transverse ¹⁵N relaxation times can be measured without determination of 2D cross peak volumes/intensities and subsequent curve fitting procedures.     

电子型掺杂高温超导体La2-xCexCuO4飞秒时间分辨动力学研究

利用飞秒时间分辨光抽运探测技术研究了电子型掺杂La_2-xCe_xCuO₄(LCCO)高温超导材料的准粒子超快动力学过程.得到低温(T<0.7T_c)、转变温度附近(0.7T_c≤T≤T_c)和高温(T>T_c)三个温区内的动力学行为.研究发 关键词： 电子型掺杂高温超导体 飞秒时间分辨 准粒子 声子瓶颈     

Saturation spectroscopy of coherent raman scattering in molecular gases

We present the theory of two-photon Raman saturation of a two-level Raman transition studied by an independent CARS process. The main goal here is to probe the saturated homogeneous Raman line shape. It is shown that there appears a saturation dip with a width determined by the relaxation timeT ₁. In the case of Doppler-broadened line the coherent Raman saturation spectroscopy may be used to determine both theT ₁ andT ₂ relaxation times.     

Influence of magnetic impurities on proton spin relaxation of water in clay

The¹H nuclear magnetic spin relaxation of water in slurry of kaolin clay was investigated in the presence of magnetite (black iron oxide, Fe₃O₄) at 0.2 T and room temperature. The water spectra at high magnetite contents showed two different resonances, presumably from surface-associated water and free interstitial water. The difference in observed resonance frequencies increased as much as 200 ppm with increasing magnetite content. The apparent nuclear magnetic resonance intensity decreased biexponentially as a function of magnetite added. The observedT ₂^* values at low magnetite contents were in accordance with the predicted values from the resonance intensities and the estimated magnetic susceptibilities. TheT ₁ relaxation was multiexponential in character, so a uniform penalty program was used for the analysis of distribution. At 0.2 T for¹H, kaolin slurry containing less than 5.5 ppm magnetite did not differ significantly from magnetite-free clay in the longitudinal relaxation rates of water. However, higher concentrations of magnetite produced features in theT ₁ distribution significantly different from those of magnetite-free clay. TheT ₂ could be approximated by monoexponential relaxation, probably because the fast-decaying components relaxed before they could be recorded. The apparent transverse relaxation ratesR ₂ increased linearly as a function of magnetite content. On the basis of the comparison of spin-echo and Carr-Purcell-Meiboom-Gill data, an empirical relation was derived to describe the signal loss due to diffusion. It can be expressed by a power function of magnetite amount, which is multiplied by the sum of volume-dependent and volume-independent terms.     

ABSORPTION INTENSITIES, PRESSURE-BROADENING AND LINE MIXING PARAMETERS OF SOME LINES OF NH3 IN THE ν4 BAND

The intensities and foreign gas broadening coefficients of 57 selected lines of the ν₄ band of NH₃ have been measured in the region of 1550 cm^-1 using a high resolution Brucker Fourier transform spectrometer. The line intensities were obtained by using the methods of absorbance at the line center and by fitting Voigt profiles to the measured shapes of the lines. The latter method also provides the collisional widths of the lines. In addition, collision cross relaxation coefficients of O₂ and air foreign gases were measured for 9 doublets of NH₃ in the ν₄ band. The J and K quantum numbers dependencies of pressure-broadening coefficients and line intensities are discussed. The observed air broadening and cross relaxation coefficients were found to be in reasonable quantitative agreement with the concentration-weighted average of the N₂ and O₂ broadening coefficients. The comparison of our present and previous results obtained for the NH₃–H₂, NH₃–air, NH₃–N₂ and NH₃–CO₂ collisions shows an increase of the pressure broadening and cross relaxation coefficients with quadrupole moment of the foreign gas. The analysis of the line intensities was based on the third-order theory of line strengths and yields effective transition moments, vibrational band strengths and correction parameters of the symmetric and antisymmetric partial bands of the fundamental ν₄ band.     

Nuclear quadrupole spin-lattice relaxation in Bi4Ge3O12 single crystals doped with atoms of d or f elements. Crystal field effects in compounds exhibiting anomalous magnetic properties

The nuclear quadrupole spin-lattice relaxation was studied in the range 4.2–300 K for single crystals of Bi₄Ge₃O₁₂ doped with minor amounts (the tenth fractions of mol%) of paramagnetic atoms of Cr, Nd, and Gd. Unusual spin dynamic features were recently found for these crystals at room temperature: a dramatic (up to 8-fold) increase in the effective nuclear quadrupole spin-spin relaxation time T ₂* occurred upon doping the pure Bi₄Ge₃O₁₂ sample. Unlike T ₂*, the effective spin-lattice relaxation time T ₁* at room temperature differs insignificantly for both doped and pure samples. But at lower temperatures, the samples exhibit considerably different behavior of the spin-lattice relaxation with temperature, which is caused by different contributions to the relaxation process of the dopant paramagnetic atoms. The distinctive maximum in the temperature dependence of the spin-lattice relaxation time for the Nd-doped crystal is shown to result from the crystal electric field effects.     

NMR and quantum chemistry study of mesoscopic structuring in 3-methylpyridine/water/NaBr mixture in the vicinity of the phase separation boundary

²³Na and ⁸¹Br NMR spin-lattice relaxation times and signal half widths (Δ_1/2) have been measured in 3-methylpyridine (3MP)/H₂O/NaBr mixture along T?=?294 and 301 K isotherms gradually increasing the mass fractions of salt (X) up to the phase separation boundary. The extreme narrowing condition and thus 1/T ₁?=?1/T ₂?=?πΔ _1/2 was found to be valid in all cases. Discontinuous changes in slope of 1/T _1,2?=?f(X) were detected, and then corresponding points on the phase diagram (X, T) were attributed to the borderline between the molecular-ionic solution and the area of enhanced mesoscopic structuring. A very strong relaxation effect was observed for ⁸¹Br nuclei reaching relaxation rates of 14,000 s^–1. ²³Na and ⁸¹Br NMR relaxation data together with calculations of quantum chemistry model of electrical field gradient tensor evidence the migration of 3-methylpyridinium at increasing X from the anions hydration shells towards cations. An interchange from migration to steady distribution regimes is observed for anions and vice versa for cations at the borderline of the structured phase.     

Changes in Porcine Muscle Water Characteristics during Growth—An in Vitro Low-Field NMR Relaxation Study

This study investigates the effects of developmental stage and muscle type on the mobility and distribution of water within skeletal muscles, using low-field ¹H-NMR transverse relaxation measurements in vitro on four different porcine muscles (M. longissimus dorsi, M. semitendinosus, M. biceps femoris, M. vastus intermedius) from a total of 48 pigs slaughtered at various weight classes between 25 kg and 150 kg. Principal component analysis (PCA) revealed effects of both slaughter weight and muscle type on the transverse relaxation decay. Independent of developmental stage and muscle type, distributed exponential analysis of the NMR T₂ relaxation data imparted the existence of three distinct water populations, T_2b, T₂₁, and T₂₂, with relaxation times of approximately 1–10, 45–120, and 200–500 ms, respectively. The most profound change during muscle growth was a shift toward faster relaxation in the intermediate time constant, T₂₁. It decreased by approx. 24% in all four muscle types during the period from 25 to 150 kg live weight. Determination of dry matter, fat, and protein content in the muscles showed that the changes in relaxation time of the intermediate time constant, T₂₁, during growth should be ascribed mainly to a change in protein content, as the protein content explained 77% of the variation in the T₂₁ time constant. Partial least squares (PLS) regression revealed validated correlations in the region of 0.58 to 0.77 between NMR transverse relaxation data and muscle development for all the four muscle types, which indicates that NMR relaxation measurements may be used in the prediction of muscle developmental stage.     

In vivo NMR T2 relaxation of experimental brain tumors in the cat: A multiparameter tissue characterization

Experimental gliomas (F98) were inoculated in cat brain for the systematic study of their in vivo T₂ relaxation time behavior. With a CPMG multi-echo imaging sequence, a train of 16 echoes was evaluated to obtain the transverse relaxation time and the magnetization M(0) at time t = 0. The magnetization decay curves were analyzed for biexponentiality. All tissues showed monoexponential T₂, only that of the ventricular fluid and part of the vital tumor tissue were biexponential. Based on these NMR relaxation parameters the tissues were characterized, their correct assignment being assured by comparison with histological slices. T₂ of normal grey and white matter was 74 ± 6 and 72 ± 6 msec, respectively. These two tissue types were distinguished through M(0) which for white matter was only 0.88 of the intensity of grey matter in full agreement with water content, determined from tissue specimens. At the time of maximal tumor growth and edema spread a tissue differentiation was possible in NMR relaxation parameter images. Separation of the three tissue groups of normal tissue, tumor and edema was based on T₂ with T_2(normal) < T_2(tumor) < T_2(edema). Using M(0) as a second parameter the differentiation was supported, in particular between white matter and tumor or edema. Animals were studied at 1–4 wk after tumor implantation to study tumor development. The magnetization M(0) of both tumor and peritumoral edema went through a maximum between the second and third week of tumor growth. T₂ of edema was maximal at the same time with 133 ± 4 msec, while the relaxation time of tumor continued to increase during the whole growth period, reaching values of 114 ± 12 msec at the fourth week. Thus, a complete characterization of pathological tissues with NMR relaxometry must include a detailed study of the developmental changes of these tissues to assure correct experimental conditions for the goal of optimal contrast between normal and pathological regions in the NMR images.     

La0.82Te0.18MnO3薄膜的输运特性和光诱导效应

用脉冲激光沉积法制备了非金属Te掺杂的钙钛矿锰氧化物La_0.82Te_0.18MnO₃单晶薄膜.该薄膜从83 K升温至373 K过程中发生金属-绝缘体相变，转变点温度为283 K.其电阻率在T<T_MI时符合电子-电子、电子-磁振子散射公式；在T>T_MI时为小极化子输运.薄膜在低温段连续激光(波长为532 nm，40 mW)作用下电阻率显著增大，电阻变化率在253 K达到最大值51.1%,该变化率远大于相同条件下的空穴掺杂材料；在高温段产生了较小的光电导，电阻变化率小于10%.这些现象主要与激光激励下自旋系统和小极化子的变化有关.La_0.82Te_0.18MnO₃薄膜在激光诱导下具有明显的与自旋相关的弛豫现象.激光开始作用时薄膜电阻率随时间的变化符合指数关系. 关键词： 0.82Te_0.18MnO₃薄膜')" href="#">La_0.82Te_0.18MnO₃薄膜 光诱导 输运特性 电子掺杂     

Air-broadening parameters in the ν3 band of NO2 using a multispectrum fitting technique

Air-broadened linewidths, pressure-induced shift coefficients and their temperature dependences were retrieved for over 1000 transitions in the ν₃ band of ¹⁴N¹⁶O₂ at 6 μm. In addition, precise line center positions and relative intensities were also determined. The results were obtained by fitting simultaneously 27 spectra recorded at high resolution (0.002-0.006 cm⁻¹) with two Fourier transform spectrometers and gas sample temperatures ranging from 206 to 298 K. It was necessary to modify the multispectrum fitting software to accommodate constraints on the retrieved parameters of closely spaced spin-split doublets in order to successfully determine their broadening and shift parameters. The variations of the widths, shifts and their temperature dependences with the quantum numbers were investigated. Subsets of the observed linewidths were reproduced to within 3% using an empirical smoothing function.     

Study of the diffusion of hydrogen in LaNi5+xH6 compounds by 1H NMR relaxation

The diffusion of hydrogen in LaNi_5+xH₆ (x=?0.2, 0.0, 0.2) has been investigated by NMR from 150 to 300 K. High-temperature data of the spin-spin relaxation time T₂ and the rotating frame spin lattice relaxation time T₁? are independent of stoichiometry but the data of the spin lattice relaxation time T₁ and low-temperature T₁? data are not, and they do not fit Torrey's relaxation model.     

Methods for correlating T1ρ and FID components in wideline H NMR studies of motionally heterogeneous polymer systems

We address the problem of correlating the observed FID and T_1ρ components in wideline ¹H relaxation measurements of motionally heterogeneous polymers, and show that different methods of data treatment can highlight different aspects of the correlations present. For a sample of polypropylene we find that the T_1ρ relaxation behaviour is driven by relaxation associated with the intermediate FID component, which strongly suggests a motionally inhomogeneous amorphous region in the sample.