Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “potential energy surface”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Stabilization of Charge-Transfer Excited States in Biological Systems: A Computational Focus on the Special Pair in Photosystem II Reaction Centers

Charge-transfer (CT) excited states play an important role in many biological processes. However, many computational approaches often inadequately address the equilibration effects of nuclear and environmental degrees of freedom on these states. One prominent example of systems in which CT states are of utmost importance is reaction centers (RC) in photosystems. Here we use a multiscale approach combined with time-dependent density functional theory to explore the lowest CT excited state of the special pair P D1 –P D2 in the Photosystem II-RC of a cyanobacterium. We find that the nonequilibrium CT excited state resides near the Soret band, making an exciton the lowest-energy excited state. However, accounting for nuclear and state-specific dielectric equilibration along the CT potential energy surface (PES), the CT state P D1 – –P D2 + stabilizes energetically below the excitonic state. Importantly, this underscores the crucial role of state-specific solvation in mapping the PES of CT states, as demonstrated in a simplified dimer model.

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

A Many-Body Perspective of Nuclear Quantum Effects in Aqueous Clusters

Nuclear quantum effects play an important role in the structure and thermodynamics of aqueous systems. Here, by performing a many-body expansion with nuclear-electronic orbital (NEO) theory, we show that proton quantization can give rise to significant energetic contributions for many-body interactions spanning several molecules in single-point energy calculations of water clusters. Although zero-point motion produces a large increase in energy at the one-body level, nuclear quantum effects serve to stabilize higher-order molecular interactions. These results are significant because they demonstrate that nuclear quantum effects play a nontrivial role in many-body interactions of aqueous systems. Our approach also provides a pathway for incorporating nuclear quantum effects into water potential energy surfaces. The NEO approach is advantageous for many-body expansion analyses because it includes nuclear quantum effects directly in the energies.

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

Toward Accurate Spin–Orbit Splittings from Relativistic Multireference Electronic Structure Theory

Most nonrelativistic electron correlation methods can be adapted to account for relativistic effects, as long as the relativistic molecular spinor integrals are available, from either a four-, two-, or one-component mean-field calculation. Furthermore, relativistic multireference correlation methods remain a relatively unexplored area, with mixed evidence regarding the improvements brought by perturbative treatments. We report, for the first time, the implementation of state-averaged four-component relativistic multireference perturbation theories to second and third order based on the driven similarity renormalization group (DSRG). With our methods, named 4c-SA-DSRG-MRPT2 and 3, we find that the dynamical correlation included on top of 4c-CASSCF references can significantly improve the spin-orbit splittings in p-block elements and potential energy surfaces when compared to 4c-CASSCF and 4c-CASPT2 results. We further show that 4c-DSRG-MRPT2 and 3 are applicable to these systems over a wide range of the flow parameter, with systematic improvement from second to third order in terms of both improved error statistics and reduced sensitivity with respect to the flow parameter.

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

Smooth Dispersion Is Physically Appropriate: Assessing and Amending the D4 Dispersion Model

The addition of dispersion corrections to density functionals is essential for accurate energy and geometry predictions. Among them, the D4 scheme is popular due to its low computational cost and high accuracy. However, due to its design, the D4 correction can occasionally lead to anomalies, such as unphysical curvature and bumps in the potential energy surface. We find these anomalies are common in the D4 model, although observable consequences are rarer than in the D3 model for reasons we explain. Nevertheless, we uncover instances of unphysical local minima and stationary points with the D4 scheme and propose two solutions that yield smoother dispersion energy as a function of nuclear position. One is trivial to implement, based on a smoother reparametrization of Gaussian weighting (D4S) to find the effective coordination number. The other replaces Gaussian weighting with soft linear interpolation (D4SL). These new approaches usually remove artificial extremum points, while maintaining accuracy.

Tkachenko, Nikolay V↗
DOE OSTI2025

Visualizing Stereodynamics in Cold Collisions through Shape Resonance Wavefunctions

