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

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

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.     

Adsorption of H<Subscript>2</Subscript>O,NH<Subscript>3</Subscript> and C<Subscript>6</Subscript>H<Subscript>6</Subscript> on alkali metal cations in internal surface of mordenite and in external surface of smectite: a DFT study

Adsorption of H₂O, NH₃ and C₆H₆ on H- and alkali metal-exchanged structures of mordenite and on corresponding cations on the smectite layer is investigated by ab initio density-functional calculations. Proton or an alkali metal cation compensates one Al/Si framework substitution and resides in the extra-framework position of zeolite or above flat smectite layer close to the Al/Si substitution. Pronounced similarities between zeolite and smectite are observed in changes of the adsorption energies and location of the external cation with changing character of the external cation. Calculated adsorption energies exhibit the following trend: E(NH₃) > E(H₂O) > E(C₆H₆). Because of looser contact with the framework, zeolitic cations are stronger adsorption centers and calculated adsorption energies of zeolites are by ~20–30% larger than cations of smectites. The highest adsorption energy is calculated for H-exchanged structures and down the group of alkali metal cations a decrease of the adsorption energy is observed. Deviations from the smooth variation of the adsorption energy are caused by: (1) formation of strong hydrogen bonds in H-exchanged structures, (2) adsorption induced migration of the external Li⁺ cation, and (3) steric hindrances of the flat C₆H₆ molecule adsorbed on the cation in the cage of zeolite.     

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.     

Pressure and temperature dependence of the viscosity of a NaAlSi<Subscript>2</Subscript>O<Subscript>6</Subscript> melt

The viscosity of a silicate melt of composition NaAlSi₂O₆ was measured at pressures from 1.6 to 5.5 GPa and at temperatures from 1,350 to 1,880°C. We employed in situ falling sphere viscometry using X-ray radiography. We found that the viscosity of the NaAlSi₂O₆ melt decreased with increasing pressure up to 2 GPa. The pressure dependence of viscosity is diminished above 2 GPa. By using the relationship between the logarithm of viscosity and the reciprocal temperature, the activation energies for viscous flow were calculated to be 3.7 ± 0.4 × 10² and 3.7 ± 0.5 × 10² kJ/mol at 2.2 and 2.9 GPa, respectively.     

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.     

High-pressure thermo-elastic properties of beryl (Al<Subscript>4</Subscript>Be<Subscript>6</Subscript>Si<Subscript>12</Subscript>O<Subscript>36</Subscript>) from ab initio calculations,and observations about the source of thermal expansion

Ab initio calculations of thermo-elastic properties of beryl (Al₄Be₆Si₁₂O₃₆) have been carried out at the hybrid HF/DFT level by using the B3LYP and WC1LYP Hamiltonians. Static geometries and vibrational frequencies were calculated at different values of the unit cell volume to get static pressure and mode-γ Grüneisen’s parameters. Zero point and thermal pressures were calculated by following a standard statistical-thermodynamics approach, within the limit of the quasi-harmonic approximation, and added to the static pressure at each volume, to get the total pressure (P) as a function of both temperature (T) and cell volume (V). The resulting P(V, T) curves were fitted by appropriate EoS’, to get bulk modulus (K ₀) and its derivative (K′), at different temperatures. The calculation successfully reproduced the available experimental data concerning compressibility at room temperature (the WC1LYP Hamiltonian provided K ₀ and K′ values of 180.2 Gpa and 4.0, respectively) and the low values observed for the thermal expansion coefficient. A zone-centre soft mode \( P6/mcc \to P\bar{1} \) phase transition was predicted to occur at a pressure of about 14 GPa; the reduction of the frequency of the soft vibrational mode, as the pressure is increased, and the similar behaviour of the majority of the low-frequency modes, provided an explanation of the thermal behaviour of the crystal, which is consistent with the RUM model (Rigid Unit Model; Dove et al. in Miner Mag 59:629–639, 1995), where the negative contribution to thermal expansion is ascribed to a geometric effect connected to the tilting of rigid polyhedra in framework silicates.     

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.     

Possible mantle phase transitions by the formation of SiO<Subscript>2</Subscript> peroxides: Implications for mantle convection

On the basis of quantum–chemical calculations of the linear to isomeric bent transition of the SiO₂ molecule, it is suggested that the bent to linear transition of SiO₂ forms can occur in melted mantle minerals of the lower mantle. This may be important for the formation of the peculiarities of mantle convection and origination of plumes.     

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.     

TiO<Subscript>2</Subscript>-based shaped catalyst for the recovery of elemental sulfur from H<Subscript>2</Subscript>S and SO<Subscript>2</Subscript> gas streams

