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.     

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     

The crystal structure of alumoklyuchevskite,K<Subscript>3</Subscript>Cu<Subscript>3</Subscript>AlO<Subscript>2</Subscript>(SO<Subscript>4</Subscript>)<Subscript>4</Subscript>

The crystal structure of the unstable mineral alumoklyuchevskite K₃Cu₃AlO₂(SO₄)₄ [monoclinic, I2, a = 18.772(7), b = 4.967(2), c = 18.468(7) Å, β = 101.66(1)°, V = 1686(1) Å] was refined to R ₁ = 0.131 for 2450 unique reflections with F ≥ 4σF _hkl. The structure is based on oxocentered tetrahedrons (OAlCu ₃ ⁷⁺ ) linked into chains via edges. Each chain is surrounded by SO₄ tetrahedrons forming a structural complex. Each complex is elongated along the b axis. This type of crystal structure was also found in other fumarole minerals of the Great Tolbachik Fissure Eruption (GTFE, Kamchatka Peninsula, Russia, 1975–1976), klyuchevskite, K₃Cu₃Fe³⁺O₂(SO₄)₄; and piypite, K₂Cu₂O(SO₄)₂.     

Tilting and distortion of CaSnO<Subscript>3</Subscript> perovskite to 7 GPa determined from single-crystal X-ray diffraction

The structural changes of CaSnO₃, a GdFeO₃-type perovskite, have been investigated to 7 GPa in a diamond-anvil cell at room temperature using single-crystal X-ray diffraction. Significant changes are observed in both the octahedral Sn–O bond lengths and tilt angles between the SnO6 octahedra. The octahedral (SnO6) site shows anisotropic compression and consequently the distortion of SnO6 increases with pressure. Increased pressure also results in a decrease of both of the inter-octahedral angles, Sn–O1–Sn and Sn–O2–Sn, indicating that octahedral tilting increases with increasing pressure, chiefly equivalent to rotation of the SnO6 octahedra about the pseudocubic <001>_p axis. The distortion in the CaO12 and SnO6 sites, along with the octahedral SnO6 tilting, is attributed to the SnO6 site being less compressible than the CaO12 site.Acknowledgments The authors acknowledge with gratitude the financial support for this work from NSF grant EAR-0105864. Ruby pressure measurements were conducted with the Raman system in the Vibrational Spectroscopy Laboratory in the Department of Geosciences at Virginia Tech with the help of Mr. Charles Farley.     

Low-temperature behaviour of ammonium ion in buddingtonite [N(D/H)<Subscript>4</Subscript>AlSi<Subscript>3</Subscript>O<Subscript>8</Subscript>] from neutron powder diffraction

The structural response of buddingtonite [N(D/H)₄AlSi₃O₈] on cooling has been studied by neutron diffraction. Data have been collected from 280 K down to 11 K, and the crystal structure refined using the Rietveld method. Rigid-body constraints were applied to the ammonium ion to explore the structural properties of ammonium in the M-site cavities at low-temperature. Low-temperature saturation is observed for almost all the lattice parameters. From the present in situ low-temperature neutron diffraction studies, there is no strong evidence of orientational order–disorder of the ammonium ions in buddingtonite.     

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.     

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.     

Thermodynamic properties,low-temperature heat-capacity anomalies,and single-crystal X-ray refinement of hydronium jarosite, (H<Subscript>3</Subscript>O)Fe<Subscript>3</Subscript>(SO<Subscript>4</Subscript>)<Subscript>2</Subscript>(OH)<Subscript>6</Subscript>

