Retrieval of stratospheric O3 and NO2 vertical profiles using zenith scattered light observations

Daily zenith scattered light intensity observations were carried out in the morning twilight hours using home-made UV-visible spectrometer over the tropical station Pune (18‡31′, 73‡51′) for the years 2000–2003. These observations are obtained in the spectral range 462–498 nm for the solar zenith angles (SZAs) varying from 87‡ to 91.5‡. An algorithm has been developed to retrieve vertical profiles of ozone (O₃) and nitrogen dioxide (NO₂) from ground-based measurements using the Chahine iteration method. This retrieval method has been checked using measured and recalculated slant column densities (SCDs) and they are found to be well matching. O₃ and NO₂ vertical profiles have been retrieved using a set of their air mass factors (AMFs) and SCDs measured over a range of 87–91.5‡ SZA during the morning. The vertical profiles obtained by this method are compared with Umkehr profiles and ozonesondes and they are found to be in good agreement. The bulk of the column density is found near layer 20–25 km. Daily total column densities (TCDs) of O₃ and NO₂ along with their stratospheric and tropospheric counterparts are derived using their vertical profiles for the period 2000–2003. The total column, stratospheric column and tropospheric column amounts of both trace gases are found to be maximum in summer and minimum in the winter season. Increasing trend is found in column density of NO₂ in stratospheric, tropospheric and surface layers, but no trend is observed in O₃ columns for above layers during the period 2000–2003     

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.     

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     

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

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

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.     

Structure and chemistry of (111) twin boundaries in MgAl<Subscript>2</Subscript>O<Subscript>4</Subscript> spinel crystals from Mogok

The atomic scale structure and chemistry of (111) twins in MgAl₂O₄ spinel crystals from the Pinpyit locality near Mogok (Myanmar, formerly Burma) were analysed using complementary methods of transmission electron microscopy (TEM). To obtain a three-dimensional information on the atomic structure, the twin boundaries were investigated in crystallographic projections and Using conventional electron diffraction and high-resolution TEM (HRTEM) analysis we have shown that (111) twins in spinel can be crystallographically described by 180° rotation of the oxygen sublattice normal to the twin composition plane. This operation generates a local hcp stacking in otherwise ccp lattice and maintains a regular sequence of kagome and mixed layers. In addition to rotation, no other translations are present in (111) twins in these spinel crystals. Chemical analysis of the twin boundary was performed by energy-dispersive X-ray spectroscopy (EDS) using a variable beam diameter (VBD) technique, which is perfectly suited for analysing chemical composition of twin boundaries on a sub-nm scale. The VBD/EDS measurements indicated that (111) twin boundary in spinel is Mg-deficient. Quantitative analyses of HRTEM (phase contrast) and HAADF-STEM (Z-contrast) images of (111) twin boundary have confirmed that Mg²⁺ ions are replaced with Be²⁺ ions in boundary tetrahedral sites. The Be-rich twin boundary structure is closely related to BeAl₂O₄ (chrysoberyl) and BeMg₃Al₈O₁₆ (taaffeite) group of intermediate polysomatic minerals. Based on these results, we conclude that the formation of (111) twins in spinel is a preparatory stage of polytype/polysome formation (taaffeite) and is a result of thermodynamically favourable formation of hcp stacking due to Be incorporation on the {111} planes of the spinel structure in the nucleation stage of crystal growth. The twin structure grows as long as the surrounding geochemical conditions allow its formation. The incorporation of Be induces a 2D-anisotropy and exaggerated growth of the crystal along the (111) twin boundary.     

Batisivite,V<Subscript>8</Subscript>Ti<Subscript>6</Subscript>[Ba(Si<Subscript>2</Subscript>O)]O<Subscript>28</Subscript>, a new mineral species from the derbylite group

