# Pressure-Induced Coordination Changes in a Pyrolitic Silicate Melt From Ab Initio Molecular Dynamics Simulations

Natalia V. Solomatova, Razvan Caracas

► To cite this version:

Natalia V. Solomatova, Razvan Caracas. Pressure-Induced Coordination Changes in a Pyrolitic Silicate Melt From Ab Initio Molecular Dynamics Simulations. Journal of Geophysical Research: Solid Earth, 2019, 124 (11), pp.11232-11250. 10.1029/2019JB018238. hal-04979713

HAL Id: hal-04979713

https://hal.science/hal-04979713v1

Submitted on 21 Mar 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# JGR Solid Earth

## RESEARCH ARTICLE

10.1029/2019JB018238

**Key Points:**

* The proportion and distribution of cation-to-oxygen coordination environments in a pyrolite melt have been determined as a function of depth

* We examine the distribution of lifetimes of cation-oxygen species at select pressure-temperature conditions

* Our results are compared to previous ab initio calculations on mafic silicate melts and to experiments on various silicate phases

**Supporting Information:**

* Supporting Information S1

* Data Set S1

**Correspondence to:**

N. V. Solomatova,
nsolomat@gmail.com

**Citation:**

Solomatova, N. V., & Caracas, R. (2019). Pressure-induced coordination changes in a pyrolitic silicate melt from ab initio molecular dynamics simulations.

*Journal of Geophysical Research: Solid Earth, 124*, 11,232–11,250. https://doi.org/10.1029/2019JB018238

Received 19 JUN 2019

Accepted 29 OCT 2019

Accepted article online 11 NOV 2019

Published online 29 NOV 2019

©2019. The Authors.
This is an open access article under the
terms of the Creative Commons
Attribution License, which permits use,
distribution and reproduction in any
medium, provided the original work is
properly cited.

<figure><img src="image_2.png" /></figure>

## Pressure-Induced Coordination Changes in a Pyrolitic Silicate Melt From Ab Initio Molecular Dynamics Simulations

N.V. Solomatova¹ ID and R. Caracas¹ ID

¹CNRS, École Normale Supérieure de Lyon, Laboratoire de Géologie de Lyon, CNRS UMR 5276, Lyon, France

**Abstract** With ab initio molecular dynamics simulations on a Na-, Ca-, Fe-, Mg-, and Al-bearing silicate melt of pyrolite composition, we examine the detailed changes in elemental coordination as a function of pressure and temperature. We consider the average coordination as well as the proportion and distribution of coordination environments at pressures and temperatures encompassing the conditions at which molten silicates may exist in present-day Earth and those of the Early Earth's magma ocean. At ambient pressure and 2,000 K, we find that the average coordination of cations with respect to oxygen is 4.0 for Si-O, 4.0 for Al-O, 3.7 for Fe-O, 4.6 for Mg-O, 5.9 for Na-O, and 6.2 for Ca-O. Although the coordination for iron with respect to oxygen may be underestimated, the coordination number for all other cations are consistent with experiments. By 15 GPa (2,000 K), the average coordination for Si-O remains at 4.0 but increases to 4.1 for Al-O, 4.2 for Fe-O, 4.9 for Mg-O, 8.0 for Na-O, and 6.8 for Ca-O. The coordination environment for Na-O remains approximately constant up to core-mantle boundary conditions (135 GPa and 4000 K) but increases to about 6 for Si-O, 6.5 for Al-O, 6.5 for Fe-O, 8 for Mg-O, and 9.5 for Ca-O. We discuss our results in the context of the metal-silicate partitioning behavior of siderophile elements and the viscosity changes of silicate melts at upper mantle conditions. Our results have implications for melt properties, such as viscosity, transport coefficients, thermal conductivities, and electrical conductivities, and will help interpret experimental results on silicate glasses.

## 1. Introduction

The Early Earth's mantle was once composed of molten silicate rock—a global magma ocean extending down to the core-mantle boundary (Stevenson, 1989; Tonks & Melosh, 1993). Today, silicate melts exist as lava flows on Earth's surface, as magmas in subterrestrial volcanic processes, and are thought to exist at various depths of the mantle, particularly at the base of the upper mantle, at the base of the transition zone, and at the core-mantle boundary. A seismically detectable low velocity layer observed at the base of the upper mantle is often attributed to the presence of melt within a narrow ~50-km thick layer (Revenaugh & Sipkin, 1994; Song et al., 2004; Tauzin et al., 2010; Vinnik & Farra, 2007; Freitas et al., 2017). It has also been proposed that dehydration melting occurs at the top of the lower mantle (Schmandt et al., 2014). Silicate melts in the mantle strongly affect seismic velocities, decreasing seismic velocities by several percent (Hammond & Humphreys, 2000). In the deepest part of the lower mantle, at the boundary between Earth's liquid iron-nickel core and solid D" layer, it is thought that pockets of silicate melts may exist, explaining some of the seismically observed ultra-low velocity zones (Vidale & Hedlin, 1998; Williams & Garnero, 1996), though it is not yet clear if they would be gravitationally stable (Bower et al., 2011).

The properties of a melt are strongly dependent on the coordination and geometric packing of its ions. Thus understanding the structure of silicate melts is essential to gain a better understanding of both the Early Earth's magma ocean and for interpreting present-day experimentally measured properties of amorphous silicates and seismically observed anomalies caused by melts.

Here we conducted simulations at pressures and temperatures relevant to conditions of both the present-day Earth and the Early Earth's magma ocean. The solidus of dry pyrolite-like compositions is roughly 1,400–1,500 K at the surface, rising rapidly to 2,000 K by 5 GPa and to 2,500 K by 25 GPa (Hirschmann, 2000; Litasov & Ohtani, 2002), while the solidus of pyrolite at midmantle pressures of 40–80 GPa is about 3,000 K rising to 3,500–4,000 K at core-mantle boundary pressures (Fiquet et al., 2010; Nomura et al., 2014; Zerr et al., 1998). Thus, our simulation temperatures of 2,000–4,000 K encompass the temperatures at which
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Pyrolite melt at 4,000 K at (a) ~0 GPa and (b) 150 GPa, the highest pressure achieved. SiOₓ and AlOₓ are shown as light gray and light blue polyhedra, respectively, highlighting the transition from predominantly fourfold coordination at low pressure to predominantly sixfold coordination at high pressure.</figcaption></figure>

melt could occur in the entire present-day mantle. In the Early Earth's magma ocean, the temperature of the molten silicate mantle may have been above the liquidus, which reaches 5,000 K at the base of the lower mantle (Fiquet et al., 2010), and so our calculations at 5,000 K are especially relevant for determining the structure of silicate melts during the early magma ocean stage of Earth.

## 2. Computational Methods

Ab initio molecular dynamics simulations were performed with density functional theory and with the Vienna ab initio simulation package (Kresse & Furthmüller, 1996). We used the projector-augmented wave method (Blöchl, 1994) to represent the core electrons and the generalized gradient approximation in the Perdew-Burke-Ernzerhof form (Perdew et al., 1996) to treat electron exchange and correlation. The electronic structure of iron is calculated at each time step and the d electrons were treated as unconstrained spin polarized at all temperatures and pressures; in other words, the spin state was not enforced to either high or low spin. A Hubbard *U* parameter, which is sensitive to the electronic and geometric environment of iron, was not applied. In previous tests on iron-bearing pyrolite melt, a *U* value of 4 eV resulted in an increase in the magnetic moment of the iron atoms with only a minor effect on pressure (<1 GPa; Caracas et al., 2019). The kinetic energy cutoffs for the plane-wave expansion of the wavefunctions and the augmentation charges were set to 550 and 800 eV, respectively. The canonical ensemble with a constant number of atoms, volume, and temperature was used with a Nose-Hoover thermostat (Nosé, 1984; Hoover, 1985). Simulations were conducted at 2,000, 3,000, 4,000, and 5,000 K and at volumes corresponding to pressures of about 0-40 GPa at 2,000 K and pressures of about 0-160 GPa at 3,000-5,000 K (negative average pressures are achieved when the melt is stretched to a volume slightly beyond the zero-pressure volume). The Brillouin zone was sampled at the gamma point. Simulations were performed with a time step of 0.5-2 fs for 10-50 picoseconds, depending on the temperature and density. The mean-square displacement as a function of time shows a ballistic regime below approximately 1,000 fs, after which the atoms reach a diffusive regime.

The composition of the bulk silicate Earth (McDonough & Sun, 1995) was modeled with a pyrolite melt with the stoichiometry NaCa₂Fe₄Mg₃₀Al₃Si₂₄O₈₉ (Figure 1; Table 1). Bond distances were determined from the pair distribution functions (Supporting Information Figure S1), which describe the probability of finding an atom type at a given distance from the reference atom. We use the first peak in the pair distribution function to approximate the average bond length; the distance at which the first minimum occurs marks the radius of the first coordination sphere of atoms that are directly bonded to the reference atom. We use this value to define the bond criterion between two atom types. The maximum acceptable bonds were then fitted with a third-order polynomial as a function of supercell lattice parameter (Figure S2) to smooth possible statistical irregularities. We chose to fit the maximum bonds as a function of supercell lattice parameter rather
<table><caption>**Table 1**

Composition of the Pyrolite Melt in Oxide Weight Percent Used for This Study Compared to the Pyrolite Model by McDonough and Sun (1995)</caption><thead><tr><th></th><th>SiO<sub>2</sub></th><th>Al<sub>2</sub>O<sub>3</sub></th><th>FeO</th><th>MgO</th><th>CaO</th><th>Na<sub>2</sub>O</th></tr></thead><tbody><tr><td>This study</td><td>44.5</td><td>4.72</td><td>9.00</td><td>37.3</td><td>3.46</td><td>0.96</td></tr><tr><td>Pyrolite model</td><td>45.0</td><td>4.45</td><td>8.05</td><td>37.8</td><td>3.55</td><td>0.36</td></tr></tbody></table>

Note. The pyrolite model also includes 0.201 wt % TiO₂, 0.384 wt % Cr₂O₃, 0.135 wt % MnO, 0.25 wt % NiO, 0.029 wt % K₂O, and 0.021 wt % P₂O₅, which are not included in our simulations.

than pressure or density due to the even distribution of data points. The fitted values were then used in all speciation analysis. Whether the maximum bonds were fitted or not resulted in a negligible difference in the speciation statistics due to the fact that, by definition, the majority of bonds exist near the peak of the pair distribution function rather than the first minimum.

## 3. Results

Using ab initio molecular dynamics simulations on pyrolite melt (i.e., a melt with approximately the composition of the bulk silicate Earth) up to core-mantle boundary pressures and temperatures of 2,000–5000 K, we compute average bond lengths, the distribution of bond lengths, and changes in cation-oxygen coordination and compare our results to previous experimental and computational studies on amorphous silicates of various compositions.

