Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “corrections”

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 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Measurement of groomed jet substructure observables in p+p collisions at $\sqrt{s} =$ 200 GeV with STAR

In this letter, measurements of the shared momentum fraction ( z g ) and the groomed jet radius ( R g ), as defined in the SoftDrop algorithm, are reported in p+p collisions at s = 200 GeV collected by the STAR experiment. These substructure observables are differentially measured for jets of varying resolution parameters from R = 0.2 – 0.6 in the transverse momentum range 15 < p T , jet < 60 GeV/ c . These studies show that, in the p T , jet range accessible at s = 200 GeV and with increasing jet resolution parameter and jet transverse momentum, the z g distribution asymptotically converges to the DGLAP splitting kernel for a quark radiating a gluon. The groomed jet radius measurements reflect a momentum-dependent narrowing of the jet structure for jets of a given resolution parameter, i.e., the larger the p T , jet , the narrower the first splitting. For the first time, these fully corrected measurements are compared to Monte Carlo generators with leading order QCD matrix elements and leading log in the parton shower, and to state-of-the-art theoretical calculations at next-to-leading-log accuracy. We observe that PYTHIA 6 with parameters tuned to reproduce RHIC measurements is able to quantitatively describe data, whereas PYTHIA 8 and HERWIG 7, tuned to reproduce LHC data, are unable to provide a simultaneous description of both z g and R g , resulting in opportunities for fine parameter tuning of these models for p+p collisions at RHIC energies. We also find that the theoretical calculations without non-perturbative corrections are able to qualitatively describe the trend in data for jets of large resolution parameters at high p T , jet , but fail at small jet resolution parameters and low jet transverse momenta.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Measurement of the $t\bar{t}$ charge asymmetry in events with highly Lorentz-boosted top quarks in pp collisions at $\sqrt{s}$ = 13 TeV

The measurement of the charge asymmetry in top quark pair events with highly Lorentz-boosted top quarks decaying to a single lepton and jets is presented. The analysis is performed using proton-proton collisions at $\sqrt{s}$ = 13 TeV with the CMS detector at the LHC and corresponding to an integrated luminosity of 138 fb -1 . The selection is optimized for top quarks produced with large Lorentz boosts, resulting in nonisolated leptons and overlapping jets. The top quark charge asymmetry is measured for events with a $t\bar{t}$ invariant mass larger than 750 GeV and corrected for detector and acceptance effects using a binned maximum likelihood fit. The measured top quark charge asymmetry of ($0.42_{-0.69}^{+0.64}$) % is in good agreement with the standard model prediction at next-to-next-to-leading order in quantum chromodynamic perturbation theory with next-to-leading-order electroweak corrections. The result is also presented for two invariant mass ranges, 750–900 and > 900 GeV .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Spontaneous Raman–LIF–CO–OH measurements of species concentration in turbulent spray flames

This paper presents new measurements of species concentrations, temperature and mixture fraction in selected regions of a turbulent ethanol spray flame. The line-Raman–LIF–COsingle bondOH setup developed at the Sandia's Combustion Research Facility is utilised to probe regions of a spray flame where laser breakdown of liquid droplets is avoided and the remaining interferences can be corrected. The spray flame is stabilised on the piloted Sydney needle spray burner, where axial translation of the liquid injecting needle in the air-blast stream can transition the spray from dilute to dense. The solution to obtaining successful measurements is found to be multifaceted and includes: the appropriate selection of flame conditions; high sensitivity of the Raman detection system permitting reduced laser energies; development of a pre-processing algorithm to reject strong droplet interferences; and application of the hybrid matrix inversion method combined with wavelet denoising to account for interference corrections and noise at the very low signal levels obtained. Unique and necessary for the successful measurements reported in this paper, a pre-processing algorithm is outlined that removes data points corrupted with strong interferences from droplets. These interferences arise from a range of sources, but the most intense are due to the laser interaction with surrounding mist or liquid fragments, such that measurements near the jet centreline are corrupted and hence discarded. Reliable measurements of mixture fraction, temperature obtained from the sum of the species number densities, and species mole fractions are reported for regions in the flames sufficiently far from the centreline. The paper demonstrates the feasibility of the judicious use of Raman scattering in turbulent spray flames, the results of which will be extremely useful for validating numerical simulations.

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

