Experimental study of the bcc-fcc phase transformations in the Fe-rich system Fe-Si at high pressures

In situ X-ray diffraction experiments at high pressure were carried out up to 8.9 GPa and 1100 °C to study phase transformations of iron and two iron-silicon alloys Fe_0.91Si_0.09 and Fe_0.83Si_0.17. For iron, the transformation from the bcc phase to the fcc phase was observed at pressures 3.8–8.2 GPa and temperatures that are consistent with previous in situ X-ray diffraction studies. Reversal of the transformation of iron was found to be sensitive to temperature; hysteresis of the transformation increased from 25 °C at 3.8 GPa to 100 °C at 7.0 GPa, primarily because the bcc-fcc phase boundary has a negative Clayperon slope. In the binary system Fe-Si, the observations of the present study indicate that the ferrite (bcc phase)-stabilizing behavior of silicon persists at high pressures and that the maximum solubility of silicon in the fcc phase increases with increasing pressure: (1) the transformation from the bcc phase to the fcc phase was observed in Fe_0.91Si_0.09 at 6.0, 7.4 and 8.9 GPa and the temperatures measured at the onset of the transformations were 300 °C higher than those in iron at similar pressures, (2) the transformation rate in Fe_0.91Si_0.09 was extremely sluggish compared to that of iron, and (3) the bcc-fcc phase transformation was not observed in Fe_0.91Si_0.09 at 4.7 GPa up to 1000 °C and in Fe_0.83Si_0.17 at 8.2 GPa and 1100 °C. Received: 1 June 1998 / Revised, accepted: 9 October 1998     

X-ray diffraction study of magnesite at high pressure and high temperature

 P–V–T measurements on magnesite MgCO₃ have been carried out at high pressure and high temperature up to 8.6 GPa and 1285 K, using a DIA-type, cubic-anvil apparatus (SAM-85) in conjunction with in situ synchrotron X-ray powder diffraction. Precise volumes are obtained by the use of data collected above 873 K on heating and in the entire cooling cycle to minimize non-hydrostatic stress. From these data, the equation-of-state parameters are derived from various approaches based on the Birch-Murnaghan equation of state and on the relevant thermodynamic relations. With K′₀ fixed at 4, we obtain K₀=103(1) GPa, α(K⁻¹)=3.15(17)×10⁻⁵ +2.32(28)×10⁻⁸ T, (∂K_T/∂T)_P=−0.021(2) GPaK⁻¹, (dα/∂P)_T=−1.81×10⁻⁶ GPa⁻¹K⁻¹ and (∂K_T/∂T)_V= −0.007(1) GPaK⁻¹; whereas the third-order Birch-Murnaghan equation of state with K′₀ as an adjustable parameter yields the following values: K₀=108(3) GPa, K′₀=2.33(94), α(K⁻¹)=3.08(16)×10⁻⁵+2.05(27) ×10⁻⁸ T, (∂K_T/∂T)_P=−0.017(1) GPaK⁻¹, (dα/∂P)_T= −1.41×10⁻⁶ GPa⁻¹K⁻¹ and (∂K_T/∂T)_V=−0.008(1) GPaK⁻¹. Within the investigated P–T range, thermal pressure for magnesite increases linearly with temperature and is pressure (or volume) dependent. The present measurements of room-temperature bulk modulus, of its pressure derivative, and of the extrapolated zero-pressure volumes at high temperatures, are in agreement with previous single-crystal study and ultrasonic measurements, whereas (∂K_T/∂T)_P, (∂α/∂P)_T and (∂K_T/∂T)_V are determined for the first time in this compound. Using this new equation of state, thermodynamic calculations for the reactions (1) magnesite=periclase+CO₂ and (2) magnesite+enstatite=forsterite+CO₂ are consistent with existing experimental phase equilibrium data. Received September 28, 1995/Revised, accepted May 22, 1996     

Melting of iron–silicon alloy up to the core–mantle boundary pressure: implications to the thermal structure of the Earth’s core

