Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cluster finding”

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
DOE OSTI2023

Constraints on S 8 from a full-scale and full-shape analysis of redshift-space clustering and galaxy–galaxy lensing in BOSS

ABSTRACT We present a novel simulation-based cosmological analysis of galaxy–galaxy lensing and galaxy redshift-space clustering. Compared to analysis methods based on perturbation theory, our simulation-based approach allows us to probe a much wider range of scales, $0.4 \, h^{-1} \, \mathrm{Mpc}$ to $63 \, h^{-1} \, \mathrm{Mpc}$, including highly non-linear scales, and marginalizes over astrophysical effects such as assembly bias. We apply this framework to data from the Baryon Oscillation Spectroscopic Survey LOWZ sample cross-correlated with state-of-the-art gravitational lensing catalogues from the Kilo Degree Survey and the Dark Energy Survey. We show that gravitational lensing and redshift-space clustering when analysed over a large range of scales place tight constraints on the growth-of-structure parameter $S_8 = \sigma _8 \sqrt{\Omega _{\rm m} / 0.3}$. Overall, we infer S8 = 0.792 ± 0.022 when analysing the combination of galaxy–galaxy lensing and projected galaxy clustering and S8 = 0.771 ± 0.027 for galaxy redshift-space clustering. These findings highlight the potential constraining power of full-scale studies over studies analysing only large scales and also showcase the benefits of analysing multiple large-scale structure surveys jointly. Our inferred values for S8 fall below the value inferred from the CMB, S8 = 0.834 ± 0.016. While this difference is not statistically significant by itself, our results mirror other findings in the literature whereby low-redshift large-scale structure probes infer lower values for S8 than the CMB, the so-called S8-tension.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Ab initio multishell valence-space Hamiltonians and the island of inversion

In the shell-model framework, valence-space Hamiltonians connecting multiple major-oscillator shells are of key interest for investigating the physics of neutron-rich nuclei, which have been the subject of intense experimental activity for decades. Here we present an extension of the ab initio valence-space in-medium similarity renormalization group, which allows the derivation of such Hamiltonians nonperturbatively. Starting from initial two- and three-nucleon forces from chiral effective field theory, we then calculate properties of nuclei in the important island-of-inversion region above oxygen, so far unexplored with ab initio methods. In this work, our results in the neon and magnesium isotopes indicate the importance of neutron excitation from the $\textit{sd}$ to $\textit{pf}$ shells and ground states dominated by intruder configurations around $\textit{N}$ = 20, consistent with the conclusions from phenomenological studies. We also benchmark the excitation spectrum of 16 O with coupled-cluster theory, finding generally good agreement, and discuss implications for ground-state energies and charge radii in oxygen and calcium isotopes. Finally we outline the proper procedure for treating the longstanding issue of center-of-mass contamination, and show that with a particular choice of valence space, these spurious states can be removed successfully.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Electromagnetic observables of open-shell nuclei from coupled-cluster theory

Here we develop a new method to describe electromagnetic observables of open-shell nuclei with two nucleons outside a closed shell. This approach combines the equation-of-motion coupled-cluster method for such systems and the Lorentz integral transform technique, expanding the applicability of coupled-cluster theory for electromagnetic observables beyond closed-(sub)shell nuclei. To validate this new approach, we compute the non-energy-weighted dipole sum rule and the dipole polarizability of 16,24 ⁢O in both the closed-(sub)shell and the new equation-of-motion coupled-cluster frameworks, finding agreement within error bars. We then analyze the evolution of the dipole polarizability along the oxygen and calcium isotopic chains. Our predictions agree well with available experimental data and other available theoretical calculations for the closed-(sub)shell 16,22 ⁢O and the open-shell 18 O. In the calcium isotopes, we observe that our dipole polarizability predictions for open-shell nuclei are lower than those of closed-(sub)shell nuclei. We expect that our predictions for 24 O and 54,56 ⁢Ca will motivate future experimental studies at the drip line.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Dark-matter-enhanced probe of relic neutrino clustering

