Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BATTELLE”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Marcus Theory of Ion-Pairing

We present a theory for ion pair dissociation and association, motivated by the concepts of the Marcus theory of electron transfer. Despite the extensive research on ion-pairing in many chemical and biological processes, much can be learned from the exploration of collective reaction coordinates. To this end, we explore two reaction coordinates, ion pair distance and coordination number. The study of the correlation between these reaction coordinates provides a new insight into the mechanism and kinetics of ion pair dissociation and association in water. The potential of mean force on these 2D-surfaces computed from molecular dynamics simulations of different monovalent ion pairs reveal a Marcus-like mechanism for ion-pairing: Water molecules rearrange forming an activated coordination state prior to ion pair dissociation or association, followed by relaxation of the coordination state due to further water rearrangement. Like Marcus theory, we find the existence of an inverted region where the transition rates are slower with increasing exergonicity. This study provides a new perspective for the future investigations of ion-pairing and transport. SR, CJM, and GKS were supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences. MDB was supported by MS3 (Materials Synthesis and Simulation Across Scales) Initiative, a Laboratory Directed Research and Development Program at Pacific Northwest National Laboratory (PNNL). The research was performed using PNNL Institutional Computing. PNNL is a multi-program national laboratory operated by Battelle for the U.S. Department of Energy.

Roy, Santanu↗
DOE OSTI2020

NWPEsSe: an Adaptive-Learning Global Optimization Algorithm for Nanosized Cluster Systems

Global optimization constitutes an important and fundamental problem in theoretical studies in many chemical fields, such as catalysis, materials or separations problems. In this paper, a novel algorithm has been developed for the global optimization of large systems including neat and ligated clusters in gas phase, and supported clusters in periodic boundary conditions. The method is based on an updated artificial bee colony (ABC) algorithm method, that allows for adaptive-learning during the search process. The new algorithm is tested against four classes of systems of diverse chemical nature: gas phase Au_55, ligated Au_8^(2+), Au_8 supported on graphene oxide and defected rutile, and a large cluster assembly ?[Co?_6 Te_8 (PEt_3 )_6][C_60 ]_n, with sizes ranging between 1 to 3 nm and containing up to 1300 atoms. Reliable global minima (GMs) are obtained for all cases, either confirming published data or reporting new lower energy structures. The algorithm and interface to other codes in the form of an independent program, Northwest Potential Energy Search Engine (NWPEsSe), is freely available and it provides a powerful and efficient approach for global optimization of nanosized cluster systems. The work described in this publication was performed at Pacific Northwest National Laboratory (PNNL), which is operated by Battelle for the United States Department of Energy (DOE) under Contract DE-AC05-76RL0180. J. Z. and V.-A. G. acknowledge support from DOE, Office of Science, Office of Basic Energy Sci-ences, Chemical, Geological and Biological Sciences Division and computing resources from PNNL’s Research Computing Facility and the National Energy Research Scientific Computing Center.

Zhang, Jun↗
DOE OSTI2021

Cryogenic Vibrationally Resolved Photoelectron Spectroscopy of OH–(H2O): Confirmation of Multidimensional Franck-Condon Simulation Results for the Transition State of the OH + H2O Reaction

We present a transition state spectroscopic study of the OH + H2O reaction using the experimental technique of cryogenic Negative Ion Photoelectron Spectroscopy (NIPES). The recorded NIPE spectrum at 193 nm exhibits multiple vibrational progressions that include excitations to the shared H atom antisymmetric stretching mode with an interval of 0.32 eV as well as other progressions, mainly involving the H bending and O…O symmetric stretching modes. The Vertical Detachment Energy (VDE) was measured at 3.53 eV, whereas an upper limit for the Adiabatic Detachment Energy (ADE) was estimated at 2.90 eV. These values are in excellent agreement with the theoretically computed values of 3.51 eV and 2.83 eV respectively, obtained at the CCSD(T)/aug-cc-pV5Z level of theory. The recorded NIPE spectrum is in very good agreement when compared to the one recently reported from four-dimensional Franck-Condon simulations, in which a similar spectral profile was predicted. Besides observing the ground state, a charge-transfer excited state in the form of [OH–(H2O)+] is identified with a relative energy of 1.39 eV, well matching the previous prediction of 1.36 eV. This work was supported by U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences, and performed using EMSL, a national scientific user facility sponsored by DOE’s Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory, which is operated by Battelle Memorial Institute for the DOE. The theoretical calculations were conducted on EMSL’s “Cascade” Supercomputer. This research also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231.

