Synthesis and crystal-chemical characterization of MgSiO3 perovskite

The orthorhombic MgSiO₃ perovskite has been synthesized with the aid of a double-stage split-sphere-type high-pressure apparatus at about 280 kbar and 1000°C. The unit cell dimensions are: a = 4.7754(3)Å, b = 4.9292(4)Å and c = 6.8969(5)Å with the probable space group Pbnm. Calculated density is 4.108 g cm^?3. Crystal structure determination has been carried out by means of both the geometrical simulation (DLS) technique and the ordinary powder X-ray analysis. The results indicate that the MgSiO₃ perovskite is closer to the ideal perovskite than ScAlO₃ perovskite.     

Elasticity of MgSiO3 in the ilmenite phase

The single-crystal elastic moduli of the ilmenite phase of MgSiO₃ have been determined from Brillouin spectroscopy. They are: C₁₁ = 472, C₁₂ = 168, C₃₃ = 382, C₁₃ = 70, C₄₄ = 106, C₁₄ = ?27, C₆₆ = 152 and C₂₅ = ?24 in GPa. These elastic properties are consistent with a structural mechanical model where the silicon octahedra are very stiff under compression and shear. This latter property yields an unexpectedly high shear modulus for the magnesium silicate ilmenite as compared with analogue compounds. The further transformation to perovskite will probably be associated with a significant increase in elastic properties since the strong silicon polyhedra form a structural network in this phase. The transformation of spinel and stishovite to ilmenite is associated with a slight density increase and a slight decrease in acoustic velocities. This transformation will probably not produce a seismic discontinuity even if it does occur in the Earth's mantle.     

High-pressure phase transformations in the system of MgSiO3-FeSiO3

Two synthetic pyroxenes (FeSiO₃, MgSiO₃) and five natural pyroxenes with compositions of about Fs₈₀En₂₀, Fs₆₀En₄₀, Fs₅₀En₅₀, Fs₄₀En₆₀, and Fs₂₀En₈₀ have been subjected to pressures up to250 ± 50kbars at a temperature of about1500 ± 200°C in a diamond anvil cell heated by an infrared laser beam. After quenching and unloading X-ray data analysis indicates that (1) those with Mg less than 50% undergo the following reactions: 2(Mg,Fe)SiO₃ (pyroxene) → (Mg,Fe)₂SiO₄ (spinel) + SiO₂ (stishovite) → 2(Mg,Fe)O (magnesiowu¨stite) + SiO₂ (stishovite) with increase of pressure, and (2) those with Mg higher than 60%, undergo the following reactions: 2(Mg,Fe)SiO₃ (pyroxene) → (Mg,Fe)₂SiO₄ (spinel) + SiO₂ (stishovite) → 2(Mg,Fe)SiO₃ (hexagonal phase) → 2(Mg,Fe)O (magnesiowu¨stite) + SiO₂ (stishovite) with increase of pressure.     

The energetics of aluminum solubility into MgSiO3 perovskite at lower mantle conditions

Experiments [T. Irifune (1994) Nature 370, 131–133; E. Ito et al. (1998) Geophys. Res. Lett. 25, 821–824; A. Kubo, M. Akaogi (2000) Phys. Earth Planet. Int. 121, 85–102] indicate that (Mg,Fe)SiO₃ perovskite, commonly believed to be the most abundant mineral in the Earth, is the preferred host phase of Al₂O₃ in the Earth’s lower mantle. Aiming to better understand the effects of Al₂O₃ on the thermoelastic properties of the lower mantle, we use atomistic models to examine the chemistry and elasticity of solid solutions within the MgSiO₃(perovskite)–Al₂O₃(corundum)–MgO(periclase) mineral assemblage under conditions pertinent to the lower mantle: low Al cation concentrations, P=25–100 GPa, and T=1000–2000 K. We assess the relative stabilities of two likely substitution mechanisms of Al into MgSiO₃ perovskite in terms of reactions involving MgSiO₃, MgO, and Al₂O₃, in a manner similar to the 0 Kelvin calculations of Brodholt [J.P. Brodholt (2000) Nature 407, 620–622] and Yamamoto et al. [T. Yamamoto et al. (2003) Earth Planet. Sci. Lett. 206, 617–625]. We determine the equilibrium composition of the assemblage by examining the chemical potentials of the Al₂O₃ and MgO components in solid solution with MgSiO₃, as functions of concentration. We find that charge coupled substitution dominates at lower mantle pressures and temperatures. Oxygen vacancy-forming substitution accounts for 3–4% of Al substitution at shallow lower mantle conditions, and less than 1% in the deep mantle. For these two pressure regimes, the corresponding adiabatic bulk moduli of aluminous perovskite are 2% and 1% lower than that of pure MgSiO₃ perovskite.     

