How to Probe Hydrated Dielectrons Experimentally: Ab Initio Simulations of the Absorption Spectra of Aqueous Dielectrons, Electron Pairs, and Hydride
Not Available
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
SEARCH · Engineering Papers
Results for “quantum mechanical calculations”
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.
A Change in C–H Activation Mechanism: Experimental and Computational Investigations of Rh-Catalyzed Disubstituted Benzene Functionalization
We report the ethenylation of 1,3- and 1,2-disubstituted benzenes using [(η 2 -C 2 H 4 ) 2 Rh(μ-OAc)] 2 as a catalyst precursor and Cu(OPiv) 2 as the oxidant. The regioselectivity of alkenylation for 1,3-disubstituted benzenes produces 3,5-disubstituted styrene products, while the alkenylation of 1,2-disubstituted benzenes produces 3,4- disubstituted styrene products. The rate of alkenylation is influenced by steric and electronic factors based on the substituents of the benzene unit. In all cases, 1,2-disubstituted benzenes react faster than 1,3-disubstituted benzenes, with a rate difference that is from 2 times up to >70 times more rapid for 1,2-disubstituted substrates. This is likely due to the difference in the number of accessible C−H bonds based on the steric protection of C−H bonds adjacent to functionality. Furthermore, the rate of alkenylation is influenced by the arene substituent electronics. The rates of alkenylation for 1,2-disubstituted benzenes follow the trend OMe > Me > CF 3 > Cl, while for 1,3-disubstituted benzenes the trend is CF 3 > Cl > Me > OMe. Using quantum mechanics DFT calculations, we found that the C−H activation step can occur by two different mechanisms. The electronic properties of substituents on the arene ring change the preferred C−H bond-breaking mechanism for 1,2-disubstituted and 1,3-disubstituted benzenes.
Aromatic compounds↗
Machine Learning Screening of Metal-Ion Battery Electrode Materials
Here, in this work we present deep neural network regression machine learning models (ML) for predicting the average voltage and the percentage change in volume of battery electrodes upon charging and discharging with metal ions. Our models exhibit good performance as measured by the average mean absolute error obtained from a 10-fold cross-validation as well as on independent test sets. We further assess the robustness our ML models by investigating their screening potential beyond the training database. We produce novel Na-ion electrodes by systematically replacing Li-ions in the original database by Na-ions, and then selecting a set of 22 electrodes that exhibit a good performance in energy density as well as small volume variations upon charging and discharging, as predicted by the machine learning model. The ML predictions for these new materials are then compared to quantum-mechanics based calculations. Our results reaffirm the significant role of machine learning techniques in the exploration of materials for battery applications.
,electrode volume change↗
Pretreatment with Sodium Methyl Mercaptide Increases Carbohydrate Yield during Kraft Pulping
Kraft pulping is the predominant technology in the pulp and paper industry for removing lignin from wood carbohydrates to produce paper, board, packaging, tissue, and specialty cellulose. However, the kraft process is energy intensive and expensive, and its yield is limited by the degradation of carbohydrates. Pretreatment can increase carbohydrate yield by limiting degradation via primary peeling of reducing end groups. However, protection of galactoglucomannan (GGM), the primary hemicellulose component of softwood, is minimal when conventional pretreatments are used. Here, we investigated the effectiveness of sodium methyl mercaptide pretreatment on southern pine wood chips under a range of experimental conditions. We found that pretreatment of biomass with 4.38% sodium methyl mercaptide at pH 12 and 105 °C for 60 min provided small but significant increases in xylan and cellulose yields relative to control conditions, but preservation of GGM was minimal. Furthermore, to provide insight into molecular-scale details of primary peeling, pretreatment, and alkaline hydrolysis, we performed classical molecular dynamics (MD) simulations under selected process conditions and quantum mechanical (QM) calculations of selected reactions. MD simulations showed that C1 of the GGM reducing end is more readily accessible by HO – and CH 3 S – ions than in cellulose. The free energy barrier for peeling calculated with QM is lower for GGM than for cellulose, indicating increased susceptibility to peeling. In addition, we found that GGM may be more susceptible to internal chain cleavage than cellulose. Thus, even though reducing end groups may be protected initially through pretreatment, new unprotected reducing end groups may be generated through alkaline hydrolysis. Taken together, these findings show the promise of methyl mercaptide as a pretreatment technology for cellulose retention and also provide molecular insight for improving its effectiveness toward GGM.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Recombination of Autodissociated Water Ions in a Nanoscale Pure Water Droplet
The recombination of water ions has diverse scientific and practical implications, ranging from acid-base chemistry and biological systems to planetary environments and applications in fuel cell and carbon conversion technologies. While spatial confinement affects the physicochemical properties of water dynamics, its impact on the recombination process has rarely been studied. In this work, we investigate the dynamics of water, the water ion distribution, and the ion recombination process in water droplets as a function of droplet size through molecular dynamics simulations and adaptive quantum mechanical/molecular mechanical calculations. We compare the dynamics of recombination in water droplet sizes ranging from 100 to 18 000 waters, both in their interiors and on their surfaces. We found that the self-diffusion of water dramatically decreases in droplets with a diameter below 2.2 nm. Using a classical RexPoN force-field, we found that the ions in 1000 H2O's spend almost 50% of the time on the surface and 0.5 nm beneath it with a slight preference for OH- ion to reside longer on the surface. We estimate that, on average, recombination in these drops occurs at 400 ps in 1000 H2O's and 1 ns in 3000 H2O's. We also found that recombination is not limited by the local structure of the surface or the size of the droplet but can be influenced by the geometry of the water wire connecting the ions as they approach each other, which can often prevent recombination. Our results provide insights to the reaction microenvironments presented by nanoscopic water droplets.
Kwon, Soonho↗
Automated discovery of a robust interatomic potential for aluminum
Machine learning, trained on quantum mechanics (QM) calculations, is a powerful tool for modeling potential energy surfaces. A critical factor is the quality and diversity of the training dataset. Here we present a highly automated approach to dataset construction and demonstrate the method by building a potential for elemental aluminum (ANI-Al). In our active learning scheme, the ML potential under development is used to drive non-equilibrium molecular dynamics simulations with time-varying applied temperatures. Whenever a configuration is reached for which the ML uncertainty is large, new QM data is collected. The ML model is periodically retrained on all available QM data. The final ANI-Al potential makes very accurate predictions of radial distribution function in melt, liquid-solid coexistence curve, and crystal properties such as defect energies and barriers. We perform a 1.3M atom shock simulation and show that ANI-Al force predictions shine in their agreement with new reference DFT calculations.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
eReaxFF force field development for BaZr 0.8 Y 0.2 O 3-δ solid oxide electrolysis cells applications
The use of solid-oxide materials in electrocatalysis applications, especially in hydrogen-evolution reactions, is promising. However, further improvements are warranted to overcome the fundamental bottlenecks to enhancing the performance of solid-oxide electrolysis cells (SOECs), which is directly linked to the more-refined fundamental understanding of complex physical and chemical phenomena and mass exchanges that take place at the surfaces and in the bulk of electrocatalysis materials. Here, we developed an eReaxFF force field for barium zirconate doped with 20 mol% of yttrium, BaZr 0.8 Y 0.2 O 3-δ (BZY20) to enable a systematic, large-length-scale, and longer-timescale atomistic simulation of solid-oxide electrocatalysis for hydrogen generation. All parameters for the eReaxFF were optimized to reproduce quantum-mechanical (QM) calculations on relevant condensed phase and cluster systems describing oxygen vacancies, vacancy migrations, electron localization, water adsorption, water splitting, and hydrogen generation on the surfaces of the BZY20 solid oxide. Using the developed force field, we performed both zero-voltage (excess electrons absent) and non-zero-voltage (excess electrons present) molecular dynamics simulations to observe water adsorption, water splitting, proton migration, oxygen-vacancy migrations, and eventual hydrogen-production reactions. Based on investigations offered in the present study, we conclude that the eReaxFF force field-based approach can enable computationally efficient simulations for electron conductivity, electron leakage, and other non-zero-voltage effects on the solid oxide materials using the explicit-electron concept. Moreover, we demonstrate how the eReaxFF force field-based atomistic-simulation approach can enhance our understanding of processes in SOEC applications and potentially other renewable-energy applications.
08 HYDROGEN↗
The ANI-1ccx and ANI-1x data sets, coupled-cluster and density functional theory properties for molecules
Abstract Maximum diversification of data is a central theme in building generalized and accurate machine learning (ML) models. In chemistry, ML has been used to develop models for predicting molecular properties, for example quantum mechanics (QM) calculated potential energy surfaces and atomic charge models. The ANI-1x and ANI-1ccx ML-based general-purpose potentials for organic molecules were developed through active learning; an automated data diversification process. Here, we describe the ANI-1x and ANI-1ccx data sets. To demonstrate data diversity, we visualize it with a dimensionality reduction scheme, and contrast against existing data sets. The ANI-1x data set contains multiple QM properties from 5 M density functional theory calculations, while the ANI-1ccx data set contains 500 k data points obtained with an accurate CCSD(T)/CBS extrapolation. Approximately 14 million CPU core-hours were expended to generate this data. Multiple QM calculated properties for the chemical elements C, H, N, and O are provided: energies, atomic forces, multipole moments, atomic charges, etc. We provide this data to the community to aid research and development of ML models for chemistry.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Bond order predictions using deep neural networks
Machine learning is an extremely powerful tool for the modern theoretical chemist since it provides a method for bypassing costly algorithms for solving the Schrödinger equation. Already, it has proven able to infer molecular and atomic properties such as charges, enthalpies, dipoles, excited state energies, and others. Most of these machine learning algorithms proceed by inferring properties of individual atoms, even breaking down total molecular energy into individual atomic contributions. In this paper, we introduce a modified version of the Hierarchically Interacting Particle Neural Network (HIP-NN) capable of making predictions on the bonds between atoms rather than on the atoms themselves. We train the modified HIP-NN to infer bond orders for a large number of small organic molecules as computed via the Natural Bond Orbital package. We demonstrate that the trained model is extensible to molecules much larger than those in the training set by studying its performance on the COMP6 dataset. This method has applications in cheminformatics and force field parameterization and opens a promising future for machine learning models to predict other quantities that are defined between atoms such as density matrix elements, Hamiltonian parameters, and molecular reactivities.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Coupling of torsion and OH-stretching in tert -butyl hydroperoxide. II. The OH-stretching fundamental and overtone spectra
We report the vibrational spectra of gas phase tert-butyl hydroperoxide have been recorded in the OH-stretching fundamental and overtone regions (Δv OH = 1–5) at room temperature using conventional Fourier transform infrared (Δv OH = 1–3) and cavity ring-down (Δv OH = 4–5) spectroscopy. In hydroperoxides, the OH-stretching and COOH torsion vibrations are strongly coupled. The double-well nature of the COOH torsion potential leads to tunneling splitting of the energy levels and, combined with the low frequency of the torsional vibration, results in spectra in the OH-stretching regions with multiple vibrational transitions. In each of the OH-stretching regions, both an OH-stretching and a stretch–torsion combination feature are observed, and we show direct evidence for the tunneling splitting in the OH-stretching fundamental region. We have developed two complementary vibrational models to describe the spectra of the OH-stretching regions, a reaction path model and a reduced dimensional local mode model, both of which describe the features of the vibrational spectra well. We also explore the torsional dependence of the OH-stretching transition dipole moment and show that a Franck–Condon treatment fails to capture the intensity in the region of the stretch–torsion combination features. The accuracy of the Franck–Condon treatment of these features improves with increasing Δv OH .
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Machine Learned Hückel Theory: Interfacing Physics and Deep Neural Networks
The Hückel Hamiltonian is an incredibly simple tight-binding model known for its ability to capture qualitative physics phenomena arising from electron interactions in molecules and materials. Part of its simplicity arises from using only two types of empirically fit physics-motivated parameters: the first describes the orbital energies on each atom and the second describes electronic interactions and bonding between atoms. By replacing these empirical parameters with machine-learned dynamic values, we vastly increase the accuracy of the extended Hückel model. The dynamic values are generated with a deep neural network, which is trained to reproduce orbital energies and densities derived from density functional theory. The resulting model retains interpretability, while the deep neural network parameterization is smooth and accurate and reproduces insightful features of the original empirical parameterization. Altogether, this work shows the promise of utilizing machine learning to formulate simple, accurate, and dynamically parameterized physics models.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Quantum mechanics based non-bonded force field functions for use in molecular dynamics simulations of materials and systems: The nitrogen and oxygen columns
Accurate Force Fields (FFs) are essential for Molecular Dynamics (MD) simulations of the dynamics of realistic materials in terms of atomic-level interactions. The FF parameters of short-range valence interactions can be derived through Quantum Mechanical (QM) calculations on model systems practical for QM (<300 atoms). Similarly, the dynamic electrostatic interactions can be described with methods such as QEq or PQEq that allow charges and polarization to adjust dynamically. However, accurately extracting long-range van der Waals (vdW) interactions from QM calculations poses challenges due to the absence of a definitive method to distinguish between the different energetic components of electrostatics, polarization, vdW, hydrogen bonding, and valence interactions. To do this we use the Perdew–Burke–Ernzerhof flavor of Density Functional Theory, including empirical D3 vdW corrections, to predict the Equation of State for each element (keeping any covalent bonds fixed), from which we obtain the two-body vdW nonbond potential. Here, we extend these calculations to include non-bonded parameters for the N and O columns of the periodic table so that we now describe columns 15 (N), 16 (O), 17 (F), and 18 (Ne) of the periodic table. For these 20 elements, we find that the two-body vdW potentials can all be mapped to a single universal two-body curve, with just three scaling parameters: Re, De, and L. We refer to this as the Universal NonBond (UNB) potential. We expect this to be useful for new MD simulations and a helpful starting point to obtain UNB parameters for the remainder of the periodic table.
Chemistry↗
Non-Hermitian quantum mechanics approach for extracting and emulating continuum physics based on bound-state-like calculations: Detailed description
Here, this work applies a reduced basis method to study the continuum physics of a finite quantum system—either few or many-body. Specifically, I develop reduced-order models, or emulators, for the underlying inhomogeneous Schrödinger equation and train the emulators against the equation's bound-state-like solutions at complex energies. The emulators rapidly and accurately interpolate and extrapolate the matrix elements of the Hamiltonian resolvent operator (Green's function) across a parameter space that includes both complex energy and other real-valued physical inputs in the Schrödinger equation. The spectra, discretized and compressed as the result of emulation, and the associated resolvent matrix elements (or amplitudes), have the defining characteristics of non-Hermitian quantum mechanics calculations, featuring complex eigenenergies with negative imaginary parts and branch cuts moved below the real axis in the complex energy plane. Therefore, one now has a method that extracts continuum physics from bound-state-like calculations and emulates those extractions in the input parameter space. Building on a prior Letter [Zhang, Phys. Rev. Lett. 135, 242501 (2025)], this article provides the full theoretical details, a comprehensive analysis of the method's performance, and a brief discussion of how it can be coupled with existing continuum approaches to perform emulations in their input parameter spaces.
ab initio calculations↗
Protonation of histidine rings using quantum-mechanical methods
Histidine can be protonated on either or both of the two N atoms of the imidazole moiety. Each of the three possible forms occurs as a result of the stereochemical environment of the histidine side chain. In an atomic model, comparing the possible protonation states in situ, looking at possible hydrogen bonding and metal coordination, it is possible to predict which is most likely to be correct. A more direct method is described that uses quantum-mechanical methods to calculate, also in situ, the minimum geometry and energy for comparison, and therefore to more accurately identify the most likely protonation state.
Histidine↗
Quantum Disordered State of Magnetic Charges in Nanoengineered Honeycomb Lattice
Abstract A quantum magnetic state due to magnetic charges is never observed, even though they are treated as quantum mechanical variables in theoretical calculations. Here, the occurrence of a novel quantum disordered state of magnetic charges in a nanoengineered magnetic honeycomb lattice of ultra‐small connecting elements is demonstrated. The experimental research, performed using spin resolved neutron scattering, reveals a massively degenerate ground state, comprised of low integer and energetically forbidden high integer magnetic charges, that manifests cooperative paramagnetism at low temperature. The system tends to preserve the degenerate configuration even under large magnetic field application. It exemplifies the robustness of disordered correlation of magnetic charges in a 2D honeycomb lattice. The realization of quantum disordered ground state elucidates the dominance of exchange energy, which is enabled due to the nanoscopic magnetic element size in nanoengineered honeycomb. Consequently, an archetypal platform is envisaged to study quantum mechanical phenomena due to emergent magnetic charges.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
A derivation of Born’s rule from symmetry
Highlights: • A very simple proof for the rule by which we calculate probabilities in quantum mechanics. We here put forward a simple argument for Born’s rule based on the requirement that transition probabilities in quantum mechanics should not be a function of the number of degrees of freedom in the system, or the dimension of the Hilbert-space, respectively. The argument presented here can be interpreted as deriving Born’s rule from a symmetry requirement.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Experimental and Computational Studies of the Production of 1,3-Butadiene from 2,3-Butanediol Using SiO2-Supported H3PO4 Derivatives
Silica-supported phosphoric acid and metal phosphate catalyzed 1,3-butadiene (BDE) production from 2,3-butanediol (2,3-BDO) was studied using experimental and computational techniques. The catalyst was initially tested in a continuous flow reactor using commercially available 2,3-BDO, leading to maximum BDE yields of 63C%. Quantum chemical mechanistic studies revealed 1,2-epoxybutane is a kinetically viable and thermodynamically stable intermediate, supported by experimental demonstration that this epoxide can be converted to BDE under standard reaction conditions. Newly proposed E2 and SN2' elementary steps were studied to rationalize the formation of BDE and all detected side-products. Additionally, using quantum mechanics/molecular mechanics (QM/MM) calculations, we modeled silica-supported phosphate catalysts to study the effect of the alkali metal center. Natural population analysis showed that phosphate oxygen atoms are more negatively charged in CsH2PO4/SiO2 than in H3PO4/SiO2. In combination with temperature-programmed desorption experiments using CO2, the results of this study suggest that the improved selectivity achieved when adding the metal center is related to an increase in the basicity of the catalyst.
butadiene↗
Understanding how chemical structure affects ignition-delay-time $\phi$-sensitivity
$\phi$-sensitivity is the change in ignition delay time (IDT) with respect to the fuel-to-air equivalence ratio ($\phi$). High $\phi$-sensitivity is a desirable fuel property for applications in advanced compression ignition and multi-mode engine designs. Understanding how $\phi$-sensitivity depends on chemical structure is essential for selecting promising biofuels from the ever-growing list of proposed candidates. Here, we investigate the effect of chemical structure on $\phi$-sensitivity with experiment, simulation, and theory. Experimental Advanced Fuel Ignition Delay Analyzer (AFIDA) measurements for 2,4-dimethylpentane and diisopropyl ether provide evidence that branching and functional groups strongly impact $\phi$-sensitivity. Further insights into this dependence are obtained with 0-D kinetic simulations with existing mechanisms for n-pentane, diethyl ether, 3-pentanone, n-heptane, 2-methylhexane, 2,4-dimethylpentane, and 2,2,3-trimethylbutane. Quantum mechanical (QM) G4 calculations of low-temperature reactions help explain the observed experimental and simulation trends. Specifically, these QM calculations provide theoretical estimates of the ketohydroperoxide (KHP) dissociation rates, the HO 2 formation rates from peroxy radical (ROO), and the “cross-over” temperatures, i.e., the temperature at which ROO dissociation is favored compared to hydroperoxyl radical (QOOH) formation. Each of these reaction rates is compared to the n-alkane reference point to determine the impact of branching and different functional groups. Although kinetic mechanisms typically assume that KHP dissociation rates are invariant of chemical environment, our QM results suggest that this rate can span a range of roughly two orders of magnitude. We also discuss the importance of including the peroxy-hydroperoxy (OO-OOH) hydrogen transfer reaction for branched ethers. Finally, the insights gained assist in proposing a highly $\phi$-sensitive compound, namely, isopropyl propyl ether.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