The kinetics of iodide oxidation by the manganese oxide mineral birnessite

The kinetics of iodide (I⁻) and molecular iodine (I₂) oxidation by the manganese oxide mineral birnessite (δ-MnO₂) was investigated over the pH range 4.5-6.25. I⁻ oxidation to iodate proceeded as a two-step reaction through an I₂ intermediate. The rate of the reaction varied with both pH and birnessite concentration, with faster oxidation occurring at lower pH and higher birnessite concentration. The disappearance of I⁻ from solution was first order with respect to I⁻ concentration, pH, and birnessite concentration, such that −d[I⁻]/dt = k[I⁻][H⁺][MnO₂], where k, the third order rate constant, is equal to 1.08 ± 0.06 × 10⁷ M⁻² h⁻¹. The data are consistent with the formation of an inner sphere I⁻ surface complex as the first step of the reaction, and the adsorption of I⁻ exhibited significant pH dependence. Both I₂, and to a lesser extent, sorbed to birnessite. The results indicate that iodine transport in mildly acidic groundwater systems may not be conservative. Because of the higher adsorption of the oxidized I species I₂ and , as well as the biophilic nature of I₂, redox transformations of iodine must be taken into account when predicting I transport in aquifers and watersheds.     

Kinetic analysis of cation exchange in birnessite using time-resolved synchrotron X-ray diffraction

In this study, we applied time-resolved synchrotron X-ray diffraction (TRXRD) to develop kinetic models that test a proposed two-stage reaction pathway for cation exchange in birnessite. These represent the first rate equations calculated for cation exchange in layered manganates. Our previous work has shown that the substitution of K, Cs, and Ba for interlayer Na in synthetic triclinic birnessite induces measurable changes in unit-cell parameters. New kinetic modeling of this crystallographic data supports our previously postulated two-stage reaction pathway for cation exchange, and we can correlate the kinetic steps with changes in crystal structure. In addition, the initial rates of cation exchange, R (ų min⁻¹), were determined from changes in unit-cell volume to follow these rate laws: R = 1.75[]^0.56, R = 41.1[]^1.10, R = 1.15[]^0.50. Thus, the exchange rates for Na in triclinic birnessite decreased in the order: Cs ? K > Ba. These results are likely a function of hydration energy differences of the cations and the preference of the solution phase for the more readily hydrated cation.     

The molecular mechanism of Mo isotope fractionation during adsorption to birnessite

Fractionation of Mo isotopes during adsorption to manganese oxides is a primary control on the global ocean Mo isotope budget. Previous attempts to explain what drives the surprisingly large isotope effect δ^97/95Mo_dissolved-δ^97/95Mo_adsorbed=1.8‰ have not successfully resolved the fractionation mechanism. New evidence from extended X-ray absorption fine structure analysis and density functional theory suggests that Mo forms a polymolybdate complex on the surfaces of experimental and natural samples. Mo in this polynuclear structure is in distorted octahedral coordination, while Mo remaining in solution is predominantly in tetrahedral coordination as . Our results indicate that the difference in coordination environment between dissolved Mo and adsorbed Mo is the cause of isotope fractionation. The molecular mechanism of metal isotope fractionation in this system should enable us to explain and possibly predict metal isotope effects in other systems where transition metals adsorb to mineral surfaces.     

Copper partitioning in a melt-vapor-brine-magnetite-pyrrhotite assemblage

