Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Oppenheimer”

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 109 records · Page 6

Proton Transfer Dynamics in the Aprotic Proton Accepting Solvent 1-Methylimidazole

The dynamics of proton transfer to the aprotic solvent 1-methylimidazole (MeIm, proton acceptor) from the photoacid 8-hydroxypyrene-1,3,6-trisulfonic acid trisodium salt (HPTS) were investigated using fast fluorescence measurements. The closely related molecule, 8- methoxypyrene-1,3,6-trisulfonic acid trisodium salt (MPTS), which is not a photoacid, was also studied for comparison. Following optical excitation, the wavelength-dependent population dynamics of HPTS in MeIm resulting from the deprotonation process were collected over the entire fluorescence emission window. Analysis of the time-dependent fluorescence spectra revealed four distinct fluorescence bands that appear and decay on different time scales. We label these four states as: protonated (P), associated I (A I ), associated II (A II ), and deprotonated (D). We find that the simple kinetic scheme of P → A I → A II → D is not consistent with the data. Instead, the kinetic scheme that describes the data has P decaying into A I , which mainly goes on to deprotonation (D), but AI can also feed into A II . A II can return to A I or decay to the ground state, but does not deprotonate within experimental error. Finally, quantum chemistry and excitedstate QM/MM Born-Oppenheimer molecular dynamics simulations indicate that A I and A II are two H-bonding conformations of MeIm to the HPTS hydroxyl, axial and equatorial, respectively.

14 SOLAR ENERGY↗

Wetting of a Hydrophobic Surface: Far-IR Action Spectroscopy and Dynamics of Microhydrated Naphthalene

The interaction of water and polycyclic aromatic hydrocarbons is of fundamental importance in areas as diverse as materials science and atmospheric and interstellar chemistry. The interplay between hydrogen bonding and dipole–π interactions results in subtle dynamics that are challenging to describe from first principles. Here, we employ far-IR action vibrational spectroscopy with the infrared free-electron laser FELIX to investigate naphthalene with one to three water molecules. We observe diffuse bands associated with intermolecular vibrational modes that serve as direct probes of the loose binding of water to the naphthalene surface. These signatures are poorly reproduced by static DFT or Møller–Plesset computations. Instead, a rationalization is achieved through Born–Oppenheimer Molecular Dynamics simulations, revealing the active mobility of water over the surface, even at low temperatures. Therefore, our work provides direct insights into the wetting interactions associated with shallow potential energy surfaces while simultaneously demonstrating a solid experimental–computational framework for their investigation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Direct Observation of the ππ* to nπ* Transition in 2-Thiouracil via Time-Resolved NEXAFS Spectroscopy

The photophysics of nucleobases has been the subject of both theoretical and experimental studies over the past decades due to the challenges posed by resolving the steps of their radiationless relaxation dynamics, which cannot be described in the framework of the Born–Oppenheimer approximation (BOA). In this context, the ultrafast dynamics of 2-thiouracil has been investigated with a time-resolved NEXAFS study at the Free Electron Laser FLASH. Near Edge X-ray Absorption Fine Structure spectroscopy (NEXAFS) can be used to observe electronic transitions in ultrafast molecular relaxation. We performed time-resolved UV-pump/X-ray probe absorption measurements at the sulfur 2s (L1) and 2p (L2/3) edges. We are able to identify absorption features corresponding to the S2 (ππ*) and S1 (nπ*) electronic states. We observe a delay of 102 ± 11 fs in the population of the nπ* state with respect to the initial optical excitation and interpret the delay as the time scale for the S2 → S1 internal conversion. We furthermore identify oscillations in the absorption signal that match a similar observation in our previous X-ray photoelectron spectroscopy study on the same molecule.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Yttrium Contenting Compositionally Complex Medium-Entropy Li-Garnet Electrolyte with Improved Ionic Conductivity

The compositionally complex medium/high-entropy design concept can greatly expand the categories and affect the properties of the materials. With such a designing concept, a medium-entropy Li-garnet electrolyte with appropriate yttrium content (formula Li 6.6 La 3 ZrNb 0.3 Ta 0.3 Hf 0.3 Y 0.1 O 12 ) shows a record-high ionic conductivity of ∼5.7 × 10 –4 S/cm, the highest reported for any single-site substituted high/medium-entropy Li-garnet. The assembled Li metal symmetric cells also show stable long-term cycling (0.1 mA/cm 2 for over 200 h). Neutron powder diffraction and Rietveld refinement results indicate that a competing conduction mechanism between (1) occupancy on high mobility of 96h sites and (2) the associated site vacancies and the bond length requires an appropriate content of Y for enhanced ionic conductivity. Li-ion hopping through the bottleneck can also contribute to the conductivity. Finally, density functional theory and Born–Oppenheimer molecular dynamics simulations also indicate the high mobility and number of hopping transitions of Li ions, contributing to the high ionic conductivity.

