SEARCH · Engineering Papers
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.
In memorium: This special issue is dedicated to Steve Attaway, who passed away on February 28, 2019. Steve Attaway worked at Sandia National Laboratories in Albuquerque, NM, for over 30 years making significant contributions in highperformance computing, shock physics, meshfree methods, the geosciences, concrete mechanics, and blast effects on structures.
The reactivity of Ru/Al2O3 catalysts in the hydrogenation of CO/CO2 gas mixtures is investigated in this work. Accordingly kinetic measurements are carried out at process conditions of industrial interest (i.e., at high COx per-pass conversion and with concentrated COx/H2 streams), and the surface species on the catalyst are investigated through transient DRIFTS-MS analysis. The catalyst shows gradual deactivation when the methanation is carried out in the presence of CO in the gas feed at low temperatures (200-300 °C). However, stable performance is observed at higher temperatures, showing CH4 yields even higher than those observed during methanation of a pure CO2 feed. The catalyst deactivation at low temperature is reversible as thermal treatments at temperatures above 350 °C are able to restore the catalytic activity. DRIFTS-MS experiments revealed that CO2 methanation involves a reaction pathway where CO2 is adsorbed as bicarbonate on Al2O3 and successively hydrogenated to methane on Ru, passing through formate and carbonyl intermediates. In the presence of CO at low temperatures, the catalyst shows a higher CO coverage of the Ru sites, a larger formate coverage of alumina sites and the presence of adsorbed carbonaceous species, identified as hydrocarbon and carboxylate species. Carboxylates remain on the surface, effectively blocking adsorption sites thus leading to catalyst deactivation. However, working above the carboxylate decomposition temperature ensures a clean catalyst surface, resulting in stable and high performance in CO/CO2 methanation. Catalyst preparation, characterization and testing activity were carried out with a joint project between the Politecnico di Milano and the Italian National Agency for New Technologies, Energy and Sustainable Economic Development (ENEA), which is gratefully acknowledged. The authors gratefully acknowledge the financial support of the DRIFTS experiments by the US Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, Chemical Sciences, Geosciences, and Biosciences Division. L.F. also acknowledges the Pacific Northwest National Laboratory (PNNL) for the financial support during his permanence as visiting student.
Forecasting high-dimensional dynamical systems is a fundamental challenge in various fields, such as geosciences and engineering. Neural Ordinary Differential Equations (NODEs), which combine the power of neural networks and numerical solvers, have emerged as a promising algorithm for forecasting complex nonlinear dynamical systems. However, classical techniques used for NODE training are ineffective for learning chaotic dynamical systems. In this work, we propose a novel NODE-training approach that allows for robust learning of chaotic dynamical systems. Here, our method addresses the challenges of non-convexity and exploding gradients associated with underlying chaotic dynamics. Training data trajectories from such systems are split into multiple, non-overlapping time windows. In addition to the deviation from the training data, the optimization loss term further penalizes the discontinuities of the predicted trajectory between the time windows. The window size is selected based on the fastest Lyapunov time scale of the system. Multi-step penalty(MP) method is first demonstrated on Lorenz equation, to illustrate how it improves the loss landscape and thereby accelerates the optimization convergence. MP method can optimize chaotic systems in a manner similar to least-squares shadowing with significantly lower computational costs. Our proposed algorithm, denoted the Multistep Penalty NODE, is applied to chaotic systems such as the Kuramoto-Sivashinsky equation, the two-dimensional Kolmogorov flow, and ERA5 reanalysis data for the atmosphere. It is observed that MP-NODE provide viable performance for such chaotic systems, not only for short-term trajectory predictions but also for invariant statistics that are hallmarks of the chaotic nature of these dynamics.
Herein we present a new modeling paradigm for optimization that we call random field optimization. Random fields are a powerful modeling abstraction that aims to capture the behavior of random variables that live on infinite-dimensional spaces (e.g., space and time) such as stochastic processes (e.g., time series, Gaussian processes, and Markov processes), random matrices, and random spatial fields. This paradigm involves sophisticated mathematical objects (e.g., stochastic differential equations and space-time kernel functions) and has been widely used in neuroscience, geoscience, physics, civil engineering, and computer graphics. Despite of this, however, random fields have seen limited use in optimization; specifically, existing optimization paradigms that involve uncertainty (e.g., stochastic programming and robust optimization) mostly focus on the use of finite random variables. This trend is rapidly changing with the advent of statistical optimization (e.g., Bayesian optimization) and multi-scale optimization (e.g., integration of molecular sciences and process engineering). Our work extends a recently-proposed abstraction for infinite-dimensional optimization problems by capturing more general uncertainty representations. Moreover, we discuss solution paradigms for this new class of problems based on finite transformations and sampling, and identify open questions and challenges.
Minerals are solid state nuclear track detectors — nuclear recoils in a mineral leave latent damage to the crystal structure. Depending on the mineral and its temperature, the damage features are retained in the material from minutes (in low-melting point materials such as salts at a few hundred °C) to timescales much larger than the 4.5 Gyr-age of the Solar System (in refractory materials at room temperature). The damage features from the O (50) MeV fission fragments left by spontaneous fission of 238U and other heavy unstable isotopes have long been used for fission track dating of geological samples. Laboratory studies have demonstrated the readout of defects caused by nuclear recoils with energies as small as O (1) keV. This whitepaper discusses a wide range of possible applications of minerals as detectors for E R ≳ O (1) keV nuclear recoils: Using natural minerals, one could use the damage features accumulated over O (10) Myr–O (1) Gyr to measure astrophysical neutrino fluxes (from the Sun, supernovae, or cosmic rays interacting with the atmosphere) as well as search for Dark Matter. Using signals accumulated over months to few-years timescales in laboratory-manufactured minerals, one could measure reactor neutrinos or use them as Dark Matter detectors, potentially with directional sensitivity. Research groups in Europe, Asia, and America have started developing microscopy techniques to read out the O (1)–O (100) damage features in crystals left by O (0.1)–O (100) keV nuclear recoils. We report on the status and plans of these programs. Furthermore, the research program towards the realization of such detectors is highly interdisciplinary, combining geoscience, material science, applied and fundamental physics with techniques from quantum information and Artificial Intelligence.
In this paper, the IUPAC-IUGS joint Task Group “Isotopes in Geosciences” recommends a value of (106.25 ± 0.38) Ga for the half-life of 147 Sm, and a corresponding decay constant λ 147 = (6.524 ± 0.024) × 10 –12 a –1 , both with a coverage factor k = 2. For the extinct radionuclide 146 Sm two very different half-lives are used in the scientific community ( c . 68 and 103 Ma), to such a degree that no consensus value can be endorsed at present by the Task Group. Pending dedicated re-investigations it is recommended that papers using the 146 Sm decay to quantify the cosmo/geological evolution of (extra)terrestrial samples perform a twin set of calculations using both proposed half-lives.
Coupled thermo-hydro-mechanical (THM) processes in fractured rock are playing a crucial role in geoscience and geoengineering applications. Diverse and conceptually distinct approaches have emerged over the past decades in both continuum and discontinuum perspectives leading to significant progress in their comprehending and modeling. This review paper offers an integrated perspective on existing modeling methodologies providing guidance for model selection based on the initial and boundary conditions. By comparing various models, one can better assess the uncertainties in predictions, particularly those related to the conceptual models. The review explores how these methodologies have significantly enhanced the fundamental understanding of how fractures respond to fluid injection and production, and improved predictive capabilities pertaining to coupled processes within fractured systems. It emphasizes the importance of utilizing advanced computational technologies and thoroughly considering fundamental theories and principles established through past experimental evidence and practical experience. The selection and calibration of model parameters should be based on typical ranges and applied to the specific conditions of applications. The challenges arising from inherent heterogeneity and uncertainties, nonlinear THM coupled processes, scale dependence, and computational limitations in representing field scale fractures are discussed. Realizing potential advances on computational capacity calls for methodical conceptualization, mathematical modeling, selection of numerical solution strategies, implementation, and calibration to foster simulation outcomes that intricately reflect the nuanced complexities of geological phenomena. Future research efforts should focus on innovative approaches to tackle the hurdles and advance the state-of-the-art in this critical field of 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.
Erodants and abradants are two types of wear-causing particles (WCP). Despite their effects on the severity and rate of wear, WCP are inadequately characterized in a surprising number of publications, especially those involving impingement erosion and three-body abrasion. That shortfall makes it difficult to correlate features of WCP with the details of worn surfaces, or to develop wear models that account for those features. It is argued that the documentation of WCP should go beyond simply reporting their composition and mean particle size. In late 1985, ASTM Committee G2 on Wear and Erosion established a task group on the characterization of WCP. At that time, image analysis was slower and less sophisticated than it is today. While that ASTM task group failed to produce a consensus standard, computerized particle characterization methods were developing in fields other than tribology, fields like geoscience, heavy sand mining, materials processing, and pharmaceuticals. A notable exception to this is ferrography, which is a widely-used diagnostic for lubricated tribosystems. In the context of dry wear, it is useful to identify which features of WCP would be beneficial to document, and to identify some techniques and scales of detail appropriate to particular tribosystems. In this paper, examples are presented for the morphological and mechanical characterization of silica sand grains, as prompted by a triboanalysis of the abrasive and erosive wear of biomass pre-processing equipment. The authors propose a minimum level of documentation for WCP, one that can enrich tribosystem analysis both in laboratory tests and field studies.
Correction was made to the funding and acknowledgment: Funding V.-A.G. was supported by the U.S. Department of Energy, Office of Science, Basic Energy Sciences, Chemical Sciences, Geosciences, and Biosciences Division. R.R. and M.-T.N. were supported by the Laboratory Directed Research and Development Program at Pacific Northwest National Laboratory (PNNL) CheMSR Agile Investment. PNNL is operated by DE-AC05-76RLO01830. J.-B.L. and J.L. were supported by the National Natural Science Foundation of China (Nos. 21433005, 91645203, and 21590792). D.C.C. was supported by Research and Innovation at the University of Nevada, Reno. ACKNOWLEDGMENTS The research was performed using resources available through Research Computing at Pacific Northwest National Laboratory (PNNL), Tsinghua National Laboratory for Information Science and Technology and the Computational Chemistry Laboratory under Tsinghua Xuetang Talents Program, and University of Nevada, Reno High Performance Computing. J.- B.L. and J.L. acknowledge discussions with Dr. Yang-Gang Wang (Southern University of Science and Technology).
This errata is to correct the acknowledgment in the original published article. G.B. (partial support), V.-A.G., M.-T.N., and R.R. were supported by the U.S. Department of Energy, Office of Science, Basic Energy Sciences, Chemical Sciences, Geosciences, and Biosciences Division. Computer resources were provided by Research Computing at Pacific Northwest National Laboratory and the National Energy Research Scientific Computing Center (NERSC), a U.S. Department of Energy Office of Science User Facility operated under Contract No. DE-AC02-05CH11231. PNNL is operated by Battelle for the U.S. Department of Energy under Contract DE-AC05-76RL01830. G.B. acknowledges the EPSRC for partial support in PhD funding. The authors acknowledge the use of the IRIDIS High Performance Computing Facility (IRIDIS 5) and associated support services at the University of Southampton in the completion of this work. We are grateful to the UK Materials and Molecular Modeling Hub (Thomas HPC) for computational resources, which is partially funded by EPSRC (EP/P020194/1).
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.
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.
Research on multiply charged anions (MCAs) in the gas phase has been intensively performed during the last decades, mainly to understand fundamental molecular physics phenomena, e.g., intramolecular Coulomb repulsion and existence of the repulsive Coulomb barrier. However, the relevance of these investigations with respect to understanding MCAs’ chemistry appears often vague. Here, we discuss how insights on the electronic structure obtained from negative ion photoelectron spectroscopy (NIPES) combined with theoretical calculations and collision induced dissociation can provide fundamental understanding on the intrinsic chemical reactivity of MCAs and their fragments. This is exemplified on our studies on polyhedral closo-borate dianions [BnXn]2- (n = 6, 10, 11, 12; X = H, F-I, CN) and their fragment ions. For example, the rational design of closo-borate dianions with specific electronic properties is described, which leads to generating highly reactive fragments. Depending on the dianionic precursor, these fragments are either tuned to bind noble gases effectively or to activate small molecules like CO and N2. The intrinsic electronic properties of closo-borate dianions are further compared to their electrochemistry in solutions, revealing solvent effects on the redox potentials. Neutral host molecules such as cyclodextrins are found to bind strongly to [BnXn]2-, and gas phase NIPES provides insights on the intrinsic host-guest interactions. Finally, outlooks including the direct NIPES of molecular fragment ions which cannot be generated in the condensed phase and their utilization in preparative mass spectrometry are discussed. J.W is grateful to the Volkswagen foundation for a Freigeist Fellowship. X. B. W acknowledges support from the U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, the Division of Chemical Sciences, Geosciences and Biosciences. We are grateful to our colleagues who support our work on gas phase ion chemistry for many years. In particular we acknowledge the members of our groups, the Asmis group (Leipzig), Laskin Group (West Lafayette), the Jenne group (Wuppertal), the Finze group (Würzburg), the Grabowsky group (Bern), and the Sun group (Shanghai). We are grateful to Edoardo Apra for the many theoretical investigations on closo-borate anions. JW acknowledges support of Harald Knorke and Markus Rohdenburg with the preparation of graphical material, and XBW thanks Qinqin Yuan for the help in organizing references.
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).
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.
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.
As a general-purpose force field for molecular simulations of layered materials and their fluid interfaces, Clayff continues to see broad usage in atomistic computational modeling for numerous geoscience and materials science applications due to its (1) success in predicting properties of bulk nanoporous materials and their interfaces, (2) transferability to a range of layered and nanoporous materials, and (3) simple functional form which facilitates incorporation into a variety of simulation codes. Here, we review applications of Clayff to model bulk phases and interfaces not included in the original parameter set and recent modifications for modeling surface terminations such as hydroxylated nanoparticle edges. Overall, we conclude with a discussion of expectations for future developments.