The melting temperature of Fe–18 wt% Si alloy was determined up to 119 GPa based on a change of laser heating efficiency and the texture of the recovered samples in the laser-heated diamond anvil cell experiments. We have also investigated the subsolidus phase relations of Fe–18 wt% Si alloy by the in-situ X-ray diffraction method and confirmed that the bcc phase is stable at least up to 57 GPa and high temperature. The melting curve of the alloy was fitted by the Simon’s equation, P(GPa)/a = (T _m(K)/T ₀) ^c , with parameters, T ₀ = 1,473 K, a = 3.5 ± 1.1 GPa, and c = 4.5 ± 0.4. The melting temperature of bcc Fe–18 wt% Si alloy is comparable with that of pure iron in the pressure range of this work. The melting temperature of Fe–18 wt% Si alloy is estimated to be 3,300–3,500 K at 135 GPa, and 4,000–4,200 K at around 330 GPa, which may provide the lower bound of the temperatures at the core–mantle boundary and the inner core–outer core boundary if the light element in the core is silicon.     

Phase equilibria for the Fe2SiO4–Fe3O4 system under high pressure

High pressure phase relation of the system Fe₂SiO₄–Fe₃O₄ was investigated by synthesis experiments using multi-anvil high pressure apparatus. A complete solid solution with spinel structure along Fe₂SiO₄–Fe₃O₄ join occurs above 9 GPa at 1200 °C. Lattice constants of the solid solution show almost linear variation with composition. A spinelloid phase is stable for intermediate compositions in the pressure range from 3 to 9 GPa. the synthesized spinelloid phase is successfully indexed assuming nickel aluminosilicate V type structure. Received: October 16, 1995 / Revised, accepted: March 19, 1997     

Solubility of water in the α, β and γ phases of (Mg,Fe)2SiO4

 The solubility of hydroxyl in the α, β and γ phases of (Mg,Fe)₂SiO₄ was investigated by hydrothermally annealing single crystals of San Carlos olivine. Experiments were performed at a temperature of 1000° or 1100 °C under a confining pressure of 2.5 to 19.5 GPa in a multianvil apparatus with the oxygen fugacity buffered by the Ni:NiO solid-state reaction. Hydroxyl solubilities were determined from infrared spectra obtained of polished thin sections in crack-free regions ≤100 μm in diameter. In the α-stability field, hydroxyl solubility increases systematically with increasing confining pressure, reaching a value of ∼20,000 H/10⁶Si (1200 wt ppm H₂O) at the α-β phase boundary near 13 GPa and 1100 °C. In the β field, the hydroxyl content is ∼400,000 H/10⁶Si (24,000 wt ppm H₂O) at 14–15 GPa and 1100 °C. In the γ field, the solubility is ∼450,000 H/10⁶Si (27,000 wt ppm H₂O) at 19.5 GPa and 1100 °C. The observed dependence of hydroxyl solubility with increasing confining pressure in the α phase reflects an increase in water fugacity with increasing pressure moderated by a molar volume term associated with the incorporation of hydroxyl ions into the olivine structure. Combined with published results on the dependence of hydroxyl solubility on water fugacity, the present results for the α phase can be summarized by the relation C _OH = A(T)f nH₂Oexp(−PΔV/RT), where A(T) = 1.1 H/10⁶Si/MPa at 1100 °C, n = 1, and ΔV = 10.6×10^–6 m³/mol. These data demonstrate that the entire present-day water content of the upper mantle could be incorporated in the mineral olivine alone; therefore, a free hydrous fluid phase cannot be stable in those regions of the upper mantle with a normal concentration of hydrogen. Free hydrous fluids are restricted to special tectonic environments, such as the mantle wedge above a subduction zone. Received: 10 February 1995 / Accepted: 23 October 1995     

Melting of Fe–Ni–Si and Fe–Ni–S alloys at megabar pressures: implications for the core–mantle boundary temperature

High pressure melting behavior of three Fe-alloys containing 5 wt% Ni and (1) 10 wt% Si, (2) 15 wt% Si or (3) 12 wt% S was investigated up to megabar pressures by in situ X-ray diffraction and laser-heated diamond anvil cell techniques. We observe a decrease in melting temperature with increasing Si content over the entire investigated pressure range. This trend is used to discuss the melting curve of pure Fe. Moreover, our measurements of eutectic melting in the Fe–Fe₃S system show a change in slope around 50 GPa concomitant with the fcc–hcp phase transition in pure solid iron. Extrapolations of our melting curve up to the core–mantle boundary pressure yield values of 3,600–3,750 K for the freezing temperature of plausible outer core compositions.     

