Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “basis sets”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Restricted Partition Function Semiclassical Transition State Theory (RPF-SCTST): Applications to Reactions Involving Hydrogen Cyanide, Hydrogen Peroxide, and Formaldehyde Oxide

We demonstrate the efficiency of the numerical calculation of thermal semiclassical transition state theory (SCTST) rates across several representative chemical reactions using the restricted partition function (RPF)─viz a function that depends only on the imaginary action associated with a reaction. Here, we treat the potential energy surfaces (PESs) through fourth order expansions around the saddle point and integrate the Hamiltonian up to second order in employing vibrational perturbation theory. We apply this formalism to uni- and bimolecular reactions with different─though relatively high─barrier heights to uncover the influence of the barrier properties, minimal energies, and separability of the rotational motion on rate constants and tunneling corrections. Although all modes are coupled within the RPF, we found in our numerical examples that the rotational component in the absence of significant rotational distortions can be separated from the vibrational contributions. Moreover, a classical treatment of the rotational contribution is adequate over the temperature range from approximately 100 to 1000 K. We also found that the choice of DFT basis set can lead to variations in rate constants of up to an order of magnitude. Lastly, different schemes for counting eligible energy levels result in rate constants that differ by no more than a factor of 2, with the remaining discrepancies attributed to the treatment of near-convergent levels in perturbation theory.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Cation, Anion, and Radical Isomers of C4H4N: Computational Characterization and Implications for Astrophysical and Planetary Environments

Nitrogen-containing ions and molecules in the gas phase have been detected in non-Earth environments such as dark molecular clouds and more recently in the atmosphere of Saturn's moon Titan. These molecules may serve as precursors to larger heterocyclic structures that provide the foundation of complex biological molecules. On Titan, molecules of m/z 66 have been detected by the Cassini mission, and species of the empirical formula C 4 H 4 N may contribute to this signature. We have characterized seven isomers of C 4 H 4 N in anionic, neutral radical, and cationic states using density functional theory. Structures were optimized using the range-separated hybrid ωB97X-V with the cc-pVTZ and aug-cc-pVTZ basis sets. Anionic and radical C4H4N favor cyclic structures with aromatic and quasi-aromatic electron arrangements, respectively. Interestingly, ionization from the radical surface to the cation induces significant changes in structural stability, and the global minimum for positively charged isomers is CH 2 CCHCNH + , a pseudo-linear species reminiscent of cyanoallene. Select formation pathways to these structures from Titan's existing or postulated gas-phase species, reactions that are also relevant for other astrophysical environments, are discussed. By characterizing C 4 H 4 N isomers, we have identified energetically stable anionic, radical, and cationic structures that may be present in Titan's atmosphere and dark molecular clouds.

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

Threshold Collision-Induced Dissociation of Hydrated Thorium(IV) Trihydroxide Cation: Experimental and Theoretical Investigation of the Binding Energies for Th(OH) 3 + (H 2 O) n Complexes ( n = 1–4)

Threshold collision-induced dissociation of Th(OH) 3 + (H 2 O) n Complexes ( n = 1–4) with xenon was performed using a guided ion beam tandem mass spectrometer. The primary dissociation pathway for all complexes is a loss of a single water molecule followed by the sequential loss of additional water molecules at higher collision energies. Here, the data were analyzed using a statistical model after accounting for lifetime effects and reactant internal and kinetic energy distributions to obtain 0 K bond dissociation energies (BDEs). These were also converted using rigid rotor/harmonic oscillator approximations to yield thermodynamic values at room temperature. The 0 K BDEs of H 2 O ligands to Th(OH) 3 + (IV) are experimentally determined for the first time as 106 ± 6, 89 ± 6, 76 ± 4, and 51 ± 4 kJ/mol for the first, second, third, and fourth water ligand added. These values agree reasonably well with values calculated at the B3LYP, B3PW91, and PBE0 levels of theory with aug-cc-pVQZ basis sets, whereas B3LYP-GD3BJ, MP2, and CCSD(T) single point energies with (without) counterpoise corrections systematically overestimate the bond energies by about 15 (20), 19 (25), and (18) kJ/mol, respectively.

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

Thermodynamics of Metal Carbonates and Bicarbonates and Their Hydrates for Mg, Ca, Fe, and Cd Relevant to Mineral Energetics

