Crystal chemistry of natural and synthetic lead oxyhalides. Part I. Crystal structure of Pb<Subscript>13</Subscript>O<Subscript>10</Subscript>Cl<Subscript>6</Subscript>

Crystals of lead oxychloride Pb₁₃O₁₀Cl₆ have been synthesized on the basis of high-temperature solid-state reactions. The Pb₁₃O₁₀Cl₆ structure was studied using X-ray single-crystal diffraction analysis. The compound is monoclinic, and the space group is C2/c; the unit-cell dimensions are a = 16.1699(14), b = 7.0086(6), c = 23.578(2) Å, β = 97.75°, and V = 2647.6(4) ų. The structure has been solved by direct methods and refined to R ₁ = 0.0505 for 2671 observed unique reflections. The structure is a 3D framework consisting of OPb₄ tetrahedrons. Chlorine atoms are located in the framework channels. The structure contains seven symmetrically independent Pb atoms, which are coordinated by 2 to 4 O^2? and 2 to 4 Cl^? anions. The synthesized compound is compared with other natural and synthetic lead oxyhalides.     

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₄.     

Crystal chemistry of selenates with mineral-like structures. III. Heteropolyhedral chains in the crystal structure of [Mg(H<Subscript>2</Subscript>O)<Subscript>4</Subscript>(SeO<Subscript>4</Subscript>)]<Subscript>2</Subscript>(H<Subscript>2</Subscript>O)

The crystal structure of a new compound [Mg(H₂O)₄(SeO₄)]₂(H₂O) (monoclinic, P2 ₁/a, a = 7.2549(12), b = 20.059(5), c = 10.3934(17) Å, β = 101.989(13), V = 1479.5(5) ų) has been solved by direct methods and refined to R ₁ = 0.059 for 2577 observed reflections with |F _hkl | ≥ 4σ|F _hkl |. The structure consists of [Mg(H₂O)₄(SeO₄)]⁰ chains formed by alternating corner-sharing Mg octahedrons and (SeO₄)^2? tetrahedrons. O atoms of Mg octahedrons that are shared with selenate tetrahedrons are in a trans orientation. The heteropoly-hedral octahedral-tetrahedral chains are parallel to the c axis and undulate within the (010) plane. The adjacent chains are linked by hydrogen bonds involving H₂O molecules not bound with M²⁺ cations.     

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₄)₂.     

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     

Pressure induced phase transition in Pb<Subscript>6</Subscript>Bi<Subscript>2</Subscript>S<Subscript>9</Subscript>

The crystal structure of Pb₆Bi₂S₉ is investigated at pressures between 0 and 5.6 GPa with X-ray diffraction on single-crystals. The pressure is applied using diamond anvil cells. Heyrovskyite (Bbmm, a = 13.719(4) Å, b = 31.393(9) Å, c = 4.1319(10) Å, Z = 4) is the stable phase of Pb₆Bi₂S₉ at ambient conditions and is built from distorted moduli of PbS-archetype structure with a low stereochemical activity of the Pb²⁺ and Bi³⁺ lone electron pairs. Heyrovskyite is stable until at least 3.9 GPa and a first-order phase transition occurs between 3.9 and 4.8 GPa. A single-crystal is retained after the reversible phase transition despite an anisotropic contraction of the unit cell and a volume decrease of 4.2%. The crystal structure of the high pressure phase, β-Pb₆Bi₂S₉, is solved in Pna2 ₁ (a = 25.302(7) Å, b = 30.819(9) Å, c = 4.0640(13) Å, Z = 8) from synchrotron data at 5.06 GPa. This structure consists of two types of moduli with SnS/TlI-archetype structure in which the Pb and Bi lone pairs are strongly expressed. The mechanism of the phase transition is described in detail and the results are compared to the closely related phase transition in Pb₃Bi₂S₆ (lillianite).     

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.     

Minerals at high pressure. Mechanics of compression from quantum mechanical calculations in a case study: the beryl (Al<Subscript>4</Subscript>Be<Subscript>6</Subscript>Si<Subscript>12</Subscript>O<Subscript>36</Subscript>)

