Dephasing and diffusion on the alveolar surface
Lukas Reinhold Buschle1, Felix Tobias Kurz1,2, Thomas Kampf3, Heinz-Peter Schlemmer1, and Christian Herbert Ziener1

1E010 Radiology, German Cancer Research Center (DKFZ), Heidelberg, Germany, 2Department of Neuroradiology, Heidelberg University Hospital, Heidelberg, Germany, 3Department of Experimental Physics 5, University of Würzburg, Würzburg, Germany

Synopsis

In lung tissue, the susceptibility difference between air-filled alveoli and surrounding tissue causes a strong dephasing of spin-bearing particles. The particles experience an averaged magnetic field due to diffusion effects. Thus, the dephasing process slows down. Both diffusion and susceptibility effects are described by the Bloch-Torrey equation that is solved for the local magnetization on the surface of alveoli. The analytical solution of the free induction decay is compared to in vivo measurements in human lung tissue.

Target audience

This work targets investigators with an interest in quantitative susceptibility weighted imaging and diffusion weighted imaging with special emphasis on theoretical aspects of MR imaging and signal formation.

Purpose

The susceptibility difference between air-filled alveoli and surrounding tissue causes microscopic magnetic field inhomogeneities that lead to a fast dephasing of the local magnetization. Yet, spin-bearing particles experience an averaged magnetic field, since they diffuse around the local field inhomogeneities. Thus, with increasing diffusion effects the dephasing process slows down. The free induction decay caused by diffusion and susceptibility effects is described by the Bloch-Torrey-equation. The first full analytic solution of this equation, solved on the surface of alveoli, is presented. It allows comparing theoretical predictions of the free induction decay with in vivo measurements in human lung tissue. Therefore, microscopic subvoxel information as e.g. the alveolar radius can be estimated. This will improve the quantitative diagnose of pulmonary emphysema or fibrosis.

Methods

A spherically shaped alveolus with radius $$$R$$$ in a static magnetic field $$$B_0$$$ is considered. Spin-bearing particles are localized on the surface of this alveolus [1]. The susceptibility difference $$$\Delta \chi$$$ between air-filled alveolus and surrounding tissue generates a local magnetic dipole field on the alveolar surface: $$$\omega(\theta) = \delta\omega [3\cos^2(\theta)-1]$$$, where $$$\theta$$$ is the angle to the static field $$$B_0$$$. The strength of the dipole field is given as $$$\delta\omega = \gamma B_0 \Delta\chi /3$$$, with the gyromagnetic ratio $$$\gamma$$$. Furthermore, the influence of diffusion on spin-bearing particles around the alveolus is characterized with diffusion coefficient $$$D$$$ and diffusion time $$$\tau = R^2/D$$$. The local transversal magnetization $$$m(\theta,t)= m_x(\theta,t)+i m_y(\theta,t)$$$ is described by the Bloch-Torrey-equation [2]: $$\frac{\partial}{\partial t} m(\theta,t) = [D\Delta – i\omega(\theta)] m(\theta,t).$$ The free induction decay $$$M(t)$$$ is obtained as a superposition of the local magnetization on the alveolar surface. Thus, solving the Bloch-Torrey-equation allows determining the free induction decay in dependence of microscopic tissue parameters.

Results

The Bloch-Torrey-equation is solved by an eigenfunction expansion: $$m(\theta,t) = m_0 e^{-2i\delta\omega t}\sum\limits_{k=0}^{\infty} a_k PS_{2k,0}(\sqrt{3i\tau\delta\omega},\cos(\theta)) e^{- \lambda_{2k,0}(\sqrt{3i\tau\delta\omega})\frac{t}{\tau}},$$ where $$$PS_{2k,0}(\sqrt{3i\tau\delta\omega},\cos(\theta))$$$ and $$$\lambda_{2k,0}(\sqrt{3i\tau\delta\omega})$$$ denote spheroidal functions and corresponding eigenvalues, respectively, according to the notation of Meixner and Schäfke [3]. The eigenvalue spectrum is shown in Fig. 1 in dependence of the parameter $$$\tau\delta\omega$$$. The local magnetization depends on the angle $$$\theta$$$ and on time $$$t$$$ and is shown in Fig. 2. The free induction decay can be obtained as a superposition of the local magnetization over all positions on the spherical shell.

