Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational Chemistry”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Atmospheric chemistry of the reaction ClO + O2 reversible reaction ClO (center dot) O2: Where it stands, what needs to be done, and why?

Possible existence and chemistry of ClO (center dot) O2 was originally proposed to explain the Norrish-Neville effect that O2 suppresses chlorine photosensitized loss of ozone. It was also thought that ClO (center dot) O2 might have some atmospheric chemistry significance. Recently, doubts have been cast on this proposal, because certain laboratory data seem to imply that the equilibrium constant of the title reaction is so small that ClO (center dot) O2 may be too unstable to matter. However, those data create only a superficial illusion to that effect, because on a closer analysis they do not disprove a moderately stable and chemically significant ClO (center dot) O2. Furthermore, our state-of-the-science accurate computational chemistry calculations also suggest that ClO (center dot) O2 may be a weakly bound ClOOO radical with a reactive (2)A ground electronic state. There is therefore a need to design and perform definitive experimental tests of the existence and chemistry of the ClO (center dot) O2 species, which we discuss and which have the potential to mediate the chlorine-catalyzed stratospheric ozone depletion.

Prasad, Sheo S.↗
DOE OSTI2024

Intrinsic Bonding and Reactivity of Actinide Clusters Poster

The effect of electronic structure on chemical bonding and reactivity is one of the most important questions in the field of chemistry. Insight into chemical bond formation involving elements that contain f electrons (i.e., the lanthanides and actinides) is critical for solving specific technical challenges, such as development of advanced nuclear fuel cycles and efficient rare earth separations, as well as a broader understanding of bonding across the entire periodic table. The primary technique for studying the fundamentals of lanthanide and actinide bonding is crystallization of metal-ligand compounds followed by structural determination with X-rays and comparison to theory with computational chemistry. While this procedure has yielded significant insights for a wide variety of elements, it is challenging to apply to the heavier trans-uranic actinides such as berkelium and californium because of limited material availability, significant radioactivity, and poorly understood chemical reactivity of these elements. We hypothesize that by forming actinide-ligand clusters in a mass spectrometer and probing their reactivity in the gas-phase, many of the challenges associated with studying the heavier trans-uranic elements can be circumvented. This project will develop a capability to address a wide range of actinide bonding and reactivity phenomena, increasing Idaho National Laboratory’s capability to address fundamental questions at the forefront of heavy element chemical and molecular science.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2011

Development of the Orion Crew Module Static Aerodynamic Database: Hypersonic - Part 1

The Orion aerodynamic database provides force and moment coefficients given the velocity, attitude, configuration, etc. of the Crew Exploration Vehicle (CEV). The database is developed and maintained by the NASA CEV Aerosciences Project team from computational and experimental aerodynamic simulations. The database is used primarily by the Guidance, Navigation, and Control (GNC) team to design vehicle trajectories and assess flight performance. The initial hypersonic re-entry portion of the Crew Module (CM) database was developed in 2006. Updates incorporating additional data and improvements to the database formulation and uncertainty methodologies have been made since then. This paper details the process used to develop the CM database, including nominal values and uncertainties, for Mach numbers greater than 8 and angles of attack between 140deg and 180deg. The primary available data are more than 1000 viscous, reacting gas chemistry computational simulations using both the Laura and Dplr codes, over a range of Mach numbers from 2 to 37 and a range of angles of attack from 147deg to 172deg. Uncertainties were based on grid convergence, laminar-turbulent solution variations, combined altitude and code-to-code variations, and expected heatshield asymmetry. A radial basis function response surface tool, NEAR-RS, was used to fit the coefficient data smoothly in a velocity-angle-of-attack space. The resulting database is presented and includes some data comparisons and a discussion of the predicted variation of trim angle of attack and lift-to-drag ratio. The database provides a variation in trim angle of attack on the order of +/-2deg, and a range in lift-to-drag ratio of +/-0.035 for typical vehicle flight conditions.

Bibb, Karen L.↗
NASA NTRS1993

