The origin of Earth's ancient magnetic field is an outstanding problem. It has recently been proposed that exsolution of MgO from the core may provide sufficient energy to drive an early geodynamo. Here we present new experiments on Mg partitioning between iron-rich liquids and silicate/oxide melts. Our results indicate that Mg partitioning depends strongly on the oxygen content in the iron-rich liquid, in contrast to previous findings that it depends only on temperature. Consequently, MgO exsolution during core cooling is drastically reduced and insufficient to drive an early geodynamo alone. Using the new experimental data, our thermal model predicts inner core nucleation at similar to 850 Ma and a nearly constant paleointensity.

The paleomagnetic record is central to our understanding of the history of the Earth. The orientation and intensity of magnetic minerals preserved in ancient rocks indicate the geodynamo has been alive since at least the Archean and possibly the Hadean. A paleomagnetic signature of the initial solidification of the inner core, arguably the singular most important event in core history, however, has remained elusive. In pursuit of this signature we investigate the assumption that the field is a geocentric axial dipole (GAD) over long time scales. We study a suite of numerical dynamo simulations from a paleomagnetic perspective to explore how long the field should be time-averaged to obtain stable paleomagnetic pole directions and intensities. We find that running averages over 20 - 40 kyr are needed to obtain stable paleomagnetic poles with alpha(95) <10 degrees, and over 40 - 120 kyr for alpha(95) <5 degrees, depending on the variability of the field. We find that models with higher heat flux and more frequent polarity reversals require longer time averages, and that obtaining stable intensities requires longer time averaging than obtaining stable directions. Running averages of local field intensity and inclination produce reliable estimates of the underlying dipole moment when reversal frequency is low. However, when heat flux and reversal frequency are increased we find that local observations tend to underestimate virtual dipole moment (VDM) by up to 50% and overestimate virtual axial dipolemoment (VADM) by up to 150%. A latitudinal dependence is found where VDM underestimates the true dipole moment more at low latitudes, while VADM overestimates the true axial dipole moment more at high latitudes. The cause for these observed intensity biases appears to be a contamination of the time averaged field by non-GAD terms, which grows with reversal frequency. We derive a scaling law connecting reversal frequency and site paleolatitude to paleointensity bias (ratio of observed to the true value). Finally we apply this adjustment to the PINT paleointensity record. These biases produce little change to the overall trend of a relatively flat but scattered intensity over the last 3.5 Ga. A more careful intensity adjustment applied during periods when the reversal frequency is known could reveal previously obscured features in the paleointensity record.

In a metal, as in Earth's core, the thermal and electrical conductivities are assumed to be correlated. In a planetary dynamo this implies a contradiction: that both electrical conductivity, which makes it easier to induce current and magnetic field, and conductive heat transport, which hinders thermal convection, should increase simultaneously. Here we show that this contradiction implies that the magnetic induction rate peaks at a particular value of electrical and thermal conductivity and derive the low- and high-conductivity limits for thermal dynamo action. A dynamo regime diagram is derived as a function of electrical conductivity and temperature for Earth's core that identifies four distinct dynamo regimes: no dynamo, thermal dynamo, compositional dynamo, and thermocompositional dynamo. Estimates for the temperature-dependent electrical conductivity of the core imply that the geodynamo may have come close to its high-conductivity "no dynamo" limit prior to inner core nucleation, consistent with recent paleomagnetic observations.