We propose heavy decaying dark matter (DM) as a new probe of the cosmic neutrino background (C⁢𝜈⁢B). Heavy DM, with mass ≳ 10 9 GeV, decaying into neutrinos can be a new source of ultrahigh-energy (UHE) neutrinos. Including this contribution along with the measured astrophysical and predicted cosmogenic neutrino fluxes, we study the scattering of UHE neutrinos with the C⁢𝜈⁢B via standard weak interactions mediated by the 𝑍 boson. We solve the complete neutrino transport equation, taking into account both absorption and reinjection effects, to calculate the expected spectrum of UHE neutrino flux at future neutrino telescopes, such as the IceCube-Gen2 radio. We argue that such observations can be used to probe the C⁢𝜈⁢B properties and, in particular, local C⁢𝜈⁢B clustering. We find that, depending on the absolute neutrino mass and the DM mass and lifetime, a local C⁢𝜈⁢B overdensity ≳ 10 6 can be probed at the IceCube-Gen2 radio within ten years of data taking.

Dark matter↗
DOE OSTI2021

Reversible disorder-order transitions in atomic crystal nucleation

Nucleation in atomic crystallization remains poorly understood, despite advances in classical nucleation theory. The nucleation process has been described to involve a nonclassical mechanism that includes a spontaneous transition from disordered to crystalline states, but a detailed understanding of dynamics requires further investigation. In situ electron microscopy of heterogeneous nucleation of individual gold nanocrystals with millisecond temporal resolution shows that the early stage of atomic crystallization proceeds through dynamic structural fluctuations between disordered and crystalline states, rather than through a single irreversible transition. Our experimental and theoretical analyses support the idea that structural fluctuations originate from size-dependent thermodynamic stability of the two states in atomic clusters. These findings, based on dynamics in a real atomic system, reshape and improve our understanding of nucleation mechanisms in atomic crystallization.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Spectral Peak Enhancement by Combining Trusted Response Elements via Machine Learning (SPECTRE-ML) v0.8.0

SPECTRE-ML (Spectral Peak Enhancement by Combining Trusted Response Elements via Machine Learning) is a machine learning based program for finding optimal clusters of radiation detector segments (i.e., pixels or voxels) in order to improve spectral performance. It provides facilities for pre-processing and analyzing training datasets, running ML algorithms, and evaluating and visualizing outputs. SPECTRE-ML outperforms simpler ad-hoc segmentation methods such as uniform depth clusters, learning detector performance trends such as dead layers, edge effects, and gain shifts. Although extensible to arbitrary highly-segmented spectroscopic radiation detectors, SPECTRE-ML currently focuses on improving spectral performance in highly-segmented CdZnTe (CZT) detectors for International Atomic Energy Agency (IAEA) non-destructive assay (NDA) safeguards tasks.

Vavrek, Jayson↗
DOE OSTI2023

DES-Y3 galaxies & ACT DR4 CMB lensing tomography

We present a measurement of the cross-correlation between the \maglim galaxies selected from the Dark Energy Survey (DES) first three years of observations (Y3) and cosmic microwave background (CMB) lensing from the Atacama Cosmology Telescope (ACT) Data Release 4 (DR4), reconstructed over $\sim 436$ $\sqdeg$ of the sky. Our galaxy sample, which covers $\sim 4143$ $\sqdeg$, is divided into six redshift bins spanning the redshift range of $0.20<z<1.05$. We adopt a blinding procedure until passing all consistency and systematics tests. After imposing scale cuts for the cross-power spectrum measurement, we reject the null hypothesis of no correlation at 9.1$\sigma$. We constrain cosmological parameters from a joint analysis of galaxy and CMB lensing-galaxy power spectra considering a flat \LCDM model, marginalized over 23 astrophysical and systematic nuisance parameters. We find the clustering amplitude $S_8\equiv \sigma_8 (\Omega_m/0.3)^{0.5} = 0.75^{+0.04}_{-0.05}$. In addition, we co nstrain the linear growth of cosmic structure as a function of redshift. Our results are consistent with recent DES Y3 analyses and suggest a preference for a lower $S_8$ compared to results from measurements of CMB anisotropies by the \textit{Planck} satellite, although at a mild level ($< 2 \sigma$) of statistical significance.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Disequilibrium Chemistry in Exoplanet Atmospheres Observed with the Hubble Space Telescope

