Insights into magma ocean dynamics from the transport properties of basaltic melt.

The viscosity of magma plays a crucial role in the dynamics of the Earth: from the crystallization of a magma ocean during its initial stages to modern-day volcanic processes. However, the pressure-dependence behavior of viscosity at high pressure remains controversial. In this study, we report the results of first-principles molecular dynamics simulations of basaltic melt to show that the melt viscosity increases upon compression along each isotherm for the entire lower mantle after showing minima at ~6 GPa. However, elevated temperatures of the magma ocean translate to a narrow range of viscosity, i.e., 0.01-0.03 Pa.s. This low viscosity implies that the crystallization of the magma ocean could be complete within a few million years. These results also suggest that the crystallization of the magma ocean is likely to be fractional, thus supporting the hypothesis that present-day mantle heterogeneities could have been generated during the early crystallization of the primitive mantle.

Behavior and properties of water in silicate melts under deep mantle conditions.

Water (H2O) as one of the most abundant fluids present in Earth plays crucial role in the generation and transport of magmas in the interior. Though hydrous silicate melts have been studied extensively, the experimental data are confined to relatively low pressures and the computational results are still rare. Moreover, these studies imply large differences in the way water influences the physical properties of silicate magmas, such as density and electrical conductivity. Here, we investigate the equation of state, speciation, and transport properties of water dissolved in Mg1-xFexSiO3 and Mg2(1-x)Fe2xSiO4 melts (for x = 0 and 0.25) as well as in its bulk (pure) fluid state over the entire mantle pressure regime at 2000-4000 K using first-principles molecular dynamics. The simulation results allow us to constrain the partial molar volume of the water component in melts along with the molar volume of pure water. The predicted volume of silicate melt + water solution is negative at low pressures and becomes almost zero above 15 GPa. Consequently, the hydrous component tends to lower the melt density to similar extent over much of the mantle pressure regime irrespective of composition. Our results also show that hydrogen diffuses fast in silicate melts and enhances the melt electrical conductivity in a way that differs from electrical conduction in the bulk water. The speciation of the water component varies considerably from the bulk water structure as well. Water is dissolved in melts mostly as hydroxyls at low pressure and as -O-H-O-, -O-H-O-H- and other extended species with increasing pressure. On the other hand, the pure water behaves as a molecular fluid below 15 GPa, gradually becoming a dissociated fluid with further compression. On the basis of modeled density and conductivity results, we suggest that partial melts containing a few percent of water may be gravitationally trapped both above and below the upper mantle-transition region. Moreover, such hydrous melts can give rise to detectable electrical conductance by means of electromagnetic sounding observations.

A magma ocean origin to divergent redox evolutions of rocky planetary bodies and early atmospheres.

Magma oceans were once ubiquitous in the early solar system, setting up the initial conditions for different evolutionary paths of planetary bodies. In particular, the redox conditions of magma oceans may have profound influence on the redox state of subsequently formed mantles and the overlying atmospheres. The relevant redox buffering reactions, however, remain poorly constrained. Using first-principles simulations combined with thermodynamic modeling, we show that magma oceans of Earth, Mars, and the Moon are likely characterized with a vertical gradient in oxygen fugacity with deeper magma oceans invoking more oxidizing surface conditions. This redox zonation may be the major cause for the Earth's upper mantle being more oxidized than Mars' and the Moon's. These contrasting redox profiles also suggest that Earth's early atmosphere was dominated by CO2 and H2O, in contrast to those enriched in H2O and H2 for Mars, and H2 and CO for the Moon.

Evidence for Fe-Si-O liquid immiscibility at deep Earth pressures.

Seismic observations suggest that the uppermost region of Earth's liquid outer core is buoyant, with slower velocities than the bulk outer core. One possible mechanism for the formation of a stably stratified layer is immiscibility in molten iron alloy systems, which has yet to be demonstrated at core pressures. We find immiscibility between liquid Fe-Si and Fe-Si-O persisting to at least 140 GPa through a combination of laser-heated diamond-anvil cell experiments and first-principles molecular dynamics simulations. High-pressure immiscibility in the Fe-Si-O system may explain a stratified layer atop the outer core, complicate differentiation and evolution of the deep Earth, and affect the structure and intensity of Earth's magnetic field. Our results support silicon and oxygen as coexisting light elements in the core and suggest that [Formula: see text] does not crystallize out of molten Fe-Si-O at the core-mantle boundary.

Transport properties of carbonated silicate melt at high pressure.

