Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximations and Expansions”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Dimensionless parameters for ballistic performance evaluation of ceramic-faced bicomponent targets against sharp-nosed projectiles

Two dimensionless parameters are proposed to give first-order approximations for the ballistic performance of ceramic-faced bicomponent armor plates impacted by conical-nosed steel projectiles. The current work collates expansive experimental data from Wilkins and Mayseless covering a broad range of material and mechanical properties, target construction, and projectile dimensions, the data is collapsed using the Hugoniot elastic limit of the ceramic strike face as a strength metric. We find that within datasets where experimental conditions are kept constant, a singular cubic curve can be generated, but variations in experimental conditions for the ballistic limit test between datasets can result in differences in the predictive curve coefficients.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Multiorbital Quantum Impurity Solver for General Interactions and Hybridizations

Here we present a numerically exact inchworm Monte Carlo method for equilibrium multiorbital quantum impurity problems with general interactions and hybridizations. We show that the method, originally developed to overcome the dynamical sign problem in certain real-time propagation problems, can also overcome the sign problem as a function of temperature for equilibrium quantum impurity models. This is shown in several cases where the current method of choice, the continuous-time hybridization expansion, fails due to the sign problem. Our method therefore enables simulations of impurity problems as they appear in embedding theories without further approximations, such as the truncation of the hybridization or interaction structure or a discretization of the impurity bath with a set of discrete energy levels, and eliminates a crucial bottleneck in the simulation of ab initio embedding problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Long-Term Dynamics of Atmospheric Sulfur Dioxide in Urban and Rural Regions of China: Urbanization and Policy Impacts

High levels of sulfur dioxide (SO 2 ) due to human activities pose a serious air pollution issue in China, especially in urban agglomerations. However, limited research has investigated the impact of anthropogenic emissions on higher SO 2 concentrations in urban regions compared to rural areas in China. Here, we analyzed the trends in SO 2 concentrations from 1980 to 2021 in China using the Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2) dataset. SO 2 column concentrations from the Copernicus Atmosphere Monitoring Service (CAMS) and the Ozone Monitoring Instrument (OMI) during the years 2007–2021 were also examined for validation and comparison purposes. Eight representative areas, including four urban regions (Pearl River Delta [PRD], Beijing-Tianjin-Hebei [BTH], Yangtze River Delta [YRD], and Sichuan Basin [SCB]) and four rural regions (Northeast Region [NER], Mongolian Region [MR], West Region [WR], and Tibetan Plateau Region [TR]) were selected for the analysis. Overall, a significant but fluctuating increase in SO 2 concentrations over China was observed during 1980–2021. During 1980–1997 and 2000–2010, there was an increase in SO 2 concentration, while during 1997–2000 and 2010–2021, a decreasing trend was observed. The average increase in SO 2 concentration was approximately 16 times higher in urban regions than in the rural background. We also found that SO 2 dynamics were highly associated with expansion of urban areas, population density, and gross domestic product. Nonetheless, since 2007, SO 2 concentrations have exhibited a downward trend, which is mainly attributed to the air pollution policies implemented by the Chinese government. Our findings highlight the need for further studies on the impact of SO 2 on regional climate change in China.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Technical Performance of Refractory Liners for Molten Chloride Salt Thermal Energy Storage Systems

