Ab-initio calculations of the proton location in topaz-OH,Al<Subscript>2</Subscript>SiO<Subscript>4</Subscript>(OH)<Subscript>2</Subscript>

The position of hydrogen in the structure of topaz-OH was determined by means of ab-initio quantum-mechanic calculations. Static lattice energy calculations predict the existence of four non-equivalent positions of protons, which are characterized by O4–H1... O1, O4–H2... O2, O4–H3... O3 and O4–H4... O4 hydrogen bonds. The distribution of the protons between positions of local equilibrium is controlled by the proton–proton avoidance rule and the strength of the hydrogen bonds. The most favourable configuration of hydrogen atoms is achieved for adjacent protons, which form O4–H3... O3 and O4–H4... O4 hydrogen bonds, respectively. The thermal excitation of atoms at a temperature of 55 K is large enough for the hydrogen atoms occasionally to change their positions to form O4–H1... O1 and O4–H2... O2 bonds. At ambient pressures and higher temperatures the protons are in a dynamic exchange between the allowed positions of local minima. As a consequence, for nearly room-temperature conditions, the dynamic change between different structural configurations leads to the violation of all possible symmetry elements and with that to space group #E5/E5#1. The flipping of the protons between different sites is achieved by simple rotation of the OH-dipole and does not produce any significant distortion of the framework of topaz, whose symmetry remains that of the space group Pbnm. Therefore, no reduction of symmetry has been observed in former X-ray studies on topaz-OH. Calculated IR absorption spectra of topaz-OH were found to be in good agreement with measured spectra. According to the calculations, the two favourable configurations of protons might correspond to the measured peak splitting within the OH-stretching range. An experimentally observed low-frequency band at 3520 cm^–1 was assigned to the OH-stretching of the O4–H3... O3 bond, while the band at 3600 cm^–1 was attributed to OH-stretching of the O4–H4... O4 hydrogen bond. The broad peak in FAR-IR frequency range at 100–150 cm^–1 is attributed to the stretching of H3... O3 and H4... O4 contacts. The rate of proton exchange at 670 K among different sites was estimated by ab-inito molecular dynamic simulations. The calculations predict that flipping of adjacent protons between O4–H3... O3 and O4–H4... O4 bonds at 670 K occur at a rate of about 1.96 THz.     

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.     

Low-temperature heat capacities of MgAl<Subscript>2</Subscript>O<Subscript>4</Subscript> and spinels of the MgCr<Subscript>2</Subscript>O<Subscript>4</Subscript>–MgAl<Subscript>2</Subscript>O<Subscript>4</Subscript> solid solution

We present results from low-temperature heat capacity measurements of spinels along the solid solution between MgAl₂O₄ and MgCr₂O₄. The data also include new low-temperature heat capacity measurements for MgAl₂O₄ spinel. Heat capacities were measured between 1.5 and 300 K, and thermochemical functions were derived from the results. No heat capacity anomaly was observed for MgAl₂O₄ spinel; however, we observe a low-temperature heat capacity anomaly for Cr-bearing spinels at temperatures below 15 K. From our data we calculate standard entropies (298.15 K) for Mg(Cr,Al)₂O₄ spinels. We suggest a standard entropy for MgAl₂O₄ of 80.9 ± 0.6 J mol⁻¹ K⁻¹. For the solid solution between MgAl₂O₄ and MgCr₂O₄, we observe a linear increase of the standard entropies from 80.9 J mol⁻¹ K⁻¹ for MgAl₂O₄ to 118.3 J mol⁻¹ K⁻¹ for MgCr₂O₄.     

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     

Shock-induced formation of kyanite (Al<Subscript>2</Subscript>SiO<Subscript>5</Subscript>) from sillimanite within a dense metamorphic rock from the Ries crater (Germany)