Viscosity and conductivity of the lower mantle; an experimental study on a MgSiO3 perovskite analogue,KZnF3

We have investigated the high-temperature rheological and electrical behaviour of single-crystal KZnF₃ fluoperovskite, an analogue of the MgSiO₃ perovskite in the lower mantle.The crystals flow by Newtonian dislocation creep (Harper-Dorn creep), predominantly on {100} planes. Below the melting point, solid-electrolyte behaviour appears, accompanied by a weakening of the mechanical properties. Geophysical implications are examined: the lower mantle can convect by Newtonian dislocation creep and an asthenosphere may exist at the bottom of the mantle. Electromagnetic interactions between the core and solid-electrolyte lower mantle may also be important.     

In situ measurements of sound velocities and densities across the orthopyroxene → high-pressure clinopyroxene transition in MgSiO3 at high pressure

Using acoustic measurement interfaced with a large volume multi-anvil apparatus in conjunction with in situ X-radiation techniques, we are able to measure the density and elastic wave velocities (V_P and V_S) for both ortho- and high-pressure clino-MgSiO₃ polymorphs in the same experimental run. The elastic bulk and shear moduli of the unquenchable high-pressure clinoenstatite phase were measured within its stability field for the first time. The measured density contrast associated with the phase transition OEN → HP-CEN is 2.6-2.9% in the pressure of 7-9 GPa, and the corresponding velocity jumps are 3-4% for P waves and 5-6% for S waves. The elastic moduli of the HP-CEN phase are K_S=156.7(8) GPa, G = 98.5(4) GPa and their pressure derivatives are K_S′=5.5(3) and G′ = 1.5(1) at a pressure of 6.5 GPa, room temperature. In addition, we observed anomalous elastic behavior in orthoenstatite at pressure above 9 GPa at room temperature. Both elastic wave velocities exhibited softening between 9 and 13-14 GPa, which we suggest is associated with a transition to a metastable phase intermediate between OEN and HP-CEN.     

Elasticity of single-crystal SmAlO3, GdAlO3 and ScAlO3 perovskites

The adiabatic single-crystal elastic moduli of SmAlO₃, GdAlO₃ and ScAlO₃, all with the orthorhombic perovskite structure, have been measured by Brillouin spectroscopy under ambient conditions. These 3 compounds display various degrees of crystallographic distortion from the ideal cubic perovskite structure. We find that longitudinal moduli in directions parallel to the axes of a pseudocubic subcell are nearly equal and insensitive to distortions of the crystal structure from cubic symmetry, whereas, the moduli C₁₁ and C₂₂, parallel to the orthorhombic axes, display pronounced anisotropy with the exception of ScAlO₃. The shear moduli also correlate with distortion from cubic symmetry, as measured by rotation, or tilt angles, of the AlO₆ octahedra. Our data support the observations of Liebermann et al. that perovskite-structure compounds define consistent elasticity trends relating bulk modulus and molar volume, and sound speed and mean atomic weight. These relationships have been used to estimate bulk and shear moduli for the high-pressure polymorphs of CaSiO₃ and MgSiO₃ with the perovskite structure.     

High-pressure phases of Co2GeO4, Ni2GeO4, Mn2GeO4 and MnGeO3: implications for the germanate-silicate modeling scheme and the Earth's mantle

