Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “INTERACTIONS”

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,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Collisional broadening of nonlinear resonant wave–particle interactions

A general procedure for understanding plasma behaviour when resonant wave–particle interactions are the sole destabilizing and transport mechanism or only heating and/or current drive source is highlighted without recourse to involved numerical or analytical treatments. These phenomena are characterized by transport that appears to be collisionless even though collisions play a central role in narrow collisional boundary layers. The order of magnitude estimates, which include nonlinear effects, are shown to provide expressions in agreement with the principal results of recent toroidal Alfvén eigenmode (TAE), toroidal magnetic field ripple, and heating and current drive treatments. More importantly, the retention of nonlinearities leads to new estimates of the alpha particle energy diffusivity at saturation for TAE modes, and the ripple threshold at which superbanana plateau evaluations of alpha particle transport are modified by nonlinear radial drift effects. In addition, the estimates indicate when quasilinear descriptions for heating and current drive will begin to fail. The phenomenological procedure demonstrates that in magnetic fusion relevant plasmas, narrow collisional boundary layers must be retained for resonant wave–particle interactions as they enhance the role of collisions, and make stochastic particle motion unlikely to be more important than other nonlinear processes.

Physics↗
DOE OSTI2025

Mechanical Roles of Polysaccharide Assembly and Interactions in Plant Cell Walls

Plants synthesize polysaccharide-based primary cell walls that possess unique microstructures and mechanical properties to accommodate plant growth and provide protection. Here, it remains challenging to assess the role of polysaccharide organization and interactions in the mechanical behavior of primary cell walls owing to their complex microstructure and highly nonlinear mechanical responses. Employing a coarse-grained molecular dynamics model developed for onion epidermal walls, this work explores the conditions under which polysaccharide assembly and interactions might play a significant role in primary cell wall mechanics. Cellulose–cellulose adhesion plays a dominant role in the wall load-bearing capacity, but when cellulose–cellulose adhesion was disrupted computationally, cellulose–xyloglucan adhesion could influence the wall load-bearing capacity. Contrary to the common concept that xyloglucans mechanically tether well-separated cellulose microfibrils, xyloglucans functioned in this case as interfibrillar adhesives capable of transmitting tensile forces between cellulose microfibrils. Our findings may inform design criteria of new materials inspired by plant cell walls.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Clathrate BaNi 2 P 4 : An Interplay of Heat and Charge Transport Due to Strong Host–Guest Interactions

Heat and charge transport properties of the metallic unconventional clathrate BaNi 2 P 4 , hosting Ba cations in oversized Ni 8 P 16 cages, are investigated. A novel method of single-crystal growth was developed, yielding 2–3 mm sized crystals of BaNi 2 P 4 . We also developed a setup to accurately measure thermal conductivity and electrical resistivity of the synthesized single crystals in a wide temperature range avoiding crystal remounting. BaNi 2 P 4 has a metallic temperature dependence of its electrical resistivity (decreasing with decreasing temperature) and manifests an unconventional T 2 power law for 50 K < T < 300 K; below 50 K, the power-law exponent increases gradually such that below 10 K the power law is T 5 , a predicted but extremely rarely experimentally observed dependence for peculiar electron–phonon interactions. Furthermore, electronic band structure calculations, consistent with measurements of de Haas–van Alphen oscillations, show large band dispersions with significant contributions of Ba orbitals to states near the Fermi level, which is atypical for clathrates. The thermal properties of BaNi 2 P 4 were probed using a combination of variable-temperature single-crystal X-ray diffraction experiments, heat capacity measurements, first-principles phonon dispersion calculations, and inelastic neutron scattering measurements. BaNi 2 P 4 exhibits significant hybridization of the Ba-guest and Ni–P-framework vibrational modes, which may be enhanced via the detected split of the Ba position, which results in strong Ba-framework interactions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Potential Skyrmion Host Fe(IO3)3: Connecting Stereoactive Lone-Pair Electron Effects to the Dzyaloshinskii-Moriya Interaction

