Acid dissociation mechanisms of Si(OH)4 and Al(H2O)6 in aqueous solution

Silicic acid and the hexa-aqua of Al³⁺ are fundamental model aqueous species of chemical importance in nature. In order to investigate their hydroxyl dissociation mechanisms, Car-Parrinello molecular dynamics (CPMD) simulations were carried out, which allow treating the solutes and solvents on the same footing. The method of constraint was employed to trigger the reactions by taking coordination number as the reaction coordinate and the thermodynamic integration was used to obtain the free-energy profiles. The approximate transition states were located and the reactant and product states were also characterized. The free-energy changes of dissociation are found about 15.0 kcal/mol and 7.7 kcal/mol for silicic acid and Al-aqua, respectively. From the simulation results, the first pKas were calculated by using two approaches, which are based on the pristine thermodynamic relation and the RDF (radial distribution function)-free energy relation, respectively. Because of more uncertainties involved in the RDF way, it is suggested that the pristine way should be favored, which shows an error margin of 1 pKa unit. This study provides an encouraging basis for applying the present methodology to predict acidity constants of those groups that are difficult to measure experimentally.     

Stable boron isotope fractionation between dissolved B(OH)3 and B(OH)4

The stable boron isotope ratio (¹¹B/¹⁰B) in marine carbonates is used as a paleo-pH recorder and is one of the most promising paleo-carbonate chemistry proxies. Understanding the thermodynamic basis of the proxy is of fundamental importance, including knowledge on the equilibrium fractionation factor between dissolved boric acid, B(OH)₃, and borate ion, B(OH)₄⁻ (, hereafter α_(B3-B4)). However, this factor has hitherto not been determined experimentally and a theoretically calculated value (Kakihana and Kotaka, 1977, hereafter KK77) has therefore been widely used. I examine the calculations underlying this value. Using the same spectroscopic data and methods as KK77, I calculate the same α_(B3−B4) = 1.0193 at 300 K. Unfortunately, it turns out that in general the result is sensitive to the experimentally determined vibrational frequencies and the theoretical methods used to calculate the molecular forces. Using analytical techniques and ab initio molecular orbital theory, the outcome for α_(B3-B4) varies between ∼1.020 and ∼1.050 at 300 K. However, several arguments suggest that α_(B3-B4) ? 1.030. Measured isotopic shifts in various ¹⁰B-, ²D-, and ¹⁸O-labeled isotopomers do not provide a constraint on stable boron isotope fractionation. I conclude that in order to anchor the fundamentals of the boron pH proxy, experimental work is required. The critics of the boron pH proxy should note, however, that uncertainties in α_(B3-B4) do not bias pH reconstructions provided that organism-specific calibrations are used.     

K-H3OFe3(SO4)2(OH)6铁矾的XRD研究

根据X射线衍射(XRD)分析发现: A Fe₃(SO₄)₂(OH)₆(A=K⁺、H₃O⁺)系列铁钒的XRD数据十分相近,难以用XRD区别,需通过能谱(EDS)辅助分析,才能区分此类铁矾。另外,此类铁矾的003和107面网间距d随K⁺含量增大而增大,且呈一元三次方程的关系;而033和220面网间距d随K⁺含量增大而减小,呈一元二次方程的关系。对该现象从铁矾晶体结构方面进行解释:K⁺、H₃O⁺离子位于较大空隙中,且沿着Z轴方向排列,当K⁺、H₃O⁺离子之间相互替换时,会导致该铁矾晶体结构在Z轴方向有较明显的变化。     

The kinetics of oxygen exchange between the GeO4Al12(OH)24(OH2)12(aq) molecule and aqueous solutions