The heats of formation of carbonate, bicarbonate and bicarbonate/hydroxide metal complexes and hydrated versions these complexes of Mg 2+ , Ca 2+ , Fe 2+ , and Cd 2+ are predicted from atomization energies using correlated molecular orbital theory at the CCSD(T) level extrapolated to the complete basis set limit following the Feller-Peterson-Dixon (FPD) approach. Using the calculated gas phase values and the available experimental solid state values, the cohesive energies of selected minerals were calculated. The gas phase decomposition energies into MO, CO 2 and H 2 O follow the order Mg ~ Ca > Cd ~ Fe, and correlate with the hardness of the metal +2 ions. Gas phase hydration energies show that the order is Mg > Fe > Ca ~ Cd. There are a number of bulk hydrated Mg and Ca complexes that occur as minerals but there are few if any for Fe and Cd, suggesting that a number of factors are important in determining the stability of the bulk mineral hydrates. The FPD heats of formation were used to benchmark a range of density functional theory exchange-correlation functionals, including those commonly used in solid state mineral calculations. None of the functionals provided chemical accuracy agreement (± 1 kcal/mol) with the FPD results The best functionals are ωB97X and ωB97X-D with mean average unsigned errors of 10 kcal/mol.

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

Calculated Reduction Potentials of Electrolyte Species in Lithium–Sulfur Batteries

In this work, reduction potentials of electrolyte molecules in the lithium–sulfur (Li/S) battery and their variations in several solvent environments are studied using the density functional theory method with Dunning’s triple-ζ correlation consistent basis set. Reliable reduction potential values are key for electrolyte additive design needed for suppressing polysulfide dissolution shuttle mechanism, resulting in poor cycle performance and severe self-discharge of the Li/S battery. Although isolated electrolyte molecules have reduction potentials outside the operating voltage range of the Li/S battery, complexation with other electrolyte species enables the electrolyte molecules to be reduced within the operating voltage range. Among the electrolyte species considered in this study, bis(fluorosulfonyl)imide (FSI – ) and fluoroethylene carbonate (FEC) yield reduction potentials within the expected range, suggesting the development of fluorine-containing additives as a promising line of research.

25 ENERGY STORAGE↗
DOE OSTI2022

Interpreting the Operando X-ray Absorption Near-Edge Structure of Supported Cu and CuPd Clusters in Conditions of Oxidative Dehydrogenation of Propane: Dynamic Changes in Composition and Size

Supported subnano-cluster catalysts are highly dynamic, developing true active sites only under the pressures and temperatures of reaction conditions. Operando X-ray absorption near-edge structure (XANES) spectroscopy can track changes in the oxidation state and the local environment of cluster atoms, providing insight into the development of these active sites. While bulk metal, oxide, and hydroxide standards are often used for fitting experimental XANES spectra to obtain average oxidation states, we recently showed that computed cluster standards of relevant compositions are a more suitable basis, producing more accurate fits. Here, we theoretically interpret the operando XANES of supported Cu 3 Pd and Cu 4 clusters during temperature-programmed reaction (TPRx) of oxidative dehydrogenation of propane. We use an expanded basis set including both globally optimized computed clusters and bulk standards. Not only can we track reversible composition/oxidation state change with temperature, but also the irreversible growth of the bulk fraction upon heating, which we attribute to cluster sintering. This has important implications for the mechanism of the catalyzed reaction and the nature of the available active sites. Here, we propose that operando XANES provides most significant insight into the nature of supported cluster catalysts in reaction conditions when interpreted using mixed computed cluster and bulk standards.

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

Integrated Experimental and Computational K-Edge X-ray Absorption Near-Edge Structure Analysis of Vanadium Catalysts

