Multipole-refined charge density study of diopside at ambient conditions

The experimental multipole electron density, ρ(r), of diopside was derived from high-resolution single-crystal diffraction at room temperature. Its topological analysis revealed predominantly ionic Si–O bonding, as found in electron density studies of other silicates. In particular, the non-bridging Si–O bonds are slightly less ionic in character than the bridging Si–O bonds. The Ca–O and Mg–O bonds are classified as pure closed-shell ionic interactions. An analysis of –∇²ρ(r) showed the presence of maxima around the oxygen atoms, associated to lone pairs domains that are involved in bonds with the surrounding ions. Calculation of atomic basins gave net charges of –1.56(12), 3.11(17), 1.79(13) and 1.88(18) e for O (averaged), Si, Ca and Mg atoms, respectively. O···O interactions between the O atoms at the vertices of the SiO₄ tetrahedron were also detected from the topological analysis of ρ(r), and indicate a cooperative interaction among the lone pairs of neighbouring oxygen atoms. All these results were also confirmed by periodic restricted Hartree–Fock (RHF) calculations. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

LiFeSi<Subscript>2</Subscript>O<Subscript>6</Subscript> and NaFeSi<Subscript>2</Subscript>O<Subscript>6</Subscript> at low temperatures: an infrared spectroscopic study

 Synthetic aegirine LiFeSi₂O₆ and NaFeSi₂O₆ were characterized using infrared spectroscopy in the frequency range 50–2000 cm⁻¹, and at temperatures between 20 and 300 K. For the C2/c phase of LiFeSi₂O₆, 25 of the 27 predicted infrared bands and 26 of 30 predicted Raman bands are recorded at room temperature. NaFeSi₂O₆ (with symmetry C2/c) shows 25 infrared and 26 Raman bands. On cooling, the C2/c–P2₁/c structural phase transition of LiFeSi₂O₆ is characterized by the appearance of 13 additional recorded peaks. This observation indicates the enlargement of the unit cell at the transition point. The appearance of an extra band near 688 cm⁻¹ in the monoclinic P2₁/c phase, which is due to the Si–O–Si vibration in the Si₂O₆ chains, indicates that there are two non-equivalent Si sites with different Si–O bond lengths. Most significant spectral changes appear in the far-infrared region, where Li–O and Fe–O vibrations are mainly located. Infrared bands between 300 and 330 cm⁻¹ show unusually dramatic changes at temperatures far below the transition. Compared with the infrared data of NaFeSi₂O₆ measured at low temperatures, the change in LiFeSi₂O₆ is interpreted as the consequence of mode crossing in the frequency region. A generalized Landau theory was used to analyze the order parameter of the C2/c–P2₁/c phase transition, and the results suggest that the transition is close to tricritical. Received: 21 January 2002 / Accepted: 22 July 2002     

An X-ray diffraction study of shock-wave-densified SiO<Subscript>2</Subscript> glasses

 The densification and structural changes in SiO₂ glass compressed up to 43.4 GPa by shock experiments are investigated quantitatively by the X-ray diffraction technique. Direct structural data (average Si–O and Si–Si distances and Si–O–Si angles, coordination number of the Si atom) of these shock-densified SiO₂ glasses have been obtained by analyzing the radial distribution function curves, RDF(r), calculated with X-ray diffraction data. The coordination number of all densified glasses is about 4 and shows almost no pressure variation. The SiO₂ glass has shown density increase of 11% at a shock compression of 26.3 GPa. This density evolution could not be explained by the coordination change. The reduction of the average Si–O–Si angle (144° at 0 GPa to 136° at 26.3 GPa) obtained from RDF(r) data may account for this density increase. This Si–O–Si angle change may be caused by shrinkage of the network structure and the increase of small rings of SiO₄ tetrahedra. For higher shock pressure, a decrease in the Si–O–Si angle to 140° was observed. This is consistent with the decrease in density at 32.0 and 43.2 GPa. This decrease in the Si–O–Si angle and density could be attributed to an annealing effect due to high after-shock residual temperature. This pressure dependence of average Si–O–Si angles in shock-densified SiO₂ glass agrees with the results of our previous Raman spectroscopic study. On the other hand, the pressure variation for the first sharp diffraction peak (FSDP) was analyzed to estimate the evolution of intermediate range structures. It is suggested that the mean d value (d _m ) obtained from the position of FSDP strongly depends on the shock and residual temperature, as well as shock pressure. Received: 29 June 2001 / Accepted: 14 November 2001     