Magnetic skyrmions, which are topologically distinct magnetic spin textures, are gaining increased attention for their unique physical properties and potential applications in spintronic devices. Here we present a design strategy for skyrmion host candidates based on combinations of magnetic spin, asymmetric building units having stereoactive lone-pair electrons, and polar lattice symmetry. To demonstrate the viability of the proposed rational design principles, we successfully synthesized a Fe(IO 3 ) 3 polycrystalline sample and single crystals by using a new simplified low-temperature pathway, which is experimentally feasible for extending materials growth of transition metal iodates. Single crystal X-ray and powder synchrotron X-ray diffraction measurements demonstrated that Fe(IO 3 ) 3 crystallizes in the polar chiral hexagonal lattice with space group P63. The combined structural features of the macroscopic electric polarization along the c-axis stemming from the coalignment of the stereoactive lone-pairs of the IO 3 – trigonal pyramid and the magnetic Fe 3+ cation residing on the 3-fold rotation axis were selected to promote asymmetric exchange coupling. We find evidence of a predicted skyrmion phase at 14 K ≤ T ≤ 16 K and 2.5 T ≤ μ 0 H ≤ 3.2 T driven by a Dzyaloshinskii–Moriya (DM) interaction, a conclusion supported by the appreciable DM exchange and the zero-field spiral antiferromagnetic ground state of Fe(IO 3 ) 3 deduced from neutron diffraction experiments. The associated magnetic modulation wavelength of the putative skyrmions is expected to be short ~18 nm, comparable to the period of the DM-driven incommensurate order. This work links stereoactive lone-pair electron effects to enhanced DM interaction, demonstrating a new approach for chemical guidelines in the search for skyrmionic states of matter.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Precise Control of Noncovalent Interactions in Semiconducting Polymers for High-Performance Organic Field-Effect Transistors

Performed through side-chain engineering or by incorporating intramolecular locking units, the directionality and dynamic nature of noncovalent interactions are particularly attractive for the design of novel semiconducting materials in a wide variety of applications. This work investigates the nature and position of hydrogen bonding (intra- versus intermolecular), with the objective of developing a rational approach to the design of new semiconducting materials with improved properties in the solid state. To control the polymer chains’ self-assembly, a π-conjugated polymer incorporating a moiety capable of generating intramolecular hydrogen bonding is evaluated against a polymer that allows for intermolecular hydrogen bonding. Characterization through various techniques, optical spectroscopies, grazing incidence wide-angle x-ray scattering, and solution small-angle neutron scattering showed that intramolecular hydrogen bonds resulted in materials with improved crystallinity and higher effective conjugation in the solid state. Additionally, the effect of the noncovalent interaction configuration on the optoelectronic properties was analyzed in organic field-effect transistor fabrication. Devices prepared from the materials with intramolecular hydrogen bonds showed significantly higher performance, with three orders of magnitude higher charge mobility than their counterparts fabricated from polymers with intermolecular hydrogen bonds. These results confirm the importance of chemical design on polymer structures and offer a novel route for the design of high-efficiency semiconducting polymers for next-generation electronics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Expanding the Design Space of Polymer–Metal Organic Framework (MOF) Gels by Understanding Polymer–MOF Interactions

The fabrication of polymer-MOF composite gels holds great potential to provide emergent properties for drug delivery, environmental remediation, and catalysis. To leverage the full potential of these composites, we investigated how the presence and chemistry of polymers impact MOF formation within the composites and, in turn, how MOFs impact polymer gelation. We show that polymers with a high density of strongly metal-binding carboxylic acids inhibit MOF formation; however, reducing the density of carboxylic acids or substituting them with weaker metal-binding hydroxyl groups permits both MOF formation and gelation within composites. Preparing composites with poly(ethylene glycol) (PEG), which does not bind MOF zirconium (Zr)-oxo clusters, and observing gelation suggests that MOFs can entrap polymer chains to create cross-links in addition to cross-linking them through polymer-Zr-oxo interactions. Both simulations and experiments show composite hydrogels formed with poly(vinyl alcohol) (PVA) to be more stable than those made with PEG, which can reptate through MOF pores upon heating. We demonstrate the generalizability of this composite formation process across different Zr-based MOFs (UiO-66, NU-901, UiO-67, and MOF-525) and by spin-coating gels into conformable films. PVA-UiO-66 composite hydrogels demonstrated high sorption and sustained release of methylene blue relative to the polymer alone (3× loading, 28× slower release), and PVA-MOF-525 composite hydrogels capably sorb the therapeutic peptide Angiotensin 1–7. By understanding the influence of polymer-MOF interactions on the structure and properties of composite gels, this work informs and expands the design space of this emerging class of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Non-Covalent Interactions and Helical Packing in Thiophene-Phenylene Copolymers: Tuning Solid-State Ordering and Charge Transport for Organic Field-Effect Transistors