Cao, Wenjin↗
DOE OSTI2021

Photoelectron Spectroscopy and Theoretical Study on Monosolvated Cyanate Analogue Clusters ECX-· Sol (ECX- = NCSe-, AsCSe- and AsCS-; Sol = H2O, CH3CN)

Six monosolvated cyanate analogue clusters ECX-?Sol (ECX– = NCSe-, AsCSe-, and AsCS-; Sol = H2O and CH3CN) were investigated using negative ion photoelectron spectroscopy (NIPES). NIPES experiments show that these clusters possess overall similar spectra compared to their respective isolated ECX- anions but shift to higher electron binding energy with solvent CH3CN stabilizing the excess electrons slightly more than H2O. For the ECX-·H2O series, vertical detachment energies and their increments relative to the bare species are measured to be 3.700 / 0.370, 3.085 / 0.425, 3.085 / 0.430 eV for NCSe-, AsCSe- and AsCS-, respectively, while the corresponding values in ECX-·CH3CN series are 3.835 / 0.505, 3.145 / 0.475, and 3.135 / 0.480 eV. Ab initio electronic structure calculations indicate the excess charges located at the terminal N and Se atoms in NCSe- and migrated to the central C atom in AsCSe- and AsCS-. For NCSe-, the solvation is driven by the interactions with the two negatively charged terminal ends, while for AsCSe- and AsCS-, the solvation revolves around the interactions with the central C atom, where all the excess negative charge is concentrated. Two nearly degenerate isomers for NCSe-·H2O are identified, one forming a single strong N…H-O hydrogen bond (HB) and the other featuring a bidentate HB with two hydroxyl H atoms pointing to N and Se ends. In contrast, the negative central C atom in AsCSe-/AsCS- allows the formation of a bifurcated HB with H2O. Similar effects are observed for the acetonitrile case, in which the three H atoms of methyl group interact with the two negatively charged terminal ends in NCSe-, while prefer binding to the central negative carbon atom in AsCSe-/AsCS-. The different binding motifs derived in this work may suggest distinctly different solvation properties in NCSe- versus AsCSe-/AsCS-, with the former anion leading to asymmetric solvation at the N end of the solute, while the latter species creating more ‘isotropic’ solvation around the central C equatorial plane. This work was supported by U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences, and performed using EMSL, a national scientific user facility sponsored by DOE’s Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory, which is operated by Battelle Memorial Institute for the DOE. The theoretical calculations were conducted on EMSL’s “Cascade” Supercomputer.

Yuan, Qinqin↗
DOE OSTI2015

Snapshots of Proton Accommodation at a Microscopic Water Surface: Understanding the Vibrational Spectral Signatures of the Charge Defect in Cryogenically Cooled H+(H2O)n=2 – 28 Clusters

In this Article, we review the role of gas-phase, size-selected protonated water clusters, H+(H2O)n, in the analysis of the microscopic mechanics responsible for the behavior of the excess proton in bulk water. We extend upon previous studies of the smaller, two-dimensional sheet-like structures to larger (n≥10) assemblies with three-dimensional cage morphologies which better mimic the bulk environment. Indeed, clusters in which a complete second solvation shell forms around a surface-embedded hydronium ion yield vibrational spectra where the signatures of the proton defect display strikingly similar positions and breadth to those observed in dilute acids. We investigate effects of the local structure and intermolecular interactions on the large red shifts observed in the proton vibrational signature upon cluster growth using various theoretical methods. We show that, in addition to sizeable anharmonic couplings, the position of the excess proton vibration can be traced to large increases in the electric field exerted on the embedded hydronium ion upon formation of the first and second solvation shells. MAJ acknowledges support from the U.S. Department of Energy under Grant No. DE-FG02- 06ER15800 as well as the facilities and staff of the Yale University Faculty of Arts and Sciences High Performance Computing Center, and by the National Science Foundation under Grant No. CNS 08-21132 that partially funded acquisition of the facilities. SMK and SSX acknowledge support from the US Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences and Biosciences. Pacific Northwest National Laboratory (PNNL) is a multiprogram national laboratory operated for DOE by Battelle. This research used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231.