Boost-phase discrimination research

The final report describes the combined work of the Computational Chemistry and Aerothermodynamics branches within the Thermosciences Division at NASA Ames Research Center directed at understanding the signatures of shock-heated air. Considerable progress was made in determining accurate transition probabilities for the important band systems of NO that account for much of the emission in the ultraviolet region. Research carried out under this project showed that in order to reproduce the observed radiation from the bow shock region of missiles in their boost phase it is necessary to include the Burnett terms in the constituent equation, account for the non-Boltzmann energy distribution, correctly model the NO formation and rotational excitation process, and use accurate transition probabilities for the NO band systems. This work resulted in significant improvements in the computer code NEQAIR that models both the radiation and fluid dynamics in the shock region.

Langhoff, Stephen R.↗
DOE OSTI2024

Distinct vibrational motions promote disparate excited-state decay pathways in cofacial perylenediimide dimers

A complex interplay of structural, electronic, and vibrational degrees of freedom underpins the fate of molecular excited states. Organic assemblies exhibit a myriad of excited-state decay processes, such as symmetry-breaking charge separation (SB-CS), excimer (EX) formation, singlet fission, and energy transfer. Recent studies of cofacial and slip-stacked perylene-3,4:9,10-bis(dicarboximide) (PDI) multimers demonstrate that slight variations in core substituents and H- or J-type aggregation can determine whether the system follows an SB-CS pathway or an EX one. However, questions regarding the relative importance of structural properties and molecular vibrations in driving the excited-state dynamics remain. Here, we use a combination of two-dimensional electronic spectroscopy, femtosecond stimulated Raman spectroscopy, and quantum chemistry computations to compare the photophysics of two PDI dimers. The dimer with 1,7-bis(pyrrolidin-1′-yl) substituents (5PDI2) undergoes ultrafast SB-CS from a photoexcited mixed state, while the dimer with bis-1,7-(3′,5′-di-t-butylphenoxy) substituents (PPDI2) rapidly forms an EX state. Examination of their quantum beating features reveals that SB-CS in 5PDI2 is driven by the collective vibronic coupling of two or more excited-state vibrations. In contrast, we observe signatures of low-frequency vibrational coherence transfer during EX formation by PPDI2, which aligns with several previous studies. We conclude that key electronic and structural differences between 5PDI2 and PPDI2 determine their markedly different photophysics.

Chemistry↗
DOE OSTI2024

Actuating Liquid Crystals Rapidly and Reversibly by Using Chemical Catalysis

Abstract Microtubules and catalytic motor proteins underlie the microscale actuation of living materials, and they have been used in reconstituted systems to harness chemical energy to drive new states of organization of soft matter (e.g., liquid crystals (LCs)). Such materials, however, are fragile and challenging to translate to technological contexts. Rapid (sub‐second) and reversible changes in the orientations of LCs at room temperature using reactions between gaseous hydrogen and oxygen that are catalyzed by Pd/Au surfaces are reported. Surface chemical analysis and computational chemistry studies confirm that dissociative adsorption of H 2 on the Pd/Au films reduces preadsorbed O and generates 1 ML of adsorbed H, driving nitrile‐containing LCs from a perpendicular to a planar orientation. Subsequent exposure to O 2 leads to oxidation of the adsorbed H, reformation of adsorbed O on the Pd/Au surface, and a return of the LC to its initial orientation. The roles of surface composition and reaction kinetics in determining the LC dynamics are described along with a proof‐of‐concept demonstration of microactuation of beads. These results provide fresh ideas for utilizing chemical energy and catalysis to reversibly actuate functional LCs on the microscale.

Chemistry↗
DOE OSTI2024

Ice Phase Classification Made Easy with Score-Based Denoising

