Electric conductivity of Fe2SiO4–Fe3O4 spinel solid solutions

 The spinel solid solution was found to exist in the whole range between Fe₃O₄ and γ-Fe₂SiO₄ at over 10 GPa. The resistivity of Fe_{3− x} Si _x O₄ (0.0<x<0.288) was measured in the temperature range of 80∼300 K by the AC impedance method. Electron hopping between Fe³⁺ and Fe²⁺ in the octahedral site of iron-rich phases gives a large electric conductivity at room temperature. The activation energy of the electron hopping becomes larger with increasing γ-Fe₂SiO₄ component. A nonlinear change in electric conductivity is not simply caused by the statistical probability of Fe³⁺–Fe²⁺ electron hopping with increasing the total Si content. This is probably because a large number of Si⁴⁺ ions occupies the octahedral site and the adjacent Fe²⁺ keeping the local electric neutrality around Si⁴⁺ makes a cluster, which generates a local deformation by Si substitution. The temperature dependence of the conductivity of solid solutions indicates the Verwey transition temperature, which decreases from 124(±2) K at x=0 (Fe₃O₄) to 102(±5) K at x=0.288, and the electric conductivity gap at the transition temperature decreases with Si⁴⁺ substitution. Received: 15 March 2000 / Accepted: 4 September 2000     

Solubility and solution mechanism of H2O in alkali silicate melts and glasses at high pressure and temperature

The solubility behavior of H₂O in melts in the system Na₂O-SiO₂-H₂O was determined by locating the univariant phase boundary, melt = melt + vapor in the 0.8-2 GPa and 1000°-1300°C pressure and temperature range, respectively. The NBO/Si-range of the melts (0.25-1) was chosen to cover that of most natural magmatic liquids. The H₂O solubility in melts in the system Na₂O-SiO₂-H₂O (X_H2O) ranges between 18 and 45 mol% (O = 1) with (∂X_H2O/∂P)_T∼14-18 mol% H₂O/GPa. The (∂X_H2O/∂P)_T is negatively correlated with NBO/Si (= Na/Si) of the melt. The (∂X_H2O/∂T)_P is in the −0.03 to +0.05 mol% H₂O/°C range, and is negatively correlated with NBO/Si. The [∂X_H2O/∂(NBO/Si)]_P,T is in the −3 to −8 mol% H₂O/(NBO/Si) range. Melts with NBO/Si similar to basaltic liquids (∼0.6-∼1.0) show (∂X_H2O/∂T)_P<0, whereas more polymerized melts exhibit (∂X_H2O/∂T)_P>0. Complete miscibility between hydrous melt and aqueous fluid occurs in the 0.8-2 GPa pressure range for melts with NBO/Si ≤0.5 at T >1100°C. Miscibility occurs at lower pressure the more polymerized the melt.     

Calorimetric study of the stability of spinelloids in the system NiAl2O4-Ni2SiO4

Enthalpies of solution in molten 2 PbO · B₂O₃ at 974 K were measured for four spinelloids, phases I (0.75 NiAl₂O₄ · 0.25 Ni₂SiO₄), II (0.60 NiAl₂O₄ · 0.40 Ni₂SiO₄), III and IV (0.50 NiAl₂O₄ · 0.50 Ni₂SiO₄) in the system NiAl₂O₄ · Ni₂SiO₄. The enthalpies (in cal per 4-oxygen mol) of formation from NiAl₂O₄ and Ni₂SiO₄ spinels are: phase I, 945±366; phase II, 1072±360; phase III, 2253±390; phase IV, 3565±544. Using these enthalpy data in combination with phase relations at high pressure at 1373 K, positive entropies of formation of the spinelloids from NiAl₂O₄ and Ni₂SiO₄ spinels were estimated (in cal mol^?1 K^?1): phase I, 1.2; phase II, 1.5; phase III, 2.0–2.3; phase IV, 3.0–3.1. The thermochemical data obtained above suggest that the spinelloids are “entropy-stabilized” phases with partially disordered cation distributions. The configurational entropies of the spinelloids were calculated based on the observed cation distribution in each spinelloid phase. The positive entropies of formation of the spinelloids from the spinel endmembers are due primarily to the configurational entropies although small positive vibrational entropy changes may also exist.     

