Spectral properties of the Bloch-Torrey operator in three dimensions

Abstract

We consider the Bloch-Torrey operator, $-\Delta +igx$, that governs the time evolution of the transverse magnetization in diffusion magnetic resonance imaging (dMRI). Using the matrix formalism, we compute numerically the eigenvalues and eigenfunctions of this non-Hermitian operator for two bounded three-dimensional domains: a sphere and a capped cylinder. We study the dependence of its eigenvalues and eigenfunctions on the parameter $g$ and on the shape of the domain (its eventual symmetries and anisotropy). In particular, we show how an eigenfunction drastically changes its shape when the associated eigenvalue crosses a branch (or exceptional) point in the spectrum. Potential implications of this behavior for dMRI are discussed.

References in corpus (8)

• Visualization of Branch Points in PT-Symmetric Waveguides
• Exceptional Points in Atomic Spectra
• Bifurcation diagram and pattern formation in superconducting wires with electric currents
• Coupling of eigenvalues of complex matrices at diabolic and exceptional points
• Coulomb analogy for nonhermitian degeneracies near quantum phase transitions
• Unfolding of eigenvalue surfaces near a diabolic point due to a complex perturbation
• Relevance of complex branch points for partial wave analysis
• Diffusion MRI/NMR at high gradients: challenges and perspectives