Fournier, Joseph A.↗
DOE OSTI2015

Ultrafast dynamics of liquid water: Energy relaxation and transfer processes of the OH stretch and the HOH bend

The vibrational energy relaxation and transfer processes of the OH stretching and the HOH bending vibrations in liquid water are investigated via the theoretical calculation of the pump-probe spectra obtained from non-equilibrium molecular dynamics simulations with the TTM3-F interaction potential. The excitation of the OH stretch induces an instantaneous response of the high frequency librational motions in the 600-1000 cm-1 range. In addition, the excess energy of the OH stretch of a water molecule quickly transfers to the OH stretches of molecules in its first hydration shell with a time constant of ~50 fs, followed by relaxation to the HOH bends of the surrounding molecules with a time constant of 230 fs. The excitation of the HOH bend also results in the ultrafast excitation of the high frequency librational motions. The energy of the excited HOH bend of a water molecule decays, with a time constant of 200 fs, mainly to the relaxation of the HOH bends of its surrounding molecules. The energies of the HOH bends were found to transfer quickly to the intermolecular motions via the coupling with the high frequency librational motions. The excess energy of the OH stretch or the HOH bend relaxes to the high frequency intermolecular librational motions and eventually to the hot ground state with a time scale of ~1 ps via the coupling with the librational and translational motions. The energy relaxation and transfer processes were found to depend on the local hydrogen bonding network; the relaxations of the excess energy of the OH stretch and the HOH bend of four- and five-coordinated molecules are faster than those of a three-coordinated molecule due to the delocalization of the vibrational motions of the former (four- and five-coordinated molecules) compared to those of the later (three-coordinated molecules). The present results highlight the importance of the high frequency intermolecular librational modes in facilitating the ultrafast energy relaxation process in liquid water via their strong nonlinear couplings with the intramolecular OH stretching and HOH bending vibrations. S.S.X. acknowledges the support of the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences. Pacific Northwest National Laboratory (PNNL) is a multiprogram national laboratory operated for DOE by Battelle. The calculation was carried out using the computing resources at the Research Center for Computational Science in Okazaki, Japan.

Imoto, Sho↗
DOE OSTI2016

Dependence of the rate of LiF ion pairing on the description of molecular interaction

We present an analysis of the dynamics of ion-pairing of Lithium Fluoride (LiF) in aqueous solvent using both detailed molecular simulation as well as reduced models within a Gener- alized Langevin Equation (GLE) framework. We explored the sensitivity of the ion-pairing phenomena to the details of descriptions of molecular interaction, comparing two empirical potentials to explicit quantum based density functional theory. We find quantitative differences in the potentials of mean force for ion-pairing as well as time dependent frictions that lead to variations in the rate constant and reactive flux correlation functions. These details reflect differences in solvent response to ion-pairing between different representations of molecular interaction and influence anharmonicity of the dynamic response. We find that the short time anharmonic response is recovered with a GLE parameterization. Recovery of the details of long time response may require extensions to the reduced model. We show that the utility of using a reduced model leads to a straight forward application of variational transition state the- ory concepts to the condensed phase system. The significance of this is reflected in the analysis of committor distributions and the variation of planar hypersurfaces, leading to an improved understanding of factors that determine the rate of LiF ion-pairing. CJM and GKS are supported by the U.S. Department of Energy‘s (DOE) Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences and Biosciences. Pacific Northwest Na- tional Laboratory (PNNL) is operated for the Department of Energy by Battelle. MDB is grateful for the support of Laboratory Directed Research and Development funding under the auspices of PNNL’s Laboratory Initiative Materials Synthesis and Simulation across Scales (MS3). Additional computing resources were generously allocated by PNNL’s Institutional Computing program. EP acknowledges support from PNNL’s Alternate Sponsored Fellowship program and IMPRS Dres- den.Support to P.J. from the Czech Science Foundation (grant P208/12/G016) and the Academy of Sciences (Praemium Academie award) is gratefully acknowledged.

