Robust Monte-Carlo Simulations in Diffusion-MRI: Effect of the substrate complexity and parameter choice on the reproducibility of results

TL;DR

This paper investigates how substrate complexity and simulation parameters affect the reproducibility of Monte-Carlo diffusion MRI simulations, proposing a framework for more realistic and consistent modeling of microstructural environments.

Contribution

It identifies key pitfalls in Monte-Carlo diffusion MRI simulations and introduces a new framework to generate complex, realistic substrates that improve simulation accuracy and reproducibility.

Findings

Simulation parameters significantly influence signal accuracy.

Simplified substrates can bias microstructure feature recovery.

The proposed framework creates complex, high-density crossing substrates.

Abstract

Monte-Carlo Diffusion Simulations (MCDS) have been used extensively as a ground truth tool for the validation of microstructure models for Diffusion-Weighted MRI. However, methodological pitfalls in the design of the biomimicking geometrical configurations and the simulation parameters can lead to approximation biases. Such pitfalls affect the reliability of the estimated signal, as well as its validity and reproducibility as ground truth data. In this work, we first present a set of experiments in order to study three critical pitfalls encountered in the design of MCDS in the literature, namely, the number of simulated particles and time steps, simplifications in the intra-axonal substrate representation, and the impact of the substrate's size on the signal stemming from the extra-axonal space. The results obtained show important changes in the simulated signals and the recovered microstructure features when changes in those parameters are introduced. Thereupon, driven by our findings from the first studies, we outline a general framework able to generate complex substrates. We show the framework's capability to overcome the aforementioned simplifications by generating a complex crossing substrate, which preserves the volume in the crossing area and achieves a high packing density. The results presented in this work,along with the simulator developed, pave the way towards more realistic and reproducible Monte-Carlo simulations for Diffusion-Weighted MRI.