Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MOLECULAR ENERGY”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

How Accurate Are Simulations and Experiments for the Lattice Energies of Molecular Crystals?

Molecular crystals play a central role in a wide range of scientific fields, including pharmaceuticals and organic semiconductor devices. However, they are challenging systems to model accurately with computational approaches because of a delicate interplay of intermolecular interactions such as hydrogen bonding and Van der Waals dispersion forces. Here, by exploiting recent algorithmic developments, we report the first set of diffusion Monte Carlo lattice energies for all 23 molecular crystals in the popular and widely used X23 dataset. Comparisons with previous state-of-the-art lattice energy predictions (on a subset of the dataset) and a careful analysis of experimental sublimation enthalpies reveals that high-accuracy computational methods are now at least as reliable as (computationally derived) experiments for the lattice energies of molecular crystals. Overall, this work demonstrates the feasibility of high-level explicitly correlated electronic structure methods for broad benchmarking studies in complex condensed phase systems, and signposts a route towards closer agreement between experiment and simulation. Published by the American Physical Society 2024

Physics↗
DOE OSTI2026

Ground and excited state gradients with end-to-end differentiable semiempirical quantum chemistry

Accurate and efficient gradients of molecular energy with respect to nuclear degrees of freedom are essential for geometry optimization and molecular dynamics, including simulations that go beyond the Born–Oppenheimer regime. A common approach involves deriving analytical formulas for new electronic structure methods, which is often conceptually difficult and requires tedious coding. Here, we implement analytical, semi-numerical, and automatic differentiation (AD)-based gradient pathways for semiempirical Hamiltonian models in the PYSEQM software package, leveraging both graphics processing unit (GPU) and central processing unit (CPU) architectures. We further extend these capabilities to excited states calculated using the configuration interaction singles and time-dependent Hartree–Fock ansätze. We benchmark wall time, peak memory usage, and accuracy across three molecular families of varying chemical complexity, including systems of up to a thousand atoms. For ground-state simulations, analytical and AD gradients achieve near-identical GPU runtimes, while semi-numerical gradients are slower on GPU but remain competitive on CPU. For excited states, both analytical and custom AD approaches using implicit differentiation show similar performance and low memory requirements, whereas gradients with full AD are memory-limited. AD gradients match analytical ones in accuracy across all tested systems, aided by a quaternion-based diatomic frame rotation for two-center quantities that ensures smooth energy surfaces. Overall, automatic differentiation emerges as a practical alternative to analytical gradients in semiempirical quantum chemistry, offering high accuracy while allowing seamless integration in AI-driven workflows and popular packages, such as PyTorch and JAX. Our results provide actionable guidance for selecting optimal gradient strategies in large-scale ground- and excited-state molecular dynamics simulations.

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

Strong adsorption of guanidinium cations to the air–water interface

Combining Deep-UV second harmonic generation spectroscopy with molecular simulations, we confirm and quantify the specific adsorption of guanidinium cations to the air–water interface. Using a Langmuir analysis of measurements at multiple concentrations, we extract the Gibbs free energy of adsorption, finding it larger than typical thermal energies. Molecular simulations clarify the role of polarizability in tuning the thermodynamics of adsorption, and establish the preferential parallel alignment of guanidinium at the air–water interface. As a polyatomic cation, guanidinium represents one of the few examples of a positively charged species to exhibit a propensity for the air-water interface. As such, these results expand on the growing body of work on specific ion adsorption.

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

Including Physics-Informed Atomization Constraints in Neural Networks for Reactive Chemistry

Machine learning interatomic potentials (MLIPs) have emerged as powerful tools for investigating atomistic systems with high accuracy and a relatively low computational cost. However, a common and unaddressed challenge with many current neural network (NN) MLIP models is their limited ability to accurately predict the relative energies of systems containing isolated or nearly isolated atoms, which appear in various reactive processes. To address this limitation, we present a mathematical technique for modifying any existing atom-centered NN architecture to account for the energies of isolated atoms. The result produces a consistent prediction of the atomization energy (AE) of a system using minimal constraints on the model. Using this technique, we build a model architecture that we call hierarchically interacting particle neural network (HIP-NN)-AE, an AE-constrained version of the HIP-NN, as well as ANI-AE, the AE-constrained version of the accurate NN engine for molecular energies (ANI). Our results demonstrate AE consistency of AE-constrained models, which drastically improves the AE predictions for the models. We compare the AE-constrained approach to unconstrained models as well as models from the literature in other scenarios, such as bond dissociation energies, bond dissociation pathways, and extensibility tests. These results show that the constraints improve the model performance in some of these tasks and do not negatively affect the performance on any tasks. The AE constraint approach thus offers a robust solution to the challenges posed by isolated atoms in energy prediction tasks.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Strong NIR II Chiroptical Response and Magnetic Anisotropy via Modular Installation of Chiral Capping Ligands on a Light-Emitting Diradicaloid Scaffold