In this study, we introduce two thiophene-phenylene-thiophene (TPT) polymers designed to leverage noncovalent intramolecular interactions to regulate main-chain conformation and enhance solid-state ordering. By incorporating unsubstituted thiophene (T) or bithiophene (2T) units, we reveal striking divergence in the thermal, morphological, and optoelectronic properties of the resulting films, facilitated by these noncovalent interactions. Using a combination of computational and experimental approaches, we show that annealing yields remarkably different polymer conformations and, consequently, charge transport properties. TPT-T undergoes a significant structural transformation, adopting a more planar backbone conformation and a highly crystalline, edge-on molecular orientation. In contrast, the introduction of a single additional thiophene unit in TPT-2T leads to a more isotropic molecular orientation with a slight preference for face-on alignment, resulting in a heterogeneous film structure that hinders charge transport despite achieving tighter molecular packing. Remarkably, despite being composed of achiral components, TPT-2T develops chirality upon annealing, indicating the formation of a helical conformation. Organic field-effect transistor measurements reveal that the well-ordered alignment in annealed TPT-T films results in higher charge carrier mobility and a narrower distribution of mobility values than in TPT-2T. These findings provide critical insights into the structure−property relationships of conjugated polymers, offering guidance for optimizing molecular design and processing strategies for highperformance organic electronic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Reactive Transport Modeling of Shale–Fluid Interactions after Imbibition of Fracturing Fluids

Injection of hydraulic fracturing fluid (HFF) into shale formations for unconventional oil/gas production results in chemical reactions in the shale matrix. Our recent experimental study determined the depths to which different types of reactions between the shale matrix and the HFF extended. In the present study, we built continuum-scale reactive transport models to understand the coupling of chemical reactions and the transport of aqueous species in these shale–HFF systems. Calibration of the model with our previous experimental results reveals that it takes hours to months for the shale matrix to completely neutralize the imbibed acids, depending primarily upon the carbonate content of the shale. Both the HFF pH and pore pH affect the location of barite precipitation, resulting in unique barite precipitation profiles extending millimeters into calcite-rich Eagle Ford shale but only tens of micrometers into low-carbonate Marcellus shale. In addition, dissolved oxygen and extracted bitumen are key to reproducing the experimental observation of Fe(III) (oxyhydr)oxide formation in the shale matrix as a result of pyrite dissolution in the shales. A comparison between the modeling results of porosity in the present study to experimentally measured permeabilities in our previous study suggests that chemical reactions occurring at a greater depth than the observable reaction zone might have impacted permeability. Our model serves as a benchmark for efficiently modeling water–rock interactions in similar systems where bulk rock samples react with a solution in batch reactors. Important reactive transport processes were ascertained via modeling, which allows for quantitative prediction of shale–HFF interactions in shale matrices given the shale and HFF compositions.

58 GEOSCIENCES↗
DOE OSTI2023

The Potential Interactions of Novec 1230 and Its Thermal Degradation Products with Actinide Metals and Oxides

The environmentally friendly and mainly chemically inert fire suppressant, Novec 1230, is being increasingly used worldwide. However, no published study has considered using Novec 1230 with fires involving radioactive material including actinides; here, therefore, this research note focuses on the possible interactions between Novec 1230 and its thermal degradation products (TDPs) with some actinide compounds (e.g., actinide metal and oxides) commonly present in a radiological facility in the event of a fire. Previous studies of halogenated fire suppressants used in actinide metal fires indicate a possibility for Novec 1230 or its TDPs to chemically interact with actinide compounds at elevated temperatures; however, these reactions are highly unlikely to cause either runaway exothermic reactions or gaseous actinide release.

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

Complexes of Platinum Group Metals with a Conformationally Locked Scorpionate in a Metal–Organic Framework: An Unusually Close Apical Interaction of Palladium(II)

