SEARCH · Engineering Papers

Results for “clusters: general”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Benchmarking an Embedded Adaptive Sampling Configuration Interaction Method for Surface Reactions: H 2 Desorption from and CH 4 Dissociation on Cu(111)

Embedded (emb-) correlated wavefunction (CW) theory enables accurate assessments of both ground- and excited-state reaction mechanisms involved in heterogeneous catalysis. Embedded multireference second-order perturbation theory (emb-MRPT2) based on reference wavefunctions generated via embedded complete active space self-consistent field (emb-CASSCF) theory is currently state-of-the-art. However, the factorial scaling of CASSCF limits the size of active space and the complexity of systems that can be studied. In this work, we assess the efficacy of an alternative CW method, adaptive sampling configuration interaction (ASCI)–which enables large active spaces to be used–for studying surface reactions. We couple ASCI with density functional embedding theory (DFET) and benchmark its performance for two reactions: H 2 desorption from and CH 4 dissociation on the Cu(111) surface. Unlike embedded complete active space second-order perturbation theory (emb-CASPT2) that accurately reproduces a measured H 2 desorption barrier, embedded ASCI, using a very large active space (though one that still comprises a small portion of the full set of orbitals) fails to do so. Adding an extra correlation term from embedded Møller–Plesset second-order perturbation theory (emb-MP2) improves the desorption barrier and endothermicity predictions. Thus, the inaccuracy of embedded ASCI comes from the missing dynamic correlation from the many other electrons and orbitals not included in the active space. For CH 4 dissociation, again embedded ASCI overestimates the dissociation barrier compared to emb-CASPT2 predictions. Adding dynamic correlation from emb-MP2 helps correct the barrier. However, this composite approach suffers from double counting of correlation within embedded ASCI followed by emb-MP2 calculations. We therefore conclude that the state-of-the-art emb-MRPT2 based on reference wavefunctions generated via emb-CASSCF remains the method of choice for studying surface reactions. emb-ASCI is useful when large active spaces beyond the limit of emb-CASSCF are essential, such as to study complex surface reactions with significant multiconfigurational character (static correlation) but weak dynamic correlation.

08 HYDROGEN↗

DOE OSTI2019

Stellar mass–halo mass relation for the brightest central galaxies of X-ray clusters since z ~ 0.65

We present the brightest cluster galaxies (BCGs) catalog for SPectroscoic IDentification of eROSITA Sources (SPIDERS) DR14 cluster program value-added catalog. We list the 416 BCGs identified as part of this process, along with their stellar mass, star formation rates (SFRs), and morphological properties. We identified the BCGs based on the available spectroscopic data from SPIDERS and photometric data from SDSS. We computed stellar masses and SFRs of the BCGs on the basis of SDSS, WISE, and GALEX photometry using spectral energy distribution fitting. Morphological properties for all BCGs were derived by Sersic profile fitting using the software package SIGMA in different optical bands ($g,r,i$). We combined this catalog with the BCGs of galaxy groups and clusters extracted from the deeper AEGIS, CDFS, COSMOS, XMM-CFHTLS, and XMM-XXL surveys to study the stellar mass–halo mass relation using the largest sample of X-ray groups and clusters known to date. This result suggests that the mass growth of the central galaxy is controlled by the hierarchical mass growth of the host halo. We find a strong correlation between the stellar mass of BCGs and the mass of their host halos. This relation shows no evolution since $z$ ~ 0.65. We measure a mean scatter of 0.21 and 0.25 for the stellar mass of BCGs in a given halo mass at low (0.1 < $z$ < 0.3) and high (0.3 < $z$ < 0.65) redshifts, respectively. We further demonstrate that the BCG mass is covariant with the richness of the host halos in the very X-ray luminous systems. We also find evidence that part of the scatter between X-ray luminosity and richness can be reduced by considering stellar mass as an additional variable.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Accurate prediction of core-level spectra of radicals at density functional theory cost via square gradient minimization and recoupling of mixed configurations

