Exact analytical PGSE signal for diffusion confined to a cylindrical surface using a spectral Laplacian formalism

Canales-Rodríguez, Erick Jorge, Tax, Chantal M. ORCID: https://orcid.org/0000-0002-7480-8817, Saez, Juan M Gorriz, Jones, Derek K. ORCID: https://orcid.org/0000-0003-4409-8049, Thiran, Jean-Philippe and Rafael-Patiño, Jonathan 2026. Exact analytical PGSE signal for diffusion confined to a cylindrical surface using a spectral Laplacian formalism. Physics in Medicine & Biology 10.1088/1361-6560/aeb2a0

PDF - Published Version Available under License Creative Commons Attribution. Download (1MB)

Abstract

Objective. To derive an exact analytical expression for the diffusion MRI signal arising from diffusion confined to a cylindrical surface under finite rectangular pulsed-gradient spin-echo (PGSE) gradients, and to develop computationally efficient strategies for repeated model evaluations. Approach. The Bloch–Torrey equation was solved using a spectral matrix formalism of the Laplace operator in the eigenbasis of the cylindrical surface. The resulting signal is expressed as a product of non-commuting matrix exponentials and is valid for arbitrary rectangular gradient durations and separations, without approximations to the diffusion propagator or spin phase distribution. A reduced real spectral basis was introduced to decrease the dimensionality of the problem. Accelerated implementations based on Strang splitting and Gauss–Legendre quadrature were developed for repeated signal evaluations and computation of the spherical mean. The analytical signal was validated against Monte Carlo diffusion simulations. Main results. The exact formulation showed excellent agreement with Monte Carlo simulations over a wide range of cylinder radii and acquisition parameters. The spectral expansion converged rapidly with mode order, and the reduced basis substantially decreased computational cost without loss of accuracy. The accelerated implementations provided large additional speed-ups while retaining controllable errors relative to the exact solution. Significance. The proposed framework extends exact finite-pulse solutions of the Bloch–Torrey equation to diffusion confined to cylindrical surfaces and provides practical implementations for directional and orientationally averaged diffusion MRI signals. It therefore offers an approximation-free reference for assessing faster signal approximations and for applications requiring large numbers of model evaluations.

Item Type:	Article
Date Type:	Published Online
Status:	In Press
Schools:	Schools > Physics and Astronomy Schools > Psychology Research Institutes & Centres > Cardiff University Brain Research Imaging Centre (CUBRIC)
Publisher:	IOP Publishing
ISSN:	0031-9155
Date of First Compliant Deposit:	13 October 2026
Last Modified:	13 Oct 2026 11:30
URI:	https://orca.cardiff.ac.uk/id/eprint/190208

Actions (repository staff only)

Edit Item