Interpretations of palaeomagnetic observations assume that naturally occurring magnetic particles can retain their primary magnetic recording over billions of years. The ability to retain a magnetic recording is inferred from laboratory measurements, where heating causes demagnetization on the order of seconds. The theoretical basis for this inference comes from previous models that assume only the existence of small, uniformly magnetized particles, whereas the carriers of palaeomagnetic signals in rocks are usually larger, non-uniformly magnetized particles, for which there is no empirically complete, thermally-activated model. This study has developed a thermally-activated numerical micromagnetic model that can quantitatively determine the energy barriers between stable states in nonuniform magnetic particles on geological time scales. We examine in detail the thermal stability characteristics of equidimensional cuboctahedral magnetite and find that contrary to previously published theories, such non-uniformly magnetized particles provide greater magnetic stability than their uniformly magnetized counterparts. Hence, non-uniformly magnetized grains, which are commonly the main remanence carrier in meteorites and rocks, can record and retain highfidelity magnetic recordings over billions of years.micromagnetics | paleomagnetism | geomagnetism S ince the 1900s magnetic recordings observed in rocks and meteorites have been studied to understand the evolution of the Earth and the Solar System. The validity of the findings from these studies depends on a theoretical understanding of rock-magnetic recordings provided by Néel (1, 2) and numerous experimental studies, for example, Strangway et. al. (3) and Evans and Wayman (4). The overwhelming evidence from these authors was that stable natural magnetic remanence (NRM) in rocks resides within ultrafine, uniformly magnetized particles, called single domain (SD) particles. Néel's theory (1, 2) for the behavior of thermally-activated SD particles describes a unique relationship between thermal and temporal stability and gave confidence that palaeomagnetic recordings that become unstable (unblocked) only at high temperatures, retain magnetic recordings from the time of their mineral crystallization, possibly as far back in time as four billion years ago.However, in the 1970s and 80s the widespread use of hysteresis parameters to characterize magnetic mineralogy (5) found that the majority of magnetic particles in rocks are not in uniform magnetic states, but are larger in size (80 -1000 nm) and contain complex magnetic states that are not described by either SD theory or the multidomain theory of micron-sized particles (2, 6). The term pseudo-single-domain (PSD) was coined for such particles and much effort was spent in determining the origin of their magnetic fidelity (Dunlop Psark, 1977, Moon and Merrill, 1985). Due to the complexity of the problem it has not been possible to determine the temporal stability of magnetisation in PSD grains on geological time scales fro...