Atomistic model of diopside-K-jadeite (CaMgSi<Subscript>2</Subscript>O<Subscript>6</Subscript>-KAlSi<Subscript>2</Subscript>O<Subscript>6</Subscript>) solid solution

Atomistic model was proposed to describe the thermodynamics of mixing in the diopside-K-jadeite solid solution (CaMgSi₂O6-KAlSi₂O₆). The simulations were based on minimization of the lattice energies of 800 structures within a 2 × 2 × 4 supercell of C2/c diopside with the compositions between CaMgSi₂O₆ and KAlSi₂O₆ and with variable degrees of order/disorder in the arrangement of Ca/K cations in M2 site and Mg/Al in Ml site. The energy minimization was performed with the help of a force-field model. The results of the calculations were used to define a generalized Ising model, which included 37 pair interaction parameters. Isotherms of the enthalpy of mixing within the range of 273–2023 K were calculated with a Monte Carlo algorithm, while the Gibbs free energies of mixing were obtained by thermodynamic integration of the enthalpies of mixing. The calculated T-X diagram for the system CaMgSi₂O₆-KAlSi₂O₆ at temperatures below 1000 K shows several miscibility gaps, which are separated by intervals of stability of intermediate ordered compounds. At temperatures above 1000 K a homogeneous solid solution is formed. The standard thermodynamic properties of K-jadeite (KAlSi₂O₆) evaluated from quantum mechanical calculations were used to determine location of several mineral reactions with the participation of the diopside-K-jadeite solid solution. The results of the simulations suggest that the low content of KalSi₂O₆ in natural clinopyroxenes is not related to crystal chemical factors preventing isomorphism, but is determined by relatively high standard enthalpy of this end member.     

Aluminum substitution mechanisms in perovskite-type MgSiO<Subscript>3</Subscript>: an investigation by Rietveld analysis

Al-containing MgSiO₃ perovskites of four different compositions were synthesized at 27 GPa and 1,873 K using a Kawai-type high-pressure apparatus: stoichiometric compositions of Mg_0.975Si_0.975Al_0.05O₃ and Mg_0.95Si_0.95Al_0.10O₃ considering only coupled substitution Mg²⁺ + Si⁴⁺ = 2Al³⁺, and nonstoichiometric compositions of Mg_0.99Si_0.96Al_0.05O_2.985 and Mg_0.97Si_0.93Al_0.10O_2.98 taking account of not only the coupled substitution but also oxygen vacancy substitution 2Si⁴⁺ = 2Al³⁺ + V_O¨. Using the X-ray diffraction profiles, Rietveld analyses were performed, and the results were compared between the stoichiometric and nonstoichiometric perovskites. Lattice parameter–composition relations, in space group Pbnm, were obtained as follows. The a parameters of both of the stoichiometric and nonstoichiometric perovskites are almost constant in the X _Al range of 0–0.05, where X _Al is Al number on the basis of total cation of two (X _Al = 2Al/(Mg + Si + Al)), and decrease with further increasing X _Al. The b and c parameters of the stoichiometric perovskites increase linearly with increasing Al content. The change in the b parameter of the nonstoichiometric perovskites with Al content is the same as that of the stoichiometric perovskites within the uncertainties. The c parameter of the nonstoichiometric perovskites is slightly smaller than that of the stoichiometric perovskites at X _Al of 0.10, though they are the same as each other at X _Al of 0.05. The Si(Al)–O1 distance, Si(Al)–O1–Si(Al) angle and minimum Mg(Al)–O distance of the nonstoichiometric perovskites keep almost constant up to X _Al of 0.05, and then the Si(Al)–O1 increases and both of the Si(Al)–O1–Si(Al) angle and minimum Mg(Al)–O decrease with further Al substitution. These results suggest that the oxygen vacancy substitution may be superior to the coupled substitution up to X _Al of about 0.05 and that more Al could be substituted only by the coupled substitution at 27 GPa. The Si(Al)–O1 distance and one of two independent Si(Al)–O2 distances in Si(Al)O₆ octahedra in the nonstoichiometric perovskites are always shorter than those in the stoichiometric perovskite at the same Al content. These results imply that oxygen defects may exist in the nonstoichiometric perovskites and distribute randomly.     