The effect of sulfur on the partitioning of Cu in a melt-vapor-brine ± magnetite ± pyrrhotite assemblage has been quantified at 800 °C, 140 MPa, f_O2 = nickel-nickel oxide (NNO), logf_S2=-3.0 (i.e., on the magnetite-pyrrhotite curve at NNO), logf_H2S=-1.3 and logf_SO2=-1. All experiments were vapor + brine saturated. Vapor and brine fluid inclusions were trapped in silicate glass and self-healed quartz fractures. Vapor and brine are dominated by NaCl, KCl and HCl in the S-free runs and NaCl, KCl and FeCl₂ in S-bearing runs. Pyrrhotite served as the source of sulfur in S-bearing experiments. The composition of fluid inclusions, glass and crystals were quantified by laser-ablation inductively coupled plasma mass spectrometry. Major element, chlorine and sulfur concentrations in glass were quantified by using electron probe microanalysis. Calculated Nernst-type partition coefficients (±2σ) for Cu between melt-vapor, melt-brine and vapor-brine are , , and , respectively, in the S-free system. The partition coefficients (±2σ) for Cu between melt-vapor, melt-brine and vapor-brine are , , and , respectively, in the S-bearing system. Apparent equilibrium constants (±1σ) describing Cu and Na exchange between vapor and melt and brine and melt were also calculated. The values of are 34 ± 21 and 128 ± 29 in the S-free and S-bearing runs, respectively. The values of are 33 ± 22 and60 ± 5 in the S-free and S-bearing runs, respectively. The data presented here indicate that the presence of sulfur increases the mass transfer of Cu into vapor from silicate melt. Further, the nearly threefold increase in suggests that Cu may be transported as both a chloride and sulfide complex in magmatic vapor, in agreement with hypotheses based on data from natural systems. Most significantly, the data demonstrate that the presence of sulfur enhances the partitioning of Cu from melt into magmatic volatile phases.     

Oxygen isotope fractionation between apatite-bound carbonate and water determined from controlled experiments with synthetic apatites precipitated at 10-37 °C

The oxygen isotope fractionation between the structural carbonate of inorganically precipitated hydroxyapatite (HAP) and water was determined in the range 10-37 °C. Values of 1000 ln α() are linearly correlated with inverse temperature (K) according to the following equation: 1000 ln α() = 25.19 (±0.53)·T⁻¹ − 56.47 (±1.81) (R² = 0.998). This fractionation equation has a slightly steeper slope than those already established between calcite and water ( [O’Neil et al., 1969] and [Kim and O’Neil, 1997]) even though measured fractionations are of comparable amplitude in the temperature range of these experimental studies. It is consequently observed that the oxygen isotope fractionation between apatite carbonate and phosphate increases from about 7.5‰ up to 9.1‰ with decreasing temperature from 37 °C to 10 °C. A compilation of δ¹⁸O values of both phosphate and carbonate from modern mammal teeth and bones confirms that both variables are linearly correlated, despite a significant scattering up to 3.5‰, with a slope close to 1 and an intercept corresponding to a 1000 ln α() value of 8.1‰. This apparent fractionation factor is slightly higher or close to the fractionation factor expected to be in the range 7-8‰ at the body temperature of mammals.     

Molecular-scale mechanisms of distribution and isotopic fractionation of molybdenum between seawater and ferromanganese oxides

The distribution of Mo between seawater and marine ferromanganese oxides has great impacts on concentration and isotopic composition of Mo in modern oxic seawater. To reveal the adsorption chemistry of Mo to ferromanganese oxides, we performed (i) detailed structural analyses of Mo surface complexes on δ-MnO₂, ferrihydrite, and hydrogenetic ferromanganese oxides by L₃- and K-edge XAFS, and (ii) adsorption experiments of Mo to δ-MnO₂ and ferrihydrite over a wide range of pHs, ionic strengths, and Mo concentrations. XAFS analyses revealed that Mo forms distorted octahedral (Oh) inner-sphere complexes on δ-MnO₂ whereas it forms a tetrahedral (Td) outer-sphere complex on ferrihydrite. In the hydrogenetic ferromanganese oxides, the dominant host phase of Mo was revealed to be δ-MnO₂. These structural information are consistent with the macroscopic behaviors of Mo in adsorption experiments, and Mo concentration in modern oxic seawater can be explained by the equilibrium adsorption reaction on δ-MnO₂. In addition, the large isotopic fractionation of Mo between seawater and ferromanganese oxides detected in previous studies can be explained by the structural difference between and adsorbed species on the δ-MnO₂ phase in ferromanganese oxides. In contrast, smaller fractionation of Mo isotopes on ferrihydrite is due to little change in the Mo local structures during its adsorption to ferrihydrite.The structures of Mo species adsorbed on crystalline Fe (oxyhydr)oxides, goethite, and hematite were also investigated at pH 8 and I = 0.70 M (NaNO₃). Our XAFS analyses revealed that Mo forms inner-sphere complexes on both minerals: Td edge-sharing (46%) and Oh double corner-sharing (54%) for goethite, and Td double corner-sharing (14%) and Oh edge-sharing (86%) for hematite. These structural information, combined with those for amorphous ferrihydrite and δ-MnO₂, show the excellent correlation with the magnitude of adsorptive isotopic fractionation of Mo reported in previous studies: the proportion of Oh species or their magnitude of distortion in Mo surface complexes become larger in the order of ferrihydrite < goethite < hematite < δ-MnO₂, a trend identical to the magnitude of isotopic fractionation.Based on the comparison with previous reports for Mo surface species on various oxides, the chemical factors that affect Mo surface complex structures were also discussed. The hydrolysis constant of cation in oxides, log K_OH (or the acidity of the oxide surfaces, PZC) is well correlated with the mode of attachment (inner- or outer-sphere) of Mo surface complexes. Furthermore, the symmetric change in Mo species from Td to Oh is suggested to be driven by the formation of inner-sphere complexes on specific sites of the oxide surfaces.     

