New force theorem.
New force theorem for determining molecular energy derivatives with respect to internuclear distance derived from Hellmann-Feynman expression
SEARCH · Engineering Papers
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.
New force theorem for determining molecular energy derivatives with respect to internuclear distance derived from Hellmann-Feynman expression
High energy and thermal energy molecular beams elastic scattering for intermolecular potentials determination
Molecular vibration absorption spectra from variation of field emitted electron energy distribution due to inelastic electron-molecule interaction at metal-vacuum interface
Electron energy-loss spectra in molecular fluorine, for energy losses from 0 to 17.0 eV, have been taken at incident electron energies of 30, 50, and 90 eV and scattering angles from 5 to 140 deg. Features in the spectra above 11.5 eV energy loss agree well with the assignments recently made from optical spectroscopy. Excitations of many of the eleven repulsive valence excited electronic states are observed and their location correlates reasonably well with recent theoretical results. Several of these excitations have been observed for the first time and four features, for which there are no identifications, appear in the spectra.
Not Available
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.
Molecular vibration spectra from inelastic interaction between electrons and absorbed molecules at metal-vacuum interface
Shock wave thickness predictions from Monte Carlo Direct Simulations, using differential scattering and the Maitland-Smith-Aziz interatomic potential, underpredict experiments as shock Mach numbers increase above about 4. Examination of several sources of data has indicated that at relatively high energies the repulsive portion of accepted potentials such as the Maitland-Smith-Aziz may be too steep. An Exponential-6 potential due to Ross, based on high energy molecular beam scattering data and shock velocity measurements in liquid argon, has been combined with the lower energy portion of the Maitland-Smith-Aziz potential. When this hybrid potential is used in Monte Carlo Direct Simulations, agreement with experiments is improved over the previous predictions using the pure Maitland-Smith-Aziz form.
The molecular constants and their adiabatic corrections have been determined for the (A 1 Sigma +) - (X 1 Sigma +) system of the isotopic lithium hydrides: (Li-6)H, (Li-7)H, (Li-6)D, and (Li-7)D. Using a fully quantum mechanical variational method, the potential energy curves (IPA potentials) are determined. Extending the variational method, we have obtained for the first time adiabatic corrections of potential energy curves from isotopic spectroscopic data. A significant difference between the potential energy curves of the lithium hydrides and the lithium deuterides has been observed. When Li-6 was replaced by Li-7, a significant difference was only observed for the (A 1 Sigma +) state, but not for the (X 1 Sigma +) state.
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.
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.
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.
Elastic scattering cross section of low energy electrons calculated by continuum Hartree-Fock equations
The purpose of this workshop was to define new lasants and transfer gases for blackbody pumped lasers. The goal was to find gases with the correct energies and lifetimes for use as lasants and/or transfer gases. The more immediate goal of the participants was to suggest molecules, currently available, (or to be synthesized) which either satisfy criteria which was developed or which deserve further study because from limited information they seem to satisfy the criteria. Thus, the workshop was interested in the identification of molecules and the rational for suggesting them. There were two basic kinds of blackbody lasers to be discussed, the cavity laser and the transfer laser. The physics of optical absorption and laser emission was addressed. The physics of the transfer laser was discussed. Molecular energy transfer from one molecule to another was added to the concepts already involved in the cavity laser. The physics of the transfer gas and the transfer-lasant gas combinations was also discussed.
Electric energy behavior for ground state of one-electron heteronuclear diatomic molecule near united atom
The rate of heating which occurs in the middle atmosphere due to four exothermic reactions involving members of the odd-hydrogen family is calculated. The following reactions are considered: O + OH yields O2 + H; H + O2 + M yields HO2 + M; H + O3 yields OH + O2; and O + HO2 yields OH + O2. It is shown that the heating rates due to these reactions rival the oxygen-related heating rates conventionally considered in middle-atmosphere models. The conversion of chemical potential energy into molecular translational energy (heat) by these odd-hydrogen reactions is shown to be a significant energy source in the middle atmosphere that has not been previously considered.
The structure of a strong plane shock wave in a monatomic rarefied perfect gas is one of the simplest problems able to be posed in kinetic theory, and one of the hardest to solve. Its simplicity lies in the absence of solid boundaries, geometrical complications, or internal molecular energy. Its difficulty arises from the great departure of the gas from equilibrium within the shock, which invalidates many of the techniques used successfully elsewhere in kinetic theory. In addition to this theoretical challenge, the modern development of ballistics and hypersonic flight has helped to stimulate extensive theoretical and experimental interest in the shock problem. The experimenters in turn have encountered great difficulties on account of the very small physical dimensions of shocks. In fact, until very recently indeed, any close comparisons of theoretical and experimental shock structure results have been rather unprofitable due to the inadequacies of both theory and experiment. During the last few years this situation has been appreciably improved by development of the Monte Carlo method. This allows idealized 'experiments' to be performed on large computers instead of in wind tunnels, using a known intermolecular force law. The most developed of these methods has been shown to be equivalent theoretically to the Boltzmann equation and to give results which agree extremely closely with measurements of high accuracy. Thus Monte Carlo results not only form the soundest basis for our present theoretical knowledge of shock wave structure, but, for purposes of developing other theories, can also be considered a very valuable experimental resource. However, such results remain very expensive to obtain. In this thesis we develop more economical kinetic theory methods for the approximate prediction of shock structure, and compare our results with those of the Monte Carlo method.
Observations made from Kitt Peak, AZ, of the 233.9 GHz emission line of O-18O-16 in the upper stratosphere and lower mesosphere are reported. A good model fit to the line is obtained by incorporation of pressure and temperature broadening effects, as well as a Zeeman splitting algorithm that uses a standard geomagnetic field model and a paramagnetic Hamiltonian description of the molecular energy states. These observations are used, along with the well known O-18O-16 mixing ratio, to establish absolute calibration for observations of other chemical species from Kitt Peak. Repeated measurements show no change in this absolute calibration between observation dates. The wide magnetic splitting (+/-1.8 MHz) exhibited by this line with only six Zeeman components provides a unique test of middle atmosphere Zeeman effect model calculations, supporting the use Of 02 lines by microwave atmospheric sounders to measure pressure and temperature.