In a diamond-anvil press coupled with YAG laser heating, the spinels of Co₂GeO₄ and Ni₂GeO₄ have been found to disproportionate into their isochemical oxide mixtures at about 250 kbar and 1400–1800°C in the same manner as their silicate analogues. At about the same P-T conditions MnGeO₃ transforms to the orthorhombic perovskite structure (space group Pbnm); the lattice parameters at room temperature and 1 bar are a₀ = 5.084 ± 0.002, b₀ = 5.214 ± 0.002, and c₀ = 7.323 ± 0.003Å with Z = 4 for the perovskite phase. The zero-pressure volume change associated with the ilmenite-perovskite phase transition in MnGeO₃ is ?6.6%. Mn₂GeO₄ disproportionates into a mixture of the perovskite phase of MnGeO₃ plus the rocksalt phase of MnO at P = 250kbar and T = 1400–1800°C. The concept of utilizing germanates as high-pressure models for silicates is valid in general. The results of this study support the previous conclusion that the lower mantle comprises predominantly the orthorhombic perovskite phase of ferromagnesian silicate.     

Pyroxene-garnet solid-solution equilibria in the systems Mg4Si4O12Mg3Al2Si3O12 and Fe4Si4O12Fe3Al2Si3O12 at high pressures and temperatures

Pyroxene-garnet solid-solution equilibria have been studied in the pressure range 41–200 kbar and over the temperature range 850–1,450°C for the system Mg₄Si₄O₁₂Mg₃Al₂Si₃O₁₂, and in the pressure range 30–105 kbar and over the temperature range 1,000–1,300°C for the system Fe₄Si₄O₁₂Fe₃Al₂Si₃O₁₂. At 1,000°C, the solid solubility of enstatite (MgSiO₃) in pyrope (Mg₃Al₂Si₃O₁₂) increases gradually to 140 kbar and then increases suddenly in the pressure range 140–175 kbar, resulting in the formation of a homogeneous garnet with composition Mg₃(Al_0.8Mg_0.6Si_0.6)Si₃O₁₂. In the MgSiO₃-rich field, the three-phase assemblage of β- or γ-Mg₂SiO₄, stishovite and a garnet solid solution is stable at pressures above 175 kbar at 1,000°C. The system Fe₄Si₄O₁₂Fe₃Al₂Si₃O₁₂ shows a similar trend of high-pressure transformations: the maximum solubility of ferrosilite (FeSiO₃) in almandine (Fe₃Al₂Si₃O₁₂) forming a homogeneous garnet solid solution is 40 mol% at 93 kbar and 1,000°C.If a pyrolite mantle is assumed, from the present results, the following transformation scheme is suggested for the pyroxene-garnet assemblage in the mantle. Pyroxenes begin to react with the already present pyrope-rich garnet at depths around 150 km. Although the pyroxene-garnet transformation is spread over more than 400 km in depth, the most effective transition to a complex garnet solid solution takes place at depths between 450 and 540 km. The complex garnet solid solution is expected to be stable at depths between 540 and 590 km. At greater depths, it will decompose to a mixture of modified spinel or spinel, stishovite and garnet solid solutions with smaller amounts of a pyroxene component in solution.     

A new high-pressure phase of Ca2Al2SiO7 and implications for the earth''s interior

Gehlenite (Ca₂Al₂SiO₇) has been found to transform to a new phase at pressures greater than 100 kbar and at about 1000°C, using a diamond-anvil pressure cell coupled with laser heating. The atoms of the new phase appear to be arranged in a perovskite-related structure similar to that described for Na₂Ti₃O₇. The structure probably consists of layers of (Al₂SiO₇)⁴⁻, which are built up from blocks of edge-sharing (Al, Si)O₆ octahedra and these blocks are joined by common octahedra corners. A small cubic unit cell with a = 3.719 ± 0.004 Å indexes completely the strong lines of the powder diffraction pattern, and a superlattice with a = 14.88 ± 0.02 Å satisfies all the observed weak lines in addition to the strong ones. However, the cell may be pseudocubic. The small cell contains a half of the gehlenite formula while the large cell contains 32 gehlenite formulae. Hence the molar volume for the new phase of Ca₂Al₂SiO₇ is calculated to be 61.96 ± 0.20 cm³ at atmospheric pressure and room temperature. The new sodium titanate-type structure is probably more closely packed than an ordinary perovskite-type structure in which all octahedral corners are shared. This view is strongly supported by the very great density of this new phase, which is about 8% denser than the equivalent mixture of CaAl₂O₄ (calcium ferrite type) plus CaSiO₃ (cubic perovskite type). The new phase is probably the most closely packed silicate known. Mg₂SiO₄ (spinel) was found to transform to an assemblage containing MgSiO₃ (perovskite) plus MgO (periclase) at P-T conditions equivalent to the upper part of the lower mantle. By reacting with MgO, the perovskite modification of both MgSiO₃ and MgSiO₃ · xAl₂O₃ may adopt the sodium titanate structure at the still greater depths of the lower mantle. If the sodium titanate structures of Mg₂(Al₂Si)O₇ and Mg₂(MgSi₂)O₇ are present in the deep part of the lower mantle, MgO does not exist as a separate phase at the mantle-core boundary. This might be an obstacle to the possibility of dissolving these oxides (specifically the FeO component) in the molten Fe in the outer core as suggested by geophysical and geochemical studies of the earth's interior. The mechanism for developing the chemical plumes in the deep mantle proposed by Anderson does not appear to be consistent with studies of phase transformations in Ca-Al-rich compounds as outlined in this paper.     