Shape resonances in cold collisions are often strongly affected by stereodynamics. These resonances can sometimes be assigned by single-channel quantization along the scattering coordinate. However, sensitive steric control of collision implies strong anisotropy of the interaction potential energy surface, which usually leads to coupling among multiple scattering channels in the strongly interacting region. Hence, the resonances might be the result of quantization in a multidimensional space. Using the cold rotationally inelastic collision between para-H 2 (v 1 = 1, j 1 = 2) and HF (v 2 = 0, j 2 = 0) as an example, we analyze four low-lying shape resonances via diagonalization of the full-dimensional Hamiltonian with a stabilization method. While some resonances can indeed be assigned with a single partial wave, others apparently involve more than one scattering channel. Furthermore, a new model based on these resonance wavefunctions is developed to better understand the stereodynamics through shape resonances in cold collisions.

Collisions↗
DOE OSTI2026

The Hidden Path to the Resonance-Stabilized Fulvenallenyl Radical (C 7 H 5 ) via the Bimolecular Reaction of Tricarbon (C 3 , X 1 Σ g + ) with 1,3-Butadiene (C 4 H 6 ; X 1 A g )

Here, we report a novel bimolecular reaction pathway forming the resonance-stabilized fulvenallenyl radical (C 7 H 5 ) via the gas-phase reaction of tricarbon (C 3 , X 1 Σ g + ) with 1,3-butadiene (C 4 H 6 , X 1 A g ). Crossed molecular beam experiments combined with high-level electronic structures and statistical calculations reveal a rich potential energy surface. The reaction proceeds via tricarbon addition to a double bond of 1,3-butadiene, overcoming a 30 kJ mol –1 barrier, followed by ring closure, isomerizations (ring opening/closures, hydrogen shifts), and eventual hydrogen atom loss, yielding the fulvenallenyl radical (p1) almost exclusively via two dominant pathways. Considering the entrance barrier, this reaction may occur in high-temperature environments, like circumstellar envelopes of carbon stars, but is suppressed in colder regions such as molecular clouds and Titan’s atmosphere. These findings highlight the unexpected reactivity of small carbon clusters in shaping the molecular complexity of our universe, from the flicker of a flame to the death of a star.

Aromatic compounds↗
DOE OSTI2026

Signatures of Antisymmetric Vibrations in the Ultrafast Dynamics of Quadrupolar Dyes

Antisymmetric molecular vibrations are central to ultrafast, nonadiabatic photophysical and photochemical processes such as conical intersection dynamics, Herzberg–Teller couplings and, potentially, singlet fission. Direct spectroscopic identification of such vibrations is, however, challenging, since they are typically Raman inactive and affect optical transitions only weakly. Here, we report experimental signatures of vibronic couplings to a high-frequency antisymmetric vibration in the excited state dynamics of a quasi-quadrupolar molecule by ultrafast two-dimensional electronic spectroscopy (2DES). The early time, sub-50 fs 2DES maps reveal an asymmetric peak pattern with characteristic low-energy cross-peaks. We show that these peaks arise from stimulated emission transitions from an anharmonic, double-minimum excited state potential energy surface formed by vibronic coupling to a high-frequency antisymmetric mode. Simulations based on a phenomenological essential state model support the results. Our findings offer a new approach for identifying antisymmetric vibrations in ultrafast 2DES and track excited state wavepacket motion before relaxation washes out the spectra.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Minimization of Disorder as a Key Design Principle for Natural Sizes of Light Harvesting 2 Complexes

The light harvesting 2 (LH2) complex of purple bacteria has excellent energy conversion efficiency. Clarifying the design principle behind such efficiency at the atomistic level is crucial for understanding its structure–function relationship and can be utilized for the design of artificial light harvesting systems. To this end, we conducted comprehensive computational investigation of the dynamical and statistical nature of electronic excited states of pigment molecules in a natural LH2 complex with 9-fold symmetry and its two non-natural in silico analogues with 6- and 12-fold symmetries. To ensure reliable and efficient all-atomistic molecular dynamics simulations, we combined a well established interpolation approach for the construction of the potential energy surface with a neural network machine learning approach. Outcomes of these calculations clarify that non-natural forms of LH2-type complexes have significantly larger quasistatic disorder than those for the natural one. In addition, non-natural systems have more disruptions of the hydrogen bonding, underscoring its crucial role for reducing the disorder. On the other hand, local environmental dynamics are relatively insensitive to the structural changes although there is moderate enhancement in the anharmonic or interatomic components for the synthetic ones. These findings based on all-atomistic simulations provide direct computational evidence that the structure and sizes of natural LH2 complexes are designed to minimize the energetic disorder. We analyze quantitative implications of these for the energy transferring capability of the LH2 complex.

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