Here, we report synthetic strategies for installing platinum group metals (PGMs: Pd, Rh, Ir, and Pt) on a scorpionate-derived linker (TpmC*) within a metal–organic framework (MOF), both by room-temperature postsynthetic metalation and by direct solvothermal synthesis, with a wide range of metal loadings relevant for fundamental studies and catalysis. In-depth studies for the palladium adduct Pd(II)@Zr-TpmC* by density-functional-theory-assisted extended X-ray absorption fine structure spectroscopy reveals that the rigid MOF lattice enforces a close Pd(II)–N apical interaction between the bidentate palladium complex and the third uncoordinated pyrazole arm of the TpmC* ligand (Pd–N apical = 2.501 ± 0.067 Å), an interaction that is wholly avoided in molecular palladium scorpionates.

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

Effect of Redox-Inactive Metal Ion–Nickel(III) Interactions on the Redox Properties and Proton-Coupled Electron Transfer Reactivity

We report Mononuclear nickel(II) and nickel(III) complexes of a bisamidate-bisalkoxide ligand, (NMe 4 ) 2 [Ni- II (HMPAB)] (1) and (NMe 4 )[Ni III (HMPAB)] (2), respectively, have been synthesized and characterized by various spectroscopic techniques including X-ray crystallography. The reaction of redox-inactive metal ions (M n+ = Ca 2+ , Mg 2+ , Zn 2+ , Y 3+ , and Sc 3+ ) with 2 resulted in 2-M n+ adducts, which was assessed by an array of spectroscopic techniques including X-ray absorption spectroscopy (XAS), electron paramagnetic resonance (EPR), and reactivity studies. The X-ray structure of Ca 2+ coordinated to Ni(III) complexes, 2Ca 2+ T, was determined and exhibited an average Ni-Ca distance of 3.1253 angstrom, close to the metal ions' covalent radius. XAS analysis of 2-Ca 2+ and 2-Y 3+ in solution further revealed an additional coordination to Ca and Y in the 2-M n+ adducts with shortened Ni-M distances of 2.15 and 2.11 Å, respectively, implying direct bonding interactions between Ni and Lewis acids (LAs). Such a short interatomic distance between Ni(III) and M is unprecedented and was not observed before. EPR analysis of 2 and 2-M n+ species, moreover, displayed rhombic signals with g av > 2.12 for all complexes, supporting the +III oxidation state of Ni. The Ni III /Ni II redox potential of 2 and 2-M n+ species was determined, and a plot of E 1/ 2 of 2-M n+ versus pK a of [M(H 2 O) n ] m+ exhibited a linear relationship, implying that the Ni III /Ni II potential of 2 can be tuned with different redox-inactive metal ions. Reactivity studies of 2 and 2-M n+ with different 4-X-2,6-ditert-butylphenol (4-XDTBP) and other phenol derivatives were performed, and based on kinetic studies, we propose the involvement of a proton-coupled electron transfer (PCET) pathway. Analysis of the reaction products after the reaction of 2 with 4-OMe-DTBP showed the formation of a Ni(II) complex (1a) where one of the alkoxide arms of the ligand is protonated. A pK a value of 24.2 was estimated for 1a. The reaction of 2-M n+ species was examined with 4-OMe-DTBP, and it was observed that the k 2 values of 2-M n+ species increase by increasing the Lewis acidity of redox-inactive metal ions. However, the obtained k 2 values for 2-M n+ species are much lower compared to the k 2 value for 2. Such a variation of PCET reactivity between 2 and 2-M n+ species may be attributed to the interactions between Ni(III) and LAs. Our findings show the significance of the secondary coordination sphere effect on the PCET reactivity of Ni(III) complexes and furnish important insights into the reaction mechanism involving high-valent nickel species, which are frequently invoked as key intermediates in Ni-mediated enzymatic reactions, solar-fuel catalysis, and biomimetic/synthetic transformation reactions.

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

Influencing Bonding Interactions of the Neptunyl (V, VI) Cations with Electron-donating and -withdrawing Groups

Neptunium makes up the largest percentage of minor actinides found in spent nuclear fuel, yet separations of this element have proven difficult due to its rich redox chemistry. Developing new reprocessing techniques should rely on understanding how to control the Np oxidation state and its interactions with different ligands. Designing new ligands for separations requires understanding how to properly tune a system toward a desired trait through functionalization. Emerging technologies for minor actinide separations focus on ligands containing carboxylate or pyridine functional groups, which are desirable due to their high degree of functionalization. So we use DFT calculations to study the interactions of carboxylate and polypyridine ligands with the neptunyl cation [Np(V/VI)O2] +/2+ . A systematic study is performed by varying the electronic properties of the carboxylate and polypyridine ligands through the inclusion of different electron-withdrawing and electron-donating R groups. We focus on how these groups can affect geometric properties, electronic structure, and bonding characterization as a function of the metal oxidation state and ligand character and discuss how these factors can play a role in neptunium ligand design principles.

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