Carbon dioxide, generally considered as the second most abundant volatile component in silicate magmas, is expected to significantly influence various melt properties. In particular, our knowledge about its dynamical effects is lacking over most of Earth's mantle pressure regime. Here, we report the first-principles molecular dynamics results on the transport properties of carbonated MgSiO3 liquid under conditions of mantle relevance. They show that dissolved CO2 systematically enhances the diffusion rates of all elements and lowers the melt viscosity on average by factors of 1.5 to 3 over the pressure range considered. It is remarkable that CO2 has very little or no influence on the electrical conductivity of the silicate melt under most conditions. Simulations also predict anomalous dynamical behavior, increasing diffusivity and conductivity and decreasing viscosity with compression in the low-pressure regime. This anomaly and the concomitant increase of pressure and temperature with depth together make these transport coefficients vary modestly over extended portions of the mantle regime. It is possible that the melt electrical conductivity under conditions corresponding to the 410- and 660-km seismic discontinuities is at a detectable level by electromagnetic sounding observation. In addition, the low melt viscosity values of 0.2 to 0.5 Pa⋅s at these depths and near the core-mantle boundary may imply high mobility of possible melts in these regions.

Carbon-bearing silicate melt at deep mantle conditions.

Knowledge about the incorporation and role of carbon in silicate magmas is crucial for our understanding of the deep mantle processes. CO2 bearing silicate melting and its relevance in the upper mantle regime have been extensively explored. Here we report first-principles molecular dynamics simulations of MgSiO3 melt containing carbon in three distinct oxidation states - CO2, CO, and C at conditions relevant for the whole mantle. Our results show that at low pressures up to 15 GPa, the carbon dioxide speciation is dominated by molecular form and carbonate ions. At higher pressures, the dominant species are silicon-polyhedral bound carbonates, tetrahedral coordination, and polymerized di-carbonates. Our results also indicate that CO2 component remains soluble in the melt at high pressures and the solution is nearly ideal. However, the elemental carbon and CO components show clustering of carbon atoms in the melt at high pressures, hinting towards possible exsolution of carbon from silicate melt at reduced oxygen contents. Although carbon lowers the melt density, the effect is modest at high pressures. Hence, it is likely that silicate melt above and below the mantle transition zone, and atop the core-mantle boundary could efficiently sequester significant amounts of carbon without being gravitationally unstable.

Solid-liquid density and spin crossovers in (Mg, Fe)O system at deep mantle conditions.

The low/ultralow-velocity zones in the Earth's mantle can be explained by the presence of partial melting, critically depending on density contrast between the melt and surrounding solid mantle. Here, first-principles molecular dynamics simulations of (Mg, Fe) O ferropericlase in the solid and liquid states show that their densities increasingly approach each other as pressure increases. The isochemical density difference between them diminishes from 0.78 (±0.7) g/cm3 at zero pressure (3000 K) to 0.16 (±0.04) g/cm3 at 135 GPa (4000 K) for pure and alloyed compositions containing up to 25% iron. The simulations also predict a high-spin to low-spin transition of iron in the liquid ferropericlase gradually occurring over a pressure interval centered at 55 GPa (4000 K) accompanied by a density increase of 0.14 (±0.02) g/cm3. Temperature tends to widen the transition to higher pressure. The estimated iron partition coefficient between the solid and liquid ferropericlase varies from 0.3 to 0.6 over the pressure range of 23 to 135 GPa. Based on these results, an excess of as low as 5% iron dissolved in the liquid could cause the solid-liquid density crossover at conditions of the lowermost mantle.

Structure and density of basaltic melts at mantle conditions from first-principles simulations.

The origin and stability of deep-mantle melts, and the magmatic processes at different times of Earth's history are controlled by the physical properties of constituent silicate liquids. Here we report density functional theory-based simulations of model basalt, hydrous model basalt and near-MORB to assess the effects of iron and water on the melt structure and density, respectively. Our results suggest that as pressure increases, all types of coordination between major cations and anions strongly increase, and the water speciation changes from isolated species to extended forms. These structural changes are responsible for rapid initial melt densification on compression thereby making these basaltic melts possibly buoyantly stable at one or more depths. Our finding that the melt-water system is ideal (nearly zero volume of mixing) and miscible (negative enthalpy of mixing) over most of the mantle conditions strengthens the idea of potential water enrichment of deep-mantle melts and early magma ocean.

Thermal conductivity of periclase (MgO) from first principles.

We combine first-principles calculations of forces with the direct nonequilibrium molecular dynamics method to determine the lattice thermal conductivity k of periclase (MgO) up to conditions representative of the Earth's core-mantle boundary (136 GPa, 4100 K). We predict the logarithmic density derivative a=(∂ln k/∂ln ρ)(T)=4.6±1.2 and that k=20±5 Wm(-1) K(-1) at the core-mantle boundary, while also finding good agreement with extant experimental data at much lower pressures.

First-principles study of enhancement of transport properties of silica melt by water.