High-pressure phase transformations in the system CaSiO3-Al2O3

Phase assemblages for five selected compositions in the system CaSiO₃-Al₂O₃ have been investigated in the pressure range 100–300 kbar and at about 1000°C in a diamond-anvil press coupled with laser heating. At pressures below about 250 kbar, the assemblage of grossularite plus corundum is stable for compositions containing more than 25 mole% Al₂O₃. Above about 250 kbar, phase assemblages for the latter compositions are truncated by those in the join CaAl₂O₄-SiO₂. Garnet solid solutions are stable between about 10 and 25 mole% Al₂O₃. Grossularite transforms to a new tetragonal form at pressures greater than about 250 kbar, but the stability field for the garnet solid solutions extends to pressures up to about 300 kbar. The perovskite modification appears to be stable at pressures above about 150 kbar, but is probably limited to nearly pure CaSiO₃ composition. Phase behaviour for calcium-bearing silicates or aluminosilicates in the lower mantle are apparently more complicated than was suggested earlier.     

High-pressure reconnaissance investigation in the system Mg3Al2Si3O12-Fe3Al2Si3O12

Two synthetic end-members and two natural solid solutions of the system Mg₃Al₂Si₃O₁₂-Fe₃Al₂Si₃O₁₂ have been found to display successive phase transformations at increasingly high pressures when they were compressed in a diamond-anvil cell and heated with a YAG laser to temperatures of approximately 1400–1800°C. X-ray diffraction studies of the quenched samples show that the iron-rich garnets apparently first transform to a garnet-related high-pressure phase, then disproportionate into a mixture of magnesiowüstite plus an unknown phase(s). The latter phase(s) may further transform to a still denser unknown phase(s). The ultimate high-pressure phase may be a perovskite-like structure as was previously found for the magnesium-rich garnets. One of the unknown phases may be the high-pressure phase of Al₂O₃ · nSiO₂ compounds. Magnesium-rich garnets display similar phase transformations as do the iron-rich garnets with the exception of the garnet-related high-pressure phase. These results disagree with a previous interpretation for the high-pressure phase of iron-silicate garnets recovered in shock-wave experiments reported by Ahrens and Graham (1972).     

High-pressure synthesis of ZnSiO3 ilmenite

High-pressure phase relations in ZnSiO₃ and Zn₂SiO₄ were investigated at about 1000°C and in the pressure range of 100–500 kbar, using a double-staged split-sphere type of high-pressure apparatus.Clinopyroxene-type ZnSiO₃ transforms directly into a polymorph with the ilmenite structure at 120 kbar. The hexagonal unit cell dimensions of the ZnSiO₃ ilmenite are determined to be $\text{a = 4.746 ± 0.001}\text{A}\text{?}\text{and}\text{c = 13.755 ± 0.002}\text{A}\text{?}$ under ambient conditions.The following reactions are also recognized at about 1000°C:

and:

The stabilities of silicate ilmenites, especially the absence of ilmenite of transition metal silicate composition, is discussed. It is pointed out that data on phase relations in zinc silicates may be informative for the consideration on those in magnesium silicates under very high pressures. It is suggested that the silicate ilmenite may be a major constituent in the lower mantle.     

