Kristian W. Paul

ATR-FTIR and density functional theory study of the structures, energetics, and vibrational spectra of phosphate adsorbed onto goethite.

Periodic plane-wave density functional theory (DFT) and molecular cluster hybrid molecular orbital-DFT (MO-DFT) calculations were performed on models of phosphate surface complexes on the (100), (010), (001), (101), and (210) surfaces of α-FeOOH (goethite). Binding energies of monodentate and bidentate HPO(4)(2-) surface complexes were compared to H(2)PO(4)(-) outer-sphere complexes. Both the average potential energies from DFT molecular dynamics (DFT-MD) simulations and energy minimizations were used to estimate adsorption energies for each configuration. Molecular clusters were extracted from the energy-minimized structures of the periodic systems and subjected to energy reminimization and frequency analysis with MO-DFT. The modeled P-O and P---Fe distances were consistent with EXAFS data for the arsenate oxyanion that is an analog of phosphate, and the interatomic distances predicted by the clusters were similar to those of the periodic models. Calculated vibrational frequencies from these clusters were then correlated with observed infrared bands. Configurations that resulted in favorable adsorption energies were also found to produce theoretical vibrational frequencies that correlated well with experiment. The relative stability of monodentate versus bidentate configurations was a function of the goethite surface under consideration. Overall, our results show that phosphate adsorption onto goethite occurs as a variety of surface complexes depending on the habit of the mineral (i.e., surfaces present) and solution pH. Previous IR spectroscopic studies may have been difficult to interpret because the observed spectra averaged the structural properties of three or more configurations on any given sample with multiple surfaces.

Periodic density functional theory calculations of bulk and the (010) surface of goethite.

BackgroundGoethite is a common and reactive mineral in the environment. The transport of contaminants and anaerobic respiration of microbes are significantly affected by adsorption and reduction reactions involving goethite. An understanding of the mineral-water interface of goethite is critical for determining the molecular-scale mechanisms of adsorption and reduction reactions. In this study, periodic density functional theory (DFT) calculations were performed on the mineral goethite and its (010) surface, using the Vienna Ab Initio Simulation Package (VASP).ResultsCalculations of the bulk mineral structure accurately reproduced the observed crystal structure and vibrational frequencies, suggesting that this computational methodology was suitable for modeling the goethite-water interface. Energy-minimized structures of bare, hydrated (one H2O layer) and solvated (three H2O layers) (010) surfaces were calculated for 1 × 1 and 3 × 3 unit cell slabs. A good correlation between the calculated and observed vibrational frequencies was found for the 1 × 1 solvated surface. However, differences between the 1 × 1 and 3 × 3 slab calculations indicated that larger models may be necessary to simulate the relaxation of water at the interface. Comparison of two hydrated surfaces with molecularly and dissociatively adsorbed H2O showed a significantly lower potential energy for the former.ConclusionSurface Fe-O and (Fe)O-H bond lengths are reported that may be useful in surface complexation models (SCM) of the goethite (010) surface. These bond lengths were found to change significantly as a function of solvation (i.e., addition of two extra H2O layers above the surface), indicating that this parameter should be carefully considered in future SCM studies of metal oxide-water interfaces.

Unimolecular decomposition of 5-aminotetrazole and its tautomer 5-iminotetrazole: new insight from isopotential searching.

Aminotetrazole compounds have become attractive ingredients in gas generating compositions, solid rocket propellants, and green pyrotechnics. Therefore, a fundamental understanding of their thermal decomposition mechanisms and thermodynamics is of great interest. In this study, the specular reflection isopotential searching method was used to investigate the unimolecular decomposition mechanisms of 5-iminotetrazole (5-ITZ), 1H-5-aminotetrazole (1H-5-ATZ), and 2H-5-aminotetrazole (2H-5-ATZ). Subsequent thermochemical analysis of the unimolecular decomposition pathways was performed at the CCSD(T)/aug-cc-pVTZ//B3LYP/6-311++G(3df,3pd) level of theory. Based upon the relative reaction barriers predicted in this study, the initial gaseous products of 5-ITZ unimolecular decomposition are HN(3) and NH(2)CN (calculated activation barrier equal to 199.5 kJ/mol). On the other hand, the initial gaseous products of 1H-5-ATZ and 2H-5-ATZ unimolecular decomposition are predicted to be N(2) and metastable CH(3)N(3) (calculated activation barriers equal to 169.2 and 153.7 kJ/mol, respectively). These predicted unimolecular decomposition products and activation barriers are in excellent agreement with thermal decomposition experiments performed by Lesnikovich et al. [Lesnikovich, A. I.; Ivashkevich, O. A.; Levchik, S. V.; Balabanovich, A. I.; Gaponik, P. N.; Kulak, A. A. Thermochim. Acta 2002, 388, 233], in which the apparent activation barriers were measured to be approximately 200 and 150 kJ/mol, respectively, for 5-ITZ and 1H-5-ATZ/2H-5-ATZ.

Environmental Fate and Transport of Energetic Materials

This project strives to improve our understanding of the environmental behavior of energetic materials (EM). A thorough understanding of how these EM interact with soil and water is expected to ultimately lead to improved remediation strategies. The immediate goals of the project are to predict a priori the chemical interactions of energetic materials with model soils; and to predict a priori the decomposition reactions of EM and resultant breakdown products. Adsorption research will involve periodic density functional theory (DFT) calculations, using plane-wave/pseudopotential codes. Decomposition studies will use the Specular Reflection Isopotential Searching algorithm of Irikura et al.