A multi-scale assessment of the impact of salinity on the desorption of chromate from hematite: Sea level rise implications.

The solubility and transport of Cr(VI) is primarily controlled by adsorption-desorption reactions at the surfaces of soil minerals such as iron oxides. Environmental properties such as pH, ionic strength, and ion competition are expected to affect the mobility and fate of Cr(VI). Sea level rise (SLR), and consequent seawater intrusion, is creating a new biogeochemical soil environment at coastal margins, potentially impacting Cr(VI) retention at contaminated sites. We employed in-situ ATR-FTIR spectroscopy and DFT calculations to investigate at the molecular level the adsorption of Cr(VI) on the hematite surface and its desorption by sulfate, as a function of pH and ionic strength. We further used a batch experiment to assess Cr(VI) desorption at varying artificial seawater (ASW) concentrations. IR results demonstrate the complexity of Cr(VI) adsorption, showing a combination of monodentate inner-sphere complexation at high pH and dichromate outer-sphere (∼75%) at low pH. The Cr(VI)-complexes exhibited desorption induced by increasing pH values (58% of desorption) and sulfate competition (∼40% desorption). ASW desorbed ∼20% more Cr(VI), even at just 1% concentration. Our findings provide insight into Cr(VI)-adsorption complexation that controls the retention and remobilization of Cr(VI) on Fe-oxide minerals. The results point to an elevated risk of Cr(VI) mobilization in contaminated soils affected by SLR.

Distinguishing different surface interactions for nucleotides adsorbed onto hematite and goethite particle surfaces through ATR-FTIR spectroscopy and DFT calculations.

Geochemical interfaces can impact the fate and transport of aqueous species in the environment including biomolecules. In this study, we investigate the surface chemistry of adsorbed nucleotides on two different minerals, hematite and goethite, using infrared spectroscopy and density functional theory (DFT) calculations. Attenuated total reflectance-Fourier transform infrared spectroscopy is used to probe the adsorption of deoxyadenosine monophosphate (dAMP), deoxyguanosine monophosphate (dGMP), deoxycytidine monophosphate (dCMP), and deoxythymidine monophosphate (dTMP) onto either hematite or goethite particle surfaces. The results show preferential adsorption of the phosphate group to either surface. Remarkably, surface adsorption of the four nucleotides onto either hematite or goethite have nearly identical experimental spectra in the phosphate region (900 to 1200 cm-1) for each mineral surface yet are distinctly different between the two minerals, suggesting differences in binding of these nucleotides to the two mineral surfaces. The experimental absorption frequencies in the phosphate region were compared to DFT calculations for nucleotides adsorbed through the phosphate group to binuclear clusters in either a monodentate or bidentate bridging coordination. Although the quality of the fits suggests that both binding modes may be present, the relative amounts differ on the two surfaces with preferential bonding suggested to be monodentate coordination on hematite and bidentate bridging on goethite. Possible reasons for these differences are discussed.

Oxide- and Silicate-Water Interfaces and Their Roles in Technology and the Environment.

Interfacial reactions drive all elemental cycling on Earth and play pivotal roles in human activities such as agriculture, water purification, energy production and storage, environmental contaminant remediation, and nuclear waste repository management. The onset of the 21st century marked the beginning of a more detailed understanding of mineral aqueous interfaces enabled by advances in techniques that use tunable high-flux focused ultrafast laser and X-ray sources to provide near-atomic measurement resolution, as well as by nanofabrication approaches that enable transmission electron microscopy in a liquid cell. This leap into atomic- and nanometer-scale measurements has uncovered scale-dependent phenomena whose reaction thermodynamics, kinetics, and pathways deviate from previous observations made on larger systems. A second key advance is new experimental evidence for what scientists hypothesized but could not test previously, namely, interfacial chemical reactions are frequently driven by "anomalies" or "non-idealities" such as defects, nanoconfinement, and other nontypical chemical structures. Third, progress in computational chemistry has yielded new insights that allow a move beyond simple schematics, leading to a molecular model of these complex interfaces. In combination with surface-sensitive measurements, we have gained knowledge of the interfacial structure and dynamics, including the underlying solid surface and the immediately adjacent water and aqueous ions, enabling a better definition of what constitutes the oxide- and silicate-water interfaces. This critical review discusses how science progresses from understanding ideal solid-water interfaces to more realistic systems, focusing on accomplishments in the last 20 years and identifying challenges and future opportunities for the community to address. We anticipate that the next 20 years will focus on understanding and predicting dynamic transient and reactive structures over greater spatial and temporal ranges as well as systems of greater structural and chemical complexity. Closer collaborations of theoretical and experimental experts across disciplines will continue to be critical to achieving this great aspiration.