Thermodynamics of mixing in diopside–jadeite,CaMgSi<Subscript>2</Subscript>O<Subscript>6</Subscript>–NaAlSi<Subscript>2</Subscript>O<Subscript>6</Subscript>, solid solution from static lattice energy calculations

Static lattice energy calculations (SLEC), based on empirical interatomic potentials, have been performed for a set of 800 different structures in a 2 × 2 × 4 supercell of C2/c diopside with compositions between diopside and jadeite, and with different states of order of the exchangeable Na/Ca and Mg/Al cations. Excess static energies of these structures have been cluster expanded in a basis set of 37 pair-interaction parameters. These parameters have been used to constrain Monte Carlo simulations of temperature-dependent properties in the range of 273–2,023 K and to calculate a temperature–composition phase diagram. The simulations predict the order–disorder transition in omphacite at 1,150 ± 20°C in good agreement with the experimental data of Carpenter (Mineral Petrol 78:433–440, 1981). The stronger ordering of Mg/Al within the M1 site than of Ca/Na in the M2 site is attributed to the shorter M1–M1 nearest-neighbor distance, and, consequently, the stronger ordering force. The comparison of the simulated relationship between the order parameters corresponding to M1 and M2 sites with the X-ray refinement data on natural omphacites (Boffa Ballaran et al. in Am Mineral 83:419–433, 1998) suggests that the cation ordering becomes kinetically ineffective at about 600°C.     

Point defects and cation tracer diffusion in (Ti<Subscript><Emphasis Type="Italic">x</Emphasis></Subscript> Fe<Subscript><Emphasis Type="Italic">1−x</Emphasis></Subscript>)<Subscript>3−δ</Subscript>O<Subscript>4</Subscript>. II. Cation tracer diffusion

 Cation tracer diffusion coefficients, D_Me ^*, for Me=Fe, Mn, Co and Ti, were measured using radioactive isotopes in the spinel solid solution (Ti _x Fe _1−x )_3−δO₄ as a function of the oxygen activity. Experiments were performed at different cationic compositions (x=0, 0.1, 0.2 and 0.3) at 1100, 1200, 1300 and 1400 °C. The oxygen activity dependence of all data for D_Me ^* at constant temperature and cationic composition can be described by equations of the type D_Me ^*=D^∘ _Me[V]. C_V·a _O2 ^2/3+D_Me[I] ^∘·a _O2 ^−2/3·D_Me[V] ^∘ and D_Me[I] ^∘ are constants and C_V is a factor of the order of unity which decreases with increasing δ. All log D_Me ^* vs. loga _O2 curves obtained for different values of x and for different temperatures go through a minimum due to a change in the type of point defects dominating the cation diffusion with oxygen activity. Cation vacancies prevail for the cation diffusion at high oxygen activities while cation interstitials become dominant at low oxygen activities. At constant values of x, D_Me[V] ^∘ decreases with increasing temperature while D_Me[I] ^∘ increases.     

In situ X-ray observations of phase transitions in MgAl<Subscript>2</Subscript>O<Subscript>4</Subscript> spinel to 40 GPa using multianvil apparatus with sintered diamond anvils

 Phase transitions in MgAl₂O₄ spinel have been studied at pressures 22–38 GPa, and at temperatures up to 1600 °C, using a combination of synchrotron radiation and a multianvil apparatus with sintered diamond anvils. Spinel dissociated into a mixture of MgO plus Al₂O₃ at pressures to 25 GPa, while it transformed to the CaFe₂O₄ (calcium ferrite) structure at higher pressures via the metastably formed oxide mixture upon increasing temperature. Neither the e-phase nor the CaTi₂O₄-type MgAl₂O₄, which were reported in earlier studies using the diamond-anvil cell, were observed in the present pressure and temperature range. The zero-pressure bulk modulus of the calcium-ferrite-type MgAl₂O₄ was calculated as K=213 (3) GPa, which is significantly lower than that reported by Yutani et al. (1997), but is consistent with a more recent result by Funamori et al. (1998) and that estimated by an ab initio calculation by Catti (2001). Received: 2 April 2002 / Accepted: 29 July 2002 Acknowledgements The authors thank Y. Higo, Y. Sueda, T.␣Ueda, Y. Tanimoto, A. Fukuyama, K. Ochi, F. Kurio and T. Kawahara for help in the in situ X-ray observations at SPring-8 (No: 2000A0061-CD-np and 2000B0093-ND-np). We also thank W.␣Utsumi, J. Ando and O. Shimomura for advice and encouragement during this study, and N. Funamori and an anonymous reviwer for comments on the article. The present study is partly supported by the grant-in-aid for Scientific Research (A) of the Ministry of Education, Science, Sport and Culture of the Japanese government (no: 11694088).     