Crystals of hydronium jarosite were synthesized by hydrothermal treatment of Fe(III)–SO₄ solutions. Single-crystal XRD refinement with R₁=0.0232 for the unique observed reflections (|F_o| > 4_F) and wR₂=0.0451 for all data gave a=7.3559(8) Å, c=17.019(3) Å, V^o=160.11(4) cm³, and fractional positions for all atoms except the H in the H₃O groups. The chemical composition of this sample is described by the formula (H₃O)_0.91Fe_2.91(SO₄)₂[(OH)_5.64(H₂O)_0.18]. The enthalpy of formation (H^o_f) is –3694.5 ± 4.6 kJ mol^–1, calculated from acid (5.0 N HCl) solution calorimetry data for hydronium jarosite, -FeOOH, MgO, H₂O, and -MgSO₄. The entropy at standard temperature and pressure (S^o) is 438.9±0.7 J mol^–1 K^–1, calculated from adiabatic and semi-adiabatic calorimetry data. The heat capacity (C_p) data between 273 and 400 K were fitted to a Maier-Kelley polynomial C_p(T in K)=280.6 + 0.6149T–3199700T^–2. The Gibbs free energy of formation is –3162.2 ± 4.6 kJ mol^–1. Speciation and activity calculations for Fe(III)–SO₄ solutions show that these new thermodynamic data reproduce the results of solubility experiments with hydronium jarosite. A spin-glass freezing transition was manifested as a broad anomaly in the C_p data, and as a broad maximum in the zero-field-cooled magnetic susceptibility data at 16.5 K. Another anomaly in C_p, below 0.7 K, has been tentatively attributed to spin cluster tunneling. A set of thermodynamic values for an ideal composition end member (H₃O)Fe₃(SO₄)₂(OH)₆ was estimated: G^o_f= –3226.4 ± 4.6 kJ mol^–1, H^o_f=–3770.2 ± 4.6 kJ mol^–1, S^o=448.2 ± 0.7 J mol^–1 K^–1, C_p (T in K)=287.2 + 0.6281T–3286000T^–2 (between 273 and 400 K).     

Single-crystal Raman spectra of YAlO<Subscript>3</Subscript> and GdAlO<Subscript>3</Subscript>: comparison to several orthorhombic ABO<Subscript>3</Subscript> perovskites

Room-temperature-polarized single-crystal Raman spectra have been measured for both GdAlO₃ and YAlO₃. Both aluminates crystallize in the orthorhombic (Pbnm) perovskite structure. Of the 24 possible Raman modes in 4 symmetries, 20 and 17 modes were observed for gadolinium and yttrium aluminates, respectively. Comparisons of the Raman spectra of these two aluminates to those of 28 other orthorhombic ABO₃ perovskites revealed remarkably similar spectral patterns, regardless of chemistry or valency of the cations. Closer examination of the effect of mass, valencies, and size of the cations on the Raman spectra versus composition revealed that for the observed modes, the A cation plays the dominant role in determining the Raman shift. In particular, the one to two lowest energy modes in each symmetry are determined by cation mass and valency no matter what the chemistry. For some perovskites with common A cations, higher energy modes were also strikingly similar. In particular, the calcium perovskites had almost all A_g modes at the same energies despite the greatly varying B cations. The second to the lowest mode in A_g and B_1g depended only on A cation mass for all perovskites. The volume plays a minor role throughout but is hard to separate from mass effects because the most massive cations are also the largest. However, if the B-cation is common, for example, aluminates or ferrites, the volume has a minor effect on the higher energy modes. These trends were not observed for all perovskites. Notable exceptions were found if a perovskite is near a phase transition or metastable, as found for three manganites. The effect of increased valency of the A cation from 2–4 to 3–3 perovskites expresses itself as relatively larger Raman shifts for the lowest energy modes. Analog studies of MgSiO₃ perovskites should be undertaken with only 2–4 perovskites. The increased understanding for the mode distributions of perovskites allows for better estimates of their thermodynamic properties through vibrational modeling.     

Magnetic properties and neutron diffraction study of two manganese sulfosalts: monoclinic MnSb<Subscript>2</Subscript>S<Subscript>4</Subscript> and benavidesite (MnPb<Subscript>4</Subscript>Sb<Subscript>6</Subscript>S<Subscript>14</Subscript>)