Literature on the theory of exoplanet atmospheric disequilibrium chemistry is rich, although its observational counterpart has yet to emerge beyond the hints provided by a few targets in dedicated studies. We report results from an uniform data reduction and analysis for a catalog of 62 Hubble Space Telescope exoplanet transit spectra where we assess the atmospheric model preference for disequilibrium chemistry (i.e., water vapor is not the dominant absorption spectral signature) over thermal equilibrium chemistry in a comparative planetology context. Where model preference assessment is possible, we find that disequilibrium occurs in about half of the atmospheres, indicating that disequilibrium processes play an important role in the composition of exoplanet atmospheres. While very hot atmospheres, over 1800 K, prefer equilibrium chemistry, we find a clustering of preference for disequilibrium in the 1200–1800 K temperature range. We suggest that UV-augmented thermochemistry may play a significant role for those atmospheres.

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

Modeling the Spectral Diversity of Quasars in the Sixteenth Data Release from the Sloan Digital Sky Survey

Abstract We present a new approach to capturing the broad diversity of emission-line and continuum properties in quasar spectra. We identify populations of spectrally similar quasars through pixel-level clustering on 12,968 high signal-to-noise ratio (S/N) spectra from the Sloan Digital Sky Survey (SDSS) in the redshift range of 1.57 < z < 2.4. Our clustering analysis finds 396 quasar spectra that are not assigned to any population, 15 misclassified spectra, and 6 quasars with incorrect redshifts. We compress the quasar populations into a library of 684 high-S/N composite spectra, anchored in redshift space by the Mg ii emission line. Principal component analysis on the library results in an eigenspectrum basis spanning 1067–4007 Å. We model independent samples of SDSS quasar spectra with the eigenbasis, allowing for a free redshift parameter. Our models achieve a median reduced χ 2 on non–broad absorption line quasar spectra that is reduced by 8.5% relative to models using the eigenspectra from the SDSS spectroscopic pipeline. A significant contribution to the relative improvement is from the ability to reconstruct the range of emission-line variation. The redshift estimates from our model are consistent with the Mg ii emission-line redshift with an average offset that displays 51.4% less redshift-dependent variation relative to the SDSS eigenspectra. Our method for developing quasar spectra models can improve automated classification and predict the intrinsic spectrum in regions affected by intervening absorbers such as Ly α , C iv , and Mg ii , thus benefiting studies of large-scale structure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

DESI EDR: Calibrating the Tully–Fisher Relationship with the DESI Peculiar Velocity Survey

We calibrate the Tully–Fisher relation (TFR) with data from the DESI Peculiar Velocity (PV) Survey taken during the Survey Validation (SV) period of the DESI galaxy redshift survey. Placing spectroscopic fibers on the centers and major axes of spatially extended spiral galaxies identified in the 2020 Siena Galaxy Atlas using the DESI Legacy Surveys, we measure the rotational velocities at 0.33R 26 for 1155 (1128 + 27 dwarf) spiral galaxies observed during SV. Using 39 spiral galaxies observed in the Coma cluster, we find a slope for the TFR of −8.32 ± 0.15 AB mag in the r band, with a scatter about the TFR of 1.12 ± 0.03 AB mag. We calibrate the zero-point of the TFR using galaxies with independent distances measured using type Ia supernovae (SNe Ia) via the cosmological distance ladder. From the SN Ia distances, we measure a zero-point of $-19.21^{+0.30}_{-0.31}$ AB mag in the r band. We produce a public catalog of the distances to these 1128 spiral galaxies observed during DESI SV as part of the DESI PV Survey with our calibrated TFR. This is, to our knowledge, the first catalog of TFR distances produced with velocities measured at a single point in the disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Dark Energy Survey Year 6 Results: Cosmological Constraints from Galaxy Clustering and Weak Lensing