X-ray Absorption Near-Edge Structure (XANES) spectroscopy is a powerful tool to reveal key structural and electronic features of isolated catalytic sites, yet insights into molecular structure and more detailed orbital analysis through a combination of experimental and computed XANES analysis are necessary for accurate interpretation of the spectra, especially when significant heterogeneity exits among the catalytic sites. Herein, we present an integrated computational and experimental strategy to determine both primary and secondary bonding interactions within the XANES pre-edge region for organovanadium complexes, which was developed using a series of well-defined molecular vanadium complexes and then applied to the characterization of a supported organovanadium olefin hydrogenation catalyst. Time-dependent density functional theory is used to predict the energy of pre-edge XANES features for a series of vanadium complexes with a variety of oxidation states and local coordination environments. Further, a calibration scheme incorporating different density functionals and basis sets is established, resulting in an optimized scheme that accurately predicts pre-edge energies with a mean absolute error of 0.40 eV. Second-shell coordination (e.g., V---V) effects within XANES are identified through the analysis of the computed dominant orbital contributions for multi-vanadium complexes. Orbital analysis also provided confirmation that the vanadium-hydride formation combined with the heterogeneity of the catalytic active species in ole-fin hydrogenation caused the energy shift and broadening of the pre-edge peak after hydrogen treatment of the silica-supported organovanadium pre-catalyst. This work further elucidates computational XANES simulations and techniques potentially guiding characterization in surface organometallic chemistry.

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

DFT Analysis of the Binding of Rare Earth Nitrates at Internal and External Surfaces of MCM-22

MCM-22-type zeolites constitute a well-characterized tunable class of aluminosilicates suitable for elucidating the fundamental aspects of the binding of rare earth elements (REEs) in layered materials. Starting from the pure silica version ITQ-1 with a unit cell of Si 72 O 144 , a model for periodic bulk crystalline MCM-22 with a finite Al concentration is provided by replacing a Si atom with an Al atom in the unit cell at suitable tetrahedral sites near an internal pore surface. Then, a H atom is added to an O atom bridging Si and Al atoms to create a Brønsted acid site (BAS). There are no internal silanol groups in this bulk model. To generate a model for an external surface, we adopt the fully hydroxylated surface structure of a layer within the ITQ-1 precursor with two silanols per lateral unit cell. A BAS on the external surface can be generated by replacing a near-surface Si atom with an Al atom and adding a H atom, as above. The strength of binding at a BAS of REE, X, taken to be present in the solution phase as nitrates, is determined from the energy change in the reaction X(NO 3 ) 3 + ≡Si–{OH}–Al≡ → ≡Si–{OX(NO 3 ) 2 }–Al≡+ HNO 3 . The strength of binding at the silanols is determined similarly. Binding energies are determined from two approaches. The first performs periodic plane-wave density functional theory (DFT) total energy analysis for an entire unit cell of MCM-22. The second utilizes cluster models capturing the local environment of REE binding sites and performs DFT analysis with localized basis sets. The two approaches yield consistent results for Nd, revealing similarly strong binding at either an external or internal BAS, but much weaker binding at a silanol site. This is consistent with the picture deduced from recent experiments. Here, we also comment on binding at Al-bridged siloxane sites, which have been suggested as alternative binding sites to BAS.

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

The Ground State Electronic Energy of Benzene

We report here on the findings of a blind challenge devoted to determining the frozen-core, full configuration interaction (FCI) ground-state energy of the benzene molecule in a standard correlation-consistent basis set of double-ζ quality. As a broad international endeavor, our suite of wave function-based correlation methods collectively represents a diverse view of the high-accuracy repertoire offered by modern electronic structure theory. In our assessment, the evaluated high-level methods are all found to qualitatively agree on a final correlation energy, with most methods yielding an estimate of the FCI value around -863 mE H . However, we find the root-mean-square deviation of the energies from the studied methods to be considerable (1.3 mE H ), which in light of the acclaimed performance of each of the methods for smaller molecular systems clearly displays the challenges faced in extending reliable, near-exact correlation methods to larger systems. While the discrepancies exposed by our study thus emphasize the fact that the current state-of-the-art approaches leave room for improvement, we still expect the present assessment to provide a valuable community resource for benchmark and calibration purposes going forward.

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

Unleashed from constrained optimization: quantum computing for quantum chemistry employing generator coordinate inspired method

Hybrid quantum-classical approaches offer potential solutions to quantum chemistry problems, yet they often manifest as constrained optimization problems. Here, we explore the interconnection between constrained optimization and generalized eigenvalue problems through the Unitary Coupled Cluster (UCC) excitation generators. Inspired by the generator coordinate method, we employ these UCC excitation generators to construct non-orthogonal, overcomplete many-body bases, projecting the system Hamiltonian into an effective Hamiltonian, which bypasses issues such as barren plateaus that heuristic numerical minimizers often encountered in standard variational quantum eigensolver (VQE). Diverging from conventional quantum subspace expansion methods, we introduce an adaptive scheme that robustly constructs the many-body basis sets from a pool of the UCC excitation generators. This scheme supports the development of a hierarchical ADAPT quantum-classical strategy, enabling a balanced interplay between subspace expansion and ansatz optimization to address complex, strongly correlated quantum chemical systems cost-effectively, setting the stage for more advanced quantum simulations in chemistry.

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