Contributions of vegetation heterogeneity within tower footprint to CO 2 flux estimations through graph neural network modeling

Net ecosystem exchange of CO 2 (Fc) measured directly by eddy covariance towers is based on various assumptions, including large, flat and homogenous land cover type. In reality, often a tower site is not large enough for flux measurements, and landscapes consist of patches of different land cover types within the flux footprint. In addition, some portions of fluxes are contributed by different cover types when a footprint exceeds the size of the target ecosystem. The contributions of non-dominant patches to Fc are often ignored. Here, in this study, we propose a novel integrated modeling framework that combines random forest (RF) and XGBoost with a residual correction module based on a deep graph convolutional network (DeeperGCN) to simulate Fc for seven flux measurement sites in southwest Michigan. High-resolution remote sensing vegetation indices, soil properties, meteorological variables, and footprint-weighted spatial features were used as model inputs at three spatial resolutions (10 m, 20 m, 30 m), and their importance in predicting Fc with DeeperGCN was assessed. We found that residual correction using DeeperGCN significantly improved prediction accuracy, with the R 2 increasing from 0.9098 to 0.9479 for RF and from 0.9235 to 0.9433 for XGBoost. At site level, the maximum improvement in R 2 reached 0.1617. Paired t-tests confirmed that these improvements were statistically significant (p < 0.05). Among all predictors, leaf area index and incoming shortwave radiation emerged as the dominant drivers of spatial residual variation, followed by precipitation, relative humidity, and selected vegetation indices. The 20 m resolution yielded the best balance between model performance and computational efficiency. In conclusion, our modeling framework effectively captures both spatial heterogeneity and nonlinear interactions, offering a robust solution for spatially explicit flux modeling in structurally diverse ecosystems beyond the study sites.

footprint model↗
DOE OSTI2025

Pointing stabilization of a 1 Hz high-power laser via machine learning

Abstract High-power lasers are vital for particle acceleration, imaging, fusion and materials processing, requiring precise control and high-energy delivery. Laser plasma accelerators (LPAs) demand laser positional stability at focus to ensure consistent electron beams in applications such as X-ray free-electron lasers and high-energy colliders. Achieving this stability is especially challenging for the low-repetition-rate lasers in current LPAs. We present a machine learning method that predicts and corrects laser pointing instabilities in real-time using a high-frequency pilot beam. By preemptively adjusting a correction mirror, this approach overcomes traditional feedback limits. Demonstrated on the BELLA petawatt laser operating at the terawatt level (30 mJ amplification), our method achieved root mean square pointing stabilization of 0.34 and 0.59 $\unicode{x3bc} \mathrm{rad}$ in the x and y directions, reducing jitter by 65% and 47%, respectively. This is the first successful application of predictive control for shot-to-shot stabilization in low-repetition-rate laser systems, paving the way for full-energy petawatt lasers and transformative advances across science, industry and security.

Amodio, Alessio↗
DOE OSTI2021

Toward a First-Principles Framework for Predicting Collective Properties of Electrolytes