We present cosmology results combining galaxy clustering and weak gravitational lensing measured in the full six years (Y6) of observations by the Dark Energy Survey (DES) covering $\sim$5000 deg$^2$. We perform a large-scale structure analysis using three two-point correlation functions (3$\times$2pt): (i) cosmic shear from 140 million source galaxy shapes, (ii) galaxy clustering of 9 million lens galaxy positions, and (iii) galaxy-galaxy lensing from their cross-correlation. We model the data in flat $Λ$CDM and $w$CDM cosmologies. The combined analysis yields $S_8\equiv σ_8 (Ω_{\rm m}/0.3)^{0.5} = 0.789^{+0.012}_{-0.012}$ and matter density $Ω_{\rm m} = 0.333^{+0.023}_{-0.028}$ in $Λ$CDM (68% CL), where $σ_8$ is the clustering amplitude. These constraints show a (full-space) parameter difference of 1.8$σ$ from a combination of cosmic microwave background (CMB) primary anisotropy datasets from Planck 2018, ACT-DR6, and SPT-3G DR1. Projected only into $S_8$ the difference is $2.6σ$. In $w$CDM the Y6 3$\times$2pt results yield $S_8 = 0.782^{+0.021}_{-0.020}$, $Ω_{\rm m} = 0.325^{+0.032}_{-0.035}$, and dark energy equation-of-state parameter $w = -1.12^{+0.26}_{-0.20}$. For the first time, we combine all DES dark-energy probes: 3$\times$2pt, SNe Ia, BAO and Clusters. In $Λ$CDM this combination yields a $2.8σ$ parameter difference from the CMB. When combining DES 3$\times$2pt with other low-redshift datasets (DESI DR2 BAO, DES SNe Ia, SPT clusters), we find a 2.3$σ$ parameter difference with CMB. A joint fit of Y6 3$\times$2pt, CMB, and those low-redshift datasets produces the tightest $Λ$CDM constraints to date: $S_8 = 0.806^{+0.006}_{-0.007}$, $Ω_{\rm m} = 0.302^{+0.003}_{-0.003}$, $h = 0.683^{+0.003}_{-0.002}$, and $\sum m_ν< 0.14$ eV (95% CL). In $w$CDM, this combination yields $w = -0.981^{+0.021}_{-0.022}$.

Abbott, T. M.C. [Cerro-Tololo InterAmerican Obs.]↗
DOE OSTI2021

Variational design principles for nonequilibrium colloidal assembly

Using large deviation theory and principles of stochastic optimal control, we show that rare molecular dynamics trajectories conditioned on assembling a specific target structure encode a set of interactions and external forces that lead to enhanced stability of that structure. Such a relationship can be formulated into a variational principle, for which we have developed an associated optimization algorithm and have used it to determine optimal forces for targeted self-assembly within nonequilibrium steady-states. We illustrate this perspective on inverse design in a model of colloidal cluster assembly within linear shear flow. We find that colloidal clusters can be assembled with high yield using specific short-range interactions of tunable complexity. Shear decreases the yields of rigid clusters, while small values of shear increase the yields of nonrigid clusters. The enhancement or suppression of the yield due to shear is rationalized with a generalized linear response theory. Furthermore, by studying 21 unique clusters made of six, seven, or eight particles, we uncover basic design principles for targeted assembly out of equilibrium.

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

Robust clustering of the local Milky Way stellar kinematic substructures with Gaia eDR3

Understanding local stellar kinematic substructures in the solar neighbourhood helps build a complete picture of the formation of the Milky Way, as well as an empirical phase space distribution of dark matter that would inform detection experiments. We apply the clustering algorithm HDBSCAN on the Gaia early third data release to identify a list of stable clusters in velocity space and action-angle space by taking into account the measurement uncertainties and studying the stability of the clustering results. We find 1405 (497) stars in 23 (6) robust clusters in velocity space (action-angle space) that are consistently not associated with noise. We discuss the kinematic properties of these structures and study whether many of the small clusters belong to a similar larger cluster based on their chemical abundances. They are attributed to the known structures: the Gaia Sausage-Enceladus, the Helmi Stream, and globular cluster NGC 3201 are found in both spaces, while NGC 104 and the thick disc (Sequoia) are identified in velocity space (action-angle space). Although we do not identify any new structures, we find that the HDBSCAN member selection of already known structures is unstable to input kinematics of the stars when resampled within their uncertainties. We therefore present the stable subset of local kinematic structures, which are consistently identified by the clustering algorithm, and emphasize the need to take into account error propagation during both the manual and automated identification of stellar structures, both for existing ones as well as future discoveries.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The SRG/eROSITA All-Sky Survey: Exploring halo assembly bias with X-ray-selected superclusters