InteraChem: Virtual Reality Visualizer for Reactive Interactive Molecular Dynamics

Interactive molecular dynamics in virtual reality (IMD-VR) simulations provide a digital molecular playground for students as an alternative or complement to traditional molecular modeling kits or 2D illustrations. Previous IMD-VR studies have used molecular mechanics to enable simulations of macromolecules such as proteins and nanostructures for the classroom setting with considerable success. Here, we present the InteraChem molecular visualizer, intended for reactive IMD-VR simulation using semiempirical and ab initio methods. InteraChem visualizes not only the molecular geometry but also (1) isosurfaces such as molecular orbitals and electrostatic potentials and (2) two-dimensional graphs of time-varying simulation quantities such as kinetic/potential energy, internal coordinates, and user-applied force. Additionally, InteraChem employs speech recognition to facilitate user interaction and introduces a novel “atom happiness” visualization using emojis to indicate the energetic feasibility of a particular bonding arrangement. Furthermore, we include a set of accompanying exercises that we have used to teach chemical reactivity in small molecular systems.

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

Improved Fast Randomized Iteration Approach to Full Configuration Interaction

In this work, we present three modifications to our recently introduced fast randomized iteration method for full configuration interaction (FCI-FRI) and investigate their effects on the method’s performance for Ne, H 2 O, and N 2 . The initiator approximation, originally developed for full configuration interaction quantum Monte Carlo, significantly reduces statistical error in FCI-FRI when few samples are used in compression operations, enabling its application to larger chemical systems. The semistochastic extension, which involves exactly preserving a fixed subset of elements in each compression, improves statistical efficiency in some cases but reduces it in others. We also developed a new approach to sampling excitations that yields consistent improvements in statistical efficiency and reductions in computational cost. We discuss possible strategies based on our findings for improving the performance of stochastic quantum chemistry methods more generally.

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

Two-Component Multireference Restricted Active Space Configuration Interaction for the Computation of L-Edge X-ray Absorption Spectra

X-ray absorption spectroscopy is a powerful probe of local electronic and nuclear structure, providing insights into chemical processes. The theoretical prediction and interpretation of metal L-edge X-ray absorption spectra is complicated by both relativistic effects, including spin-orbit coupling, and by the multiconfigurational nature of the states involved. This work details an exact two-component multireference restricted active space configuration interaction (X2C-MRRASCI) scheme that uses an exact two-component state average complete active space self-consistent-field method (X2C-SA-CASSCF), which includes the spin-orbit coupling in a variational manner, for the accurate description of the electronic structure before using a restricted active space configuration interaction method to describe the core excited states of the X-ray spectrum. Benchmark calculations are presented for a series of iron containing complexes, with results showing key features of the spectrum being reproduced, including ligand to metal charge transfer and shake-up excitations.

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

Molecular Modeling of Surfactant Interaction on Phospholipid Bilayers Mimicking Corneal Epithelium

Surfactants found in consumer products can compromise eye corneal membrane integrity upon accidental exposure. Traditional in vitro and in vivo approaches to evaluate membrane–surfactant interaction pose experimental limitations such as species variability, reproducibility, and most often do not provide the overall picture. These limitations motivate the use of in silico models to study phenomena like cellular disruption assays caused by surfactants at the molecular scale. In this work, coarse-grained molecular dynamics simulations have been employed to investigate how nonionic alcohol ethoxylate (AE) and anionic surfactant alcohol ethoxy sulfate (AES) interact with lipid bilayer liposomes that mimic corneal epithelial cell membranes. The spherical liposome is composed of 1,2-dihexadecanoyl-sn-glycero-3-phosphocholine (DPPC), 1,2-di(9Z-octadecenoyl)-sn-glycero-3-phosphoethanolamine (DOPE), 1,2-di(9Z-octadecenoyl)-sn-glycero-3-phospho-l-serine (DOPS), and cholesterol, resembling the composition of the corneal epithelial cells’ membrane bilayer. The simulation consisted of varying degrees of representative surfactant compositions and two initial types of surfactant configurations within or outside the liposome. Our results reveal that both surfactants induce outer leaflet bulging, agreeing with membrane solubilization models. The more highly ethoxylated surfactant, AE, caused more consistent inner leaflet disruption than AES, resulting in significantly more water permeation and membrane thinning. In addition, both surfactants increase the lateral diffusion of lipids within the membrane layers, with higher ethoxylated AE showing a stronger effect than AES. This study demonstrates how surfactant structure and localization influence bilayer membrane integrity, offering mechanistic insights into the irritation potential, thus guiding the rational design of effective surfactant-based formulations.