Electric conductivity of Fe2SiO4–Fe3O4 spinel solid solutions

 The spinel solid solution was found to exist in the whole range between Fe₃O₄ and γ-Fe₂SiO₄ at over 10 GPa. The resistivity of Fe_{3− x} Si _x O₄ (0.0<x<0.288) was measured in the temperature range of 80∼300 K by the AC impedance method. Electron hopping between Fe³⁺ and Fe²⁺ in the octahedral site of iron-rich phases gives a large electric conductivity at room temperature. The activation energy of the electron hopping becomes larger with increasing γ-Fe₂SiO₄ component. A nonlinear change in electric conductivity is not simply caused by the statistical probability of Fe³⁺–Fe²⁺ electron hopping with increasing the total Si content. This is probably because a large number of Si⁴⁺ ions occupies the octahedral site and the adjacent Fe²⁺ keeping the local electric neutrality around Si⁴⁺ makes a cluster, which generates a local deformation by Si substitution. The temperature dependence of the conductivity of solid solutions indicates the Verwey transition temperature, which decreases from 124(±2) K at x=0 (Fe₃O₄) to 102(±5) K at x=0.288, and the electric conductivity gap at the transition temperature decreases with Si⁴⁺ substitution. Received: 15 March 2000 / Accepted: 4 September 2000     

Eutectic temperatures and melting relations in the Fe–O–S system at high pressures and temperatures

We have investigated melting relations in the Fe–O–S ternary system in the pressure range of 15–27 GPa and 1873 K. Subsolidus phase relations are Fe, Fe₃S₂, and FeO up to 17 GPa and Fe, Fe₃S, and FeO above this pressure. The eutectic temperature slightly decreases from ambient pressure to 17 GPa, whereas increases above this pressure. The eutectic temperature in this study is 100 K lower than that in the Fe–S binary system. The oxygen content in the Fe–O–S eutectic liquid drops when the coexisting solid phases changes from FeS to Fe₃S₂. The cotectic lines in the ternary phase diagram lie close to the Fe–FeS binary axis. The isothermal sections indicate that oxygen solubility in the Fe–O–S liquid increases with increasing temperature, and with increasing sulfur content. The solubility of sulfur in the solid Fe has a maximum value at the eutectic temperature, and decreases with increasing temperature. Our results could have important implications for formation and composition of the Martian core.     

Fe–Mg partitioning between post-perovskite and ferropericlase in the lowermost mantle

Fe–Mg partitioning between post-perovskite and ferropericlase has been studied using a laser-heated diamond anvil cell at pressures up to 154 GPa and 2,010 K which corresponds to the conditions in the lowermost mantle. The composition of the phases in the recovered samples was determined using analytical transmission electron microscopy. Our results reveal that the Fe–Mg partition coefficient between post-perovskite and ferropericlase (K _D^PPv/Fp) increases with decreasing bulk iron content. The compositional dependence of K _D^PPv/Fp on the bulk iron content explains the inconsistency in previous studies, and the effect of the bulk iron content is the most dominant factor compared to other factors, such as temperature and aluminum content. Iron prefers ferropericlase compared to post-perovskite over a wide compositional range, whereas the iron content of post-perovskite (X _Fe^PPv, the mole fraction) does not exceed a value of 0.10. The iron-rich ferropericlase phase may have significant influence on the physical properties, such as the seismic velocity and electrical conductivity at the core–mantle boundary region.     

Density measurements of liquid Fe–Si alloys at high pressure using the sink–float method