Diffractive Imaging of Conical Intersections Amplified by Resonant Infrared Fields

The fate of virtually all photochemical reactions is determined by conical intersections. These are energetically degenerate regions of molecular potential energy surfaces that strongly couple electronic states, thereby enabling fast relaxation channels. Their direct spectroscopic detection relies on weak features that are often buried beneath stronger, less interesting contributions. For azobenzene photoisomerization, a textbook photochemical reaction, we demonstrate how a resonant infrared field can be employed during the conical intersection passage to significantly enhance its coherence signatures in time-resolved Xray diffraction while leaving the product yield intact. Furthermore, this transition-state amplification holds promise to bring signals of conical intersections above the detection threshold.

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

Control of Photoswitching Kinetics with Strong Light–Matter Coupling in a Cavity

Most photochemistry occurs in the regime of weak light–matter coupling, in which a molecule absorbs a photon and then performs photochemistry from its excited state. In the strong coupling regime, enhanced light–matter interactions between an optical field and multiple molecules lead to collective hybrid light–matter states called polaritons. This strong coupling leads to fundamental changes in the nature of the excited states including multi-molecule delocalized excitations, modified potential energy surfaces, and dramatically altered energy levels relative to non-coupled molecules. Here, the effect of strong light–matter coupling on covalent photochemistry has not been well explored. Photoswitches undergo reversible intramolecular photoreactions that can be readily monitored spectroscopically. In this work, we study the effect of strong light–matter coupling on the kinetics of photoswitching within optical cavities. Reproducing prior experiments, photoswitching of spiropyran/merocyanine photoswitches is decelerated in a cavity. Fulgide photoswitches, however, show the opposite effect, with strong coupling accelerating photoswitching. While modified merocyanine switching can be explained by changes in radiative decay rates or the amount of light in the cavity, modified fulgide switching kinetics suggest direct changes to excited-state reaction kinetics.

Cavities↗
DOE OSTI2023

Combining Force Fields and Neural Networks for an Accurate Representation of Chemically Diverse Molecular Interactions

A key goal of molecular modeling is the accurate reproduction of the true quantum mechanical potential energy of arbitrary molecular ensembles with a tractable classical approximation. The challenges are that analytical expressions found in general purpose force fields struggle to faithfully represent the intermolecular quantum potential energy surface at close distances and in strong interaction regimes; that the more accurate neural network approximations do not capture crucial physics concepts, e.g., nonadditive inductive contributions and application of electric fields; and that the ultra-accurate narrowly targeted models have difficulty generalizing to the entire chemical space. We therefore designed a hybrid wide-coverage intermolecular interaction model consisting of an analytically polarizable force field combined with a short-range neural network correction for the total intermolecular interaction energy. Here, we describe the methodology and apply the model to accurately determine the properties of water, the free energy of solvation of neutral and charged molecules, and the binding free energy of ligands to proteins. The correction is subtyped for distinct chemical species to match the underlying force field, to segment and reduce the amount of quantum training data, and to increase accuracy and computational speed. For the systems considered, the hybrid ab initio parametrized Hamiltonian reproduces the two-body dimer quantum mechanics (QM) energies to within 0.03 kcal/mol and the nonadditive many-molecule contributions to within 2%. Simulations of molecular systems using this interaction model run at speeds of several nanoseconds per day.

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

Chemistry of Sugar Formation in the Gas Phase: Following the Activated Aldehyde