Conspectus Given the universal importance of electrolyte solutions, it is natural to expect that we have a nearly complete understanding of the fundamental properties of these solutions (e.g., the chemical potential) and that we can therefore explain, predict, and control the phenomena occurring in them. In fact, reality falls short of these expectations. But, recent advances in the simulation and modeling of electrolyte solutions indicate that it should soon be possible to make progress toward these goals. In this Account, we will discuss the use of first-principles interaction potentials based in quantum mechanics (QM) to enhance our understanding of electrolyte solutions. Specifically, we will focus on the use of quantum density functional theory (DFT) combined with molecular dynamics simulation (DFT-MD) as the foundation for our approach. The overarching concept is to understand and accurately reproduce the balance between local or short-ranged (SR) structural details and long-range (LR) correlations, allowing the prediction of the thermodynamics of both single ions in solution as well as the collective interactions characterized by activity/osmotic coefficients. In doing so, relevant collective motions and driving forces characterized by chemical potentials can be determined. Here, in this Account, we will make the case that understanding electrolyte solutions requires a faithful QM representation of the SR nature of the ion–ion, ion–water, and water–water interactions. However, the number of molecules that is required for collective behavior makes the direct application of high-level QM methods that contain the best SR physics untenable, making methods that balance accuracy and efficiency a practical goal. Alternatives such as continuum solvent models (CSMs) and empirically based classical molecular dynamics have been extensively employed to resolve this problem but without yet overcoming the fundamental issue of SR accuracy. We will demonstrate that accurately describing the SR interaction is imperative for predicting both intrinsic properties, namely, at infinite dilution, and collective properties of electrolyte solutions. DFT has played an important role in our understanding of condensed phase systems, e.g., bulk liquid water, the air–water interface, ions in bulk, and at the air–water interface. This approach holds huge promise to provide benchmark calculations of electrolyte solution properties that will allow for the development and improvement of more efficient methods, as well as an enhanced understanding of fundamental phenomena. However, the standard protocol using the generalized gradient approximation with van der Waals (vdW) correction requires improvement in order to achieve a high level of quantitative accuracy. Simply simulating with higher level DFT functionals may not be the best route considering the significant computational cost. Alternative methods of incorporating information from higher levels of QM should be explored; e.g., using force matching techniques on small clusters, where high level benchmark calculations are possible, to develop ideal correction terms to the DFT functional is a promising possibility. We argue that DFT with statistical mechanics is becoming an increasingly useful framework enabling the prediction of collective electrolyte properties.

Duignan, Timothy T.↗
DOE OSTI2020

Structural Tolerance Factor Approach to Defect-Resistant I 2 -II-IV-X 4 Semiconductor Design

Recent work on quaternary semiconductors Cu 2 BaSn(S,Se) 4 and Ag 2 BaSnSe 4 for photovoltaic and thermoelectric applications, respectively, has shown the promise of exploring the broader family of defect-resistant I 2 -II-IV-X 4 materials (where I, II, and IV refer to the formal oxidation state of the metal cations and X is a chalcogen anion) with tetrahedrally coordinated I/IV cations and larger II cations (i.e., Sr, Ba, Pb, and Eu) for optoelectronic and energy-related applications. Chemical dissimilarity among the II and I/IV atoms represents an important design motivation because it presents a barrier to antisite formation, which otherwise may act as electronically harmful defects. We herein show how all 31 experimentally reported I 2 -II-IV-X 4 examples (with large II cations and tetrahedrally coordinated smaller I/IV cations), which form within five crystal structure types, are structurally linked. Based on these structural similarities, we derive a set of tolerance factors that serve as descriptors for phase stability within this family. Despite common usage in the well-studied perovskite system, Shannon ionic radii are found to be insufficient for predicting metal–chalcogen bond lengths, pointing to the need for experimentally derived correction factors as part of an empirically driven learning approach to structure prediction. Here, we use the tolerance factors as a predictive tool and demonstrate that four new I 2 -II-IV-X 4 compounds, Ag 2 BaSiS 4 , Ag 2 PbSiS 4 , Cu 2 PbGeS 4 , and Cu 2 SrSiS 4 , can be synthesized in correctly predicted phases. One of these compounds, Ag 2 PbSiS 4 , shows potentially promising optoelectronic properties for photovoltaic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Insights into Controls of Mineralogy and Pore Structure on the Density of Methane Adsorption Phase in Shales under Supercritical Conditions