Batisivite has been found as an accessory mineral in the Cr-V-bearing quartz-diopside metamorphic rocks of the Slyudyanka Complex in the southern Baikal region, Russia. A new mineral was named after the major cations in its ideal formula (Ba, Ti, Si, V). Associated minerals are quartz, Cr-V-bearing diopside and tremolite; calcite; schreyerite; berdesinskiite; ankangite; V-bearing titanite; minerals of the chromite-coulsonite, eskolaite-karelianite, dravite-vanadiumdravite, and chernykhite-roscoelite series; uraninite; Cr-bearing goldmanite; albite; barite; zircon; and unnamed U-Ti-V-Cr phases. Batisivite occurs as anhedral grains up to 0.15–0.20 mm in size, without visible cleavage and parting. The new mineral is brittle, with conchoidal fracture. Observed by the naked eye, the mineral is black and opaque, with a black streak and resinous luster. Batisivite is white in reflected light. The microhardness (VHN) is 1220–1470 kg/mm² (load is 30 g), the mean value is 1330 kg/mm². The Mohs hardness is near 7. The calculated density is 4.62 g/cm³. The new mineral is weakly anisotropic and bireflected. The measured values of reflectance are as follows (λ, nm—R _max ^′ /R _min ^′ ): 440—17.5/17.0; 460—17.3/16.7; 480—17.1/16.5; 500—17.2/16.6; 520—17.3/16.7; 540—17.4/16.8; 560—17.5/16.8; 580—17.6/16.9; 600—17.7/17.1; 620—17.7/17.1; 640—17.8/17.1; 660—17.9/17.2; 680—18.0/17.3; 700—18.1/17.4. Batisivite is triclinic, space group P \(\overline 1\); the unit-cell dimensions are: a = 7.521(1) Å, b = 7.643(1) Å, c = 9.572(1) Å, α = 110.20°(1), β = 103.34°(1), γ = 98.28°(1), V = 487.14(7) ų, Z = 1. The strongest reflections in the X-ray powder diffraction pattern [d, Å (I, %)(hkl)] are: 3.09(8)(12\(\overline 2\)); 2.84, 2.85(10)(021, 120); 2.64(8)(21\(\overline 3\)); 2.12(8)(31\(\overline 3\)); 1.785(8)(32\(\overline 4\)), 1.581(10)(24\(\overline 2\)); 1.432, 1.433(10)(322, 124). The chemical composition (electron microprobe, average of 237 point analyses, wt %) is: 0.26 Nb₂O₅, 6.16 SiO₂, 31.76 TiO₂, 1.81 Al₂O₃, 8.20 VO₂, 26.27 V₂O₃, 12.29 Cr₂O₃, 1.48 Fe₂O₃, 0.08 MgO, 11.42 BaO; the total is 99.73. The VO₂/V₂O₃ ratio has been calculated. The simplified empirical formula is (V _4.8 ³⁺ Cr_2.2V _0.7 ⁴⁺ Fe_0.3)_8.0(Ti_5.4V _0.6 ⁴⁺ )_6.0[Ba(Si_1.4Al_0.5O_0.9)]O₂₈. An alternative to the title formula could be a variety (with the diorthogroup Si₂O₇) V₈Ti₆[Ba(Si₂O₇)]O₂₂. Batisivite probably pertains to the V ₈ ³⁺ Ti ₆ ⁴⁺ [Ba(Si₂O)]O₂₈-Cr ₈ ³⁺ Ti ₆ ⁴⁺ [Ba(Si₂O)]O₂₈ solid solution series. The type material of batisivite has been deposited in the Fersman Mineralogical Museum, Russian Academy of Sciences, Moscow.     

The H<Subscript>2</Subscript>O solubility of alkali basaltic melts: an experimental study