### 3.1. Bond Lengths

The first maximum in the pair distribution function corresponds to the most frequent bond distance that occurs in the simulation and approximately to the average bond length. The maximum acceptable bond distances were determined from the first minimum in the pair distribution function (Figure S2). Using the aforementioned definition, average bond lengths at ambient pressure and 2,000 K are 1.64 Å for Si-O, 1.77 Å for Al-O, 1.89 Å for Fe-O, 1.98 Å for Mg-O, 2.30 Å for Na-O, and 2.30 Å for Ca-O. At upper mantle conditions, the average cation-oxygen bond lengths for Si-O, Al-O, and Mg-O gently increase with increasing pressure while the average bond lengths for Fe-O and Ca-O increase with increasing pressure more steeply (Figure 2). In the same pressure range, the average bond lengths for Na-O generally decrease with increasing pressure. At transition-zone and lower-mantle pressures (i.e., above roughly 15 GPa), the bond lengths for all cation-oxygen pairs decrease monotonically with increasing pressure (Figure 2). The slope of the change in bond distance with depth is dependent on a balance between the increase in pressure and increase in coordination. At low pressures, bond distances increase slightly with increasing pressure because the

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Average bond lengths as a function of pressure, as approximated by the position of the first peak of the corresponding cation-oxygen pair distribution functions.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Average coordination numbers of cations with respect to oxygen as a function of pressure for cation-oxygen pairs: (a) Si-O, (b) Al-O, (c) Mg-O, (d) Fe-O, (e) Ca-O, and (f) Na-O. Solid lines with solid symbols are the results of this study where blue is 2,000 K, green is 3,000 K, yellow-orange is 4,000 K and red is 5,000 K. Ab initio molecular dynamics simulations on midocean ridge basalt (MORB) melt are shown as a gray dashed line, as reported averaged over temperatures of 1,800–4,000 K and pressures at each volume (Bajgain et al., 2015). Ab initio molecular dynamics simulations at 3,000 K on CaMgSi₂O₆ and Mg₂SiO₄ melts are shown in green dashed-dotted and dotted curves, respectively (Sun et al., 2011; de Koker et al., 2008). The average coordination of silicon from ab initio calculations on MgSiO₃ melt at 3,000 K (Stixrude & Karki, 2005) overlaps nearly perfectly with the silicon's average coordination in CaMgSi₂O₆ (Sun et al., 2011) and so is not shown here for clarity.</figcaption></figure>

coordination number and polyhedra size increase very rapidly (section 3.2) while at high pressures, bond distances decrease with increasing pressure as the ions get closer together and the polyhedra shrink in size.

## 3.2. Atomic Coordinations

The average coordination of all cations with respect to oxygen generally increases with increasing pressure. At upper mantle conditions, calcium and sodium undergo sharp increases in coordination while silicon and aluminum undergo more gentle increases in coordination and at lower mantle pressures, the coordination environments of calcium and sodium remain roughly constant. The average coordination of all cations with respect to oxygen at all pressures and temperatures, as well as the errors associated with pressure are provided in Tables S1–S4.

### 3.2.1. Silicon-Oxygen

As the most abundant cation in silicate melts, the chemical environment of silicon is most strongly linked to the properties of the melt, such as density and viscosity, and it is therefore essential to understand changes in coordination of silicon as a function of depth. At ambient pressure and 2,000 K, we find that silicon has an
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Atomic coordination of (a) silicon, (b) aluminum, (c) magnesium, (d) and iron to oxygen at 4,000 K as a function of density (corresponding pressures are provided on top). Coordination numbers are labeled next to each curve. See Supporting Information Figure S4 for the atomic coordination of silicon, aluminum, and magnesium at 2,000 K and 3,000 K compared to 4,000 K.</figcaption></figure>

average coordination of 4, existing almost entirely as SiO₄ tetrahedra (Figure 3a). With increasing temperature and constant ambient pressure, the average coordination approaches 3 as the proportion of fourfold-coordinated silicon decreases at the expense of threefold-coordinated silicon. The proportions of the most common coordination environments of silicon with respect to oxygen are plotted in Figure 4a at 4,000 K and at all temperatures in Figures S3 and S4. Only about 50% of silicon has fourfold coordination at 5,000 K compared to nearly 100% at 2,000 K. Unlike at 2,000 K, where the average coordination of silicon remains constant at pressures of 0–5 GPa, at 5,000 K, the average coordination increases rapidly from threefold to fourfold within the first 5 GPa. Overall, with increasing temperature, the distribution of coordination environments increases while the average coordination decreases (Figure 5a). Fitting a Gaussian curve to the coordination numbers at 2,000 K produces a full width at half maximum (FWHM) of 0.8, ranging from 3.6 to 4.4, whereas at 5,000 K, the full width at half maximum increases to 1.5, ranging from 2.9 to 4.4.

With increasing pressure, the average coordination number of silicon increases with a general trend from fourfold at low pressure to fivefold at 40–50 GPa to sixfold at core-mantle boundary pressures. At 5,000 K, the proportion of fourfold-coordinated silicon increases with increasing pressure up until 5 GPa, after which, it steadily decreases. Meanwhile, the proportion of threefold-coordinated silicon decreases with pressure, reaching 0% by about 50 GPa, and the proportions of fivefold-coordinated silicon steeply increases with pressure. Sixfold-coordinated silicon remains under 1% abundance until 10 GPa and under 5% abundance until 20 GPa for all temperatures, after which, its abundance increases rapidly. At core-mantle boundary pressures, the proportion of sixfold silicon is about 80% at 3,000 K, 70% at 4,000 K, and 55% at 5,000 K.

### 3.2.2. Aluminum-Oxygen

At ambient pressure, despite having a larger average bond length compared to silicon, aluminum displays a similar coordination environment to silicon, existing predominantly in fourfold coordination at 2,000 K and decreasing slightly in coordination and increasing in distribution of coordination environments at higher
<figure><img src="image_7.png" /><figcaption>**Figure 5.** The distribution of coordination environments for (a) silicon and (b) aluminum with respect to oxygen at pressures of 0 GPa (ambient pressure), 14 GPa (base of the upper mantle, UM), 24 GPa (base of the transition zone, TZ), and 135 GPa (core-mantle boundary, CMB) and temperatures of 2,000–5,000 K. Calculations were performed at constant volume and temperature (see computational methods in section 2 for more details); to determine the coordination numbers at exactly the same pressures for all temperatures, the coordination numbers were locally fitted along isotherms as a function of pressure with a third-order polynomial.</figcaption></figure>

temperatures. However, with increasing pressure, the proportion of fourfold-coordinated aluminum decreases more rapidly than the proportion of fourfold-coordinated silicon (Figure 3b). By 40–50 GPa, there is no fourfold-coordinated aluminum left in the melt, whereas a small fraction of fourfold-coordinated silicon persists up to lowermost mantle pressures. The proportion of fivefold-coordinated aluminum increases with increasing pressure reaching a maximum of about 50% at 10–20 GPa (depending on the temperature), after which, the population of fivefold-coordinated aluminum decreases to 5% by core-mantle-boundary pressures (Figure 4b). The proportions of sixfold and sevenfold aluminum increase with increasing pressure. Between 50 and 100 GPa, roughly 50% of aluminum is sixfold coordinated, 25% is fivefold coordinated, and 25% is sevenfold coordinated, where differences in temperature appear to be mostly a result of scatter due to the small quantity of aluminum atoms in the melt. The proportions of the most common coordination environments for aluminum with respect to oxygen are plotted as a function of pressure in Figure S4, while the distribution of the coordination environments at select pressures are plotted in Figure 5b.

### 3.2.3. Magnesium-Oxygen

We find that at ambient pressure and 2,000 K, the average coordination of magnesium in pyrolite melt is 4.6 (Figure 3c), existing as a mixture of threefold (10%), fourfold (40%), fivefold (40%), and sixfold (10%) coordination (Figures S4 and S6). With increasing pressure, the average coordination of magnesium increases to sixfold by 20 GPa, sevenfold by 45 GPa, and nearly eightfold by 135 GPa. There is generally a negative correlation between temperature and average coordination number. At ambient pressure, average coordination of magnesium decreases from 4.6 at 2,000 K to 4 at 5,000 K while at 40 GPa, it decreases from about 7 at 2,000 K down to 6.6 at 5,000 K (Figure 3c). The distribution of environments is also slightly and positively correlated with temperature; interestingly though, at core-mantle boundary pressures, temperature no longer affects the distribution (Figure 6a). The proportion of fourfold-coordinated magnesium decreases rapidly at low pressure to under 5% by 20 GPa, while the concentrations of fivefold, sixfold, and sevenfold coordination increase with pressure and reach a maximum at about 3, 20, and 50 GPa, respectively, after which, their concentrations decrease with increasing pressure as higher coordination begin to dominate (Figures 4c and S6).
<figure><img src="image_8.png" /></figure>

**Figure 6.** The distribution of coordination environments for (a) magnesium (b) and iron with respect to oxygen at pressures of 0 GPa (ambient pressure), 14 GPa (base of the upper mantle, UM), 24 GPa (base of the transition zone, TZ), and 135 GPa (core-mantle boundary, CMB) and temperatures of 2,000–5,000 K. Calculations were performed at constant volume and temperature (see computational methods in section 2 for more details); to determine the coordination numbers at exactly the same pressures for all temperatures, the coordination numbers were locally fitted along isotherms as a function of pressure with a third-order polynomial.

### 3.2.4. Iron-Oxygen

Although it might be expected that the average coordination of iron would be higher than the average coordination of magnesium, ab initio molecular dynamics simulations of silicate melts seem to predict a lower average coordination (Bajgain et al., 2015; this study). Non-spin-polarized calculations on midocean ridge basalt (MORB) melt produced an average iron to oxygen coordination of 3.5 at 0 GPa increasing to 4.5 by about 14 GPa (Bajgain et al., 2015). In comparison, in our spin-polarized calculations on pyrolite melt, the average coordination of iron at 2,000 K increases from 3.7 at 0 GPa to 5.5 at 15 GPa; by 25 GPa, iron is on average sixfold coordinated (Figure 3d). At core-mantle boundary conditions, iron is approximately sixfold to sevenfold coordinated. The proportions of the most common coordination environments for iron with respect to oxygen are plotted as a function of pressure in Figure S7, and the distribution of the coordination environments at select pressures are plotted in Figure 5b. At ambient pressure, iron has a magnetic moment of about 3.5 µB where the iron atoms remain entirely in the high-spin state. With increasing pressure, the magnetic moment of iron gradually decreases. An example of the distribution of local magnetic moments of iron as a function of pressure is shown in Figure S8.