The partitioning of molybdenum(VI) between aqueous liquid and vapour at temperatures up to 370 °C

We have conducted experiments to evaluate the vapour-liquid fractionation of Mo(VI) in the system MoO₃-NH₃-H₂O at 300-370 °C and saturated vapour pressure, using a two-chamber autoclave that allows separate trapping of the vapour and liquid. The measured total Mo concentrations in each phase were used to calculate a distribution coefficient, , which increases as the density of the vapour approaches that of the liquid, and is greater than one for pH ? 4. Molybdenum speciation in the vapour is described by a single complex, MoO₃H₂O. By contrast, thermodynamic modeling of the distribution of Mo species in the liquid indicates that bimolybdate (HMoO₄⁻) is the dominant aqueous species at the conditions of our experiments, and that molybdate (MoO₄²⁻) and molybdic acid (H₂MoO₄⁰) are present in smaller quantities. As vapour-liquid fractionation occurs between neutral species, it is governed by the reaction H₂MoO₄⁰(aq) = MoO₃ · H₂O(g). Fractionation is therefore controlled by the concentration of H₂MoO₄⁰ in the liquid, which increases with increasing temperature and decreasing pH. Owing to the pH dependence of , it cannot be used to describe Mo fractionation in aqueous vapour-liquid systems with compositions different than those of this study. We have therefore calculated a composition-independent (Henry’s Law) constant, , for each experimental point, using the measured total Mo concentration in the vapour and the modeled concentration of H₂MoO₄⁰ in the liquid. This constant may be applied to aqueous vapour-liquid systems of known liquid composition to estimate the concentration of Mo in a vapour for which little chemical information is available, and thereby supplement the available fractionation data for natural porphyry-forming systems. The results of this study demonstrate that at conditions typical of natural porphyry ore-forming systems, a significant amount of molybdenum fractionates into the vapour over the liquid, and the vapour may transport quantities of Mo in excess of that in the liquid at pH conditions below those of the muscovite-microcline reaction boundary.     

Experimental constraints on the mobility of Rhenium in silicate liquids

The volatization of Rhenium (Re) from melts of natural basalt, dacite and a synthetic composition in the CaO-MgO-Al₂O₃-SiO₂ system has been investigated at 0.1 MPa and 1250-1350 °C over a range of fO₂ conditions from log fO₂ = −10 to −0.68. Experiments were conducted using open top Pt crucibles doped with Re and Yb. Analysis of quenched glasses by laser ablation-inductively coupled plasma mass spectrometry (LA-ICP-MS) normal to the melt/gas interface showed concentration profiles for Re, to which a semi-infinite one-dimensional diffusion model could be applied to extract diffusion coefficients (D). The results show Re diffusivity in basalt at 1300 °C in air is log D_Re = −7.2 ± 0.3 cm²/s and increases to log D_Re = −6.6 ± 0.3 cm²/s when trace amounts of Cl were added to the starting material. At fO₂ conditions below the nickel-nickel oxide (NNO) buffer Re diffusivity decreases to and to in dacitic melt. In the CMAS composition, . The diffusivity of Re is comparable to Ar and CO₂ in basalt at 500 MPa favoring its release as a volatile. Our results support the contention that subaerial degassing is the cause of lower Re concentrations in arc-type and ocean island basalts compared to mid-ocean ridge basalts.     

