skip to main content

Title: Thermal Conductivity of Silicate Liquid Determined by Machine Learning Potentials

Silicate liquids are important agents of thermal evolution, yet their thermal conductivity is largely unknown. Here, we determine the thermal conductivity of a silicate liquid by combining the Green‐Kubo method with a machine learning potential ofab initioquality over the entire pressure regime of the mantle. We find that the thermal conductivity of MgSiO3liquid is 1.1 W m−1 K−1at the 1 bar melting point, and 4.0 W m−1 K−1at core‐mantle boundary conditions. The thermal conductivity increases with compression, while remaining nearly constant on isochoric heating. The pressure dependence arises from the increasing bulk modulus on compression, and the weak temperature dependence arises from the saturation of the phonon mean free path due to structural disorder. The thermal conductivity of silicate liquids is less than that of ambient mantle, a contrast that may be important for understanding melt generation, and heat flux from the core.

more » « less
Award ID(s):
Author(s) / Creator(s):
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Seismic and mineralogical studies have suggested regions at Earth’s core-mantle boundary may be highly enriched in FeO, reported to exhibit metallic behavior at extreme pressure-temperature (PT) conditions. However, underlying electronic processes in FeO remain poorly understood. Here we explore the electronic structure ofB1-FeO at extreme conditions with large-scale theoretical modeling using state-of-the-art embedded dynamical mean field theory (eDMFT). Fine sampling of the phase diagram reveals that, instead of sharp metallization, compression of FeO at high temperatures induces a gradual orbitally selective insulator-metal transition. Specifically, atPTconditions of the lower mantle, FeO exists in an intermediate quantum critical state, characteristic of strongly correlated electronic matter. Transport in this regime, distinct from insulating or metallic behavior, is marked by incoherent diffusion of electrons in the conductingt2gorbital and a band gap in theegorbital, resulting in moderate electrical conductivity (~105S/m) with modestPTdependence as observed in experiments. Enrichment of solid FeO can thus provide a unifying explanation for independent observations of low seismic velocities and elevated electrical conductivities in heterogeneities at Earth’s mantle base.

    more » « less
  2. Abstract

    The interplay between crystal–melt and grain boundary interfaces in partially melted polycrystalline aggregates controls many physical properties of mantle rocks. To understand this process at the fundamental level requires improved knowledge about the interfacial structures and energetics. Here, we report the results of first-principles molecular dynamics simulations of two grain boundaries of (0l1)/[100] type for tilt angles of 30.4° and 49.6° and the corresponding solid–liquid interfaces in Mg2SiO4forsterite at the conditions of the upper mantle. Our analysis of the simulated position time series shows that structural distortions at the solid–liquid interfacial region are stronger than intergranular interfacial distortions. The calculated formation enthalpy of the solid–solid interfaces increases nearly linearly from 1.0 to 1.4 J/m2for the 30.4° tilt and from 0.8 to 1.0 J/m2for the 49.6° tilt with pressure from 0 to 16 GPa at 1500 K, being consistent with the experimental data. The solid–liquid interfacial enthalpy takes comparable values in the range 0.9 to 1.5 J/m2over similar pressure interval. The dihedral angle of the forsterite–melt system estimated using these interfacial enthalpies takes values in the range of 67° to 146°, showing a decreasing trend with pressure. The predicted dihedral angle is found to be generally larger than the measured data for silicate systems, probably caused by compositional differences between the simulation and the measurements.

    more » « less
  3. Abstract

    In Landau’s Fermi liquid picture, transport is governed by scattering between quasi-particles. The normal liquid3He conforms to this picture but only at very low temperature. Here, we show that the deviation from the standard behavior is concomitant with the fermion-fermion scattering time falling below the Planckian time,$$\frac{\hslash }{{k}_{{{{{{{{\rm{B}}}}}}}}}T}$$kBTand the thermal diffusivity of this quantum liquid is bounded by a minimum set by fundamental physical constants and observed in classical liquids. This points to collective excitations (a sound mode) as carriers of heat. We propose that this mode has a wavevector of 2kFand a mean free path equal to the de Broglie thermal length. This would provide an additional conducting channel with aT1/2temperature dependence, matching what is observed by experiments. The experimental data from 0.007 K to 3 K can be accounted for, with a margin of 10%, if thermal conductivity is the sum of two contributions: one by quasi-particles (varying as the inverse of temperature) and another by sound (following the square root of temperature).

    more » « less
  4. Abstract

    Ultrafast time‐domain thermoreflectance (TDTR) is utilized to extract the through‐plane thermal conductivity (ΛLSCO) of epitaxial La0.5Sr0.5CoO3−δ(LSCO) of varying thickness (<20 nm) on LaAlO3and SrTiO3substrates. These LSCO films possess ordered oxygen vacancies as the primary means of lattice mismatch accommodation with the substrate, which induces compressive/tensile strain and thus controls the orientation of the oxygen vacancy ordering (OVO). TDTR results demonstrate that the room‐temperatureΛLSCOof LSCO on both substrates (1.7 W m−1K−1) are nearly a factor of four lower than that of bulk single‐crystal LSCO (6.2 W m−1K−1). Remarkably, this approaches the lower limit of amorphous oxides (e.g., 1.3 W m−1K−1for glass), with no dependence on the OVO orientation. Through theoretical simulations, origins of the glass‐like thermal conductivity of LSCO are revealed as a combined effect resulting from oxygen vacancies (the dominant factor), Sr substitution, size effects, and the weak electron/phonon coupling within the LSCO film. The absence of OVO dependence in the measuredΛLSCOis rationalized by two main effects: (1) the nearly isotropic phononic thermal conductivity resulting from the imperfect OVO planes when δ is small; (2) the missing electronic contribution toΛLSCOalong the through‐plane direction for these ultrathin LSCO films on insulating substrates.

    more » « less
  5. 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. 
    more » « less