State-specific orbital optimized approaches are more accurate at predicting core-level spectra than traditional linear-response protocols, but their utility had been restricted due to the risk of "variational collapse" down to the ground state. We employ the recently developed square gradient minimization algorithm to reliably avoid variational collapse and study the effectiveness of orbital optimized density functional theory (DFT) at predicting second period element 1s core-level spectra of open-shell systems. Several density functionals (including SCAN, B3LYP, and ωB97X-D3) are found to predict excitation energies from the core to singly occupied levels with high accuracy (≤0.3 eV RMS error) against available experimental data. Higher excited states are, however, more challenging by virtue of being intrinsically multiconfigurational. We thus present a configuration interaction inspired route to self-consistently recouple single determinant mixed configurations obtained from DFT, in order to obtain approximate doublet states. This recoupling scheme is used to predict the C K-edge spectra of the allyl radical, the O K-edge spectra of CO + , and the N K-edge of NO 2 with high accuracy relative to experiment, indicating substantial promise in using this approach for the computation of core-level spectra for doublet species [vs more traditional time dependent DFT, equation of motion coupled cluster singles and doubles (EOM-CCSD), or using unrecoupled mixed configurations]. Finally, we also present general guidelines for computing core-excited states from orbital optimized DFT.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Completion Effects on Diagnosing Multistage Facture Treatments with Distributed Temperature Sensing

Summary Distributed temperature sensing (DTS) is a valuable tool to diagnose multistage hydraulic fracture treatments. When a stage interval is shut in, the clusters that take more fluid during pumping warm up more slowly. Therefore, the fluid volume injected into each cluster can be quantitatively interpreted by numerical inversion of the warm-back temperature behavior. This general concept assumes that the different warm-back behavior is controlled by only the injected fluid volume; however, recent observations of DTS data indicate that completion configurations significantly influence the warm-back behavior. This paper investigates the completion effects on the DTS interpretation. In ideal conditions, when a stage is fractured, the upstream stage intervals should show an almost uniform temperature that is close to the injected fluid temperature. This is due to the high fluid velocity of injected fluid in the wellbore, and the upstream intervals have not been perforated (noncommunicating intervals), so the only heat transfer is heat conduction between the wellbore fluid and the surrounding reservoir. But the field DTS data show considerably irregular variations in temperature along the upstream stage intervals. These variations are caused by the completion effects. The nonuniform temperature profile is caused by different heat transfer behavior induced by completion hardware along the production casing string, such as joints, clamps, and blast protectors, and by the sensing cable location in the cement, as well as the cement quality. Because the varying heat transfer behavior impacts the warm-back behavior as well as the temperature profile, the completion effects need to be considered in DTS interpretation. A method of DTS interpretation considering the completion effects to diagnose multistage fracture treatments was developed. Because the heat transfer between a wellbore and a reservoir depends on the overall heat transfer coefficient describing heat conduction through the completion in a forward model, this parameter needs to be tuned all along the wellbore. To calibrate the completion effect, the temperature inversion is conducted using the temperature measured at a stage interval that is upstream of a stage interval currently being treated. Because the interpreted stage interval is not perforated at that time, the thermal behavior at the noncommunicating interval is governed by only the heat conduction through the completion environment. Once the effective values of the overall heat transfer coefficient are estimated along the interpreted stage interval, they can be assumed to be constant physical parameters. Then, the fluid volume distribution is interpreted by using the effective overall heat transfer coefficient profile along each interval. This study provides a field application of the developed interpretation method. The new interpretation method provides more accurate diagnosis of fracture treatments by DTS interpretation.

Engineering↗

DOE OSTI2023

First Constraints on Growth Rate from Redshift-space Ellipticity Correlations of SDSS Galaxies at 0.16 < z < 0.70

We report the first constraints on the growth rate of the universe, f(z)σ 8 (z), with intrinsic alignments (IAs) of galaxies. We measure the galaxy density-intrinsic ellipticity cross-correlation and intrinsic ellipticity autocorrelation functions over 0.16 < z < 0.7 from luminous red galaxies (LRGs) and LOWZ and CMASS galaxy samples in the Sloan Digital Sky Survey (SDSS) and SDSS-III BOSS survey. We detect clear anisotropic signals of IA due to redshift-space distortions. By combining measured IA statistics with the conventional galaxy clustering statistics, we obtain tighter constraints on the growth rate. The improvement is particularly prominent for the LRG, which is the brightest galaxy sample and known to be strongly aligned with underlying dark matter distribution; using the measurements on scales above 10 h -1 Mpc, we obtain $f{\sigma }_{8}={0.5196}_{-0.0354}^{+0.0352}$ (68% confidence level) from the clustering-only analysis and $f{\sigma }_{8}={0.5322}_{-0.0291}^{+0.0293}$ with clustering and IA, meaning 19% improvement. The constraint is in good agreement with the prediction of general relativity, f σ 8 = 0.4937 at z = 0.34. For LOWZ and CMASS samples, the improvement of constraints on f σ 8 is found to be 10% and 3.5%, respectively. Our results indicate that the contribution from IA statistics for cosmological constraints can be further enhanced by carefully selecting galaxies for a shape sample.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