Mn²⁺Sb₂S₄, a monoclinic dimorph of clerite, and benavidesite (Mn²⁺Pb₄Sb₆S₁₄) show well-individualized single chains of manganese atoms in octahedral coordination. Their magnetic structures are presented and compared with those of iron derivatives, berthierite (Fe²⁺Sb₂S₄) and jamesonite (Fe²⁺Pb₄Sb₆S₁₄). Within chains, interactions are antiferromagnetic. Like berthierite, MnSb₂S₄ shows a spiral magnetic structure with an incommensurate 1D propagation vector [0, 0.369, 0], unchanged with temperature. In berthierite, the interactions between identical chains are antiferromagnetic, whereas in MnSb₂S₄ interactions between chains are ferromagnetic along c-axis. Below 6 K, jamesonite and benavidesite have commensurate magnetic structures with the same propagation vector [0.5, 0, 0]: jamesonite is a canted ferromagnet and iron magnetic moments are mainly oriented along the a-axis, whereas for benavidesite, no angle of canting is detected, and manganese magnetic moments are oriented along b-axis. Below 30 K, for both compounds, one-dimensional magnetic ordering or correlations are visible in the neutron diagrams and persist down to 1.4 K.     

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.     

Solid solution behaviour of CaSiO<Subscript>3</Subscript> and MgSiO<Subscript>3</Subscript> perovskites

Using density functional simulations within the generalized gradient approximation and projector-augmented wave method together with thermodynamic modelling, the reciprocal solubilities of MgSiO₃ and CaSiO₃ perovskites were calculated for pressures and temperatures of the Earth’s lower mantle from 25 to 100 GPa and 0 to 6,000 K, respectively. The solubility of Ca in MgSiO₃ at conditions along a mantle adiabat is found to be less than 0.02 atoms per formula unit. The solubility of Mg in CaSiO₃ is even lower, and most important, the extent of solid solution decreases with pressure. To dissolve CaSiO₃ perovskite completely in MgSiO₃ perovskite, a solubility of 7.8 or 2.3 mol% would be necessary for a fertile pyrolytic or depleted harzburgitic mantle, respectively. Thus, for any reasonable geotherm, two separate perovskites will be present in fertile mantle, suggesting that Ca-perovskite will be residual to low degree melting throughout the entire mantle. At the solidus, CaSiO₃ perovskite might completely dissolve in MgSiO₃ perovskite only in a depleted mantle with <1.25 wt% CaO. These implications may be modified if Ca solubility in MgSiO₃ is increased by other major mantle constituents such as Fe and Al.     

Elasticity of CaIrO<Subscript>3</Subscript> with perovskite and post-perovskite structure

Compression behaviors of CaIrO₃ with perovskite (Pv) and post-perovskite (pPv) structures have been investigated up to 31.0(1.0) and 35.3(1) GPa at room temperature, respectively, in a diamond-anvil cell with hydrostatic pressure media. CaIrO₃ Pv and pPv phases were compressed with the axial compressibility of β _a > β _c > β _b and β _b > β _a > β _c, respectively and no phase transition was observed in both phases up to the highest pressure in the present study. The order of axial compressibility for pPv phase is consistent with the crystallographic consideration for layer structured materials and previous experimental results. On the other hand, Pv phase shows anomalous compression behavior in b axis, which exhibit constant or slightly expanded above 13 GPa, although the applied pressure remained hydrostatic. Volume difference between Pv and pPv phases was gradually decreased with increasing pressure and this is consistent with the results of theoretical study based on the ab initio calculation. Present results, combined with theoretical study, suggest that these complicate compression behaviors in CaIrO₃ under high pressure might be caused by the partially filled electron of Ir⁴⁺. Special attention must be paid in case of using CaIrO₃ as analog materials to MgSiO₃, although CaIrO₃ exhibits interesting physical properties under high pressure.     

Equation of state and crystal structure of Sb<Subscript>2</Subscript>S<Subscript>3</Subscript> between 0 and 10 GPa

