Hydrolysis of zeolite framework aluminum and its impact on acid catalyzed alkane reactions
Not Available
Xue, Nianhua↗
SEARCH · Engineering Papers
Results for “MULTIPROGRAMMING”
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.
Counting surface redox sites in carbon-supported electrocatalysts by cathodic stripping of O deposited from N2O
Not Available
Egbert, Jonathan D.↗
Thermodynamic and Kinetic Roles of H2 in Structure Evolution of Urchin-like Co: A density functional theory study
Small gas molecules working as good capping agents play important roles in controlling the morphologies and surface structures of the metal nanocrystals. In present work, the thermodynamic and kinetic roles of H2 molecule played on the morphology of Co nanocrystal were systematically studied based on density functional theory (DFT). The Gibbs surface free energies (?) of Co(100), Co(110) and Co(111) under different hydrogen surface coverages were investigated by ab initio thermodynamics, based on which i) the phase diagrams of stable H coverage on each plane were obtained; ii) the morphology evolutions of the Co nanocrystals with various surface hydrogen coverages were further predicted by Wulff construction. It was revealed that addition of H2 could change the facet stability generating diverse morphological Co nucleus. The kinetic role of H2 was further explored by DFT during adatom Co surface diffusions under different H coverages (?H), which suggested that except for Co(100) at ?H of 0.56 ML, the surface H would hinder Co surface diffusions. The projected density of state (PDOS) gave a deeper insight that the electronic structures of the Co adatom could be alerted by addition of the surface H, which thereby affected its surface diffusion ability. X. Wang, N. Liu, Q. Zhang, X. Liang, and B. Chen were financially supported by the National Natural Science Foundation of China (No. 21476012, 21571012, and 91534201). D. Mei was supported by the US Department of Energy (DOE), 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. Computing time was granted by the facility user proposal grand challenge of computational catalysis of the William R. Wiley Environmental Molecular Sciences Laboratory (EMSL). EMSL is a national scientific user facility located at Pacific Northwest National Laboratory (PNNL) and sponsored by DOE’s Office of Biological and Environmental Research.
Nanocrystal Co, Surface free energy, Surface diffu↗
Quenching of electron transfer reactions through coadsorption: A study of oxygen photodesorption from TiO2(110)
Using temperature programmed desorption (TPD) and photon stimulated desorption (PSD), we show that coadsorbates of varying binding energies on the rutile TiO2(110) surface exert a commensurate inhibiting influence on the hole-mediated photodesorption of adsorbed O2. A variety of coadsorbates (Ar, Kr, Xe, N2, CO, CO2, CH4, N2O, acetone, methanol or water) were shown to quench O2 photoactivity, with the extent correlating with the coadsorbate’s gas phase basicity, which in turn determines the strength of the coadsorbate-Ti4+ bond. Coadsorbed rare gases inhibited the photodesorption of O2 by ~10-25%, whereas strongly bound species (water, methanol and acetone) nearly completely inhibited O2 PSD. We suggest that coadsorption of these molecules inhibit the arrival probability of holes to the surface. Band bending effects, which vary with the extent of charge transfer between the coadsorbate and the TiO2(110) surface, are not expected to be significant in the cases of the rare gases and physisorbed species. These results indicate that neutral coadsorbates can exert a significant influence on charge transfer events by altering the interfacial dipole in the vicinity of the target molecule. This work was supported by the US Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences & Biosciences. The work 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 (PNNL). PNNL is a multiprogram national laboratory operated for DOE by Battelle under Contract DE-AC05-76RL01830.
Environmental Molecular Sciences Laboratory↗
A New, Dispersion-Driven Intermolecular Arrangement for the Benzene–Water Octamer Complex: Isomers and Analysis of their Vibrational Spectra
The anharmonic spectra of the gas phase cubic water octamer (W8, D2d and S4 isomers) have been theoretically calculated at the second order Møller-Plesset perturbation (MP2) and the Coupled Cluster with Single, Double and a perturbative estimate of Triple replacements [CCSD(T)] theories. The CCSD(T) harmonic frequencies are the first ones reported for this cluster. An additional band at ~3500 cm-1 is reported, in a spectral area that was not previously accessible experimentally due to technical reasons. The IR spectra of the S4 isomer have a larger number of fundamental bands than the spectra of the D2d isomer, in accordance with the presence of lower symmetry in the former compared to the latter. When W8 interacts with benzene (BZ), the dispersion interaction plays a major role in determining the resulting intermolecular arrangement. Calculations at the MP2, DFT (with the B97XD functional which includes dispersion corrections) and CCSD(T) levels of theory suggest an optimal arrangement in which BZ is almost parallel to one of the faces of the W8 cube. This is in sharp contrast with the previously reported structure in which one of the “free” OH bonds of the W8 cube pointed towards the center of BZ (Science 276, 1678 (1997)) that was determined at the DFT level with the dispersionless B3LYP functional and was used to assign the experimentally measured IR spectra. Five low-lying isomers, three of the S4-like and two of the D2d-like type, were determined and their spectra were assigned. The perturbation of BZ to the W8 bands amounts to (i) the localization of the normal modes of W8 thus resulting in more IR active bands and (ii) the lowering of the overall symmetry of the complex that results in the splitting of the doubly degenerate bands on the bare W8. Our results further suggest that a future recording of the IR spectra in the HOH bending region can definitively aid in the assignment of the various isomers of both the BZ and the BZ-W8 complexes. EM 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 the US 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.
Miliordos, Evangelos↗
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.↗
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↗
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.↗
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.↗
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↗
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↗
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↗
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↗
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 ↗
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↗
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↗
Observation of molecular hydrogen produced from bridging hydroxyls on anatase TiO2(101)
Anatase TiO2 is used extensively in a wide range of catalytic and photocatalytic processes and is a promising catalyst for hydrogen production. Here, we show that molecular hydrogen was produced from bridging hydroxyls, HOb, on the (101) surface of single crystal anatase (a-TiO2(101)). Electron bombardment at 30 K produced bridging oxygen vacancies in the surface. Deuterated bridging hydroxyls, DOb, were subsequently formed via dissociation of adsorbed D2O and confirmed by infrared reflection-absorption spectroscopy. During temperature programmed desorption (TPD) spectroscopy, D2 desorption was observed at 520 K. Density functional theory calculations show that both H2 and H2O production from HOb are endothermic at 0 K on a-TiO2(101), but H2 (H2O) desorption is entropically driven above 230 K (800 K). The calculated activation barrier for H2 desorption is 1.40 eV, which is similar to the desorption energy obtained from analysis of the D2 TPD spectra. The H2 desorption likely proceeds in two steps: H atom diffusion on the surface and then recombination. The calculations indicate that diffusion of the hydrogen atoms into the subsurface could also be an important process. We thank Drs. Z. Dohnalek and R. Rousseau for stimulating discussions. This work was supported by the US Department of Energy, Office of Science, 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 DOE by Battelle. We also acknowledge the WPI computational facilities for providing resources to perform the DFT studies.
Deskins, N A.↗
The melting temperature of liquid water with the effective fragment potential
Direct simulation of the solid-liquid water interface with the effective fragment potential (EFP) via the constant enthalpy and pressure (NPH) ensemble was used to estimate the melting temperature (Tm) of ice-Ih. Initial configurations and velocities, taken from equilibrated constant pressure and temperature (NPT) simulations at T = 300 K, 350 K and 400 K, respectively, yielded corresponding Tm values of 378±16 K, 382±14 K and 384±15 K. These estimates are consistently higher than experiment, albeit to the same degree with previously reported estimates using density functional theory (DFT)-based Born-Oppenheimer simulations with the Becke-Lee-Yang-Parr functional plus dispersion corrections (BLYP-D). KRB was supported by a Computational Science Graduate Fellowship from the Department of Energy. MSG was supported by a U.S. National Science Foundation Software Infrastructure (SI2) grant (ACI – 1047772). SSX acknowledges 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.
Brorsen, Kurt R.↗