The compositional dependence on the density of liquid Fe alloys under high pressure is important for estimating the amount of light elements in the Earth’s outer core. Here, we report on the density of liquid Fe–Si at 4 GPa and 1,923 K measured using the sink–float method and our investigation on the effect of the Si content on the density of the liquid. Our experiments show that the density of liquid Fe–Si decreases from 7.43 to 2.71 g/cm³ non-linearly with increasing Si content (0–100 at%). The molar volume of liquid Fe–Si calculated from the measured density gradually decreases in the compositional range 0–50 at% Si, and increases in the range 50–100 at% Si. It should be noted that the estimated molar volume of the alloys shows a negative volume of mixing between Fe and Si. This behaviour is similar to Fe–S liquid (Nishida et al. in Phys Chem Miner 35:417–423, 2008). However, the excess molar volume of mixing for the liquid Fe–Si is smaller than that of liquid Fe–S. The light element contents in the outer core estimated previously may be an underestimation if we take into account the possible negative value of the excess mixing volume of iron–light element alloys in the outer core.     

Elasticity of TiO2 rutile to 1800 K

The ambient pressure elastic properties of single-crystal TiO₂ rutile are reported from room temperature (RT) to 1800 K, extending by more than 1200 oK the maximum temperature for which rutile elasticity data are available. The magnitudes of the temperature derivatives decrease with increasing temperature for five of the six adiabatic elastic moduli (C _ij ). At RT, we find (units, GPa): C ₁₁=268(1); C ₃₃=484(2); C ₄₄=123.8(2); C ₆₆=190.2(5); C ₂₃=147(1); and C ₁₂=175(1). The temperature derivatives (units, GPa K⁻¹) at RT are: (∂C ₁₁/∂T) _P =−0.042(5); (∂C ₃₃/∂T) _P =−0.087(6); (∂C ₄₄/∂T) _P =−0.0187(2); (∂C ₆₆/∂T) _P =−0.067(2); (∂C ₂₃/∂T) _P =−0.025; and (∂C ₁₂/∂T) _P −0.048(5). The values for K _S (adiabatic bulk modulus) and μ (isotropic shear modulus) and their temperature derivatives are K _S =212(1) GPa; μ=113(1) GPa; (∂K _S /∂T) _P =−0.040(4) GPa K⁻¹; and (∂μ/∂T) _P =−0.018(1) GPa K⁻¹. We calculate several dimensionless parameters over a large temperature range using our new data. The unusually high values for the Anderson-Gròneisen parameters at room temperature decrease with increasing temperature. At high T, however, these parameters are still well above those for most other oxides. We also find that for TiO₂, anharmonicity, as evidenced by a non-zero value of [∂ln (K _T )/∂lnV] _T , is insignificant at high T, implying that for the TiO₂ analogue of stishovite, thermal pressure is independent of volume (or pressure). Systematic relations indicate that ∂² K _S /∂T∂P is as high as 7×10⁻⁴ K⁻¹ for rutile, whereas ∂²μ/∂T∂P is an order of magnitude less. Received: 19 September 1997 / Revised, accepted: 27 February 1998     

Fulgurite morphology: a classification scheme and clues to formation

Fulgurites are natural glasses formed by cloud-to-ground lightning. Several different morphologies of fulgurites have been reported in previous studies, including sand fulgurites, rock fulgurites, and clay fulgurites. Herein, we examine sand, clay, and caliche fulgurites and demonstrate that these differ systematically in their morphology. We further use morphological features to constrain properties of fulgurite-forming lightning strikes. We classify fulgurites into four types of morphologies with an additional minor type. Type I fulgurites are sand fulgurites consisting of thin, glass walls; type II fulgurites are clay fulgurites, consisting of thick, melt rich walls; type III fulgurites are caliche fulgurites, consisting of thick, glass poor walls; and type IV fulgurites are rock fulgurites, consisting of glasses with walls consisting of surrounding, unmelted rock. Fulgurite morphology shows that the energy of fulgurite-forming strikes is between 1 and 30 MJ/m of fulgurite formed, suggests heating rates in the order of 1,000 K/s, and lightning channel thicknesses of about 1 mm diameter. Lightning generates mixtures of at least two components in most fulgurites: an SiO₂ glass identified as lechatelierite and a groundmass of more varied composition. In addition to these four primary types, a fifth type—droplet fulgurites—is morphologically dissimilar from the other types, but is compositionally related to the type II or IV fulgurites. Additionally, two fulgurites, both from York County, Pennsylvania, USA, showed the reduction of iron to iron metal with an assortment of Fe–Ti and Si–P compounds with stoichiometry that ranges from nearly pure Fe metal to FeSi. These metal silicides include stoichiometric Fe₃Si, Fe₂Si, and Fe₅Si₃, and possibly Fe₈Si₃ and Fe₇Si₃, and provide a terrestrial source for these phases, which are typically associated with extraterrestrial material.     

