“…Modeling the complex gray matter microstructure is remarkably challenging, but, even in white matter, the multi-component simulations of phase evolution could only partially explain the variability of measured susceptibility across TEs ( Cronin et al, 2017 ). Even though QSM provides reproducible measurements of tissue susceptibility for matching TEs at the same field strength ( Lancione et al, 2019 , Rua et al, 2020 , Spincemaille et al, 2020 ) when scanning parameters such as coverage ( Elkady et al, 2016 , Karsa et al, 2018 ) and spatial resolution ( Karsa et al, 2018 , Zhou et al, 2017 ) are fixed and the subject’s head is accurately positioned along the same orientation ( Cronin and Bowtell, 2018 , Lancione et al, 2017 , Li et al, 2012 ), the choice of these experimental factors can affect differently the susceptibility values measured in distinct populations. In fact, the net magnitude and phase signal from a voxel with inhomogeneous composition sampled at a certain TE contain information on the subcompartments whose relaxation time T2* is not much shorter than TE.…”