Sugars are produced by living organisms, and are required building blocks for life as we know it, which raises the foundational question of how sugars formed in a prebiotic environment. The abiotic formose reaction produces sugars from formaldehyde, but our understanding of its initiation step remains murky, with chemists invoking the concept of an “activated aldehyde” to seed this reaction. Singlet hydroxycarbenes, high-energy isomers of aldehydes, were recently reported to facilitate sugar formation under cold, nonaqueous conditions relevant to interstellar environments. Here, we generate singlet methylhydroxycarbene ( 1 CH 3 –C̈–OH) from the photodissociation of pyruvic acid and experimentally measure its gas-phase reaction with d 4 -acetaldehyde using multiplexed photoionization mass spectrometry. The C 4 H 4 D 4 O 2 isomer d 4 -acetoin is the sole product, which we kinetically link to the reactant CH 3 –C̈–OH, and attribute to a carbonyl-ene formation mechanism. We see no evidence of 3-hydroxybutanal, the C–H insertion product expected in carbene chemistry. Using automated exploration we calculate stationary points on the potential energy surface and report master equation rate coefficients from T = 20–600 K, providing quantitative kinetics of this fast reaction for use in chemical models. The prereactive complex in this reaction is stabilized by both hydrogen bonding and electrophilic carbene C═O interactions. These effects create a short-range dynamical bottleneck for the reaction besides the long- and midrange barrierless bottlenecks. Combined with recent reports of 1 HC̈OH production from methanol photodissociation and pyruvic acid production in cold irradiated ices, this work provides evidence that singlet hydroxycarbene + aldehyde chemistry is a feasible path to prebiotic sugar formation.

aldehydes↗
DOE OSTI2025

Acid-Catalyst is Required for Hydrogenolysis of a Late Metal Hydroxo Complex

The cationic AuIII pincer complex [(tBuPCP)AuIII -OH]OTf (AuOH, tBuPCP = 2,6-(CH2PtBu2)2C6H3) was prepared and characterized. In contrast to hydrogenolysis reactions of other late metal-hydroxides, the hydrogenolysis of AuOH to form the corresponding AuIII -H and water is shown to require an acid catalyst. The kinetic data, collected in acetone-d6, are consistent with a reaction mechanism involving initial protonation to form the AuIII aquo complex, [(tBuPCP)AuIII -OH2](OTf)2 (AuOH2), which was isolated and characterized. Following substitution of the water by coordinating solvent, the dicationic AuIII center undergoes electrophilic activation of H2 to generate the [(tBuPCP)AuIII -H]OTf product (AuH). The mechanism of hydrogenolysis is supported by experimental and computational (DFT) results. A AuIII dihydrogen complex (Au(H2)) was located along the potential energy surface for the reaction and computations indicate that deprotonation of Au(H2) by acetone solvent is barrierless. The bonding and properties of the proposed AuIII -H2 complex are consistent with predictions made for high-valent, electrophilic metal-dihydrogen complexes.