High-precision unit-cell volume data of stibnite, collected in the pressure range of 0–10 GPa, was used for fitting a third-order Birch–Murnaghan equation of state. The zero-pressure volume, bulk modulus and its pressure derivative were found to be 487.73(6) ų, 26.91(14) GPa and 7.9(1), respectively. A series of X-ray intensity data was collected in the same pressure range using a CCD-equipped Bruker diffractometer. The high-pressure structures were all refined to R1(|F₀|>4) values of approximately 0.03. Crystal-chemical parameters as polyhedron volume, centroid and eccentricity were calculated for the seven coordinated cation positions using the software IVTON. The cation eccentricity appears to be a very useful tool for quantification of the lone electron pair activity. U₂S₃, Dy₂S₃ and Nd₂Te₃ are all isostructural with stibnite, but the cations in these materials have no lone electron pair. Their eccentricity is much smaller than that of Sb, and close to zero. This confirms that the stibnite structure type alone does not force eccentricity upon the cations involved and it is the lone electron pairs of Sb that generate the eccentricity of cation positions in the structures of stibnite. At increasing pressure the eccentricity of Sb is decreasing. It is therefore reasonable to conclude that the lone electron pair activity is decreasing with increasing pressure.     

Elastic behavior of vanadinite,Pb<Subscript>10</Subscript>(VO<Subscript>4</Subscript>)<Subscript>6</Subscript>Cl<Subscript>2</Subscript>, a microporous non-zeolitic mineral

The high-pressure behavior of a vanadinite (Pb₁₀(VO₄)₆Cl₂, a = b = 10.3254(5), c = 7.3450(4) Å, space group P6₃/m), a natural microporous mineral, has been investigated using in-situ HP-synchrotron X-ray powder diffraction up to 7.67 GPa with a diamond anvil cell under hydrostatic conditions. No phase transition has been observed within the pressure range investigated. Axial and volume isothermal Equations of State (EoS) of vanadinite were determined. Fitting the P–V data with a third-order Birch-Murnaghan (BM) EoS, using the data weighted by the uncertainties in P and V, we obtained: V ₀ = 681(1) Å³, K ₀ = 41(5) GPa, and K′ = 12.5(2.5). The evolution of the lattice constants with P shows a strong anisotropic compression pattern. The axial bulk moduli were calculated with a third-order “linearized” BM-EoS. The EoS parameters are: a ₀ = 10.3302(2) Å, K ₀(a) = 35(2) GPa and K′(a) = 10(1) for the a-axis; c ₀ = 7.3520(3) Å, K ₀(c) = 98(4) GPa, and K′(c) = 9(2) for the c-axis (K ₀(a):K ₀(c) = 1:2.80). Axial and volume Eulerian-finite strain (fe) at different normalized stress (Fe) were calculated. The weighted linear regression through the data points yields the following intercept values: Fe _a (0) = 35(2) GPa for the a-axis, Fe _c (0) = 98(4) GPa for the c-axis and Fe _V (0) = 45(2) GPa for the unit-cell volume. The slope of the regression lines gives rise to K′ values of 10(1) for the a-axis, 9(2) for the c-axis and 11(1) for the unit cell-volume. A comparison between the HP-elastic response of vanadinite and the iso-structural apatite is carried out. The possible reasons of the elastic anisotropy are discussed.     

Peierls dislocation modelling in perovskite (CaTiO<Subscript>3</Subscript>): comparison with tausonite (SrTiO<Subscript>3</Subscript>) and MgSiO<Subscript>3</Subscript> perovskite

We present here a numerical modelling study of dislocations in perovskite CaTiO₃. The dislocation core structures and properties are calculated through the Peierls–Nabarro model using the generalized stacking fault (GSF) results as a starting model. The GSF are determined from first-principles calculations using the VASP code. The dislocation properties such as collinear, planar core spreading and Peierls stresses are determined for the following slip systems: [100](010), [100](001), [010](100), [010](001), [001](100), [001](010), and All dislocations exhibit lattice friction, but glide appears to be easier for [100](010) and [010](100). [001](010) and [001](100) exhibit collinear dissociation. Comparing Peierls stresses among tausonite (SrTiO₃), perovskite (CaTiO₃) and MgSiO₃ perovskite demonstrates the strong influence of orthorhombic distortions on lattice friction. However, and despite some quantitative differences, CaTiO₃ appears to be a satisfactory analogue material for MgSiO₃ perovskite as far as dislocation glide is concerned.     

Hydroxylborite,Mg<Subscript>3</Subscript>(BO<Subscript>3</Subscript>)(OH)<Subscript>3</Subscript>, a new mineral species and isomorphous series fluoborite-hydroxylborite