Proxima Centauri b provides an unprecedented opportunity to understand the evolution and nature of terrestrial planets orbiting M dwarfs. Although Proxima Cen b orbits within its star's habitable zone, multiple plausible evolutionary paths could have generated different environments that may or may not be habitable. Here, we use 1-D coupled climate-photochemical models to generate self-consistent atmospheres for several evolutionary scenarios, including high-O-2, high-CO2, and more Earth-like atmospheres, with both oxic and anoxic compositions. We show that these modeled environments can be habitable or uninhabitable at Proxima Cen b's position in the habitable zone. We use radiative transfer models to generate synthetic spectra and thermal phase curves for these simulated environments, and use instrument models to explore our ability to discriminate between possible planetary states. These results are applicable not only to Proxima Cen b but to other terrestrial planets orbiting M dwarfs. Thermal phase curves may provide the first constraint on the existence of an atmosphere. We find that James Webb Space Telescope (JWST) observations longward of 10 mu m could characterize atmospheric heat transport and molecular composition. Detection of ocean glint is unlikely with JWST but may be within the reach of larger-aperture telescopes. Direct imaging spectra may detect O-4 absorption, which is diagnostic of massive water loss and O-2 retention, rather than a photosynthetic biosphere. Similarly, strong CO2 and CO bands at wavelengths shortward of 2.5 mu m would indicate a CO2-dominated atmosphere. If the planet is habitable and volatile-rich, direct imaging will be the best means of detecting habitability. Earth-like planets with microbial biospheres may be identified by the presence of CH4-which has a longer atmospheric lifetime under Proxima Centauri's incident UV-and either photosynthetically produced O-2 or a hydrocarbon haze layer.

We describe a software package called VPLanet that simulates fundamental aspects of planetary system evolution over Gyr timescales, with a focus on investigating habitable worlds. In this initial release, eleven physics modules are included that model internal, atmospheric, rotational, orbital, stellar, and galactic processes. Many of these modules can be coupled to simultaneously simulate the evolution of terrestrial planets, gaseous planets, and stars. The code is validated by reproducing a selection of observations and past results. VPLanet is written in C and designed so that the user can choose the physics modules to apply to an individual object at runtime without recompiling, i.e., a single executable can simulate the diverse phenomena that are relevant to a wide range of planetary and stellar systems. This feature is enabled by matrices and vectors of function pointers that are dynamically allocated and populated based on user input. The speed and modularity of VPLanet enables large parameter sweeps and the versatility to add/remove physical phenomena to assess their importance. VPLanet is publicly available from a repository that contains extensive documentation, numerous examples, Python scripts for plotting and data management, and infrastructure for community input and future development.

The presence of a planetary magnetic field is an important ingredient for habitability. The coexistence of a solid and a liquid core can facilitate the maintenance of a compositionally driven dynamo; however, the likelihood of such a configuration in super-Earths is unknown. Recently, shock experiments and ab initio calculations have constrained the stability, equations of state, and melting properties of ultrahigh pressure core and mantle phases. Here, we investigate the internal structure of super-Earth exoplanets with a range of total masses and core mass fractions ranging from that of Mars (0.2) to Mercury (0.68), including an Earth-like bulk composition. We examine the effect of the initial core-mantle boundary temperature (T-CMB) on their internal structure and identify regimes with coexisting solid and liquid cores, and deep mantle melting. We find that the range of T-CMB for which an inner core is growing increases with the total planet mass and even more with the core mass fraction. Therefore, our results suggest that super-Earths should have a crystallizing core over a large temperature range. We also find that the presence of a growing inner core is likely to be accompanied by a partially liquid lower mantle, which will likely influence planetary thermal evolution. We estimate the initial CMB temperature after super-Earth accretion by assuming an accretional heat retention efficiency similar to Earth. We find that massive super-Earths are expected to have an initial internal temperature consistent with a partially liquid core, allowing for the possibility of thermal and compositional dynamo action.

MgO exsolution has been proposed to drive an early geodynamo. Experimental studies, however, have drawn different conclusions regarding the applicability of MgO exsolution. While many studies suggest that significant Mg can dissolve into the Earth's core, the amount of MgO exsolved out of the core, which hinges on the temperature dependence of MgO solubility, remains unclear. Here we present new high-temperature experiments to better constrain the temperature and compositional dependence of Mg partitioning between Fe alloys and silicate liquids. Our experiments show that Mg partitioning is weakly dependent on temperature, while confirming its strong dependence on oxygen content in Fe alloys. This implies that MgO exolution is limited as the core cools but can help drive an early geodyanamo if the core heat loss is slightly subadiabatic. If an exosolution-driven geodynamo did occur, it was likely over a limited time span that depends on the core thermal history and conductivity.