### 3.2.5. Calcium-Oxygen

Having a much larger ionic radius than magnesium or iron, calcium exhibits a higher coordination environment, which also increases in coordination more rapidly with pressure (Figure 3f). At ambient pressure, calcium exhibits a very broad distribution of coordination environments, ranging from fourfold to eightfold at 2,000 K and threefold to sevenfold at 5,000 K (Figure 7a). At the base of the upper mantle, 20% of calcium is sevenfold, 40% is eightfold, 25% is ninefold, and 10% is tenfold (the remaining 5% is fivefold and elevenfold). With increasing pressure, the distribution of coordination environments narrows. For example, at core-mantle boundary pressure, calcium exists mostly in ninefold (30%) and tenfold (45%) coordinations with lesser proportions of eightfold (10%) and elevenfold (15%) coordinations (Figure S9). At 2,000 K, the average coordination of calcium by oxygen increases steeply from 6.3 to 8 with increasing pressure from 0 to 10 GPa while with increasing temperature, the average coordination drops from 6.3 to about 5 at 2,000 to 5,000 K (Figure 3f). The average coordination of calcium increases to 7 by 6 GPa and 8 by 15 GPa, independent of temperature. Above 40 GPa, the average coordination number flattens out at around 9–10 up to core-mantle boundary conditions.
<figure><img src="image_9.png" /></figure>

**Figure 7.** The distribution of coordination environments for (a) calcium (b) and sodium with respect to oxygen at pressures of 0 GPa (ambient pressure), 14 GPa (base of the upper mantle, UM), 24 GPa (base of the transition zone, TZ), and 135 GPa (core-mantle boundary, CMB) and temperatures of 2,000–5,000 K. Calculations were performed at constant volume and temperature (see computational methods in section 2 for more details); to determine the coordination numbers at exactly the same pressures for all temperatures, the coordination numbers were locally fitted along isotherms as a function of pressure with a third-order polynomial.

### 3.2.6. Sodium-Oxygen

Due to the low concentration of sodium (one atom per 153-atom supercell), the statistics for sodium are poor and the scatter in the data is subsequently large, so results should be evaluated cautiously. In general, sodium behaves similarly to calcium. At low pressure, the average coordination ranges from 6 at 2,000 K down to 3.5 at 5,000 K. The average coordination of sodium increases rapidly to 8–9 by 15 GPa, after which, the average coordination plateaus at around 8 (Figure 3e). Similar to calcium, the distribution of coordination environments is very large at low pressure and sharpens at high pressure as the coordination of sodium by oxygen atoms is maximized (Figure 7b). The proportions of the most common coordination environments for sodium with respect to oxygen are plotted as a function of pressure in Figure S10.

### 3.2.7. Lifetimes of Species

In a molecular dynamics simulation, bonds are formed and broken as the ions move through the melt, forming various types of species. Here we focus on cation-anion species, MOₓ (where M is one of the six cations, O is oxygen, and x is the number of oxygen atoms surrounding the cation, as defined by the first coordination shell in the pair distribution function). Each MOₓ species lives for some amount of time before the coordination changes, and we refer to the time they exist for as the “lifetime” of the species. Some species, such as SiO₄ tend to live for very long times (thousands of femtoseconds), while others are much more short lived (tens of femtoseconds). We performed a careful analysis of the lifetime of MOₓ species. In Figure 8, we show the lifetimes of species at ambient pressure (~0 GPa) and 2,000 K. The longest surviving SiO₄ species survives the entirety of the simulation (20,646 fs) with an average SiO₄ lifetime of 2,868 fs. Although aluminum has a similar coordination environment to silicon at these conditions, AlO₄ tetrahedra are not trapped within as deep of a potential minimum as SiO₄, and so the longest surviving AlO₄ species only lives for about a fifth of the simulation time lives (4.805 fs) and has an average lifetime of 331 fs. Behaving more like network modifiers even in low coordination environments, iron and magnesium polyhedra have much shorter lifetimes, ranging from 20 to 65 fs for the more populous threefold to sixfold coordinations while sodium and calcium have even shorter lifetimes of 10–25 fs for the threefold to eightfold coordinations. Na-O and Ca-O bonds break and reform frequently as the relatively weakly bonded oxygen atoms migrate through the melt much more easily than the oxygen atoms bonded to silicon tetrahedra. As temperature increases, the
<figure><img src="image_10.png" /><figcaption>**Figure 8.** Lifetimes of each species at ambient pressure (~0 GPa) and 2,000 K. The x axis represents each individual MOₓ species sorted in decreasing order of lifetimes, and the y axis is the lifetime (the amount of time each species lived for) in femtoseconds. The average length of time that each type of species existed for in femtoseconds is labeled within each figure. The total simulation length was 20,600 fs. For example, within the simulation, there were about 160 SiO₄ species that existed for 0–20,600 fs with an average lifetime of about 2,900 fs. In contrast, there were 5,500 MgO₄ species that existed for only 0–280 fs with average lifetime of about 50 fs.</figcaption></figure>

lifetimes of species decreases. For example, at 0 GPa and 4,000 K, the average lifetime of SiO₄ tetrahedra is 244 fs, more than an order of magnitude lower than at 2,000 K (Figure S11). At midmantle solidus conditions, many different types of coordination environments exist for silicon with respect to oxygen. At 80 GPa and 3,000 K, for example, silica tetrahedra exist for an average of 30 fs while silica octahedra exist for an average of 119 fs (Figure 9), decreasing down to 16 and 43 fs, respectively, at 150 GPa and 5,000 K (Figure S12). Aluminum-oxygen polyhedra, ranging from fourfold to eightfold coordination exist on average for about 20–50 fs at 80 GPa and 3,000 K, decreasing down to 15–20 fs at 150 GPa and 5,000 K. Iron-oxygen and magnesium-oxygen polyhedra with fivefold to eightfold coordination exist on average for 10–30 fs and similarly, sodium-oxygen and calcium-oxygen polyhedra with eightfold to elevenfold coordination exist on average for 10–25 fs at 80 GPa and 3,000 K (Figure 9), decreasing further by 150 GPa and 5,000 K (Figure S10).

## 4. Discussion

We compare our results to ab initio molecular dynamics simulations on a MORB melt (Bajgain et al., 2015), Mg₂SiO₄ melt (de Koker et al., 2008), CaMgSi₂O₆ melt (Sun et al., 2011), and MgSiO₃ melt (Stixrude & Karki, 2005). We note that Bajgain et al. (2015) averaged their coordination numbers over pressures corresponding to each volume and temperatures of 1,800–4,000 K (resulting in an average temperature of ~3,000 K). For
<figure><img src="image_11.png" /><figcaption>**Figure 9.** Lifetimes of each species at 80 GPa and 3,000 K. The x axis represents each individual MOₓ species sorted in decreasing order of lifetimes, and the y axis is the lifetime (the amount of time each species lived for) in femtoseconds. The average length of time that each type of species existed for in femtoseconds is labeled within each figure. The total simulation length was 23.9 picoseconds.</figcaption></figure>

Mg₂SiO₄ (de Koker et al., 2008) and CaMgSi₂O₆ (Sun et al., 2011), we compare results at 3,000 K. The average cation-anion coordinations are systematically higher in MORB melt (Bajgain et al., 2015) compared to pyrolite melt for all cations except for iron, which has a lower coordination number at all pressures, likely due to the difference in spin. Our simulations are spin polarized (section 2 for more details), whereas the MORB simulations are non-spin polarized, resulting in a smaller ionic radius and consequently, smaller coordination. It has been shown in previous calculations on pyrolite that iron remains in a magnetized state for most of the pressures and temperatures of the magma ocean (Caracas et al., 2019).

We also compare our results to experiments on silicate glasses. The coordination environment of cations with respect to oxygen has been determined experimentally using a wide range of techniques, including infrared spectroscopy, nuclear magnetic resonance spectroscopy (NMR), small-angle X-ray scattering, and inelastic X-ray Raman scattering.

**4.1. Coordination Environments**

**4.1.1. Silicon-Oxygen**

In crystalline materials, silicon is generally fourfold coordinated by oxygen at upper mantle pressures (e.g., forsterite and enstatite) (Smith, 1959; Smyth & Hazen, 1973) and increases to sixfold in the lower mantle
(e.g., periclase and bridgmanite) (Hazen, 1976; Ito & Matsui, 1978). Unlike crystalline materials, melts have a distribution of coordination environments, which smoothly increase with increasing pressure. Silicon exists in mostly fourfold coordination at ambient pressure and temperatures of 2,000–4,000 K in pyrolite melt (Figure 3a), in agreement with ab initio simulations on MORB (Bajgain et al., 2015; Karki et al., 2018), Mg₂SiO₄ (de Koker et al., 2008), CaMgSi₂O₆ (Sun et al., 2011), and MgSiO₃ melts (Stixrude & Karki, 2005). At 2,000 K and ambient pressure, 1% of silicon in pyrolite melt exists in threefold coordination, 1% exists in fivefold coordination, and the remaining 98% exists in fourfold coordination. At higher temperatures, the amount of threefold-coordinated silicon increases, which has never been experimentally investigated, mostly due to the fact that experiments are limited in temperature. Although Karki et al. (2018) do not explicitly report threefold silicon in MORB melt, at 3,000 K and ambient pressure, it is reported that 91% of silicon is in fourfold coordination and 5% of silicon is in fivefold coordination with an average coordination of 4, indicating that 4% of silicon is in threefold coordination. At equivalent conditions, we predict that 5% of silicon is threefold, 93% is fourfold, and 1% is fivefold, in agreement with the results of Karki et al. (2018) despite the different melt chemistries. The nontrivial amount of threefold-coordinated silicon at >3,000 K suggests that it is more relevant than previously thought, particularly in the Early Earth's magma ocean where the surface temperatures were much above the liquidus, especially if an atmosphere was forming.

