A. Pommier (2025, GRL).
The liquid outer core of Mercury is fundamental to the sustainability of its dynamo, which depends strongly on the composition and related physical properties of the iron alloy. Core composition is influenced by inner core crystallization and interactions with the adjacent reduced mantle. Chemical transport between core and mantle analogs was investigated using impedance spectroscopy at 5 GPa and up to 1973 K and electron microscopy. Results support an outer core that contains significant amounts of alloying agents (>~20 at.% Si, Ni, C, S). Bulk diffusion coefficients across the silicate-core interface are 1-5.10-8 m2/s at 1973 K. Electrical resistivity of Fe-Ni-Si-C/S liquids was measured at 5 GPa and up to 2123 K to estimate the thermal conductivity of the outer core, which ranges from 17 to 31 W/m.K. These low thermal conductivity values, compared to Fe-Si liquids, could increase the power available to the dynamo during core cooling. An Fe-Ni-Si-C/S outer core could enhance dynamo efficiency until carbides and Fe-Ni crystallize. Diamond is not expected in Mercury's core.
Support: the Carnegie Endowment.
Liquidus projections of phase diagrams for the Fe-Ni-C system at 6 GPa from the NIST Phase Equilibria Diagrams Online Database (2024) and references therein. The red area corresponds to a region of realistic core compositions ( 17 at.% (20 wt.%) Ni and 16 at.% (4 wt.%) C). Any liquid from this region would need to reach eutectic point E (i.e., particularly low temperature) in order to crystallize diamond. However, Teutectic is much lower than the temperature of the core, and therefore, core crystallization products will be carbides and Fe-Ni solids, but not diamond.
A. Pommier, M.J. Tauber, H. Pirotte, G. D. Cody, A. Steele, E. S. Bullock, B. Charlier, and B. O. Mysen (2023, GCA).
Elucidating sulfur speciation and bonding at the atomic scale is required to understand transport properties in S-bearing melts such as diffusivity, viscosity, and electrical conductivity at high temperature and pressure. These properties are fundamental to modeling the evolution of terrestrial planets and moons. Despite several investigations of sulfur speciation in glasses formed by temperature-quenching of their melts, questions remain regarding the structural role of S and its effect on transport properties under highly reducing conditions such as S-rich lava on Mercury.
We studied the role of sulfur in silicate glasses formed by quenching of melts at 1573-1673 K synthesized in evacuated silica tubes under highly reducing conditions (Delta IW=-5.8 and -6.4, with IW the iron-wüstite oxygen fugacity buffer). The compositions reproduce the silicate portion of enstatite chondrites, representative of the northern volcanic plains at the surface of Mercury. Major cations in the silicate glasses included Al, Mg, Ca, Na, and K; S was varied from 0 to ~5 wt.%. The samples were characterized with impedance spectroscopy performed at 2 and 4 GPa and from 475 up to 1738 K using a multi-anvil press and the 4-electrode technique, 29Si MAS NMR spectroscopy, Raman spectroscopy, and electron microscopy.
We observe that conductivity generally increases with the S content of the glass, though no systematic correlation is observed. Electrical activation energy Ea below the glass transition temperature ranges from 0.56 to 1.10 eV, in agreement with sodium being the main charge carrier in all samples. The glass transition is located at 650-750 K based on impedance measurements. Above Tg, Ea decreases (0.35-0.68 eV) and the conductivities of the samples are comparable (~5-8.10-3S/m) until 973 K. At T>1600 K, the melt fraction is 50-70% and melt conductivity varies from 0.8 to 2.4 S/m, with the melt containing 5 wt.% S being the most conductive. 29Si NMR results reveal that a portion of S bonds with silicon by substitution for oxygen, as previously observed in Na2S-SiO2 glasses (Asahi et al., 1998), affecting slightly the degree of polymerization of the glass. This result is in accord with a study on viscosity of S-bearing melts (Mouser et al., 2022). Raman spectra show that some S is isolated from the glass-network structure and combines with divalent cations, Ca2+and Mg2+, to form sulfide clusters in the glass, consistent with previous work (Namur et al., 2016). Our results also indicate that CaS and MgS are not as strong network modifiers as CaO and MgO. Our spectroscopic measurements do not reveal direct interactions between sulfur and sodium, which is consistent with only a moderate effect of S on conductivity. This study not only illuminates the role of sulfur in highly reduced complex silicate glasses and melts, but also exemplifies a multi-disciplinary approach that would be useful for the investigation of other geomaterials.
The multi-anvil cell assemblies used for electrical measurements are available to the scientific community via ASU.
Support: the Carnegie Endowment, NSF-CAREER award.
Back-scattered electron image of an S-bearing glass sample containing spherical sulfides retrieved from an electrical experiment at 2 GPa and quenched at 1738 K.
A. Pommier, K. Leinenweber, T. Tran (2019, EPSL); C. Davies, A. Pommier, S. Greenwood, A. Wilson (2024, EPSL).
The thermal state and chemistry of the metallic core of terrestrial planets and moons governs their evolution and in particular, the intrinsic magnetic field. The weak intrinsic magnetic field of Mercury is intimately tied to the structure and cooling history of its metallic core. Recent constraints about the planet’s internal structure are consistent with the presence of a FeS layer overlying a silicon-bearing core. The planet is characterised by unusual chemical characteristics and a weak magnetic field generated in a large metallic core, and its early evolution was also marked by the presence of a magnetic field, widespread volcanism, and global contraction.
In Pommier et al. (2019), we performed 4-electrode resistivity experiments on core analogues up to 10 GPa and over wide temperature ranges in order to investigate the insulating properties of core materials. Our results show that the FeS layer is liquid and insulating, and that the electrical resistivity of a miscible Fe-Si(-S) core is comparable to the one of an immiscible Fe-S, Fe-Si core. The difference in electrical resistivity between the FeS-rich layer and the underlying Fe-Si or Fe-Si-S core is at least 1 log unit at pressure and temperature conditions relevant to Mercury’s interior. If present, a thick (>40 km) FeS-rich shell is expected to maintain high temperatures across the core, and if the temperature in this layer departs from an adiabat, then this might affect the core cooling rate. The presence of a liquid and insulating shell is not inconsistent with a thermally stratified core in Mercury and might impact the generation and sustainability of a magnetic field. To test the latter, modeling the planet's evolution is required.
Therefore, in Davies et al. (2024), we developed a parameterized model of coupled core-mantle thermal and magnetic evolution considering a layered Fe-Si(-S) core structure with chemical and physical properties of the mantle and the core based on previous laboratory studies. We seek successful solutions that are consistent with observations of Mercury’s long-lived dynamo, total global contraction, present-day crustal thickness, and present-day interior structure. Successful solutions have a mantle reference viscosity > 10^21 Pa s, a silicon concentration in the core > 13 wt.%, a present inner core radius of ∼ 1000 − 1200 km and a thermally stable layer ∼ 500 − 800 km thick below the core-mantle boundary. Our results show that if present, a molten FeS layer atop the core has minimal effect on Mercury’s long-term thermal and magnetic evolution. Predictions from our models can be tested with upcoming Bepi-Colombo observations.
Support: NSF CAREER, NSF NSF-COMPRES IV EOID, and NSF-NERC grant.
The interior of present-day Mercury: the large metallic core (about 85 % of the planet's radius) has a fundamental effect on the cooling of the planet, including the magnetic field. Image credit: Navid Marvi/Carnegie Science.
A. Pommier, M. J. Tauber, C. Renggli, C. J. Davies, A. Wilson (2025, JGR-Planets).
Alkaline earth sulfides could be abundant in Mercury’s mantle, and knowledge of their melting and transport properties, such as electrical conductivity, is needed to study the inner structure. Defining the conductivity of the interior is a scientific objective of the ESA-JAXA BepiColombo mission. We conducted electrical experiments at pressures in the range 2 – 5 GPa and at temperatures up to ~2400 K on analogs of natural sulfides, i.e., Ca1-xMgxS with minor impurities (CaS, Ca0.8Mg0.2S, Ca0.6Mg0.4S, Ca0.4Mg0.6S, Ca0.2Mg0.8S, and MgS).
Electrical conductivity increases nonuniformly with temperature, with no systematic dependence on cation composition. At relatively low temperatures (near 1100 K, the conductivities span a wide range, whereas at higher temperatures the values converge within ~0.5 – 7 S/m at 1800 K and 5 GPa. The conductivity trends are complex, and likely reflect contributions from divalent cations, alkali metal and carbon impurities, which would similarly contribute to the conductivity of Mercury's crust and mantle. Melting is identified by a sharp increase in conductivity between ~1850 and 2100 K at 5 GPa. These transition temperatures are consistent with the presence of impurities.
Using electrical studies on relevant silicate minerals and petrological observations, we developed electrical conductivity-depth profiles of Mercury's mantle. Depending on the interconnectivity of the sulfide phase, the conductivity at the base of the mantle containing 8 vol.% sulfide ranges from ~0.2 to > 8 S/m. Our results can be tested with future observations from the ESA-JAXA BepiColombo mission.
Support: the Carnegie Endowment.
Thermal profile (left) and corresponding electrical conductivity-pressure profiles (middle and right) across the silicate portion of Mercury. Left: Temperature from Davies et al. (2024). Middle: Effect of sulfide content on the bulk conductivity of a silicate mantle made of 70 vol.% olivine (forsterite) and 30 vol.% orthopyroxene (enstatite). The conductivity of the sulfide phase is from this study (Mg0.8Ca0.2S), and the bulk conductivity is calculated using the geometric mean. Right: Effect of sulfide interconnectivity on the bulk conductivity of the silicate mantle. 8 vol.% sulfide is considered, and bulk conductivity is calculated using the modified Archie's law. m corresponds to the cementation exponent. A value of 1 for cementation exponent m corresponds to an interconnected sulfide phase.