Pluharova, Eva↗
DOE OSTI2016

Local Aqueous Solvation Structure Around Ca2+ During Ca2+---Cl– Pair Formation

The molecular details of single ion solvation around Ca2+ and ion-pairing of Ca2--Cl- are investigated using ab initio molecular dynamics. The use of empirical dispersion corrections to the BLYP functional are investigated by comparison to experimentally available extended X-ray absorption fine structure (EXAFS) measurements, which probes the first solvation shell in great detail. Besides finding differences in the free-energy for both ion-pairing and the coordination number of ion solvation between the quantum and classical descriptions of interaction, there were important differences found between dispersion corrected and uncorrected density functional theory (DFT). Specifically, we show significantly different free-energy landscapes for both coordination number of Ca2+ and its ion-pairing with Cl- depending on the DFT simulation protocol. Our findings produce a self-consistent treatment of short-range solvent response to the ion and the intermediate to long-range collective response of the electrostatics of the ion-ion interaction to produce a detailed picture of ion-pairing that is consistent with experiment. MDB is supported by MS3 (Materials Synthesis and Simulation Across Scales) Initiative at Pacific Northwest National Laboratory. It was conducted under the Laboratory Directed Research and Development Program at PNNL, a multiprogram national laboratory operated by Battelle for the U.S. Department of Energy. CJM acknowledges support from US Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences & Biosciences. This research used resources of the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. Additional computing resources were generously allocated by PNNL's Institutional Computing program. The authors thank Prof. Tom Beck for discussions regarding QCT, and Drs. Greg Schenter and Shawn Kathmann for insightful comments.

Baer, Marcel D.↗
DOE OSTI2016

Tribute to Bruce C. Garrett

This special issue of The Journal of Physical Chemistry B honors the many and wide-ranging contributions of Bruce C. Garrett to theoretical chemistry. Bruce’s contributions have spanned the scope of theoretical chemistry from undergraduate trajectory studies with Don Bunker at UC Irvine, to graduate work on semiclassical reactive dynamics at UC Berkeley with Professor William H. Miller, through postdoctoral studies mainly on chemical reaction theory at the University of Minnesota with one of us, and then on to a long independent career studying a wide array of complex problems, mainly in kinetics and liquid-phase chemistry at Battelle Columbus, at Chemical Dynamics Corporation, and at Pacific Northwest National Laboratory (PNNL). Let's start by recalling some of his research. This work was supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences. The calculations were carried out using computer resources provided by the Office of Basic Energy Sciences.

Dang, Liem X.↗
DOE OSTI2016

Divalent Ion Parameterization Strongly Affects Conformation and Interactions of an Anionic Biomimetic Polymer

The description of peptides and the use of molecular dynamics simulations to refine structures and investigate the dynamics on an atomistic scale are well developed. Through a consensus in this community over multiple decades, parameters were developed for molecular interactions that only require the sequence of amino-acids and an initial guess for the three-dimensional structure. The recent discovery of peptoids will require a retooling of the currently available interaction potentials in order to have the same level of confidence in the predicted structures and pathways as there is presently in the peptide counterparts. Here we present modeling of peptoids using a combination of ab initio molecular dynamics (AIMD) and atomistic resolution classical forcefield (FF) to span the relevant time and length scales. To properly account for the dominant forces that stabilize ordered structures of peptoids, namely steric-, electrostatic, and hydrophobic interactions mediated through sidechain-sidechain interactions in the FF model, those have to be first mapped out using high fidelity atomistic representations. A key feature here is not only to use gas phase quantum chemistry tools, but also account for solvation effects in the condensed phase through AIMD. One major challenge is to elucidate ion binding to charged or polar regions of the peptoid and its concomitant role in the creation of local order. Here, similar to proteins, a specific ion effect is observed suggesting that both the net charge and the precise chemical nature of the ion will need to be described. MDD was supported by MS3 (Materials Synthesis and Simulation Across Scales) Initiative at Pacific Northwest National Laboratory. Research was funded by the Laboratory Directed Research and Development program at Pacific Northwest National Laboratory. MDB acknowledges support from US Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Material & Engineering. CJM acknowledges support from US Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences & Biosciences. PNNL is a multiprogram national laboratory operated by Battelle for the U.S. Department of Energy.