Experiments were conducted to determine the water solubility of alkali basalts from Etna, Stromboli and Vesuvius volcanoes, Italy. The basaltic melts were equilibrated at 1,200°C with pure water, under oxidized conditions, and at pressures ranging from 163 to 3,842 bars. Our results show that at pressures above 1 kbar, alkali basalts dissolve more water than typical mid-ocean ridge basalts (MORB). Combination of our data with those from previous studies allows the following simple empirical model for the water solubility of basalts of varying alkalinity and fO₂ to be derived: \textH ₂ \textO( \textwt% ) = \text H ₂ \textO_\textMORB ( \textwt% ) + ( 5.84 ×10 ^{- 5} *\textP - 2.29 ×10 ^{- 2} ) ×( \textNa₂ \textO + \textK₂ \textO )( \textwt% ) + 4.67 ×10 ^{- 2} ×\Updelta \textNNO - 2.29 ×10 ^{- 1} {\text{H}}_{ 2} {\text{O}}\left( {{\text{wt}}\% } \right) = {\text{ H}}_{ 2} {\text{O}}_{\text{MORB}} \left( {{\text{wt}}\% } \right) + \left( {5.84 \times 10^{ - 5} *{\text{P}} - 2.29 \times 10^{ - 2} } \right) \times \left( {{\text{Na}}_{2} {\text{O}} + {\text{K}}_{2} {\text{O}}} \right)\left( {{\text{wt}}\% } \right) + 4.67 \times 10^{ - 2} \times \Updelta {\text{NNO}} - 2.29 \times 10^{ - 1} where H₂O_MORB is the water solubility at the calculated P, using the model of Dixon et al. (1995). This equation reproduces the existing database on water solubilities in basaltic melts to within 5%. Interpretation of the speciation data in the context of the glass transition theory shows that water speciation in basalt melts is severely modified during quench. At magmatic temperatures, more than 90% of dissolved water forms hydroxyl groups at all water contents, whilst in natural or synthetic glasses, the amount of molecular water is much larger. A regular solution model with an explicit temperature dependence reproduces well-observed water species. Derivation of the partial molar volume of molecular water using standard thermodynamic considerations yields values close to previous findings if room temperature water species are used. When high temperature species proportions are used, a negative partial molar volume is obtained for molecular water. Calculation of the partial molar volume of total water using H₂O solubility data on basaltic melts at pressures above 1 kbar yields a value of 19 cm³/mol in reasonable agreement with estimates obtained from density measurements.     

Calorimetric study on high-pressure transitions in KAlSi<Subscript>3</Subscript>O<Subscript>8</Subscript>

KAlSi₃O₈ sanidine dissociates into a mixture of K₂Si₄O₉ wadeite, Al₂SiO₅ kyanite and SiO₂ coesite, which further recombine into KAlSi₃O₈ hollandite with increasing pressure. Enthalpies of KAlSi₃O₈ sanidine and hollandite, K₂Si₄O₉ wadeite and Al₂SiO₅ kyanite were measured by high-temperature solution calorimetry. Using the data, enthalpies of transitions at 298 K were obtained as 65.1 ± 7.4 kJ mol^–1 for sanidine wadeite + kyanite + coesite and 99.3 ± 3.6 kJ mol^–1 for wadeite + kyanite + coesite hollandite. The isobaric heat capacity of KAlSi₃O₈ hollandite was measured at 160–700 K by differential scanning calorimetry, and was also calculated using the Kieffer model. Combination of both the results yielded a heat-capacity equation of KAlSi₃O₈ hollandite above 298 K as C_p=3.896 × 10²–1.823 × 10³T^–0.5–1.293 × 10⁷T^–2+1.631 × 10⁹T^–3 (C_p in J mol^–1 K^–1, T in K). The equilibrium transition boundaries were calculated using these new data on the transition enthalpies and heat capacity. The calculated transition boundaries are in general agreement with the phase relations experimentally determined previously. The calculated boundary for wadeite + kyanite + coesite hollandite intersects with the coesite–stishovite transition boundary, resulting in a stability field of the assemblage of wadeite + kyanite + stishovite below about 1273 K at about 8 GPa. Some phase–equilibrium experiments in the present study confirmed that sanidine transforms directly to wadeite + kyanite + coesite at 1373 K at about 6.3 GPa, without an intervening stability field of KAlSiO₄ kalsilite + coesite which was previously suggested. The transition boundaries in KAlSi₃O₈ determined in this study put some constraints on the stability range of KAlSi₃O₈ hollandite in the mantle and that of sanidine inclusions in kimberlitic diamonds.     

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.     

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.     

Experimental multipole-refined and theoretical charge density study of LiGaSi<Subscript>2</Subscript>O<Subscript>6</Subscript> clinopyroxene at ambient conditions