Accurate identification of ice phases is essential for understanding various physicochemical phenomena. However, such classification for structures simulated with molecular dynamics is complicated by the complex symmetries of ice polymorphs and thermal fluctuations. For this purpose, both traditional order parameters and data-driven machine learning approaches have been employed, but they often rely on expert intuition, specific geometric information, or large training data sets. In this work, we present an unsupervised phase classification framework that combines a score-based denoiser model with a subsequent model-free classification method to accurately identify ice phases. Further, the denoiser model is trained on perturbed synthetic data of ideal reference structures, eliminating the need for large data sets and labeling efforts. The classification step utilizes the smooth overlap of atomic position (SOAP) descriptors as the atomic fingerprint, ensuring Euclidean symmetries and transferability to various structural systems. Our approach achieves a remarkable 100% accuracy in distinguishing ice phases of test trajectories using only seven ideal reference structures of ice phases as model inputs. This demonstrates the generalizability of the score-based denoiser model in facilitating phase identification for complex molecular systems. The proposed classification strategy can be broadly applied to investigate structural evolution and phase identification for a wide range of materials, offering new insights into the fundamental understanding of water and other complex systems.

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

Improving Bond Dissociations of Reactive Machine Learning Potentials through Physics-Constrained Data Augmentation

In the field of computational chemistry, predicting bond dissociation energies (BDEs) presents well-known challenges, particularly due to the multireference character of reactive systems. Many chemical reactions involve configurations where single-reference methods fall short, as the electronic structure can significantly change during bond breaking. As generating training data for partially broken bonds is a challenging task, even state-of-the-art reactive machine learning interatomic potentials (MLIPs) often fail to predict reliable BDEs and smooth dissociation curves. By contrast, simple and inexpensive physics-based models, such as the well-established Morse potential, do not suffer from any such limitations. This work leverages the Morse potential to improve reactive MLIPs by augmenting the training data set with inexpensive Morse data along the dissociation pathways. Further, this physics-constrained data augmentation (PCDA) approach results in MLIPs with smooth bond dissociation curves as well as near coupled-cluster level BDEs, all without requiring any expensive multireference quantum mechanical calculations. A case study for methane combustion demonstrates how the PCDA approach can improve an existing reactive MLIP, namely, ANI-1xnr. In conclusion, not only are the BDEs and bond dissociation curves for all radicals and molecules significantly improved compared to ANI-1xnr but the PCDA-trained MLIP retains the reliability of ANI-1xnr when performing reactive molecular dynamics simulations.

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

DONKEY: A Flexible and Accurate Algorithm for Clustering

We propose an accurate clustering algorithm suitable for the varied and multidimensional data sets that correspond to temporal snapshots from on-the-fly nonadiabatic trajectory-based simulations of photoexcited dynamics. The algorithm approximates the underlying probability density function using variable kernel density estimation, with local maxima corresponding to cluster centers. Each data point is then assigned to one of the maxima by employing a maximization procedure. Finally, clusters artificially separated by minor fluctuations in the probability density are merged. The algorithm does not require parameter tuning, which ensures flexibility and reduces the risk of bias. It is tested on several synthetic data sets, where it consistently outperforms conventional clustering algorithms. As a final example, the algorithm is applied to the excited dynamics of the norbornadiene ⇌ quadricyclane (C 7 H 8 ) molecular photoswitch, demonstrating how distinct reaction pathways can be identified.

algorithms↗
DOE OSTI2025

LeaPP: Learning Pathways to Polymorphs through Machine Learning Analysis of Atomic Trajectories

Understanding the mechanisms underlying crystal nucleation and growth is crucial for many technological applications. Due to the short length and time scales involved, crystal nucleation is often studied using molecular simulations. Most existing approaches to extract the nucleation mechanism from simulations focus on the analysis of static snapshots of the configurations, potentially overlooking subtle local fluctuations and the history of the particles involved in the formation of solid nuclei. Here, in this work, we propose a novel methodology called LeaPP that categorizes nucleation trajectories based on the temporal information of their constituent particles. We leverage the time evolution of the local environment of the crystallizing particles to encapsulate the relationship between the structure and dynamics and distinguish between different evolving particle paths. Identification of the distinct particle paths further enables characterizing the nucleation trajectories into different pathways. Collectively, LeaPP provides a more nuanced understanding of nucleation through an unsupervised approach with lesser dependence on traditional order parameters. Furthermore, the pathways identified by LeaPP are predictive of the resulting polymorph. We demonstrate LeaPP on three different systems─Lennard-Jones-like particles, Ni 3 Al, and water on surfaces. The general methodology underlying LeaPP─considering the time evolution of the building blocks─applies to a wide range of self-assembly problems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Simulating Crystallization in a Colloidal System Using State Predictive Information Bottleneck Based Enhanced Sampling