We report rates of oxygen exchange with bulk solution for an aqueous complex, ^IVGeO₄Al₁₂(OH)₂₄(OH₂)₁₂⁸⁺(aq) (GeAl₁₂), that is similar in structure to both the ^IVAlO₄Al₁₂(OH)₂₄(OH₂)₁₂⁷⁺(aq) (Al₁₃) and ^IVGaO₄Al₁₂(OH)₂₄(OH₂)₁₂⁷⁺(aq) (GaAl₁₂) molecules studied previously. All of these molecules have ε-Keggin-like structures, but in the GeAl₁₂ molecule, occupancy of the central tetrahedral metal site by Ge(IV) results in a molecular charge of +8, rather than +7, as in the Al₁₃ and GaAl₁₂. Rates of exchange between oxygen sites in this molecule and bulk solution were measured over a temperature range of 274.5 to 289.5 K and 2.95 < pH < 4.58 using ¹⁷O-NMR.Apparent rate parameters for exchange of the bound water molecules (η-OH₂) are k_ex²⁹⁸ = 200 (±100) s⁻¹, ΔH^‡ = 46 (±8) kJ · mol⁻¹, and ΔS^‡ = −46 (±24) J · mol⁻¹ K⁻¹ and are similar to those we measured previously for the GaAl₁₂ and Al₁₃ complexes. In contrast to the Al₁₃ and GaAl₁₂ molecules, we observe a small but significant pH dependence on rates of solvolysis that is not yet fully constrained and that indicates a contribution from the partly deprotonated GeAl₁₂ species.The two topologically distinct μ₂-OH sites in the GeAl₁₂ molecule exchange at greatly differing rates. The more labile set of μ₂-OH sites in the GeAl₁₂ molecule exchange at a rate that is faster than can be measured by the ¹⁷O-NMR isotopic-equilibration technique. The second set of μ₂-OH sites have rate parameters of k_ex²⁹⁸ = 6.6 (±0.2) · 10⁻⁴ s⁻¹, ΔH^‡ = 82 (±2) kJ · mol⁻¹, and ΔS^‡ = −29 (±7) J · mol⁻¹ · K⁻¹, corresponding to exchanges ≈40 and ≈1550 times, respectively, more rapid than the less labile μ₂-OH sites in the Al₁₃ and GaAl₁₂ molecules. We find evidence of nearly first-order pH dependence on the rate of exchange of this μ₂-OH site with bulk solution for the GeAl₁₂ molecule, which contrasts with Al₁₃ and GaAl₁₂ molecules.     

Free energy of formation for green rust sodium sulphate (NaFe6Fe3(OH)18(SO4)2(s))

In a recent study, sulphate-bearing green rust (GR_SO4) was shown to incorporate Na⁺ in its structure (NaFe^II₆Fe^III₃(OH)₁₈(SO₄)_2(s); GR_Na,SO4). The compound was synthesised by aerial oxidation of Fe(OH)_2(s) in the presence of NaOH. This paper reports on its free energy of formation .Freshly synthesised GR_Na,SO4 was titrated with 0.5 M H₂SO₄ in an inert atmosphere at 25 °C, producing dissolved Fe²⁺ and magnetite or goethite. Solution concentrations, PHREEQC and the MINTEQ database were used to calculate reaction constants for the reactions:

     

Solubility of Fe-ettringite (Ca6[Fe(OH)6]2(SO4)3 · 26H2O)

The solubility of Fe-ettringite (Ca₆[Fe(OH)₆]₂(SO₄)₃ · 26H₂O) was measured in a series of precipitation and dissolution experiments at 20 °C and at pH-values between 11.0 and 14.0 using synthesised material. A time-series study showed that equilibrium was reached within 180 days of ageing. After equilibrating, the solid phases were analysed by XRD and TGA while the aqueous solutions were analysed by ICP-OES (calcium, sulphur) and ICP-MS (iron). Fe-ettringite was found to be stable up to pH 13.0. At higher pH-values Fe-monosulphate (Ca₄[Fe(OH)₆]₂(SO₄) · 6H₂O) and Fe-monocarbonate (Ca₄[Fe(OH)₆]₂(CO₃) · 6H₂O) are formed. The solubilities of these hydrates at 25 °C are:      

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.     

Theoretical exploration of the water exchange mechanism of the polyoxocation GaO4Al12(OH)24(H2O)12 in aqueous solution