Daily, Michael D.↗
DOE OSTI2019

The Diverse Nature of Ion Speciation at the Nanoscale Hydrophobic/Water Interface

Herein we investigate ion pairing at weakly charged hydrophobic/water interfaces. In particular, the speciation of the thiocyanate (SCN-) ion in the vicinity of the nano aqueous (NaSCN solution) - hydrophobic (dodecane) interface in the presence of tetraalkylammonium groups attached to dodecyl chains (DTA+). This combined simulation and experimental approach finds that the SCN- interfacial ion pairs already form at concentrations as low as 5 mM and exist in three distinct, detectable structural motifs controlled by the SCN- concentration in the aqueous phase. The first is the association with the DTA+, followed by hydrophobic adsorption to the interface with the sulfur end of the SCN-, followed by an oppositely aligned SCN- absorption with a solvated nitrogen end of the SCN-. These different structures are formed due to subtle shifts in interactions and their elucidation requires multiple surface sensitive measurements in conjunction with the description of the interface using a representation of the interactions with quantum density functional theory (DFT). Together, this study reveals that we can both simultaneously control and detect subtle changes in the speciation of ions in the vicinity of soft interfaces. This work is supported by the Julia Jacobi Foundation and the European Research Council (grant number 616305). Christopher J. Mundy is supported by the US Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences (BES) Division of Chemical Sciences, Geosciences, and Biosciences. Marcel D. Baer is supported within the Materials Synthesis and Simulation Across Scales (MS3) Initiative through the Laboratory Research and Development Program at Pacific Northwest National Laboratory (PNNL). The calculations utilized resources of the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility supported by the Office of Science of the US DOE under Contract No. DE-AC02-05CH11231. PNNL is a multiprogram national laboratory operated by Battelle for the US DOE.

Zdrali, Evangelia↗
DOE OSTI2020

Experimental and DFT Calculated IR Spectra of Guests in Zeolites: Acyclic Olefins and Host-Guest Interactions

We performed experimental and periodic density functional theory (DFT) IR spectroscopy to investigate the adsorption of acyclic olefins over both acidic and non-acidic zeolites. Two conjugated polyenes, 2,4-dimethyl-1,3-pentadiene (I) and 2,6-dimethyl-2,4,6-octatriene (II) were studied to probe organic intermediates that can form during methanol conversion, and that can lead to deactivating species known collectively as “coke.” We computed vibrational spectra using zeolite-adsorbed and gas-phase models for both neutral and protonated forms of I and II, and compared these DFT results to diffuse reflectance IR Fourier transform (DRIFT) spectra of zeolite-guest systems. Our experimental and computational results are precise enough to pinpoint that the gauche s-cis conformation of species I predominates during adsorption over de-aluminated zeolite beta. Computed zeolite-adsorbed spectra of the protonated species I and II best represent the DRIFT spectra obtained after the adsorption of the olefins on HMOR at 20 ?, with computed bands at 1543 and 1562 cm–1 for molecules I+ and II+, respectively, attributed to the allylic stretching mode, ?(C=C–C+). These computed band frequencies are within 6 cm–1 of experimental data and confirm that the interaction between neutral acyclic olefins and acidic zeolites leads to protonation of the olefin. Comparison of computed spectra of the protonated species in the gas phase to those in the zeolite indicates that the electrostatic interaction between alkenyl and alkadienyl cations and negative zeolite framework does not significantly impact the position of the allylic stretching bands. These results highlight that computed spectroscopy and thermodynamics coupled with experimental spectra can be used for elucidating complex mixtures in zeolites, and certain spectral features of adsorbed olefins can be accurately modeled by gas-phase calculations. ACKNOWLEDGMENTS BM and SMA were supported by the National Science Foundation under Award # CBET-1512442. Acknowledgment is made to the Donors of the American Chemical Society Petroleum Research Fund for support (or partial support) of this research (#589978-ND5, EDH and FCJ). MDB was supported by the US Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences (BES), Division of Materials Science and Engineering. CJM was supported by the DOE, Office of Science, BES, Division of Chemical Sciences, Geosciences, and Biosciences. Pacific Northwest National Laboratory (PNNL) is a multiprogram national laboratory operated by Battelle for the U.S. DOE. Computing resources were generously provided through PNNL’s institutional computing (PIC).