For experimental verification of the theoretical results the free induction decay in human lung tissue was measured with a PRESS-sequence in a cubic $$$(15mm)^3$$$ voxel. In Fig. 3 the experimental free induction decay is compared to the theoretical free induction decay assuming typical microscopic parameters. Both curves agree very well, especially for small times.

Discussion

This work treats dephasing of the local magnetization under susceptibility and diffusion effects on a spherical shell. The local magnetization is governed by the Bloch-Torrey-equation, where the influence of susceptibility and diffusion effects is taken into account. As expected, the local magnetization exhibits the same symmetry as the local Larmor frequency. Due to this angular dependence of the local Larmor frequency, the eigenvalues are in general complex-valued. Therefore, the local magnetization as well as the free induction decay are dependent on both transversal components.

The free induction decay is monoexponential only for a large influence of diffusion (small values of $$$\tau\delta\omega$$$), since mainly the first summand with $$$k=0$$$ in the infinite sum contributes. With decreasing diffusion effects, more and more terms contribute and the free induction decay shows significant differences from a monoexponential decay.

The considered model is an appropiate description of human lung tissue in inspiration, since water molecues are mainly localized on the surface of the alveolus. Furthermore the averaged contribution of all neighboring alveoli vansihes in inspiration.

Since the free induction decay depends only on the parameter $$$\tau\delta\omega \propto \gamma B_0 \Delta\chi R^2/D$$$, the alveolar radius $$$R$$$ can be extracted from measured free induction decays [4]. Indeed, a calculated free induction decay and in vivo measurements in human lung tissue show an excellent agreement.

With the provided method it is possible to quantify the local mean alveolar size in an imaging voxel. This is highly relevant for the clinical application of lung MRI, since the acquisition of a free induction decay is a robust and fast imaging method.

Acknowledgements

No acknowledgement found.

References

[1] A. G. Cutillo. Application of magnetic resonance to the study of lung. Futura Publishing Company, Inc., Armonk, NY, 1996.

[2] H. C. Torrey. Bloch equations with diffusion terms. Phys Rev, 104:563–565, 1956.

[3] J. Meixner and F. W. Schäfke. Mathieusche Funktionen und Sphäroidfunktionen. Springer-Verlag, Berlin Heidelberg, 1954.

[4] R. Mulkern, S. Haker, H. Mamata, E. Lee, D. Mitsouras, K. Oshio, M. Balasubramanian, and H. Hatabu. Lung parenchymal signal intensity in MRI: A technical review with educational aspirations regarding reversible versus irreversible transverse relaxation effects in common pulse sequences. Concepts Magn Reson Part A, 43A:29–53, 2014.

Figures

Real and imaginary part of the eigenvalues $$$\lambda_{2k,0}(\sqrt{3i\tau\delta\omega})$$$. In the motional narrowing limit where the influence of diffusion dominates, the eigenvalues become purely real: $$$\lambda_{2k,0}(0) = 2k[2k+1]$$$.

Local magnetization $$$m_x(\theta, t)$$$ in dependence of the normalized time $$$t/\tau$$$ and angle $$$\theta$$$ for $$$\tau\delta\omega=10$$$. With increasing time, the x-component of the local magnetization decreases. However, it stays symmetric around $$$\theta = \pi/2$$$ for all times.

Comparison of in vivo measurements in human lung tissue with theoretical predictions. The measurements were performed with a PRESS-sequence at a Siemens Magnetom Aera 1.5T scanner in inspiration. The theoretical prediction and experimental measurement agree very well, especially for small times.



Proc. Intl. Soc. Mag. Reson. Med. 24 (2016)
1615