With compression at upper mantle conditions, the coordination of silicon in pyrolite melt increases at a slower pace than of other cations, remaining close to fourfold over a large pressure range due to the deep potential minimum of silica tetrahedra. We find that with further increasing depth, its average coordination number increases following a general trend from fourfold at shallow depths (<10 GPa and 2,000 K) to fivefold at midmantle depths (40–50 GPa and 3,000 K) to sixfold at core-mantle boundary conditions (135 GPa and 4,000 K). At core-mantle boundary pressures, the proportion of sixfold silicon decreases with increasing temperature from 80% at 3,000 K to 70% at 4,000 K to 55% at 5,000 K. The proportions of silicon in fourfold, fivefold, and sixfold coordination in pyrolite melt are essentially identical to the proportions in CaMgSi₂O₆ melt at all pressures (Figure S3). Compared to simulations on MORB melt (Bajgain et al., 2015; Karki et al., 2018), the proportions of fourfold and fivefold silicon are slightly larger in pyrolite melt and CaMgSi₂O₆ melt (Sun et al., 2011) while the proportion of sixfold silicon is a bit lower because of the different degrees of polymerization or M: Si ratios (where M is Mg, Ca, Fe, and other network-modifier cations), which cumulatively explains the larger average coordination of silicon in MORB melt. The average coordination of silicon at high pressure is likely dependent on the degree of polymerization of the melt. At ambient pressure, the ratio, non-bridging oxygen per tetrahedrally coordinated cation (NBO/T) is often used to quantify the degree of polymerization such that olivine is completely depolymerized with an NBO/T ratio of 4, pyrolite is moderately depolymerized with an NBO/T ratio of 2.6, MgSiO₃ has an NBO/T ratio of 2, MORB is relatively polymerized with an NBO/T ratio of about 0.7, and SiO₂ and feldspars are fully polymerized with an NBO/T ratio of 0. Although at high pressure the NBO/T ratio is no longer valid, the relative polymerization of the melts at ambient pressure gives insight to their relative behavior at high pressure.

Using the first peak of the radial distribution function for MgSiO₃ glass, which corresponds to the combined influence of the coordination environments of silicon and magnesium with respect to oxygen, Kono et al. (2018) suggest that the coordination of silicon increases from fourfold to sixfold to about 70 GPa, after which, it remains constant at sixfold up to about 110 GPa. In our simulations on pyrolite melt, the coordination of silicon increases to 5.6 by 75 GPa and up to 6 by 140 GPa at 3,000 K without a plateau, although the density of data is not high enough to eliminate the possibility of a plateau between data points. Sanloup, Drewitt, Crépisson, et al. (2013) conducted X-ray diffraction experiments on basalt melt up to 60 GPa, finding that the coordination of silicon increases to sixfold by 35 GPa and plateauing at an average coordination number of 6.2 up to 60 GPa. In contrast, silicon is still fivefold coordinated in our simulations on pyrolite melt at 35 GPa and 2,000–3,000 K, likely due to the lower degree of polymerization in pyrolite melt.

Using mid-infrared absorption spectroscopy on fully polymerized glasses (silica [SiO₂], anorthite [CaAl₂Si₂O₈], and diopside [CaMgSi₂O₆]), Williams and Jeanloz (1988) observed a gradual increase in the coordination of silicon from entirely fourfold at 0 GPa to increasingly sixfold at higher pressures. At 40 GPa, the highest pressure achieved, the silica and diopside glasses contained coordination environments ranging from fourfold to sixfold, whereas the anorthite glass appeared to contain no fourfold silicon by 40 GPa. With inelastic X-ray Raman scattering, Lin et al. (2007) demonstrated that silicon in SiO₂ glass
undergoes a coordination change from fourfold like to sixfold like between 10 and 22 GPa; however, the spectra for the glass are much broader than the spectra for the crystalline quartz and stishovite standards, which is explained by the possible presence of fivefold coordination. The combined effect of broadening and increasing of the coordination environments may also yield the same effect. Compared to these results on silica glass, the fourfold to sixfold transition in silicon in pyrolite melt occurs at much higher pressures.

Several studies have calculated the average coordination of silicon with respect to oxygen in SiO₂ glass using X-ray diffraction, generally finding that the coordination number of silicon increases from fourfold to sixfold by about 40–50 GPa, followed by a more gradual increase in coordination (Maede et al., 1992; Sato & Funamori, 2010; Brazhkin et al., 2011; Prescher et al., 2017); however, slight differences in the trends were observed. Using average Si-O distances in SiO₂ glass measured with X-ray diffraction and comparing to crystalline materials, Meade et al. (1992) found that silicon exists almost entirely in fourfold coordination between 0 and 10 GPa, increasing rapidly towards sixfold between 10 and 40 GPa. Prescher et al. (2017) also find that the coordination of silicon in SiO₂ plateaus at fourfold between 0 to 14 GPa, after which, it increases steeply to sixfold by about 50 GPa. Similarly, in our simulated pyrolite melt at 2,000 K, we observe a plateau at fourfold coordination between 0 and 10 GPa, after which, it rises sharply to fivefold by 40 GPa, the highest pressure achieved. In contrast, Brazhkin et al. (2011) observe a steep continuous increase in coordination from 4 to 5.7 between 0 and 40 GPa without a plateau in both the compression and decompression regimes. Above 50 GPa, Brazhkin et al. and Prescher et al. observe a shallow increase in coordination up to about 6.75 by core-mantle boundary pressures while Sato and Funamori (2010) observe a plateau in coordination at sixfold up to about 100 GPa. Although in our simulations, we also observe a change of slope at around 50 GPa, the overall coordination numbers at 3,000–5,000 K are lower, reaching a maximum of sixfold coordination by core-mantle boundary pressures. The difference is likely due to the different chemistries (e.g., the presence of network-modifying cations in pyrolite melt and thus much lower degree of polymerization) and to a lesser extent, the higher temperatures in our simulations.

While the pressure-induced coordination change of silicon by oxygen from fourfold to sixfold is well documented in both silicate melts and silicate minerals (Liu & Bassett, 1986), fivefold silicon rarely occurs in crystalline phases and has consequently received less attention in silicate melts. Fivefold silicon is rare even in high-pressure crystalline phases, with the most notable example in CaSi₂O₅ (Kanzaki et al., 1991; Angel, 1997); instead, silicon in minerals typically shows either fourfold or sixfold coordination or a mixture of the two (Finger & Hazen, 2000). In minerals, the pressure-induced progression from fourfold to sixfold is either direct or occurs via a mixed-coordination phase. For example, silicon was observed to transition from entirely fourfold to mixed fourfold and sixfold to entirely sixfold coordination in forsterite as it transforms to metastable high-pressure polymorphs (Finkelstein et al., 2014). Unlike crystalline materials, melts and glasses have a continuous distribution of coordination environments, which gradually increase in coordination with increasing pressure, and thus include coordinations that are rarely observed in minerals, such as fivefold and sevenfold coordinations.

Even at ambient pressure, a small proportion of fivefold-coordinated silicon was experimentally observed in various alkali silicate glasses, with the proportion of fivefold-coordinated silicon increasing with increasing pressure and temperature. For example, NMR experiments on K₂Si₄O₉ and Na₂Si₄O₉ glasses found that 0–0.16% of silicon was fivefold in glasses quenched from ambient pressure and 3–4% of silicon was found to be fivefold in glasses quenched from 6 GPa (Stebbins & McMillan, 1993) while NMR experiments on sodium aluminosilicate glasses observe that 0.3–3% of silicon is fivefold in ambient-pressure glasses and 0–3% is sixfold in glasses synthesized at 6 GPa (Kelsey et al., 2009). Consistent with the experiments, we predict that at 2,000 K, about 1% of Si is fivefold at 0 GPa and 4% is fivefold at 6 GPa (from extrapolation between the data points at 3 and 9 GPa). Our simulations show that the proportion of fivefold-coordinated silicon continues to increase with a maximum proportion of fivefold coordination at about 40–50 GPa, accounting for nearly 50% of SiOₓ species at all temperatures (Figure 5a).

### 4.1.2. Aluminum-Oxygen

Aluminum typically substitutes for tetrahedrally coordinated silicon in silicate crystals (e.g., in feldspars) but may behave as both a network former and a network modifier in silicate melts, even at ambient pressure. We find that at ambient pressure, aluminum displays a similar coordination environment to silicon, existing predominantly in fourfold coordination at ambient pressure, but displays a slightly broader distribution
environment. For example, at 2,000 K and ambient pressure, 5%, 87%, and 8% of aluminum exists in three-fold, fourfold, and fivefold coordinations, respectively, compared to 1%, 98%, and 1% for silicon, respectively. At higher temperatures, the amount of threefold and fivefold aluminum increases to 10% and 23%, respectively. Although Karki et al. (2018) do not report threefold aluminum in MORB melt, at 3,000 K and ambient pressure, it is reported that 66% of aluminum is in fourfold coordination and 24% of aluminum is in fivefold coordination with an average coordination of 4.2, indicating that 8% of aluminum is in threefold coordination, in agreement with our predictions despite the different chemistries.

Compared to silicon, the coordination of aluminum increases much more rapidly with depth. By 500 km, the average coordination is fivefold and by 1,200 km, it is sixfold with no fourfold aluminum remaining (in contrast to silicon, which exists in fourfold coordination even at lowermost mantle pressures, albeit in small proportions). A similar trend was observed in ab initio simulations of MORB melt (Bajgain et al., 2015; Karki et al., 2018). Compared to pyrolite melt, the proportions of sixfold and sevenfold aluminum in MORB melt increase slightly more rapidly with pressure at the expense of fourfold and fivefold aluminum, resulting in an overall larger average coordination, likely as a result of a larger alumina concentration compared to pyrolite (16 wt % compared to 4.7 wt % Al₂O₃).

Experiments on aluminosilicate and tectosilicate glasses have demonstrated that the coordination of aluminum is mostly fourfold with the proportion of fivefold-coordinated aluminum ranging from 0% to 30%, depending strongly on composition, but generally increasing with increasing alumina content (Risbud et al., 1987; Sato et al., 1991; Poe et al., 1992; Toplis et al., 2000; Sen & Youngman, 2004; Neuville, Cormier, Flank, et al., 2004; Neuville, Cormier, Massiot, 2004; Neuville et al., 2006; Thompson & Stebbins, 2011; Malfait et al., 2012), increasing non-bridging oxygen (NBO) content (Lee, 2004; Lee et al., 2004; Lee, 2010; Park & Lee, 2018) and increasing field strength of the modifying cations (Allwardt et al., 2005; Kelsey et al., 2009; Morin et al., 2014). Kelsey et al. (2009) conducted a detailed NMR study on a series of sodium aluminosilicate glasses synthesized at 1 atm and 6 GPa. In all the glasses synthesized at 1 atm, a small fraction of fivefold-coordinated aluminum was present (0.8–2.2%) without any sixfold-coordinated aluminum whereas in the glasses synthesized at 6 GPa, both fivefold and sixfold aluminum were observed. The average coordination of aluminum was found to be fourfold for all glasses at 0 GPa but ranged from 4.6 to 5.5 at 6 GPa. We predict that at 2,000 K, aluminum's average coordination is 4 at 0 GPa and increases to about 4.4 by 6 GPa where 55% is fourfold, 35% is fivefold, and 5% is sixfold. In comparison, Kelsey et al. (2009) observed proportions of 58% fourfold, 28% fivefold, and 14% sixfold in the (Al/Al+Si)=0.2 glass with an average coordination of 4.6, and Bajgain et al. (2005) predicted roughly 45% fourfold, 35% fivefold, and 15% sixfold in MORB melt at 6 GPa from extrapolation.