Manookian, Babgen↗
DOE OSTI2020

Binding of Formic Acid on Anatase TiO2(101)

The adsorption of formic acid (FA) and the formation of formate species on metal oxide surfaces are of great interest in catalysis. Formic acid is used to probe the adsorption properties of surface sites and formates are common intermediates in many catalytic reactions. Here we focus on the interaction of FA with a prototypical anatase TiO2(101) surface by using a combination of scanning tunneling microscopy (STM), infrared reflection absorption spectroscopy (IRAS), electron stimulated desorption (ESD), and density functional theory (DFT) to assess the coverage dependent evolution of different FA-derived surface species at low temperatures (80–240 K). We find that isolated FA adsorbs at 80 K in both monodentate (MD) and bidentate (BD) configurations on top of undercoordinated Ti sites (Ti5c). The MD form is likely deprotonated and readily converts upon annealing to bridging BD on two neighboring Ti5c sites. DFT calculations show that molecularly-bound MD FA is metastable and readily converts into a deprotonated state in the proximity of subsurface oxygen vacancy. The stability of the MD species increases as the availability of the paired Ti5c sites for BD formation diminishes at high coverages. Upon surface saturation, a mixed configuration of alternating bidentate and monodentate species with coverage of 2/3 FA/Ti5c represents the most favorable configuration. This is in contrast with the adsorption of FA on rutile TiO2(110) and many other oxides, where bidentate formate species dominate. Y.W., A.D., G.A.K., R.R., N.G.P., and Z.D. were supported by the U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences (BES), Division of Chemical Sciences, Geosciences and Biosciences (CSGB) under FWP 47319. B.W. and A.S. acknowledge the support of DOE BES, CSGB Division under Award DESC0007347. The experimental studies were performed in EMSL, a national scientific user facility sponsored by the Department of Energy’s Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory (PNNL). PNNL is a multiprogram national laboratory operated for DOE by Battelle. The authors also acknowledge computational resources from the TIGRESS high-performance computer center at Princeton University.

Wang, Yang↗
DOE OSTI2021

Formation of Gas-Phase Allyl Radicals from Glycerol on Rutile TiO2(110)

The reaction pathways of glycerol on partially reduced TiO2(110) have been studied by a combination of molecular beam dosing and temperature-programmed desorption (TPD). Glycerol is found to bind more strongly to TiO2(110) than glycols due to an additional hydroxyl group in the molecule. When the Ti surface sites are saturated by glycerol, the majority of glycerol molecules (~ 95%) undergo further reactions to yield carbon-free (water and hydrogen) and carbon-containing products that are desorbed during the TPD ramp. In contrast to glycols for which both water and hydrogen arise exclusively from the cleavage of glycol O-H bonds, we identify two distinct reaction channels resulting from the cleavage of O-H and C-H bonds, respectively, for both water and hydrogen evolution from glycerol on TiO2(110). Quantitative determination of the desorption yields of water and hydrogen reveals that the C-H bond scission channel which occurs at higher temperature than the O-H bond scission channel, is a minor reaction channel for water formation, while it is a major reaction channel responsible for hydrogen evolution. For the carbon-containing products, we exclude both carbon deposition and carbon oxides (CO and CO2) formation, and identify the allyl radical as the major carbon-containing product. The formation of the allyl radical is accompanied by the formation of various minor carbon-containing products including propylene, acrolein, allyl alcohol, acetone, formaldehyde and ethylene. The different surface chemistries observed between glycerol and simpler alcohols and glycols on TiO2(110), suggest that smaller oxygenates cannot serve as models for larger oxygenates in probing their reaction pathways on oxide surfaces. The new reaction channels observed for glycerol on TiO2(110) may open up a new way to convert bio-renewable glycerol directly to pure hydrogen and hydrocarbons on TiO2. This work was supported by the US Department of Energy, Office of Sciences, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences & Biosciences, and performed in EMSL, a national scientific user facility sponsored by the Department of Energy’s Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory (PNNL). PNNL is a multiprogram national laboratory operated for the DOE by Battelle.