Lipids↗
DOE OSTI2020

Modern Approaches to Exact Diagonalization and Selected Configuration Interaction with the Adaptive Sampling CI Method

Recent advances in selected configuration interaction methods have made them competitive with the most accurate techniques available and, hence, creating an increasingly powerful tool for solving quantum Hamiltonians. In this work, we build on recent advances from the adaptive sampling configuration interaction (ASCI) algorithm. We show that a useful paradigm for generating efficient selected CI/exact diagonalization algorithms is driven by fast sorting algorithms, much in the same way iterative diagonalization is based on the paradigm of matrix vector multiplication. We present several new algorithms for all parts of performing a selected CI, which includes new ASCI search, dynamic bit masking, fast orbital rotations, fast diagonal matrix elements, and residue arrays. The ASCI search algorithm can be used in several different modes, which includes an integral driven search and a coefficient driven search. The algorithms presented here are fast and scalable, and we find that because they are built on fast sorting algorithms they are more efficient than all other approaches we considered. After introducing these techniques, we present ASCI results applied to a large range of systems and basis sets to demonstrate the types of simulations that can be practically treated at the full-CI level with modern methods and hardware, presenting double- and triple-ζ benchmark data for the G1 data set. The largest of these calculations is Si$_2$H$_6$ which is a simulation of 34 electrons in 152 orbitals. We also present some preliminary results for fast deterministic perturbation theory simulations that use hash functions to maintain high efficiency for treating large basis sets.

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

The many-body expansion for aqueous systems revisited: 1. Water-water interactions

We revisit the Many-Body Expansion (MBE) for water-water interactions by examining the effects of the basis set, including those resulting from the Basis Set Superposition Error (BSSE) correction, on the various terms for selected sizes of water clusters up to n = 21. The analysis is performed at the second order Møller-Plesset (MP2) perturbation theory with the family of augmented correlation consistent basis sets up to five zeta quality for the (H2O)n, n = 7, 10, 13, 16 and 21 clusters, for which we report either the complete MBE (n = 7, 10) or the one through the 6-body (n = 13) and the 5-body terms (n = 16, 21). Our results suggest that any sizeable contributions to the total cluster binding energy arising from the 5-body and larger terms are solely an artifact of the size of the basis set. Indeed, all terms above the 4-body converge monotonically to practically zero at the Complete Basis Set (CBS) limit. The magnitudes of the 5-body and above terms are accurately reproduced even with the smaller basis set of the series (aug-cc-pVDZ) once the BSSE correction is taken into account. The same level of theory (MP2/aug-cc-pVDZ, BSSEcorrected) also accurately reproduces the MP2/CBS values of the 3- and 4-body terms. The contribution of electron correlation to the 3- and 4-body terms is quite small, so that neglecting the correlation contribution in all terms above the 3-body results in an error of ??(0.1%). The BSSE correction to the largest 2-body term in the MBE was accurately estimated from the function ?? [1 + erf(-?? ??)], which is proportional to the common (overlapping) area between two Gaussian distributions whose centers are separated by R, with the constants ??, ?? fitted to the calculated BSSE corrections for the individual 2-body terms of the clusters with each basis set. Our results unquestionably demonstrate that the MBE for water-water interactions converges to the 4-body term, since any finite terms above the 4-body are artifacts of the size of the basis set and can be safely truncated at the 4-body term when either a very large basis set is used or BSSE corrections are taken into account even with the aug-cc-pVDZ basis set. We expect these findings to have important consequences in the pursuit of accurate many-body molecular dynamics simulations for aqueous systems.

Heindel, Joseph↗