Several studies have observed a minimum in the magnitude of the viscosity and/or plateauing of the viscosity between 0.5 and 6 GPa (Kushiro 1986; Poe et al., 1997; Suzuki et al., 2005; Allwardt et al., 2007). It has been hypothesized that the viscosity minimum may correspond to the maximum abundance of fivefold aluminum (Poe et al., 1997; Poe & Rubie, 2000). However, we observe a maximum in fivefold aluminum at 10–20 GPa (where the maximum increases with increasing temperature; Figures S4 and S5). The viscosity minimum may instead be related to the stagnation of the formation of SiO₆ (Figure S3), which in our study remains at ~0% abundance between 0 and 9 GPa, after which, it increases steeply to 36% abundance by 36 GPa.

### 4.1.3. Magnesium-Oxygen

Despite the different melt compositions, our prediction that the average coordination of magnesium by oxygen increases from about 4.5 at ambient pressure to nearly 8 by 135 GPa is in agreement with ab initio simulations on MORB melt (Bajgain et al., 2015), Mg₂SiO₄ melt (de Koker et al., 2008), and CaMgSi₂O₆ melt (Sun et al., 2011; Figure 3c). In fact, at 3,000 K and pressures of about 40–120 GPa, the average coordination of magnesium in pyrolite melt generally falls between the average coordination of magnesium in Mg₂SiO₄ and CaMgSi₂O₆ melts such that there is a general decrease in magnesium's average coordination with increasing magnesium content at approximately lower mantle pressures.

X-ray scattering, NMR, and Mg K-edge X-ray absorption near edge structure (XANES) experiments have found that the average coordination of magnesium in ambient-pressure silicate glasses varies from fourfold to sixfold, depending on the chemistry (Shimoda et al., 2008; Trcera et al., 2009). The average coordination of magnesium was found to be 4.1–4.5 in MgSiO₃ enstatite glass (Waseda & Toguri, 1977; Yin et al., 1983), 4.8 in Mg₃Al₂Si₃O₁₂ pyrope glass (Okuno & Marumo, 1993), and 5 in diopside and albite glasses (Ildefonse et al.,
1995; Dien et al., 1999). As a magnesium-rich glass, the calculated average coordination of 4.6 in pyrolite melt is consistent with the experimental studies on the enstatite and pyrope glass systems.

In crystalline phases, magnesium has a much larger coordination number. At low pressure, magnesium has a coordination number of 6 in pyroxenes (Warren & Bragg, 1929; Smith, 1959) and 8 in garnets (Novak & Gibbs, 1971) while in bridgmanite of the lower mantle, magnesium exists in a mixture of octahedral and highly distorted dodecahedral sites (Dobson & Jacobsen, 2004; Jephcoat et al., 1999), yet we find that in mafic silicate melts, the average coordination of magnesium remains below eightfold in most of the lower mantle, indicating that the packing of silicate melts is less efficient than crystalline materials. Even though the coordination of silicon in melts is similar to silicon in crystalline phases and the coordination increases at a comparable rate, network-modifying cations are persistently surrounded by fewer oxygen atoms than in crystalline phases at equivalent pressure-temperature conditions.

### 4.1.4. Iron-Oxygen

In ambient-pressure minerals, ferrous iron is typically in sixfold coordination (Smyth & Hazen, 1973; Cameron et al., 1973) while in silicate melts and glasses, ferrous iron exists in fourfold, fivefold, and sixfold coordinations, depending on composition as well as interpretation of experimental results. Using radial distribution functions from high-temperature X-ray diffraction experiments, Waseda et al. (1980) calculated coordination numbers of iron ranging from 4.3 to 5.7 in FeO-Fe₂O₃-SiO₂ melts (increasing with decreasing silica content) and with Mössbauer spectroscopy and Fe K-edge XANES on a wide variety of ferrosilicate glasses, Jackson et al. (2005) found that the average coordination of ferrous iron is as low as 3.9 in sodium-iron silicate glass and as high as 5.2 in basaltic glass. With similar results, a detailed Mössbauer spectroscopy study on aluminosilicate glasses suggests that ferrous iron is mostly fourfold in sodium aluminosilicate glasses and fivefold to sixfold in calcium and magnesium aluminosilicate glasses (Mysen, 2006). Ferrous iron was also found to be in fourfold coordination in Na₂FeSi₃O₈ and K₂FeSi₃O₈ glasses (Waychunas et al., 1988), Fe₀.₄Mg₀.₈Ca₀.₈Si₂O₆ glass (Calas & Petiau, 1983), and Fe₂SiO₄ glass (Cooney & Sharma, 1990) with X-ray absorption spectroscopy; however, in a careful X-ray diffraction study on fayalite melts, Sanloup, Drewitt, Konôpková, et al. (2013) found that ferrous iron has an average coordination number of 4.8 at 0 GPa and 7.2 at 7.5 GPa.

With XANES and multiple-scattering simulations on submarine basaltic glass, Wu et al. (1999) observed iron in predominantly fourfold coordination with lesser amounts of fivefold-coordinated iron and no indication of sixfold coordination. Rossano et al. (2000) found that ferrous iron exists in fourfold and fivefold coordinations in CaFeSi₂O₆ glass using a combination of extended X-ray absorption fine structure and molecular dynamics simulations. Comparing Na-, K-, and Al-bearing melts to glasses with X-ray absorption spectroscopy, Wilke et al. (2007) found an average coordination close to 5, with a slightly higher concentration of fourfold coordinated iron in melts with respect to glasses. Iron was interpreted to be predominantly in sixfold coordination in albite-diopside glasses (Keppler, 1992), basalt glass (Bell & Mao, 1974), and a variety of Fe-Ti glasses (Nolet et al., 1979) with optical absorption spectroscopy. Using time-domain Mössbauer and optical absorption spectroscopy on reduced basaltic glass, Solomatova et al. (2017) found that less than 15–20% of ferrous iron is in fourfold coordination and the rest of iron is consistent with fivefold and/or sixfold coordination. Thus, although the proportions of the site environments depend on composition, the reported coordination environments also depend strongly on the experimental technique and method of analysis and interpretation.

Ab initio calculations on silicate melts have found that iron has an average coordination of about 3.5–3.7 at ambient pressure and 2,000–3,000 K (Bajgain et al., 2015; this study), below experimentally determined values. Two main possibilities may explain the discrepancy: the coordination of iron in glasses is lower than in melts (e.g., Wilke et al., 2007) or density functional theory fails to accurately predict the geometric environment of iron likely due to its failure to accurately treat the d electrons. Indeed, it is difficult to imagine that iron has a lower coordination compared to magnesium, especially considering that the difference is not trivial (a difference of one oxygen atom at ambient pressure in our simulations).

### 4.1.5. Calcium-Oxygen

XANES and X-ray diffraction experiments on calcium-bearing silicate glasses have found that at ambient pressure, calcium bonds to oxygen with a wide range of bond distances, ranging from 2.29 to 2.85 Å with an average coordination of about 7, ranging from 6 to 7.4, depending on the experimental technique and
composition of the glass (Taylor & Brown, 1979; Okuno & Marumo, 1993; Eckersley et al., 1988; Cormier et al., 2000; Cormier & Neuville, 2004, Neuville, Cormier, Massiot, 2004). However, the proportions of calcium coordination environments in amorphous silicates have not yet been determined experimentally. Our results demonstrate that calcium reaches an average coordination of 9–10 by roughly 40 GPa and stays roughly constant up to core-mantle boundary conditions, whereas calcium in majorite garnet is already twelvefold within the upper mantle (Hazen et al., 1994).

### 4.1.6. Sodium-Oxygen

At ambient pressure and temperature of 2,000 K, the coordination of sodium is 6 but increases very rapidly with pressure. By the base of the upper mantle, it reaches a maximum of 8–9. The coordination of sodium in crystalline low albite is fivefold or sevenfold (Phillips et al., 1988), bounding our determined coordination of sixfold in ambient-pressure pyrolite. In amorphous silicates, experimentally determining sodium-oxygen bond distances can be arbitrary as the first and second coordination spheres are continuous rather than distinct. NMR experiments on sodium-rich silicate glasses and melts have estimated that the average coordination of sodium is about 6–8 without a dependence on composition (Maekawa et al., 1997; Xue & Stebbins, 1993).

### 4.1.7. Partitioning of Minor Elements Between the Molten Iron and Silicate Melt

The structure of silicate melts likely affects the partitioning behavior of minor elements between silicate melt and iron metal (Hillgren et al., 1996; Keppler & Rubie, 1993; Mann et al., 2012). The core-to-mantle partition coefficient of highly and moderately siderophile elements decreases with increasing pressure more steeply between 0 and ~6 GPa and more shallowly above ~6 GPa (Kegler et al., 2008; Mann et al., 2012). In pyrolite melt, we observe several sharp changes in the atomic coordination environments within the first 10 GPa. Up to about 10 GPa, the concentration of sixfold-coordinated silicon remains below 1%, after which, its concentration increases steeply. At 2,000 K, the concentration of sixfold-coordinated silicon reaches 4% at 20 GPa and 36% at 36 GPa while at 3,000 K, the concentration of sixfold-coordinated silicon is 7% at 23 GPa and 41% at 27 GPa (Figure S3). Meanwhile, the concentrations of fivefold-coordinated magnesium and sevenfold-coordinated calcium reach a local maximum such that their concentrations increase with increasing pressure up to 3–5 GPa at 2,000–3,000 K, accounting for roughly 45% and 35% of their corresponding species, and decrease with increasing pressure up to 40 GPa, at which point, their relative concentrations reach 0% (Figures S6 and S9). These coordination environments may provide a favorable environment for siderophile elements within the silicate melt; and so, the steep increase in their abundance at low pressure may be responsible for the steep slope observed in their core-to-mantle partition coefficients below ~6 GPa. Furthermore, experiments on carbon-bearing iron melt have shown a compressibility increase at about 5 GPa (Sanloup et al., 2011). Thus, we believe that the combined structural changes in silicate and iron melt may be responsible for the kink in the partitioning behavior between the two melts (Keppler & Rubie, 1993; Mann et al., 2012).

## 5. Conclusions