The location of Fe3+ ions in forsterite (Mg2SiO4)

The electron spin resonance spectrum of Fe³⁺ in a single crystal of forsterite was studied. Two distinct patterns of about equal intensities were observed which are due to Fe³⁺ at two distinct positions with 4a (M1)and4c (M2or Si) symmetry of Pbnm. The assignment of the 4c pattern to Fe³⁺ ions at the Si position cannot be excluded by symmetry but it is unlikely. The Hamiltonian parameters A and E/D are consistent with the conclusion that Fe³⁺ in this crystal is disordered over two distinct octahedral positions.     

Phase transition in CaSiO3 perovskite

A phase transition in pure CaSiO₃ perovskite was investigated at 27 to 72 GPa and 300 to 819 K by in-situ X-ray diffraction experiments in an externally-heated diamond-anvil cell. The results show that CaSiO₃ perovskite takes a tetragonal form at 300 K and undergoes phase transition to a cubic structure above 490–580 K in a pressure range studied here. The transition boundary is strongly temperature-dependent with a slightly positive dT / dP slope of 1.1 (± 1.3) K/GPa. It is known that the transition temperature depends on Al₂O₃ content dissolved in CaSiO₃ perovskite [Kurashina et al., Phys. Earth Planet. Inter. 145 (2004) 67–74]. The phase transition in CaSiO₃(+ 3 wt.% Al₂O₃) perovskite therefore could occur in a cold subducted mid-oceanic ridge basalt (MORB) crust at about 1200 K in the upper- to mid-lower mantle. This phase transition is possibly ferroelastic-type and may cause large seismic anomalies in a wide depth range.     

Radiation-induced decomposition of U(VI) phases to nanocrystals of UO2

U⁶⁺-phases are common alteration products, under oxidizing conditions, of uraninite and the UO₂ in spent nuclear fuel. These U⁶⁺-phases are subjected to a radiation field caused by the α-decay of U, or in the case of spent nuclear fuel, incorporated actinides, such as ²³⁹Pu and ²³⁷Np. In order to evaluate the effects of α-decay events on the stability of the U⁶⁺-phases, we report, for the first time, the results of ion beam irradiations (1.0 MeV Kr²⁺) of U⁶⁺-phases. The heavy-particle irradiations are used to simulate the ballistic interactions of the recoil-nucleus of an α-decay event with the surrounding structure. The Kr²⁺-irradiation decomposed the U⁶⁺-phases to UO₂ nanocrystals at doses as low as 0.006 displacements per atom (dpa). U⁶⁺-phases accumulate substantial radiation doses (∼1.0 displacement per atom) within 100,000 yr if the concentration of incorporated ²³⁹Pu is as high as 1 wt.%. Similar nanocrystals of UO₂ were observed in samples from the natural fission reactors at Oklo, Gabon. Multiple cycles of radiation-induced decomposition to UO₂ followed by alteration to U⁶⁺-phases provide a mechanism for the remobilization of incorporated radionuclides.     

Phosphate in Angra dos Reis: Structure and composition of the Ca3(PO4)2 minerals

Extraterrestrial calcium phosphates (“whitlockites”) have the anhydrous β-Ca₃(PO₄)₂ structure, which is different from that of hydrous terrestrial whitlockite. This has been confirmed by X-ray refinement of the structure of a phosphate from the achondrite Angra dos Reis. In the β-Ca₃(PO₄)₂ structure, there is one crystallographic site, Ca(IIA), which is half-occupied by calcium, and which seems to have an energetically unfavorable configuration; natural phosphates with this configuration (including Angra dos Reis) have composition Ca₁₉(Mg,Fe)₂(PO₄)₁₄. Stability of the structure is probably increased by substitution of Na for Ca in Ca(IIA) giving composition Ca₁₈ (Mg,Fe)₂Na₂(PO₄)₁₄, which occurs in chondrites; by vacancy of Ca(IIA), with rare earths and yttrium substituting for calcium in other sites for charge balance, giving composition Ca₁₆(Y,RE)₂(Mg,Fe)₂(PO₄)₁₄, found in lunar rocks; or by replacing Ca with hydrogen, giving composition Ca₁₈(Mg,Fe)₂H₂(PO₄)₁₄, which is terrestrial whitlockite. Lack of the favorable substitutions of Na, (Y, RE) or H in Angra dos Reis phosphate implies that these elements were relatively scarce in its environment of formation.     