With gradually increasing enhancing population and industrialization, exploitation of shale gas is enhancing in the world to satisfy the growing demand for energy worldwide. To exploit shale gas efficiently, accurate assessment of shale gas in place (GIP) is necessary for determining production strategies. Unlike free and dissolved gas, adsorbed gas contributes to the shale GIP up to 85% due to the well-developed surface area and micropores. Measured excess adsorption of methane (up to over 95% of subsurface shale gas) is often corrected to absolute adsorption to obtain the actual amount of adsorption under supercritical (or geological) conditions. During the correction, adsorption phase density (APD) is critical. However, the APD of subsurface shale gas and the effects of shale properties (e.g., mineralogy and pore structure) on APD remain poorly understood. A series of high-pressure methane adsorption isotherms (HPMAI) on Caney Shales were collected and analyzed in conjunction with other United States and Chinese shales from the literature at 35–125 °C and up to 15 MPa. Here, a three-layer Ono–Kondo (OK3) model is utilized to derive the temperature- and pressure-dependent APD coupling low-pressure nitrogen adsorption isotherms (LPNAI) and HPMAI. X-ray diffraction and organic geochemistry are combined to reveal the mineralogy. Brunauer–Emmett–Teller, Barrett–Joyner–Halenda, and Horvath–Kawazoe analyses via LPNAI are used to investigate the pore structures. Results show that APD increases with organic matter (OM) proxied by total organic carbon and decreases with clay minerals and the sum of quartz and feldspar. OM dramatically contributes to the APD as multiple-layer adsorption exists, and the APD for OM could be 1.4–8.5 times that for clay minerals. Other inorganic minerals contribute less to APD. The properties and constitution of the surface area instead of the volume fraction contribute to the APD in shales. APD does not show an obvious correlation with micropore volume, likely related to the ratio of micropore volume to the total pore volume. Here, we provide a significant and comprehensive study of petrological factors that impact the APD of subsurface shale gas, which will improve the estimation of supercritical adsorption and shale GIP under reservoir conditions. Also, the findings in this work can provide applications for subsurface carbon dioxide adsorption and storage.

58 GEOSCIENCES↗
DOE OSTI2022

Ecological Traits of Fish for Mercury Biomonitoring: Insights from Compound-Specific Nitrogen and Stable Mercury Isotopes

Here we coupled compound-specific isotopic analyses of nitrogen (N) in amino acids (δ 15 N Glu , δ 15 N Phe ) and mercury stable isotopes (δ 202 Hg, Δ 199 Hg) to quantify ecological traits governing the concentration, variability, and source of Hg in largemouth bass (LB) and pike gudgeon (PG) across four rivers, South Korea. PG displayed uniform Hg concentration (56–137 ng/g), trophic position (TP corrected ; 2.6–3.0, n = 9), and N isotopes in the source amino acid (δ 15 N Phe ; 7–13‰), consistent with their specialist feeding on benthic insects. LB showed wide ranges in Hg concentration (45–693 ng/g), TP corrected (2.8–3.8, n = 14), and δ 15 N Phe (1.3–16‰), reflecting their opportunistic feeding behavior. Hg sources assessed using Hg isotopes reveal low and uniform Δ 199 Hg in PG (0.20–0.49‰), similar to Δ 199 Hg reported in sediments. LB displayed site-specific δ 202 Hg (-0.61 to -0.04‰) and Δ 199 Hg (0.53–1.09‰). At the Yeongsan River, LB displayed elevated Δ 199 Hg and low δ 15 N Phe , consistent with Hg and N sourced from the atmosphere. LB at the Geum River displayed low Δ 199 Hg and high δ 15 N Phe , both similar to the isotope values of anthropogenic sources. Our results suggest that a specialist fish (PG) with consistent ecological traits and Hg concentration is an effective bioindicator species for Hg. When accounting for Hg sources, however, LB better captures site-specific Hg sources.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