Grazing-incidence diffraction reveals cellulose and pectin organization in hydrated plant primary cell wall.

The primary cell wall is highly hydrated in its native state, yet many structural studies have been conducted on dried samples. Here, we use grazing-incidence wide-angle X-ray scattering (GIWAXS) with a humidity chamber, which enhances scattering and the signal-to-noise ratio while keeping outer onion epidermal peels hydrated, to examine cell wall properties. GIWAXS of hydrated and dried onion reveals that the cellulose ([Formula: see text]) lattice spacing decreases slightly upon drying, while the (200) lattice parameters are unchanged. Additionally, the ([Formula: see text]) diffraction intensity increases relative to (200). Density functional theory models of hydrated and dry cellulose microfibrils corroborate changes in crystalline properties upon drying. GIWAXS also reveals a peak that we attribute to pectin chain aggregation. We speculate that dehydration perturbs the hydrogen bonding network within cellulose crystals and collapses the pectin network without affecting the lateral distribution of pectin chain aggregates.

A perspective on iron (Fe) in the atmosphere: air quality, climate, and the ocean.

As scientists engage in research motivated by climate change and the impacts of pollution on air, water, and human health, we increasingly recognize the need for the scientific community to improve communication and knowledge exchange across disciplines to address pressing and outstanding research questions holistically. Our professional paths have crossed because our research activities focus on the chemical reactivity of Fe-containing minerals in air and water, and at the air-sea interface. (Photo)chemical reactions driven by Fe can take place at the surface of the particles/droplets or within the condensed phase. The extent and rates of these reactions are influenced by water content and biogeochemical activity ubiquitous in these systems. One of these reactions is the production of reactive oxygen species (ROS) that cause damage to respiratory organs. Another is that the reactivity of Fe and organics in aerosol particles alter surficial physicochemical properties that impact aerosol-radiation and aerosol-cloud interactions. Also, upon deposition, aerosol particles influence ocean biogeochemical processes because micronutrients such as Fe or toxic elements such as copper become bioavailable. We provide a perspective on these topics and future research directions on the reactivity of Fe in atmospheric aerosol systems, from sources to short- and long-term impacts at the sinks with emphasis on needs to enhance the predictive power of atmospheric and ocean models.

Microbial Denitrification: Active Site and Reaction Path Models Predict New Isotopic Fingerprints.

The study of isotopic fingerprints in nitrate (δ15N, δ18O, Δ17O) has enabled pivotal insights into the global nitrogen cycle and revealed new knowledge gaps. Measuring populations of isotopic homologs of intact NO3 - ions (isotopologues) shows promise to advance the understanding of nitrogen cycling processes; however, we need new theory and predictions to guide laboratory experiments and field studies. We investigated the hypothesis that the isotopic composition of the residual nitrate pool is controlled by the N-O bond-breaking step in Nar dissimilatory nitrate reductase using molecular models of the enzyme active sites and associated kinetic isotope effects (KIEs). We integrated the molecular model results into reaction path models representing the reduction of nitrate under either closed-system or steady-state conditions. The predicted intrinsic KIE (15ε and 18ε) of the Nar active site matches observed fractionations in both culture and environmental studies. This is what would be expected if the isotopic composition of marine nitrate were controlled by dissimilatory nitrate reduction by Nar. For a closed system, the molecular models predict a pronounced negative 15N-18O clumping anomaly in residual nitrate. This signal could encode information about the amount of nitrate consumed in a closed system and thus constrain initial nitrate concentration and its isotopic composition. Similar clumped isotope anomalies can potentially be used to distinguish whether a system is open or closed to new nitrate addition. These mechanistic predictions can be tested and refined in combination with emerging ESI-Orbitrap measurements.

