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...

Complex magnetic domain structures and the energy barriers between them are responsible for pseudo‐single‐domain phenomena in rock magnetism and contribute significantly to the magnetic remanence of paleomagnetic samples. This article introduces MERRILL, an open source software package for three‐dimensional micromagnetics optimized and designed for the calculation of such complex structures. MERRILL has a simple scripting user interface that requires little computational knowledge to use but provides research strength algorithms to model complex, inhomogeneous domain structures in magnetic materials. It uses a finite element/boundary element numerical method, optimally suited for calculating magnetization structures of local energy minima (LEM) in irregular grain geometries that are of interest to the rock and paleomagnetic community. MERRILL is able to simulate the magnetic characteristics of LEM states in both single grains, and small assemblies of interacting grains, including saddle‐point paths between nearby LEMs. Here the numerical model is briefly described, and an overview of the scripting language and available commands is provided. The open source nature of the code encourages future development of the model by the scientific community.

Magnetite (Fe3O4) is an important magnetic mineral to Earth scientists, as it carries the dominant magnetic signature in rocks, and the understanding of its magnetic recording fidelity provides a critical tool in the field of palaeomagnetism. However, reliable interpretation of the recording fidelity of Fe3O4 particles is greatly diminished over time by progressive oxidation to less magnetic iron oxides, such as maghemite (γ-Fe2O3), with consequent alteration of remanent magnetization potentially having important geological significance. Here we use the complementary techniques of environmental transmission electron microscopy and off-axis electron holography to induce and visualize the effects of oxidation on the magnetization of individual nanoscale Fe3O4 particles as they transform towards γ-Fe2O3. Magnetic induction maps demonstrate a change in both strength and direction of remanent magnetization within Fe3O4 particles in the size range dominant in rocks, confirming that oxidation can modify the original stored magnetic information.

Zircon crystals from the Jack Hills, Western Australia, are one of the few surviving mineralogical records of Earth’s first 500 million years and have been proposed to contain a paleomagnetic record of the Hadean geodynamo. A prerequisite for the preservation of Hadean magnetization is the presence of primary magnetic inclusions within pristine igneous zircon. To date no images of the magnetic recorders within ancient zircon have been presented. Here we use high-resolution transmission electron microscopy to demonstrate that all observed inclusions are secondary features formed via two distinct mechanisms. Magnetite is produced via a pipe-diffusion mechanism whereby iron diffuses into radiation-damaged zircon along the cores of dislocations and is precipitated inside nanopores and also during low-temperature recrystallization of radiation-damaged zircon in the presence of an aqueous fluid. Although these magnetites can be recognized as secondary using transmission electron microscopy, they otherwise occur in regions that are indistinguishable from pristine igneous zircon and carry remanent magnetization that postdates the crystallization age by at least several hundred million years. Without microscopic evidence ruling out secondary magnetite, the paleomagnetic case for a Hadean–Eoarchean geodynamo cannot yet been made.

Recordings of magnetic fields, thought to be crucial to our solar system’s rapid accretion, are potentially retained in unaltered nanometric low-Ni kamacite (~ metallic Fe) grains encased within dusty olivine crystals, found in the chondrules of unequilibrated chondrites. However, most of these kamacite grains are magnetically non-uniform, so their ability to retain four-billion-year-old magnetic recordings cannot be estimated by previous theories, which assume only uniform magnetization. Here, we demonstrate that non-uniformly magnetized nanometric kamacite grains are stable over solar system timescales and likely the primary carrier of remanence in dusty olivine. By performing in-situ temperature-dependent nanometric magnetic measurements using off-axis electron holography, we demonstrate the thermal stability of multi-vortex kamacite grains from the chondritic Bishunpur meteorite. Combined with numerical micromagnetic modeling, we determine the stability of the magnetization of these grains. Our study shows that dusty olivine kamacite grains are capable of retaining magnetic recordings from the accreting solar system.