Realistic Microstructure Simulator (RMS): Monte Carlo simulations of diffusion in three-dimensional cell segmentations of microscopy images

TL;DR

This paper introduces RMS, a Monte Carlo simulation pipeline that efficiently models diffusion in realistic 3D cell microstructures from microscopy images, enabling accurate validation of biophysical MRI models.

Contribution

The RMS pipeline is the first to simulate MR signals in complex, realistic cell geometries with permeable and irregular membranes, utilizing GPU acceleration and the corner reflector concept.

Findings

Validated diffusion in mouse brain white matter axons.

Demonstrated computational efficiency with GPU parallelization.

Provided theoretical framework for discretized diffusion and membrane interactions.

Abstract

Background: Monte Carlo simulations of diffusion are commonly used as a model validation tool as they are especially suitable for generating the diffusion MRI signal in complicated tissue microgeometries. New method: Here we describe the details of implementing Monte Carlo simulations in three-dimensional (3d) voxelized segmentations of cells in microscopy images. Using the concept of the corner reflector, we largely reduce the computational load of simulating diffusion within and exchange between multiple cells. Precision is further achieved by GPU-based parallel computations. Results: Our simulation of diffusion in white matter axons segmented from a mouse brain demonstrates its value in validating biophysical models. Furthermore, we provide the theoretical background for implementing a discretized diffusion process, and consider the finite-step effects of the particle-membrane reflection and permeation events, needed for efficient simulation of interactions with irregular boundaries, spatially variable diffusion coefficient, and exchange. Comparison with existing methods: To our knowledge, this is the first Monte Carlo pipeline for MR signal simulations in a substrate composed of numerous realistic cells, accounting for their permeable and irregularly-shaped membranes. Conclusions: The proposed RMS pipeline makes it possible to achieve fast and accurate simulations of diffusion in realistic tissue microgeometry, as well as the interplay with other MR contrasts. Presently, RMS focuses on simulations of diffusion, exchange, and T1 and T2 NMR relaxation in static tissues, with a possibility to straightforwardly account for susceptibility-induced T2* effects and flow.