A chloride-based molten-salt system that uses a ternary blend of MgCl2/KCl/NaCl is investigated to provide higher temperature thermal energy storage capability. Despite higher thermal stability, molten chlorides present several unique challenges, including the design of internal refractory-ceramic liners to prevent the corrosion and thermal stress of alloy tank shells. This work discusses issues and potential solutions related to containment of molten chloride salt, specifically the optimization of the refractory material at the molten salt interface (hot face). The down-selected hot face candidate limits permeation of salt through the material and forms a highly stable secondary surface phase in equilibrium with the molten salt. A mortar is created using the corrosion resistant hot face brick. Brick and mortar composite are subjected to mechanical stress/strain analysis, in order to calculate composite material properties and better inform thermomechanical models. The U.S. Department of Energy Generation 3 (DOE Gen3) program seeks to develop higher efficiency CSP plants that can provide cost-competitive, flexible power in the U.S. electric grid. The proposed Gen3 Liquid Pathway CSP plant closely resembles the configuration of current nitrate salt power towers with two-tank storage (Gen2). The differences between Gen2 and Gen3 include the types of compatible materials used in salt storage tank construction. Stainless steel loses strength at Gen3 temperatures, and although nickel superalloys would be capable of withstanding sustained high temperatures, these materials are prohibitively expensive at scale. Uninsulated tank shells also pose a significant risk as common steels are highly susceptible to chemical attack from molten chloride salt. To address these concerns, refractory-ceramic based containment materials are proposed to line the inside of the hot and cold storage tanks. In doing so, stainless or carbon steel shells may be used in construction depending on the level of insulation provided. The composition of the internal liner requires careful consideration to maximize the efficacy of multiple parameters including corrosion resistance, strength at operating temperature, durability, and cost. This is particularly true for the material at the interface with the salt, known as the "hot face", which is responsible for protecting the insulating layers between the tank shell and the hot face brick layer. The molten salt in this system is superheated over 300 °C above its freeze temperature. Therefore, unlike other industrial processes which use refractory-lined vessels, it is not expected that a freeze plane will develop in the hot face. Therefore, the hot face must be designed to withstand chemical corrosion and inhibit permeation of molten salt into the insulating layers. A down selection was performed to identify a hot face candidate best equipped to maintain thermal, mechanical, and chemical integrity when exposed to molten salt over extended periods of time. Long-duration chemical capability experiments were conducted with the down selected hot face refractory fully immersed in molten chloride salt for up to 3000 hours. The average salt penetration does not exceed 100 microns. When extrapolated to 20 and 30 years of continuous exposure, the expected salt penetration depth is approximately 2.0 mm and 2.9 mm, respectively. A magnesium-rich secondary phases develops at the salt/refractory interface. X-ray diffraction identifies the material as forsterite (Mg2SiO4), which is reported to form synthetically in molten chloride salt solutions. These results suggest the selected hot face will adequately inhibit salt permeation. While there is optimism that the hot face brick will inhibit salt permeation, mortar joints are typically the weakest point of a refractory brick liner. From a thermochemical perspective, differing thermal expansion coefficients may result in the mortar and brick to grow independent of each other, creating gaps through which molten salt can penetrate. To address this issue, NREL has developed an in-house mortar composed of the down selected hot face brick that has been shown to be compatible with the salt. Compressive stress/strain analyses have been performed on the brick/ mortar composites to generate stress/strain curves. Modulus of elasticity and Poisson's ratio of the composite may be calculated from the stress/strain curves, in order to provide more representative data to finite element mechanical models for accurate approximation of stress on the tank shell and the amount of thermal expansion expected within the tank liner.

41 EE - Solar Energy Technologies Office (EE-4S)↗
DOE OSTI2021

Assessing the future of global energy-for-water

This study incorporates the energy demands of water abstraction, treatment, distribution, and post-use (wastewater) treatment into the Global Change Assessment Model (GCAM), an integrated human-Earth systems model, and analyzes a range of scenarios that estimate the future evolution of this demand of energy. The study builds on prior research on the historical demands of energy-for-water, fitting these demands of energy into existing nation-level energy balances, where demands for energy-for-water have heretofore been aggregated with other energy use in the commercial and public services, industrial, agricultural, and/or electric power sectors. This development allows for improved projections of future energy demands in general, and in this study, allows for assessment of the energy implications of different scenarios of improvements in water access and water quality that are consistent with the Sustainable Development Goals. In our baseline scenario, energy-for-water approximately doubles from 2010 to 2050, and in the scenario with the greatest level of water access, standards for treatment of wastewater, and irrigation expansion, the consequent energy-for-water demands increase three-fold from 2010 to 2050. Note also that these increases do not consider future scarcity-driven increases in energy requirements per unit of water abstraction, treatment, and distribution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Quantum Joule expansion of one-dimensional systems

Here, we investigate the Joule expansion of nonintegrable quantum systems that contain bosons or spinless fermions in one-dimensional lattices. A barrier initially confines the particles to be in half of the system in a thermal state described by the canonical ensemble and is removed at time t = 0. We investigate the properties of the time-evolved density matrix, the diagonal ensemble density matrix, and the corresponding canonical ensemble density matrix with an effective temperature determined by the total energy conservation using exact diagonalization. The weights for the diagonal ensemble and the canonical ensemble match well for high initial temperatures that correspond to negative effective final temperatures after the expansion. At long times after the barrier is removed, the time-evolved Rényi entropy of subsystems bigger than half can equilibrate to the thermal entropy with exponentially small fluctuations. The time-evolved reduced density matrix at long times can be approximated by a thermal density matrix for small subsystems. Few-body observables, like the momentum distribution function, can be approximated by a thermal expectation of the canonical ensemble with strongly suppressed fluctuations. The negative effective temperatures for finite systems go to non-negative temperatures in the thermodynamic limit for bosons, but there is a true thermodynamic effect for fermions, which is confirmed by finite-temperature density matrix renormalization group calculations. We propose the Joule expansion as a way to dynamically create negative temperature states for fermion systems with repulsive interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Interpolating between small- and large-$\mathscr{g}$ expansions using Bayesian model mixing

