14 research outputs found

    Ewald methods for polarizable surfaces with application to hydroxylation and hydrogen bonding on the (012) and (001) surfaces of alpha-Fe2O3

    Full text link
    We present a clear and rigorous derivation of the Ewald-like method for calculation of the electrostatic energy of the systems infinitely periodic in two-dimensions and of finite size in the third dimension (slabs) which is significantly faster than existing methods. Molecular dynamics simulations using the transferable/polarizable model by Rustad et al. were applied to study the surface relaxation of the nonhydroxylated, hydroxylated, and solvated surfaces of alpha-Fe2O3 (hematite). We find that our nonhydroxylated structures and energies are in good agreement with previous LDA calculations on alpha-alumina by Manassidis et al. [Surf. Sci. Lett. 285, L517, 1993]. Using the results of molecular dynamics simulations of solvated interfaces, we define end-member hydroxylated-hydrated states for the surfaces which are used in energy minimization calculations. We find that hydration has a small effect on the surface structure, but that hydroxylation has a significant effect. Our calculations, both for gas-phase and solution-phase adsorption, predict a greater amount of hydroxylation for the (012) surface than for the (001) surface. Our simulations also indicate the presence of four-fold coordinated iron ions on the (001) surface.Comment: 23 pages, REVTeX (LaTeX), 8 figures not included, e-mail to [email protected], paper accepted in Surface Scienc

    GMIN: A computerized chemical equilibrium model using a constrained minimization of the Gibbs free energy

    No full text
    This document is a user's manual and technical reference for the computerized chemical equilibrium model GMIN. GMIN calculates the chemical composition of systems composed of pure solid phases, solid-solution phases, gas phases, adsorbed phases, and the aqueous phase. In the aqueous phase model, the excess solution free energy is modeled by using the equations developed by PITZER and his coworkers, which are valid to high ionic strengths. The Davies equation can also be used. Activity coefficients for nonideal soild-solution phases are calculated using parameters of polynomial expansion in mole fraction of the excess free energy of mixing. The free energy of adsorbed phase species is described by the triple-layer site-binding model. The mathematical algorithm incorporated into GMIN is based upon a constrained minimization of the Gibbs free energy. This algorithm is numerically stable and reliably converges to a free energy minimum. The data base for GMIN contains all standard chemical potentials and Pitzer ion-interaction parameters necessary to model the system Na-K-Ca-Mg-H-Cl-SO{sub 4}-CO{sub 2}-B(OH){sub 4}-H{sub 2}0 at 25{degrees}C

    Interfacial reactivity of radionuclides: emerging paradigms from molecular-level observations

    Get PDF
    Over the past few decades an increasing array of molecular-level analytical probes has provided new detailed insight into mineral and radionuclide interfacial reactivity in subsurface environments. This capability has not only helped change the way mineral surface reactivity is studied but also how field scale contaminant migration problems are addressed and ultimately resolved. Here we review relatively new interfacial reactivity paradigms and assess their implications for future research directions. Specific examples include understanding the following: the role of site-to-site electron conduction at mineral surfaces and through bulk mineral phases and the effects of local chemical environment on the stability of intermediate species in oxidation-reduction reactions and the importance of mechanistic reaction pathways for defining possible reaction products and thermodynamic driving force. The discussion also includes examples of how detailed molecular/microscopic characterization of field samples has changed the way complex contaminant migration problems are conceptualized and modeled

    Predictive geochemical modeling of interactions between uranium-mill-tailings solutions and sediments in a flow-through system: model formulations and preliminary results

    No full text
    An equilibrium thermodynamic conceptual model consisting of minerals and solid phases was developed to represent a soil column. A computer program was used as a tool to solve the system of mathematical equations imposed by the conceptual chemical model. The combined conceptual model and computer program were used to predict aqueous phase compositions of effluent solutions from permeability cells packed with geologic materials and percolated with uranium mill tailings solutions. Initial calculations of ion speciation and mineral solubility and our understanding of the chemical processes occurring in the modeled system were used to select solid phases for inclusion in the conceptual model. The modeling predictions were compared to the analytically determined column effluent concentrations. Hypotheses were formed, based on modeling predictions and laboratory evaluations, as to the probable mechanisms controlling the migration of selected contaminants. An assemblage of minerals and other solid phases could be used to predict the concentrations of several of the macro constituents (e.g., Ca, SO/sub 4/, Al, Fe, and Mn) but could not be used to predict trace element concentrations. These modeling conclusions are applicable to situations where uranium mill tailings solutions of low pH and high total dissolved solids encounter either clay liners or natural geologic materials that contain inherent acid neutralizing capacities. 116 references, 22 figures, 6 tables

    MININR: a geochemical computer program for inclusion in water flow models - an application study

    No full text
    MININR is a reduced form of the computer program MINTEQ which calculates equilibrium precipitation/dissolution of solid phases, aqueous speciation, adsorption, and gas phase equilibrium. The user-oriented features in MINTEQ were removed to reduce the size and increase the computational speed. MININR closely resembles the MINEQL computer program developed by Westall (1976). The main differences between MININR and MINEQL involve modifications to accept an initial starting mass of solid and necessary changes for linking with a water flow model. MININR in combination with a simple water flow model which considers only dilution was applied to a laboratory column packed with retorted oil shale and percolated with distilled water. Experimental and preliminary model simulation results are presented for the constituents K/sup +/, Na/sup +/, SO/sub 4//sup 2 -/, Mg/sup 2 +/, Ca/sup 2 +/, CO/sub 3//sup 2 -/ and pH
    corecore