How Beneficial Is Pretraining on a Narrow Domain-Specific Corpus for Information Extraction about Photocatalytic Water Splitting?

Language models trained on domain-specific corpora have been employed to increase the performance in specialized tasks. However, little previous work has been reported on how specific a “domain-specific” corpus should be. Here, we test a number of language models trained on varyingly specific corpora by employing them in the task of extracting information from photocatalytic water splitting. We find that more specific corpora can benefit performance on downstream tasks. Furthermore, PhotocatalysisBERT, a pretrained model from scratch on scientific papers on photocatalytic water splitting, demonstrates improved performance over previous work in associating the correct photocatalyst with the correct photocatalytic activity during information extraction, achieving a precision of 60.8(+11.5)% and a recall of 37.2(+4.5)%.

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

KBKit: A Python Toolkit for Kirkwood–Buff Theory from Molecular Dynamics

Thermodynamic properties of liquid mixtures govern processes that range from drug delivery to energy storage, yet extracting these properties from molecular simulations remains challenging. Kirkwood–Buff (KB) theory offers a rigorous route by linking microscopic pair distribution functions to macroscopic free energies, but practical use of the theory has been hindered by two obstacles: (i) the long simulations needed to obtain well-converged Kirkwood-Buff integrals (KBIs) and (ii) the specialized corrections required to translate finite-size data to the thermodynamic limit. $\texttt{KBKit}$ is an open-source Python package that removes these barriers. It automatically computes KBIs and derived thermodynamic quantities from GROMACS input files, applies state-of-the-art finite-size corrections, and provides built-in diagnostic tools to quantify statistical uncertainty. Written with modern software-engineering practices—continuous integration, extensive unit testing, and thorough documentation—$\texttt{KBKit}$ is both reliable and easy to extend. By condensing complex KBI analysis into a few intuitive commands, $\texttt{KBKit}$ enables researchers to incorporate KB theory into routine simulation workflows and accelerate the discovery of solution-phase thermodynamics.

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

Revisiting Competing Paths in Electrochemical CO 2 Reduction on Copper via Embedded Correlated Wavefunction Theory

We re-evaluate two key steps in the mechanism of CO 2 reduction on copper at a higher level of theory capable of correcting inherent errors in density functional theory (DFT) approximations, namely, embedded correlated wavefunction (ECW) theory. Here, we consider the CO reduction step on Cu(111), which is critical to understanding reaction selectivity. Here, we optimize embedding potentials at the periodic plane-wave DFT level using density functional embedding theory (DFET). All possible adsorption sites (adsites) for each adsorbate then are screened with ECW theory at the catalytically active site to refine the local electronic structure. Unsurprisingly, DFT and ECW theory predict different adsite preferences, largely because of DFT’s inability to properly situate the CO 2π* level. Differing preferred adsites suggest that different reaction pathways could emerge from DFT versus ECW theory. Starting from these preferred ECW theory adsites, we then obtain reaction pathways at the plane-wave DFT level using the climbing-image nudged elastic band method to determine minimum energy paths. Thereafter, we perform ECW calculations at the catalytically active site to correct the energetics at each interpolated structure (image) along the reaction pathways. Via this approach, we confirm that the first step in CO reduction via hydrogen transfer on Cu(111) is to form hydroxymethylidyne (*COH) instead of formyl (*CHO). Although the prediction to preferentially form *COH is consistent with that of DFT, the two theories predict quite different structural and mechanistic behaviors, suggesting that verification is needed for other parts of the mechanism of CO 2 reduction, which is the subject of ongoing work.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

The many-body expansion for aqueous systems revisited: II. Alkali metal and halide ion – water interactions

