Structural evolution of liquid silicates under conditions in Super-Earth interiors

Molten silicates at depth are crucial for planetary evolution, yet their local structure and physical properties under extreme conditions remain elusive due to experimental challenges. In this study, we utilize in situ X-ray diffraction (XRD) at the Matter in Extreme Conditions (MEC) end-station of the Linear Coherent Linac Source (LCLS) at SLAC National Accelerator Laboratory to investigate liquid silicates. Using an ultrabright X-ray source and a high-power optical laser, we probed the local atomic arrangement of shock-compressed liquid (Mg,Fe)SiO3 with varying Fe content, at pressures from 81(9) to 385(40) GPa. We compared these findings to ab initio molecular dynamics simulations under similar conditions. Results indicate continuous densification of the O-O and Mg-Si networks beyond Earth's interior pressure range, potentially altering melt properties at extreme conditions. This could have significant implications for early planetary evolution, leading to notable differences in differentiation processes between smaller rocky planets, such as Earth and Venus, and super-Earths, which are exoplanets with masses nearly three times that of Earth.

Ultrafast x-ray detection of low-spin iron in molten silicate under deep planetary interior conditions

The spin state of Fe can alter the key physical properties of silicate melts, affecting the early differentiation and the dynamic stability of the melts in the deep rocky planets. The low-spin state of Fe can increase the affinity of Fe for the melt over the solid phases and the electrical conductivity of melt at high pressures. However, the spin state of Fe has never been measured in dense silicate melts due to experimental challenges. We report detection of dominantly low-spin Fe in dynamically compressed olivine melt at 150 to 256 gigapascals and 3000 to 6000 kelvin using laser-driven shock wave compression combined with femtosecond x-ray diffraction and x-ray emission spectroscopy using an x-ray free electron laser. The observation of dominantly low-spin Fe supports gravitationally stable melt in the deep mantle and generation of a dynamo from the silicate melt portion of rocky planets.

The Critical Point and the Supercritical State of Alkali Feldspars: Implications for the Behavior of the Crust During Impacts

The position of the vapor-liquid dome and of the critical point determine the evolution of the outermost parts of the protolunar disk during cooling and condensation after the Giant Impact. The parts of the disk in supercritical or liquid state evolve as a single thermodynamic phase; when the thermal trajectory of the disk reaches the liquid-vapor dome, gas and melt separate leading to heterogeneous convection and phase separation due to friction. Different layers of the proto-Earth behaved differently during the Giant Impact depending on their constituent materials and initial thermodynamic conditions. Here we use first-principles molecular dynamics to determine the position of the critical point for NaAlSi3O8 and KAlSi3O8 feldspars, major minerals of the Earth and Moon crusts. The variations of the pressure calculated at various volumes along isotherms yield the position of the critical points: 0.5-0.8 g cm-3 and 5500-6000 K range for the Na-feldspar, 0.5-0.9 g cm-3 and 5000-5500 K range for the K-feldspar. The simulations suggest that the vaporization is incongruent, with a degassing of O2 starting at 4000 K and gas component made mostly of free Na and K cations, O2, SiO and SiO2 species for densities below 1.5 g cm-3. The Hugoniot equations of state imply that low-velocity impactors (<8.3 km s-1) would at most melt a cold feldspathic crust, whereas large impacts in molten crust would see temperatures raise up to 30000 K.

In situ X-ray diffraction of silicate liquids and glasses under dynamic and static compression to megabar pressures

Properties of liquid silicates under high-pressure and high-temperature conditions are critical for modeling the dynamics and solidification mechanisms of the magma ocean in the early Earth, as well as for constraining entrainment of melts in the mantle and in the present-day core-mantle boundary. Here we present in situ structural measurements by X-ray diffraction of selected amorphous silicates compressed statically in diamond anvil cells (up to 157 GPa at room temperature) or dynamically by laser-generated shock compression (up to 130 GPa and 6,000 K along the MgSiO3 glass Hugoniot). The X-ray diffraction patterns of silicate glasses and liquids reveal similar characteristics over a wide pressure and temperature range. Beyond the increase in Si coordination observed at 20 GPa, we find no evidence for major structural changes occurring in the silicate melts studied up to pressures and temperatures exceeding Earth's core mantle boundary conditions. This result is supported by molecular dynamics calculations. Our findings reinforce the widely used assumption that the silicate glasses studies are appropriate structural analogs for understanding the atomic arrangement of silicate liquids at these high pressures.