Reaction pathways, solvent effects and energy barriers have been investigated for the water exchange of the polyoxocation GaO₄Al₁₂(OH)₂₄(H₂O)₁₂⁷⁺ (K-GaAl₁₂) in aqueous solution by means of supermolecule density functional theory calculations. In the proposed reaction pathway, the supermolecular reactant K-GaAl₁₂15H₂O first loses a water ligand to form an intermediate with a five-coordinated aluminum atom, and then the incoming water molecule in the second coordination sphere attacks the intermediate with a five-coordinated aluminum atom to produce the reaction product. Our calculated results indicate that the water exchange of K-GaAl₁₂ proceeds via a dissociative mechanism, and that the reverse reaction of Step II is the most favorable dissociative pathway, with a barrier height of 31.3 kJ mol⁻¹. The calculated transition-state rate for the favorable dissociative pathway is much larger than the experimental rate constant, but is close to the data calculated for Al₃₀ by molecular dynamics. The transmission coefficient was also predicted on the basis of both the calculated transition-state rate and the experimental rate. Our calculated results also indicate that both the explicit solvent effect and the bulk solvent effect have obvious effects on the barrier heights of the water exchange reaction of K-GaAl₁₂. By comparison, the water exchange mechanism for K-GaAl₁₂ was found to be more similar to that for mineral surfaces than that for monomeric aluminum species.     

Microscope-photometric methods for non-destructive Fe2+ – Fe3+ determination in chloritoids (Fe2+, Mn2+, Mg)2 (Al,Fe3+)4 Si2O10(OH)4

Six natural chloritoid crystals with Fe²⁺ and Fe³⁺ contents ranging from 4.15 to 12.81 and from 0.411 to 0.849g-atoms/l, respectively, as determined by means of microprobe and Mössbauer techniques, served as reference material to develop non-destructive microscope-spectrophotometric methods for quantitative Fe²⁺ – Fe³⁺ determinations in chloritoids from unpolarized spectra of (001) platelets. Fe²⁺ concentrations in g-atom/l can be obtained from [ [Fe³⁺]=C₁xD₁/t where D₁ = log₁₀(I_0/I at 28,000 cm^-1 and t=crystal thickness in cm; C₁ is a conttant that may be influenced somewhat by experimental conditions and is found to be 0.002289 with the experimental set-up used in this study. Fe²⁺ concentrations in g-atom/l can be obtained from [Fe²⁺]=C₁xD₁/D₁-C₃ with D₂=log₁₀(I₀/I) at 16,300 cm^?1 and constants C₄ = 45.36 and C₅ = 3.540. Due to the uncertainties in absorbance measurements, D₁ and D₂ and the thickness measurements, the accuracies are ±0.05 and ±0.15 g-atom/l for [Fe³⁺] and [Fe²⁺], respectively. The determinations may be carried out on chloritoid grains in normal thin sections with an areal resolution of ～10 μm.     

Structural change and mineralogical transformation mechanism of aluminum hydroxide gels from forced hydrolysis Al（Ⅲ） solutions containing AlO4Al12（OH）24（H2O）12^7＋ polyoxycation during aging

The structural change and mineralogy of Al gel during aging time were investigated by using spectroscopy techniques. The results indicated that: 1) the aggregation extent of solution-gel system increases with aging time, and the structure of amorphous gel becomes more short-ordered; 2) after six months, the gel formats nordstrandite and little gibbsite; 3) a marked decrease in the number of (Al-OH)_oh bands occurring at 610 cm⁻¹ and increase in the number of (Al-OH₂)_oh bands occurring at 555 cm⁻¹ indicate that the gel undergoes rearrangement-like process during aging; 4) the gel constantly contains Al-O tetrahedron of Keggin structure, but the signal peak occurring at ≈61×10⁻⁶ of ²⁷Al MAS NMR have a slight shift to downfield with aging time. A mineralogical transformation mechanism for hydrolysis Al(III) solution was proposed.     

Density functional studies of the structural characteristics, Al NMR chemical shifts and water-exchange reactions of Al30O8(OH)56(H2O)26(Al30) in aqueous solution