H<Subscript>2</Subscript>O storage capacity of MgSiO<Subscript>3</Subscript> clinoenstatite at 8–13 GPa, 1,100–1,400°C

We present H₂O analyses of MgSiO₃ pyroxene crystals quenched from hydrous conditions in the presence of olivine or wadsleyite at 8–13.4 GPa and 1,100–1,400°C. Raman spectroscopy shows that all pyroxenes have low clinoenstatite structure, which we infer to indicate that the crystals were high clinoenstatite (C2/c) during conditions of synthesis. H₂O analyses were performed by secondary ion mass spectrometry and confirmed by unpolarized Fourier transform infrared spectroscopy on randomly oriented crystals. Measured H₂O concentrations increase with pressure and range from 0.08 wt.% H₂O at 8 GPa and 1,300°C up to 0.67 wt.% at 13.4 GPa and 1,300°C. At fixed pressure, H₂O storage capacity diminishes with increasing temperature and the magnitude of this effect increases with pressure. This trend, which we attribute to diminishing activity of H₂O in coexisting fluids as the proportion of dissolved silicate increases, is opposite to that observed previously at low pressure. We observe clinoenstatite 1.4 GPa below the pressure stability of clinoenstatite under nominally dry conditions. This stabilization of clinoenstatite relative to orthoenstatite under hydrous conditions is likely owing to preferential substitution of H₂O into the high clinoenstatite polymorph. At 8–11 GPa and 1,200–1,400°C, observed H₂O partitioning between olivine and clinoenstatite gives values of D ^ol/CEn between 0.65 and 0.87. At 13 GPa and 1,300°C, partitioning between wadsleyite and clinoenstatite, D ^wd/CEn, gives a value of 2.8 ± 0.4.     

EPR study of coordination of Ag and Pb cations in BaB<Subscript>2</Subscript>O<Subscript>4</Subscript> crystals and barium borate glasses

It is shown the possibility to determine the coordination of paramagnetic ions in disordered solid structures, e.g., in barium borate glasses. For this purpose the electron paramagnetic resonance (EPR) method was used to study α-and β-BaB₂O₄ crystals and glasses of 45·BaO × 55·B₂O₃ and 40·BaO × 60·B₂O₃ (mol%) composition activated by Ag⁺ and Pb²⁺ ions. After the samples were exposed to X-rays at 77 K, different EPR centers were observed in them. In α-and β-BaB₂O₄ crystals and glasses the EPR centers Ag²⁺, Ag⁰, Pb⁺, Pb³⁺, and hole centers of O⁻ type were studied. The EPR parameters of these centers and their arrangement in crystal structure were determined. It is shown that Pb³⁺ ions in β-BaB₂O₄ crystals occupy Ba²⁺ position in an irregular polyhedron from the eight oxygen, whereas in α-BaB₂O₄ crystals they occupy Bа₂ position in a sixfold coordination. Pb⁺ ions in α-BaB₂O₄ crystals occupy Bа₁ position in a ninefold coordination from oxygen. In barium borate glasses, Pb³⁺ ions were studied in coordination polyhedron from six oxygen atoms and in a polyhedron from nine to ten oxygen atoms. It is assumed that the established difference in the structural position of Pb³⁺ ions in glasses is due to their previous incorporation in associative cation–anion complexes (AC) and “free” structure-forming cations (FC). Computer simulations have been performed to analyze the stability of specific associative complexes and to compare their bond lengths with experimental data.     