Low-energy molecular lumiphores have seen increased interest due to potential imaging and communications applications. Specifically, molecules that emit in the near-infrared (NIR, 700–1700 nm) or telecom (∼1260–1625 nm) regions, where attenuation is minimized in biological tissue and optical fibers, respectively, can drastically improve image resolution and depth penetration; however, bright low-energy emission is rare due to exponentially decreasing quantum yields in this region. Chiral molecules exhibiting strong NIR or telecom absorption/emission would be of particular interest due to advanced security and spintronics applications, but these compounds remain scarce and are currently restricted to lanthanide or nanoparticle-based systems. Here, we report the synthesis of a chiral organic-based NIR emitter, enabled by simple peripheral installation of chiral capping ligands. These ligands twist the achiral NIR-emissive core, breaking inversion symmetry. Furthermore, this twisting enables strong chiroptical responses observed via static and transient circular dichroism and increased magnetic anisotropy through electron paramagnetic resonance (EPR) measurements, establishing this strategy as a promising method for the development of new chiral emitters and sensors.

Electron paramagnetic resonance spectroscopy↗
DOE OSTI2024

Informing solar blind radioluminescence imaging through a calibrated spectrum

While direct radiation detection methods offer great insight into the origin of ionizing particles and photons, their use to locate contaminated areas or concealed radioactive sources can lead to undue exposure of personnel and equipment to ionizing radiation or the potential for contamination. These same sources induce ultraviolet (UV) optical photon fluorescence in air – a process referred to as radioluminescence – that may be imaged from low dose regions over larger attenuation lengths than ionizing radiation. However, most optical detection methods are limited to low lighting conditions to image the more abundant ultraviolet-A (UV-A) photons. To extend this capability to room light or daytime conditions, the solar blind region (Ultraviolet-C (UV-C), <280 nm) can be tapped. Though the emission yield of UV-C photons is roughly two orders of magnitude lower than that in the UV-A regime, the UV-C offers dramatic improvements in signal-to-noise ratios under bright lighting conditions due to decreased background interferences. The yield of specific UV-C lines, if present in the literature at all, varies widely, which has a large impact in modeling and analyzing standoff UV-C measurement scenarios. Thus, we have captured improved radioluminescence spectra over 250–400 nm and identified observed emission peaks in ambient air. Many of the UV-C photons produced by ionizing radiation excitation result from high-energy, molecular nitrogen Gaydon-Herman transitions which have had limited study to date for this application. Relating these findings to published UV-A yields, we estimate emissions between 0.15–0.19 photons/MeV over 250–280 nm. Additionally, we also use a commercial corona-discharge imaging camera to demonstrate outdoor UV-C radiation mapping of alpha and gamma emitters from 50 and 75 m standoffs, respectively. The imaged “counts” are compared to optically modelled values and show the same trend over distance.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Suppressed Cosmic-Ray Energy Densities in Molecular Clouds from Streaming Instability-regulated Transport

Abstract Cosmic rays (CRs) are the primary driver of ionization in star-forming molecular clouds (MCs). Despite their potential impacts on gas dynamics and chemistry, no simulations of star cluster formation following the creation of individual stars have included explicit cosmic-ray transport (CRT) to date. We conduct the first numerical simulations following the collapse of a 2000 M ⊙ MC and the subsequent star formation including CRT using the STAR FORmation in Gaseous Environments framework implemented in the GIZMO code. We show that when CRT is streaming-dominated, the CR energy in the cloud is strongly attenuated due to energy losses from the streaming instability. Consequently, in a Milky Way–like environment the median CR ionization rate in the cloud is low ( ζ ≲ 2 × 10 −19 s −1 ) during the main star-forming epoch of the calculation and the impact of CRs on the star formation in the cloud is limited. However, in high-CR environments, the CR distribution in the cloud is elevated ( ζ ≲ 6 × 10 −18 ), and the relatively higher CR pressure outside the cloud causes slightly earlier cloud collapse and increases the star formation efficiency by 50% to ∼13%. The initial mass function is similar in all cases except with possible variations in a high-CR environment. Further studies are needed to explain the range of ionization rates observed in MCs and explore star formation in extreme CR environments.

