Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density gradient function”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Small-scale density irregularities in the nightside Venus ionosphere - Comparison of theory and observations

The theory of the lower-hybrid-drift instability is compared with observations of plasma density irregularities. Marginal stability boundaries (gamma = 0) for the lower-hybrid-drift instability and the occurrence of small-scale density fluctuations are presented as a function of magnetic field B and density n. For plasma density gradient scale lengths in the range 2-10 km, 80-85 percent of the density fluctuations lie in the unstable B/n parameter regime (gamma greater than 0). Stability boundaries for the onset of instability at a wavelength for which the Pioneer Venus Orbiter could measure a Doppler-shifted frequency of 100 Hz are presented. The conditions on B and n for instability in this situation are more stringent than those for marginal stability, especially at low densities (n less than or equal to 5000/cu cm). In general, the instability is most likely to be active in regions of low beta: high magnetic field strength and low density, as found in ionospheric holes or troughs.

Huba, J. D.↗
NASA NTRS1995

Properties of the solar wind electrons between 1 and 3.3 AU from Ulysses thermal noise measurements

In order to describe the distribution function f(v) of the solar wind electrons, the simplest model which is commonly used consists of the sum of two Maxwellians representing two distinct populations: a core (density n(sub c), temperature T(sub c)) and a halo (density n(sub h), temperature T(sub h)). It is possible, with the latter assumptions on the electron f(v), to determine the quasi-thermal noise (QTN) induced on an antenna by the motion of the ambient electrons in the solar wind. Using this distribution and the spectroscopy of thermal noise measurements from the radio receiver on Ulysses in the ecliptic plane, we deduce the total electron density N(sub e), the core temperature T(sub c), and the core and halo kinetic pressures N(sub c)T(sub c) and N(sub h)T(sub h). From these electron parameters, we can define a 'global' electron temperature as T(sub e) = (N(sub c)T(sub c) + N(sub h)T(sub h))/N(sub e). Here we present different radial gradients of T(sub e), between 1 and 3.3 AU, as a function of three classes of N(sub e) at 1 AU: low, intermediate, and high densities. In general all these gradients are found to be positive with different polytrope power law indexes between N(sub e) and T(sub e), which are in general lower than unity. We also show different behaviors of the ratio N(sub h)T(sub h)/N(sub c)T(sub c) for each density class considered. Some possible interpretations for these observations are discussed.

Maksimovic, M.↗
DOE OSTI2022

Ab initio investigation of elastic properties of dilute Cu alloys for high-gradient accelerating structures

The elastic properties and solid solution strengthening of several Cu alloys that may be promising for high gradient applications were studied using ab initio calculations based on density functional theory. Specifically, the temperature and solute concentration dependence of elastic properties including bulk and shear moduli, thermal expansion coefficient, and heat capacity of these alloys were calculated. We found that the bulk and shear moduli are linearly dependent on the solute concentration in the dilute regime and change significantly for different Cu alloys. The thermal expansion coefficient and heat capacity are comparatively less sensitive to solute concentrations. Solid solution strengthening is observed to be linearly correlated with the change of lattice parameter of Cu. This work provides an in-depth understanding of the elastic properties of Cu alloys, which is key for predicting their performance at high gradient in normal conducting accelerating structures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A Meta-Generalized Gradient Approximation for the Cavity-Dependent Exchange-Correlation Interaction in Strongly Coupled Light–Matter Systems

Strong light–matter coupling in optical cavities enables the manipulation of chemical and physical properties without altering molecular composition. Theoretical modeling of such phenomena requires exchange-correlation (XC) functionals that account for both electron–electron and electron–photon (ep) interactions within quantum electrodynamical density functional theory (QEDFT). In this work, we develop a meta-generalized gradient approximation (meta-GGA) specifically targeting the cavity-dependent XC interaction in strongly coupled light–matter systems. This novel approximation is built upon a new semilocal polarizability approximation, which draws from the jellium-with-a-gap model, and can be extended to a “global hybrid” variant that goes beyond the isotropic model from previous approximations. The polarizability model yields significantly improved dispersion coefficients and benchmark calculations with the cavity-dependent XC functional demonstrate improved agreement with QED Hartree–Fock (QED-HF) reference energies. Application to the regioselectivity of brominated nitrobenzene intermediates reveals the functional’s capacity to capture cavity-induced energetic shifts. In conclusion, our results advance the Jacob’s ladder of functionals for QEDFT and provide a practical tool for modeling polaritonic chemistry.