The synthetic LiGaSi₂O₆ clinopyroxene is monoclinic C2/c at room-T. Its experimental electron density, ρ(r), has been derived starting from accurate room-T single-crystal diffraction data. Topological analysis confirms an intermediate ionic-covalent character for Si–O bonding, as found by previous electron-density studies on other silicates such as diopside, coesite and stishovite. The non-bridging Si–O bonds have more covalent character than the bridging ones. The Ga–O bonds have different bonding characters, the Ga–O2 bond being more covalent than the two Ga–O1 bonds. Li–O bonds are classified as pure closed-shell ionic interactions. Similar to spodumene (LiAlSi₂O₆), Li has sixfold coordination, but the bond critical points associated to the two longest bonds are characterized by very low electron density values. Similar to what previously found in spodumene and diopside, O···O interactions were detected from the topological analysis of ρ(r), and indicate a cooperative interaction among the lone pairs of neighbouring oxygen atoms. In particular, this kind of interaction has been obtained for the O1···O1 edge shared between two Ga octahedra. Integration over the atomic basins gives net charges of −1.39(10), 2.82(10), 1.91(10) and 0.82(8) e for O (averaged), Si, Ga and Li atoms, respectively. Periodic Hartree–Fock and DFT calculations confirm the results obtained by multipole refinement of the experimental data. Moreover, the theoretical topological properties of the electron density distribution on the Si₂O₆ group are very similar to those calculated for spodumene. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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

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

Synthesis,TEM characterization and thermal behaviour of LiNiSi<Subscript>2</Subscript>O<Subscript>6</Subscript> pyroxene

A pyroxene with composition LiNiSi₂O₆ was synthesized at T = 1,473 K and P = 2.0 GPa; the cell parameters at T = 298 K are a = 9.4169(6) Å, b = 8.4465(7) Å, c = 5.2464(3) Å, β = 110.534(6)°, V = 390.78(3) Å³. TEM examination of the LiNiSi₂O₆ pyroxene showed the presence of h + k odd reflections indicative of a primitive lattice, and of antiphase domains obtained by dark field imaging of the h + k odd reflections. A HT in situ investigation was performed by examining TEM selected area diffraction patterns collected at high temperature and synchrotron radiation powder diffraction. In HTTEM the LiNiSi₂O₆ was examined together with LiCrSi₂O₆ pyroxene. In LiCrSi₂O₆ the h + k odd critical reflections disappear at about 340 K; they are sharp up to the transition temperature and do not change their shape until they disappear. In LiNiSi₂O₆ the h + k odd reflections are present up to sample deterioration at 650 K. A high temperature synchrotron radiation powder diffraction investigation was performed on LiNiSi₂O₆ between 298 and 773 K. The analysis of critical reflections and of changes in cell parameters shows that the space group is P-centred up to the highest temperature. The comparative analysis of the thermal and spontaneous strain contributions in P2₁/c and C2/c pyroxenes indicates that the high temperature strain in P-LiNiSi₂O₆ is very similar to that due to thermal strain only in C2/c spodumene and that a spontaneous strain contribution related to pre-transition features is not apparent in LiNiSi₂O₆. A different high-temperature behaviour in LiNiSi₂O₆ with respect to other pyroxenes is suggested, possibly in relation with the presence of Jahn–Teller distortion of the M1 polyhedron centred by low-spin Ni³⁺.     

High-pressure synthesis,long-term stability of single crystals of diboron trioxide,B<Subscript>2</Subscript>O<Subscript>3</Subscript>, and an empirical electronic polarizability of <Superscript>[3]</Superscript>B<Superscript>3+</Superscript>

Single crystals of B₂O₃ are needed for the precise determination of the refractive indices used to calculate the electronic polarizability α of 3-coordinated boron. The α(B) values in turn are used to predict mean refractive indices of borate minerals. Since the contribution of boron to the total polarizability of a mineral is very low, the synthetic compound B₂O₃ represents an ideal model system because of its high molar content of boron. Millimeter-sized crystals were synthesized at 1 GPa in a piston-cylinder apparatus. The samples were heated above the liquidus (800 °C), subsequently cooled at 15 °C/h to 500 °C and finally quenched. The refractive indices were determined by the immersion method using a microrefractometer spindle stage. The refractive indices n _o = 1.653 (3) and n _e = 1.632 (3) correspond to a total polarizability for B₂O₃ of α = 4.877 Å³. These values were used to determine the electronic polarizability of boron of α(B) = 0.16 Å³. Although the surface of the B₂O₃ crystals was coated with a hydrous film immediately after being exposed to air, its bulk crystallinity is retained for a period of at least 2 months.     