Precipitation and mixing properties of the “disordered” (Mn,Ca)CO3 solid solution

The enthalpy of mixing of the calcite-rhodochrosite (Ca,Mn)CO₃ solid solution was determined at 25 °C from calorimetric measurements of the enthalpy of precipitation of solids with different compositions. A detailed study of the broadening of powder X-ray diffraction peaks shows that most of the precipitates are compositionally homogeneous. All the experimental enthalpy of mixing (ΔH^m) values are positive and fit reasonably well (R² = 0.86) to a Guggenheim function of three terms:

     

The origin of Zn isotope fractionation in sulfides

Isotope fractionation of Zn between aqueous sulfide, chloride, and carbonate species (Zn²⁺, Zn(HS)₂, , , ZnS(HS)⁻, ZnCl⁺, ZnCl₂, , and ZnCO₃) was investigated using ab initio methods. Only little fractionation is found between the sulfide species, whereas carbonates are up to 1‰ heavier than the parent solution. At pH > 3 and under atmospheric-like CO₂ pressures, isotope fractionation of Zn sulfides precipitated from sulfidic solutions is affected by aqueous sulfide species and the δ⁶⁶Zn of sulfides reflect these in the parent solutions. Under high P_CO2 conditions, carbonate species become abundant. In high P_CO2 conditions of hydrothermal solutions, Zn precipitated as sulfides is isotopically nearly unfractionated with respect to a low-pH parent fluid. In contrast, negative δ⁶⁶Zn down to at least −0.6‰ can be expected in sulfides precipitated from solutions with pH > 9. Zinc isotopes in sulfides and rocks therefore represent a potential indicator of mid to high pH in ancient hydrothermal fluids.     

Arsenate uptake by calcite: Macroscopic and spectroscopic characterization of adsorption and incorporation mechanisms

Batch uptake experiments and X-ray element mapping and spectroscopic techniques were used to investigate As(V) (arsenate) uptake mechanisms by calcite, including adsorption and coprecipitation. Batch sorption experiments in calcite-equilibrated suspensions (pH 8.3; PCO₂ = 10^−3.5 atm) reveal rapid initial sorption to calcite, with sorption rate gradually decreasing with time as available sorption sites decrease. An As(V)-calcite sorption isotherm determined after 24 h equilibration exhibits Langmuir-like behavior up to As concentrations of 300 μM. Maximum distribution coefficient values (K_d), derived from a best fit to a Langmuir model, are ∼190 L kg⁻¹.Calcite single crystals grown in the presence of As(V) show well-developed rhombohedral morphology with characteristic growth hillocks on surfaces at low As(V) concentrations (?5 μM), but habit modification is evident at As(V) concentrations ?30 μM in the form of macrostep development preferentially on the − vicinal surfaces of growth hillocks. Micro-X-ray fluorescence element mapping of surfaces shows preferential incorporation of As in the − vicinal faces relative to + vicinals. EXAFS fit results for both adsorption and coprecipitation samples confirm that As occurs in the 5+ oxidation state in tetrahedral coordination with oxygen, i.e., as arsenate. For adsorption samples, As(V) forms inner-sphere surface complexes via corner-sharing with Ca octahedra. As(V) coprecipitated with calcite substitutes in carbonate sites but with As off-centered, as indicated by two Ca shells, and with likely disruption of local structure. The results indicate that As(V) interacts strongly with the calcite surface, similar to often-cited analog phosphate, and uptake can occur via both adsorption and coprecipitation reactions. Therefore, calcite may be effective for partial removal of dissolved arsenate from aquatic and soil systems.     

Formation of todorokite from vernadite in Ni-rich hemipelagic sediments