Li-garnet↗

Highly Active Hydrogen Evolution Reaction (HER) Catalysts Formed by Energetic Pt n Cluster Deposition: Deposition Dynamics and the HER Mechanism

Mass-selected Pt n + (n ≤ 7) were deposited at variable energies on highly oriented pyrolytic graphite (HOPG), creating highly active hydrogen evolution reaction (HER) electrocatalysts. HER mass activities were ~2 to >10 times higher than those for the surface atoms in bulk Pt and for Pt n deposited on several other supports. Thus, high activity reflects the Pt-C structures formed by energetic Pt n -HOPG impacts, in addition to high Pt surface availability. The Pt n /HOPG electrodes were probed by X-ray photoelectron spectroscopy, low energy ion scattering, and electron microscopy. Born-Oppenheimer molecular dynamics (BOMD) was used to simulate Pt n - HOPG impacts, revealing the types of structures formed at different energies, then DFT was used to probe their most important HER pathways. For low deposition energies, the Pt n deposit onto the HOPG surface with sub-unit sticking probability, aggregating at defects. With increasing deposition energy, the sticking probability initially decreases, then rises to unity as subplantation and defect creation allow formation of strongly bonded platinum-carbon structures. Barriers for HER on these structures were found to be low and weakly dependent on Pt n size, consistent with experiment. The activities were highest for small covalently-bonded Pt-C structures created at high deposition energies. The larger aggregated structures formed at low energies were less active, but still substantially better than the bulk Pt surface monolayer. The catalysts were stable in repeated potential cycling at reducing potentials, but electrodes containing subplanted Pt became more active when scanned to oxidizing potentials, due to emergence of subplanted Pt onto the surface.

08 HYDROGEN↗

Time-resolved relaxation and fragmentation of polycyclic aromatic hydrocarbons investigated in the ultrafast $\mathrm{XUV}$-$\mathrm{IR}$ regime

Polycyclic aromatic hydrocarbons (PAHs) play an important role in interstellar chemistry and are subject to high energy photons that can induce excitation, ionization, and fragmentation. Previous studies have demonstrated electronic relaxation of parent PAH monocations over 10–100 femtoseconds as a result of beyond-Born-Oppenheimer coupling between the electronic and nuclear dynamics. Here, we investigate three PAH molecules: fluorene, phenanthrene, and pyrene, using ultrafast XUV and IR laser pulses. Simultaneous measurements of the ion yields, ion momenta, and electron momenta as a function of laser pulse delay allow a detailed insight into the various molecular processes. We report relaxation times for the electronically excited PAH * , PAH +* and PAH 2+* states, and show the time-dependent conversion between fragmentation pathways. Additionally, using recoil-frame covariance analysis between ion images, we demonstrate that the dissociation of the PAH 2+ ions favors reaction pathways involving two-body breakup and/or loss of neutral fragments totaling an even number of carbon atoms.

79 ASTRONOMY AND ASTROPHYSICS↗

Ultrafast lattice disordering can be accelerated by electronic collisional forces

In the prevalent picture of ultrafast structural phase transitions, the atomic motion occurs in a slowly varying potential energy surface determined adiabatically by the fast electrons. However, this ignores non-conservative forces caused by electron-lattice collisions, which can significantly influence atomic motion. Most ultrafast techniques only probe the average structure and are less sensitive to random displacements, and therefore do not detect the role played by non-conservative forces in phase transitions. Here we show that the lattice dynamics of the prototypical insulator-to-metal transition of VO 2 cannot be described by a potential energy alone. We use the sample temperature to control the preexisting lattice disorder before ultrafast photoexcitation across the phase transition and our ultrafast diffuse scattering experiments show that the fluctuations characteristic of the rutile metal develop equally fast (120 fs) at initial temperatures of 100 K and 300 K. This indicates that additional non-conservative forces are responsible for the increased lattice disorder. In conclusion, these results highlight the need for more sophisticated descriptions of ultrafast phenomena beyond the Born-Oppenheimer approximation as well as ultrafast probes of spatial fluctuations beyond the average unit cell measured by diffraction.

36 MATERIALS SCIENCE↗

Attosecond response of molecules to impulsive ionization