Approximation↗
DOE OSTI2022

Understanding the fundamental driver of semiconductor radiation tolerance with experiment and theory

Space is the operating environment of a multitude of systems that our society is heavily reliant on today, however, maintaining operability there necessitates special consideration of the electronic systems' tolerance of space radiation. Electronic systems are critically dependent on the electronic properties of their semiconductor components, which are modified by space radiation with an adverse impact on the space system performance. What innate property allows some semiconductors to sustain little damage while others accumulate defects rapidly with dose is poorly understood, which limits the extent to which radiation tolerance can be implemented as a design criterion. Here, to gain insight into what properties are drivers of semiconductor radiation tolerance, the first step is to generate a dataset of the relative radiation tolerance of a broad sampling of semiconductors. To accomplish this, Rutherford backscatter channeling experiments are used to compare the displaced lattice atom buildup in InAs, InP, GaP, GaN, ZnO, MgO, and Si as a function of stepwise alpha particle dose. With this experimental information on radiation-induced incorporation of interstitial defects in hand, hybrid density functional theory electron densities (and their derived quantities) are calculated and their gradient and Laplacian are evaluated to obtain key fundamental information about the interactions in each material. It is shown that simple, undifferentiated values (which are typically used to describe bond strength) are insufficient to predict radiation tolerance. Instead, the curvature of the electron density at bond critical points provides a measure of radiation tolerance consistent with the experimental results obtained. This curvature and associated forces surrounding bond critical points have the potential to disfavor the localization of displaced lattice atoms at these points, favoring their diffusion toward perfect lattice positions. With this criterion to predict radiation tolerance, simple density functional theory simulations can be conducted on potential new materials to gain insight into how they may operate in demanding high radiation environments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Slow ferromagnetic fluctuations in the kagome metal Sc 3 ⁢ Mn 3 ⁢Al 7 ⁢ Si 5 revealed by 27 Al NMR

Static and dynamical magnetic and electronic properties of the kagome metal Sc 3 ⁢Mn 3 ⁢ Al 7⁢ Si 5 have been investigated by 27 Al nuclear magnetic resonance (NMR) measurements. Two distinct 27 Al -NMR signals with two different values of quadrupolar frequencies of 𝜈 Q = 1.55(2) and 1.07(2) MHz are observed, which are assigned to Al(1) and Al(2), respectively. From the detailed NMR spectrum measurements under three different magnetic field directions and the density functional theory calculations, the principal axes of the electric field gradient for each Al site have been determined. The temperature dependence of Knight shift (𝐾) shows a similar temperature dependence of the DC magnetic susceptibility 𝜒 except for the low-temperature region below ∼50 K where 𝐾 is almost constant while 𝜒 keeps increasing, which suggests that the increase in 𝜒 at low temperatures is not intrinsic. 27 Al spin-lattice relaxation rate divided by temperature (1/𝑇 1 ⁢𝑇) is found to be constant, confirming the metallic state of Sc 3 ⁢Mn 3 ⁢ Al 7⁢ Si 5 from a microscopic point of view. Based on a Korringa ratio analysis using the 𝑇 1 and 𝐾 data, ferromagnetic fluctuations are found to dominate in Sc 3 ⁢Mn 3 ⁢ Al 7⁢ Si 5 . In conclusion, these fluctuations are suggested to be very slow with frequencies on the order of kilohertz or lower.