The Claus process has been used for the conversion of H₂S and SO₂ to elemental sulfur. These two sulfur compounds need special attention because they are very poisonous with negative impact on both the environment and human health. Here, highly active Fe–Ni/TiO₂ catalyst has been prepared and shaped by three different binders (bentonite, polyethylene glycol and carboxymethyl cellulose) into extrudes. Comparing the mechanical strength and surface area of prepared extrudes, the optimal shaped catalyst was selected with 20% of bentonite, 2% of PEG and 2% of CMC. The optimal catalyst was characterized by X-ray powder diffraction, temperature-programmed reduction, Brunauer–Emmett–Teller specific surface area, Barrett–Joyner–Halenda, scanning electron microscopy and energy-dispersive X-ray techniques and used for sulfur recovery process. The performance of this product for sulfur recovery via Claus process was excellent with the conversion of hydrogen sulfide of 76.77% and sulfur dioxide of 97.83%. The catalyst also provides high hydrolysis activity of CS₂ (83.06%). Therefore, a highly active TiO₂-supported shaped catalyst with 85.62% of conversion efficiency has been prepared successfully to convert the small amounts of H₂S, SO₂ and CS₂ to elemental sulfur.     

Cocrystallization of isomorphic components in solutions and crystal zoning: An example of the (Ba,Pb)(NO<Subscript>3</Subscript>)<Subscript>2</Subscript> series

The phase equilibrium and growth of mixed (Ba,Pb)(NO₃)₂ crystals in aqueous solutions were investigated. The microcrystallization method was adapted to systems with isomorphic components, where crystal composition deviates from the thermodynamic equilibrium even at a low supercooling. The solid phase is characterized by continuous miscibility of the components within the analyzed ranges of temperature and solution composition. The diagram is characterized by substantial nonlinearity of the solubility isotherms and nonuniformity of the solid isocomposite location. Quasiequilibrium paths change their slope depending on solution enrichment in Ba(NO₃)₂ and depletion in Pb(NO₃)₂. The area of curvature of quasiequilibrium paths is coordinated with the area of the changing isotherm slope, i.e., the area of “remarkable” points, where the behavior of the system changes drastically. Examples of theoretical zoning of a crystal approximated to a sphere were calculated for paths at a temperature decreasing from 50 to 15°C. The content of the Ba component decreases toward the periphery with the consecutive overgrowing of zones. Each zone corresponds to one gram of (Ba,Pb)(NO₃)₂. Crystals grown from different initial solutions consist of different numbers of zones owing to the nonlinearity of solubility isotherms. Specific features of mixed crystal formation should be taken into account in genetic interpretation of natural minerals of mixed composition.     

The Role of CaCO<Subscript>3</Subscript> Reactions in the Contemporary Oceanic CO<Subscript>2</Subscript> Cycle

The present analysis adjusts previous estimates of global ocean CaCO₃ production rates substantially upward, to 133 × 10¹² mol yr^?1 plankton production and 42 × 10¹² mol yr^?1 shelf benthos production. The plankton adjustment is consistent with recent satellite-based estimates; the benthos adjustment includes primarily an upward adjustment of CaCO₃ production on so-called carbonate-poor sedimentary shelves and secondarily pays greater attention to high CaCO₃ mass (calcimass) and turnover of shelf communities on temperate and polar shelves. Estimated CaCO₃ sediment accumulation rates remain about the same as they have been for some years: ~20 × 10¹² mol yr^?1 on shelves and 11 × 10¹² mol yr^?1 in the deep ocean. The differences between production and accumulation of calcareous materials call for dissolution of ~22 × 10¹² mol yr^?1 (~50 %) of shelf benthonic carbonate production and 122 × 10¹² mol yr^?1 (>90 %) of planktonic production. Most CaCO₃ production, whether planktonic or benthonic, is assumed to take place in water depths of <100 m, while most dissolution is assumed to occur below this depth. The molar ratio of CO₂ release to CaCO₃ precipitation (CO₂↑/CaCO₃↓) is <1.0 and varies with depth. This ratio, Ψ, is presently about 0.66 in surface seawater and 0.85 in ocean waters deeper than about 1000 m. The net flux of CO₂ associated with CaCO₃ reactions in the global ocean in late preindustrial time is estimated to be an apparent influx from the atmosphere to the ocean, of +7 × 10¹² mol C yr^?1, at a time scale of 10²–10³ years. The CaCO₃-mediated influx of CO₂ is approximately offset by CO₂ release from organic C oxidation in the water column. Continuing ocean acidification will have effects on CaCO₃ and organic C metabolic responses to the oceanic inorganic C cycle, although those responses remain poorly quantified.     