First-principles molecular dynamics simulations show that water (8.25 wt%) dramatically affects the transport properties of SiO2 liquid increasing the diffusivity and decreasing the viscosity by an order of magnitude. At 3000 K, the diffusivity of Si, O, and H, and the viscosity vary anomalously with pressure. Highly mobile protons make the hydrous liquid a potential superionic conductor. The predicted dynamical changes are associated with structural depolymerization and water speciation, which changes from being dominated by hydroxyls at low pressure to extended structures at high pressure.

Viscosity of MgSiO3 liquid at Earth's mantle conditions: implications for an early magma ocean.

Understanding the chemical and thermal evolution of Earth requires knowledge of transport properties of silicate melts at high pressure and high temperature. Here, first-principles molecular dynamics simulations show that the viscosity of MgSiO3 liquid varies by two orders of magnitude over the mantle pressure regime. Addition of water systematically lowers the viscosity, consistent with enhanced structural depolymerization. The combined effects of pressure and temperature along model geotherms lead to a 10-fold increase in viscosity with depth from the surface to the base of the mantle. Based on these calculations, efficient heat flux from a deep magma ocean may have exceeded the incoming solar flux early in Earth's history.

A framework for detecting glaucomatous progression in the optic nerve head of an eye using proper orthogonal decomposition.

Glaucoma is the second leading cause of blindness worldwide. Often, the optic nerve head (ONH) glaucomatous damage and ONH changes occur prior to visual field loss and are observable in vivo. Thus, digital image analysis is a promising choice for detecting the onset and/or progression of glaucoma. In this paper, we present a new framework for detecting glaucomatous changes in the ONH of an eye using the method of proper orthogonal decomposition (POD). A baseline topograph subspace was constructed for each eye to describe the structure of the ONH of the eye at a reference/baseline condition using POD. Any glaucomatous changes in the ONH of the eye present during a follow-up exam were estimated by comparing the follow-up ONH topography with its baseline topograph subspace representation. Image correspondence measures of L1-norm and L2 -norm, correlation, and image Euclidean distance (IMED) were used to quantify the ONH changes. An ONH topographic library built from the Louisiana State University Experimental Glaucoma study was used to evaluate the performance of the proposed method. The area under the receiver operating characteristic curves (AUCs) was used to compare the diagnostic performance of the POD-induced parameters with the parameters of the topographic change analysis (TCA) method. The IMED and L2-norm parameters in the POD framework provided the highest AUC of 0.94 at 10 degrees field of imaging and 0.91 at 15 degrees field of imaging compared to the TCA parameters with an AUC of 0.86 and 0.88, respectively. The proposed POD framework captures the instrument measurement variability and inherent structure variability and shows promise for improving our ability to detect glaucomatous change over time in glaucoma management.

Atomistic visualization: space-time multiresolution integration of data analysis and rendering.

Time-varying three-dimensional scattered data representing snapshots of atomic configurations produced by molecular dynamics simulations are not illuminating by themselves; gaining insight into them poses a tremendous challenge. In order to take the advantage of maximal information offered by these simulations, we have proposed an efficient scheme, which integrates various analysis and rendering tasks together in order to support interactive visualization of the data at space-time multiresolution. Additional data produced by various analytical techniques on the fly represent the atomic system under consideration at diverse length- (e.g., nearest neighbor, next-nearest neighbor or beyond) and time- (e.g., instantaneous, finite intervals or overall averages) scales. In particular, the radial distribution functions, coordination environments, clusters and rings are computed and visualized to understand the structural behavior whereas a variety of displacement data and covariance matrices are explored to understand the dynamical behavior. While the spatial distributions of atoms need to be reproduced correctly during rendering, we take the advantage of high flexibility in rendering other attributes because of the lack of their direct physical relevance. A combination of different techniques including animation, color maps, pathlines, different types of glyphs, and graphics hardware accelerated approach is exploited to render the original and extracted data. First-principles molecular dynamics simulation data for liquid systems are used to justify the effectiveness and usefulness of the proposed scheme.

Large-scale molecular dynamics simulations of alkanethiol self-assembled monolayers.

Large-scale molecular dynamics simulations of self-assembled alkanethiol monolayer systems have been carried out using an all-atom model involving a million atoms to investigate their structural properties as a function of temperature, lattice spacing, and molecular chain length. Our simulations show that the alkanethiol chains of 13-carbons tilt from the surface normal by a collective angle of 25 degrees along next-nearest-neighbor direction at 300 K. The tilt structure of 13-carbon alkanethiol system is found to depend strongly on temperature and exhibits hysteresis. At 350 K the 13-carbon alkanethiol system transforms to a disordered phase characterized by small collective tilt angle, flexible tilt direction, and random distribution of backbone planes. The tilt structure also depends on lattice spacing: With increasing lattice spacing a the tilt angle increases rapidly from a nearly zero value at a=4.7 A to as high as 34 degrees at a=5.3 A at 300 K for 13-carbon alkanethiol system. Finally, the effects of the molecular chain length on the tilt structure are significant at high temperatures.