High-pressure and high-temperature in situ X-ray diffraction study of iron and corundum to 68 GPa using an internally heated diamond anvil cell

Using powder X-ray diffraction of heated solids to pressures reaching 68 GPa, the pressure-volume-temperature (PVT) data on corundum Al₂O₃ and ɛ-Fe were determined with the following results: *Corundum,*Iron, *Al₂O₃*ɛ-Fe Isothermal bulk*258 (2)*164 (3)  modulus K'_{300, 1} (GPa) Pressure derivative K_{300, 1}*4.88 (4)*5.36 (16) Temperature derivative*–0.020 (2)*–0.043 (3)  (∂K _T,1 /∂T) _P (GPa/K) Molar volume V_300,1*25.59 (2)*6.76 (2)  (cm³/mol) Isobaric thermal expansion at 1 atm (0.101 MPa) is given by (K^–1): α _T =2.6 (2) 10^–5+1.81 (9) 10^–9 T–0.67 (6)/T ² for corundum, and α _T =5.7 (4) 10^–5+4.2 (4) 10^–9 T–0.17 (7)/T ² for iron ɛ-Fe. Received: 1 March 1997 / Revised, accepted: 21 August 1997     

Hydrostatic compression and crystal structure of pyrope to 33 GPa

The equation of state and crystal structure of pyrope were determined by single crystal X-ray diffraction under hydrostatic conditions to 33 GPa, a pressure that corresponds to a depth of about 900 km in the lower mantle. The bulk modulus K _T0 and its pressure derivative K ^' _T0 were determined simultaneously from an unweighted fit of the volume data at different pressures to a third order Birch-Murnaghan equation of state. They are 171(2) GPa and 4.4(2), respectively. Over the whole pressure range, MgO₈ polyhedra showed the largest compression of 18.10(8)%, followed by AlO₆ and SiO₄ polyhedra, with compression of 11.7(1)% and 4.6(1)%, respectively. The polyhedral bulk moduli for MgO₈, AlO₆ and SiO₄ are 107(1), 211(11) and 580(24) GPa, respectively, with K ^' _T0 fixed to 4. Significant compression of up to 1.8(1)% in the very rigid Si−O bonding in pyrope could be detected to 33 GPa. Changes in the degree of polyhedral distortion for all three types of polyhedra could also be observed. These changes could be found for the first time for AlO₆ and SiO₄ in pyrope. It seems that the compression of pyrope crystal structure is governed by the kinking of the Al−O−Si angle between the octahedra and tetrahedra. No phase transition could be detected to 33 GPa. Received: 24 March 1997 / Revised, accepted: 29 July 1997     

Phase relations in Fe–Ni–C system at high pressures and temperatures

We performed comparative study of phase relations in Fe_1−x Ni _x (0.10 ≤ x ≤ 0.22 atomic fraction) and Fe_0.90Ni_0.10−x C _x (0.1 ≤ x ≤ 0.5 atomic fraction) systems at pressures to 45 GPa and temperatures to 2,600 K using laser-heated diamond anvil cell and large-volume press (LVP) techniques. We show that laser heating of Fe,Ni alloys in DAC even to relatively low temperatures can lead to the contamination of the sample with the carbon coming from diamond anvils, which results in the decomposition of the alloy into iron- and nickel-rich phases. Based on the results of LVP experiments with Fe–Ni–C system (at pressures up to 20 GPa and temperatures to 2,300 K) we demonstrate decrease of carbon solubility in Fe,Ni alloy with pressure.     

Si in the core? New high-pressure and high-temperature experimental data