Todorokite is considered to form from vernadite in nature and commonly concentrates nickel. However, this mineralogical transformation has never been imaged nor explained mechanistically, and its effect on the uptake of nickel has never been quantified at the molecular-level. We have characterized these reactions at the macroscopic, microscopic, nanoscopic and atomic scales in a marine manganese concretion by combining transmission electron microscopy, electron and X-ray microprobe analysis, powder and micro X-ray diffraction, and Mn and Ni K-edge EXAFS spectroscopy. The concretion was collected during the Ticoflux II expedition near the Nicoya Peninsula, Costa Rica, and is representative of Mn deposits in hemipelagic sediments. It consists of 5 to 25 μm aggregates, shaped like sea-urchins, with a core of 7Å-vernadite (1.0 wt% Ni), a rim of 10Å-vernadite (3.8 wt% Ni), and an outermost region of todorokite fibers (1.9 wt% Ni) that extend outwards. The crystallites of 7Å-vernadite are single- to bi-layered, with hexagonal layer symmetry (a = b = 2.83 Å), and an average structural formula of . The crystallites of 10Å-vernadite contain 10 to 20 layers semi-coherently stacked in the ab plane and uniformly separated in the [0 0 1] direction by ∼9 Å due to the intercalation of hydrated Mg²⁺ cations. The average structural formula of 10Å-vernadite is if the layers contain vacancy sites, or alternately , if they contain Mn³⁺. The average formula of todorokite is .A genetic model is proposed based on combining these new data with previously published results. The thermodynamically unstable 7Å-vernadite transforms via dissolution-recrystallization to semi-ordered Mg-rich 10Å-vernadite. Nickel is released from dissolved biogenic silica or reduced organic matter, and taken up mainly in the Mn layer of 10Å-vernadite. Interlayer magnesium serves as a template to the further topotactic transformation of 10Å-vernadite to todorokite. The dimension of the todorokite tunnels in the [0 0 1] direction is uniform and determined by the size of the hydrated Mg²⁺ ion (8.6 Å). The tunnel dimension in the [1 0 0] direction depends on the density of Mg²⁺ in the interlayer and the superstructure of the phyllomanganate layer. If the parent phyllomanganate contains high amounts of Mg²⁺ (i.e., high layer charge), or Mn³⁺ and Mn⁴⁺ cations ordered following the Mn³⁺-Mn⁴⁺-Mn⁴⁺ sequence as in synthetic triclinic birnessite, then the tunnel dimension is ideally 3 × 3 octahedral chain widths in both crystallographic directions. Otherwise, the tunnel dimension is incoherent in the [1 0 0] direction (i.e., T(3,n) tunnel structure), as has been observed in all natural todorokites. Natural todorokite is defective because the precursor natural phyllomanganates either have a layer charge deficit below 0.33e per octahedral site, or rarely are triclinic birnessite. The abundance of Mg in seawater and its key role in converting phyllomanganate to tectomanganate with T(3,n) tunnel structure explain why todorokite is common in marine ferromanganese oxides, and seldom present in terrestrial environments. The topotactic phase transformation described here is the only known route to todorokite crystallization. This implies that all natural todorokites may be authigenic because they are formed in situ from a phyllomanganate precursor.     

Modeling the effects of bond environment on equilibrium iron isotope fractionation in ferric aquo-chloro complexes