Chen, Long↗
DOE OSTI2015

Nanocrystalline Anatase Titania Supported Vanadia Catalysts: Facet-dependent Structure of Vanadia

Titania supported vanadia, a classic heterogeneous catalyst for redox reactions, typically has nonhomogeneous vanadia species on various titania facets, making it challenging not only to determine and quantify each species but also to decouple their catalytic contributions. We prepared truncated tetragonal bipyramidal (TiO2-TTB) and rod-like (TiO2-Rod) anatase titania with only {101} and {001} facets at ratios of about 80:20 and 93:7, respectively, and used them as supports of sub-monolayer vanadia. The structure and redox properties of supported vanadia were determined by XRD, TEM, XPS, EPR, Raman, FTIR and TPR, etc. It was found that vanadia preferentially occupy TiO2 {001} facets and form isolated O=V4+(O-Ti)2 species, and with further increase in vanadia surface coverage, isolated O=V5+(O-Ti)3 and oligomerized O=V5+(O-M)3 (M = Ti or V) species form on TiO2 {101} facets. The discovery on support facet-dependent structure of vanadia on anatase titania is expected to enable the elucidation of structure-function correlations on high surface area TiO2 supported vanadia catalysts. This work was supported by U. S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Biosciences and Geosciences. The research was performed in the Environmental Molecular Sciences Laboratory (EMSL), a national scientific user facility sponsored by the DOE Office of Biological and Environmental Research, and located at Pacific Northwest National Laboratory (PNNL). PNNL is operated for DOE by Battelle.

catalyst characterization↗
DOE OSTI2015

Photochemical Outcomes of Adsorbed Oxygen: Desorption, Dissociation and Passivation by Coadsorbed Water

A mixed Fe+Cr oxide surface was used to explore the photochemical fate of adsorbed O2 under ultrahigh vacuum (UHV) conditions. This mixed oxide surface possessed a magnetite-like (111) structure based on low energy electron diffraction (LEED), with its chemical behavior resembling that of Fe3O4(111). Oxygen adsorption at 40 K resulted in two chemisorption states, a strongly bound form that desorbed in temperature programmed desorption (TPD) at 230 K and a weakly bound form that evolved at 100 K. The former was assigned to charge transfer adsorption at Fe2+ sites and the latter to electrostatic binding at Fe3+ sites. A minority state of O2 was also detected in TPD at ~160 K, and tentatively assigned to adsorption at Cr3+ sites. The 230 K O2 state was the focus of photochemical studies employing four wavelengths of light from the red to the UV. Irrespective of wavelength, O2 molecules in the 230 K state preferentially photodesorbed when irradiated, with comparable rates across the visible and an order of magnitude greater in the UV. Also, approximately 10% of adsorbed O2 irreversibly photodissociated, irrespective of wavelength, with the resulting fragments blocking access to cation sites for subsequent O2 adsorption. Preadsorbed water also blocked O2 adsorption, but post-adsorbed water stabilized O2 at Fe2+ sites in TPD to 285 K. However, water-stabilized O2 molecules were insensitive to photodesorption. O2 photodissociation and photo-passivation both represent potentially adverse outcomes in the release of O2 during the heterogeneous water photooxidation reaction. The author thanks Drs. Sara Chamberlin and Scott Chambers for supplying the film used in this work. This work was supported by the US Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences & Biosciences. Pacific Northwest National Laboratory (PNNL) is a multiprogram national laboratory operated for DOE by Battelle. The research reported here was performed using EMSL, a national scientific user facility sponsored by the Department of Energy's Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory.

