Insights into core-mantle differentiation from bulk Earth melt simulations.

The earth is thought to have gone through complex physicochemical changes during the accretion and magma ocean stages. To better understand this evolution process at the fundamental level, we investigate the behavior of a bulk earth melt system by simulating the composition FeMgSiO (in wt%) at high pressure. A deep neural network potential trained by first-principles data can enable accurate molecular dynamics simulation of large supercells that greatly enhances sampling for reliable evaluation of elemental partitioning.

Optimal surface-tension isotropy in the Rothman-Keller color-gradient lattice Boltzmann method for multiphase flow.

The Rothman-Keller color-gradient (CG) lattice Boltzmann method is a popular method to simulate two-phase flow because of its ability to deal with fluids with large viscosity contrasts and a wide range of interfacial tensions. Two fluids are labeled red and blue, and the gradient in the color difference is used to compute the effect of interfacial tension. It is well known that finite-difference errors in the color-gradient calculation lead to anisotropy of interfacial tension and errors such as spurious currents.

Evidence of lower-mantle slab penetration phases in plate motions.

It is well accepted that subduction of the cold lithosphere is a crucial component of the Earth's plate tectonic style of mantle convection. But whether and how subducting plates penetrate into the lower mantle is the subject of continuing debate, which has substantial implications for the chemical and thermal evolution of the mantle. Here we identify lower-mantle slab penetration events by comparing Cenozoic plate motions at the Earth's main subduction zones with motions predicted by fully dynamic models of the upper-mantle phase of subduction, driven solely by downgoing plate density.