Progress on Grizzly Development for Reactor Pressure Vessels and Reinforced Concrete Structures

The Grizzly code has been under development for the U.S. Department of Energy’s Light Water Reactor Sustainability program to provide predictive tools for the evolution of material degradation in critical light water reactor structural components due to long-term exposure to the environmental conditions of normal reactor operation, and for the effects of this degradation on the safety of these components. This development has primarily targeted reactor pressure vessels and reinforced concrete structures. Work performed during Fiscal Year 2019 added important features for both of these types of structures. For probabilistic fracture mechanics analysis of reactor pressure vessels, capabilities were added to account for the effects of warm prestressing and residual stresses, as well as to model crack initiation, growth, and arrest. For modeling microstructure evolution in reactor pressure vessels, a cluster dynamics model for evolution of Mn-Ni-Si phases has been implemented. For modeling degraded reinforced concrete structures, capabilities to consider the combined effects of damage and creep in the concrete constitutive model and model nonlinear behavior of reinforcing bars has been added. These newly-developed features fill important gaps in Grizzly’s feature set for engineering-scale analysis of degraded reactor pressure vessels and reinforced concrete structures. Grizzly development in both of these areas has reached a point where it is ready for extensive testing on relevant engineering problems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Can assembly bias explain the lensing amplitude of the BOSS CMASS sample in a Planck cosmology?

ABSTRACT In this paper, we investigate whether galaxy assembly bias can reconcile the 20–40 ${{\ \rm per\ cent}}$ disagreement between the observed galaxy projected clustering signal and the galaxy–galaxy lensing signal in the Baryon Oscillation Spectroscopic Survey CMASS galaxy sample. We use the suite of abacuscosmos lambda cold dark matter simulations at Planck best-fitting cosmology and two flexible implementations of extended halo occupation distribution (HOD) models that incorporate galaxy assembly bias to build forward models and produce joint fits of the observed galaxy clustering signal and the galaxy–galaxy lensing signal. We find that our models using the standard HODs without any assembly bias generalizations continue to show a 20–40 ${{\ \rm per\ cent}}$ overprediction of the observed galaxy–galaxy lensing signal. We find that our implementations of galaxy assembly bias do not reconcile the two measurements at Planck best-fitting cosmology. In fact, despite incorporating galaxy assembly bias, the satellite distribution parameter, and the satellite velocity bias parameter into our extended HOD model, our fits still strongly suggest a $\sim \! 34{{\ \rm per\ cent}}$ discrepancy between the observed projected clustering and galaxy–galaxy lensing measurements. It remains to be seen whether a combination of other galaxy assembly bias models, alternative cosmological parameters, or baryonic effects can explain the amplitude difference between the two signals.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

DESI 2024 II: sample definitions, characteristics, and two-point clustering statistics

We present the samples of galaxies and quasars used for DESI 2024 cosmological analyses, drawn from the DESI Data Release 1 (DR1). We describe the construction of largescale structure (LSS) catalogs from these samples, which include matched sets of synthetic reference ‘randoms’ and weights that account for variations in the observed density of the samples due to experimental design and varying instrument performance. We detail how we correct for variations in observational completeness, the input ‘target’ densities due to imaging systematics, and the ability to confidently measure redshifts from DESI spectra. We then summarize how remaining uncertainties in the corrections can be translated to systematic uncertainties for particular analyses. We describe the weights added to maximize the signalto-noise of DESI DR1 2-point clustering measurements. We detail measurement pipelines applied to the LSS catalogs that obtain 2-point clustering measurements in configuration and Fourier space. The resulting 2-point measurements depend on window functions and normalization constraints particular to each sample, and we present the corrections required to match models to the data. We compare the configuration- and Fourier-space 2-point clustering of the data samples to that recovered from simulations of DESI DR1 and find they are, generally, in statistical agreement to within 2% in the inferred real-space over-density field. The LSS catalogs, 2-point measurements, and their covariance matrices will be released publicly with DESI DR1.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Thermofield Theory for Finite-Temperature Electronic Structure