Universal relationship of compression shocks in two-dimensional Yukawa systems

Using molecular dynamical simulations, compressional shocks in two-dimensional (2D) dusty plasmas are quantitatively investigated under various conditions. A universal relationship between the thermal and the drift velocities after shocks is discovered in 2D Yukawa systems. Using the equation of state of 2D Yukawa liquids, and the obtained pressure from the Rankine-Hugoniot relation, an analytical relation between the thermal and the drift velocities is derived, which well agrees with the discovered universal relationship for various conditions.

Electrical conductivity and magnetic dynamos in magma oceans of Super-Earths

Super-Earths are extremely common among the numerous exoplanets that have been discovered. The high pressures and temperatures in their interiors are likely to lead to long-lived magma oceans. If their electrical conductivity is sufficiently high, the mantles of Super-Earth would generate their own magnetic fields. With ab initio simulations, we show that upon melting, the behavior of typical mantle silicates changes from semi-conducting to semi-metallic. The electrical conductivity increases and the optical properties are substantially modified. Melting could thus be detected with high-precision reflectivity measurements during the short time scales of shock experiments. We estimate the electrical conductivity of mantle silicates to be of the order of 100 Ω⁻¹ cm⁻¹, which implies that a magnetic dynamo process would develop in the magma oceans of Super-Earths if their convective velocities have typical values of 1 mm/s or higher. We predict exoplanets with rotation periods longer than 2 days to have multipolar magnetic fields.

With the discovery of large rocky exoplanets called Super-Earths, questions have arisen regarding the properties of their interiors and their ability to produce a magnetic field. Here, the authors show that under high pressure, molten silicates are semi-metallic and that magma oceans would host a dynamo process.

Continuous Sound Velocity Measurements along the Shock Hugoniot Curve of Quartz

We report continuous measurements of the sound velocity along the principal Hugoniot curve of α quartz between 0.25 and 1.45 TPa, as determined from lateral release waves intersecting the shock front as a function of time in decaying-shock experiments. The measured sound velocities are lower than predicted by prior models, based on the properties of stishovite at densities below ∼7  g/cm^{3}, but agree with density functional theory molecular dynamics calculations and an empirical wide-regime equation of state presented here. The Grüneisen parameter calculated from the sound velocity decreases from γ∼1.3 at 0.25 TPa to 0.66 at 1.45 TPa. In combination with evidence for increased (configurational) specific heat and decreased bulk modulus, the values of γ suggest a high thermal expansion coefficient at ∼0.25-0.65  TPa, where SiO_{2} is thought to be a bonded liquid. From our measurements, dissociation of the molecular bonds persists to ∼0.65-1.0  TPa, consistent with estimates by other methods. At higher densities, the sound velocity is close to predictions from previous models, and the Grüneisen parameter approaches the ideal gas value.

Evidence for a phase transition in silicate melt at extreme pressure and temperature conditions

Laser-driven shock compression experiments reveal the presence of a phase transition in MgSiO(3) over the pressure-temperature range 300-400 GPa and 10 000-16 000 K, with a positive Clapeyron slope and a volume change of ∼6.3 (±2.0) percent. The observations are most readily interpreted as an abrupt liquid-liquid transition in a silicate composition representative of terrestrial planetary mantles, implying potentially significant consequences for the thermal-chemical evolution of extrasolar planetary interiors. In addition, the present results extend the Hugoniot equation of state of MgSiO(3) single crystal and glass to 950 GPa.

Vibrational density of states and Lindemann melting law

We examine the Lindemann melting law at different pressures using the vibrational density of states (DOS), equilibrium melting curve, and Lindemann parameter delta(L) (fractional root-mean-squared displacement, rmsd, at equilibrium melting) calculated independently from molecular dynamics simulations of the Lennard-Jones system. The DOS is obtained using spectra analysis of atomic velocities and accounts for anharmonicity. The increase of delta(L) with pressure is non-negligible: delta(L) is about 0.116 and 0.145 at ambient and extreme pressures, respectively. If the component of rmsd normal to a reflecting plane as in the Debye-Waller-factor-type measurements using x rays is adopted for delta(L), these values are about 0.067 (+/-0.002) and 0.084 (+/-0.003), and are comparable with experimental and calculated values for face-centered-cubic elements. We find that the Lindemann relation holds accurately at ambient and high pressures. The non-negligible pressure dependence of delta(L) suggests that caution should be exerted in applying the Lindemann law to obtaining the high pressure melting curve anchored at ambient pressure.