skip to main content

Title: Structure and Properties of Liquid Fe‐C Alloys at High Pressures by Experiments and First‐Principles Calculations
It is believed that the core formation processes sequestered a large majority of Earth’s carbon into its metallic core. Incorporation of carbon to liquid iron may significantly influence its properties under physicochemical conditions pertinent to the deep magma ocean and thus the chemical evolution of terrestrial planets and moons. Compared to available experimental data on the physical properties of crystalline iron alloys under pressure, there is a remarkable lack of data on the properties of liquid iron‐rich alloys, due to experimental challenges. Here we review experimental and computational results on the structure and properties of iron or iron‐nickel liquids alloyed with carbon upon compression. These laboratory data provide an important foundation on which the interpretation of ultrahigh pressure laboratory data and the verification of theoretical data will have to be based. The low‐pressure data can be used to validate results from theoretical calculations at the same conditions, and high‐pressure calculations can be used to estimate and predict liquid properties under core conditions. Availability of the liquid properties of Fe‐C liquids will provide essential data for stringent tests of carbon‐rich core composition models for the outer core.
Award ID(s):
Publication Date:
Journal Name:
Geophysical monograph
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. Light elements in Earth’s core play a key role in driving convection and influencing geodynamics, both of which are crucial to the geodynamo. However, the thermal transport properties of iron alloys at high-pressure and -temperature conditions remain uncertain. Here we investigate the transport properties of solid hexagonal close-packed and liquid Fe-Si alloys with 4.3 and 9.0 wt % Si at high pressure and temperature using laser-heated diamond anvil cell experiments and first-principles molecular dynamics and dynamical mean field theory calculations. In contrast to the case of Fe, Si impurity scattering gradually dominates the total scattering in Fe-Si alloys with increasing Si concentration, leading to temperature independence of the resistivity and less electron–electron contribution to the conductivity in Fe-9Si. Our results show a thermal conductivity of ∼100 to 110 W⋅m −1 ⋅K −1 for liquid Fe-9Si near the topmost outer core. If Earth’s core consists of a large amount of silicon (e.g., > 4.3 wt %) with such a high thermal conductivity, a subadiabatic heat flow across the core–mantle boundary is likely, leaving a 400- to 500-km-deep thermally stratified layer below the core–mantle boundary, and challenges proposed thermal convection in Fe-Si liquid outer core.
  2. In this study a high-pressure membrane reactor (MR) was employed to carry-out the methanol synthesis (MeS) reaction. Syngas was fed into the MR shell-side where a commercial MeS catalyst was used, while the tube-side is swept with a high boiling point liquid with good solubility towards methanol. A mesoporous alumina ceramic membrane was utilized, after its surface had been modified to be rendered more hydrophobic. The efficiency of the MR was investigated under a variety of experimental conditions (different pressures, temperatures, sweep liquid flow rates, and types of sweep liquids). The results reveal improved per single-pass carbon conversions when compared to the conventional packed-bed reactor. An ionic liquid (IL), 1-ethyl-3-methylimidazolium tetrafluoroborate ([EMIM][BF4]) was utilized in the MR as the sweep liquid. The experimental results are compared to those previously reported by our Group (Li and Tsotsis, J. Membrane Sci., 2019) while using a conventional petroleum-derived solvent as sweep liquid, tetraethylene glycol dimethyl ether (TGDE). Enhanced carbon conversion (over the petroleum-derived solvent) was obtained using the IL.
  3. Core formation in small planetary bodies likely involves percolation of immiscible liquids (e.g. S- and C- rich iron alloys) through pore spaces in a silicate medium. The manner in which this phenomenon occurs is not fully understood. Furthermore, it is unknown whether the metallic melts can physically segregate during percolation. To improve our understanding of core formation in small planetesimals, we performed analog experiments. We used an emulsion of oil and water to simulate an emulsion of S-rich and C-rich iron alloys, respectively. The experiments were performed in a Hele-Shaw cell, a thin “channel” made of two acrylic plates (51 cm x 15 cmx 1.3 cm) kept apart with a thin aluminum plate (0.27 mm). A U-shaped cut out of the aluminum plate formed the channel. We used a syringe pump to inject the emulsion into the channel through a hole in the top plate. We investigated the effect of injection rate and droplet size on the percolation behavior of the emulsion. We observed that droplet velocity was size dependent. The smallest droplet size detected was 0.0133 mm2 with a velocity of 0.67 mm/s. Medium size droplets ranged from 0.03mm2 – ~10 mm2 with average velocity of ~0.43 mm/s. Largermore »droplets moved faster: the largest droplet, with an area of 91.4 mm2, had a velocity of 7.95 mm/s. We suggest that (1) suspended droplets slow down when they begin to touch the Hele-Shaw plates (medium size droplets), and (2) droplets flow faster when they become large enough to deform with the flow. We also tested percolation through a channel filled with polydisperse acrylic particles of diameter < 50 µm. When injected into the granular matrix, the oil formed a wetting front while the water advanced in “pulses”. These pulses may represent the faster flow of larger water droplets. In conclusion, the size of the droplets affects their velocity and possibly their ability to migrate through pore networks. The results suggest that immiscible liquids could potentially segregate due to different percolation efficiencies of the non-wetting/wetting phases. Consequently, this would affect the distribution of the metallic components within differentiated planetesimals.« less
  4. Huber, Julie A. (Ed.)
    ABSTRACT Wind-driven upwelling followed by relaxation results in cycles of cold nutrient-rich water fueling intense phytoplankton blooms followed by nutrient depletion, bloom decline, and sinking of cells. Surviving cells at depth can then be vertically transported back to the surface with upwelled waters to seed another bloom. As a result of these cycles, phytoplankton communities in upwelling regions are transported through a wide range of light and nutrient conditions. Diatoms appear to be well suited for these cycles, but their responses to them remain understudied. To investigate the bases for diatoms’ ecological success in upwelling environments, we employed laboratory simulations of a complete upwelling cycle with a common diatom, Chaetoceros decipiens , and coccolithophore, Emiliania huxleyi . We show that while both organisms exhibited physiological and transcriptomic plasticity, the diatom displayed a distinct response enabling it to rapidly shift-up growth rates and nitrate assimilation when returned to light and available nutrients following dark nutrient-deplete conditions. As observed in natural diatom communities, C. decipiens highly expresses before upwelling, or frontloads, key transcriptional and nitrate assimilation genes, coordinating its rapid response to upwelling conditions. Low-iron simulations showed that C. decipiens is capable of maintaining this response when iron is limiting to growth,more »whereas E. huxleyi is not. Differential expression between iron treatments further revealed specific genes used by each organism under low iron availability. Overall, these results highlight the responses of two dominant phytoplankton groups to upwelling cycles, providing insight into the mechanisms fueling diatom blooms during upwelling events. IMPORTANCE Coastal upwelling regions are among the most biologically productive ecosystems. During upwelling events, nutrient-rich water is delivered from depth resulting in intense phytoplankton blooms typically dominated by diatoms. Along with nutrients, phytoplankton may also be transported from depth to seed these blooms then return to depth as upwelling subsides creating a cycle with varied conditions. To investigate diatoms’ success in upwelling regions, we compare the responses of a common diatom and coccolithophore throughout simulated upwelling cycles under iron-replete and iron-limiting conditions. The diatom exhibited a distinct rapid response to upwelling irrespective of iron status, whereas the coccolithophore’s response was either delayed or suppressed depending on iron availability. Concurrently, the diatom highly expresses, or frontloads, nitrate assimilation genes prior to upwelling, potentially enabling this rapid response. These results provide insight into the molecular mechanisms underlying diatom blooms and ecological success in upwelling regions.« less
  5. Volatiles including carbon and hydrogen are generally considered to be more soluble in silicate melts than in mantle rocks. How these melts contribute to the storage and distribution of key volatiles in Earth's interior today and during its early evolution, however, remains largely unknown. It is essential to improve our knowledge about volatiles-bearing silicate magmas over the entire mantle pressure regime. Here we investigate molten Mg FexSiO3 ( , 0.25) containing both carbon and hydrogen using first-principles molecular dynamics simulations. Our results show that the dissolution mechanism of the binary volatiles in melts varies considerably under different conditions of pressure and redox. When incorporated as CO2 and H2O components (corresponding to oxidizing conditions) almost all carbon and hydrogen form bonds with oxygen. Their speciation at low pressure consists of predominantly isolated molecular CO2, carbonates, and hydroxyls. More oxygenated species, including tetrahedrally coordinated carbons, hydrogen (O-H-O) bridges, various oxygen-joined complexes appear as melt is further compressed. When two volatiles are incorporated as hydrocarbons CH4 and C2H6 (corresponding to reducing conditions), hydroxyls are prevalent with notable presence of molecular hydrogen. Carbon-oxygen bonding is almost completely suppressed. Instead carbon is directly correlated with itself, hydrogen, and silicon. Both volatiles also show strong affinitymore »to iron. Reduced volatile speciation thus involves polymerized complexes comprising of carbon, hydrogen, silicon, and iron, which can be mostly represented by two forms: C1−4H1−5Si0−5O0−2 (iron-free) and C5−8H1−8Si0−6Fe5−8O0−2. The calculated partial molar volumes of binary volatiles in their oxidized and reduced incorporation decrease rapidly initially with pressure and then gradually at higher pressures, thereby systematically lowering silicate melt density. Our assessment of the calculated opposite effects of the volatile components and iron on melt density indicates that melt-crystal density crossovers are possible in the present-day mantle and also could have occurred in early magma ocean environments. Melts at upper mantle and transition zone conditions likely dissolve carbon and hydrogen in a wide variety of oxidized and non-oxygenated forms. Deep-seated partial melts and magma ocean remnants at lower mantle conditions may exsolve carbon as complex reduced species possibly to the core during core-mantle differentiation while retaining a majority of hydrogen as hydroxyls-associated species.« less