Regularized machine learning on molecular graph model explains systematic error in DFT enthalpies

Abstract A major goal of materials research is the discovery of novel and efficient heterogeneous catalysts for various chemical processes. In such studies, the candidate catalyst material is modeled using tens to thousands of chemical species and elementary reactions. Density Functional Theory (DFT) is widely used to calculate the thermochemistry of these species which might be surface species or gas-phase molecules. The use of an approximate exchange correlation functional in the DFT framework introduces an important source of error in such models. This is especially true in the calculation of gas phase molecules whose thermochemistry is calculated using the same planewave basis set as the rest of the surface mechanism. Unfortunately, the nature and magnitude of these errors is unknown for most practical molecules. Here, we investigate the error in the enthalpy of formation for 1676 gaseous species using two different DFT levels of theory and the ‘ground truth values’ obtained from the NIST database. We featurize molecules using graph theory. We use a regularized algorithm to discover a sparse model of the error and identify important molecular fragments that drive this error. The model is robust to rigorous statistical tests and is used to correct DFT thermochemistry, achieving more than an order of magnitude improvement.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Computing molecular excited states on a D-Wave quantum annealer

Abstract The possibility of using quantum computers for electronic structure calculations has opened up a promising avenue for computational chemistry. Towards this direction, numerous algorithmic advances have been made in the last five years. The potential of quantum annealers, which are the prototypes of adiabatic quantum computers, is yet to be fully explored. In this work, we demonstrate the use of a D-Wave quantum annealer for the calculation of excited electronic states of molecular systems. These simulations play an important role in a number of areas, such as photovoltaics, semiconductor technology and nanoscience. The excited states are treated using two methods, time-dependent Hartree–Fock (TDHF) and time-dependent density-functional theory (TDDFT), both within a commonly used Tamm–Dancoff approximation (TDA). The resulting TDA eigenvalue equations are solved on a D-Wave quantum annealer using the Quantum Annealer Eigensolver (QAE), developed previously. The method is shown to reproduce a typical basis set convergence on the example $$\hbox {H}_2$$ H 2 molecule and is also applied to several other molecular species. Characteristic properties such as transition dipole moments and oscillator strengths are computed as well. Three potential energy profiles for excited states are computed for $$\hbox {NH}_3$$ NH 3 as a function of the molecular geometry. Similar to previous studies, the accuracy of the method is dependent on the accuracy of the intermediate meta-heuristic software called qbsolv.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

The water trimer reaction OH + (H 2 O) 3 → (H 2 O) 2 OH + H 2 O

All important stationary points on the potential energy surface (PES) for the reaction OH + (H 2 O) 3 → (H 2 O) 2 OH + H 2 O have been fully optimized using the “gold standard” CCSD(T) method with the large Dunning correlation-consistent cc-pVQZ basis sets. Three types of pathways were found. For the pathway without hydrogen abstraction, the barrier height of the transition state (TS1) is predicted to lie 5.9 kcal mol -1 below the reactants. The two major complexes (H 2 O) 3 ∙∙∙OH (CP1 and CP2a) are found to lie 6.3 and 11.0 kcal mol -1 , respectively, below the reactants [OH + (H 2 O) 3 ]. For one of the H-abstraction pathways the lowest classical barrier height is predicted to be much higher, 6.1 kcal mol -1 (TS2a) above the reactants. For the other H-abstraction pathway the barrier height is even higher, 15.0 (TS3) kcal mol -1 . Vibrational frequencies and the zero-point vibrational energies connected to the PES are also reported. Fianlly, the energy barriers for the H-abstraction pathways are compared with those for the OH + (H 2 O) 2 and OH + H 2 O reactions, and the effects of the third water on the energetics are usually minor (0.2 kcal mol -1 ).

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

The many-body expansion for aqueous systems revisited: III. Hofmeister ion – water interactions