Hydrogen isotope exchange kinetics between H2O and H4SiO4 from ab initio calculations

Hydrogen isotope exchange between water and orthosilicic acid (H₄SiO₄) was modeled using B3LYP calculations and classical transition-state theory. Configurations of 1, 2, 3 and 7 water molecules and H₄SiO₄ were used to investigate energetically viable reaction pathways. An upper-bound of 71 kJ/mol was assumed for the zero-point energy corrected barrier (ZPECB) because this is the experimentally determined activation energy for Si-O bond breaking (Rimstidt and Barnes, 1980) and ZPECB is expected to be close to this value. Long range solvation forces were accounted for using the integral equation formalism polarized continuum model (IEFPCM; Cancès et al., 1997). Primary and secondary isotope effects were computed by exchanging hydrogen atoms with deuterium. Results show that reaction mechanisms involving 3 and 7 water molecules have ZPECB of 34 to 38 kJ/mol, whereas those involving 1 and 2 water molecules have ZPECB in excess of the set upper-bound. The lower range of ZPECB with 3 or 7 water molecules is reasonable to explain rapid hydrogen isotope exchange with silicates. Rate constant calculations accounting for tunneling, anharmonicity and scaling factors indicate that the reaction is fast and equilibrium can be assumed under most geologic conditions.     

Molecular H2O in armenite, BaCa2Al6Si9O30·2H2O, and epididymite, Na2Be2Si6O15·H2O: Heat capacity, entropy and local-bonding behavior of confined H2O in microporous silicates

Armenite, ideal formula BaCa₂Al₆Si₉O₃₀·2H₂O, and its dehydrated analog BaCa₂Al₆Si₉O₃₀ and epididymite, ideal formula Na₂Be₂Si₆O₁₅·H₂O, and its dehydrated analog Na₂Be₂Si₆O₁₅ were studied by low-temperature relaxation calorimetry between 5 and 300 K to determine the heat capacity, C_p, behavior of their confined H₂O. Differential thermal analysis and thermogravimetry measurements, FTIR spectroscopy, electron microprobe analysis and powder Rietveld refinements were undertaken to characterize the phases and the local environment around the H₂O molecule.The determined structural formula for armenite is Ba_0.88(0.01)Ca_1.99(0.02)Na_0.04(0.01)Al_5.89(0.03)Si_9.12(0.02)O₃₀·2H₂O and for epididymite Na_1.88(0.03)K_0.05(0.004)Na_0.01(0.004)Be_2.02(0.008)Si_6.00(0.01)O₁₅·H₂O. The infrared (IR) spectra give information on the nature of the H₂O molecules in the natural phases via their H₂O stretching and bending vibrations, which in the case of epididymite only could be assigned. The powder X-ray diffraction data show that armenite and its dehydrated analog have similar structures, whereas in the case of epididymite there are structural differences between the natural and dehydrated phases. This is also reflected in the lattice IR mode behavior, as observed for the natural phases and the H₂O-free phases. The standard entropy at 298 K for armenite is S° = 795.7 ± 6.2 J/mol K and its dehydrated analog is S° = 737.0 ± 6.2 J/mol K. For epididymite S° = 425.7 ± 4.1 J/mol K was obtained and its dehydrated analog has S° = 372.5 ± 5.0 J/mol K. The heat capacity and entropy of dehydration at 298 K are Δ = 3.4 J/mol K and ΔS^rxn = 319.1 J/mol K and Δ = −14.3 J/mol K and ΔS^rxn = 135.7 J/mol K for armenite and epididymite, respectively. The H₂O molecules in both phases appear to be ordered. They are held in place via an ion-dipole interaction between the H₂O molecule and a Ca cation in the case of armenite and a Na cation in epididymite and through hydrogen-bonding between the H₂O molecule and oxygen atoms of the respective silicate frameworks. Of the three different H₂O phases ice, liquid water and steam, the C_p behavior of confined H₂O in both armenite and epididymite is most similar to that of ice, but there are differences between the two silicates and from the C_p behavior of ice. Hydrogen-bonding behavior and its relation to the entropy of confined H₂O at 298 K is analyzed for various microporous silicates.The entropy of confined H₂O at 298 K in various silicates increases approximately linearly with increasing average wavenumber of the OH-stretching vibrations. The interpretation is that decreased hydrogen-bonding strength between a H₂O molecule and the silicate framework, as well as weak ion-dipole interactions, results in increased entropy of H₂O. This results in increased amplitudes of external H₂O vibrations, especially translations of the molecule, and they contribute strongly to the entropy of confined H₂O at T < 298 K.     

