A simulation environment for diffusion weighted MR experiments in complex media

Simulations of diffusion in neural tissues have traditionally been limited to analytical solutions, to grid‐based solvers, or to small‐scale Monte Carlo simulations. None of these approaches has had the capability to simulate realistic complex neural tissues on the scale of even a single voxel of reasonable (i.e., clinical) size. An approach is described that combines a Monte Carlo Brownian dynamics simulator capable of simulating diffusion in arbitrarily complex polygonal geometries with a signal integrator flexible enough to handle a variety of pulse sequences. Taken together, this package provides a complete and general simulation environment for diffusion MRI experiments. The simulator is validated against analytical solutions for unbounded diffusion and diffusion between parallel plates. Further results are shown for aligned fibers, varying packing density and permeability, and for crossing straight fibers. Magn Reson Med, 2009. © 2009 Wiley‐Liss, Inc.

[1]  A. Anderson.  Measurement of fiber orientation distributions using high angular resolution diffusion imaging , 2005, Magnetic resonance in medicine.

[2]  Daniel C. Alexander,et al.  Camino: Open-Source Diffusion-MRI Reconstruction and Processing , 2006 .

[3]  C. Beaulieu,et al.  The basis of anisotropic water diffusion in the nervous system – a technical review , 2002, NMR in biomedicine.

[4]  William Schroeder,et al.  The Visualization Toolkit: An Object-Oriented Approach to 3-D Graphics , 1997 .

[5]  L. Frank.  Characterization of anisotropy in high angular resolution diffusion‐weighted MRI , 2002, Magnetic resonance in medicine.

[6]  Daniel C. Alexander,et al.  Finite pulse widths improve fibre orientation estimates in diffusion tensor MRI , 2006 .

[7]  Wolfgang Dreher,et al.  Diffusion in compartmental systems. I. A comparison of an analytical model with simulations , 2003, Magnetic resonance in medicine.

[8]  A. Szafer,et al.  An analytical model of restricted diffusion in bovine optic nerve , 1997, Magnetic resonance in medicine.

[9]  J. E. Tanner,et al.  Restricted Self‐Diffusion of Protons in Colloidal Systems by the Pulsed‐Gradient, Spin‐Echo Method , 1968 .

[10]  Erik De Schutter,et al.  Computational neuroscience : realistic modeling for experimentalists , 2000 .

[11]  J. Gore,et al.  Theoretical Model for Water Diffusion in Tissues , 1995, Magnetic resonance in medicine.

[12]  R W Cox,et al.  AFNI: software for analysis and visualization of functional magnetic resonance neuroimages. , 1996, Computers and biomedical research, an international journal.

[13]  Paul T. Callaghan,et al.  Pulsed-Gradient Spin-Echo NMR for Planar, Cylindrical, and Spherical Pores under Conditions of Wall Relaxation , 1995 .

[14]  Scott B. Baden,et al.  A large scale Monte Carlo simulator for cellular microphysiology , 2004, 18th International Parallel and Distributed Processing Symposium, 2004. Proceedings..

[15]  P. Basser,et al.  Estimation of the effective self-diffusion tensor from the NMR spin echo. , 1994, Journal of magnetic resonance. Series B.

[16]  L. Frank.  Anisotropy in high angular resolution diffusion‐weighted MRI , 2001, Magnetic resonance in medicine.

[17]  A Mohoric,et al.  Computer simulation of the spin-echo spatial distribution in the case of restricted self-diffusion. , 2001, Journal of magnetic resonance.

[18]  Joel R. Stiles,et al.  Miniature Endplate Current Rise Times <100 mu s from Improved Dual Recordings Can be Modeled with Passive Acetylcholine Diffusion from a Synaptic Vesicle , 1996 .

[19]  P. Basser,et al.  A model for diffusion in white matter in the brain. , 2005, Biophysical journal.

[20]  Kuchel,et al.  Analytical Solutions and Simulations for Spin-Echo Measurements of Diffusion of Spins in a Sphere with Surface and Bulk Relaxation , 1996, Journal of magnetic resonance. Series B.

[21]  M. Baslow,et al.  Evidence supporting a role for N-acetyl-l-aspartate as a molecular water pump in myelinated neurons in the central nervous system An analytical review , 2002, Neurochemistry International.

[22]  Scott N. Hwang,et al.  An image‐based finite difference model for simulating restricted diffusion , 2003, Magnetic resonance in medicine.

[23]  Erik De Schutter,et al.  Monte Carlo Methods for Simulating Realistic Synaptic Microphysiology Using MCell , 2000 .

[24]  V. Wedeen,et al.  Diffusion MRI of Complex Neural Architecture , 2003, Neuron.

[25]  T. Bartol,et al.  Miniature endplate current rise times less than 100 microseconds from improved dual recordings can be modeled with passive acetylcholine diffusion from a synaptic vesicle. , 1996, Proceedings of the National Academy of Sciences of the United States of America.