High-pressure high-temperature experiments have been carried out up to 25 GPa and 2200°C in a multianvil press on assemblages made of silicates and iron-silicon alloys. At 20 GPa, silicon is extracted from the metal phase, forming stishovite reaction rims around metal grains. The silicon content in metal has been measured by analytical electron microscopy and electron microprobe. In contrast with earlier experiments, the present data were obtained by using silicon-rich metal alloys as starting materials instead of studying incorporation of silicon in initially silicon-free metal. As in most of previous studies carried out below 25 GPa, the silicon content in liquid metal increases with increasing pressure and with decreasing oxygen fugacity. The oxygen fugacity in most experiments was calculated by using two independent buffers: iron/?stite (IW) and SiO₂/Si, allowing to link consistently the Fe contents in silicates, the Si contents in metal and the temperatures of the experiments. At oxygen fugacities 4 log units below IW, silicates are in equilibrium with Si-rich metallic alloys (up to 17 wt% of Si in metal at 20 GPa and 2200°C). Extrapolation to 2 log units below IW leads to less than 0.1 wt% Si in the metal phase. Presence of several wt% of silicon in the Earth’s core thus requires highly reduced initial materials that, if equilibrated at conditions relevant to small planets, should already contain significant amount of silicon dissolved in metal.     

Precambrian banded iron formations in the North China Craton: Silicon and oxygen isotopes and genetic implications

Banded iron formations (BIFs) are Precambrian chemical marine sedimentary formations that record major transitions of chemical composition, and oxidation–reduction state of oceans at the time of their deposition. In this paper, we report silicon and oxygen isotope compositions of a variety of BIFs from the North China Craton (NCC) in order to deduce the mechanism of their formation. Quartz in the various types of BIFs from the NCC are generally depleted in ³⁰Si, where δ³⁰Si_NBS-28 values range from − 2.0‰ to − 0.3‰ (average, − 0.8‰), similar to δ³⁰Si_NBS-28 values measured from modern submarine black chimneys and sinters. The δ¹⁸O_V-SMOW values of quartz in the BIFs are relatively high (8.1‰–21.5‰; average, 13.1‰), similar to those of siliceous rock formed by hydrothermal activities. The δ³⁰Si_NBS-28 values of quartz in magnetite bands are commonly lower than those of quartz in adjacent siliceous bands within the same sample, whereas δ¹⁸O_V-SMOW values are higher in the magnetite bands. A negative correlation is observed between δ³⁰Si_NBS-28 and δ¹⁸O_V-SMOW values of quartz from siliceous and magnetite bands in BIF from Fuping, Hebei Province. The isotopic compositions of silicon and oxygen of quartz in BIFs provide insights for the formation mechanisms of silicon–iron cyclothems in BIFs. After the silicon- and iron-rich hydrothermal solution was injected onto the seabed, the abrupt temperature drop caused oversaturation of silicic acid, resulting in rapid precipitation of SiO₂ and deposition of siliceous layers. Ferric hydroxide was precipitated later than SiO₂ because of low free-oxygen concentration in the ocean bottom. Progressive mixing of hydrothermal solution with seawater caused a continuous drop in temperature and an increase in Eh values, resulting in gradual oxidation of hydrothermal Fe^2 + and deposition of iron-rich layers. In summary, each silicon–iron cyclothem marks a large-scale submarine hydrothermal exhalation. The periodic nature of these exhalations resulted in the formation of regular silicon–iron cyclothems. The widespread distribution of BIFs indicates that volcanism and submarine hydrothermal exhalation were extensive; the low δ³⁰Si_NBS-28 and high δ¹⁸O _V-SMOW values of the BIFs indicate that the temperature of seawater was relatively high at the time of BIF formation, and that concentrations of Fe^2 + and H₄SiO₄ in seawater were saturated.     