In this paper, the structure of the Al₃₀O₈(OH)₅₆(H₂O)₂₆¹⁸⁺(Al₃₀) polyoxocation in aqueous solution is investigated, including an exploration of its water-exchange reaction using a supramolecular model. Thirty-one solvent water molecules were explicitly included in the supramolecular model to approximate the influence of the solvent. The calculated results indicated that both the gas-phase and the supramolecular models could correctly reproduce the structure of the Al₃₀ polyoxocation, but the supramolecular model described the structure more accurately. Using the supramolecular model, we calculated the ²⁷Al NMR chemical shifts of various aluminum atoms using HF and GIAO methods, and they compared well to the chemical shifts determined experimentally. The water-exchange reaction of the Al₃₀ polyoxocation could not be simulated with the gas phase model because of a proton-transfer reaction that is induced by the highly positive charge of the Al₃₀ polyoxocation. However, the inclusion of an explicit second solvation sphere lowered the acidity of the coordinated water molecules and allowed simulation of the water exchange reaction.     

Solid-liquid equilibria of Mg(OH)2(cr) and Mg2(OH)3Cl·4H2O(cr) in the system Mg-Na-H-OH-Cl-H2O at 25°C

The solubility of crystalline Mg(OH)₂(cr) was determined by measuring the equilibrium H⁺ concentration in water, 0.01-2.7 m MgCl₂, 0.1-5.6 m NaCl, and in mixtures of 0.5 and 5.0 m NaCl containing 0.01-0.05 m MgCl₂. In MgCl₂ solutions above 2 molal, magnesium hydroxide converted into hydrated magnesium oxychloride. The solid-liquid equilibrium of Mg₂(OH)₃Cl·4H₂O(cr) was studied in 2.1-5.2 m MgCl₂. Using known ion interaction Pitzer coefficients for the system Mg-Na-H-OH-Cl-H₂O (25°C), the following equilibrium constants at I = 0 are calculated:

     

Dissolution of jarosite [KFe3(SO4)2(OH)6] at pH 2 and 8: Insights from batch experiments and computational modelling

Jarosite [KFe₃(SO₄)₂(OH)₆] is a mineral that is common in acidic, sulphate-rich environments, such as acid sulphate soils derived from pyrite-bearing sediments, weathering zones of sulphide ore deposits and acid mine or acid rock drainage (ARD/AMD) sites. The structure of jarosite is based on linear tetrahedral-octahedral-tetrahedral (T-O-T) sheets, made up from slightly distorted FeO₆ octahedra and SO₄ tetrahedra. Batch dissolution experiments carried out on synthetic jarosite at pH 2, to mimic environments affected by ARD/AMD, and at pH 8, to simulate ARD/AMD environments recently remediated with slaked lime (Ca(OH)₂), suggest first order dissolution kinetics. Both dissolution reactions are incongruent, as revealed by non-ideal dissolution of the parent solids and, in the case of the pH 8 dissolution, because a secondary goethite precipitate forms on the surface of the dissolving jarosite grains. The pH 2 dissolution yields only aqueous K, Fe, and SO₄. Aqueous, residual solid, and computational modelling of the jarosite structure and surfaces using the GULP and MARVIN codes, respectively, show for the first time that there is selective dissolution of the A- and T-sites, which contain K and SO₄, respectively, relative to Fe, which is located deep within the T-O-T jarosite structure. These results have implications for the chemistry of ARD/AMD waters, and for understanding reaction pathways of ARD/AMD mineral dissolution.     

The kinetics of dissolution of zircon (ZrSiO4)

The dissolution of amorphous zircon obeys first order kinetics, and, over long periods of time decreases due to precipitation of new phases from solution. The dissolution reaction has an activation energy of 12.8 K cal/mole. Crystalline zircon has low solubility in aqueous solution and methods capable of detecting Zr and Si in solution at concentrations of the order of 0.001 ppb are needed to successfully evaluate the kinetics of dissolution of crystalline zircon.