Numerical simulations indicate that the clustering of dark matter halos is not only dependent on the halo masses but has a secondary dependence on other properties, such as the assembly history of the halo. This phenomenon, known as the halo assembly bias (HAB), has been found mostly on galaxy scales; observational evidence on larger scales is scarce. In this work, we propose a novel method for exploring HAB on cluster scales using large samples of superclusters. Leveraging the largest-ever X-ray galaxy cluster and supercluster samples obtained from the first SRG/eROSITA all-sky survey, we constructed two subsamples of galaxy clusters that consist of supercluster members and isolated clusters, respectively. After correcting for the selection effects on redshift, mass, and survey depth, we computed the excess in the concentration of the intracluster gas of isolated clusters with respect to supercluster members, defined as δc gas ≡ c gas , ISO /c gas , SC –1 , to investigate the environmental effect on the concentration of clusters, a sign of HAB on cluster scales. We find that the average gas mass concentration of isolated clusters is a few percent higher than that of supercluster members, with a maximum significance of 2.8σ. The result for δc gas varies with the overdensity ratio, f, in supercluster identification, cluster mass proxies, and mass and redshift ranges but remains positive in almost all the measurements. We measure slightly larger δc gas when adopting a higher f for supercluster identification. The δc gas is also higher for low-mass and low-redshift clusters. We performed weak lensing analyses to compare the total mass concentration of the two classes and find a similar trend in total mass concentration as obtained from the gas mass concentration. Our results are consistent with the prediction of HAB on cluster scales, where halos located in denser environments are less concentrated; this trend is stronger for halos with lower masses and at lower redshifts. These phenomena can be explained by the fact that clusters in denser environments, such as superclusters, have experienced more mergers than isolated clusters in their assembling history. This work paves the way to explore HAB with X-ray superclusters and demonstrates that large samples of superclusters with X-ray and weak-lensing data can advance our understanding of the evolution of the large-scale structure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Possibilities and Limitations of Kinematically Identifying Stars from Accreted Ultra-faint Dwarf Galaxies

Abstract The Milky Way has accreted many ultra-faint dwarf galaxies (UFDs), and stars from these galaxies can be found throughout our Galaxy today. Studying these stars provides insight into galaxy formation and early chemical enrichment, but identifying them is difficult. Clustering stellar dynamics in 4D phase space ( E , L z , J r , J z ) is one method of identifying accreted structure that is currently being utilized in the search for accreted UFDs. We produce 32 simulated stellar halos using particle tagging with the Caterpillar simulation suite and thoroughly test the abilities of different clustering algorithms to recover tidally disrupted UFD remnants. We perform over 10,000 clustering runs, testing seven clustering algorithms, roughly twenty hyperparameter choices per algorithm, and six different types of data sets each with up to 32 simulated samples. Of the seven algorithms, HDBSCAN most consistently balances UFD recovery rates and cluster realness rates. We find that, even in highly idealized cases, the vast majority of clusters found by clustering algorithms do not correspond to real accreted UFD remnants and we can generally only recover 6% of UFDs remnants at best. These results focus exclusively on groups of stars from UFDs, which have weak dynamic signatures compared to the background of other stars. The recoverable UFD remnants are those that accreted recently, z accretion ≲ 0.5. Based on these results, we make recommendations to help guide the search for dynamically linked clusters of UFD stars in observational data. We find that real clusters generally have higher median energy and J r , providing a way to help identify real versus fake clusters. We also recommend incorporating chemical tagging as a way to improve clustering results.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Dark Energy Survey Year 6 results: Intra-cluster light from redshift 0.2 to 0.5