We report a Many Body Energy (MBE) analysis of aqueous ionic clusters containing anions and cations at the two opposite ends of the Hofmeister series, viz. the kosmotropes Ca2+, SO42- and chaotropes NH4+ and ClO4- with 9 water molecules to quantify the how these ions in altering the interaction between the water molecules in their immediate surrounding. The current results are contrasted to the ones reported earlier for water clusters as well as for alkali metal and halide ion aqueous clusters of the same size, which lie in the middle of the Hofmeister series. Through this analysis, noteworthy differences between the MBE of kosmotropes and chaotropes were identified. The MBE of kosmotropes is dominated by ion-water interactions that extends beyond the 4-body term, the point at which the MBE of pure water converges. The percentage contribution of the 2- B to the total cluster binding energy is noticeably larger. The disruption due to the dominant ion results in weak, unfavorable water-water interactions. The MBE for chaotropes, on the other hand, was found to converge more quickly as it more closely resembles that of pure water clusters. Chaotropes exhibit weaker overall binding energies and ion-water interactions with more favorable water-water interactions, somewhat recovering the pattern of the 2-4 body terms exemplified by pure water clusters. More importantly, both kosmotropic and chaotropic ions exhibit an anticorrelation between the 2-B ion-water (I-W) and water-water (W-W) interactions as well as between the 3-B (I-W-W) and (I-W) interactions. The consideration of two different structural arrangements (ion inside and outside of a water cluster) suggests that fully solvated (ion inside) chaotropes disrupt the hydrogen bonding network in a similar manner as partially solvated (ion outside) kosmotropes and offer useful insights into the modeling requirements of bulk vs. an interface. Finally, the 2-B contribution to the total Basis Set Superposition Error (BSSE) correction for the kosmotropic and chaotropic ions follows the previously reported erf profile vs. intermolecular distance. When scaled for the corresponding dimer energies and distances, a single profile fits the current results together with all previously reported ones for the pure water and halide water clusters.

Herman, Kristina M.↗
DOE OSTI2022

Diffusion Monte Carlo evaluation of disiloxane linearisation barrier

The disiloxane molecule is a prime example of silicate compounds containing the Si-O-Si bridge. The molecule is of significant interest within the field of quantum chemistry, owing to the difficulty in theoretically predicting its properties. Herein, the linearisation barrier of disiloxane is investigated using a fixed-node diffusion Monte Carlo (FNDMC) approach, which is one of the most reliable ab initio methods in accounting for the electronic correlation. Calculations utilizing the density functional theory (DFT) and the coupled cluster method with single and double substitutions, including noniterative triples (CCSD(T)) are carried out alongside FNDMC for comparison. It is concluded that FNDMC successfully predicts the disiloxane linearisation barrier and does not depend on the completeness of the basis-set as much as DFT or CCSD(T), thus establishing its suitability.

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

Potential energy profile for the Cl + (H 2 O) 3 → HCl + (H 2 O) 2 OH reaction. A CCSD(T) study

Four different reaction pathways are initially located for the reaction of Cl atom plus water trimer Cl + (H 2 O) 3 → HCl + (H 2 O) 2 OH using a standard DFT method. As found for the analogous fluorine reaction, the geometrical and energetic results for the four chlorine pathways are closely related. However, the energetics for the Cl reaction are very different from those for fluorine. Here in this paper, we investigate the lowest-energy chlorine pathway using the “gold standard” CCSD(T) method in conjunction with correlation-consistent basis sets up to cc-pVQZ. Structurally, the stationary points for the water trimer reaction Cl + (H 2 O) 3 may be compared to those for the water monomer reaction Cl + H 2 O and water dimer reaction Cl + (H 2 O) 2 . Based on the CCSD(T) energies, the title reaction is endothermic by 19.3 kcal mol -1 , with a classical barrier height of 16.7 kcal mol -1 between the reactants and the exit complex. There is no barrier for the reverse reaction. The Cl … (H 2 O) 3 entrance complex lies 5.3 kcal mol -1 below the separated reactants. The HCl … (H 2 O) 2 OH exit complex is bound by 8.6 kcal mol -1 relative to the separated products. The Cl + (H 2 O) 3 reaction is somewhat similar to the analogous Cl + (H 2 O) 2 reaction, but qualitatively different from the Cl + H 2 O reaction. It is reasonable to expect that the reactions between the chlorine atom and larger water clusters may be similar to the Cl + (H 2 O) 3 reaction. The potential energy profile for the Cl + (H 2 O) 3 reaction is radically different from that for the valence isoelectronic F + (H 2 O) 3 system, which may be related to the different bond energies between HCl and HF.

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