Using ab initio molecular dynamics calculations, we performed a detailed analysis of the local coordination environments of the major cations in pyrolite melt, a good approximation of the bulk silicate Earth composition (McDonough & Sun, 1995). We considered the pressure and temperature conditions spanning the entire Earth's mantle as well as the conditions of the magma ocean that were prevalent during the Early Earth, comparing our results to previous experimental and computational studies on amorphous silicates.

We find a broad distribution of the cation-oxygen coordination environments. The sharpest distribution is found in silicon at uppermost mantle conditions, occupying predominantly fourfold tetrahedral coordination while the broadest distribution of site environments is found in calcium and sodium at upper mantle conditions, sharpening slightly at higher pressures as the coordination is effectively maximized. The average coordination numbers for the cations with respect to oxygen generally increase with increasing pressure more steeply in the upper mantle and more gradually in the lower mantle. We computed the lifetimes of all the species present in our simulations. At 2,000 K and ambient-pressure conditions, we find $SiO_4$ tetrahedra as the longest-living species, with an average close to 3 picoseconds while $AlO_4$ tetrahedra exist for an average of 0.3 picoseconds. In contrast, the lifetimes of most of the coordination polyhedra of the network-modifying cations, that is Fe, Mg, Ca, and Na, are on the order of a few tens of femtoseconds at
most. As temperature increases, the lifetimes of the species decrease, as more the atoms have more kinetic energy, and bonds are more easily broken and reformed.

We observe a non-negligible contribution from exotic coordination environments (relative to crystalline Earth materials), such as threefold-coordinated cations at ambient pressure and temperatures of 3,000 to 5,000 K, predominantly fivefold-coordinated aluminum at midmantle conditions, and significant amounts of sevenfold coordination of all cations at lowermost mantle conditions. In fact, this study is the first computational study to highlight the presence of threefold-coordinated cations, which likely played a significant role in the physical properties of the Early Earth's magma ocean. Sharp changes in coordination behavior of fivefold-coordinated magnesium and sevenfold-coordinated calcium at about 3–5 GPa and the sharp increase in sixfold-coordinated silicon above 10 GPa may explain a kink found in experiments in the pressure dependence of the partitioning coefficient of minor elements between the magma ocean and liquid iron. Our detailed structural results may also be used to improve geodynamic modeling of planetary magma oceans, as the density, viscosity, and compressibility behavior of melts are directly linked to the geometric packing of its ions.

## Acknowledgments

This research was supported by the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program (Grant Agreement n°681818-IMPACT to RC). We acknowledge access to the GENCI supercomputers (Occigen, Ada, and Curie) through the stl2816 series of eDARI computing grants. All speciation data presented in this manuscript are publicly available at https://doi.org/10.5281/zenodo.3527985.

## References

* Allwardt, J. R., Stebbins, J. F., Schmidt, B. C., Frost, D. J., Withers, A. C., & Hirschmann, M. M. (2005). Aluminum coordination and the densification of high-pressure aluminosilicate glasses. *American Mineralogist*, 90(7), 1218–1222. https://doi.org/10.2138/am.2005.1836

* Allwardt, J. R., Stebbins, J. F., Terasaki, H., Du, L.-S., Frost, D. J., Withers, A. C., et al. (2007). Effect of structural transitions on properties of high-pressure silicate melts: 27Al NMR, glass densities, and melt viscosities. *American Mineralogist*, 92(7), 1093–1104. https://doi.org/10.2138/am.2007.2530

* Angel, R. J. (1997). Transformation of fivefold-coordinated silicon to octahedral silicon in calcium silicate, CaSi₂O₅. *American Mineralogist: Journal of Earth and Planetary Materials*, 82(7–8), 836–839.

* Bajgain, S., Ghosh, D. B., & Karki, B. B. (2015). Structure and density of basaltic melts at mantle conditions from first-principles simulations. *Nature Communications*, 6, 8578. https://doi.org/10.1038/ncomms9578

* Bell, P. M., & Mao, H. K. (1974). Crystal-field spectra of Fe²⁺ and Fe³⁺ in synthetic basalt glass as a function of oxygen fugacity. *Carnegie Institute of Washington Year Book*, 73, 496–497.

* Blöchl, P. E. (1994). Projector augmented-wave method. *Physical Review B*, 50(24), 17953.

* Bower, D. J., Wicks, J. K., Gurnis, M., & Jackson, J. M. (2011). A geodynamic and mineral physics model of a solid-state ultralow-velocity zone. *Earth and Planetary Science Letters*, 303(3–4), 193–202. https://doi.org/10.1016/j.epsl.2010.12.035

* Brazhkin, V. V., Lyapin, A. G., & Trachenko, K. (2011). Atomistic modeling of multiple amorphous-amorphous transitions in SiO₂ and GeO₂ glasses at megabar pressures. *Physical Review B*, 83(13), 132103. https://doi.org/10.1103/PhysRevB.83.132103

* Calas, G., & Petiau, J. (1983). Coordination of iron in oxide glasses through high-resolution K-edge spectra: Information from the pre-edge. *Solid State Communications*, 48(7), 625–629.

* Cameron, M., Sueno, S., Prewitt, C. T., & Papike, J. J. (1973). High-temperature crystal chemistry of acmite, diopside, hedenbergite jadeite, spodumene and ureyite. *American Mineralogist: Journal of Earth and Planetary Materials*, 58(7–8), 594–618.

* Caracas, R., Hirose, K., Nomura, R., & Ballmer, M. D. (2019). Melt-crystal density crossover in a deep magma ocean. *Earth and Planetary Science Letters*, 516, 202–211. https://doi.org/10.1016/j.epsl.2019.03.031

* Cooney, T. F., & Sharma, S. K. (1990). Structure of glasses in the systems Mg₂SiO₄-Fe₂SiO₄, Mn₂SiO₄-Fe₂SiO₄, Mg₂SiO₄-CaMgSiO₄, and Mn₂SiO₄-CaMnSiO₄. *Journal of Non-Crystalline Solids*, 122(1), 10–32. https://doi.org/10.1016/0022-3093(90)90220-G

* Cormier, L., & Neuville, D. R. (2004). Ca and Na environments in Na₂O-CaO-Al₂O₃-SiO₂ glasses: Influence of cation mixing and cation-network interactions. *Chemical Geology*, 213(1–3), 103–113, https://doi.org/10.1016/j.chemgeo.2004.08.049

* Cormier, Laurent, Neuville, D. R., & Calas, G. (2000). Structure and properties of low-silica calcium aluminosilicate glasses. *Journal of Non-Crystalline Solids*, 274(1–3), 110–114. https://doi.org/10.1016/S0022-3093(00)00209-X

* Dien, L., Mingsheng, P., & Murata, T. (1999). Coordination and local structure of magnesium in silicate minerals and glasses; Mg K-edge XANES study. *The Canadian Mineralogist*, 37(1), 199–206.

* Dobson, D. P., & Jacobsen, S. D. (2004). The flux growth of magnesium silicate perovskite single crystals. *American Mineralogist*, 89(5–6), 807–811. https://doi.org/10.2138/am-2004-5-615

* Eckersley, M. C., Gaskell, P. H., Barnes, A. C., & Chieux, P. (1988). Structural ordering in a calcium silicate glass. *Nature*, 335(6190), 525.

* Finger, L. W., & Hazen, R. M. (2000). Systematics of high-pressure silicate structures. *Reviews in Mineralogy and Geochemistry*, 41(1), 123–155.

* Finkelstein, G. J., Dera, P. K., Jahn, S., Oganov, A. R., Holl, C. M., Meng, Y., & Duffy, T. S. (2014). Phase transitions and equation of state of forsterite to 90 GPa from single-crystal X-ray diffraction and molecular modeling. *American Mineralogist*, 99(1), 35–43. https://doi.org/10.2138/am.2014.4526

* Fiquet, G., Auzende, A. L., Siebert, J., Corgne, A., Bureau, H., Ozawa, H., & Garbarino, G. (2010). Melting of peridotite to 140 gigapascals. *Science*, 329(5998), 1516–1518. https://doi.org/https://doi.org/10.1126/science.1192448

* Freitas, D., Manthilake, G., Schiavi, F., Chantel, J., Bolfan-Casanova, N., Bouhifd, M. A., & Andrault, D. (2017). Experimental evidence supporting a global melt layer at the base of the Earth's upper mantle. *Nature Communications*, 8(1), 2186, https://doi.org/10.1038/s41467-017-02275-9

* Hammond, W. C., & Humphreys, E. D. (2000). Upper mantle seismic wave velocity: Effects of realistic partial melt geometries. *Journal of Geophysical Research Solid Earth*, 105(B5), 10975–10986. https://doi.org/10.1029/2000JB900042

* Hazen, R. M. (1976). Effects of temperature and pressure on the cell dimension and X-ray temperature factors of periclase. *American Mineralogist*, 61(3–4), 266–271.

* Hazen, R. M., Downs, R. T., Finger, L. W., Conrad, P. G., & Gasparik, T. (1994). Crystal chemistry of Ca-bearing majorite. *American Mineralogist: Journal of Earth and Planetary Materials*, 79(5–6), 581–584.
Hillgren, V. J., Drake, M. J., & Rubie, D. C. (1996). High pressure and high temperature metal-silicate partitioning of siderophile elements: The importance of silicate liquid composition. *Geochimica et Cosmochimica Acta*, **60**(12), 2257–2263. https://doi.org/10.1016/0016-7037(96)00079-8

Hirschmann, M. M. (2000). Mantle solidus: Experimental constraints and the effects of peridotite composition. *Geochemistry, Geophysics, Geosystems*, **1**(10). https://doi.org/10.1029/2000GC000070

Hoover, W. G. (1985). Canonical dynamics: Equilibrium phase-space distributions. *Physical Review A*, **31**(3), 1695.

Ildefonse, P., Calas, G., Flank, A. M., & Lagarde, P. (1995). Low Z elements (Mg, Al, and Si) K-edge X-ray absorption spectroscopy in minerals and disordered systems. *Nuclear Instruments and Methods in Physics Research Section B: Beam Interactions with Materials and Atoms*, **97**(1–4), 172–175.

Ito, E., & Matsui, Y. (1978). Synthesis and crystal-chemical characterization of MgSiO3 perovskite. *Earth and Planetary Science Letters*, **38**(2), 443–450.

Jackson, W. E., Farges, F., Yeager, M., Mabrouk, P. A., Rossano, S., Waychunas, G. A., et al. (2005). Multi-spectroscopic study of Fe (II) in silicate glasses: Implications for the coordination environment of Fe (II) in silicate melts. *Geochimica et Cosmochimica Acta*, **69**(17), 4315–4332. https://doi.org/10.1016/j.gca.2005.01.008