Simulation of the molecular mass distributions of Si anions on the liquidus of the Na<Subscript>2</Subscript>O-SiO<Subscript>2</Subscript> system with regard for the disproportionation of <Emphasis Type="Italic">Q</Emphasis> structons

A new version of the STRUCTON (2009) computer model is proposed for the simulation of the molecular mass distributions (MMD) characterizing the diversity of anions in silicate melts depending on their polymerization and temperature. In contrast to earlier versions, the new version of the model accounts for disproportionation reactions of Q ⁿ species and makes use of their proportions in the statistical simulations of the origin of real Si-O complexes. The new potentialities of the STRUCTON program package are illustrated by its application to studying the structural-chemical characteristics of melts in the Na₂O-SiO₂ system along its liquidus line, including the points of eutectics and phase transitions at 0.333 ≤ $ N_{SiO_2 } $ N_{SiO_2 } < 0.500. This problem is solved with the use of a temperature-composition dependence of polymerization constants K _p ^Na in the Toop-Samis approximation. The variations in K _p ^Na were proved to be as large as three orders of magnitude due to both the temperature effect at a constant composition and the composition effect at a constant temperature. The results of the MMD simulations on the liquidus show that the concentration of the SiO₄⁴⁻ ion strongly decreases, and the proportion of chain species increases compared to those at a stochastic distribution. The concentration of the Si₂O₇⁶⁻ anion reaches its maximum (∼42%) at 40 mol % in the liquid, i.e., the composition of Na₆Si₂O₇. At $ N_{SiO_2 } $ N_{SiO_2 } > 0.40, this ion dominates over the SiO₄⁴⁻ monomer. More silicic melts with $ N_{SiO_2 } $ N_{SiO_2 } ≥ 0.45, are dominated by (Si _n O_3n )³ⁿ⁻ ring species, and the concentrations of these species are related as (Si₃O₉)⁶⁻ > (Si₄O₁₂)⁸⁻ > (Si₅O₁₅)¹⁰⁻. The maximum concentration of these flat rings also occurs near the composition of stoichiometric metasilicate with Si/O = 0.333. The comparison of the dependence of the average size of anions i _av and the average number of their species on depolymerization indicates that a change in the proportion of Q ⁿ species in melt at decreasing temperature results in structural restyling and an increase in the average size of Si-O complexes. The average number of anion species thereby decreases compared to that in a stochastic MMD. The results presented in this publication direct the progress in the thermodynamic theory of silicate melts to a new avenue that makes use of the capabilities and advantages of the ion-polymer model, the theory of associated solutions, spectroscopic data, and the experimental study of variations in oxide activities depending on composition and temperature.     

Equation of state,structural behaviour and phase diagram of synthetic MgFe<Subscript>2</Subscript>O<Subscript>4</Subscript>, as a function of pressure and temperature

The behaviour of synthetic Mg-ferrite (MgFe₂O₄) has been investigated at high pressure (in situ high-pressure synchrotron radiation powder diffraction at ESRF) and at high temperature (in situ high-temperature X-ray powder diffraction) conditions. The elastic properties determined by the third-order Birch–Murnaghan equation of state result in K₀=181.5(± 1.3) GPa, K=6.32(± 0.14) and K= –0.0638 GPa^–1. The symmetry-independent coordinate of oxygen does not show significant sensitivity to pressure, and the structure shrinking is mainly attributable to the shortening of the cell edge (homogeneous strain). The lattice parameter thermal expansion is described by _a0+_a1*(T–298)+_a2/(T–298)², where _a0=9.1(1) 10^–6 K^–1, _a1=4.9(2) 10^–9 K^–2 and _a2= 5.1(5) 10^–2 K. The high-temperature cation-ordering reaction which MgFe-spinel undergoes has been interpreted by the ONeill model, whose parameters are = 22.2(± 1.8) kJ mol^–1 and =–17.6(± 1.2) kJ mol^–1. The elastic and thermal properties measured have then been used to model the phase diagram of MgFe₂O₄, which shows that the high-pressure transition from spinel to orthorombic CaMn₂O₄-like structure at T < 1700 K is preceded by a decomposition into MgO and Fe₂O₃.