Potential protonation sites in the Al<Subscript>2</Subscript>SiO<Subscript>5</Subscript> polymorphs based on polarized FTIR spectroscopy and properties of the electron density distribution

Potential protonation sites for, kyanite, sillimanite, and andalusite, located in a mapping of the (3, −3) critical points displayed by their L(r) = −∇²ρ(r) distributions, are compared with polarized single-crystal FTIR spectra of kyanite and sillimanite determined earlier and with andalusite measured in this study. For andalusite, seven peaks were observed when the electric vector, E, is parallel to [100]: four intense ones at 3,440, 3,460, 3,526, and 3,597 cm⁻¹ and three weaker ones at 3,480, 3,520, and 3,653 cm⁻¹. Six peaks, three intense ones at 3,440, 3,460, and 3,526 cm⁻¹ and three weaker ones at 3,480, 3,520, and 3,653 cm⁻¹ when E parallels [010]. No peaks were observed when E is parallel to [001]. The concentration of water in andalusite varies between 110 and 168 ppm by weight % H₂O. Polarized FTIR spectra indicate that the OH vector is parallel to (001) in andalusite and sillimanite and in kyanite. Examination of the L(r) (3, −3) critical points in comparison with the polarized FTIR indicates that H prefers to bond to the oxygen atoms O1 and O2 in andalusite and O2 and O4 in sillimanite which correspond to the underbonded oxygen atoms and those with the largest L(r) maxima. In kyanite, comparison of the FTIR spectrum and the critical points indicates that H will preferentially bond to the two 4-coordinated O2 and O6 atoms.     

High-pressure Raman studies and heat capacity measurements on the MgSiO<Subscript>3</Subscript> analogue CaIr<Subscript>0.5</Subscript>Pt<Subscript>0.5</Subscript>O<Subscript>3</Subscript>

Raman spectroscopy and heat capacity measurements have been used to study the post-perovskite phase of CaIr_0.5Pt_0.5O₃, recovered from synthesis at a pressure of 15 GPa. Laser heating CaIr_0.5Pt_0.5O₃ to 1,900 K at 60 GPa produces a new perovskite phase which is not recoverable and reverts to the post-perovskite polymorph between 20 and 9 GPa on decompression. This implies that Pt-rich CaIr_1−xPt_xO₃ perovskites including the end member CaPtO₃ cannot easily be recovered to ambient pressure from high P–T synthesis. We estimate an increase in the thermodynamic Grüneisen parameter across the post-perovskite to perovskite transition of 34%, of similar magnitude to those for (Mg,Fe)SiO₃ and MgGeO₃, suggesting that CaIr_0.5Pt_0.5O₃ is a promising analogue for experimental studies of the competition in energetics between perovskite and post-perovskite phases of magnesium silicates in Earth’s lowermost mantle. Low-temperature heat capacity measurements show that CaIrO₃ has a significant Sommerfeld coefficient of 11.7 mJ/mol K² and an entropy change of only 1.1% of Rln2 at the 108 K Curie transition, consistent with the near-itinerant electron magnetism. Heat capacity results for post-perovskite CaIr_0.5Rh_0.5O₃ are also reported.     

Electron density distributions and bonded interactions for the fibrous zeolites natrolite, mesolite and scolecite and related materials

 For the fibrous zeolites natrolite, Na₂[Al₂Si₃O₁₀]·2H₂O, mesolite, Na₂Ca₂[Al₂Si₃O₁₀]₃·8H₂O, and scolecite, Ca[Al₂Si₃O₁₀]·3H₂O, with topologically identical aluminosilicate framework structures, accurate single-crystal X-ray diffraction data have been analyzed by least-squares refinements using generalized scattering factor (GSF) models. The final agreement indices were R(F ) = 0.0061, 0.0165, and 0.0073, respectively. Ensuing calculations of static deformation [Δρ(r)], and total, [ρ(r)], model electron density distributions served to study chemical bonding, in particular by topological electron density analyses yielding bond critical point (bcp) properties and in situ cation electronegativities. The results for 32 SiO, 24 AlO, 14 CaO, and 12 NaO unique bonds are compiled and analyzed in terms of both mean values and correlations between bond lengths, bonded oxygen radii, bcp densities, curvatures at the bcps, and electronegativities. Comparison with recent literature data obtained from both experimental electron density studies on minerals and model calculations for geometry-optimized molecules shows that the majority of the present findings conforms well with chemical expectation and with the trends observed from molecular modeling. For the SiO bond, the shared interaction is indicated to increase with decreasing bond length, whereas the AlO bond is of distinctly more polar nature, as is the NaO bond compared to CaO. Also, the observed ranges of the Si and Al in situ electronegativities and their mean electronegativities agree well with both Pauling's values and model calculation results, and statistically significant correlations are obtained which are consistent with trends described for oxide and nitride molecules. Received: 10 May 1999 / Revised, accepted: 14 September 1999     