Jephcoat, A. P., Hriljac, J. A., Cammon, C. A., O'Neill, H. S. C., Rubie, D. C., & Finger, L. W. (1999). High-resolution synchrotron X-ray powder diffraction and Rietveld structure refinement of two (Mg<sub>0.95</sub>,Fe<sub>0.05</sub>)SiO<sub>3</sub> perovskite samples synthesized under different oxygen fugacity conditions. *American Mineralogist*, **84**(3), 214–220.

Kanzaki, M., Stebbins, J. F., & Xue, X. (1991). Characterization of quenched high pressure phases in CaSiO<sub>3</sub> system by XRD and <sup>29</sup>Si NMR. *Geophysical Research Letters*, **18**(3), 463–466.

Karki, B. B., Ghosh, D. B., & Bajgain, S. K. (2018). Simulation of silicate melts under pressure. *Magmas Under Pressure*, 419–453. https://doi.org/10.1016/b978-0-12-811301-1.00016-2

Kegler, P., Holzheid, A., Frost, D. J., Rubie, D. C., Dohmen, R., & Palme, H. (2008). New Ni and Co metal-silicate partitioning data and their relevance for an early terrestrial magma ocean. *Earth and Planetary Science Letters*, **268**(1), 28–40. https://doi.org/10.1016/j.epsl.2007.12.020

Kelsey, K. E., Stebbins, J. F., Mosenfelder, J. L., & Asimow, P. D. (2009). Simultaneous aluminum, silicon, and sodium coordination changes in 6 GPa sodium aluminosilicate glasses. *American Mineralogist*, **94**(8–9), 1205–1215. https://doi.org/10.2138/am.2009.3177.

Keppler, H. (1992). Crystal field spectra and geochemistry of transition metal ions in silicate melts and glasses. *American Mineralogist*, **77**(1–2), 62–75.

Keppler, H., & Rubie, D. C. (1993). Pressure-induced coordination changes of transition-metal ions in silicate melts. *Nature*, **364**(6432), 54–56. https://doi.org/10.1038/364054a0

de Koker, N. P., Stixrude, L., & Karki, B. B. (2008). Thermodynamics, structure, dynamics, and freezing of Mg<sub>2</sub>SiO<sub>4</sub> liquid at high pressure. *Geochimica et Cosmochimica Acta*, **72**(5), 1427–1441. https://doi.org/10.1016/j.gca.2007.12.019

Kono, Y., Shibazaki, Y., Kenney-Benson, C., Wang, Y., & Shen, G. (2018). Pressure-induced structural change in MgSiO<sub>3</sub> glass at pressures near the Earth's core-mantle boundary. *Proceedings of the National Academy of Sciences*, **115**(8), 1742–1747. https://doi.org/10.1073/pnas.1716748115

Kresse, G., & Furthmüller, J. (1996). Software VASP, vienna (1999). *Phys. Rev. B*, **54**(11), 169.

Kushiro, I. (1986). *Viscosity of partial melts in the upper mantle*. *Journal of Geophysical Research*: Solid Earth, 9343–9350. https://doi.org/10.1029/JB091iB09p09343@10.1002/(ISSN)2169-9356.PARTMELT1

Lee, S. K. (2010). Effect of pressure on structure of oxide glasses at high pressure: Insights from solid-state NMR of quadrupolar nuclides. *Solid State Nuclear Magnetic Resonance*, **38**(2-3), 45–57. https://doi.org/10.1016/j.ssnmr.2010.10.002

Lee, S. K. (2004). Structure of silicate glasses and melts at high pressure: Quantum chemical calculations and solid-state NMR. *The Journal of Physical Chemistry B*, **108**(19), 5889–5900. https://doi.org/10.1021/jp037575d

Lee, S. K., Cody, G. D., Fei, Y., & Mysen, B. O. (2004). Nature of polymerization and properties of silicate melts and glasses at high pressure. *Geochimica et Cosmochimica Acta*, **68**(20), 4189–4200. https://doi.org/10.1016/j.gca.2004.04.002

Lin, J.-F., Fukui, H., Prendergast, D., Okuchi, T., Cai, Y. Q., Hiraoka, N., et al. (2007). Electronic bonding transition in compressed SiO<sub>2</sub> glass. *Physical Review B*, **75**(1), 012201. https://doi.org/10.1103/PhysRevB.75.012201

Litasov, K., & Ohtani, E. (2002). Phase relations and melt compositions in CMAS-pyrrolite–H<sub>2</sub>O system up to 25 GPa. *Physics of the Earth and Planetary Interiors*, **134**(1–2), 105–127. https://doi.org/10.1016/S0031-9201(02)00152-8

Maekawa, H., Nakao, T., Shimokawa, S., & Yokokawa, T. (1997). Coordination of sodium ions in NaAlO<sub>2</sub>–SiO<sub>2</sub> melts: A high temperature <sup>23</sup>Na NMR study. *Physics and Chemistry of Minerals*, **24**(1), 53–65, https://doi.org/10.1007/s002690050

Malfait, W. J., Verel, R., Ardia, P., & Sanchez-Valle, C. (2012). Aluminum coordination in rhyolite and andesite glasses and melts: Effect of temperature, pressure, composition and water content. *Geochimica et Cosmochimica Acta*, **77**, 11–26. https://doi.org/10.1016/j.gca.2011.11.011

Mann, U., Frost, D. J., Rubie, D. C., Becker, H., & Audôtat, A. (2012). Partitioning of Ru, Rh, Pd, Re, Ir and Pt between liquid metal and silicate at high pressures and high temperatures-Implications for the origin of highly siderophile element concentrations in the Earth's mantle. *Geochimica et Cosmochimica Acta*, **84**, 593–613. https://doi.org/10.1016/j.gca.2012.01.026

McDonough, W. F., & Sun, S.-S. (1995). The composition of the Earth. *Chemical Geology*, **120**(3-4), 223–253. https://doi.org/10.1016/0009-2541(94)00140-4

Meade, C., Hemley, R. J., & Mao, H. K. (1992). High-pressure X-ray diffraction of SiO<sub>2</sub> glass. *Physical Review Letters*, **69**(9), 1387. https://doi.org/10.1103/PhysRevLett.69.1387, DOI: https://doi.org/10.1103/PhysRevLett.69.1387

Morin, E. I., Wu, J., & Stebbins, J. F. (2014). Modifier cation (Ba, Ca, La, Y) field strength effects on aluminum and boron coordination in aluminoborosilicate glasses: The roles of fictive temperature and boron content. *Applied Physics A*, **116**(2), 479–490. https://doi.org/10.1007/s00339-014-8369-4

Mysen, B. O. (2006). The structural behavior of ferric and ferrous iron in aluminosilicate glass near meta-aluminosilicate joins. *Geochimica et Cosmochimica Acta*, **70**(9), 2337–2353. https://doi.org/10.1016/j.gca.2006.01.026

Neuville, D. R., Cormier, L., Flank, A.-M., Briois, V., & Massiot, D. (2004). Al speciation and Ca environment in calcium aluminosilicate glasses and crystals by Al and Ca K-edge X-ray absorption spectroscopy. *Chemical Geology*, **213**(1–3), 153–163. https://doi.org/10.1016/j.chemgeo.2004.08.039

Neuville, D. R., Cormier, L., & Massiot, D. (2004). Al environment in tectosilicate and peraluminous glasses: A <sup>27</sup>Al MQ-MAS NMR, Raman, and XANES investigation. *Geochimica et Cosmochimica Acta*, **68**(24), 5071–5079. https://doi.org/10.1016/j.gca.2004.05.048
Neuville, D. R., Cormier, L., & Massiot, D. (2006). Al coordination and speciation in calcium aluminosilicate glasses: Effects of composition determined by ²⁷Al MQ-MAS NMR and Raman spectroscopy. *Chemical Geology*, **229**(1–3), 173–185. https://doi.org/10.1016/j.chemgeo.2006.01.019

Nolet, D. A., Burns, R. G., Flamm, S. L., & Besancon, J. R. (1979). Spectra of Fe-Ti silicate glasses-Implications to remote-sensing of planetary surfaces. In *Lunar and Planetary Science Conference Proceedings* (Vol. 10, pp. 1775-1786).

Nomura, R., Hirose, K., Uesugi, K., Ohishi, Y., Tsuchiyama, A., Miyake, A., & Ueno, Y. (2014). Low core-mantle boundary temperature inferred from the solidus of pyrolite. *Science*, **343**(6170), 522–525. https://doi.org/10.1126/science.1248186

Nosé, S. (1984). A unified formulation of the constant temperature molecular dynamics methods. *The Journal of Chemical Physics*, **81**(1), 511–519.

Novak, G. A., & Gibbs, G. V. (1971). The crystal chemistry of the silicate garnets. *American Mineralogist: Journal of Earth and Planetary Materials*, **56**(5–6), 791–825.

Okuno, M., & Marumo, F. (1993). The structure analyses of pyrope (Mg₃Al₂Si₃O₁₂) and grossular (Ca₃Al₂Si₃O₁₂) glasses by X-ray diffraction method. *Mineralogical Journal*, **16**(8), 407–415.

Park, S. Y., & Lee, S. K. (2018). Probing the structure of Fe-free model basaltic glasses: A view from a solid-state ²⁷Al and ¹⁷O NMR study of Na-Mg silicate glasses, Na₂O-MgO-Al₂O₃-SiO₂ glasses, and synthetic Fe-free KLB-1 basaltic glasses. *Geochimica et Cosmochimica Acta*, **238**, 563–579. https://doi.org/10.1016/j.gca.2018.07.032

Perdew, J. P., Burke, K., & Ernzerhof, M. (1996). Generalized gradient approximation made simple. *Physical Review Letters*, **77**(18), 3865. https://doi.org/10.1103/PhysRevLett.77.3865, 3868, DOI: https://doi.org/10.1103/PhysRevLett.77.3865

Phillips, B. L., Kirkpatrick, R. J., & Hovis, G. L. (1988). ²⁷Al, ²⁹Si, and ²³Na MAS NMR study of an Al, Si ordered alkali feldspar solid solution series. *Physics and Chemistry of Minerals*, **16**(3), 262–275.

Poe, B. T., McMillan, P. F., Angell, C. A., & Sato, R. K. (1992). Al and Si coordination in SiO₂-Al₂O₃ glasses and liquids: A study by NMR and IR spectroscopy and MD simulations. *Chemical Geology*, **96**(3–4), 333–349. https://doi.org/10.1016/0009-2541(92)90063-B

Poe, B. T., McMillan, P. F., Rubie, D. C., Chakraborty, S., Yarger, J., & Diefenbacher, J. (1997). Silicon and oxygen self-diffusivities in silicate liquids measured to 15 Gigapascals and 2800 Kelvin. *Science*, **276**(5316), 1245–1248. https://doi.org/10.1126/science.276.5316.1245