Four or five sets of ab initio models, including Unrestricted Hartree Fock (UHF) and hybrid Density Functional Theory (DFT) are calculated for each species in a series of aqueous ferric aquo-chloro complexes: , , , FeCl₃(H₂O)₃, FeCl₃(H₂O)₂, , FeCl₅H₂O²⁻, , ) in order to determine the relative isotopic fractionation among the complexes, to compare the results of different models for the same complexes, to examine factors that influence the magnitude of the isotopic fractionation, and to compare bond-partner-driven fractionation with redox-driven fractionation.Relative to , all models show a nearly linear decrease in ⁵⁶Fe/⁵⁴Fe as the number of Cl⁻ ions per Fe³⁺ ion increases, with slopes of −0.8‰ to −1.0‰ per Cl⁻ at 20 °C. At 20 °C, 1000 ln β (β = ⁵⁶Fe/⁵⁴Fe reduced partition function ratio relative to a dissociated Fe atom) values range from 8.93‰ to 9.73‰ for , 8.04-9.12‰ for , 7.61-8.73‰ for , 7.14-8.25‰ for , and 3.09-4.41‰ for . The fractionation between and ranges from 1.5‰ to 2.6‰, depending on the model; this is comparable in magnitude to fractionation effects due to Fe³⁺/Fe²⁺ redox reactions. β values from the UHF models are consistently higher than those from the hybrid DFT models.Isotopic fractionation is shown to be sensitive to differences in ligand bond stiffness (above), coordination number, bond length, and the frequency of the asymmetric Fe-X stretching vibrational mode, as predicted by previous theoretical studies. Complexes with smaller coordination numbers have higher 1000 ln β (7.46‰, 5.25‰, and 3.48‰ for , ,, respectively, from the B3LYP/6-31G(d) model). Species with the same number of chlorides but fewer waters also show the effect of coordination number on 1000 ln β: (7.46‰ vs. 7.05‰ for FeCl₃(H₂O)₂ vs. FeCl₃(H₂O)₃ and 5.25‰ vs. 4.94‰ for vs. FeCl₅H₂O²⁻ with the B3LYP/6-31G(d) model). As more Fe-Cl bonds substitute for Fe-OH₂ bonds (with a resulting decrease in β), the lengths of the Fe-Cl bonds and the Fe-O bonds increase.Preliminary modeling of shows an Fe³⁺/Fe²⁺ fractionation of 3.2‰ for the B3LYP/6-31G(d) model, in agreement with previous studies. The addition of an explicit outer hydration sphere of 12 H₂O molecules to models of improves agreement with measured vibrational frequencies and bond lengths; 1000 ln β increases by 0.8-1.0‰. An additional hydration sphere around increases 1000 ln β by only 0.1‰.Isotopic fractionations predicted for this simple system imply that ligands present in an aqueous iron environment are potentially important drivers of fractionation, and suggest that significant fractionation effects are likely in other aqueous systems containing sulfides or organic ligands. Fractionation effects due to both speciation and redox must be considered when interpreting iron isotope fractionations in the geological record.     

Silicon isotopic fractionation during adsorption of aqueous monosilicic acid onto iron oxide

The quantification of silicon isotopic fractionation by biotic and abiotic processes contributes to the understanding of the Si continental cycle. In soils, light Si isotopes are selectively taken up by plants, and concentrate in secondary clay-sized minerals. Si can readily be retrieved from soil solution through the specific adsorption of monosilicic acid () by iron oxides. Here, we report on the Si-isotopic fractionation during adsorption on synthesized ferrihydrite and goethite in batch experiment series designed as function of time (0-504 h) and initial concentration (ic) of Si in solution (0.21-1.80 mM), at 20 °C, constant pH (5.5) and ionic strength (1 mM). At various contact times, the δ²⁹Si vs. NBS28 compositions were determined in selected solutions (ic = 0.64 and 1.06 mM Si) by MC-ICP-MS in dry plasma mode with external Mg doping with an average precision of ±0.08‰ (±2σ_SEM). Per oxide mass, ferrihydrite (74-86% of initial Si loading) adsorbed more Si than goethite (37-69%) after 504 h of contact over the range of initial Si concentration 0.42-1.80 mM. Measured against its initial composition (δ²⁹Si = +0.01 ± 0.04‰ (±2σ_SD)), the remaining solution was systematically enriched in ²⁹Si, reaching maximum δ²⁹Si values of +0.70 ± 0.07‰ for ferrihydrite and +0.50 ± 0.08‰ for goethite for ic 1.06 mM. The progressive ²⁹Si enrichment of the solution fitted better a Rayleigh distillation path than a steady state model. The fractionation factor ²⁹ε (±1σ_SD) was estimated at −0.54 ± 0.03‰ for ferrihydrite and −0.81 ± 0.12‰ for goethite. Our data imply that the sorption of onto synthetic iron oxides produced a distinct Si-isotopic fractionation for the two types of oxide but in the same order than that generated by Si uptake by plants and diatoms. They further suggest that the concentration of light Si isotopes in the clay fraction of soils is partly due to sorption onto secondary clay-sized iron oxides.     

