Ab initio lattice thermal conductivity of MgSiO3 perovskite as found in Earth's lower mantle

Recent progress in the study on phonon heat transport property of Earth's lower mantle minerals

The lattice thermal conductivities (κlat) of Earth's lower mantle (LM) minerals is a crucial parameter in the study of deep Earth dynamics and its determination is also one of the grand challenges in condensed matter physics. Here, we review recent progress on theoretical and experimental studies for theκlatunder high pressure (P) and high temperature (T) condition up to 150 GPa and 4000 K. After the critical parameters necessary to obtain converged values of theκlatare summarized, the theoreticalκlatof the LM minerals, determined through various computational methodologies, is compiled along with experimental findings. Although significant scattering is found in the experimental results at LMP,T, the quantum anharmonic lattice dynamics theory combined with the phonon Boltzmann transport theory demonstrates a clear relationship in theκlatof the end-member LM phases, MgO, MgSiO3bridgmanite (Brg) and post-perovskite (PPv),κlatMgO>>κlatPPv>κlatBrg, and a discontinuous change in theκlatby ∼20%-50% expected across the Brg-PPv transition. Knowledge on the additional but geophysically important factors, such as the effects of iron solid solution, isotopic mass difference, and higher order crystal anharmonicity are also summarized in detail. Current problems and future perspectives are finally mentioned.

Unexpected thermal transport properties of MgSiO3monolayer at extreme conditions

The thermal transport properties of mantle minerals are of paramount importance to understand the thermal evolution processes of the Earth. Here, we perform extensively structural searches of two-dimensional MgSiO3monolayer by CALYPSO method and first-principles calculations. A stable MgSiO3monolayer withPmm2 symmetry is uncovered, which possesses a wide indirect band gap of 4.39 eV. The calculations indicate the lattice thermal conductivities of MgSiO3monolayer are 49.86 W (mK)-1and 9.09 W (mK)-1inxandydirections at room temperature. Our findings suggest that MgSiO3monolayer is an excellent low-dimensional thermoelectric material with highZTvalue of 4.58 from n-type doping in theydirection at 2000 K. The unexpected anisotropic thermal transport of MgSiO3monolayer is due to the puckered crystal structure and the asymmetric phonon dispersion as well as the distinct electron states around the Fermi level. These results offer a detailed description of structural and thermal transport properties of MgSiO3monolayer at extreme conditions.

Non-monotonic pressure dependence of the thermal conductivity of boron arsenide

Recent experiments demonstrate that boron arsenide (BAs) is a showcase material to study the role of higher-order four-phonon interactions in affecting heat conduction in semiconductors. Here we use first-principles calculations to identify a phenomenon in BAs and a related material - boron antimonide, that has never been predicted or experimentally observed for any other material: competing responses of three-phonon and four-phonon interactions to pressure rise cause a non-monotonic pressure dependence of thermal conductivity, κ, which first increases similar to most materials and then decreases. The resulting peak in κ shows a strong temperature dependence from rapid strengthening of four-phonon interactions relative to three-phonon processes with temperature. Our results reveal pressure as a knob to tune the interplay between the competing phonon scattering mechanisms in BAs and similar compounds, and provide clear experimental guidelines for observation in a readily accessible measurement regime.

Thermal properties of materials are driven by complex many-body interactions among thermal atomic vibrations called phonons. Here the authors show, from first-principles, that a full description of the unusual thermal behaviour in boron arsenide requires considering often-neglected four-phonon interactions.

Pressure Impact on the Stability and Distortion of the Crystal Structure of CeScO3

The effects of high pressure on the crystal structure of orthorhombic (Pnma) perovskite-type cerium scandate were studied in situ under high pressure by means of synchrotron X-ray powder diffraction, using a diamond-anvil cell. We found that the perovskite-type crystal structure remains stable up to 40 GPa, the highest pressure reached in the experiments. The evolution of unit-cell parameters with pressure indicated an anisotropic compression. The room-temperature pressure-volume equation of state (EOS) obtained from the experiments indicated the EOS parameters V₀ = 262.5(3) ų, B₀ = 165(7) GPa, and B₀' = 6.3(5). From the evolution of microscopic structural parameters like bond distances and coordination polyhedra of cerium and scandium, the macroscopic behavior of CeScO₃ under compression was explained and reasoned for its large pressure stability. The reported results are discussed in comparison with high-pressure results from other perovskites.

Lattice Thermal Conductivity of MgSiO3 Perovskite from First Principles

