Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles”

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 1,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

First-Principles Study of Interaction between Molecules and Lewis Acid Zeolites Manipulated by Injection of Energized Charge Carriers

Metal dopants in Lewis acid zeolites can be applied to manipulate the interaction between the active sites and molecular reactants. However, only a few elements can be introduced into the zeolite framework to generate strong Lewis acid sites. In this study, we propose an alternative approach to tune the interaction between zeolite and molecules by introducing energized charge carriers that are photoexcited in a complex of plasmonic metal and zeolites. Here, through density functional theory (DFT) calculations, we show that although the energized electrons do not impact the adsorption of ammonia, which is used as a probe molecule, on Zr- and Si-BEA, ammonia’s heat of adsorption can be reduced up to 60% on Sn-BEA and 20% on Ti-BEA, respectively, because of the localized population of the electrons at the Sn and Ti center. This most pronounced effect in Sn-BEA is attributed to the Sn-induced midgap states and the strong electronegativity of Sn atoms. In contrast, the energized holes can strengthen the interaction between zeolites and ammonia by at least 10%. Particularly, if an energized hole is present in the system, the heat of adsorption of ammonia within Si-BEA can increase 3-fold, shifting from weak to moderate interaction. We find that, similar to closed sites, the effect of energized charge carriers can also be extended to open sites of Lewis acid zeolites. Our study provides fundamental insights for unlocking the potentials of integrating plasmonic metals with zeolites and brings an additional degree of control to tune the interfacial interaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

First-Principles Investigation of Adsorbate–Adsorbate Interactions on Ni(111), Ni(211), and Ni(100) Surfaces

Periodic density functional theory calculations were employed to investigate the coadsorption patterns between CO and key water gas shift reaction (WGSR) intermediates, i.e., H2O, H, OH, O, and COOH, on Ni(111), Ni(100), and Ni(211) single-crystal surfaces. It has been shown that although the nature of these adsorbate pair interactions are predominantly repulsive, as commonly assumed in the literature, the interaction pattern for each adsorbate pair can be complicated by facet lattice structures and potential intermolecular hydrogen bonds, which play an unsubtle role in influencing the WGSR redox and carboxyl pathways. This investigation can enrich our descriptions of the adsorbate lateral interactions for surface chemistry modeling and will also provide theoretical insights into tuning the properties of nanoscale nickel catalysts by manipulating the surface structure and morphology for hydrogen fuel production.

water-gas shift reaction↗
DOE OSTI2019

A first-principles investigation of gas phase ring opening reaction of furan over HZSM-5 and Ga-substituted ZSM-5

The furan ring-opening reaction using three catalytic models, HZSM-5, the extra-framework sites [GaO]/ZSM-5, and [Ga(OH)(2)]/ZSM-5, was investigated using periodic density functional theory as a model reaction for mechanistic understanding of catalytic vapor phase biomass upgrading. The formation of formyl allene from the furan was investigated in detail by computing the energy profiles and reaction barriers. Based on the computed free energy profiles, the HZSM-5 and [Ga(OH)(2)]/ZSM-5 are active catalysts, whereas [GaO]/ZSM-5 is not active for the furan ring-opening reaction. In HZSM-5, the likely rate-controlling step is furan CH2-O bond cleavage upon protonation (1.70 eV). The computed energy difference between the highest and lowest points in the free energy profile is 1.89 eV. The [Ga(OH)(2)](+) is likely the catalytic center for the gallium-incorporated ZSM-5 catalyst. The rate-controlling step is a hydrogen transfer reaction, which requires a reaction barrier of 2.48 eV. The computed energy difference between the highest and lowest points in the Ga(OH)(2) catalyzed free energy profile is 2.61 eV. These computational studies provide new mechanistic understanding of the role of catalytic extra-framework sites in the vapor phase upgrading reactions.

Zhou, Mingxia↗
DOE OSTI2026

First-Principles Insights into Proton-Coupled Electron Transfer versus Hydrogen Evolution Reaction Selectivity from a Base-Appended Cobaltocene Mediator

