Exact analytical PGSE signal for diffusion confined to a cylindrical surface using a spectral Laplacian formalism
This provides a computationally efficient framework for diffusion MRI applications requiring many model evaluations, though it is incremental as it builds on existing spectral methods for a specific geometry.
The authors derived an exact analytical solution for diffusion MRI signals when diffusion is confined to a cylindrical surface, overcoming inaccuracies in existing models at high diffusion weightings, and validated it against Monte Carlo simulations with efficient numerical implementations.
Pulsed-gradient spin-echo (PGSE) MRI experiments probe molecular self-diffusion through spin phase accumulation under time-dependent magnetic field gradients. For diffusion confined to cylindrical surfaces, existing analytical signal models typically rely on the narrow-pulse limit, approximate treatments of finite gradient durations, or the Gaussian phase approximation, which become increasingly inaccurate at high diffusion weightings. Here, we derive an exact analytical solution of the Bloch-Torrey equation for diffusion confined to a cylindrical surface under finite PGSE gradients and obtain the corresponding diffusion MRI signal expression valid for arbitrary gradient durations and separations. The derivation is based on a spectral matrix formalism of the Laplace operator in the eigenbasis of the confining geometry. The signal is expressed as a product of non-commuting matrix exponentials, without approximations to the diffusion propagator or the spin phase distribution. We further introduce a reduced real spectral basis exploiting the symmetry of the cylindrical surface, substantially improving computational efficiency. Building on this exact formulation, we develop efficient numerical strategies for repeated signal evaluations, including a Strang splitting approximation of the matrix exponentials and an efficient computation of the spherical mean signal using Gauss-Legendre quadrature. The analytical signal is validated against Monte Carlo simulations over a wide range of cylinder radii and experimental parameters. The accelerated implementations are benchmarked against the exact formulation to quantify accuracy-runtime trade-offs. These results establish a computationally efficient framework for evaluating directional and orientationally averaged diffusion MRI signals in applications requiring large numbers of model evaluations.