Wave function methods have offered a robust, systematically improvable means to study ground-state properties in quantum many-body systems. Theories like coupled cluster and their derivatives provide highly accurate approximations to the energy landscape at a reasonable computational cost. Analogues of such methods to study thermal properties, though highly desirable, have been lacking because evaluating thermal properties involve a trace over the entire Hilbert space, which is a formidable task. Besides, excited-state theories are generally not as well studied as ground-state ones. In this mini-review, we present an overview of a finite-temperature wave function formalism based on thermofield dynamics to overcome these difficulties. Thermofield dynamics allows us to map the equilibrium thermal density matrix to a pure state, i.e., a single wave function, albeit in an expanded Hilbert space. Ensemble averages become expectation values over this so-called thermal state. Around this thermal state, we have developed a procedure to generalize ground-state wave function theories to finite temperatures. As explicit examples, we highlight formulations of mean-field, configuration interaction, and coupled cluster theories for thermal properties of Fermions in the grand-canonical ensemble. To assess the quality of these approximations, we also show benchmark studies for the one-dimensional Hubbard model, while comparing against exact results. We will see that the thermal methods perform similarly to their ground-state counterparts, while merely adding a prefactor to the asymptotic computational cost. Furthermore, they also inherit all the properties, good or bad, from the ground-state methods, signifying the robustness of our formalism and the scope for future development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

The history of metal enrichment traced by X-ray observations of high-redshift galaxy clusters

ABSTRACT We present the analysis of deep X-ray observations of 10 massive galaxy clusters at redshifts 1.05 < z < 1.71, with the primary goal of measuring the metallicity of the intracluster medium (ICM) at intermediate radii, to better constrain models of the metal enrichment of the intergalactic medium. The targets were selected from X-ray and Sunyaev–Zel’dovich effect surveys, and observed with both the XMM–Newton and Chandra satellites. For each cluster, a precise gas mass profile was extracted, from which the value of r500 could be estimated. This allows us to define consistent radial ranges over which the metallicity measurements can be compared. In general, the data are of sufficient quality to extract meaningful metallicity measurements in two radial bins, r < 0.3r500 and 0.3 < r/r500 < 1.0. For the outer bin, the combined measurement for all 10 clusters, Z/Z⊙ = 0.21 ± 0.09, represents a substantial improvement in precision over previous results. This measurement is consistent with, but slightly lower than, the average metallicity of 0.315 solar measured at intermediate-to-large radii in low-redshift clusters. Combining our new high-redshift data with the previous low-redshift results allows us to place the tightest constraints to date on models of the evolution of cluster metallicity at intermediate radii. Adopting a power-law model of the form Z ∝ (1 + z)γ, we measure a slope $\gamma = -0.5^{+0.4}_{-0.3}$, consistent with the majority of the enrichment of the ICM having occurred at very early times and before massive clusters formed, but leaving open the possibility that some additional enrichment in these regions may have occurred since a redshift of 2.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Effective interactions between nuclear clusters

In this work, the effective interactions between two nuclear clusters, d + d, t + t, and α + α, are investigated within a cluster model using local nucleon-nucleon (NN) forces. It is shown that the interaction in the spin-aligned d + d system is repulsive for all intercluster distances, whereas the α + α and spin-aligned t + t systems are attractive at intermediate distances. The Pauli blocking between identical-nucleon pairs is responsible for the cluster-cluster repulsion and becomes dominant in the shallow binding limit. We demonstrate that two d clusters could be bound if the NN force has nonzero range and is strong enough to form a deeply bound d cluster, or if the NN force has both even-parity and odd-parity attraction. Effective dimer-dimer interactions for general quantum systems of two-component fermions are also discussed in heavy-light mass limit, where one component is much heavier than the other, and their relation to intercluster interactions in nuclear systems are discussed. Our findings provide a conceptual foundation for conclusions obtained numerically in the literature, that increasing the range or strength of the local part of the attractive nucleon-nucleon interaction results in a more attractive cluster-cluster interaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

A theoretical investigation of the hydrolysis of uranium hexafluoride: the initiation mechanism and vibrational spectroscopy