Bayesian model mixing (BMM) is a statistical technique that can be used to combine models that are predictive in different input domains into a composite distribution that has improved predictive power over the entire input space. We explore the application of BMM to the mixing of two expansions of a function of a coupling constant $\mathscr{g}$ that are valid at small and large values of $\mathscr{g}$ respectively. This type of problem is quite common in nuclear physics, where physical properties are straightforwardly calculable in strong and weak interaction limits or at low and high densities or momentum transfers, but difficult to calculate in between. Interpolation between these limits is often accomplished by a suitable interpolating function, e.g., Padé approximants, but it is then unclear how to quantify the uncertainty of the interpolant. We address this problem in the simple context of the partition function of zero-dimensional Φ 4 theory, for which the (asymptotic) expansion at small $\mathscr{g}$ and the (convergent) expansion at large $\mathscr{g}$ are both known. We consider three mixing methods: linear mixture BMM, localized bivariate BMM, and localized multivariate BMM with Gaussian processes. Here, we find that employing a Gaussian process in the intermediate region between the two predictive models leads to the best results of the three methods. The methods and validation strategies we present here should be generalizable to other nuclear physics settings.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

A mesoscale 3D model of irradiated concrete informed via a 2.5 U-Net semantic segmentation

The concrete biological shield in light-water reactors is exposed to neutron and gamma irradiation, which deteriorates the concrete’s mechanical properties in the long term. To assess the irradiation-induced damage, predictive mechanical models are developed and used in parallel with the characterization of irradiated concrete samples. Realistic 3D simulation domains can drastically improve a model’s prediction. In this work, we utilized x-ray computed tomography (XCT) data of a concrete specimen to reconstruct its 3D microstructure. The XCT data shows low contrast between the concrete’s aggregates and cement paste, resulting in poor image segmentation when using traditional unsupervised techniques. To address this issue, we developed and trained a 2.5D U-Net model on only 24 pre-labeled XCT layers to segment 651 layers of the XCT data. The overall F1-score of the model is approximately 96%. Then, we created a 3D finite element (FE) mesh based on the stack of segmented images. The FE model contains radiation-induced expansion, damage, and creep. The constitutive equations are adapted to each phase (aggregates and cement paste). Here, we simulated the effects of neutron irradiation in the concrete specimen as well as the specimen’s mechanical response to uniaxial compression. Finally, model validation was performed using experimental data on similar concrete specimens in the literature.

2.5D U-Net↗
DOE OSTI2026

Air Emissions Impacts of Fuel Substitution in the United States Cement Industry

High-heat industrial processes cannot be easily electrified, often rely on comparatively high-emitting solid fuels, and are essential for economic growth. Cement production in particular is critical to supporting construction and the expansion and maintenance of infrastructure systems in countries around the world. Fuels for cement production include petcoke, coal, and waste solvents, and their combustion contributes approximately 40% of the CO 2 emissions from global cement production. Broadening the fuel resources to include agricultural residues and other wastes can reduce emissions and divert waste from landfills. We assess the potential to use 16 different alternative fuels across the 88 integrated cement plants in the United States. The share of alternative fuels in the cement production sector could increase from the current level of 17% to as high as 73%. Fuel substitution can reduce the current fuel-related life-cycle emissions (28.6 ± 5.7 Mton CO 2 /year) by 19–45%. The greatest national-scale human health cobenefits occur when alternative fuels are prioritized at facilities currently using coal and petcoke. NO x emissions outcomes from this strategy are mixed, but SO x and particulate matter emissions reductions drive net human health benefits. These results suggest that targeted fuel substitution in heavy industries can yield outsized emissions benefits.

alternative fuels↗
DOE OSTI2021

Annual Status Report (FY 2020): Performance Assessment for the Environmental Restoration Disposal Facility