Thermodynamic properties of H4SiO4 in the ideal gas state as evaluated from experimental data

Solid phases of silicon dioxide react with water vapor with the formation of hydroxides and oxyhydroxides of silica. Recent transpiration and mass-spectrometric studies convincingly demonstrate that H₄SiO₄ is the predominant form of silica in vapor phase at water pressure in excess of 10⁻² MPa. Available literature transpiration and solubility data for the reactions of solid SiO₂ phases and low-density water, extending from 424 to 1661 K, are employed for the determination of Δ_fG⁰, Δ_fH⁰ and S⁰ of H₄SiO₄ in the ideal gas state at 298.15 K, 0.1 MPa. In total, there are 102 data points from seven literature sources. The resulting values of the thermodynamic functions of H₄SiO₄(g) are: Δ_fG⁰ = −1238.51 ± 3.0 kJ mol⁻¹, Δ_fH⁰ = −1340.68 ± 3.5 kJ mol⁻¹ and S⁰ = 347.78 ± 6.2 J K⁻¹ mol⁻¹. These values agree quantitatively with one set of ab initio calculations. The relatively large uncertainties are mainly due to conflicting data for H₄SiO₄(g) from various sources, and new determinations of would be helpful. The thermodynamic properties of this species, H₄SiO₄(g), are necessary for realistic modeling of silica transport in a low-density water phase. Applications of this analysis may include the processes of silicates condensation in the primordial solar nebula, the precipitation of silica in steam-rich geothermal systems and the corrosion of SiO₂-containing alloys and ceramics in moist environments.     

Residence times for protons bound to three oxygen sites in the AlO4Al12(OH)24(H2O)12 polyoxocation

Although, the kinetic reactivity of a mineral surface is determined, in part, by the rates of exchange of surface-bound oxygens and protons with bulk solution, there are no elementary rate data for minerals. However, such kinetic measurements can be made on dissolved polynuclear clusters, and here we report lifetimes for protons bound to three oxygen sites on the AlO₄Al₁₂(OH)₂₄(H₂O)₁₂⁷⁺ (Al₁₃) molecule, which is a model for aluminum-hydroxide solids in water. Proton lifetimes were measured using ¹H NMR at pH ∼ 5 in both aqueous and mixed solvents. The ¹H NMR peak for protons on bound waters (η-H₂O) lies near 8 ppm in a 2.5:1 mixture of H₂O/acetone-d₆ and broadens over the temperature range −20 to −5 °C. Extrapolated to 298 K, the lifetime of a proton on a η-H₂O is τ²⁹⁸ ∼ 0.0002 s, which is surprisingly close to the lifetime of an oxygen in the η-H₂O (∼0.0009 s), but in the same general range as lifetimes for protons on fully protonated monomer ions of trivalent metals (e.g., Al(H₂O)₆³⁺). The lifetime is reduced somewhat by acid addition, indicating that there is a contribution from the partly deprotonated Al₁₃ molecule in addition to the fully protonated Al₁₃ at self-buffered pH conditions. Proton lifetimes on the two distinct sets of hydroxyls bridging two Al(III) (μ₂-OH) differ substantially and are much shorter than the lifetime of an oxygen at these sites. The average lifetimes for hydroxyl protons were measured in a 2:1 mixture of H₂O/dmso-d₆ over the temperature range 3.7-95.2 °C. The lifetime of a hydrogen on one of the μ₂-OH was also measured in D₂O. The τ²⁹⁸ values are ∼0.013 and ∼0.2 s in the H₂O/dmso-d₆ solution and the τ²⁹⁸ value for the μ₂-OH detectable in D₂O is τ²⁹⁸ ∼ 0.013 s. The ¹H NMR peak for the more reactive μ₂-OH broadens slightly with acid addition, indicating a contribution from an exchange pathway that involves a proton or hydronium ion. These data indicate that surface protons on minerals will equilibrate with near-surface waters on the diffusional time scale.     