Elasticity and equation of state of Li<Subscript>2</Subscript>B<Subscript>4</Subscript>O<Subscript>7</Subscript>

A single crystal X-ray diffraction study on lithium tetraborate Li₂B₄O₇ (diomignite, space group I4₁ cd) has been performed under pressure up to 8.3 GPa. No phase transitions were found in the pressure range investigated, and hence the pressure evolution of the unit-cell volume of the I4₁ cd structure has been described using a third-order Birch–Murnaghan equation of state (BM-EoS) with the following parameters: V ₀  = 923.21(6) ų, K ₀  = 45.6(6) GPa, and K′ = 7.3(3). A linearized BM-EoS was fitted to the axial compressibilities resulting in the following parameters a ₀  = 9.4747(3) Å, K _0a  = 73.3(9) GPa, K′ _a  = 5.1(3) and c ₀  = 10.2838(4) Å, K _0c  = 24.6(3) GPa, K′ _c  = 7.5(2) for the a and c axes, respectively. The elastic anisotropy of Li₂B₄O₇ is very large with the zero-pressure compressibility ratio β _0c /β _0a  = 3.0(1). The large elastic anisotropy is consistent with the crystal structure: A three-dimensional arrangement of relatively rigid tetraborate groups [B₄O₇]²⁻ forms channels occupied by lithium along the polar c–axis, and hence compression along the c axis requires the shrinkage of the lithium channels, whereas compression in the a direction depends mainly on the contraction of the most rigid [B₄O₇]²⁻ units. Finally, the isothermal bulk modulus obtained in this work is in general agreement with that derived from ultrasonic (Adachi et al. in Proceedings-IEEE Ultrasonic Symposium, 228–232, 1985; Shorrocks et al. in Proceedings-IEEE Ultrasonic Symposium, 337–340, 1981) and Brillouin scattering measurements (Takagi et al. in Ferroelectrics, 137:337–342, 1992).     

High-pressure phase behavior of MnTiO<Subscript>3</Subscript>: decomposition of perovskite into MnO and MnTi<Subscript>2</Subscript>O<Subscript>5</Subscript>

The phase relations and compression behavior of MnTiO₃ perovskite were examined using a laser-heated diamond-anvil cell, X-ray diffraction, and analytical transmission electron microscopy. The results show that MnTiO₃ perovskite becomes unstable and decomposes into MnO and orthorhombic MnTi₂O₅ phases at above 38 GPa and high temperature. This is the first example of ABO₃ perovskite decomposing into AO + AB₂O₅ phases at high pressure. The compression behavior of volume, axes, and the tilting angle of TiO₆ octahedron of MnTiO₃ perovskite are consistent with those of other A²⁺B⁴⁺O₃ perovskites, although no such decomposition was observed in other perovskites. FeTiO₃ is also known to decompose into two phases, instead of transforming into the CaIrO₃-type post-perovskite phase and we argue that one of the reasons for the peculiar behavior of titanate is the weak covalency of the Ti–O chemical bonds.     

H<Subscript>2</Subscript>O storage capacity of olivine and low-Ca pyroxene from 10 to 13 GPa: consequences for dehydration melting above the transition zone