Hydroxylborite, a new mineral species, an analogue of fluoborite with OH > F, has been found at the Titovsky deposit (57°41′N, 125°22′E), the Chersky Range, Dogdo Basin, Sakha-Yakutia Republic, Russia. Prismatic crystals of the new mineral are dominated by the {10\(\overline 1 \)0} faces without distinct end forms and reach (1?1.5) × (0.1?0.2) mm in size. Radial aggregates of such crystals occur in the mineralized marble adjacent to the boron ore (suanite-kotoite-ludwigite). Calcite, dolomite, Mg-rich ludwigite, kotoite, szaibelyite, clinohumite, magnetite, serpentine, and chlorite are associated minerals. Hydroxylborite is transparent colorless, with a white streak and vitreous luster. The new mineral is brittle. The Mohs’ hardness is 3.5. The cleavage is imperfect on {0001}. The density measured with equilibration in heavy liquids is 2.89(1) g/cm³; the calculated density is 2.872 g/cm³. The wave numbers of the absorption bands in the IR spectrum of hydroxylborite are (cm^?1; sh is shoulder): 3668, 1233, 824, 742, 630sh, 555sh, 450sh, and 407. The new mineral is optically uniaxial, negative, ω = 1.566(1), and ε = 1.531(1). The chemical composition (electron microprobe, H₂O measured with the Penfield method, wt %) is 18.43 B₂O₃, 65.71 MgO, 10.23 F, 9.73 H₂O, 4.31-O = F₂, where the total is 99.79. The empirical formula calculated on the basis of 6 anions pfu is as follows: Mg_3.03B_0.98[(OH)_2.00F_1.00]O_3.00. Hydroxylborite is hexagonal, and the space group is P6₃/m. The unit-cell dimensions are: a = 8.912(8) Å, c = 3.112(4) Å, V = 214.05(26) ų, and Z = 2. The strongest reflections in the X-ray powder pattern [d, Å (I, %)(hkil)] are: 7.69(52)(01\(\overline 1 \)0), 4.45(82)(11\(\overline 2 \)0), 2.573(65)(03\(\overline 3 \)0), 2.422(100)(02\(\overline 2 \)1), and 2.128(60)(12\(\overline 3 \)1). The compatibility index 1 ? (K _p/K _c) is 0.038 (excellent) for the calculated density and 0.044 (good) for the measured density. The type material of hydroxylborite is deposited in the Fersman Mineralogical Museum, Russian Academy of Sciences, Moscow (inventory number 91968) and the Geological Museum of the All-Russia Institute of Mineral Resources, Moscow (inventory number M-1663).     

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.     

An infrared and <Superscript>1</Superscript>H MAS NMR investigation of strong hydrogen bonding in ussingite,Na<Subscript>2</Subscript>AlSi<Subscript>3</Subscript>O<Subscript>8</Subscript>(OH)

The mineral ussingite, Na₂AlSi₃O₈(OH), an interrupted tectosilicate, has strong hydrogen bonding between OH and the other nonbridging oxygen atom in the structure. Infrared spectra contain a strongly polarized, very broad OH-stretching band with an ill-defined maximum between 1500 and 1800 cm^–1, and a possible OH librational bending mode at 1295 cm^–1. The IR spectra confirm the orientation of the OH vector within the triclinic unit cell as determined from X-ray refinement (Rossi et al. 1974). There are three distinct bands in the ¹H NMR spectrum of ussingite: a predominant band at 13.5 ppm (TMS) representing 90% of the structural hydrogen, a second band at 15.9 ppm corresponding to 8% of the protons, and a third band at 11.0 ppm accounting for the remaining 2% of structural hydrogen. From the correlation between hydrogen bond length and ¹H NMR chemical shift (Sternberg and Brunner 1994), the predominant hydrogen bond length (H...O) was calculated to be 1.49 Å, in comparison to the hydrogen bond length determined from X-ray refinement (1.54 Å). The population of protons at 15.9 ppm is consistent with 5–8% Al–Si disorder. Although the ussingite crystal structure and composition are similar to those of low albite, the bonding environment of OH in low albite and other feldspars, as characterized through IR and ¹H NMR, is fundamentally different from the strong hydrogen bonding found in ussingite.