Effect of oxygen fugacity on the H2O storage capacity of forsterite in the carbon-saturated systems

High pressure experiments have been performed in the systems Mg₂SiO₄-C-O-H and Mg₂SiO₄-K₂CO₃-C at 6.3 GPa and 1200 to 1600 °C using a split-sphere multi-anvil apparatus. In the Mg₂SiO₄-C-O-H system the composition of fluid was modeled by adding different amounts of water and stearic acid. The fO₂ was controlled by the Mo-MoO₂ or Fe-FeO oxygen buffers. Several experiments in the Mg₂SiO₄-C-O-H system and all experiments in the Mg₂SiO₄-K₂CO₃-C system have been conducted without buffering the fO₂. Forsterite in the system Mg₂SiO₄-K₂CO₃-C does not reveal OH absorption bands in the IR spectra, while forsterite coexisting with carbon-bearing fluid and silicate melt at logfO₂ from FMQ-2 to FMQ-5 (from 2 to 5 log units below fayalite-magnetite-quartz oxygen buffer) contains 800-1850 wt. ppm H₂O. The maximum concentrations were detected at 1400 °C and FMQ-3.5. We observed an increase in the solidus temperature in the system Mg₂SiO₄-C-O-H from 1200 to above 1600 °C with log fO₂ decreasing from FMQ-2 to FMQ-5. The increase of the solidus temperature and the broadening of the stability field of the H₂O-H₂-CH₄ subsolidus fluid phase at 1400-1600 °C explain the high H₂O storage capacity of forsterite relative to that crystallized from carbon-free, oxidized, hydrous, silicic melt. At temperatures above 1400 °C liquidus forsterite precipitated along with diamond from oxidized (FMQ-1) carbonate-silicate melt and from silicate melt dissolving the moderately reduced C-O-H fluid (from FMQ-2 to FMQ-3.5). Formation of diamond was not detected under ultra-reduced conditions (FMQ-5) at 1200-1600 °C. Olivine co-precipitating with diamond from dry carbonate-silicate or hydrous-silicic fluid/melt can provide information on the H₂O contents and speciation of the diamond-forming media in the mantle. The conditions for minimum post-crystallization alteration of olivine and its hydrogen content are discussed.     

花岗岩—H2O—HF体系相关系及氟对花岗质熔体结构的影响

通过在0.1GPa压力下钠长花岗岩-H₂O-HF体系相关系实验获得,随体系F含量的增加,石英的温度稳定域上限升高,长石的温度稳定域上限降低;石英、碱性长石的晶-液平衡热力学计算表明,F导致花岗质熔体中组分SiO₂的活度增加,组分NaAlSi₃O₈和KAlSi₃O₈的活度减小,且NaAlSi₃O₈活度较KAlSi₃O₈活度减小幅度大。这些结果显示了F在花岗质熔体中与Si以外的阳离子Al、Na、K等产生了结合,且F与Na结合的优先性大于K,破坏了具有电荷平衡离子Na、K的AlO₂^-四面体,使熔体架状网格中Si/(Si+Al)和K/Na比值增大。通过F与H₂O对花岗岩体系相平衡的影响比较,作者认为F不与Si结合是它与OH^-在干扰熔体结构方面的最大区别。     