Using the full 6 years of imaging data from the Dark Energy Survey, we study the surface brightness profiles of galaxy cluster central galaxies and intra-cluster light. We apply a ‘stacking’ method to over 4000 galaxy clusters identified by the redMaPPer cluster finder in the redshift range of 0.2−0.5. This yields high-signal-to-noise circularly averaged profile measurements of the central galaxy and intra-cluster light out to 1 Mpc from the cluster centre. Using redMaPPer richness as a cluster mass indicator, we find that the brightness of the intra-cluster light has a strong mass dependence throughout the 0.2−0.5 redshift range, and this dependence grows stronger at a larger radius. In terms of redshift evolution, we find some evidence that the central galaxy, as well as the diffuse light within the transition region between the cluster central galaxy and intra-cluster light within 80 kpc from the centre, may be growing over time. At larger radii, more than 80 kpc away from the cluster centre, we do not detect evidence of additional redshift evolution beyond the cluster mass dependence, consistent with the findings from the IllustrisTNG hydrodynamic simulation. We speculate that the major driver of intra-cluster light growth, especially at large radii, is associated with cluster mass growth. Finally, we find that the colour of the cluster central galaxy and intra-cluster light displays a radial gradient that becomes bluer at a larger radius, which is consistent with a stellar stripping and disruption origin of intra-cluster light as suggested by simulation studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

ADDGALS: Simulated Sky Catalogs for Wide Field Galaxy Surveys

Abstract We present a method for creating simulated galaxy catalogs with realistic galaxy luminosities, broadband colors, and projected clustering over large cosmic volumes. The technique, denoted Addgals (Adding Density Dependent GAlaxies to Lightcone Simulations), uses an empirical approach to place galaxies within lightcone outputs of cosmological simulations. It can be applied to significantly lower-resolution simulations than those required for commonly used methods such as halo occupation distributions, subhalo abundance matching, and semi-analytic models, while still accurately reproducing projected galaxy clustering statistics down to scales of r ∼ 100 h −1 kpc . We show that Addgals catalogs reproduce several statistical properties of the galaxy distribution as measured by the Sloan Digital Sky Survey (SDSS) main galaxy sample, including galaxy number densities, observed magnitude and color distributions, as well as luminosity- and color-dependent clustering. We also compare to cluster–galaxy cross correlations, where we find significant discrepancies with measurements from SDSS that are likely linked to artificial subhalo disruption in the simulations. Applications of this model to simulations of deep wide-area photometric surveys, including modeling weak-lensing statistics, photometric redshifts, and galaxy cluster finding, are presented in DeRose et al., and an application to a full cosmology analysis of Dark Energy Survey (DES) Year 3 like data is presented in DeRose et al. We plan to publicly release a 10,313 square degree catalog constructed using Addgals with magnitudes appropriate for several existing and planned surveys, including SDSS, DES, VISTA, Wide-field Infrared Survey Explorer, and Rubin Observatory’s Legacy Survey of Space and Time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Density functional theory of water with the machine-learned DM21 functional

The delicate interplay between functional-driven and density-driven errors in density functional theory (DFT) has hindered traditional density functional approximations (DFAs) from providing an accurate description of water for over 30 years. Recently, the deep-learned DeepMind 21 (DM21) functional has been shown to overcome the limitations of traditional DFAs as it is free of delocalization error. To determine if DM21 can enable a molecular-level description of the physical properties of aqueous systems within Kohn–Sham DFT, we assess the accuracy of the DM21 functional for neutral, protonated, and deprotonated water clusters. Furthermore, we find that the ability of DM21 to accurately predict the energetics of aqueous clusters varies significantly with cluster size. Additionally, we introduce the many-body MB-DM21 potential derived from DM21 data within the many-body expansion of the energy and use it in simulations of liquid water as a function of temperature at ambient pressure. We find that size-dependent functional-driven errors identified in the analysis of the energetics of small clusters calculated with the DM21 functional result in the MB-DM21 potential systematically overestimating the hydrogen-bond strength and, consequently, predicting a more ice-like local structure of water at room temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