Connecting Thermodynamics of Alkali Ion Exchange on the Quartz (101) Surface with Density Functional Theory Calculations.

Periodic plane-wave density functional theory (DFT) calculations were performed on the α-quartz (SiO2) (101) surface to model exchange of adsorbed Li+ and either Na+, K+, or Rb+ in inner- and outer-sphere adsorbed, and aqueous configurations, which are charge-balanced with 2 Cl-. SiO- or SiOH groups represented the adsorption surface sites. The SiO- models included 58 H2O and 2 H3O+ molecules to approximate an aqueous environment, whereas the SiOH models had 59 H2O and 1 H3O+ molecules. The goal of this work is to calculate the heats of exchange for these alkali ions and to compare the results with those measured by flow microcalorimetry to ascertain the most probable mechanisms for these cations exchanging on the α-quartz (101) surface. Energy minimizations of each alkali ion adsorbed as outer-sphere complexes on SiOH surface sites, and as inner- and outer-sphere complexes on SiO- surface sites, were used to determine the energy of exchange (ΔEex) with Li+ for comparison with experimentally determined ΔHex values. Here, we present a novel method for calculating ΔEex using the difference in energies of geometry-optimized end member models. The aqueous and surface structures produced are similar to those observed experimentally. Although the trend for the calculated ΔEex values is consistent with those from the heats of exchange measured experimentally, the magnitude of our modeled ΔEex results is significantly larger than select experimental data from the literature [Peng, L. Zeta-Potentials and Enthalpy Changes in the Process of Electrostatic Self-Assembly of Cations on Silica Surface. Powder Technol. 2009, 193(1), 46-49]; we discuss the reasons for this discrepancy herein. The relative energy differences of the various configurations modeled have implications for the measurements of the surface charge via potentiometric titrations due to the more active role of alkali cations in quartz surface chemistry that have been previously considered as inert background electrolytes.

Adsorption of Organic Acids and Phosphate to an Iron (Oxyhydr)oxide Mineral: A Combined Experimental and Density Functional Theory Study.

The interaction of soil organic matter with mineral surfaces is a critical reaction involved in many ecosystem services, including stabilization of organic matter in the terrestrial carbon pool and bioavailability of plant nutrients. Using model organic acids typically present in soil solutions, this study couples laboratory adsorption studies with density functional theory (DFT) to provide physical insights into the nature of the chemical bonding between carboxylate functional groups and a model FeOOH cluster. Topological determination of electron density at bond critical points using quantum theory of atoms in molecules (QTAIM) analysis revealed that the presence of multiple bonding paths between the organic acid and the FeOOH cluster is essential in determining the competitive adsorption of organic acids and phosphate for FeOOH surface adsorption sites. The electron density and Laplacian parameter values from QTAIM indicated that the primary carboxylate-FeOOH bond was more ionic than covalent in nature. The experimental and computational results provide molecular-level evidence of the important role of electrostatic forces in the bonding between carboxylic acids and Fe-hydroxides. This knowledge may assist in the formulation of management studies to meet the challenges of maintaining ecosystems services in the face of a changing climate.

In Situ and Real-Time ATR-FTIR Temperature-Dependent Adsorption Kinetics Coupled with DFT Calculations of Dimethylarsinate and Arsenate on Hematite Nanoparticles.