H2O diffusion in dacitic and andesitic melts

The diffusion of water in dacitic and andesitic melts was investigated at temperatures of 1458 to 1858 K and pressures between 0.5 and 1.5 GPa using the diffusion couple technique. Pairs of nominally dry glasses and hydrous glasses containing between 1.5 and 6.3 wt.% dissolved H₂O were heated for 60 to 480 s in a piston cylinder apparatus. Concentration profiles of hydrous species (OH groups and H₂O molecules) and total water (C_H2Ot = sum of OH and H₂O) were measured along the cylindrical axis of the diffusion sample using IR microspectroscopy. Electron microprobe traverses show no significant change in relative proportions of anhydrous components along H₂O profiles, indicating that our data can be treated as effective binary interdiffusion between H₂O and the rest of the silicate melt. Bulk water diffusivity (D_H2Ot) was derived from profiles of total water using a modified Boltzmann-Matano method as well as using fittings assuming a functional relationship between D_H2Ot and C_H2Ot. In dacitic melts D_H2Ot is proportional to C_H2Ot up to 6 wt.%. In andesitic melts the dependence of D_H2Ot on C_H2Ot is less pronounced. A pressure effect on water diffusivity could not be resolved for either dacitic or andesitic melt in the range 0.5 to 1.5 GPa. Combining our results with previous studies on water diffusion in rhyolite and basalt show that for a given water content D_H2Ot increases monotonically with increasing melt depolymerization at temperatures >1500 K. Assuming an Arrhenian behavior in the whole compositional range, the following formulation was derived to estimate D_H2Ot (m²/s) at 1 wt.% H₂O_t in melts with rhyolitic to andesitic composition as a function of T (K), P (MPa) and S (wt.% SiO₂):

     

Thermodynamics of open networks: Ordering and entropy in NaAlSiO4 glass,liquid, and polymorphs

The thermodynamic properties of carnegieite and NaAlSiO₄ glass and liquid have been investigated through C _p determinations from 10 to 1800 K and solution-calorimetry measurements. The relative entropies S ₂₉₈-S₀ of carnegieite and NaAlSiO₄ glass are 118.7 and 124.8 J/mol K, respectively. The low-high carnegieite transition has been observed at 966 K with an enthalpy of transition of 8.1±0.3 kJ/mol, and the enthalpy of fusion of carnegieite at the congruent melting point of 1799 K is 21.7±3 kJ/mol. These results are consistent with the reported temperature of the nepheline-carnegieite transition and available thermodynamic data for nepheline. The entropy of quenched NaAlSiO₄ glass at 0 K is 9.7±2 J/mol K and indicates considerable ordering among AlO₄ and SiO₄ tetrahedra. In the liquid state, progressive, temperature-induced Si, Al disordering could account for the high configurational heat capacity. Finally, the differences between the entropies and heat capacities of nepheline and carnegieite do not seem to conform to current polyhedral modeling of these properties     

Low-temperature heat capacities,entropies and enthalpies of Mg<Subscript>2</Subscript>SiO<Subscript>4</Subscript> polymorphs,and α−β−γ and post-spinel phase relations at high pressure