We investigate lattice thermal conductivity κ of MgSiO₃ perovskite (pv) by ab initio lattice dynamics calculations combined with exact solution of linearized phonon Boltzmann equation. At room temperature, κ of pristine MgSiO₃ pv is found to be 10.7 W/(m · K) at 0 GPa. It increases linearly with pressure and reaches 59.2 W/(m · K) at 100 GPa. These values are close to multi-anvil press measurements whereas about twice as large as those from diamond anvil cell experiments. The increase of k with pressure is attributed to the squeeze of weighted phase-spaces phonons get emitted or absorbed. Moreover, we find κ exhibits noticeable anisotropy, with κ _zz being the largest component and [Formula: see text] being about 25%. Such extent of anisotropy is comparable to those of upper mantle minerals such as olivine and enstatite. By analyzing phonon mean free paths and lifetimes, we further show that the weak temperature dependence of κ observed in experiments should not be caused by phonons reaching 'minimum' mean free paths. These results clarify the microscopic mechanism of thermal transport in MgSiO₃ pv, and provide reference data for understanding heat conduction in the Earth's deep interior.

Exotic thermoelectric behavior in nitrogenated holey graphene

Due to a high power factor and low thermal conductivity, azTexceeding 5.0 is reached in nitrogenated holey graphene.

Giant Phonon Anharmonicity and Anomalous Pressure Dependence of Lattice Thermal Conductivity in Y2Si2O7 silicate

Modification of lattice thermal conductivity (κL) of a solid by means of hydrostatic pressure (P) has been a crucially interesting approach that targets a broad range of advanced materials from thermoelectrics and thermal insulators to minerals in mantle. Although it is well documented knowledge that thermal conductivity of bulk materials normally increase upon hydrostatic pressure, such positive relationship is seriously challenged when it comes to ceramics with complex crystal structure and heterogeneous chemical bonds. In this paper, we predict an abnormally negative trend dκL/dP < 0 in Y2Si2O7 silicate using density functional theoretical calculations. The mechanism is disclosed as combined effects of slightly decreased group velocity and significantly augmented scattering of heat-carrying acoustic phonons in pressured lattice, which is originated from pressure-induced downward shift of low-lying optic and acoustic phonons. The structural origin of low-lying optic phonons as well as the induced phonon anharmonicity is also qualitatively elucidated with respect to intrinsic bonding heterogeneity of Y2Si2O7. The present results are expected to bring deeper insights for phonon engineering and modulation of thermal conductivity in complex solids with diverging structural flexibility, enormous bonding heterogeneity, and giant phonon anharmonicity.

Anharmonic force constants extracted from first-principles molecular dynamics: applications to heat transfer simulations

A systematic method to calculate anharmonic force constants of crystals is presented. The method employs the direct-method approach, where anharmonic force constants are extracted from the trajectory of first-principles molecular dynamics simulations at high temperature. The method is applied to Si where accurate cubic and quartic force constants are obtained. We observe that higher-order correction is crucial to obtain accurate force constants from the trajectory with large atomic displacements. The calculated harmonic and anharmonic force constants are, then, combined with the Boltzmann transport equation (BTE) and non-equilibrium molecular dynamics (NEMD) methods in calculating the thermal conductivity. The BTE approach successfully predicts the lattice thermal conductivity of bulk Si, whereas NEMD shows considerable underestimates. To evaluate the linear extrapolation method employed in NEMD to estimate bulk values, we analyze the size dependence in NEMD based on BTE calculations. We observe strong nonlinearity in the size dependence of NEMD in Si, which can be ascribed to acoustic phonons having long mean-free-paths and carrying considerable heat. Subsequently, we also apply the whole method to a thermoelectric material Mg2Si and demonstrate the reliability of the NEMD method for systems with low thermal conductivities.

Effect of mass disorder on the lattice thermal conductivity of MgO periclase under pressure

Thermal conductivity of mantle materials controlling the heat balance and thermal evolution of the Earth remains poorly constrained as the available experimental and theoretical techniques are limited in probing minerals under the relevant conditions. We report measurements of thermal conductivity of MgO at high pressure up to 60 GPa and 300 K via diamond anvil cells using the time-domain thermoreflectance technique. These measurements are complemented by model calculations which take into account the effect of temperature and mass disorder of materials within the Earth. Our model calculations agree with the experimental pressure dependencies at 300 and 2000 K for MgO. Furthermore, they predict substantially smaller pressure dependence for mass disordered materials as the mechanism of scattering changes. The calculated thermal conductivity at the core-mantle boundary is smaller than the majority of previous predictions resulting in an estimated total heat flux of 10.4 TW, which is consistent with modern geomodeling estimates.