Surface complexation of Cu on birnessite (δ-MnO2): Controls on Cu in the deep ocean

Hexagonal birnessite (δ-MnO₂) is a close analogue to the dominant phase in hydrogenetic marine ferromanganese crusts and nodules. These deposits contain ∼0.25 wt.% Cu which is believed to be scavenged from the overlying water column where Cu concentrations are near 0.1 μg/L. Here, we measured the sorption of Cu on δ-MnO₂ as a function of pH and surface loading. We characterized the nature of the Cu sorption complex at pH 4 and 8 using EXAFS spectroscopy and find that, at pH 4, Cu sorbs to birnessite by inner-sphere complexation on the {0 0 1} surface at sites above Mn vacancies to give a three to fourfold coordinated complex with 6 Mn neighbors at ∼3.4 Å. At pH 8, however, we find that some Cu has become structurally incorporated into the MnO₂ layer by occupying the vacancy sites to give 6 Mn neighbors at ∼2.91 Å. Density functional calculations on and clusters predict a threefold coordinated surface complex and show that the change from surface complexation to structural incorporation is a response to protonation of oxygens surrounding the vacancy site. Consequently, we propose that the transformation between sorption via surface complex and vacancy site occupancy should be reversible. By fitting the Cu sorption as a function of surface loading and pH to the formation of the observed and predicted surface complex, we developed a surface complexation model (in the basic Stern approximation) for the sorption of Cu onto birnessite. Using this model, we demonstrate that the concentration of inorganic Cu in the deep ocean should be several orders of magnitude lower than the observed total dissolved Cu. We propose that the observed total dissolved Cu concentration in the oceans reflects solubilization of Cu by microbially generated ligands.     

Influence of the enzyme dissimilatory sulfite reductase on stable isotope fractionation during sulfate reduction

The stable isotopes of sulfate are often used as a tool to assess bacterial sulfate reduction on the macro scale. However, the mechanisms of stable isotope fractionation of sulfur and oxygen at the enzymatic level are not yet fully understood. In batch experiments with water enriched in ¹⁸O we investigated the effect of different nitrite concentrations on sulfur isotope fractionation by Desulfovibrio desulfuricans.With increasing nitrite concentrations, we found sulfur isotope enrichment factors ranging from −11.2 ± 1.8‰ to −22.5 ± 3.2‰. Furthermore, the δ¹⁸O values in the remaining sulfate increased from approximately 50-120‰ when ¹⁸O-enriched water was supplied. Since ¹⁸O-exchange with ambient water does not take place in sulfate, but rather in intermediates of the sulfate reduction pathway (e.g. ), we suggest that nitrite affects the steady-state concentration and the extent of reoxidation of the metabolic intermediate sulfite to sulfate during sulfate reduction. Given that nitrite is known to inhibit the production of the enzyme dissimilatory sulfite reductase, our results suggest that the activity of the dissimilatory sulfite reductase regulates the kinetic isotope fractionation of sulfur and oxygen during bacterial sulfate reduction. Our novel results also imply that isotope fractionation during bacterial sulfate reduction strongly depends on the cell internal enzymatic regulation rather than on the physico-chemical features of the individual enzymes.     

Isotopic fractionations associated with phosphoric acid digestion of carbonate minerals: Insights from first-principles theoretical modeling and clumped isotope measurements