The Bader topological analysis has been applied to ab initio computed electron densities of beryl, in order to clarify its mechanism of compression. Full structural optimization and total energy (E) calculations were performed at different cell volumes (V _c). The pressure at each volume and the equation of state were estimated from the first and second derivatives of the resultant E(V _c) curve. The total (negative) potential energy of the crystal, sum of both attractive and repulsive electrostatic terms, was found to systematically decrease (i.e., it moved to more negative values) up to the highest pressure considered (28.4 GPa), indicating that interelectronic and internuclear repulsions are not the only terms controlling the compressibility, at least in the pressure range investigated. Electronic kinetic energy increases as the cell volume is reduced, leading to a parallel increase of the total energy. Both structure at equilibrium and compressibility are therefore due to the balance between the opposing kinetic and potential energy terms. The Bader theory has been used to identify the topological atoms within the structure and to calculate their properties, with particular attention to the forces driving the structural relaxation at high pressure. On a qualitative basis, the obtained results are expected to be transferable to the discussion of compressibility of other mineral phases.     

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).     

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.     

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.     

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.     

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.     

Calcite and barite precipitation in CaCO<Subscript>3</Subscript>-BaSO<Subscript>4</Subscript>-NaCl and BaSO<Subscript>4</Subscript>-NaCl-CaCl<Subscript>2</Subscript> aqueous systems: kinetic and microstructural study

During the production of hydrocarbons from subterranean reservoirs, scaling with calcium carbonate and barium sulfate causes flux decline and dangerous problems in production facilities. This work is intended to study the effect of calcium ions on the precipitation of barium sulfate (barite); then, the effect of the formed barite on calcium carbonate crystallization. The conductometric and pH methods were used to follow the progress of the precipitation reaction in aqueous medium. The obtained precipitates were characterized by FTIR, RAMAN, SEM, and XRD. It was shown that Ca²⁺ in the reaction media does not affect the microstructure of barite even for higher calcium–barium molar ratio. It influences the precipitation kinetics and the solubility of barite by the formation of CaSO₄° ion pairing as a predominant role of complex formation (CaSO4) and the increase of the ionic strength. In Ca(HCO₃)₂-BaSO₄-NaCl aqueous system, experiments have showed that added or formed barite in the reaction media accelerates calcite precipitation. No effect on the microstructure of heterogeneous formed calcite which remain calcite shape. However the presence of carbonate ions affects slightly the microstructure of barite.     

Compressibility of strontium orthophosphate Sr<Subscript>3</Subscript>(PO<Subscript>4</Subscript>)<Subscript>2</Subscript> at high pressure

High pressure in situ synchrotron X-ray diffraction experiment of strontium orthophosphate Sr₃(PO₄)₂ has been carried out to 20.0 GPa at room temperature using multianvil apparatus. Fitting a third-order Birch–Murnaghan equation of state to the P–V data yields a volume of V ₀ = 498.0 ± 0.1 Å³, an isothermal bulk modulus of K _T  = 89.5 ± 1.7 GPa, and first pressure derivative of K _T ′ = 6.57 ± 0.34. If K _T ′ is fixed at 4, K _T is obtained as 104.4 ± 1.2 GPa. Analysis of axial compressible modulus shows that the a-axis (K _a  = 79.6 ± 3.2 GPa) is more compressible than the c-axis (K _c  = 116.4 ± 4.3 GPa). Based on the high pressure Raman spectroscopic results, the mode Grüneisen parameters are determined and the average mode Grüneisen parameter of PO₄ vibrations of Sr₃(PO₄)₂ is calculated to be 0.30(2).     

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.     

Temperature-jump induced cation exchange kinetics in (Co<Subscript>0.1</Subscript>Mg<Subscript>0.9</Subscript>)<Subscript>2</Subscript>SiO<Subscript>4</Subscript> olivine: an in situ optical spectroscopic study

In the olivine crystal structure, cations are distributed over two inequivalent octahedral sites, M1 and M2. Kinetics of cation exchange between the two octahedral sites in (Co_0.1Mg_0.9)₂SiO₄ single crystal have been studied in the temperature range from 600 to 800°C by monitoring the time evolution of the absorbance of Co²⁺ ions in M1 or M2 sites using optical spectroscopy after rapid temperature jumps. It was found from such temperature-jump induced relaxation experiments that with increasing temperature the absorbance of Co²⁺ ions in the M1 site decreases while that in the M2 site increases. This indicates a tendency of Co²⁺ cations to populate the M2 site with increasing temperatures and vice versa. The experimental relaxation data can be modeled using a triple exponential equation based on theoretical analysis. Activation energies of 221 ± 4 and 213 ± 10 kJ/mol were derived from relaxation experiments on the M2 site and M1 site, respectively, for the cation exchange processes in (Co_0.1Mg_0.9)₂SiO₄ olivine. Implications for cation diffusion at low temperatures are discussed.