When matter interacts with energetic radiation it can undergo sudden, or impulsive, ionization. This process can drive chemical change and occurs widely in space and planetary atmospheres, yet its comprehensive description challenges our current theoretical and computational capabilities as it requires advanced treatment of electron correlation and nonadiabatic dynamics beyond the Born–Oppenheimer approximation. Here, in this study, we measure the response of the para-aminophenol molecule to sudden ionization. Using attosecond X-ray absorption spectroscopy, we resolve the ultrafast dynamics of the ionized molecule with atomic precision. A subfemtosecond decay corresponds to states undergoing non-radiative decay, whereas few-femtosecond oscillatory signatures are associated with electronic wavepacket motion in stable cation states that later couple to nuclear motion. We compare our measurement with state-of-the-art computational modelling, qualitatively reproducing the observed response across multiple timescales. These results provide a benchmark for computational models of sudden ionization and ultrafast charge motion in matter.

Driver, Taran [SLAC National Accelerator Laborator↗

Real-time degradation dynamics of hydrated per- and polyfluoroalkyl substances (PFASs) in the presence of excess electrons

Per- and polyfluoroalkyl substances (PFASs) are synthetic chemicals that are harmful to both the environment and human health. In this work, using self-interaction-corrected Born–Oppenheimer molecular dynamics simulations, we provide the first real-time assessment of PFAS degradation in the presence of excess electrons. In particular, we show that the initial phase of the degradation involves the transformation of an alkane-type C–C bond into an alkene-type C[double bond, length as m-dash]C bond in the PFAS molecule, which is initiated by the trans elimination of fluorine atoms bonded to these adjacent carbon atoms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-fidelity first principles nonadiabaticity: diabatization, analytic representation of global diabatic potential energy matrices, and quantum dynamics

We report nonadiabatic dynamics, which goes beyond the Born–Oppenheimer approximation, has increasingly been shown to play an important role in chemical processes, particularly those involving electronically excited states. Understanding multistate dynamics requires rigorous quantum characterization of both electronic and nuclear motion. However, such first principles treatments of multi-dimensional systems have so far been rather limited due to the lack of accurate coupled potential energy surfaces and difficulties associated with quantum dynamics. In this Perspective, we review recent advances in developing high-fidelity analytical diabatic potential energy matrices for quantum dynamical investigations of polyatomic uni- and bi-molecular nonadiabatic processes, by machine learning of high-level ab initio data. Special attention is paid to methods of diabatization, high fidelity construction of multi-state coupled potential energy surfaces and property surfaces, as well as quantum mechanical characterization of nonadiabatic nuclear dynamics. To illustrate the tremendous progress made by these new developments, several examples are discussed, in which direct comparison with quantum state resolved measurements led to either confirmation of the observation or sometimes reinterpretation of the experimental data. The insights gained in these prototypical systems greatly advance our understanding of nonadiabatic dynamics in chemical systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Polaron absorption in aligned conjugated polymer films: breakdown of adiabatic treatments and going beyond the conventional mid-gap state model

This study provides the first experimental polarized intermolecular and intramolecular optical absorption components of field-induced polarons in regioregular poly(3-hexylthiophene-2,5-diyl), rr-P3HT, a polymer semiconductor. Highly aligned rr-P3HT thin films were prepared by a high temperature shear-alignment process that orients polymer backbones along the shearing direction. rr-P3HT in-plane molecular orientation was measured by electron diffraction, and out-of-plane orientation was measured through series of synchrotron X-ray scattering techniques. Then, with molecular orientation quantified, polarized charge modulation spectroscopy was used to probe mid-IR polaron absorption in the ℏω = 0.075 — 0.75 eV range and unambiguously assign intermolecular and intramolecular optical absorption components of hole polarons in rr-P3HT. This data represents the first experimental quantification of these polarized components and allowed long-standing theoretical predictions to be compared to experimental results. The experimental data is discrepant with predictions of polaron absorption based on an adiabatic framework that works under the Born–Oppenheimer approximation, but the data is entirely consistent with a more recent nonadiabatic treatment of absorption based on a modified Holstein Hamiltonian. This nonadiabatic treatment was used to show that both intermolecular and intramolecular polaron coherence break down at length scales significantly smaller than estimated structural coherence in either direction. Importantly, this strongly suggests that polaron delocalization is fundamentally limited by energetic disorder in rr-P3HT.

36 MATERIALS SCIENCE↗

Exploring non-adiabaticity to CO reduction reaction through ab initio molecular dynamics simulation

Non-adiabatic chemical reaction refers to the electronic excitation during reactions. This effect cannot be modeled by the ground-state Born-Oppenheimer molecular dynamics (BO-MD), where the electronic structure is at the ground state for every step of ions' movement. Although the non-adiabatic effect has been explored extensively in gas phase reactions, its role in electrochemical reactions, such as water splitting and CO 2 reduction, in electrolyte has been rarely explored. On the other hand, electrochemical reactions usually involve electron transport; thus, a non-adiabatic process can naturally play a significant role. In this work, using one-step CO 2 reduction as an example, we investigated the role of the non-adiabatic effect in the reaction. The reaction barriers were computed by adiabatic BO-MD and non-adiabatic real-time time dependent density functional theory (rt-TDDFT). We found that by including the non-adiabatic effect, rt-TDDFT could increase the reaction barrier up to 6% compared to the BO-MD calculated barrier when the solvent model is used to represent water. Simulations were carried out using explicit water molecules around the reaction site under different overpotentials, and similar non-adiabatic effects were found.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

How good are polarizable and flexible models for water: Insights from a many-body perspective

We present a systematic analysis of state-of-the-art polarizable and flexible water models from a many-body perspective, with a specific focus on their ability to represent the Born–Oppenheimer potential energy surface of water from the gas to the liquid phase. Using coupled cluster data in the completed basis set limit as a reference, we examine the accuracy of the polarizable models in reproducing individual many-body contributions to interaction energies and harmonic frequencies of water clusters and compare their performance with that of MB-pol, an explicit many-body model that has been shown to correctly predict the properties of water across the entire phase diagram. Based on these comparisons, we use MB-pol as a reference to analyze the ability of the polarizable models to reproduce the energy landscape of liquid water under ambient conditions. We find that, while correctly reproducing the energetics of minimum-energy structures, the polarizable models examined in this study suffer from inadequate representations of many-body effects for distorted configurations. To investigate the role played by geometry-dependent representations of 1-body charge distributions in reproducing coupled cluster data for both interaction and many-body energies, we introduce a simplified version of MB-pol that adopts fixed atomic charges and demonstrate that the new model retains the same accuracy as the original MB-pol model. Based on the analyses presented in this study, we believe that future developments of both polarizable and explicit many-body models should continue in parallel and would benefit from synergistic efforts aimed at integrating the best aspects of the two theoretical/computational frameworks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental and theoretical study of energy transfer in a chromophore triad: What makes modeling dynamics successful?

Simulation of electronic dynamics in realistically large molecular systems is a demanding task that has not yet achieved the same level of quantitative prediction already realized for its static counterpart. This is particularly true for processes occurring beyond the Born–Oppenheimer regime. Non-adiabatic molecular dynamics (NAMD) simulations suffer from two convoluted sources of error: numerical algorithms for dynamics and electronic structure calculations. While the former has gained increasing attention, particularly addressing the validity of ad hoc methodologies, the effect of the latter remains relatively unexplored. Indeed, the required accuracy for electronic structure calculations to reach quantitative agreement with experiment in dynamics may be even more strict than that required for static simulations. In this study, we address this issue by modeling the electronic energy transfer in a donor–acceptor–donor (D–A–D) molecular light harvesting system using fewest switches surface hopping NAMD simulations. In the studied system, time-resolved experimental measurements deliver complete information on spectra and energy transfer rates. Subsequent modeling shows that the calculated electronic transition energies are “sufficiently good” to reproduce experimental spectra but produce over an order of magnitude error in simulated dynamical rates. We further perform simulations using artificially shifted energy gaps to investigate the complex relationship between transition energies and modeled dynamics to understand factors affecting non-radiative relaxation and energy transfer rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Comparing (stochastic-selection) ab initio multiple spawning with trajectory surface hopping for the photodynamics of cyclopropanone, fulvene, and dithiane

Ab Initio Multiple Spawning (AIMS) simulates the excited-state dynamics of molecular systems by representing nuclear wavepackets in a basis of coupled traveling Gaussian functions, called trajectory basis functions (TBFs). New TBFs are spawned when nuclear wavepackets enter regions of strong nonadiabaticity, permitting the description of non-Born–Oppenheimer processes. The spawning algorithm is simultaneously the blessing and the curse of the AIMS method: it allows for an accurate description of the transfer of nuclear amplitude between different electronic states, but it also dramatically increases the computational cost of the AIMS dynamics as all TBFs are coupled. Recently, a strategy coined stochastic-selection AIMS (SSAIMS) was devised to limit the ever-growing number of TBFs and tested on simple molecules. In this work, we use the photodynamics of three different molecules—cyclopropanone, fulvene, and 1,2-dithiane—to investigate (i) the potential of SSAIMS to reproduce reference AIMS results for challenging nonadiabatic dynamics, (ii) the compromise achieved by SSAIMS in obtaining accurate results while using the smallest average number of TBFs as possible, and (iii) the performance of SSAIMS in comparison to the mixed quantum/classical method trajectory surface hopping (TSH)—both in terms of its accuracy and computational cost. We show that SSAIMS can accurately reproduce the AIMS results for the three molecules considered at a much cheaper computational cost, often close to that of TSH. We deduce from these tests that an overlap-based criterion for the stochastic-selection process leads to the best agreement with the reference AIMS dynamics for the smallest average number of TBFs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nuclear–electronic orbital approach to quantization of protons in periodic electronic structure calculations

The nuclear–electronic orbital (NEO) method is a well-established approach for treating nuclei quantum mechanically in molecular systems beyond the usual Born–Oppenheimer approximation. In this work, we present a strategy to implement the NEO method for periodic electronic structure calculations, particularly focused on multicomponent density functional theory (DFT). The NEO-DFT method is implemented in an all-electron electronic structure code, FHI-aims, using a combination of analytical and numerical integration techniques as well as a resolution of the identity scheme to enhance computational efficiency. After validating this implementation, proof-of-concept applications are presented to illustrate the effects of quantized protons on the physical properties of extended systems, such as two-dimensional materials and liquid–semiconductor interfaces. Specifically, periodic NEO-DFT calculations are performed for a trans-polyacetylene chain, a hydrogen boride sheet, and a titanium oxide–water interface. The zero-point energy effects of the protons as well as electron–proton correlation are shown to noticeably impact the density of states and band structures for these systems. These developments provide a foundation for the application of multicomponent DFT to a wide range of other extended condensed matter systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The multichannel i -propyl + O2 reaction system: A model of secondary alkyl radical oxidation

The i-propyl + O2 reaction mechanism has been investigated by definitive quantum chemical methods to establish this system as a benchmark for the combustion of secondary alkyl radicals. Focal point analyses extrapolating to the ab initio limit were performed based on explicit computations with electron correlation treatments through coupled cluster single, double, triple, and quadruple excitations and basis sets up to cc-pV5Z. The rigorous coupled cluster single, double, and triple excitations/cc-pVTZ level of theory was used to fully optimize all reaction species and transition states, thus, removing some substantial flaws in reference geometries existing in the literature. The vital i-propylperoxy radical (MIN1) and its concerted elimination transition state (TS1) were found 34.8 and 4.4 kcal mol−1 below the reactants, respectively. Two β-hydrogen transfer transition states (TS2, TS2′) lie above the reactants by (1.4, 2.5) kcal mol−1 and display large Born–Oppenheimer diagonal corrections indicative of nearby surface crossings. An α-hydrogen transfer transition state (TS5) is discovered 5.7 kcal mol−1 above the reactants that bifurcates into equivalent α-peroxy radical hanging wells (MIN3) prior to a highly exothermic dissociation into acetone + OH. The reverse TS5 → MIN1 intrinsic reaction path also displays fascinating features, including another bifurcation and a conical intersection of potential energy surfaces. An exhaustive conformational search of two hydroperoxypropyl (QOOH) intermediates (MIN2 and MIN3) of the i-propyl + O2 system located nine rotamers within 0.9 kcal mol−1 of the corresponding lowest-energy minima.

Chemistry↗

Constrained nuclear–electronic orbital method for periodic density functional theory: Application to H 2 chemisorption on Si(001) surfaces

The nuclear–electronic orbital (NEO) method provides a powerful computational framework for incorporating nuclear quantum effects (NQE) in electronic structure calculations beyond the Born–Oppenheimer approximation. By incorporating additional constraints to the position operator on quantum particles like protons, the NEO method enables calculation of effective potential that accounts for NQE. Here, in this work, we present a new constrained NEO (cNEO) formulation for density functional theory (cNEO-DFT) calculations in the context of extended periodic systems. Using the nudged elastic band method, we discuss an application of the cNEO-DFT approach to studying the adsorption of a hydrogen molecule on the Si(001) surfaces. The calculation shows how NQE impacts the reaction energetics. The proton density changes are computed along the reaction pathways. This work demonstrates the capability of the new cNEO-DFT method to study a wide range of chemical processes, such as surface reactions where the quantum nature of light atoms like protons is non-negligible.

Chemical processes↗