An effective phase shift diffusion equation method for analysis of PFG normal and fractional diffusions

Revealing the effect of phase and time coupling on NMR relaxation rate by random walks in phase space

Phase-time coupling is a natural process in the phase random walks of a spin system; however, its effect on the nuclear magnetic resonance (NMR) relaxation is a challenge to the established theories such as the second-order quantum perturbation theory. This paper extends the recently developed phase diffusion method to treat the phase-time coupling effect, based on uncoupled phase diffusions, and coupled random walks. The instantaneous projection of the rotating random field is employed to get the accumulated phase of the NMR observable. In the static frame and the rotating frame, the phase diffusion coefficients are derived. The obtained theoretical results show that the phase-time coupling has a significant impact on the NMR relaxation rate: The angular frequency ω in the spectral density is modified to an apparent angular frequency ηω, where η is the phase-time coupling constant. The strongest coupling has η equaling 2, while η equaling 1 corresponds to the traditional results. As an example, the modified relaxation time expressions based on both monoexponential and nonmonoexponential functions can successfully fit the previously reported ^{13}C T_{1} NMR experimental data of polyisobutylene (PIB) in the blend of PIB and head-to-head poly(propylene). In contrast, the traditional relaxation rate expression based on the monoexponential time correlation function cannot fit such experimental data. With the phase-time coupling, the obtained characteristic time of the segmental motion is faster than that from conventional results.

Pulsed field gradient signal attenuation of restricted anomalous diffusions in plate, sphere, and cylinder with wall relaxation

The effect of boundary relaxation on pulsed field gradient (PFG) anomalous restricted diffusion is investigated in this paper. The PFG signal attenuation expressions of anomalous diffusion in plate, sphere, and cylinder are derived based on fractional calculus. In addition, approximate expressions for boundary relaxation induced short time signal attenuation under zero gradient field and boundary relaxation affected short time apparent diffusion coefficients are given in this paper. Unlike the exponential signal attenuation in normal diffusion, the PFG signal attenuation in anomalous diffusion with boundary relaxation is either a Mittag-Leffler-function-based attenuation or a stretched-exponential-function-based attenuation. The stretched exponential attenuations of all three structures clearly show the diffractive pattern. In contrast, only in the plate structure does the Mittag-Leffler-function-based attenuation display an obvious diffractive pattern. Additionally, anomalous diffusion with smaller time derivative order α has a weaker diffractive pattern and less signal attenuation. Moreover, the results demonstrate that boundary relaxation induced signal attenuation is significantly affected by the anomalous diffusion when no gradient field is applied. Meanwhile, the boundary relaxation significantly affects PFG signal attenuation of anomalous diffusion in the following ways: The boundary relaxation results in reduced radius from the minimum of the diffractive patterns, and it results in an increased apparent diffusion coefficient and decreased surfaces to volume ratio in varying the diffusion time experiment; the boundary relaxation also substantially affects the apparent diffusion coefficient of sphere structure in the variation of gradient experiment.

Age-related microstructural and physiological changes in normal brain measured by MRI γ-metrics derived from anomalous diffusion signal representation

Nowadays, increasing longevity associated with declining cerebral nervous system functions, suggests the need for continued development of new imaging contrast mechanisms to support the differential diagnosis of age-related decline. In our previous papers, we developed a new imaging contrast metrics derived from anomalous diffusion signal representation and obtained from diffusion-weighted (DW) data collected by varying diffusion gradient strengths. Recently, we highlighted that the new metrics, named γ-metrics, depended on the local inhomogeneity due to differences in magnetic susceptibility between tissues and diffusion compartments in young healthy subjects, thus providing information about myelin orientation and iron content within cerebral regions. The major structural modifications occurring in brain aging are myelinated fibers damage in nerve fibers and iron accumulation in gray matter nuclei. Therefore, we investigated the potential of γ-metrics in relation to other conventional diffusion metrics such as DTI, DKI and NODDI in detecting age-related structural changes in white matter (WM) and subcortical gray matter (scGM). DW-images were acquired in 32 healthy subjects, adults and elderly (age range 20-77 years) using 3.0T and 12 b-values up to 5000 s/mm². Association between diffusion metrics and subjects' age was assessed using linear regression. A decline in mean γ (Mγ) in the scGM and a complementary increase in radial γ (γ_⊥) in frontal WM, genu of corpus callosum and anterior corona radiata with advancing age were found. We suggested that the increase in γ⊥ might reflect declined myelin density, and Mγ decrease might mirror iron accumulation. An increase in D_// and a decrease in the orientation dispersion index (ODI) were associated with axonal loss in the pyramidal tracts, while their inverted trends within the thalamus were thought to be linked to reduced architectural complexity of nerve fibers. γ-metrics together with conventional diffusion-metrics can more comprehensively characterize the complex mechanisms underlining age-related changes than conventional diffusion techniques alone.

The fractal derivative wave equation: Application to clinical amplitude/velocity reconstruction imaging

This paper proposes a dissipative acoustic wave equation in which the fractal derivative is employed to represent dissipation. The proposed model is derived from the viscoelastic constitutive relationship via the fractal derivative. It is noted that the fractal derivative is a local operator and avoids the expensive computational costs of non-local fractional derivative, which is popular in recent decades to describe frequency-dependent dissipation in acoustic wave propagation in soft materials. The proposed model is tested to simulate the clinical amplitude/velocity reconstruction imaging of breast tumors, where the reflecting plate is imaged as an elevated line in correspondence to tumor. Numerical experiments show that the present model is capable of indicating the size, position and quantity of tumors. The comparative study confirms that the fractal derivative acoustic wave equation has an advantage over the fractional derivative model regarding computational costs.

Signal attenuation of PFG restricted anomalous diffusions in plate, sphere, and cylinder

