Geological materials and archeological artifacts containing magnetic particles record the direction and intensity of the past geomagnetic field as they cool. These thermoremanent magnetizations are our primary source of information on the behavior of the Earth's magnetic field. Obtaining reliable paleointensities and paleodirections from samples with a large variation in grain sizes is a challenge due to differences in magnetic behavior between grains that differ in size, shape, and chemistry (e.g., Dunlop & Özdemir, 1997;Tauxe & Yamazaki, 2015). In paleomagnetic measurement techniques that rely on bulk measurements, the contributions of individual grains are measured collectively, that is, the signals of many millions of grains result in a single magnetic moment for the entire sample. This possibly obscures information from grains that record the paleofield well by the signal of nonperfect recorders in the sample. Especially, the presence of large (>>1 μm), multidomain (MD) grains often prevents a reliable interpretation of a magnetic signal from a bulk sample. Samples consisting of predominantly single-domain (SD) grains or slightly larger (<1 μm) pseudo-single domain (PSD) grains with complex domain structures such as vortices or "flower states" generally produce more reliable paleomagnetic data (e.g., Nagy et al., 2017Nagy et al., , 2019.Over the past decades, a number of studies have focused on high-end magnetometry techniques to assess the magnetic state of magnetic recorders and micromagnetic processes in them on a (sub) micrometer scale (e.g.,