Temperature-dependent kinetic studies of the adsorption of critical pollutants onto reactive components in soils and removal technologies provide invaluable rate information and mechanistic insight. Using attenuated total internal reflection Fourier transform infrared spectroscopy, we collected in situ spectra as a function of time, concentration, and temperature in the range of 5-50 °C (278-323 K) for the adsorption of arsenate (iAs) and dimethylarsinate (DMA) on hematite nanoparticles at pH 7. These experimental data were modeled with density functional theory (DFT) calculations on the energy barriers between surface complexes. The Langmuir adsorption kinetic model was used to extract values of the fast (<5 min) and slow (6-10 min) observed adsorption rate, initial rate constants of adsorption and desorption, Arrhenius parameters, effective activation energies (ΔEa), and pre-exponential factors (A). The trend in the kinetic parameters correlated with the type of surface complexes that iAs and DMA form, which are mostly bidentate binuclear compared to a mix of outer sphere and monodentate, respectively. The observed initial adsorption rates were found to be more sensitive to changes in the aqueous concentration of the arsenicals than slow rates. On average, iAs adsorbs 2.5× faster and desorbs 4× slower than dimethylarsinate (DMA). The ΔEa and A values for the adsorption of iAs bidentate complexes are statistically higher than those extracted for outer-sphere DMA by a factor of 3. The DFT results on adsorption energies and ΔEa barriers are consistent with the experimental data and provide a mechanistic explanation for the low ΔEa values observed. The presence of defect sites with under-coordinated Fe atoms or exchangeable surface water (i.e., Fe-OH2 groups) lowers activation barriers of adsorption. These results suggest that increasing organic substitutions on arsenate at the expense of As-O bonds decreases the effective energy barrier for complex formation and lowers the number of collisional orientations that result in binding to the hematite surface.

Evaluating Computational Chemistry Methods for Isotopic Fractionation between CO2(g) and H2O(g).

Quantum mechanical calculations can be useful in predicting equilibrium isotopic fractionations of geochemical reactions. However, these computational chemistry methods vary widely in their effectiveness in the prediction of various physical observables. Most studies employing the approach known as density functional theory (DFT) to model these observable quantities focus on predictive accuracy for energetics and geometries. In this study, several density functionals are evaluated against experimental bond lengths, harmonic vibrational frequencies, frequency shifts upon isotopic substitution, and 18O/16O isotopic fractionation between CO2(g) and H2O(g). Successful prediction of harmonic vibrational frequencies strongly correlates with successful prediction of isotopic fractionation, despite the possible introduction of errors by the harmonic approximation. Harmonic experimental frequencies, not anharmonic ones, must be used when comparing spectra and when predicting isotope fractionation. The B3LYP and X3LYP functionals perform more accurately in the evaluation of both harmonic vibrational frequencies and isotopic fractionation factors using the 6-311+G(d,p) and 6-311++G(2d,p) basis sets, achieving fractionation factor errors of 0.2-0.6‰ at 25 °C out of a total fractionation of 51‰. Error cancellation between vibrational frequencies and the harmonic approximation is crucial to their success. The above combination of exchange-correlation functionals and basis sets also well predicts the vibrational properties of interacting CO2 and H2O molecules, suggesting that they may be applicable to more complex geochemical reactions involving C and O isotopic fractionations.

Gibbsite (100) and Kaolinite (100) Sorption of Cadmium(II): A Density Functional Theory and XANES Study of Structures and Energies.

Due to the potential toxicity of cadmium (Cd2+) and its presence in various waste products found in the environment, it is necessary to develop methods to attenuate and remediate Cd2+ waste. Sorption of Cd2+ to mineral surfaces is a potential route to accomplish this goal. This work focused on improving our molecular-scale understanding of the chemistry of Cd2+ interactions with gibbsite and kaolinite mineral surfaces. Plane-wave density functional theory (DFT) energy minimization calculations and molecular dynamics simulations were used to study the adsorption energies and the nature of the bonds between Cd2+ and the mineral surfaces for possible inner- and outer-sphere surface complexes. Models resulting from the DFT calculations were used to calculate theoretical XANES spectra that were compared with experimental Cd LIII XANES of aqueous Cd2+ as a proxy for outer-sphere Cd2+ hydrated complexes associated with the mineral surfaces. These studies suggest that Cd2+ would favorably bond to the (100) surfaces of both kaolinite and gibbsite through a bidentate mononuclear interaction. However, the results indicate that mixtures of surface complexes on these minerals are likely.

Simulations of Cellulose Synthesis Initiation and Termination in Bacteria.

The processivity of cellulose synthesis in bacterial cellulose synthase (CESA) was investigated using molecular dynamics simulations and the hybrid quantum mechanics and molecular mechanics approach. Our results suggested that cellulose synthesis in bacterial CESA can be initiated with H2O molecules. The chain length or degree of polymerization (DOP) of the product cellulose is related to the affinity of the cellulose chain to the transmembrane tunnel of the enzyme. This opens up the possibility of generating mutants that would produce cellulose chains with desired chain lengths that could have applications in the biofuel and textile fields that depend on the DOP of cellulose chains.