A dense (~3.34 g cm^–3) garnet–sillimanite-rich metamorphic rock from the suevite breccia of the Ries impact crater was studied by scanning-electron microscopy and Raman microprobe spectroscopy. In the strongly shocked rock clast kyanite was formed from sillimanite under momentary high pressures of natural shock waves. Kyanite aggregates were found as thin (~0.3–2.0 m) seams along grain boundaries between, and fractures within, sillimanite grains. Within these seams kyanite c-axes are oriented perpendicular to original grain boundaries and fractures. In addition, larger (up to 10 m) isolated kyanite grains were rarely found within host sillimanite. Filamentary kyanite aggregates and isolated crystals typically show shrinkage cracks due to volume decrease (~10%). Locally, broad interstices between sillimanite crystals are filled with aluminosilicate glass containing a high volume fraction of sub-micrometer-sized euhedral crystals. The silica-rich glass suggests incongruent melting of sillimanite at local post-shock temperatures significantly higher than 1,300°C. The edges of adjacent sillimanite grains are thermally and chemically altered. The local generation of temperature spikes is attributed to strong shock wave interactions due to very high shock impedance contrasts.     

In situ high-temperature X-ray and neutron powder diffraction study of cation partitioning in synthetic Mg(Fe<Subscript>0.5</Subscript>Al<Subscript>0.5</Subscript>)<Subscript>2</Subscript>O<Subscript>4</Subscript> spinel

The intra-crystalline cation partitioning over T- and M-sites in a synthetic Mg(Fe,Al)₂O₄ spinel sample has been determined as a function of temperature by Rietveld structure refinements from powder diffraction data, combining in situ high-temperature neutron powder diffraction (NPD; POLARIS diffractometer, at ISIS, Rutherford Appleton Laboratory, UK), to determine the Mg and Al occupancy factors, with in situ high-temperature X-ray powder diffraction, to fix the Fe³⁺ distribution. The results obtained agree with a two-stage reaction, in which an initial exchange between Fe³⁺ and Mg, the former leaving and the latter entering tetrahedral sites, is successively followed by a rearrangement involving also Al. The measured cation distribution has then been compared and discussed with that calculated by the Maximum Configuration Entropy principle, for which only NPD patterns have been used. The cation partitioning has finally been interpreted in the light of the configuration model of O’Neill and Navrotsky.     

Pressure-induced structural evolution of pyrite-type SiO<Subscript>2</Subscript>

The equation of state and pressure-induced structural evolution of pyrite-type SiO₂ have been investigated based on synchrotron X-ray diffraction measurements in a diamond-anvil cell. The Rietveld refinement revealed that the oxygen coordinate x of pyrite-type SiO₂ increases with increasing pressure. The SiO₆ coordination polyhedra of pyrite-type SiO₂ is less compressible than the unit cell, and the increase in x induces a rotation of the SiO₆ coordination polyhedra to fill the blank space in the unit cell. Thus, the volume reduction in pyrite-type SiO₂ is achieved mainly by the rotation of the SiO₆ polyhedra, rather than by the compression of the SiO₆ polyhedra. In addition, the increase in x with increasing pressure enhances a distortion of the coordination polyhedra of pyrite-type SiO₂, implying that pyrite-type SiO₂ is not likely to transform into a fluorite-type structure at higher pressures.     

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.     

Thermoelastic properties of chromium oxide Cr<Subscript>2</Subscript>O<Subscript>3</Subscript> (eskolaite) at high pressures and temperatures

A new synchrotron X-ray diffraction study of chromium oxide Cr₂O₃ (eskolaite) with the corundum-type structure has been carried out in a Kawai-type multi-anvil apparatus to pressure of 15 GPa and temperatures of 1873 K. Fitting the Birch–Murnaghan equation of state (EoS) with the present data up to 15 GPa yielded: bulk modulus (K _0,T0), 206 ± 4 GPa; its pressure derivative K′_0,T , 4.4 ± 0.8; (?K _0,T /?T) = ?0.037 ± 0.006 GPa K^?1; a = 2.98 ± 0.14 × 10^?5 K^?1 and b = 0.47 ± 0.28 × 10^?8 K^?2, where α _0,T  = a + bT is the volumetric thermal expansion coefficient. The thermal expansion of Cr₂O₃ was additionally measured at the high-temperature powder diffraction experiment at ambient pressure and α _0,T0 was determined to be 2.95 × 10^?5 K^?1. The results indicate that coefficient of the thermal expansion calculated from the EoS appeared to be high-precision because it is consistent with the data obtained at 1 atm. However, our results contradict α ₀ value suggested by Rigby et al. (Brit Ceram Trans J 45:137–148, 1946) widely used in many physical and geological databases. Fitting the Mie–Grüneisen–Debye EoS with the present ambient and high-pressure data yielded the following parameters: K _0,T0 = 205 ± 3 GPa, K′_0,T  = 4.0, Grüneisen parameter (γ ₀) = 1.42 ± 0.80, q = 1.82 ± 0.56. The thermoelastic parameters indicate that Cr₂O₃ undergoes near isotropic compression at room and high temperatures up to 15 GPa. Cr₂O₃ is shown to be stable in this pressure range and adopts the corundum-type structure. Using obtained thermoelastic parameters, we calculated the reaction boundary of knorringite formation from enstatite and eskolaite. The Clapeyron slope (with \({\text{d}}P/{\text{d}}T = - 0.014\) GPa/K) was found to be consistent with experimental data.     