Pulsed field gradient (PFG) NMR is a noninvasive tool to study anomalous diffusion, which exists widely in many systems such as in polymer or biological systems, in porous material, in single file structures and in fractal geometries. In a real system, the diffusion could be a restricted or a tortuous anomalous diffusion, rather than a free diffusion as the domains for fast and slow transport could coexist. Though there are signal attenuation expressions for free anomalous diffusion in literature, the signal attenuation formalisms for restricted anomalous diffusion is very limited, except for a restricted time-fractional diffusion within a plate reported recently. To better understand the PFG restricted fractional diffusion, in this paper, the PFG signal attenuation expressions were derived for three typical structures (plate, sphere, and cylinder) based on two models: fractal derivative model and fractional derivative model. These signal attenuation expressions include two parts, the time part T_n(t) and the space part X_n(r). Unlike normal diffusion, the time part T_n(t) in time-fractional diffusion can be either a Mittag-Leffler function from the fractional derivative model or a stretched exponential function from the fractal derivative model. However, provided the restricted normal diffusion and the restricted time-fractional diffusion are in an identical structure, they will have the same space part X_n(r) as both diffusions have the same space derivative parameter β equaling 2, therefore, they should have similar diffractive patterns. The restricted general fractional diffusion within a plate is also investigated, which indicates that at a long time limit, the diffusion type is insignificant to the diffractive pattern that depends only on the structure and the gradient pulses. The expressions describing the time-dependent behaviors of apparent diffusion coefficient D_f,app for restricted anomalous diffusion are also proposed in this paper. Both the short and long time-dependent behaviors of D_f,app are distinct from that of normal diffusion. The general expressions for PFG restricted curvilinear diffusion of tube model were derived in a conventional way and its result agree with that obtained from the fractional derivative model with α equaling 1/2. Additionally, continuous-time random walk simulation was performed to give good support to the theoretical results. These theoretical results reported here will be valuable for researchers in analyzing PFG anomalous diffusion.

Instantaneous signal attenuation method for analysis of PFG fractional diffusions

An instantaneous signal attenuation (ISA) method for analyzing pulsed field gradient (PFG) fractional diffusion (FD) has been developed, which is modified from the propagator approach developed in 2001 by Lin et al. for analyzing PFG normal diffusion. Both, the current ISA method and the propagator method have the same fundamental basis that the total signal attenuation (SA) is the accumulation of all the ISA, and the ISA is the average SA of the whole diffusion system at each moment. However, the manner of calculating ISA is different. Unlike the use of the instantaneous propagator in the propagator method, the current method directly calculates ISA as A(K(t'),t'+dt')/A(K(t'),t'), where A(K(t'),t'+dt') and A(K(t'),t') are the SA. This modification makes the current method applicable to PFG FD as the instantaneous propagator may not be obtainable in FD. The ISA method was applied to study PFG SA including the effect of finite gradient pulse widths (FGPW) for free FD, restricted FD and the FD affected by a non-homogeneous gradient field. The SA expressions were successfully obtained for all three types of free FDs while other current methods still have difficulty in obtaining all of them. The results from this method agree with reported results such as that obtained by the effective phase shift diffusion equation (EPSDE) method. The M-Wright phase distribution approximation was also used to derive an SA expression for time FD as a comparison, which agrees with ISA method. Additionally, the continuous-time random walk (CTRW) simulation was performed to simulate the SA of PFG FD, and the simulation results agree with the analytical results. Particularly, the CTRW simulation results give good support to the analytical results including FGPW effect for free FD and restricted time FD based on a fractional derivative model where there have been no corresponding theoretical reports to date. The theoretical SA expressions including FGPW obtained here such as [Formula: see text] may be applied to analyze PFG FD in polymer or biological systems with improved accuracy where SGP approximation cannot be satisfied. The method can perhaps provide new insight to FD MRI and hence benefit the development of diffusion biomarkers based on fractional derivative.

Analyzing the special PFG signal attenuation behavior of intermolecular MQC via the effective phase shift diffusion equation method

Inter-molecular multiple quantum coherence (iMQC) has important applications in NMR and MRI. However, the current theoretical methods still have some difficulties in analyzing the behavior of iMQC signal attenuation of pulsed field gradient diffusion experiments. In this paper, the iMQC diffusion experiments were analyzed by an effective phase shift diffusion equation (EPSDE) method, which is based on the idea that the accumulating phase shift (APS) can be viewed as the result of a diffusion process in virtual phase space (VPS) with effective diffusion coefficient K(2)(t) D (rad(2)/s) where K(t)=∫0 (t)γg(t')dt' is a wavenumber and D is the physical diffusion coefficient of the spin carrier in the real space. The term K(t(tot)) z1 needs to be added to the APS when K(t(tot)) is not zero. Most of the time, K(t(tot)) equals zero. However, in iMQC experiments, the condition K(t(tot)) equaling zero or being non-zero for each spin depends on the gradient pulse setting. The signal attenuations of these two types of iMQC, zero or non-zero K(t(tot)), were analyzed in detail for free and restricted diffusions, which shows that there are significant differences between these two types of iMQC. Particularly, if an apparent diffusion coefficient D(app) is used to analyze the signal attenuation, it equals nD for zero K(t(tot)) which agrees with current theoretical and experimental reports, while for non-zero K(t(tot)), it equals (2n - 1) D which agrees with experimental results from the literature; there are no similar theoretical results reported for comparison. The result that D(app) equals (2n - 1) D is important because the higher value of D(app) means that non-zero K(t(tot)) iMQC can potentially provide more contrast and measure slower diffusion rates than zero K(t(tot)) iMQC. The EPSDE method provides a new way to analyze iMQC diffusion experiments.