mixed oxide, photodesorption, oxygen, temperature ↗
DOE OSTI2016

Reaction Rate Theory in Coordination Number Space: An Application to Ion Solvation

Understanding reaction mechanisms in many chemical and biological processes require application of rare event theories. In these theories, an effective choice of a reaction coordinate to describe a reaction pathway is essential. To this end, we study ion solvation in water using molecular dynamics simulations and explore the utility of coordination number (n = number of water molecules in the first solvation shell) as the reaction coordinate. Here we compute the potential of mean force (W(n)) using umbrella sampling, predicting multiple metastable n-states for both cations and anions. We find with increasing ionic size, these states become more stable and structured for cations when compared to anions. We have extended transition state theory (TST) to calculate transition rates between n-states. TST overestimates the rate constant due to solvent-induced barrier recrossings that are not accounted for. We correct the TST rates by calculating transmission coefficients using the reactive flux method. This approach enables a new way of understanding rare events involving coordination complexes. We gratefully acknowledge Liem Dang and Panos Stinis for useful discussion. This research used resources of the National Energy Research Scientific Computing Center, a DOE Office of Science User Facility supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. SR, CJM, and GKS were supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences. MDB was supported by MS3 (Materials Synthesis and Simulation Across Scales) Initiative, a Laboratory Directed Research and Development Program at Pacific Northwest National Laboratory (PNNL). PNNL is a multiprogram national laboratory operated by Battelle for the U.S. Department of Energy.

Roy, Santanu↗
DOE OSTI2016

Anharmonicity and confinement in zeolites: Structure, spectroscopy, and adsorption free energy of ethanol in H-ZSM-5

Here, to account for thermal and entropic effects caused by the dynamics of the motion of the reaction intermediates, ethanol adsorption on the Brønsted acid site of the H-ZSM-5 catalyst has been studied at different temperatures and ethanol loadings using ab initio molecular dynamics (AIMD) simulations, infrared (IR) spectroscopy and calorimetric measurements. At low temperatures (T ≤ 400 K) and ethanol loading, a single ethanol molecule adsorbed in H-ZSM-5 forms a Zundel-like structure where the proton is equally shared between the oxygen of the zeolite and the oxygen of the alcohol. At higher ethanol loading, a second ethanol molecule helps to stabilize the protonated ethanol at all temperatures by acting as a solvating agent. The vibrational density of states (VDOS), as calculated from the AIMD simulations, are in excellent agreement with measured IR spectra for C 2 H 5 OH, C 2 H 5 OD and C 2 D 5 OH isotopomers and support the existence of both monomers and dimers. A quasi-harmonic approximation (QHA), applied to the VDOS obtained from the AIMD simulations, provides estimates of adsorption free energy within ~10 kJ/mol of the experimentally determined quantities, whereas the traditional approach, employing harmonic frequencies from a single ground state minimum, strongly overestimates the adsorption free energy by at least ~30 kJ/mol. This discrepancy is traced back to the inability of the harmonic approximation to represent the contributions to the vibrational motions of the ethanol molecule upon confinement in the zeolite. KA, MFR, GBM were supported by the Long Term Structural Methusalem Funding by the Flemish Government – grant number BOF09/01M00409. MSL, VAG, RR and JAL were supported by the US Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences & Biosciences. PNNL is a multiprogram national laboratory operated for DOE by Battelle. Computational resources were provided at W. R. Wiley Environmental Molecular Science Laboratory (EMSL), a national scientific user facility sponsored by the Department of Energy’s Office of Biological and Environmental Research located at PNNL, the National Energy Research Scientific Computing Center (NERSC) at Lawrence Berkeley National Laboratory and the Stevin Supercomputer Infrastructure at Ghent University.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