DOE O 435.1 and DOE M 435.1-1 require that a determination of continued adequacy of the performance assessment (PA) (CP-60089), composite analysis, and disposal authorization statement (DAS) be made annually, and these guidelines must be used to consider the results of data collection and analysis from research, field studies, and monitoring as well as provide the need to update any radioactive waste management basis documents. Beginning in 1996, the Environmental Restoration Disposal Facility (ERDF) started accepting low-level radioactive, hazardous, and mixed wastes generated during cleanup activities at the Hanford Site. ERDF is composed of a series of cells or disposal areas and can accommodate future design expansions as needed. Currently, there are 10 cells. During this reporting period (fiscal year 2020, which extended from October 1, 2019, through September 30, 2020), approximately 3.39E+04 U.S. tons (3.07E+04 metric tons) of waste was disposed at ERDF. From ERDF inception through September 30, 2020, approximately 18.5 million U.S. tons of waste has been disposed at ERDF, which equates to the consumption of approximately 88% of the disposal volume. As a condition of the DAS, disposal operations within ERDF must be in accordance with the waste acceptance criteria (ERDF-00011) that provide specific radionuclide disposal limits, waste form restrictions, and descriptions of acceptable waste packages in compliance with DOE M 435.1-1 requirements. The ERDF waste acceptance criteria stipulate that waste destined for disposal at ERDF be controlled based on source, physical form, and contaminant concentration and activity levels. There have been no changes to the physical configuration of ERDF or to the waste forms (source, physical form, etc.). No new Unreviewed Disposal Question Screenings or Evaluations have been generated during this reporting period. Therefore, there are no noted impacts to the PA, composite analysis, DAS, or radioactive waste management basis documents resulting from the evaluations and screenings. Sum of fraction analysis shows that the disposed inventory meets both the concentration and inventory threshold requirements. A sum of fractions value is computed for ERDF sensitive radionuclides contributing to the all pathways and air pathway inventory limits. Computed values were 8.85E-02 and 1.78E-01, respectively. The disposed waste inventory remained well under the PA imposed limits. Required monitoring was satisfactorily completed during the fiscal year reporting period (fiscal year 2020). Compliance with performance objectives were met as each of the reported values were well below the established limit. Overall, there are no substantive changes to primary PA assumptions nor changes to the PA analysis conclusion; therefore, compliance with DOE O 435.1 and the DAS is maintained.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Fourier Neural Networks as Function Approximators and Differential Equation Solvers

We present a Fourier neural network (FNN) that can be mapped directly to the Fourier decomposition. The choice of activation and loss function yields results that replicate a Fourier series expansion closely while preserving a straightforward architecture with a single hidden layer. The simplicity of this network architecture facilitates the integration with any other higher-complexity networks, at a data pre- or postprocessing stage. We validate this FNN on naturally periodic smooth functions and on piecewise continuous periodic functions. We showcase the use of this FNN for modeling or solving partial differential equations with periodic boundary conditions. The main advantages of the current approach are the validity of the solution outside the training region, interpretability of the trained model, and simplicity of use.

Fourier decomposition↗
DOE OSTI2021

Calculation and interpretation of classical turning surfaces in solids

Abstract Classical turning surfaces of Kohn–Sham potentials separate classically allowed regions (CARs) from classically forbidden regions (CFRs). They are useful for understanding many chemical properties of molecules but need not exist in solids, where the density never decays to zero. At equilibrium geometries, we find that CFRs are absent in perfect metals, rare in covalent semiconductors at equilibrium, but common in ionic and molecular crystals. In all materials, CFRs appear or grow as the internuclear distances are uniformly expanded. They can also appear at a monovacancy in a metal. Calculations with several approximate density functionals and codes confirm these behaviors. A classical picture of conduction suggests that CARs should be connected in metals, and disconnected in wide-gap insulators, and is confirmed in the limits of extreme compression and expansion. Surprisingly, many semiconductors have no CFR at equilibrium, a key finding for density functional construction. Nonetheless, a strong correlation with insulating behavior can still be inferred. Moreover, equilibrium bond lengths for all cases can be estimated from the bond type and the sum of the classical turning radii of the free atoms or ions.

Chemistry↗
DOE OSTI2021

B-PINNs: Bayesian Physics-informal Neural Networks for Forward and Inverse PDE Problems with Noisy Data

We propose a Bayesian physics-informed neural network (B-PINN) to solve both forward and inverse nonlinear problems described by partial differential equations (PDEs) and noisy data. In this Bayesian framework, the Bayesian neural network (BNN) combined with a PINN for PDEs serves as the prior while the Hamiltonian Monte Carlo (HMC) or the variational inference (VI) could serve as an estimator of the posterior. B-PINNs make use of both physical laws and scattered noisy measurements to provide predictions and quantify the aleatoric uncertainty arising from the noisy data in the Bayesian framework. Compared with PINNs, in addition to uncertainty quantification, B-PINNs obtain more accurate predictions in scenarios with large noise due to their capability of avoiding overfitting. We conduct a systematic comparison between the two different approaches for the B-PINNs posterior estimation (i.e., HMC or VI), along with dropout used for quantifying uncertainty in deep neural networks. Our experiments show that HMC is more suitable than VI with mean field Gaussian approximation for the B-PINNs posterior estimation, while dropout employed in PINNs can hardly provide accurate predictions with reasonable uncertainty. Finally, we replace the BNN in the prior with a truncated Karhunen-Loève (KL) expansion combined with HMC or a deep normalizing flow (DNF) model as posterior estimators. The KL is as accurate as BNN and much faster but this framework cannot be easily extended to high-dimensional problems unlike the BNN based framework.