Astronomy & Astrophysics↗
DOE OSTI2025

Tetracene Functionalized Si(111) Achieves Enhanced Solar-to-Chemical Energy Conversion via Molecular Acceptor States

The properties of semiconductor|liquid interfaces play a critical role in determining the efficiency of solar-to-hydrogen (STH) conversion. Here, we investigate how molecular functionalization of Si(111) and Si(111)|TiO 2 surfaces impacts photoelectrochemical (PEC) hydrogen production efficiency. We find that functionalization of ∼3% of the atop sites of Si(111) with either 9-anthracene (Anth) or 5-tetracene (Tet), with the remaining sites passivated by methyl groups, provides substrates with high electronic quality and low surface oxide densities, as determined by X-ray photoelectron spectroscopy (XPS) measurements. Surface photovoltage (SPV) spectroscopy shows that surfaces modified with Anth or Tet exhibit an increased photovoltage, with Tet-functionalized surfaces yielding an additional 192 meV relative to methyl-terminated Si(111), indicating improved charge separation for Si-Tet. Further improvement in onset potential was achieved by replacing a nitrogen-containing TiO 2 atomic layer deposition (ALD) precursor (TDMAT) with a precursor lacking nitrogen (TTIP), which eliminates the parasitic defect band in the TiO 2 overlayer (p-Si(111)-Tet|TTIP-TiO 2 |Pt: V OC = +0.283 ± 0.041 V vs RHE). Density functional theory (DFT) analysis demonstrates that compared with Anth-modified Si(111), the Tet-modified surface exhibits more hybridized Si(111)-Tet states closer to the silicon band edges. Mercury contact current–voltage (I–V, dark) measurements quantified the relative interfacial density of states of Si-Tet, Si-Anth and Si-Me surfaces─revealing that the interfacial state density was highest for Si-Tet. This suggests that such hybridized interfaces serve to capture better photoexcited charge, which enables facile electron transfer to molecular acceptors in solution. Altogether, the data indicate that beneficial hybrid molecular LUMO surface states interacting with the Si conduction band edge results in improved hydrogen evolution (HER) performance for p-Si devices.

Group theory↗
DOE OSTI2024

Energy conversion and transport in molecular-scale junctions

Molecular-scale junctions (MSJs) have been considered the ideal testbed for probing physical and chemical processes at the molecular scale. Due to nanometric confinement, charge and energy transport in MSJs are governed by quantum mechanically dictated energy profiles, which can be tuned chemically or physically with atomic precision, offering rich possibilities beyond conventional semiconductor devices. While charge transport in MSJs has been extensively studied over the past two decades, understanding energy conversion and transport in MSJs has only become experimentally attainable in recent years. As demonstrated recently, by tuning the quantum interplay between the electrodes, the molecular core, and the contact interfaces, energy processes can be manipulated to achieve desired functionalities, opening new avenues for molecular electronics, energy harvesting, and sensing applications. This Review provides a comprehensive overview and critical analysis of various forms of energy conversion and transport processes in MSJs and their associated applications. We elaborate on energy-related processes mediated by the interaction between the core molecular structure in MSJs and different external stimuli, such as light, heat, electric field, magnetic field, force, and other environmental cues. Key topics covered include photovoltaics, electroluminescence, thermoelectricity, heat conduction, catalysis, spin-mediated phenomena, and vibrational effects. Furthermore, the review concludes with a discussion of existing challenges and future opportunities, aiming to facilitate in-depth future investigation of promising experimental platforms, molecular design principles, control strategies, and new application scenarios.

Charge transport↗
DOE OSTI2025

Energy-efficient, Large-scale Molecular Dynamics Simulations via Hardware- and Algorithm-level Optimization

This work aims to develop a framework for energy-efficient computing that will enable molecular dynamics (MD) simulations of large-scale phenomena with atomic precision and simultaneously remove computational bottlenecks limiting the speed of MD simulations. We seek to implement such an approach through the development of surrogate models for the interatomic force calculation combined with the use of mixed numerical precision formats. For a model system of neutral atoms (only pairwise interactions), significant force calculation efficiency improvements were achieved, without detrimental effects on atomic structures or average energies, using single precision, by developing a surrogate model (deep neural network), and by quantizing this surrogate model. For a model system of charged atoms, the reciprocal-space calculation of electrostatic interactions was identified as the main bottleneck, and the development of a surrogate model should be pursued to achieve an estimated one-order-of-magnitude additional speedup.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Neuromorphic heat transport effects in a molecular junction