The low-temperature isobaric heat capacities (C _p) of β- and γ-Mg₂SiO₄ were measured at the range of 1.8–304.7 K with a thermal relaxation method using the Physical Property Measurement System. The obtained standard entropies (S°₂₉₈) of β- and γ-Mg₂SiO₄ are 86.4 ± 0.4 and 82.7 ± 0.5 J/mol K, respectively. Enthalpies of transitions among α-, β- and γ-Mg₂SiO₄ were measured by high-temperature drop-solution calorimetry with gas-bubbling technique. The enthalpies of the α−β and β−γ transitions at 298 K (ΔH°₂₉₈) in Mg₂SiO₄ are 27.2 ± 3.6 and 12.9 ± 3.3 kJ/mol, respectively. Calculated α−β and β−γ transition boundaries were generally consistent with those determined by high-pressure experiments within the errors. Combining the measured ΔH°₂₉₈ and ΔS°₂₉₈ with selected data of in situ X-ray diffraction experiments at high pressure, the ΔH°₂₉₈ and ΔS°₂₉₈ of the α−β and β−γ transitions were optimized. Calculation using the optimized data tightly constrained the α−β and β−γ transition boundaries in the P, T space. The slope of α−β transition boundary is 3.1 MPa/K at 13.4 GPa and 1,400 K, and that of β−γ boundary 5.2 MPa/K at 18.7 GPa and 1,600 K. The post-spinel transition boundary of γ-Mg₂SiO₄ to MgSiO₃ perovskite plus MgO was also calculated, using the optimized data on γ-Mg₂SiO₄ and available enthalpy and entropy data on MgSiO₃ perovskite and MgO. The calculated post-spinel boundary with a Clapeyron slope of −2.6 ± 0.2 MPa/K is located at pressure consistent with the 660 km discontinuity, considering the error of the thermodynamic data.     

Ab initio calculated geometries and charge distributions for H4SiO4 and H6Si2O7 compared with experimental values for silicates and siloxanes

Ab initio STO-3G molecular orbital theory has been used to calculate energy-optimized Si-O bond lengths and angles for molecular orthosilicic and pyrosilicic acids. The resulting bond length for orthosilicic acid and the nonbridging bonds for pyrosilicic acid compare well with Si-OH bonds observed for a number of hydrated silicate minerals. Minimum energy Si-O bond lengths to the bridging oxygen of the pyrosilicic molecule show a close correspondence with bridging bond length data observed for the silica polymorphs and for gas phase and molecular crystal siloxanes when plotted against the SiOSi angle. In addition, the calculations show that the mean Si-O bond length of a silicate tetrahedron increases slightly as the SiOSi angle narrows. The close correspondence between the Si-O bond length and angle variations calculated for pyrosilicic acid and those observed for the silica polymorphs and siloxanes substantiates the suggestion that local bonding forces in solids are not very different from those in molecules and clusters consisting of the same atoms with the same coordination numbers. An extended basis calculation for H₄SiO₄ implies that there are about 0.6 electrons in the 3d-orbitals on Si. An analysis of bond overlap populations obtained from STO-3G* calculations for H₆Si₂O₇ indicates that Si-O bond length and SiOSi angle correlations may be ascribed to changes in the hybridization state of the bridging oxygen and (d – p) π-bonding involving all five of the 3d AO's of Si and the lone-pair AO's of the oxygen. Theoretical density difference maps calculated for H₆Si₂O₇ show a build-up of charge density between Si and O, with the peak-height charge densities of the nonbridging bonds exceeding those of the bridging bonds by about 0.05 e Å^?3. In addition, atomic charges (+1.3 and ?0.65) calculated for Si and O in a SiO₂ moiety of the low quartz structure conform reasonably well with the electroneutrality postulate and with experimental charges obtained from monopole and radial refinements of diffraction data recorded for low quartz and coesite.     

Heat capacity of γ-Fe2SiO4 between 5 and 303 K and derived thermodynamic properties

A multi-anvil device was used to synthesize 24 mg of pure γ-Fe₂SiO₄ crystals at 8.5 GPa and 1,273 K. The low-temperature heat capacity (C _p) of γ-Fe₂SiO₄ was measured between 5 and 303 K using the heat capacity option of a physical properties measurement system. The measured heat capacity data show a broad λ-transition at 11.8 K. The difference in the C _p between fayalite and γ-Fe₂SiO₄ is reduced as the temperature increases in the range of 50–300 K. The gap in C _p data between 300 and 350 K of γ-Fe₂SiO₄ is an impediment to calculation of a precise C _p equation above 298 K that can be used for phase equilibrium calculations at high temperatures and high pressures. The C _p and entropy of γ-Fe₂SiO₄ at standard temperature and pressure (S°₂₉₈) are 131.1 ± 0.6 and 140.2 ± 0.4 J mol⁻¹ K⁻¹, respectively. The Gibbs free energy at standard pressure and temperature (ΔG° _f,298) is calculated to be −1,369.3 ± 2.7 J mol⁻¹ based on the new entropy data. The phase boundary for the fayalite–γ-Fe₂SiO₄ transition at 298 K based on current thermodynamic data is located at 2.4 ± 0.6 GPa with a slope of 25.4 bars/K, consistent with extrapolated results of previous experimental studies.     