The low-pressure stability of phase A,Mg7Si2O8(OH)6

 Phase A, Mg₇Si₂O₈(OH)₆, is a dense hydrous magnesium silicate whose importance as a host of H₂O in the Earth’s mantle is a subject of debate. We have investigated the low-pressure stability of phase A in experiments on the reaction phase A=brucite+forsterite. Experiments were conducted in piston-cylinder and multi-anvil apparatus, using mixtures of synthetic phase A, brucite and forsterite. The reaction was bracketed between 2.60 and 2.75 GPa at 500° C, between 3.25 and 3.48 GPa at 600° C and between 3.75 and 3.95 GPa at 650^° C. These pressures are much lower than observed in the synthesis experiments of Yamamoto and Akimoto (1977). At 750^° C the stability field of brucite + chondrodite was entered. The enthalpy of formation and entropy of phase A at 1 bar (10⁵ Pa), 298 K, were derived from the experimental brackets on the reaction phase A=brucite+forsterite using a modified version of the thermodynamic dataset THERMOCALC of Holland and Powell (1990), which includes a new equation of state of H₂O derived from the molecular dynamics simulations of Brodholt and Wood (1993). The data for phase A are: ΔH ^o _f =−7126±8 kJ mol^-1, S ^o=351 J K^-1 mol^-1. Incorporating these data into THERMOCALC allows the positions of other reactions involving phase A to be calculated, for example the reaction phase A + enstatite=forsterite+vapour, which limits the stability of phase A in equilibrium with enstatite. The calculated position of this reaction (753° C at 7 GPa to 937° C at 10 GPa) is in excellent agreement with the experimental brackets of Luth (1995) between 7 and 10 GPa, supporting the choice of equation of state of H₂O used in THERMOCALC. Comparison of our results with calculated P-T paths of subducting slabs (Peacock et al. 1994) suggests that, in the system MgO–SiO₂–H₂O, phase A could crystallise in compositions with Mg/Si>2 at pressures as low as 3 GPa. In less Mg rich compositions phase A could crystallise at pressures above approximately 6 GPa. Received: 3 July 1995/Accepted: 14 December 1995     

P-V equation of state of portlandite, Ca(OH)2, from powder neutron diffraction data

 A high pressure neutron powder diffraction study of portlandite [Ca(OH)₂] has been performed at ISIS facility (U.K.); nine spectra have been collected increasing the pressure by steps, up to 10.9 GPa, by means of a Paris-Edinburgh cell installed on the POLARIS diffractometer. The tensorial formalism of the lagrangian finite strain theory and the Birch-Murnaghan equation of state have been used to determine, independently, two values of the bulk modulus of portlandite, obtaining K ₀=38.3(±1.1) GPa [linear incompressibilities: K _0a=188.4(±9.9), K _0c=64.5(±2.5) GPa] and K ₀=34.2(±1.4) GPa, respectively. The present results comply with values from previous measurements by X-ray diffraction [K ₀=37.8(±1.8) GPa] and Brillouin spectroscopy [K ₀=31.7(±2.5) GPa]. Reasonably, Ca(OH)₂ has revealed to be bulkly softer than Mg(OH)₂ [K ₀=41(±2), K _0a=313, K _0c=57 GPa]. The Ca(OH)₂ linear incompressibility values reflect the nature of forces acting to stabilize the (001) layer structure and, further, prove that the replacement Ca/Mg mainly affects the elastic properties in the (001) plane, rather than along the [001] direction. Data from a full refinement of the structure at room pressure are reported. Received January 12, 1996/Revised, accepted June 15, 1996     

Isothermal compression behavior of (Mg,Fe)O using neon as a pressure medium

We present isothermal volume compression behavior of two polycrystalline (Mg,Fe)O samples with FeO = 39 and 78 mol% up to ~90 GPa at 300 K using synchrotron X-ray diffraction and neon as a pressure-transmitting medium. For the iron-rich (Mg_0.22Fe_0.78)O sample, a structural transition from the B1 structure to a rhombohedral structure was observed at 41.6 GPa, with no further indication of changes in structural or compression behavior changes up to 93 GPa. In contrast, a change in the compression behavior of (Mg_0.61Fe_0.39)O was observed during compression at P ≥ 71 GPa and is indicative of a spin crossover occurring in the Fe²⁺ component of (Mg_0.61Fe_0.39)O. The low-spin state exhibited a volume collapse of ~3.5%, which is a larger value than what was observed for a similar composition in a laser-heated NaCl medium. Upon decompression, the volume of the high-spin state was recovered at approximately 65 GPa. We therefore bracket the spin crossover at 65 ≤ P (GPa) ≤ 77 at 300 K (Mg_0.61Fe_0.39)O. We observed no deviation from the B1 structure in (Mg_0.61Fe_0.39)O throughout the pressure range investigated.