Search

International audience; The complex transverse water proton magnetization subject to diffusion-encoding magnetic field gradient pulses in a heterogeneous medium such as brain tissue can be modeled by the Bloch-Torrey partial differential equation. The spatial integral of the solution of this equation in realistic geometry provides a gold-standard reference model for the diffusion MRI signal arising from different tissue micro-structures of interest. A closed form representation of this reference diffusion MRI signal has been derived twenty years ago, called Matrix Formalism that makes explicit the link between the Laplace eigenvalues and eigenfunctions of the biological cell and its diffusion MRI signal. In addition, once the Laplace eigendecomposition has been computed and saved, the diffusion MRI signal can be calculated for arbitrary diffusion-encoding sequences and b-values at negligible additional cost.Up to now, this representation, though mathematically elegant, has not been often used as a practical model of the diffusion MRI signal, due to the difficulties of calculating the Laplace eigendecomposition in complicated geometries. In this paper, we present a simulation framework that we have implemented inside the MATLAB-based diffusion MRI simulator SpinDoctor that efficiently computes the Matrix Formalism representation forrealistic neurons using the finite elements method. We show the Matrix Formalism representation requires around a few hundred eigenmodes to match the reference signal computed by solving the Bloch-Torrey equation when the cell geometry comes from realistic neurons. As expected, the number of required eigenmodes to match the reference signal increases with smaller diffusion time and higher b-values. We also converted the eigenvalues to alength scale and illustrated the link between the length scale and the oscillation frequency of the eigenmode in the cell geometry. We gave the transformation that links the Laplace eigenfunctions to the eigenfunctions of the Bloch-Torrey operator and computed the Bloch-Torrey eigenfunctions and eigenvalues. This work is another step in bringing advanced mathematical tools and numerical method development to the simulation and modeling ofdiffusion MRI.

The complex transverse water proton magnetization subject to diffusion-encoding magnetic field gradient pulses in a heterogeneous medium can be modeled by the multiple compartment Bloch-Torrey partial differential equation (BTPDE). A mathematical model for the time-dependent apparent diffusion coefficient (ADC), called the H-ADC model, was obtained recently using homogenization techniques on the BTPDE. Under the assumption of negligible water exchange between compartments, the H-ADC model produces the ADC of a diffusion medium from the solution of a diffusion equation (DE) subject to a time-dependent Neumann boundary condition. This paper describes a publicly available Matlab toolbox called SpinDoctor that can be used 1) to solve the BTPDE to obtain the dMRI signal (the toolbox provides a way of robustly fitting the dMRI signal to obtain the fitted ADC); 2) to solve the DE of the H-ADC model to obtain the ADC; 3) a short-time approximation formula for the ADC is also included in the toolbox for comparison with the simulated ADC. The PDEs are solved by P 1 finite elements combined with build-in Matlab routines for solving ordinary differential equations. The finite element mesh generation is performed using an external package called Tetgen that is included in the toolbox. SpinDoctor provides built-in options of including 1) spherical cells with a nucleus; 2) cylindrical cells with a myelin layer; 3) an extra-cellular space (ECS) enclosed either a) in a box or b) in a tight wrapping around the cells; 4) deformation of canonical cells by bending and twisting. 5) permeable membranes for the BT-PDE (the H-ADC assumes negligible permeability). Built-in diffusion-encoding pulse sequences include the Pulsed Gradient Spin Echo and the Oscilating Gradient Spin Echo., 49 pages, 18 figures