Theoretical reaction pathways for the formation of [Si(OH)5]1− and the deprotonation of orthosilicic acid in basic solution

Ab initio, molecular orbital calculations at the 6-31G¹ level including second-order Møller-Plesset electron correlation predict that the species [Si(OH)₅]¹⁻ is dynamically stable in a distorted trigonal bipyramid configuration. Reaction pathways for Si(OH)₄ + (OH)⁻ → [Si(OH)₅]¹⁻ → [(OH)₃SiO]¹⁻—H₂O are also calculated. The first reaction represents the formation of pentacoordinate Si from orthosilicic acid and hydroxide. The activation energy for adding a fifth Si-(OH)⁻ bond to the Si(OH)₄ molecule is ≈0.1 eV /molec (≈10kJ/mol). The second reaction is the deprotonation of the Si(OH)₄ which forms as a hydroxyl group leaves the [Si(OH)₅]¹⁻ molecule. Removal of a bond from this complex requires 0.9 eV/molecule (≈85 kJ/mol). Lengthening the Si—OH₂ distance results in the isolated molecules [(OH)₃SiO]¹⁻ + H₂O. This represents dehydration of the deprotonated orthosilicic acid.[Si(OH)₅]¹⁻ and [(OH)₃SiO]¹⁻- H₂O have the same energetic stability within the accuracy of these calculations. The potential energies of the isolated molecular systems [(OH)₃SiO]¹⁻+ H₂O and Si(OH)₄ + (OH)⁻ are considerably higher. These results suggest that [Si(OH)₅]¹⁻ may be a stable species or reaction intermediate in dissolution of silicate minerals in basic aqueous solutions.     

Thermochemical study of glasses in the system NaAlSi3O8-KAlSi3O8-Si4O8 and the join Na1.6Al1.6Si2.4O8-K1.6Al1.6Si2.4O8

Enthalpies of solution in 2PbO · B₂O₃ at 981 K have been measured for glasses in the system albite-orthoclase-silica and along the join Na_1.6Al_1.6Si_2.4O₈-K_1.6Al_1.6Si_2.4O₈. The join KAlSi₃O₈-Si₄O₈ shows zero heat of mixing similar to that found previously for NaAlSi₃O₈-Si₄O₈ glasses. Albite-orthoclase glasses show negative heats of mixing symmetric about Ab₅₀Or₅₀ (W_n = ? 2.4 ± 0.8 kcal). Negative heats of (Na, K) mixing are also found at $\text{Si}\text{(Si + Al)}\text{= 0.6.}$ Ternary excess enthalpies of mixing in the glassy system Ab-Or-4Q are positive but rarely exceed 1 kcal mol^?1.Using earlier studies of the thermodynamic properties of the crystals, the present calorimetric data and the “two-lattice” entropy model, the albite-orthoclase phase diagram is calculated in good agreement with experimental data. Attempts to calculate albite-silica and orthoclase-silica phase diagrams reveal complexities probably related to significant (but unknown) mutual solid solubility between cristobalite and alkali feldspar and to the very small heat and entropy of fusion of SiO₂.     

Ab initio valence force field calculations for quartz