Quantum Calculations on Plant Cell Wall Component Interactions.

Density functional theory calculations were performed to assess the relative interaction energies of plant cell wall components: cellulose, xylan, lignin and pectin. Monomeric and tetramer linear molecules were allowed to interact in four different configurations for each pair of compounds. The M05-2X exchange-correlation functional which implicitly accounts for short- and mid-range dispersion was compared against MP2 and RI-MP2 to assess the reliability of the former for modeling van der Waals forces between these PCW components. Solvation effects were examined by modeling the interactions in the gas phase, in explicit H2O, and in polarized continuum models (PCM) of solvation. PCMs were used to represent water, methanol, and chloroform. The results predict the relative ranges of each type of interaction and when specific configurations will be strongly preferred. Structures and energies are useful as a basis for testing classical force fields and as guidance for coarse-grained models of PCWs.

Cyclodextrin-enhanced 1,4-dioxane treatment kinetics with TCE and 1,1,1-TCA using aqueous ozone.

Enhanced reactivity of aqueous ozone (O3) with hydroxypropyl-ß-cyclodextrin (HPßCD) and its impact on relative reactivity of O3 with contaminants were evaluated herein. Oxidation kinetics of 1,4-dioxane, trichloroethylene (TCE), and 1,1,1-trichloroethane (TCA) using O3 in single and multiple contaminant systems, with and without HPßCD, were quantified. 1,4-Dioxane decay rate constants for O3 in the presence of HPßCD increased compared to those without HPßCD. Density functional theory molecular modeling confirmed that formation of ternary complexes with HPßCD, O3, and contaminant increased reactivity by increasing reactant proximity and through additional reactivity within the HPßCD cavity. In the presence of chlorinated co-contaminants, the oxidation rate constant of 1,4-dioxane was enhanced. Use of HPßCD enabled O3 reactivity within the HPßCD cavity and enhanced 1,4-dioxane treatment rates without inhibition in the presence of TCE, TCA, and radical scavengers including NaCl and bicarbonate. Micro-environmental chemistry within HPßCD inclusion cavities mediated contaminant oxidation reactions with increased reaction specificity.

Kinetic analysis of cellulose synthase of Gluconacetobacter hansenii in whole cells and in purified form.

The Gram-negative bacterium, Gluconacetobacter hansenii, has been long studied and is a model for cellulose synthesis. It produces cellulose, using the enzyme AcsA-AcsB, of exceptionally high crystallinity in comparison to the cellulose of higher plants. We determined the rate of cellulose synthesis in whole cells measured as moles of glucose incorporated into cellulose per second per mole of enzyme. This was determined by quantifying the rate of cellulose synthesis (over a short time span, such that the enzyme concentration is not changing due to cell growth) and the amount of enzyme in the whole cell by quantitative western blotting. We found that the whole cell rate of 24 s-1 is much faster than the kcat, measured from steady-state kinetic analysis, of 1.7 s-1. Our whole cell rates are consistent with previous studies using microscopy. We postulate that the rationale for this difference is the presence of an alternative in vivo priming mechanism. This in turn can increase the rate of initiation, which we previously postulated to be the rate-limiting step in catalysis.

The Shape of Native Plant Cellulose Microfibrils.

Determining the shape of plant cellulose microfibrils is critical for understanding plant cell wall molecular architecture and conversion of cellulose into biofuels. Only recently has it been determined that these cellulose microfibrils are composed of 18 cellulose chains rather than 36 polymers arranged in a diamond-shaped pattern. This study uses density functional theory calculations to model three possible habits for the 18-chain microfibril and compares the calculated energies, structures, 13C NMR chemical shifts and WAXS diffractograms of each to evaluate which shape is most probable. Each model is capable of reproducing experimentally-observed data to some extent, but based on relative theoretical energies and reasonable reproduction of all variables considered, a microfibril based on 5 layers in a 34443 arrangement is predicted to be the most probable. A habit based on a 234432 arrangement is slightly less favored, and a 6 × 3 arrangement is considered improbable.