Accurate core excitation and ionization energies from a state-specific coupled-cluster singles and doubles approach

Here, we investigate the use of orbital-optimized references in conjunction with single-reference coupled-cluster theory with single and double substitutions (CCSD) for the study of core excitations and ionizations of 18 small organic molecules, without the use of response theory or equation-of-motion (EOM) formalisms. Three schemes are employed to successfully address the convergence difficulties associated with the coupled-cluster equations, and the spin contamination resulting from the use of a spin symmetry-broken reference, in the case of excitations. In order to gauge the inherent potential of the methods studied, an effort is made to provide reasonable basis set limit estimates for the transition energies. Overall, we find that the two best-performing schemes studied here for ΔCCSD are capable of predicting excitation and ionization energies with errors comparable to experimental accuracies. The proposed ΔCCSD schemes reduces statistical errors against experimental excitation energies by more than a factor of two when compared to the frozen-core core–valence separated (FC-CVS) EOM-CCSD approach – a successful variant of EOM-CCSD tailored towards core excitations.

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

An extensive assessment of the performance of pairwise and many-body interaction potentials in reproducing ab initio benchmark binding energies for water clusters n = 2–25

We assess the performance of 7 pairwise additive (TIP3P, TIP4P, TIP4P-ice, TIP5P, OPC, SPC, SPC/E) and 8 families of many-body potentials (q-AQUA, HIPPO, AMOEBA, EFP, TTM, WHBB, MB-pol, MB-UCB) in reproducing high-level ab initio benchmark values, CCSD(T) or MP2 at the complete basis set (CBS) limit for the binding energy and the many-body expansion (MBE) of water clusters n = 2–11, 16–17, 20, 25. By including a large range of cluster sizes having dissimilar hydrogen bonding networks, we obtain an understanding of how these potentials perform for different hydrogen bonding arrangements that are mostly outside of their parameterization range. While it is appropriate to compare the results of ab initio based many-body potentials directly to the electronic binding energies (D e 's), the pairwise additive ones are compared to the enthalpies at T = 298 K, ΔH(298 K), as the latter class of force fields are parametrized to reproduce enthalpies (implicitly accounting for zero-point energy corrections) rather than binding energies. We find that all pairwise additive potentials considered overestimate the reference ΔH values for the n = 2–25 clusters by >13%. For the water dimer (n = 2) in particular, the errors are in the range 83–119% for the pairwise additive potentials studied since these are based on an effective rather than the true 2-body interaction specifically designed as a means of partially accounting for the missing many-body terms. This stronger 2-body interaction is achieved by an enhanced monomer dipole moment that mimics its increase from the gas phase monomer to the condensed phase value. Indeed, for cluster sizes n ≥ 4 the percent deviations become slightly smaller (albeit all exceeding 13%). In contrast, we find that the many-body potentials perform more accurately in reproducing the electronic binding energies (D e 's) throughout the entire cluster range (n = 2–25), all reproducing the ab initio benchmark binding energies within ±7% of the respective CBS values. We further assess the ability of a subset of the many-body potentials (MB-UCB, q-AQUA, MB-pol, and TTM2.1-F) to also reproduce the magnitude of the ab initio many-body energy terms for water cluster sizes n = 7, 10, 16 and 17. The potentials show an overall good agreement with the available benchmark values. However, we identify characteristic differences upon comparing the many-body terms at both the ab initio-optimized geometries and the respective potential-optimized geometries to the reference ab initio values. Additionally, by applying this analysis to a wide range of cluster sizes, trends in the MBE of the potentials with increasing cluster size can be identified. Finally, in an attempt to draw a parallel between the pairwise additive and many-body potentials, we report the analysis of the individual molecular dipole moments for water clusters with 1 to ~4 solvation shells with the TTM2.1-F potential. We find that the internally solvated water molecules have in general a larger molecular dipole moment ranging from 2.6–3.0 D. This justifies the use of an enhanced, with respect to the gas-phase value, molecular dipole moment for the pairwise additive potentials, which is intended to fold in the many body terms into an effective (enhanced) pairwise interaction through the choice of the charges. These results have important implications for the development of future generations of efficient, transferable, and highly accurate classical interaction potentials in both the pairwise additive and many-body categories.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