Poe, B. T., & Rubie, D. C. (2000). Transport properties of silicate melts at high pressure. *Physics Meets Mineralogy*, 340–353.

Prescher, C., Prakapenka, V. B., Stefanski, J., Jahn, S., Skinner, L. B., & Wang, Y. (2017). Beyond sixfold coordinated Si in SiO₂ glass at ultrahigh pressures. *Proceedings of the National Academy of Sciences of the United States of America*, **114**(38), 10041–10046. https://doi.org/10.1073/pnas.1708882114

Revenaugh, J., & Sipkin, S. A. (1994). Seismic evidence for silicate melt atop the 410-km mantle discontinuity. *Nature*, **369**(6480), 474–476. https://doi.org/10.1038/369474a0

Risbud, S. H., Kirkpatrick, R. J., Taglialavore, A. P., & Montez, B. (1987). Solid-state NMR Evidence of 4-, 5-, and 6-fold aluminum sites in roller-quenched SiO₂-Al₂O₃ glasses. *Journal of the American Ceramic Society*, **70**(1), C-10-C-12. https://doi.org/10.1111/j.1151-2916.1987.tb04859.x

Rossano, S., Ramos, A. Y., & Delaye, J.-M. (2000). Environment of ferrous iron in CaFeSi₂O₆ glass; Contributions of EXAFS and molecular dynamics. *Journal of Non-Crystalline Solids*, **273**(1–3), 48–52. https://doi.org/10.1016/S0022-3093(00)00124-1

Sanloup, C., Drewitt, J. W. E., Crépisson, C., Kono, Y., Park, C., Cammon, C., et al. (2013). Structure and density of molten fayalite at high pressure. *Geochimica et Cosmochimica Acta*, **118**, 118–128. https://doi.org/10.1016/j.gca.2013.05.012

Sanloup, C., van Westrenen, W., Dasgupta, R., Maynard-Casely, H., & Perrillat, J.-P. (2011). Compressibility change in iron-rich melt and implications for core formation models. *Earth and Planetary Science Letters*, **306**(1), 118–122. https://doi.org/10.1016/j.epsl.2011.03.039

Sanloup, C., Drewitt, J. W. E., Konópková, Z., Dalladay-Simpson, P., Morton, D. M., Rai, N., et al. (2013). Structural change in molten basalt at deep mantle conditions. *Nature*, **503**(7474), 104–107. https://doi.org/10.1038/nature12668

Sato, R. K., McMillan, P. F., Dennison, P., & Dupree, R. (1991). High-resolution ²⁷Al and ²⁹Si MAS NMR investigation of SiO₂-Al₂O₃ glasses. *Journal of Physical Chemistry*, **95**(11), 4483–4489.

Sato, T., & Funamori, N. (2010). High-pressure structural transformation of SiO₂ glass up to 100 GPa. *Physical Review B*, **82**(18), 184102. https://doi.org/10.1103/PhysRevB.82.184102

Schmandt, B., Jacobsen, S. D., Becker, T. W., Liu, Z., & Dueker, K. G. (2014). Dehydration melting at the top of the lower mantle. *Science*, **344**(6189), 1265–1268. https://doi.org/https://doi.org/10.1126/science.1253358

Sen, S., & Youngman, R. E. (2004). High-resolution multinuclear NMR structural study of binary aluminosilicate and other related glasses. *The Journal of Physical Chemistry B*, **108**(23), 7557–7564. https://doi.org/10.1021/jp031348u

Shimoda, K., Nemoto, T., & Saito, K. (2008). Local structure of magnesium in silicate glasses: a ²⁵Mg 3QMAS NMR study. *The Journal of Physical Chemistry B*, **112**(22), 6747–6752. https://doi.org/10.1021/jp711417t

Smith, J. V. (1959). The crystal structure of proto-enstatite, MgSiO₃. *Acta Crystallographica*, **12**(7), 515–519, https://doi.org/10.1021/jp711417t

Smyth, J. R., & Hazen, R. M. (1973). The crystal structures of forsterite and hortonolite at several temperatures up to 900 C. *American Mineralogist: Journal of Earth and Planetary Materials*, **58**(7–8), 588–593.

Solomatova, N. V., Jackson, J. M., Sturhahn, W., Rossman, G. R., & Roskosz, M. (2017). Electronic environments of ferrous iron in rhyolitic and basaltic glasses at high pressure. *Journal of Geophysical Research: Solid Earth*, **122**(8), 6306–6322. https://doi.org/10.1002/2017JB014363

Song, T. R. A., Helmberger, D. V., & Grand, S. P. (2004). Low-velocity zone atop the 410-km seismic discontinuity in the northwestern United States. *Nature*, **427**(6974), 530.

Stebbins, J. F., & McMillan, P. (1993). Compositional and temperature effects on five-coordinated silicon in ambient pressure silicate glasses. *Journal of Non-Crystalline Solids*, **160**(1–2), 116–125. https://doi.org/10.1016/0022-3093(93)90292-6

Stixrude, L., & Karki, B. (2005). Structure and freezing of MgSiO₃ liquid in Earth's lower mantle. *Science*, **310**(5746), 297–299. https://doi.org/10.1126/science.1116952

Sun, N., Stixrude, L., de Koker, N., & Karki, B. B. (2011). First principles molecular dynamics simulations of diopside (CaMgSi₂O₆) liquid to high pressure. *Geochimica et Cosmochimica Acta*, **75**(13), 3792–3802. https://doi.org/10.1016/j.gca.2011.04.004

Suzuki, A., Ohtani, E., Terasaki, H., & Funakoshi, K. (2005). Viscosity of silicate melts in CaMgSi₂O₆-NaAlSi₂O₆ system at high pressure. *Physics and Chemistry of Minerals*, **32**(2), 140–145. https://doi.org/10.1007/s00269-005-0452-0

Tauzin, B., Debayle, E., & Wittlinger, G. (2010). Seismic evidence for a global low-velocity layer within the Earth's upper mantle. *Nature Geoscience*, **3**, 718. https://doi.org/10.1038/ngeo0969

Taylor, M., & Brown, G. E. Jr. (1979). Structure of mineral glasses—I. The feldspar glasses NaAlSi₃O₈, KAlSi₃O₈, CaAl₂Si₂O₈. *Geochimica et Cosmochimica Acta*, **43**(1), 61–75.
Thompson, L. M., & Stebbins, J. F. (2011). Non-bridging oxygen and high-coordinated aluminum in metaluminous and peraluminous calcium and potassium aluminosilicate glasses: High-resolution ¹⁷O and ²⁷Al MAS NMR results. *American Mineralogist*, *96*(5-6), 841-853. https://doi.org.10.2138/am.2011.3680

Tonks, W. B., & Melosh, H. J. (1993). Magma ocean formation due to giant impacts. *Journal of Geophysical Research Planets*, 98(E3), 5319-5333. https://doi.org.10.1029/92JE02726.

Toplis, M. J., Kohn, S. C., Smith, M. E., & Poplett, I. J. (2000). Fivefold-coordinated aluminum in tectosilicate glasses observed by triple quantum MAS NMR. *American Mineralogist*, 85(10), 1556-1560. https://doi.org.10.2138/am-2000-1031

Trcera, N., Cabaret, D., Rossano, S., Farges, F., Flank, A.-M., & Lagarde, P. (2009). Experimental and theoretical study of the structural environment of magnesium in minerals and silicate glasses using X-ray absorption near-edge structure. *Physics and Chemistry of Minerals*, 36(5), 241-257. https://doi.org.10.1007/s00269-008-0273-z

Vidale, J. E., & Hedlin, M. A. (1998). Evidence for partial melt at the core-mantle boundary north of Tonga from the strong scattering of seismic waves. *Nature*, 391(6668), 682.

Vinnik, L., & Farra, V. (2007). Low S velocity atop the 410-km discontinuity and mantle plumes. *Earth and Planetary Science Letters*, 262(3-4), 398-412. https://doi.org.10.1038/35601

Warren, B., & Bragg, W. L. (1929). XII. The structure of diopside, CaMg (SiO₃)₂. *Zeitschrift Für Kristallographie-Crystalline Materials*, 69(1-6), 168-193.

Waseda, Y., Shiraishi, Y., & Toguri, J. M. (1980). The structure of the molten FeO-Fe₂O₃-SiO₂ system by X-ray diffraction. *Transactions of the Japan Institute of Metals*, 21(1), 51-62. https://doi.org/10.2320/matertrans1960.21.51

Waseda, Y., & Toguri, J. M. (1977). The structure of molten binary silicate systems CaO-SiO₂ and MgO-SiO₂. *Metallurgical Transactions B*, 8(4), 563-568.

Waychunas, G. A., Brown, G. E., Ponader, C. W., & Jackson, W. E. (1988). Evidence from X-ray absorption for network-forming Fe²⁺ in molten alkali silicates. *Nature*, 332(6161), 251.

Wilke, M., Farges, F., Partzsch, G. M., Schmidt, C., & Behrens, H. (2007). Speciation of Fe in silicate glasses and melts by in-situ XANES spectroscopy. *American Mineralogist*, 92(1), 44-56. https://doi.org.10.2138/am.2007.1976

Williams, Q., & Garnero, E. J. (1996). Seismic evidence for partial melt at the base of Earth's mantle. *Science*, 273(5281), 1528-1530. https://doi.org.10.1126/science.273.5281.1528

Williams, Q., & Jeanloz, R. (1988). Spectroscopic evidence for pressure-induced coordination changes in silicate glasses and melts. *Science*, 239(4842), 902-905. https://doi.org/10.1126/science.239.4842.902

Wu, Z., Bonnin-Mosbah, M., Duraud, J. P., Métrich, N., & Delaney, J. S. (1999). XANES studies of Fe-bearing glasses. *Journal of Synchrotron Radiation*, 6(3), 344-346, https://doi.org.10.1107/S0909049598015969

Xue, X., & Stebbins, J. F. (1993). ²³Na NMR chemical shifts and local Na coordination environments in silicate crystals, melts and glasses. *Physics and Chemistry of Minerals*, 20(5), 297-307, https://doi.org.10.1007/BF00215100

Yin, C. D., Okuno, M., Morikawa, H., & Marumo, F. (1983). Structure analysis of MgSiO₃ glass. *Journal of Non-Crystalline Solids*, 55(1), 131-141. https://doi.org.10.1016/0022-3093(83)90013-3

Zerr, A., Diegeler, A., & Boehler, R. (1998). Solidus of Earth's deep mantle. *Science*, 281(5374), 243-246. https://doi.org.10.1126/science.281.5374.243