We present a detailed study of the Many-Body Expansion (MBE) for alkali metal and halide ion-water interactions and quantify the effect of these ions on the strength of the surrounding aqueous hydrogen bonding environment. Building on our previous work on neutral water clusters [J. P. Heindel and S. S. Xantheas, J. Chem. Theor. Comput. 16 (11), 6843–6855 (2020)], we carry out the complete MBE for ion-water clusters, Z+/-(H2O)9, where Z = Li+, K+, Cs+, Cl-, Br-, I- and compare with the results for (H2O)10. The 2-B ion-water (I-W) interaction represents a larger percentage of the total cluster binding energy compared to a pure water cluster of the same size with the total 3-B term being smaller and of opposite sign (repulsive) whereas higher order terms are essentially negligible. The same oscillating behavior around zero for MBE terms higher than the 5-B with basis set that was reported for water clusters is also observed for the ion-water clusters considered with Basis Set Superposition Error (BSSE) corrections amending this like in the water cluster case. A remarkable, linear anti-correlation between the total 2-B ion-water and the total 2-B and 3-B water-water interactions is found, quantifying the effect of the different ions in disrupting and altering (weakening) the neighboring hydrogen bonded water network. Our results suggest a universal behavior of the two different families of ions (alkali metals and halides) for both the correlations of the various components of the total binding energies as well as the estimate of the 2-B BSSE correction which is reported to follow a common profile for ion-water and water-water interactions when cast in terms of reduced distances and energies. We expect our insights into the nature of BSSE will be important for future ab initio based, many-body molecular dynamics studies.

ion-water interactions, water-water interaction, m↗
DOE OSTI2021

Phase Equilibrium of Water with Hexagonal and Cubic Ice Using the SCAN Functional

Machine learning models are rapidly becoming widely used to simulate complex physicochemical phenomena with ab initio accuracy. Here, we use one such model as well as direct density functional theory (DFT) calculations to investigate the phase equilibrium of water, hexagonal ice (Ih), and cubic ice (Ic), with an eye toward studying ice nucleation. The machine learning model is based on deep neural networks and has been trained on DFT data obtained using the SCAN exchange and correlation functional. We use this model to drive enhanced sampling simulations aimed at calculating a number of complex properties that are out of reach of DFT-driven simulations and then employ an appropriate reweighting procedure to compute the corresponding properties for the SCAN functional. This approach allows us to calculate the melting temperature of both ice polymorphs, the driving force for nucleation, the heat of fusion, the densities at the melting temperature, the relative stability of ices Ih and Ic, and other properties. We find a correct qualitative prediction of all properties of interest. In some cases, quantitative agreement with experiment is better than for state-of-the-art semiempirical potentials for water. Furthermore, our results also show that SCAN correctly predicts that ice Ih is more stable than ice Ic.

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

A Massively Parallel Implementation of the CCSD(T) Method Using the Resolution-of-the-Identity Approximation and a Hybrid Distributed/Shared Memory Parallelization Model

In this work, a parallel algorithm is described for the coupled-cluster singles and doubles method augmented with a perturbative correction for triple excitations [CCSD(T)] using the resolution-of-the-identity (RI) approximation for two-electron repulsion integrals (ERIs). The algorithm bypasses the storage of four-center ERIs by adopting an integral-direct strategy. The CCSD amplitude equations are given in a compact quasi-linear form by factorizing them in terms of amplitude-dressed three-center intermediates. A hybrid MPI/OpenMP parallelization scheme is employed, which uses the OpenMP-based shared memory model for intranode parallelization and the MPI-based distributed memory model for internode parallelization. Parallel efficiency has been optimized for all terms in the CCSD amplitude equations. Two different algorithms have been implemented for the rate-limiting terms in the CCSD amplitude equations that entail and -scaling computational costs, where N O and N V denote the number of correlated occupied and virtual orbitals, respectively. One of the algorithms assembles the four-center ERIs requiring N V 4 and N O 2 N V 2 -scaling memory costs in a distributed manner on a number of MPI ranks, while the other algorithm completely bypasses the assembling of quartic memory-scaling ERIs and thus largely reduces the memory demand. It is demonstrated that the former memory-expensive algorithm is faster on a few hundred cores, while the latter memory-economic algorithm shows a better strong scaling in the limit of a few thousand cores. The program is shown to exhibit a near-linear scaling, in particular for the compute-intensive triples correction step, on up to 8000 cores. The performance of the program is demonstrated via calculations involving molecules with 24–51 atoms and up to 1624 atomic basis functions. As the first application, the complete basis set (CBS) limit for the interaction energy of the π-stacked uracil dimer from the S66 data set has been investigated. This work reports the first calculation of the interaction energy at the CCSD(T)/aug-cc-pVQZ level without local orbital approximation. The CBS limit for the CCSD correlation contribution to the interaction energy was found to be -8.01 kcal/mol, which agrees very well with the value -7.99 kcal/mol reported by Schmitz, Hättig, and Tew [ Phys. Chem. Chem. Phys. 2014 , 16 , 22167-22178]. The CBS limit for the total interaction energy was estimated to be -9.64 kcal/mol.

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