Here, we investigate crystal nucleation in supersaturated colloid suspensions using enhanced molecular dynamics simulations augmented with machine learning techniques. The simulations reveal that crystallization in the model colloidal system studied here, with particles interacting through a repulsive screened Coulomb Yukawa potential, proceeds from vapor to dense liquid droplet to crystalline phases across multiple high barriers. Employing a one-dimensional reaction coordinate derived from the State Predictive Information Bottleneck framework, our simulations capture back-and-forth phase transitions across multiple barriers effectively in biased metadynamics simulations. We obtain relative free energy differences between different phases and also quantify the roles of different molecular level features in driving the phase changes.

Chemistry↗
DOE OSTI2025

Crystal Nucleation Kinetics and Mechanism: Influence of Interaction Potential

Modulating liquid-to-solid transitions and the resulting crystalline structure for tailored properties is much desired. Colloidal systems are exemplary to this end, but the fundamental knowledge gaps in relating the influence of intermolecular interactions to crystallization behavior continue to hinder progress. In this study, we address this knowledge gap by studying nucleation and growth in systems with modified Lennard-Jones potential. Specifically, we study the commonly used 12-6 potential and a softer 7-6 potential. The thermodynamic state point for the study is chosen such that both systems are investigated at the same level of supercooling and pressure. Under these conditions, we find that the nucleation rate for both systems is comparable. Interestingly, the nucleation pathways and resulting crystal structures are different. In the 12-6 system, nucleation and growth occur predominantly through the FCC structure. Softening the potential alters the critical nucleus composition and introduces two distinct nucleation pathways. One pathway predominantly leads to the nucleus with a body-centered cubic (BCC) structure, while the other favors the face-centered cubic (FCC) arrangement. Our study illustrates that polymorph selection can be achieved through modifications to intermolecular interactions without impacting nucleation kinetics. The results have significant implications in designing approaches for polymorph selection and modulating self-assembly mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2024

UV-Induced Reaction Pathways in Bromoform Probed with Ultrafast Electron Diffraction

For many chemical reactions, it remains notoriously difficult to predict and experimentally determine the rates and branching ratios between different reaction channels. This is particularly the case for reactions involving short-lived intermediates, whose observation requires ultrafast methods. The UV photochemistry of bromoform (CHBr 3 ) is among the most intensely studied photoreactions. Yet, a detailed understanding of the chemical pathways leading to the production of atomic Br and molecular Br 2 fragments has proven challenging. In particular, the role of isomerization and/or roaming and their competition with direct C–Br bond scission has been a matter of continued debate. Here, in this work, gas-phase ultrafast megaelectronvolt electron diffraction (MeV-UED) is used to directly study structural dynamics in bromoform after single 267 nm photon excitation with femtosecond temporal resolution. The results show unambiguously that isomerization contributes significantly to the early stages of the UV photochemistry of bromoform. In addition to direct C–Br bond breaking within <200 fs, formation of iso-CHBr 3 (Br-CH-Br-Br) is observed on the same time scale and with an isomer lifetime of >1.1 ps. The branching ratio between direct dissociation and isomerization is determined to be 0.4 ± 0.2:0.6 ± 0.2, i.e., approximately 60% of molecules undergo isomerization within the first few hundred femtoseconds after UV excitation. The structure and time of formation of iso-CHBr 3 compare favorably with the results of an ab initio molecular dynamics simulation. The lifetime and interatomic distances of the isomer are consistent with the involvement of a roaming reaction mechanism.

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