Understanding energy transport at the nanoscale is an open and fundamental challenge in the molecular sciences with direct implications for the design of new electronics, computing devices, and materials. While nanoscale energy transport under steady-state conditions has been studied extensively, there is much less known about energy transport under time-dependent driving forces, particularly in the far-from-equilibrium regime. In this work, we use nonequilibrium molecular dynamics simulations and stochastic thermodynamics to investigate energy transport in a well-studied nanoscale system—a molecular junction—subjected to a time-periodic temperature gradient. The primary observation is that molecular junctions can exhibit heat transport hysteresis, a phenomenon in which the heat flux through a system depends not only on the instantaneous value of a time-dependent temperature bias but also on the temporal history of that bias. The presented findings illustrate that molecular junctions can exhibit the specific memory effect—heat transport hysteresis—that is essential for the design of thermal neuromorphic computers. This work elucidates a potential pathway toward the realization of such devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Hybrid Quantum Mechanical, Molecular Mechanical, and Machine Learning Potential for Computing Aqueous-Phase Adsorption Free Energies on Metal Surfaces

Performing reliable computer simulations of elementary processes occurring at metal–water interfaces is pivotal for novel catalyst design in sustainable energy applications. Computational catalyst design hinges on the ability to reliably and efficiently compute the potential energy surface (PES) of the system. Here, due to the large system sizes needed for studying processes at liquid water–metal interfaces, these systems can currently not be described using density functional theory (DFT). In this work, we used a hybrid quantum mechanical, molecular mechanical, and machine learning potential for studying the adsorption behavior of phenol, atomic hydrogen, 2-butanol, and 2-butanone on the (0001) facet of Ru under reducing conditions when Ru is not oxidized. Specifically, we describe the adsorbate and the surrounding metal atoms at the DFT level of theory. Here, we also considered the electrostatic field effect of the water molecules on adsorbate–metal interactions. Next, for the water–water and water–adsorbate interactions, we used established classical force fields. Finally, for the water–Ru surface interaction, for which no reliable force fields have been published, we used Behler–Parrinello high-dimensional neural network potentials (HDNNPs). Employing this setup, we used our explicit solvation for metal surface (eSMS) approach to compute the aqueous-phase effect on the low-coverage adsorption of selected molecules and atoms on the (0001) facet of Ru. In agreement with previous experimental and computational studies of oxygenated molecules over transition metal facets, we found that liquid water destabilizes the tested adsorbates on Ru(0001). Interestingly, our findings indicate that adsorbates on Ru are less affected by the presence of an aqueous phase than on other transition metals (e.g., Pt), highlighting the necessity of experimental investigations of Ru-based catalytic systems in liquid water.

Adsorption↗
DOE OSTI2025

Stimulus-Responsive Modulation of Solvation Environments in Solid Catalysts

Liquid environments play a crucial role in the biological processes occurring in living organisms as well as in many human-made processes involving electrochemistry, photo-, and thermocatalysis. In the majority of these systems, aqueous phases are ubiquitous due to water’s natural abundance. Water molecules, however, can exert large changes in the chemical environment of catalytically active sites, altering the reaction rates, selectivity, and catalyst stability. These solvation effects induced by water molecules near catalytic sites can drastically change the energy landscape and unlock unique reaction pathways with far more favorable kinetics. In nature, living organisms couple these complex interactions with detection, communication, and actuation mechanisms to induce self-regulatory behavior, ensuring stability of the system and thus long-term durability. Extrapolating this behavior to heterogeneous catalysis is desirable because the resulting “smart materials” can potentially unlock new chemical conversion processes with higher atom efficiency, rates, and stability. The combination of polymer chemistry and heterogeneous catalysis has introduced versatile approaches for creating materials that can respond to cues in the reaction medium that alter the accessibility, intrinsic activity, and selectivity of the catalyst. To achieve this, one could combine stimulus-responsive polymers, which undergo a large volumetric phase transition in response to an external stimulus, with a solid catalyst. This chemo-mechanical response has been employed to create a variety of nanoreactor vessels with stimulus-responsive character that turn on- and off- depending on the reaction conditions. In this Account, we focus on the impact of these polymer coatings on the solvation environment around the active site and the implications of these effects on the reaction energy landscape, molecular arrangement of the solvent, electric fields at the catalyst–liquid interface, binding energy, and mobility of surface reaction intermediates. These seemingly subtle changes in solvent molecules induced by the presence of polymers can have a tremendous impact on the development of bioinspired heterogeneous catalysts, reliable chemical clocks, micro/nanoreactors, and robots. The large library of polymer chemistries offers a plethora of combinations of stimulus-responsive mechanisms (e.g., temperature, pH, light, magnetic field, solvent composition), providing the possibility of creating homeostatic catalysts à la carte.