Ardon, Yotam [University of Pennsylvania] (ORCID:0↗
DOE OSTI2026

Ion-Electron Coupling-Driven Redox Behavior in Metal–Organic Frameworks

Redox-active metal–organic frameworks (MOFs) have long been proposed as electronic transport platforms, yet the microscopic origin of their conductivity remains debated. A theoretical demonstration reveals charge transport in a Zn(pyrazole–naphthalene diimide (NDI)) MOF arising not from delocalized band-like states but from redox hopping between discrete linker sites. Using ab initio molecular dynamics simulations combined with electronic structure analysis, we established a direct link among electron injection, structural reorganization, and transport. Electron accumulation proceeds sequentially and site-selectively from imide and carbonyl groups of the NDI core progressively involving pyrazole N atoms at higher reduction states, through a hierarchy of redox-active sites. In contrast, Zn nodes remain essentially redox-inactive, which confirms their structural role. Density-of-states analysis corroborates a transport regime dominated by linker-centered states with evolving p-character upon reduction, resulting in dynamically reconfigured conduction networks. Real-time trajectories reveal anisotropic linker-to-linker electron transfer modulated by counterion coordination. This cooperative ion–electron regime emerges from potential energy surface collapse into a single low-barrier transition (ΔG ‡ ≈ 45 meV), where ionic and electronic motions evolve adiabatically on the same free-energy landscape. Elucidating redox conductivity in Zn(pyrazole–NDI) MOFs provides a theoretical framework for use in neuromorphic computing and related technologies.

Charge transfer↗
DOE OSTI2021

Nonequilibrium charge-density-wave order beyond the thermal limit

The interaction of many-body systems with intense light pulses may lead to novel emergent phenomena far from equilibrium. Recent discoveries, such as the optical enhancement of the critical temperature in certain superconductors and the photo-stabilization of hidden phases, have turned this field into an important research frontier. Here, we demonstrate nonthermal charge-density-wave (CDW) order at electronic temperatures far greater than the thermodynamic transition temperature. Using time- and angle-resolved photoemission spectroscopy and time-resolved X-ray diffraction, we investigate the electronic and structural order parameters of an ultrafast photoinduced CDW-to-metal transition. Tracking the dynamical CDW recovery as a function of electronic temperature reveals a behaviour markedly different from equilibrium, which we attribute to the suppression of lattice fluctuations in the transient nonthermal phonon distribution. A complete description of the system’s coherent and incoherent order-parameter dynamics is given by a time-dependent Ginzburg-Landau framework, providing access to the transient potential energy surfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Chiral photochemistry of achiral molecules

Chirality is a molecular property governed by the topography of the potential energy surface (PES). Thermally achiral molecules interconvert rapidly when the interconversion barrier between the two enantiomers is comparable to or lower than the thermal energy, in contrast to thermally stable chiral configurations. In principle, a change in the PES topography on the excited electronic state may diminish interconversion, leading to electronically prochiral molecules that can be converted from achiral to chiral by electronic excitation. Here we report that this is the case for two prototypical examples – cis-stilbene and cis-stiff stilbene. Both systems exhibit unidirectional photoisomerization for each enantiomer as a result of their electronic prochirality. We simulate an experiment to demonstrate this effect in cis-stilbene based on its interaction with circularly polarized light. Our results highlight the drastic change in chiral behavior upon electronic excitation, opening up the possibility for asymmetric photochemistry from an effectively nonchiral starting point.

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

Step-by-step state-selective tracking of fragmentation dynamics of water dications by momentum imaging

Abstract The double photoionization of a molecule by one photon ejects two electrons and typically creates an unstable dication. Observing the subsequent fragmentation products in coincidence can reveal a surprisingly detailed picture of the dynamics. Determining the time evolution and quantum mechanical states involved leads to deeper understanding of molecular dynamics. Here in a combined experimental and theoretical study, we unambiguously separate the sequential breakup via D + + OD + intermediates, from other processes leading to the same D + + D + + O final products of double ionization of water by a single photon. Moreover, we experimentally identify, separate, and follow step by step, two pathways involving the b 1 Σ + and a 1 Δ electronic states of the intermediate OD + ion. Our classical trajectory calculations on the relevant potential energy surfaces reproduce well the measured data and, combined with the experiment, enable the determination of the internal energy and angular momentum distribution of the OD + intermediate.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Architector for high-throughput cross-periodic table 3D complex building

Abstract Rare-earth and actinide complexes are critical for a wealth of clean-energy applications. Three-dimensional (3D) structural generation and prediction for these organometallic systems remains a challenge, limiting opportunities for computational chemical discovery. Here, we introduce Architector , a high-throughput in-silico synthesis code for s-, p-, d-, and f-block mononuclear organometallic complexes capable of capturing nearly the full diversity of the known experimental chemical space. Beyond known chemical space, Architector performs in-silico design of new complexes including any chemically accessible metal-ligand combinations. Architector leverages metal-center symmetry, interatomic force fields, and tight binding methods to build many possible 3D conformers from minimal 2D inputs including metal oxidation and spin state. Over a set of more than 6,000 x-ray diffraction (XRD)-determined complexes spanning the periodic table, we demonstrate quantitative agreement between Architector-predicted and experimentally observed structures. Further, we demonstrate out-of-the box conformer generation and energetic rankings of non-minimum energy conformers produced from Architector , which are critical for exploring potential energy surfaces and training force fields. Overall, Architector represents a transformative step towards cross-periodic table computational design of metal complex chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