Depleted uranium hexafluoride (UF 6 ), a stockpiled byproduct of the nuclear fuel cycle, reacts readily with atmospheric humidity, but the mechanism is poorly understood. Here we compare several potential initiation steps at a consistent level of theory, generating underlying structures and vibrational modes using hybrid density functional theory (DFT) and computing relative energies of stationary points with double-hybrid (DH) DFT. A benchmark comparison is performed to assess the quality of DH-DFT data using reference energy differences obtained using a complete-basis-limit coupled-cluster (CC) composite method. The associated large-basis CC computations were enabled by a new general-purpose pseudopotential capability implemented as part of this work. Dispersion-corrected parameter-free DH-DFT methods, namely PBE0-DH-D3(BJ) and PBE-QIDH-D3(BJ), provided mean unsigned errors within chemical accuracy (1 kcal mol -1 ) for a set of barrier heights corresponding to the most energetically favorable initiation steps. The hydrolysis mechanism is found to proceed via intermolecular hydrogen transfer within van der Waals complexes involving UF 6 , UF 5 OH, and UOF 4 , in agreement with previous studies, followed by the formation of a previously unappreciated dihydroxide intermediate, UF 4 (OH) 2 . The dihydroxide is predicted to form under both kinetic and thermodynamic control, and, unlike the alternate pathway leading to the UO 2 F 2 monomer, its reaction energy is exothermic, in agreement with observation. Finally, harmonic and anharmonic vibrational simulations are performed to reinterpret literature infrared spectroscopy in light of this newly identified species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Residual mean field model of valence quarks in the nucleon

Abstract We develop a non-perturbative model for valence parton distribution functions (PDFs) based on the mean field interactions of valence quarks in the nucleonic interior. The main motivation for the model is to obtain a mean field description of the valence quarks as a baseline to study the short range quark–quark interactions that generate the high x tail of PDFs. The model is based on the separation of the valence three-quark cluster and residual system in the nucleon. Then the nucleon structure function is calculated within the effective light-front diagrammatic approach introducing nonperturbative light-front valence quark and residual wave functions. Within the model a new relation is obtained between the position, $$x_p$$ x p , of the peak of $$xq_V(x)$$ x q V ( x ) distribution of the valence quark and the effective mass of the residual system, $$m_R$$ m R , in the form: $$x_{p} \approx {1\over 4} (1-{m_R\over m_N})$$ x p ≈ 1 4 ( 1 - m R m N ) at starting $$Q^2$$ Q 2 . This relation explains the difference in the peak positions for d- and u-quarks through the expected difference of residual masses for valence d- and u-quark distributions. The parameters of the model are fixed by fitting the calculated valence quark distributions to the phenomenological PDFs. This allowed us to estimate the overall mean field contribution in baryonic and momentum sum rules for valence d- and u-quarks. Finally, the evaluated parameters of the non-perturbative wave functions of valence 3q-cluster and residual system can be used in calculation of other quantities such as nucleon form factors, generalized partonic and transverse momentum distributions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Stability and mobility of tungsten clusters on tungsten (110) surface: Ab initio and atomistic simulations

Quantifying the surface transport properties of tungsten (W) is of prime importance to understand the formation of nano-fuzz in fusion plasma-facing conditions. The stability and mobility of W adatom clusters (W n , n = 2-9) on the W(110) surface has been investigated by computer simulations, including ab initio calculations using density functional theory (DFT) and molecular statics (MS) simulations with multiple W interatomic potentials. The DFT results demonstrate that the sequential binding energy generally increases with number of W adatoms, except for the 5 th and 7 th W adatoms. The most common elemental migration steps of W n (n>2) clusters are observed to consist of monomer and dimer hops, while larger W n clusters can also diffuse by dissociation and recombination of smaller clusters. The threshold migration energy of W 9 is the highest, then followed by W 8 , W 4 , and W 6 , while W 3 , W 5 , and W 7 have similar migration energies. Compared to DFT, each interatomic potential evaluated overestimates the binding energies of W n clusters. Additionally, the embedded-atom potential developed by Juslin and Wirth adequately predicts the threshold migration energy of W n (n>2) clusters, although it predicts different underlying migration mechanisms. The results show that interaction mechanism between W adatoms controls the stability and mobility of W n clusters on the W(110) surface.

36 MATERIALS SCIENCE↗

DOE OSTI2025

The SRG/eROSITA All-Sky Survey: Constraints on ultralight axion dark matter through galaxy cluster number counts