Stability of (Cs,K)Al<Subscript>4</Subscript>Be<Subscript>5</Subscript>B<Subscript>11</Subscript>O<Subscript>28</Subscript> (londonite) at high pressure and high temperature: a potential neutron absorber material

The stability and the thermo-elastic behaviour of a natural londonite

Structure and phase transition of CaGe<Subscript>2</Subscript>O<Subscript>5</Subscript> revisited

The structure of CaGe₂O₅ between room temperature and 923 K has been determined by X-ray powder diffraction. A continuous phase transition from triclinic C1¯ to monoclinic C2/c symmetry at T_c=714±3 K is observed. The transition is accompanied by a weak heat capacity anomaly. This anomaly and the strain analysis based on the measured lattice parameters indicate a classical second-order phase transition. The order parameter, as measured by the strain component e₂₃, is associated with the displacement of the Ca cation. Electronic structure optimization by density functional methods is used to verify the centric space group of the low-temperature structure of CaGe₂O₅.     

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.     

Adsorption behaviors of Cd^2＋ on Fe2O3/MnO2 and the effects of coexisting ions under alkaline conditions

This study describes the adsorption features of cadmium on Fe2O3 and MnO2 in alkaline saline conditions. The adsorption reached equilibrium in 6 hours under alkaline conditions. The absorption of cadmium on Fe2O3 and MnO2 was consistent with Freundlich absorption isotherms, and the corresponding adsorption capacities were 16.3 and 16.7 mg·g-1, respectively. Moreover, the adsorption quantity of cadmium on Fe2O3 and MnO2 rose with increasing pH from acidic to neutral, and reached the maximum at pH= 9. The coexisting chlorides reduced the adsorption capacity of Fe2O3 and MnO2. The influence intensities of different cations follow the order of CaCl2>>KCl>NaCl. However, the influence of sodium salts on the capacities of Fe2O3 and MnO2 to adsorb cadmium appeared more complicated: the relatively low concentrations of sodium salts could reduce the adsorption capacity; with increasing concentrations of sodium salts, e.g. NaCl and NaNO3. The adsorption capacity decreased continually. Moreover, due to the competition adsorption and precipitation effects, the adsorption capabilities of Na2CO3, NaH2PO4 and Na2HSO4 could also be reduced and cadmium concentrations in the solution were reduced as well.     

Atomistic simulation of the structure and segregation to the (0001) and surfaces of Fe<Subscript>2</Subscript>O<Subscript>3</Subscript>

The segregation of ten isovalent impurities (Al³⁺ Cr³⁺, Eu³⁺, Gd³⁺, Ho³⁺, La³⁺, Lu³⁺, Nd³⁺, Tb³⁺, Y³⁺) to the and the (0001) surfaces of haematite (-Fe₂O₃) have been studied using atomistic simulation where the forces between the atoms are modelled using the Born model of solids. Segregation is found to be energetically favoured in virtually every case. The results for the surface show that the most favourable impurity surface concentration is 33.33%. The (0001) surface has two possible terminations, one terminated by iron atoms and the other by oxygen. No minimum is calculated for the Fe termination of the (0001) surface at low temperatures, but when the effect of raising the temperature is considered, an energy minimum is found, also at 33.33% impurity coverage. In contrast, the O terminated (0001) surface has a minimum in the segregation energy for between 16.67 and 33.33% depending on the cation being considered.