catalysts↗
DOE OSTI2025

Machine learning-accelerated path integral molecular dynamics simulations of reactive organic electrolytes

Hydrogen bonded electrolytes that exhibit accelerated proton transport via sequential reactive hops have drawn interest for their promise in clean energy applications. Molecular dynamics simulations of these electrolytes offer the opportunity to uncover microscopic mechanistic details that could be used to design and tune the properties of candidate electrolyte technologies. However, accurately modeling the proton transfer reactions and transport properties that give rise to high charge conductivites in these electrolytes proves computationally challenging because of the need to perform lengthy condensed phase simulations, treating both the electronic and nuclear degrees of freedom quantum mechanically. In this paper, we demonstrate that such a modeling task can be efficiently achieved with the use of density functional theory (DFT)-trained machine learning potentials (MLP) to accelerate path integral molecular dynamics (PIMD) simulations. We highlight the practical utility of this approach by using it to benchmark how closely PIMD simulations employing different DFT exchange–correlation functionals reproduce the composition-dependent densities, diffusion coefficients, and electrical conductivities of mixtures consisting of imidazole and levulinic acid. Even with the speedup afforded by our MLPs, PIMD simulations remain quite expensive. Furthermore, in order to render PIMD more computationally tractable, we introduce and benchmark the accuracy of a ring polymer contraction approach that leverages a computationally efficient short-range MLP to accelerate our PIMD simulations by an additional factor of four.

Chemical bonding↗
DOE OSTI2025

Methods in PES-Learn: Direct-Fit Machine Learning of Born–Oppenheimer Potential Energy Surfaces

The release of PES-L EARN version 1.0 as an open-source software package for the automatic construction of machine learning models of semi-global molecular potential energy surfaces (PESs) is presented. Improvements to PES-L EARN ’s interoperability are stressed with new Python API that simplifies workflows for PES construction via interaction with QCSchema input and output infrastructure. In addition, a new machine learning method is introduced to PES-L EARN : kernel ridge regression (KRR). The capabilities of KRR are emphasized with examination of select semi-global PESs. All machine learning methods available in PES-L EARN are benchmarked with benzene and ethanol datasets from the rMD17 database to illustrate PES-L EARN ’s performance ability. Fitting performance and timings are assessed for both systems. Finally, the ability to predict gradients with neural network models is presented and benchmarked with ethanol and benzene. PES-L EARN is an active project and welcomes community suggestions and contributions.

kernel ridge regression↗
DOE OSTI2024

Controlling Interfacial Charge Separation Energetics and Kinetics

The overall goal is to probe electron transfer kinetics of the Marcus inverted region at interfaces. This project specifically uses a dye‐semiconductor interface with a reversible electron transfer reagent in solution where the semiconductor conduction band and electron transfer reagent ideally have a potential energy difference firmly within the Marcus inverted region. The goals of this project are: (Goal 1) Identify dye design elements that give dyes with low energy ground‐state oxidation potentials while retaining desirable electron transfer kinetics at the semiconductor interface. (Goal 2) Identify dye designs with excited‐state oxidation potentials set to fixed energies through molecular design while the dye ground‐state oxidation potentials are varied to allow for the evaluation of varied electron transfer reagents/redox shuttles (RS) with different free energies for electron transfer. (Goal 3) Within the context of the functionality identified in goals 1 and 2, goal 3 seeks to design dyes with features maximizing charge carrier lifetimes by minimizing recombination losses at metal oxide interfaces with high voltage systems. (Goal 4) Probe electron transfer kinetics in a high voltage‐dye sensitized solar cell with both the oxidized dye and RS firmly within the Marcus inverted region with regard to the TiO2 CB. (Goal 5) Identify tailored sensitizer/redox shuttle pairs with desirable electron transfer kinetics in defined spectral regions capable of “banding off” incident solar irradiation to give high voltage multijunction systems capable of exceeding the Shockley‐Queisser limit and powering catalytic processes.

14 SOLAR ENERGY↗