The onset of hydrous partial melting in the mantle above the transition zone is dictated by the H₂O storage capacity of peridotite, which is defined as the maximum concentration that the solid assemblage can store at P and T without stabilizing a hydrous fluid or melt. H₂O storage capacities of minerals in simple systems do not adequately constrain the peridotite water storage capacity because simpler systems do not account for enhanced hydrous melt stability and reduced H₂O activity facilitated by the additional components of multiply saturated peridotite. In this study, we determine peridotite-saturated olivine and pyroxene water storage capacities at 10–13 GPa and 1,350–1,450°C by employing layered experiments, in which the bottom ~2/3 of the capsule consists of hydrated KLB-1 oxide analog peridotite and the top ~1/3 of the capsule is a nearly monomineralic layer of hydrated Mg# 89.6 olivine. This method facilitates the growth of ~200-μm olivine crystals, as well as accessory low-Ca pyroxenes up to ~50 μm in diameter. The presence of small amounts of hydrous melt ensures that crystalline phases have maximal H₂O contents possible, while in equilibrium with the full peridotite assemblage (melt + ol + pyx + gt). At 12 GPa, olivine and pyroxene water storage capacities decrease from ~1,000 to 650 ppm, and ~1,400 to 1,100 ppm, respectively, as temperature increases from 1,350 to 1,450°C. Combining our results with those from a companion study at 5–8 GPa (Ardia et al., in prep.) at 1,450°C, the olivine water storage capacity increases linearly with increasing pressure and is defined by the relation C_{\textH2 \textO}^\textolivine ( \textppm ) = 57.6( ±16 ) ×P( \textGPa ) - 169( ±18 ). C_{{{\text{H}}_{2} {\text{O}}}}^{\text{olivine}} \left( {\text{ppm}} \right) = 57.6\left( { \pm 16} \right) \times P\left( {\text{GPa}} \right) - 169\left( { \pm 18} \right). Adjustment of this trend for small increases in temperature along the mantle geotherm, combined with experimental determinations of D_{\textH2 \textO}^{\textpyx/olivine} D_{{{\text{H}}_{2} {\text{O}}}}^{\text{pyx/olivine}} from this study and estimates of D_{\textH2 \textO}^{\textgt/\textolivine} D_{{{\text{H}}_{2} {\text{O}}}}^{{{\text{gt}}/{\text{olivine}}}} , allows for estimation of peridotite H₂O storage capacity, which is 440 ± 200 ppm at 400 km. This suggests that MORB source upper mantle, which contains 50–200 ppm bulk H₂O, is not wet enough to incite a global melt layer above the 410-km discontinuity. However, OIB source mantle and residues of subducted slabs, which contain 300–1,000 ppm bulk H₂O, can exceed the peridotite H₂O storage capacity and incite localized hydrous partial melting in the deep upper mantle. Experimentally determined values of D_{\textH2 \textO}^{\textpyx/\textolivine} D_{{{\text{H}}_{2} {\text{O}}}}^{{{\text{pyx}}/{\text{olivine}}}} at 10–13 GPa have a narrow range of 1.35 ± 0.13, meaning that olivine is probably the most important host of H₂O in the deep upper mantle. The increase in hydration of olivine with depth in the upper mantle may have significant influence on viscosity and other transport properties.     

First-principles-based calculations of the CaCO<Subscript>3</Subscript>–MgCO<Subscript>3</Subscript> and CdCO<Subscript>3</Subscript>–MgCO<Subscript>3</Subscript> subsolidus phase diagrams

 Planewave pseudopotential calculations of supercell total energies were used as bases for first-principles calculations of the CaCO₃–MgCO₃ and CdCO₃–MgCO₃ phase diagrams. Calculated phase diagrams are in qualitative to semiquantitative agreement with experiment. Two unobserved phases, Cd₃Mg (CO₃)₄ and CdMg₃(CO₃)₄, are predicted. No new phases are predicted in the CaCO₃–MgCO₃ system, but a low-lying metastable Ca₃Mg(CO₃)₄ state, analogous to the Cd₃Mg(CO₃)₄ phase is predicted. All of the predicted lowest-lying metastable states, except for huntite CaMg₃(CO₃)₄, have dolomite-related structures, i.e. they are layer structures in which A _m B _n cation layers lie perpendicular to the rhombohedral [111] vector. Received: 6 May 2002 / Accepted: 23 October 2002 Acknowledgements This work was partially supported by NSF contract DMR-0080766 and NIST.     