Shock-induced CO2 loss from CaCO3; implications for early planetary atmospheres

We report new results of shock recovery experiments on single crystal calcite. Recovered samples are subjected to thermogravimetric analysis. This yields the maximum amount of post-shock CO₂, the decarbonization interval, ΔT, and the energy of association (or vaporization), ΔEV, for the removal of remaining CO₂ in shock-loaded calcite. Comparison of post-shock CO₂ with that initially present determines shock-induced CO₂ loss as a function of shock pressure. Incipient to complete CO₂ loss occurs over a pressure range of 10to 70GPa. The latter pressure should be considered a lower bound. Comparable to results on hydrous minerals, ΔT and ΔEV decrease systematically with increasing shock pressure. This indicates that shock loading leads to both the removal of structural volatiles and weakening of bonds between the volatile species and remainder of the crystal lattice.Optical and scanning electron microscopy (SEM) reveal structural changes, which are related to the shock-loading. Comparable to previous findings on shocked antigorite is the occurrence of dark, diffuse areas, which can be resolved as highly vesicular areas as observed with a scanning electron microscope. These areas are interpreted as representing quenched partial melts, into which shock-released CO₂ has been injected.The experimental results are used to place bonds on models of impact production of CO₂ during accretion of the terrestrial planets.     

Synthesis of a perovskite-type polymorph of CaSiO3

Synthetic crystalline (wollastonite) and glass forms of CaSiO₃ have been compressed to loading pressures above 160 kbar and heated to about 1500° C by a laser in a diamond-anvil cell. After cooling, an X-ray diffraction study carried out whilst the sample was maintained at high pressure revealed that it had transformed to a cubic perovskite-type 3olymorph with a = 3.485 ± 0.008A?. After release of pressure, however, the sample showed a mixture of glass plus a few weak lines corresponding to ε-CaSiO₃ which is thus interpreted as a retrogressive transition product. The density of the perovskite polymorph of CaSiO₃ is about 9.2% greater than that of an isochemical mixture of CaO + SiO₂ (stishovite) at about 160 kbar.     

The join grossularite-calcite through the system CaO—Al2O3—SiO2—CO2 at 30 kilobars: Crystallization range of silicates and carbonates on the liquidus

At 30 kbar, calcite melts congruently at 1615°C, and grossularite melts incongruently to liquid + gehlenite (tentative identification) at 1535°C. The assemblage calcite + grossularite melts at 1450°C to produce liquid + vapor, with piercing point at about 49 wt.% CaCO₃. Vapor phase is present in all hypersolidus phase fields except for those with less than about 7% CaCO₃ or 8% Ca₃Al₂Si₃O₁₂. These results, together with known liquidus data for CaO—SiO₂—CO₂ and inferred results for CaO—Al₂O₃—CO₂ and Al₂O₃—SiO₂—CO₂, permit construction of the position of the CO₂- saturated liquidus surface in the quaternary system, and estimation of the positions of liquidus field boundaries separating some of the primary crystallization fields on this surface. The field of calcite is separated from those for grossularite and quartz by a field boundary with about 50% dissolved CaCO₃. Crystallization paths of silicate liquids in the range Ca₂SiO₄—Ca₃Al₂Si₃O₁₂—SiO₂, with some dissolved CO₂, will terminate at a quaternary eutectic on this field boundary, with the precipitation of calcite together with grossularite and quartz, at a temperature below 1450°C. Addition of Al₂O₃ to CaO—SiO₂—CO₂ in amounts sufficient to stabilize garnet thus causes little change in the general liquidus pattern as far as carbonates and silicates are concerned. With addition of MgO, we anticipate that silicate liquids with dissolved CO₂ will also follow liquidus paths to fields for the precipitation of carbonates; we conclude that similar paths link kimberlite and some carnbonatite magmas.