Phosphoric acid digestion has been used for oxygen- and carbon-isotope analysis of carbonate minerals since 1950, and was recently established as a method for carbonate ‘clumped isotope’ analysis. The CO₂ recovered from this reaction has an oxygen isotope composition substantially different from reactant carbonate, by an amount that varies with temperature of reaction and carbonate chemistry. Here, we present a theoretical model of the kinetic isotope effects associated with phosphoric acid digestion of carbonates, based on structural arguments that the key step in the reaction is disproportionation of H₂CO₃ reaction intermediary. We test that model against previous experimental constraints on the magnitudes and temperature dependences of these oxygen isotope fractionations, and against new experimental determinations of the fractionation of ¹³C-¹⁸O-containing isotopologues (‘clumped’ isotopic species). Our model predicts that the isotope fractionations associated with phosphoric acid digestion of carbonates at 25 °C are 10.72‰, 0.220‰, 0.137‰, 0.593‰ for, respectively, ¹⁸O/¹⁶O ratios (1000 lnα^∗) and three indices that measure proportions of multiply-substituted isotopologues . We also predict that oxygen isotope fractionations follow the mass dependence exponent, λ of 0.5281 (where ). These predictions compare favorably to independent experimental constraints for phosphoric acid digestion of calcite, including our new data for fractionations of ¹³C-¹⁸O bonds (the measured change in Δ₄₇ = 0.23‰) during phosphoric acid digestion of calcite at 25 °C.We have also attempted to evaluate the effect of carbonate cation compositions on phosphoric acid digestion fractionations using cluster models in which disproportionating H₂CO₃ interacts with adjacent cations. These models underestimate the magnitude of isotope fractionations and so must be regarded as unsucsessful, but do reproduce the general trend of variations and temperature dependences of oxygen isotope acid digestion fractionations among different carbonate minerals. We suggest these results present a useful starting point for future, more sophisticated models of the reacting carbonate/acid interface. Examinations of these theoretical predictions and available experimental data suggest cation radius is the most important factor governing the variations of isotope fractionation among different carbonate minerals. We predict a negative correlation between acid digestion fractionation of oxygen isotopes and of ¹³C-¹⁸O doubly-substituted isotopologues, and use this relationship to estimate the acid digestion fractionation of for different carbonate minerals. Combined with previous theoretical evaluations of ¹³C-¹⁸O clumping effects in carbonate minerals, this enables us to predict the temperature calibration relationship for different carbonate clumped isotope thermometers (witherite, calcite, aragonite, dolomite and magnesite), and to compare these predictions with available experimental determinations. The success of our models in capturing several of the features of isotope fractionation during acid digestion supports our hypothesis that phosphoric acid digestion of carbonate minerals involves disproportionation of transition state structures containing H₂CO₃.     

Oxygen and hydrogen isotope fractionation in serpentine-water and talc-water systems from 250 to 450 °C, 50 MPa

Oxygen and hydrogen isotope fractionation factors in the talc-water and serpentine-water systems have been determined by laboratory experiment from 250 to 450 °C at 50 MPa using the partial exchange technique. Talc was synthesized from brucite + quartz, resulting in nearly 100% exchange during reaction at 350 and 450 °C. For serpentine, D-H exchange was much more rapid than ¹⁸O-¹⁶O exchange when natural chrysotile fibers were employed in the initial charge. In experiments with lizardite as the starting charge, recrystallization to chrysotile enhanced the rate of ¹⁸O-¹⁶O exchange with the coexisting aqueous phase. Oxygen isotope fractionation factors in both the talc-water and serpentine-water systems decrease with increasing temperature and can be described from 250 to 450 °C by the relationships: 1000 ln  = 11.70 × 10⁶/T² − 25.49 × 10³/T + 12.48 and 1000 ln  = 3.49 × 10⁶/T² − 9.48 where T is temperature in Kelvin. Over the same temperature interval at 50 MPa, talc-water D-H fractionation is only weakly dependent on temperature, similar to brucite and chlorite, and can be described by the equation: 1000 ln = 10.88 × 10⁶/T² − 41.52 × 10³/T + 5.61 where T is temperature in Kelvin. Our D-H serpentine-water fractionation factors calibrated by experiment decrease with temperature and form a consistent trend with fractionation factors derived from lower temperature field calibrations. By regression of these data, we have refined and extended the D-H fractionation curve from 25 to 450 °C, 50 MPa as follows: 1000 ln  = 3.436 × 10⁶/T² − 34.736 × 10³/T + 21.67 where T is temperature in Kelvin. These new data should improve the application of D-H and ¹⁸O-¹⁶O isotopes to constrain the temperature and origin of hydrothermal fluids responsible for serpentine formation in a variety of geologic settings.