IR calibrations for water determination in olivine,r-GeO<Subscript>2</Subscript>, and SiO<Subscript>2</Subscript> polymorphs

Mineral-specific IR absorption coefficients were calculated for natural and synthetic olivine, SiO₂ polymorphs, and GeO₂ with specific isolated OH point defects using quantitative data from independent techniques such as proton–proton scattering, confocal Raman spectroscopy, and secondary ion mass spectrometry. Moreover, we present a routine to detect OH traces in anisotropic minerals using Raman spectroscopy combined with the “Comparator Technique”. In case of olivine and the SiO₂ system, it turns out that the magnitude of ε for one structure is independent of the type of OH point defect and therewith the peak position (quartz ε = 89,000 ± 15,000  \textl \textmol_{\textH2\textO}^-1 \textcm^-2\text{l}\,\text{mol}_{{\text{H}_2}\text{O}}^{-1}\,\text{cm}^{-2}), but it varies as a function of structure (coesite ε = 214,000 ± 14,000  \textl \textmol_{\textH2\textO}^-1 \textcm^-2\text{l}\,\text{mol}_{{\text{H}_2}\text{O}}^{-1}\,\text{cm}^{-2}; stishovite ε = 485,000 ± 109,000  \textl \textmol_{\textH2\textO}^-1 \textcm^-2\text{l}\,\text{mol}_{{\text{H}_2}\text{O}}^{-1}\,\text{cm}^{-2}). Evaluation of data from this study confirms that not using mineral-specific IR calibrations for the OH quantification in nominally anhydrous minerals leads to inaccurate estimations of OH concentrations, which constitute the basis for modeling the Earth’s deep water cycle.     

Crystal chemistry and OH defect concentrations in spodumene from different granitic pegmatites

Thirty spodumene samples of distinct paragenetic types (primary magmatic, secondary after petalite and hydrothermal) from variety of granitic pegmatites were characterized by electron microprobe, polarized FTIR spectroscopy and Mössbauer spectroscopy. The FTIR spectra of OH (weak sharp pleochroic bands at 3,425, 3,410, 3,395 cm⁻¹ and in the 3,500–3,470 spectral region) are strongly polarized with maximum absorption parallel to n_γ. The majority of OH dipoles are presumably generated by a partial replacement of O2 oxygen atoms with an orientation pointing above the Li vacancy site. The separation of the bands probably resulted from a replacement of the coordinating Al by Fe and Si by Al. Homogeneous spodumene mostly close to its ideal formula LiAlSi₂O₆ shows Fe (0.00–0.10 apfu as Fe³⁺; Fe³⁺ >> Fe²⁺) and Na (0.00–0.04 apfu) as the only minor cations and Fe³⁺Al₋₁ substitution up to 10 mol% of the LiFe³⁺Si₂O₆ component. Hydrogen concentrations (from 0.1 up to <5 ppm H₂O by weight) vary as a function of genetic type with the highest amounts in high-temperature magmatic spodumene. Differences among particular genetic types of spodumene are related to maximum solubility of OH in spodumene structure at given P–T conditions and at actual chemical composition of spodumene. OH defect concentrations in spodumene follow a trend, LT/LP pyroxenes containing lower hydrogen contents compared to HT/HP ones. The hydrogen contents in particular genetic types of spodumene and their decrease with decreasing T and P are consistent with petrologic models of the pegmatite (sub)types formations.     

Vibrational states of H<Subscript>2</Subscript>O in beryl: physical aspects

Single-crystal polarized Raman spectra (3,000–4,000 cm⁻¹ at 3 ≤ T ≤ 300 K) were measured for synthetic alkali-free and natural beryl, Be₂Al₃Si₆O₁₈·xH₂O, to determine the behavior of H₂O molecules of both Type I and Type II in the cavities. At low temperature, the H₂O molecules of Type I displace from the center of cavity and give rise to very weak hydrogen bonding with the host lattice. The H₂O Type I translational motion is characterized by substantial anharmonicity and looks like a motion of “a particle in the box” with a frequency of 6.3 cm⁻¹. Water Type II is characterized by a free rotation with respect to the C ₂ molecule axis, and it makes possible the water nuclear isomers (i.e. ortho- and para-) to be observed at low temperature.