Ding, Qing-Ping [Ames Laboratory (AMES), Ames, IA ↗
NASA NTRS1985

Vapor Growth of Alloy-Type Semiconductor Crystals

The present effort is part of a continuing research program directed towards the investigation of basic vapor transport phenomena and of crystal growth properties of electronic materials. The primary purpose of ground-based studies is the development and definition of optimum experimental parameters for flight experiments. The ground-based effort includes the investigation of gravity-driven convection effects on mass transport rates and on crystal morphology for different orientations of the density gradient with respect to the gravity vector, and as a function of pressure and of temperature. In addition to the experimental tasks, theoretical efforts involve the quantitative thermodynamic analysis of the systems under investigation, the computation of fluid dynamic parameters, and the consideration of other possible effects on fluid flow under vertical, stabilizing and microgravity conditions. The specific experiments to be performed in a microgravity environment include the investigation of vapor transport and crystal growth phenomena of the GeSe-Xenon system and of the mass flux and growth of bulk and layer-type crystals of Hg sub (1-x) Cd sub x Te using HgI as a transport agent.

Wiedemeier, H.↗
DOE OSTI2024

Vacancy-mediated transport and segregation tendencies of solutes in fcc nickel under diffusional creep: A density functional theory study

The Nabarro-Herring (NH) diffusional creep theory postulates the vacancy-mediated transport of atoms under a stress gradient as the creep mechanism under low-stress and high-temperature conditions. In multicomponent alloys, we premise that this stress-assisted flow of vacancies to and from grain boundaries will produce elemental segregation. An observation of such segregation, validated with theoretical predictions, can provide the necessary experimental evidence for the occurrence of NH creep. Theoretical calculations of the segregation tendencies via analyzing the dominant solute diffusion mechanisms and the difference in diffusivities of the elements are therefore essential. To this end, this study applies density functional theory calculations of migration barriers and solute-vacancy binding energies as input to the self-consistent mean-field theory to assess the vacancy-mediated diffusion mechanisms, transport coefficients, and segregation tendencies of Co, Cr, Mo, Re, Ta, and W solutes in face-centered-cubic Ni. We find Co, Re, and W to be slow diffusers at high temperatures and Cr, Mo, and Ta to be fast diffusers. Further analysis shows that the slow diffusers tend to always enrich at vacancy sinks over a wide range of temperatures. In contrast, the fast diffusers show a transition from depletion to enrichment as the temperature lowers. Furthermore, our analysis of the segregation tendencies under tensile hydrostatic strains shows that slow diffusers are largely unaffected by the strain and favor enrichment. On the other hand, the fast diffusers exhibit high sensitivity to strain and their segregation tendency can transition from depletion to enrichment at a given temperature. Finally, the transport coefficients calculated in this work are expected to serve as input to mesoscale microstructure models to provide a more rigorous assessment of solute segregation under NH creep conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Symmetry Breaking with the SCAN Density Functional Describes Strong Correlation in the Singlet Carbon Dimer

The SCAN (strongly constrained and appropriately normed) meta-generalized gradient approximation (meta-GGA), which satisfies all 17 exact constraints that a meta-GGA can satisfy, accurately describes equilibrium bonds that are normally correlated. With symmetry breaking, it also accurately describes some sd equilibrium bonds that are strongly correlated. While sp equilibrium bonds are nearly always normally correlated, the C 2 singlet ground state is known from correlated wave function theory to be a rare case of strong correlation in an sp equilibrium bond. Earlier work that calculated atomization energies of the molecular sequence B 2 , C 2 , O 2 , and F 2 in the local spin density approximation (LSDA), the Perdew–Burke–Ernzerhof (PBE) GGA, and the SCAN meta-GGA, without symmetry breaking in the molecule, found that only SCAN was accurate enough to reveal an anomalous under-binding for C 2 . Here, this work shows that spin symmetry breaking in singlet C 2 , which involves the appearance of net up- and down-spin densities on opposite sides (not ends) of the bond, corrects that underbinding, with a small SCAN atomization-energy error more like that of the other three molecules, suggesting that symmetry breaking with an advanced density functional might reliably describe strong correlation. This article also discusses some general aspects of symmetry breaking and the insights into strong correlation that symmetry breaking can bring. The normally correlated low-lying triplet excited state has the right vertical excitation energy in SCAN but not in LSDA or PBE, where the triplet is a false ground state. Fractional occupation numbers are found only for the symmetry-unbroken singlet and only in LSDA and PBE GGA.

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

Study of self-interaction errors in density functional predictions of dipole polarizabilities and ionization energies of water clusters using Perdew-Zunger and locally scaled self-interaction corrected methods

Here, we studied the effects of self-interaction error (SIE) on the static dipole polarizabilities of water clusters modelled with three increasingly sophisticated, non-empirical density functional approximations (DFAs), viz. the local spin density approximation (LDA), the Perdew-Burke-Ernzherof (PBE) generalized-gradient approximation (GGA), and the strongly constrained and appropriately normed (SCAN) meta-GGA, using the Perdew-Zunger self-interaction-correction (PZ-SIC) energy functional in the Fermi-Lowdin orbital SIC (FLO-SIC) framework. Our results show that while all three DFAs overestimate the cluster polarizabilities, the description systematically improves from LDA to PBE to SCAN. The self-correlation free SCAN predicts polarizabilities quite accurately with a mean absolute error (MAE) of 0.58 Bohr 3 with respect to coupled cluster singles and doubles (CCSD) values. Removing SIE using PZ-SIC correctly reduces the DFA polarizabilities, but over-corrects, resulting in underestimated polarizabilities in SIC-LDA, -PBE, and -SCAN. Finally, we applied a recently proposed local-scaling SIC (LSIC) method using a quasi self-consistent scheme and using the kinetic energy density ratio as an iso-orbital indicator. The results show that the LSIC polarizabilities are in excellent agreement with mean absolute error of 0.08 Bohr 3 for LSIC-LDA and 0.06 Bohr 3 for LSIC-PBE with most recent CCSD polarizabilities. Likewise, the ionization energy estimates as an absolute of highest occupied energy eigenvalue predicted by LSIC are also in excellent agreement with CCSD(T) ionization energies with MAE of 0.4 eV for LSIC-LDA and 0.04 eV for LSIC-PBE. The LSIC-LDA predictions of ionization energies are comparable to the reported GW ionization energies while the LSIC-PBE ionization energies are more accurate than reported GW results.

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

A data-driven method for modelling dissipation rates in stratified turbulence

We present a deep probabilistic convolutional neural network (PCNN) model for predicting local values of small-scale mixing properties in stratified turbulent flows, namely the dissipation rates of turbulent kinetic energy and density variance, $\varepsilon$ and $\chi$ . Inputs to the PCNN are vertical columns of velocity and density gradients, motivated by data typically available from microstructure profilers in the ocean. The architecture is designed to enable the model to capture several characteristic features of stratified turbulence, in particular the dependence of small-scale isotropy on the buoyancy Reynolds number $Re_b:=\varepsilon /(\nu N^2)$ , where $\nu$ is the kinematic viscosity and $N$ is the background buoyancy frequency, the correlation between suitably locally averaged density gradients and turbulence intensity and the importance of capturing the tails of the probability distribution functions of values of dissipation. Empirically modified versions of commonly used isotropic models for $\varepsilon$ and $\chi$ that depend only on vertical derivatives of density and velocity are proposed based on the asymptotic regimes $Re_b\ll 1$ and $Re_b\gg 1$ , and serve as an instructive benchmark for comparison with the data-driven approach. When trained and tested on a simulation of stratified decaying turbulence which accesses a range of turbulent regimes (associated with differing values of $Re_b$ ), the PCNN outperforms assumptions of isotropy significantly as $Re_b$ decreases, and additionally demonstrates improvements over the fitted empirical models. A differential sensitivity analysis of the PCNN facilitates a comparison with the theoretical models and provides a physical interpretation of the features enabling it to make improved predictions.

42 ENGINEERING↗
NASA NTRS1984

Plasma and magnetospheric research

Research activities on the following topics were summarized: (1) software for the Space Plasma computer Analysis Network (SPAN), (2) plasmaspheric field-aligned temperature gradients, (3) the shift in spacecraft potential as a function of plasma density, (4) plasma flow, (5) the Fabry-Perot interferometer, and (6) the Differential Ion Flux Probe (DIFP).

Comfort, R. H.↗
NASA NTRS1998

Seasonally Frozen Soil Monitoring Using Passive Microwave Satellite Data and Simulation Modeling

Satellite data and simulation modeling were used to assess seasonally frozen soils in the central US - Canada borders area (46-53 degrees N and 96-108 degrees). We used Scanning Multichannel Microwave Radiometer (SMMR) satellite data to delineate the top layer of frozen soils. SMMR is a passive microwave sensor having five channels (6.6, 10, 18, 21 and 37 GHz) with a horizontal and vertical polarization. SMRR data are available between 1978-1987 with noon and midnight overpass and footprint sizes between 25 km and 150 km. SMMR data were processed from resampled 1/4 degree grid cells during fall freeze-up and spring thaw (fall 1985 - spring 1987). The dielectric properties of a target may directly affect the satellite signal. The dielectric value is an order of magnitude smaller for frozen soil water. There are other significant changes to the emitted microwave signal from changes to the surface physical temperature, attenuation of the soil signal from plant water and soil moisture. We further characterized the temporal and spatial dynamic of frozen soils using the FroST (Frozen Soil Temperature) simulation model. The FroST model was used to further predict soil water and ice content, and soil temperature. SMMR results were compared versus 5-cm soil temperature data from available weather stations (14 in Canada and 11 for available months in the US). SMMR data were analyzed as a function of frequency, polarization, polarization difference, and "frequency gradient". In addition, vegetation density, physical temperature and snow depth were also considered. Preliminary analysis of SMMR derived frozen soil/thaw classification using a simple threshold classification indicates a mean overall classification accuracy by season of 85 percent. A sensitivity analysis for different soils with varying amounts of snow was conducted with FroST, which showed that the amount of snow, and the time of snow fall and melt affected the ice and water content, and depth of thaw. These results indicate a potential source of flooding and erosion under conditions when melting snow and spring rains provide a source of infiltrating water.

Toll, D. L.↗
DOE OSTI2024

Self‐Consistent Convolutional Density Functional Approximations: Application to Adsorption at Metal Surfaces

The exchange-correlation (XC) functional in density functional theory is used to approximate multi-electron interactions. A plethora of different functionals are available, but nearly all are based on the hierarchy of inputs commonly referred to as “Jacob's ladder.” This paper introduces an approach to construct XC functionals with inputs from convolutions of arbitrary kernels with the electron density, providing a route to move beyond Jacob's ladder. We derive the variational derivative of these functionals, showing consistency with the generalized gradient approximation (GGA), and provide equations for variational derivatives based on multipole features from convolutional kernels. A proof-of-concept functional, PBEq, which generalizes the PBEα framework with mathematical equation being a spatially-resolved function of the monopole of the electron density, is presented and implemented. It allows a single functional to use different GGAs at different spatial points in a system, while obeying PBE constraints. Analysis of the results underlines the importance of error cancellation and the XC potential in data-driven functional design. After testing on small molecules, bulk metals, and surface catalysts, the results indicate that this approach is a promising route to simultaneously optimize multiple properties of interest.

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

Structure of pressure-gradient-driven current singularity in ideal magnetohydrodynamic equilibrium

Abstract Singular currents typically appear on rational surfaces in non-axisymmetric ideal magnetohydrodynamic (MHD) equilibria with a continuum of nested flux surfaces and a continuous rotational transition. These currents have two components: a surface current (Dirac δ -function in flux surface labeling) that prevents the formation of magnetic islands, and an algebraically divergent Pfirsch–Schlüter current density when a pressure gradient is present across the rational surface. On flux surfaces adjacent to the rational surface, the traditional treatment gives the Pfirsch–Schlüter current density scaling as J ∼ 1 / Δ ι , where Δ ι is the difference of the rotational transform relative to the rational surface. If the distance s between flux surfaces is proportional to Δ ι , the scaling relation J ∼ 1 / Δ ι ∼ 1 / s will lead to a paradox that the Pfirsch–Schlüter current is not integrable. In this work, we investigate this issue by considering the pressure-gradient-driven singular current in the Hahm–Kulsrud–Taylor problem, which is a prototype for singular currents arising from resonant magnetic perturbations. We show that not only the Pfirsch–Schlüter current density but also the diamagnetic current density are divergent as ∼ 1 / Δ ι . However, due to the formation of a Dirac δ -function current sheet at the rational surface, the neighboring flux surfaces are strongly packed with s ∼ ( Δ ι ) 2 . Consequently, the singular current density J ∼ 1 / s , making the total current finite, thus resolving the paradox. Furthermore, the strong packing of flux surfaces causes a steepening of the pressure gradient near the rational surface, with ∇ p ∼ d p / d s ∼ 1 / s . In general non-axisymmetric MHD equilibrium, contrary to Grad’s conjecture that the pressure profile is flat around densely distributed rational surfaces, our result suggests a pressure profile that densely steepens around them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Low-pressure inductively coupled plasmas in hydrogen: impact of gas heating on the spatial distribution of atomic hydrogen and vibrationally excited states

Non-equilibrium inductively coupled plasmas (ICPs) operating in hydrogen are of significant interest for applications including large-area materials processing. Increasing control of spatial gas heating, which drives the formation of neutral species density gradients and the rate of gas-temperature-dependent reactions, is critical. In this study, we use 2D fluid-kinetic simulations with the Hybrid Plasma Equipment Model to investigate the spatially resolved production of atomic hydrogen in a low-pressure planar ICP operating in pure hydrogen (10–20 Pa or 0.075–0.15 Torr, 300 W). The reaction set incorporates self-consistent calculation of the spatially resolved gas temperature and 14 vibrationally excited states. We find that the formation of neutral-gas density gradients, which result from spatially non-uniform electrical power deposition at constant pressure, can drive significant variations in the vibrational distribution function and density of atomic hydrogen when gas heating is spatially resolved. This highlights the significance of spatial gas heating on the production of reactive species in relatively high-power-density plasma processing sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Heterogeneity of neutrophils in arterial hypertension

Highlights: • Normal-density and low-density neutrophils from patients with hypertension influx more Na{sup +} compared to normotensive individuals. • Normal-density neutrophils from patients with hypertension generate low levels of reactive oxygen species. • Normal-density and low-density neutrophils from patients with hypertension exhibit lower levels of apoptosis. • Normal-density neutrophils from patients with hypertension maintain an activated and degranulated phenotype. • Low-density neutrophils from patients with hypertension show atypical mitochondrial morphology. Cellular heterogeneity and diversity are recognized to contribute to the functions of neutrophils under homeostatic and pathological conditions. We previously suggested that the chronic inflammatory responses associated with hypertension (HTN) are related to the participation of different subpopulations of neutrophils. Two populations of neutrophils can be obtained by density gradient centrifugation: normal-density neutrophils (NDN) and low-density neutrophils (LDN). However, the lack of standardized functional protocols has limited phenotypic characterization and functional comparisons of LDN and NDN. Based on their capability to incorporate Na{sup +}, maturity and activation stage, we characterized NDN and LDN in blood samples from ten patients with HTN and ten healthy individuals (HI) using flow cytometry. We compared the levels of reactive oxygen species (ROS), generation of neutrophil extracellular traps (NETs) and levels of apoptosis in NDN and LDN. In general, the NDN and LDN subpopulations from patients with HTN exhibited higher levels of sodium influx and ROS, and lower levels of apoptosis than the corresponding NDN and LDN subsets from HI. Transmission electron microscopy revealed NDN and LDN from patients with HTN exhibited alterations to mitochondrial morphology and fewer cytoplasmic granules than the corresponding HI subpopulations. Our results indicate both the NDN and LDN subpopulations enhance the effects of inflammation that contribute to the pathophysiology of HTN. Further detailed studies are required to characterize the events during ontogeny of the myeloid lineage that result in the diverse phenotypic characteristics of each subpopulation of LDN and NDN.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Delocalization error poisons the density-functional many-body expansion

The many-body expansion is a fragment-based approach to large-scale quantum chemistry that partitions a single monolithic calculation into manageable subsystems. This technique is increasingly being used as a basis for fitting classical force fields to electronic structure data, especially for water and aqueous ions, and for machine learning. Here, we show that the many-body expansion based on semilocal density functional theory affords wild oscillations and runaway error accumulation for ion–water interactions, typified by F − (H 2 O) N with N ≳ 15. We attribute these oscillations to self-interaction error in the density-functional approximation. The effect is minor or negligible in small water clusters, explaining why it has not been noticed previously, but grows to catastrophic proportion in clusters that are only moderately larger. This behavior can be counteracted with hybrid functionals but only if the fraction of exact exchange is ≳50%, whereas modern meta-generalized gradient approximations including ωB97X-V, SCAN, and SCAN0 are insufficient to eliminate divergent behavior. Other mitigation strategies including counterpoise correction, density correction (i.e., exchange–correlation functionals evaluated atop Hartree–Fock densities), and dielectric continuum boundary conditions do little to curtail the problematic oscillations. In contrast, energy-based screening to cull unimportant subsystems can successfully forestall divergent behavior. These results suggest that extreme caution is warranted when the many-body expansion is combined with density functional theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