Initiation, Elongation, and Termination of Bacterial Cellulose Synthesis.

Cellulose is the major component of the plant cell wall and composed of ß-linked glucose units. Use of cellulose is greatly impacted by its physical properties, which are dominated by the number of individual cellulose strand within each fiber and the average length of each strand. Our work described herein provides a complete mechanism for cellulose synthase accounting for its processivity and mechanism of initiation. Using ionic liquids and gel permeation chromatography, we obtain kinetic constants for initiation, elongation, and termination (release of the cellulose strand from the enzyme) for two bacterial cellulose synthases (Gluconacetobacter hansenii and Rhodobacter sphaeroides). Our results show that initiation of synthesis is primer-independent. After initiation, the enzyme undergoes multiple cycles of elongation until the strand is released. The rate of elongation is much faster than that of steady-state turnover. Elongation requires cyclic addition of glucose (from uridine diphosphate-glucose) and then strand translocation by one glucose unit. Translocations greater than one glucose unit result in termination requiring reinitiation. The rate of the strand release, relative to the rate of elongation, determines the processivity of the enzyme. This mechanism and the measured rate constants were supported by kinetic simulation. With the experimentally determined rate constants, we are able to simulate steady-state kinetics and mimic the size distribution of the product. Thus, our results provide for the first time a mechanism for cellulose synthase that accounts for initiation, elongation, and termination.

Density functional theory modeling of chromate adsorption onto ferrihydrite nanoparticles.

Density functional theory (DFT) calculations were performed on a model of a ferrihydrite nanoparticle interacting with chromate ([Formula: see text]) in water. Two configurations each of monodentate and bidentate adsorbed chromate as well as an outer-sphere and a dissolved bichromate ([Formula: see text]) were simulated. In addition to the 3-D periodic planewave DFT models, molecular clusters were extracted from the energy-minimized structures. Calculated interatomic distances from the periodic and cluster models compare favorably with Extended X-ray Absorption Fine Structure spectroscopy values, with larger discrepancies seen for the clusters due to over-relaxation of the model substrate. Relative potential energies were derived from the periodic models and Gibbs free energies from the cluster models. A key result is that the bidentate binuclear configuration is the lowest in potential energy in the periodic models followed by the outer-sphere complex. This result is consistent with observations of the predominance of bidentate chromate adsorption on ferrihydrite under conditions of high surface coverage (Johnston Environ Sci Technol 46:5851-5858, 2012). Cluster models were also used to perform frequency analyses for comparison with observed ATR FTIR spectra. Calculated frequencies on monodentate, bidentate binuclear, and outer-sphere complexes each have infrared (IR)-active modes consistent with experiment. Inconsistencies between the thermodynamic predictions and the IR-frequency analysis suggest that the 3-D periodic models are not capturing key components of the system that influence the adsorption equilibria under varying conditions of pH, ionic strength and electrolyte composition. Model equilibration via molecular dynamics (MD) simulations is necessary to escape metastable states created during DFT energy minimizations based on the initial classical force field MD-derived starting configurations.

Fourier-transform infrared spectroscopy (FTIR) analysis of triclinic and hexagonal birnessites.

The characterization of birnessite structures is particularly challenging for poorly crystalline materials of biogenic origin, and a determination of the relative concentrations of triclinic and hexagonal birnessite in a mixed assemblage has typically required synchrotron-based spectroscopy and diffraction approaches. In this study, Fourier-transform infrared spectroscopy (FTIR) is demonstrated to be capable of differentiating synthetic triclinic Na-birnessite and synthetic hexagonal H-birnessite. Furthermore, IR spectral deconvolution of peaks resulting from MnO lattice vibrations between 400 and 750cm-1 yield results comparable to those obtained by linear combination fitting of synchrotron X-ray absorption fine structure (EXAFS) data when applied to known mixtures of triclinic and hexagonal birnessites. Density functional theory (DFT) calculations suggest that an infrared absorbance peak at ~1628cm-1 may be related to OH vibrations near vacancy sites. The integrated intensity of this peak may show sensitivity to vacancy concentrations in the Mn octahedral sheet for different birnessites.