Cascade CO 2 Insertion in Carbanion Ionic Liquids Driven by Structure Rearrangement

The CO 2 chemisorption in state-of-the-art sorbents based on oxide/hydroxide/amine moieties is driven by strong chemical bonding formation in the carbonate/bicarbonate/carbamate products, which in turn leads to high energy input in sorbent regeneration. In addition, the CO 2 uptake capacity was limited by the active sites’ utilization efficiency, with each active site incorporating one CO 2 molecule or less. In this work, a new concept and generation of sorbent was developed to achieve cascade insertion of multiple CO 2 molecules by leveraging structure rearrangement as the driving force, leading to in situ generation of extra CO 2 -binding sites and significantly reduced energy input for CO 2 release. The designed ionic liquids (ILs) containing carbanions with conjugated and asymmetric structure, deprotonated (methylsulfonyl)acetonitrile ([MSA]) anion, allowed the cascade insertion of two CO 2 molecules via consecutive C–C and O–C bond formations. The proton transfer and structure rearrangement of the carboxylic acid intermediates played critical roles in stabilizing the first integrated CO 2 and generating extra electron-rich oxygen sites for the insertion of the second CO 2 . The structure variation and reaction pathway were confirmed by operando spectroscopy, magnetic resonance spectroscopy (NMR), mass spectroscopy, and computational chemistry. The energy input in sorbent regeneration could be further reduced by harnessing the phase-changing behavior of the carbanion salts in ether solutions upon reacting with CO 2 , avoiding the energy consumption in heating the solvent. In conclusion, the fundamental insights obtained herein provide a promising approach to greatly improve the CO 2 sorption performance via sophisticated molecular-scale structural engineering of the sorbents.

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

Three rate-determining protein roles in photosynthetic O 2 -evolution addressed by time-resolved experiments on genetically modified photosystems

Light-driven water splitting by plants, algae and cyanobacteria is pivotal for global bioenergetics and biomass formation. A manganese cluster bound to the photosystem II proteins catalyzes the complex reaction at high rate, but the rate-determining factors are insufficiently understood. Here we trace the oxygen-evolution transition by time-resolved polarography and infrared spectroscopy for cyanobacterial photosystems genetically modified at two strategic sites, complemented by computational chemistry. Our results highlight three rate-determining roles of the protein environment of the metal cluster: acceleration of proton-coupled electron transfer, acceleration of substrate-water insertion after O 2 -formation, and balancing of rate-determining enthalpic and entropic contributions. Whereas in general the substrate-water insertion step may be unresolvable in time-resolved experiments, here it likely becomes traceable because of deceleration by genetic modification. Our results may stimulate new time-resolved experiments on substrate-water insertion in photosynthesis, clarification of enthalpy-entropy compensation in enzyme catalysis, and knowledge-guided development of inorganic catalyst materials.

Bioenergetics↗
DOE OSTI2025

Accurate Dehydrogenation Enthalpies Dataset for Liquid Organic Hydrogen Carriers

This contribution presents a comprehensive extension of the QM9 dataset (originally at 133 K molecules) with the calculation of G4MP2 enthalpies for 9,841 molecules, featuring up to nine heavy atoms. We present QM9-LOHC, a (de)hydrogenation dataset of 10,373 reactions, including a minimum of 5.5% weight hydrogen storage capacity in line with the Department of Energy standards for Liquid Organic Hydrogen Carriers (LOHC). By utilizing the accurate quantum chemical method G4MP2 we expand the QM9 database and explore new avenues for the exploration of hydrogen storage technologies (electrochemical LOHCs, alkali metal-LOHCs, and mixtures of LOHCs). The QM9-LOHC dataset, with its focus on reactions that vary only by hydrogen saturation levels, provides a needed data resource for advancing the design and optimization of both conventional and innovative LOHC systems, and high-fidelity data for molecular discovery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