Non linear PDEs, Noisy data, Bysian physics inform↗
DOE OSTI2023

Physics-based adaptivity of a spectral method for the Vlasov–Poisson equations based on the asymmetrically-weighted Hermite expansion in velocity space

We propose a spectral method for the 1D-1V Vlasov–Poisson system where the discretization in velocity space is based on asymmetrically-weighted Hermite functions, dynamically adapted via a scaling α and shifting u of the velocity variable. Specifically, at each time instant an adaptivity criterion selects new values of α and u based on the numerical solution of the discrete Vlasov–Poisson system obtained at that time step. Once the new values of the Hermite parameters α and u are fixed, the Hermite expansion is updated and the discrete system is further evolved for the next time step. The procedure is applied iteratively over the desired temporal interval. The key aspects of the adaptive algorithm are: the map between approximation spaces associated with different values of the Hermite parameters that preserves total mass, momentum and energy; and the adaptivity criterion to update α and u based on physics considerations relating the Hermite parameters to the average velocity and temperature of each plasma species. For the discretization of the spatial coordinate, we rely on Fourier functions and use the implicit midpoint rule for time stepping. The resulting numerical method possesses intrinsically the property of fluid-kinetic coupling, where the low-order terms of the expansion are akin to the fluid moments of a macroscopic description of the plasma, while kinetic physics is retained by adding more spectral terms. Moreover, the scheme features conservation of total mass, momentum and energy associated in the discrete, for periodic boundary conditions. A set of numerical experiments confirms that the adaptive method outperforms the non-adaptive one in terms of accuracy and stability of the numerical solution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Back-to-back dijet production in DIS at arbitrary Bjorken x: TMD gluon distributions to twist-3 accuracy

We derive the gluon transverse-momentum-dependent (TMD) operator structure of back-to-back\\\\r\\\\nquark–antiquark dijet production in deep inelastic scattering at arbitrary Bjorken-x to twist-3 ac\\\\r\\\\ncuracy. Working at leading order in the strong coupling and in the kinematic regime where the\\\\r\\\\ntransverse momentum imbalance of the jets is much smaller than their individual transverse mo\\\\r\\\\nmenta, we perform a systematic gradient expansion of the quark propagator in a background gluon\\\\r\\\\nfield. This expansion organizes multiple interactions with the target in terms of longitudinal Wilson\\\\r\\\\nlines and gauge-invariant field-strength insertions, yielding a TMD description valid beyond the\\\\r\\\\nstrict high-energy eikonal (x → 0) approximation. We obtain explicit cross sections for longitudi\\\\r\\\\nnally and transversely polarized virtual photons, identifying all contributing gluon TMD operators\\\\r\\\\nup to twist-3, including structures involving F+−, Fij, and three-gluon correlators. The full lon\\\\r\\\\ngitudinal phase eixP+z− associated with Bjorken-x is retained throughout. In the small-x limit,\\\\r\\\\nour results reproduce the known sub-eikonal expressions obtained in the Color Glass Condensate\\\\r\\\\nframework, establishing a direct connection between the general-x TMD expansion and high-energy\\\\r\\\\nfactorization. We further reduce the operator basis using equations of motion, minimizing the num\\\\r\\\\nber of independent nonperturbative matrix elements entering the cross section. This work provides\\\\r\\\\na systematic foundation for extending TMD analyses of dijet production beyond leading twist, es\\\\r\\\\ntablishing a unified operator framework valid at arbitrary Bjorken-x that smoothly interpolates\\\\r\\\\nbetween moderate- and small-x descriptions of gluon TMDs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Dynamic high pressure phase transformation of ZrW 2 O 8

Phase transformations under high strain rates (dynamic compression) are examined in situ on ZrW 2 O 8 , a negative thermal expansion ternary ceramic displaying polymorphism. Amorphization, consistent with prior quasi-static measurements, was observed at a peak pressure of 3.0 GPa under dynamic conditions, which approximate those expected during fabrication. Evidence of partial amorphization was observed at lower pressure (1.8 GPa) that may be kinetically restrained by the short (<~150 ns) time scale of the applied high pressure. The impact of kinetics of pressure-induced amorphization from material fabrication methods is briefly discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