Ultralight axions are hypothetical scalar particles that influence the evolution of large-scale structures of the Universe. Depending on their mass, they can potentially be part of the dark matter component of the Universe as candidates commonly referred to as fuzzy dark matter. While strong constraints have been established for pure fuzzy dark matter models, the more general scenario where ultralight axions constitute only a fraction of the dark matter has been limited to only a few observational probes. In this work, we use the galaxy cluster number counts obtained from the first All-Sky Survey (eRASS1) of the SRG/eROSITA mission together with gravitational weak lensing data from the Dark Energy Survey, the Kilo-Degree Survey, and the Hyper Suprime-Cam to constrain the fraction of ultralight axions in the mass range 10 −32 eV to 10 −24 eV. We put upper bounds on the ultralight axion relic density Ω a in independent logarithmic axion mass bins by performing a full cosmological parameter inference. We find an exclusion region in the intermediate ultralight axion mass regime with the tightest bounds reported so far in the mass bins around m a = 10 −27 eV with Ω a < 0.0035 and m a = 10 −26 eV with Ω a < 0.0079; both are at a 95% confidence level. When combined with cosmic microwave background probes, these bounds are tightened to Ω a < 0.0030 in the m a = 10 −27 eV mass bin and Ω a < 0.0058 in the m a = 10 −26 eV mass bin, with both at a 95% confidence level. This is the first time that constraints on ultralight axions have been obtained using the growth of structure measured by galaxy cluster number counts. These results pave the way for large surveys, which can be utilized to obtain tight constraints on the mass and relic density of ultralight axions with better theoretical modeling of the abundance of halos.

cosmological parameters↗

DOE OSTI2020

The Generality of the GUGA MRCI Approach in COLUMBUS for Treating Complex Quantum Chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Water Network Shape-Dependence of Local Interactions with the Microhydrated –NO 2 – and –CO 2 – Anionic Head Groups by Cold Ion Vibrational Spectroscopy

We report the structural evolutions of water networks and solvatochromic response of the CH 3 NO 2 – radical anion in the OH and CH stretching regions by analysis of the vibrational spectra displayed by cryogenically cooled CH 3 NO 2 – ·(H 2 O) n=1–6 clusters. The OH stretching bands evolve with a surprisingly large discontinuity at n = 6, which features the emergence of an intense, strongly red-shifted band along with a weaker feature that appears in the region assigned to a free OH fundamental. Very similar behavior is displayed by the perdeuterated carboxylate clusters, RCO 2 – ·(H 2 O) n=5–7 (R = CD 3 CD 2 ), indicating that this behavior is a general feature in the microhydration of the triatomic anionic domain and not associated with CH oscillators. Electronic structure calculations trace this behavior to the formation of a “book” isomer of the water hexamer that adopts a configuration in which one of the water molecules resides in an acceptor–acceptor–donor (AAD) (A = acceptor, D = donor) H-bonding site. Excitation of the bound OH in the AAD site explores the local network topology best suited to stabilize an incipient –XO 2 H–OH–(H 2 O) 2 intracluster proton-transfer reaction. These systems thus provide particularly clear examples where the network shape controls the potential energy landscape that governs water network-mediated, intracluster proton transfer. The CH stretching bands of the CH 3 NO 2 – ·(H 2 O) n=1–6 clusters also exhibit strong solvatochromic shifts, but in this case, they smoothly blue-shift with increasing hydration with no discontinuity at n = 6. Furthermore, this behavior is analyzed in the context of the solute-ion polarizability response and partial charge transfer to the water networks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Uncertainty-aware molecular dynamics from Bayesian active learning for phase transformations and thermal transport in SiC

Abstract Machine learning interatomic force fields are promising for combining high computational efficiency and accuracy in modeling quantum interactions and simulating atomistic dynamics. Active learning methods have been recently developed to train force fields efficiently and automatically. Among them, Bayesian active learning utilizes principled uncertainty quantification to make data acquisition decisions. In this work, we present a general Bayesian active learning workflow, where the force field is constructed from a sparse Gaussian process regression model based on atomic cluster expansion descriptors. To circumvent the high computational cost of the sparse Gaussian process uncertainty calculation, we formulate a high-performance approximate mapping of the uncertainty and demonstrate a speedup of several orders of magnitude. We demonstrate the autonomous active learning workflow by training a Bayesian force field model for silicon carbide (SiC) polymorphs in only a few days of computer time and show that pressure-induced phase transformations are accurately captured. The resulting model exhibits close agreement with both ab initio calculations and experimental measurements, and outperforms existing empirical models on vibrational and thermal properties. The active learning workflow readily generalizes to a wide range of material systems and accelerates their computational understanding.

36 MATERIALS SCIENCE↗