Performing selective proton-coupled electron transfer (PCET) to substrates such as N 2 , CO 2 , and unsaturated organic molecules under electrochemical conditions requires the suppression of the competing hydrogen evolution reaction (HER). To address this challenge, our laboratory previously demonstrated a PCET mediator strategy using a dimethylaniline-appended cobaltocene complex, [(CpCoCp NMe2 )H] + , which performs selective reductive chemistry while suppressing the HER. However, the origin of the suppressed, yet still observable, HER has not been thoroughly established. In this work, we perform density functional theory (DFT) calculations to elucidate the HER mechanism involving this redox mediator and to provide atomistic insights into the bifurcation between the PCET and HER pathways. We find that protonation of the aniline moiety to form [CpCoCp NMe2H ] + is more favorable, both kinetically and thermodynamically, than formation of the ring-protonated species [(CpCo(Cp-H) NMe2 )] + . Furthermore, PCET to acetophenone is energetically more favorable via [CpCoCp NMe2H ] + than via [(CpCo(Cp-H) NMe2 )] +1/0 . In contrast, the most favorable HER pathway involves the ring-protonated Co(I) species. These results offer mechanistic insights into HER versus PCET bifurcation and establish guiding principles for designing PCET mediators for selective electroreductive transformations.

evolution reactions↗
DOE OSTI2020

First Principles Nonadiabatic Excited-State Molecular Dynamics in NWChem

