jats:titleAbstract</jats:title>jats:pThree‐dimensional geometries of silicate‐hosted magnetic inclusions from the Harcus intrusion, South Australia have been determined using focused‐ion‐beam nanotomography (FIB‐nt). By developing an effective workflow, the geometries were reconstructed for magnetic particles in a plagioclase (162) and a pyroxene (282), respectively. For each inclusion, micromagnetic modeling using Micromagnetic Earth Related Rapid Interpreted Language Laboratory provided averaged hysteresis loops and backfield remanence curves of 20 equidistributed field directions together with average jats:italicM</jats:italic>jats:subjats:italics</jats:italic></jats:sub>, jats:italicM</jats:italic>jats:subjats:italicrs</jats:italic></jats:sub>, jats:italicH</jats:italic>jats:subjats:italicc</jats:italic></jats:sub>, and jats:italicH</jats:italic>jats:subjats:italiccr</jats:italic></jats:sub>. The micromagnetic structures within each silicate are single‐domain (SD), single‐vortex (SV), multivortex (MV) and multidomain states. They have been analyzed using domain‐state diagnostic plots, such as the Day plot and the Néel plot. SD particles can be subdivided into groups with dominant uniaxial anisotropy (jats:italicM</jats:italic>jats:subjats:italicrs</jats:italic></jats:sub>/jats:italicM</jats:italic>jats:subjats:italics</jats:italic></jats:sub> ∼0.5 and 10 < jats:italicH</jats:italic>jats:subjats:italicc</jats:italic></jats:sub> < 100 mT) and mixed uniaxial/multiaxial anisotropy (jats:italicM</jats:italic>jats:subjats:italicrs</jats:italic></jats:sub>/jats:italicM</jats:italic>jats:subjats:italics</jats:italic></jats:sub> ∼0.7 and 10 < jats:italicH</jats:italic>jats:subjats:italicc</jats:italic></jats:sub> < 30 mT). Most SV particles lie on a trend with 0 < jats:italicM</jats:italic>jats:subjats:italicrs</jats:italic></jats:sub>/jats:italicM</jats:italic>jats:subjats:italics</jats:italic></jats:sub> < 0.1 and 0 < jats:italicH</jats:italic>jats:subjats:italicc</jats:italic></jats:sub> < 10 mT, while others display a broad range of intermediate jats:italicM</jats:italic>jats:subjats:italicrs</jats:italic></jats:sub>/jats:italicM</jats:italic>jats:subjats:italics</jats:italic></jats:sub> and jats:italicH</jats:italic>jats:subjats:italicc</jats:italic></jats:sub> values. SV and MV states do not plot on systematic grain‐size trends. Instead, the multicomponent mixture of domain states within each silicate spans the entire range of natural variability seen in bulk samples. This questions the interpretation of bulk average hysteresis parameters in terms of grain size alone. FIB‐nt combined with large‐scale micromagnetic simulations provides a more complete characterization of silicate‐hosted carriers of stable magnetic remanence. This approach will improve the understanding of single‐crystal paleomagnetism and enable primary paleomagnetic data to be extracted from ancient rocks.</jats:p>