AbstractPurposeTo create a realisticin-silicohead phantom for the second Quantitative Susceptibility Mapping (QSM) Reconstruction Challenge and for future evaluations of processing algorithms for QSM.MethodsWe created a digital whole-head tissue property phantom by segmenting and post-processing high-resolution (0.64mm isotropic), multi-parametric MRI data acquired at 7T from a healthy volunteer. We simulated the steady-state magnetization at 7T using a Bloch simulator and mimicked a Cartesian sampling scheme through Fourier-based processing. Computer code for generating the phantom and performing the MR simulation was designed to facilitate flexible modifications of the phantom in the future, such as the inclusion of pathologies, as well as the simulation of a wide range of acquisition protocols. Specifically, the following parameters and effects were implemented: repetition time and echo time, voxel size, background fields, and RF phase biases. Diffusion weighted imaging phantom data is provided allowing future investigations of tissue microstructure effects in phase and QSM algorithms.ResultsThe brain-part of the phantom featured realistic morphology with spatial variations in relaxation and susceptibility values similar to the in vivo setting. We demonstrated some of the phantom’s properties, including the possibility of generating phase data with non-linear evolution over echo time due to partial volume effects or complex distributions of frequency shifts within the voxel.ConclusionThe presented phantom and computer programs are publicly available and may serve as a ground truth in future assessments of the faithfulness of quantitative susceptibility reconstruction algorithms.