Computational simulation of non-adiabatic molecular dynamics is an indispensable tool for understanding complex photoinduced processes such as internal conversion, energy transfer, charge separation, and spatial localization of excitons, to name a few. Here, we report an implementation of the fewest-switches surface hopping algorithm in the NWChem computational chemistry program. Here the surface hopping method is combined with linear-response time-dependent density functional theory calculations of adiabatic excited state potential energy surfaces. To treat quantum transitions between arbitrary electronic Born{Oppenheimer states, we have implemented both numerical and analytical differentiation schemes for derivative non-adiabatic couplings. A numerical approach for the time-derivative non-adiabatic couplings together with an analytical method for calculating non-adiabatic coupling vectors is an efficient combination for surface hopping approaches. Additionally, electronic decoherence schemes and a state reassigned unavoided crossings algorithm are also implemented to improve the accuracy of the simulated dynamics and to handle trivial unavoided crossings. We apply our code to study the ultrafast decay of photoexcited benzene, including a detailed analysis of the potential energy surface, population decay time scales, and vibrational coordinates coupled to the excitation dynamics. The development presented in this work is a baseline for future implementations of more sophisticated frameworks for simulating of non-adiabatic molecular dynamics in NWChem.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Combined First-Principles Calculations of Electron–Electron and Electron–Phonon Self-Energies in Condensed Systems

In this work, we present a method to efficiently combine the computation of electron–electron and electron–phonon self-energies, which enables the evaluation of electron–phonon coupling at the G 0 W 0 level of theory for systems with hundreds of atoms. In addition, our approach, which is a generalization of a method recently proposed for molecules [J. Chem. Theory Comput.2018, 14, 6269–6275], enables the inclusion of nonadiabatic and temperature effects at no additional computational cost. We present results for diamond and defects in diamond and discuss the importance of numerically accurate G 0 W 0 band structures to obtain robust predictions of zero point renormalization (ZPR) of band gaps, and of the inclusion of nonadiabatic effects to accurately compute the ZPR of defect states in the band gap.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Protein–Ligand Binding Free-Energy Calculations with ARROW–A Purely First-Principles Parameterized Polarizable Force Field

Protein–ligand binding free-energy calculations using molecular dynamics (MD) simulations have emerged as a powerful tool for in silico drug design. Here, we present results obtained with the ARROW force field (FF)–a multipolar polarizable and physics-based model with all parameters fitted entirely to high-level ab initio quantum mechanical (QM) calculations. ARROW has already proven its ability to determine solvation free energy of arbitrary neutral compounds with unprecedented accuracy. The ARROW FF parameterization is now extended to include coverage of all amino acids including charged groups, allowing molecular simulations of a series of protein–ligand systems and prediction of their relative binding free energies. We ensure adequate sampling by applying a novel technique that is based on coupling the Hamiltonian Replica exchange (HREX) with a conformation reservoir generated via potential softening and nonequilibrium MD. ARROW provides predictions with near chemical accuracy (mean absolute error of ~0.5 kcal/mol) for two of the three protein systems studied here (MCL1 and Thrombin). The third protein system (CDK2) reveals the difficulty in accurately describing dimer interaction energies involving polar and charged species. Overall, for all of the three protein systems studied here, ARROW FF predicts relative binding free energies of ligands with a similar accuracy level as leading nonpolarizable force fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

A Lattice Kinetic Monte Carlo Solver for First-Principles Microkinetic Trend Studies

Here, mean-field microkinetic models in combination with Brønsted–Evans–Polanyi like scaling relations have proven highly successful in identifying catalyst materials with good or promising reactivity and selectivity. Analysis of the microkinetic model by means of lattice kinetic Monte Carlo promises a faithful description of a range of atomistic features involving short-range ordering of species in the vicinity of an active site. In this paper, we use the “fruit fly” example reaction of CO oxidation on fcc(111) transition and coinage metals to motivate and develop a lattice kinetic Monte Carlo solver suitable for the numerically challenging case of vastly disparate rate constants. As a result, we show that for the case of infinitely fast diffusion and absence of adsorbate-adsorbate interaction it is, in fact, possible to match the prediction of the mean-field-theory method and the lattice kinetic Monte Carlo method. As a corollary, we conclude that lattice kinetic Monte Carlo simulations of surface chemical reactions are most likely to provide additional insight over mean-field simulations if diffusion limitations or adsorbate–adsorbate interactions have a significant influence on the mixing of the adsorbates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Heterogeneous Reduction Pathways for Hg(II) Species on Dry Aerosols: A First-Principles Computational Study

Here, the atmospheric lifetime of mercury is greatly impacted by redox chemistry resulting from the high deposition rate of reactive mercury (Hg(II)) compared to elemental mercury (Hg 0 ). Recent laboratory and field studies have observed the reduction of Hg(II) but the chemical mechanism for this reaction has not been identified. Recent laboratory studies have shown that the reduction reaction is heterogeneous and can occur on iron and sodium chloride aerosol surfaces. This study explores the use of density functional theory calculations to discern the reduction pathways of HgCl 2 , HgBr 2 , Hg(NO 3 ) 2 , and HgSO 4 on clean Fe(110), NaCl(100), and NaCl(111) Na surfaces. In doing so, potential energy surfaces have been prepared for the various reduction pathways, indicating that the reduction pathway leading to the production of gas-phase elemental mercury is highly favorable on Fe(110) and NaCl(111) Na . Moreover, the Fe(110) surface requires an external energy source of approximately 0.5 eV to desorb the reduced mercury, whereas the NaCl(111) Na surface requires no energy input. The results indicate that a number of mercury species can be reduced on metallic iron and sodium chloride surfaces, which are known aerosol components, and that a photochemical reaction involving the aerosol surface is likely needed for the reaction to be catalytic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Can Excited State Electronic Coherence Be Tuned via Molecular Structural Modification? A First-Principles Quantum Electronic Dynamics Study of Pyrazolate-Bridged Pt(II) Dimers

Materials and molecular systems exhibiting long-lived electronic coherence can facilitate coherent transport, opening the door to efficient charge and energy transport beyond traditional methods. Recently, signatures of a possible coherent, recurrent electronic motion were identified in femtosecond pump-probe spectroscopy experiments on a binuclear platinum complex, where a persistent periodic beating in the transient absorption signal’s anisotropy was observed. In this study, we investigate the excitonic dynamics that underlie the suspected electronic coherence for a series of binuclear platinum complexes exhibiting a range of interplatinum distances. Results suggest that the long-lived coherence can only result when competitive electronic couplings are in balance. At longer Pt-Pt distances, the electronic couplings between the two halves of the binuclear system weaken, and exciton localization and recombination is favored on short time scales. For short Pt-Pt distances, electronic couplings between the states in the coherent superposition are stronger than the coupling with other excitonic states, leading to long-lived coherence.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

The CP-PAW Code Package for First-Principles Calculations from a User’s Perspective

CP-PAW is a combined electronic structure and ab initio molecular dynamics code to perform mixed quantum and classical simulations of atomistic condensed phase systems, such as solids, liquids, and molecular systems. As the name suggests, the CP-PAW code unifies the all-electron projector augmented-wave (PAW) method with the Car–Parrinello (CP) approach to determine not only the electronic and nuclear ground states of condensed matter but also to study their properties and dynamics. In addition to briefly outlining the underlying theory, the focus will be on the unique aspects of CP-PAW and how to correctly employ them as a user. How to install CP-PAW using the new build system will also be briefly mentioned.

Blöchl, Peter E [Institute for Theoretical Physic↗