Comparing the Expense and Accuracy of Methods to Simulate Atomic Vibrations in Rubrene

Atomic vibrations can inform about materials properties from hole transport in organic semiconductors to correlated disorder in metal–organic frameworks. Currently, there are several methods for predicting these vibrations using simulations, but the accuracy–efficiency tradeoffs have not been examined in depth. In this work, rubrene is used as a model system to predict atomic vibrational properties using six different simulation methods: density functional theory, density functional tight binding, density functional tight binding with a Chebyshev polynomial-based correction, a trained machine learning model, a pretrained machine learning model called ANI-1, and a classical forcefield model. The accuracy of each method is evaluated by comparison to the experimental inelastic neutron scattering spectrum. All methods discussed here show some accuracy across a wide energy region, though the Chebyshev-corrected tight-binding method showed the optimal combination of high accuracy with low expense. We then offer broad simulation guidelines to yield efficient, accurate results for inelastic neutron scattering spectrum prediction.

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

Single-Point Extrapolation to the Complete Basis Set Limit through Deep Learning

Machine learning (ML) offers an attractive method for making predictions about molecular systems while circumventing the need to run expensive electronic structure calculations. Once trained on ab initio data, the promise of ML is to deliver accurate predictions of molecular properties that were previously computationally infeasible. In this work, we develop and train a graph neural network model to correct the basis set incompleteness error (BSIE) between a small and large basis set at the RHF and B3LYP levels of theory. Our results show that, when compared to fitting to the total potential, an ML model fitted to correct the BSIE is better at generalizing to systems not seen during training. We test this ability by training on single molecules while evaluating on molecular complexes. We also show that ensemble models yield better behaved potentials in situations where the training data is insufficient. However, even when only fitting to the BSIE, acceptable performance is only achieved when the training data sufficiently resemble the systems one wants to make predictions on. As a result, the test error of the final model trained to predict the difference between the cc-pVDZ and cc-pV5Z potential is 0.184 kcal/mol for the B3LYP density functional, and the ensemble model accurately reproduces the large basis set interaction energy curves on the S66x8 dataset.

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

Symmetrized Drude Oscillator Force Fields Improve Numerical Performance of Polarizable Molecular Dynamics

Drude oscillator potentials are a popular and computationally efficient class of polarizable models that represent each polarizable atom as a positively charged Drude core harmonically bound to a negatively charged Drude shell. Here we show that existing force fields that place all non-Coulomb forces on the Drude core and none on the shell inadvertently couple the dipole to non-Coulombic forces. This introduces errors where interactions with neutral particles can erroneously induce atomic polarization, leading to spurious polarizations in the absence of an electric field, exacerbating violations of equipartition in the employed Carr–Parinello scheme. A suitable symmetrization of the interaction potential that correctly splits the force between the Drude core and shell can correct this shortcoming, improving the stability and numerical performance of Drude oscillator-based simulations. The symmetrization procedure is straightforward and only requires the rescaling of a few force field parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