We have derived valence force constants for the tetrahedral SiO₄ unit and the inter-tetrahedral SiOSi linkage from previous ab initio molecular orbital calculations on H₄SiO₄ and H₆Si₂O₇ using a split-valence polarized Gaussian basis set (6-31G^*), and used these to calculate the infrared and Raman active vibrational modes of α-quartz. The calculation gives frequencies approximately 15% greater than experiment, as expected from harmonic force constants obtained at this level of Hartree-Fock theory, but the calculation gives the correct distribution of modes within each frequency range. Calculated 28–30 Si and 16–18 O isotope shifts and pressure shifts to 6 GPa are also in reasonable agreement with experiment. We have also used our ab initio force field to calculate the vibrational spectrum for β-quartz. The results suggest either that inclusion of a torsional force constant is important for determining the stability of this high temperature polymorph, or that the β-quartz has a disordered structure with lower symmetry (P6₂) domains, as suggested by earlier diffraction studies.     

Thermal expansion of hydrated six-coordinate silicon in thaumasite,Ca<Subscript>3</Subscript>Si(OH)<Subscript>6</Subscript>(CO<Subscript>3</Subscript>)(SO<Subscript>4</Subscript>)·12H<Subscript>2</Subscript>O

Thaumasite, Ca₃Si(OH)₆(CO₃)(SO₄)12H₂O, occurs as a low-temperature secondary alteration phase in mafic igneous and metamorphic rocks, and is recognized as a product and indicator of sulfate attack in Portland cement. It is also the only mineral known to contain silicon in six-coordination with hydroxyl (OH)^? that is stable at ambient P–T conditions. Thermal expansion of the various components of this unusual structure has been determined from single-crystal X-ray structure refinements of natural thaumasite at 130 and 298 K. No phase transitions were observed over this temperature range. Cell parameters at room temperature are: a= 11.0538(6) Å, c=10.4111(8) Å and V=1101.67(10) ų, and were measured at intervals of about 50 K between 130 and 298 K, resulting in mean axial and volumetric coefficients of thermal expansion (×10^?5K^?1); α _a =1.7(1), α _c =2.1(2), and α _V =5.6(2). Although the unit cell and ^VIIICaO₈ polyhedra show significant positive thermal expansion over this temperature range, the silicate octahedron, sulfate tetrahedron, and carbonate group show zero or negative thermal expansion, with α _V (^VISiO₆) = ?0.6 ± 1.1, α _V (^IVSO₄)=?5.8 ± 1.4, and α _R (C–O)= 0.0 ± 1.8 (×10^?5 K^?1). Most of the thermal expansion is accommodated by lengthening of the R(O...O) hydrogen bond distances by on average 5σ, although the hydrogen bonds involving hydroxyl sites on ^VISi expand twice as much as those on molecular water, causing the [Ca₃Si(OH)₆(H₂O)₁₂]⁴⁺ columns to expand in diameter more than they move apart over this temperature range. The average Si–OH bond length of the six-coordinated Si atom 〈R(^VISi–OH)〉 in thaumasite is 1.783(1) Å, being about 0.02 Å (?20σ) shorter than ^VISi–OH in the dense hydrous magnesium silicate, phase D, MgSi₂H₂O₆.     

H2O diffusion models in rhyolitic melt with new high pressure data

Water diffusion in silicate melts is important for understanding bubble growth in magma, magma degassing and eruption dynamics of volcanos. Previous studies have made significant progress on water diffusion in silicate melts, especially rhyolitic melt. However, the pressure dependence of H₂O diffusion is not constrained satisfactorily. We investigated H₂O diffusion in rhyolitic melt at 0.95–1.9 GPa and 407–1629 °C, and 0.2–5.2 wt.% total water (H₂O_t) content with the diffusion-couple method in a piston-cylinder apparatus. Compared to previous data at 0.1–500 MPa, H₂O diffusivity is smaller at higher pressures, indicating a negative pressure effect. This pressure effect is more pronounced at low temperatures. Assuming H₂O diffusion in rhyolitic melt is controlled by the mobility of molecular H₂O (H₂O_m), the diffusivity of H₂O_m (D_H2Om) at H₂O_t ≤ 7.7 wt.%, 403–1629 °C, and ≤ 1.9 GPa is given by

D_H2Om=D₀exp(aX),