Contents of Fixed-Ammonium （ NH4^＋ ） in Lamprophyres in the Zhenyuan Gold Orefield, Yunnan Province, China： Implications for Its Characteristics of the Source Region

This paper determined the fixed-ammonium (NH₄ ⁺) contents of lamprophyres in the Zhenyuan gold orefield, Yunnan Province, China. The results show that the NH₄ ⁺ contents of minettes in the orefield range from 120 × 10⁻⁶ to 469 × 10⁻⁶ and those of kersantites from 74.3 × 10⁻⁶ to 136 × 10⁻⁶. These values are higher than those of other mantle-derived rocks (less than 50 × 10⁻⁶), but lower than those of carbonaceous wall rocks in the orefield (from 1200 × 10⁻⁶ to 1343 × 10⁻⁶). Combining with the Sr isotopic composition, this paper suggested that lamprophyres in the orefield with high NH₄ ⁺ contents relative to other mantle-derived rocks would not have resulted from the primary magma contaminated by crustal materials in the process of rising or in the magma chamber, but from mantel metasomatism. This research project was financially supported jointly by the State Key Basic Research Program of the People’s Republic of China (No. G1999043203) and the Chinese Academy of Sciences Knowledge-Innovation Program (No. KZCX3-SW-125).     

Barioferrite BaFe<Subscript>12</Subscript>O<Subscript>19</Subscript>: A new mineral species of the magnetoplumbite group from the Haturim Formation in Israel

A new mineral barioferrite—a natural analogue of synthetic barium ferrite Ba Fe ₁₂ ³⁺ O₁₉—has been identified in the central part of a metamorphosed barite nodule in the rock of the Haturim Formation (Mottled Zone) on the southern slope of Mount Ye’elim in Israel. The mineral is associated with barite, calcite, magnetite, and maghemite and occurs as tiny platy crystals up to 3 × 15 × 15 μm and their irregular aggregates. Barioferrite is black with streaks of brown, and its luster is submetallic. Its Calculated density is 5.31 g/cm³. The mineral is brittle; cleavage is absent. IR absorption bands (cm^?1) are observed at 635 (shoulder), 582, 544, 433, and 405 (shoulder). Barioferrite is characterized by ferrimagnetic behavior. Under a microscope in reflected light, barioferrite is grayish white with brownish red internal reflections, the pleochroism is weak (from gray-white on R _o to gray-white with a brown tint on R _e), and the bireflectance is weak with distinct anisotropy. The reflectance values of R _o/R _e, % (λ, nm) are 24.51/22.80 (470), 24.17/22.25 (546), 23.65/21.68 (589), and 22.67/20.85 (650). The chemical composition (electron microprobe, wt %; the ranges are given in parentheses) is BaO 13.13 (12.5–13.8), Fe₂O₃ 86.47 (85.5–87.5), and 99.60 in total. The empirical formula is Ba_0.95Fe _12.03 ³⁺ O₁₉. Barioferrite is hexagonal with space group P6₃/mmc, a = 5.875 (3) Å, c = 23.137 (19) Å, V = 691.6 (5) ų, and Z = 2. The strongest lines of the X-ray powder diffraction pattern [d, Å, (I, 5) (hkl)] are 2.938(46) (110), 2.770(100) (107), 2.624 (84) (114, 200), 2.420(44) (203), 2.225(40) (205), and 1.627(56) (304, 2.0.11). The holotype specimen of barioferrite is deposited at the Mineralogical Museum of St. Petersburg State University; its catalogue number is 1/19436.     

Crystal chemistry of selenates with mineral-like structures: VIII. Butlerite chains in the structure of K(UO<Subscript>2</Subscript>)(SeO<Subscript>4</Subscript>)(OH)(H<Subscript>2</Subscript>O)

A new potassium uranyl selenate compound K(UO₂)(SeO₄)(OH)(H₂O) has been synthesized for the first time using the technique of evaporation from water solution. Its crystal structure has been solved by direct methods (monoclinic, P2₁/c,a = 8.0413(9) Å, b = 8.0362(9) Å, c = 11.6032(14) Å, β = 106.925(2)°, V = 717.34(14) ų) and refined to R ₁ = 0.0319 (wR ₂ = 0.0824) for 1285 reflections with |F ₀| > 4σ _F . The structure consists of [(UO₂(SeO₄)(OH)(H₂O)]^? chains extending along axis b. In the chains, the uranyl pentagonal bipyramids are linked via bridged hydroxyl anions and tetrahedral oxoanions [SeO₄]^2?. Potassium ions are situated between these chains. No chains of that type have been observed in uranyl compounds earlier, but they had been detected in the structures of butlerite, parabutlerite, uklonskovite, fibroferrite, and a number of synthetic compounds.