Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density correlation”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The Bonding Nature and Adhesion of Polyacrylic Acid Coating on Li-Metal for Li Dendrite Prevention

The success of polyacrylic acid (PAA) to suppress Li dendrite growth suggests that the mechanical properties of polymer-based coatings, including the modulus, toughness, and interfacial adhesion are important design criteria. However, the measurement of the adhesion of thin PAA, as well as other polymer coatings to the reactive Li-metal anode surface is limited experimentally and challenging computationally. In this paper, a strategy was proposed to estimate the adhesion and delamination of the PAA(polymer)/Li interface, based on the bonding nature at the simpler PAA (oligomer)/Li interfaces using density functional theory calculations. It has been shown that the carboxylic acid groups in PAA reacted strongly with metallic Li, which significantly enhances the interfacial adhesion through the Li–O bonds formation, Li ionization and its incorporation into PAA, and –H or –OH termination of Li after decomposition of the COOH functional group. During delamination, it was found that the most likely PAA delamination route involved breaking partial Li–O bonds and lifting some ionized Li atoms from the Li-metal, especially for the Li atoms that showed a charge closer to +1 or are bonded with two O atoms from PAA. Based on the average bonding energies from PAA(oligomer)/Li interface delamination calculations, the work of separation, W sep , of the PAA(polymer)/Li interface was estimated to be ~1.0 (J/m 2 ). The high W sep of PAA (polymer)/Li was comparable with the Li 2 O/Li interface and higher than Li 2 CO 3 /Li and LiF/Li interfaces. This order correlated well with the areal density of Li–O bonds, which can serve as a descriptor for the interfacial adhesion. Furthermore, this computational approach can be applied to other interfaces with polymer-based coatings.

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

California-Kepler Survey. IX. Revisiting the Minimum-mass Extrasolar Nebula with Precise Stellar Parameters

We investigate a possible correlation between the solid surface density Σ of the minimum-mass extrasolar nebula (MMEN) and the host star mass M {sub ⋆} and metallicity [Fe/H]. Leveraging on the precise host star properties from the California-Kepler Survey (CKS), we found that Σ=50{sub −20}{sup +33} g cm{sup −2} (a/1 au){sup −1.75±0.07} (M {sub ⋆}/M {sub ⊙}){sup 1.04±0.22} 10{sup 0.22±0.05[Fe/H]} for Kepler-like systems (1–4R {sub ⊕}; a < 1 au). The strong M {sub ⋆} dependence is reminiscent of previous dust continuum results that the solid disk mass scales with M {sub ⋆}. The weaker [Fe/H] dependence shows that sub-Neptune planets, unlike giant planets, form readily in lower metallicity environment. The innermost region (a < 0.1 au) of an MMEN maintains a smooth profile despite a steep decline of planet occurrence rate: a result that favors the truncation of disks by corotating magnetospheres with a range of rotation periods, rather than the sublimation of dust. The Σ of Kepler multitransiting systems shows a much stronger correlation with M {sub ⋆} and [Fe/H] than singles. This suggests that the dynamically hot evolution that produced single systems also partially removed the memory of formation in disks. Radial-velocity planets yielded a MMEN very similar to CKS planets; transit-timing-variation planets’ postulated convergent migration history is supported by their poorly constrained MMEN. We found that lower mass stars have a higher efficiency of forming/retaining planets: for Sun-like stars, about 20% of the solid mass within ∼1 au are converted/preserved as sub-Neptunes, compared to 70% for late-K to early-M stars. This may be due to the lower binary fraction, lower giant-planet occurrence, or the longer disk lifetime of lower mass stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Enhanced superconductivity and electron correlations in intercalated ZrTe 3

We report charge density waves (CDWs) with superconductivity, competing Fermi surface instabilities, and collective orders have captured much interest in two-dimensional van der Waals (vdW) materials. Understanding the CDW suppression mechanism, its connection to the emerging superconducting state, and electronic correlations provides opportunities for engineering the electronic properties of vdW heterostructures and thin-film devices. Using a combination of the thermal transport, x-ray photoemission spectroscopy, Raman measurements, and first-principles calculations, we observe an increase in electronic correlations of the conducting states as the CDW is suppressed in ZrT e3 with 5% Cu and Ni intercalation in the vdW gap. As superconductivity emerges, intercalation brings not only decoupling of quasi-one-dimensional conduction electrons with phonons as a consequence of intercalation-induced lattice expansion but also a drastic increase in Zr 2+ at the expense of Zr 4+ metal atoms. These observations not only demonstrate the potential of atomic intercalates in the vdW gap for ground-state tuning but also illustrate the crucial role of the Zr metal valence in the formation of collective electronic orders.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Ab initio calculation of atomic solid hydrogen phases based on Gutzwiller many-body wave functions

We apply two ab initio many-body methods based on Gutzwiller wave functions, i.e., correlation matrix renormalization theory (CMRT) and Gutzwiller conjugate gradient minimization (GCGM), to the study of crystalline phases of atomic hydrogen. Both methods avoid empirical Hubbard U parameters and are free from double-counting issues. CMRT employs a Gutzwiller-type approximation that enables efficient calculations, while GCGM goes beyond this approximation to achieve higher accuracy at higher computational cost. By benchmarking against available quantum Monte Carlo (QMC) results, we demonstrate that while both methods are more accurate than the widely used density-functional theory, GCGM systematically captures additional correlation energy missing in CMRT, leading to significantly improved total energy predictions. We also show that by including the correlation energy Ec from local density approximation in the CMRT calculation, CMRT + E c produces energy in better agreement with the QMC results in these hydrogen lattice systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Core electron temperature turbulence and transport during sawtooth oscillations in the DIII-D tokamak

Abstract Sawteeth are one of the concerning instabilities in ITER and future burning plasma experiments. Sawtooth dynamics and its interaction with broadband plasma turbulence has been a challenge for predictive simulations of core transport in future fusion devices. This study provides new observations of core turbulence behavior during sawtooth oscillations in DIII-D hydrogen L-mode neutral beam injection heated plasmas in an inner wall limited configuration. A strong correlation of electron temperature and density turbulence levels with the sawtooth oscillation phase has been observed at locations inside the T e inversion radius and/or safety factor q = 1 magnetic surface. The T e turbulence amplitude in the core during the sawtooth ramp exhibits a critical T e gradient behavior inside but not near the T e inversion radius/ q = 1 magnetic surface. The most unstable mode calculated from the trapped gyro-landau fluid turbulence simulations reveal a change from low- k ion-type to low-k electron-type modes from pre- to post- sawtooth crash time periods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The Predictive Power of Exact Constraints and Appropriate Norms in Density Functional Theory

Ground-state Kohn-Sham density functional theory provides, in principle, the exact ground-state energy and electronic spin densities of real interacting electrons in a static external potential. In practice, the exact density functional for the exchange-correlation (xc) energy must be approximated in a computationally efficient way. About 20 mathematical properties of the exact xc functional are known. In this work, we review and discuss these known constraints on the xc energy and hole. By analyzing a sequence of increasingly sophisticated density functional approximations (DFAs), we argue that (a) the satisfaction of more exact constraints and appropriate norms makes a functional more predictive over the immense space of many-electron systems and (b) fitting to bonded systems yields an interpolative DFA that may not extrapolate well to systems unlike those in the fitting set. We discuss both how the class of well-described systems has grown along with constraint satisfaction and the possibilities for future functional development.

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

Monte Carlo Explicitly Correlated Second-Order Many-Body Green’s Function Calculations of Semiconductor Band Gaps

A systematically converging series of ab initio, post-density-functional, size-consistent, electron-correlated approximations is desired for predictive computing of felectronic band structures of insulating, semiconducting, and metallic solids. A series that meets all of these desiderata (except the applicability to metals) is ab initio many-body Green's function theory based on Gaussian-type-orbital (GTO) basis sets. Here, its leading-order approximation, the second-order Green's function (GF2) method in the diagonal and frequency-independent approximations with the aug-cc-pVDZ basis set, is applied to the fundamental band gaps of three semiconductors (diamond, silicon, and silicon carbide in the zincblende structure) using cluster models. Corrections are made to the basis-set-incompleteness errors by the explicit-correlation (F12) ansatz (GF2-F12) for the valence band edges. The crystals are modeled as surface-passivated clusters of increasing sizes, whose wave functions are expanded by up to 2709 GTO basis functions. Immense computational costs of these calculations are overcome by the highly scalable stochastic algorithm of the Monte Carlo GF2-F12 method, whose operation cost per state increases only as a cubic power of system size, which has a tiny memory footprint and easily achieves near-perfect parallel efficiency on thousands of CPUs or on hundreds of GPUs. The correlated, F12-corrected highest-occupied and lowest-unoccupied molecular-orbital energy (HOMO-LUMO) gap is 5.78 ± 0.07 eV for C 87 H 76 as compared with the experimental value of the fundamental (indirect) band gap of bulk diamond at 5.48 eV. The correlated, F12-corrected HOMO-LUMO gaps for Si 75 H 76 and Si 32 C 43 H 76 are 2.56 ± 0.15 eV and 3.50 ± 0.12 eV, respectively, which are expected to decrease further with increasing cluster sizes. As a result, the experimental fundamental (indirect) band gaps of bulk silicon and silicon carbide are 1.17 eV and 2.42 eV, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

SUNBIRD : a simulation-based model for full-shape density-split clustering

Combining galaxy clustering information from regions of different environmental densities can help break cosmological parameter degeneracies and access non-Gaussian information from the density field that is not readily captured by the standard two-point correlation function (2PCF) analyses. However, modelling these density-dependent statistics down to the non-linear regime has so far remained challenging. We present a simulation-based model that is able to capture the cosmological dependence of the full shape of the density-split clustering (DSC) statistics down to intra-halo scales. Our models are based on neural-network emulators that are trained on high-fidelity mock galaxy catalogues within an extended-ΛCDM framework, incorporating the effects of redshift-space, Alcock–Paczynski distortions, and models of the halo–galaxy connection. Our models reach sub-percent level accuracy down to $1 \, h^{-1}\text{Mpc}$ and are robust against different choices of galaxy–halo connection modelling. When combined with the galaxy 2PCF, DSC can tighten the constraints on ω cdm , σ 8 , and n s by factors of 2.9, 1.9, and 2.1, respectively, compared to a 2PCF-only analysis. DSC additionally puts strong constraints on environment-based assembly bias parameters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Interacting spin and charge density waves in the kagome metal FeGe

Unveiling the interplay between spin density wave (SDW) and charge density wave (CDW) orders in correlated electron materials is important in obtaining a comprehensive understanding of their electronic, structural, and magnetic properties. Kagome lattice materials are interesting because their flat electronic bands, Dirac points, and Van Hove singularities can enable a variety of exotic electronic and magnetic phenomena. The kagome metal FeGe (the B35 phase), which exhibits a CDW order deep within an A-type antiferromagnetic (AFM) phase, was found to respond dramatically to postgrowth annealing—with the ability to tune the CDW repeatedly from long-range order to negligible order. Additionally, neutron scattering studies suggest that incommensurate magnetic peaks that onset at 𝑇 Canting = 𝑇 SDW ≈ 60 K in the system arise from a SDW order instead of the AFM double-cone structure. Here, in this study, we use inelastic neutron scattering to show that two distinct spin excitations exist below 𝑇 Canting corresponding to two coexisting magnetic orders in the system in both sets of annealed samples with and without CDW. While CDW order or negligible order can dramatically affect the onset temperature of 𝑇 Canting and elastic incommensurate magnetic scattering, its impact on low-energy spin fluctuations is more limited. In both samples, a pair of gapless incommensurate spin excitations arising from the SDW order wave vector coexist with gapped commensurate spin waves from the A-type AFM order across 𝑇 Canting . The low-energy spin excitations for both samples couple dynamically to the lattice through enhanced magnetic scattering intensity on cooling below 𝑇 CDW , regardless of the status of the static long-range CDW order. The incommensurate SDW order in the long-range CDW ordered sample also induces a tiny in-plane lattice distortion of the kagome lattice that is absent in the negligible CDW ordered sample, in a way that is different from the previously known SDW and CDW ordering materials.

charge density waves↗
DOE OSTI2025

Achieving Multimodal and Multicolor Luminescence in LaAlO 3 :Pr 3+ , Gd 3+ via Trap Engineering and Energy Transfer

Achieving multimodal luminescence within a single phosphor is vital for multifunctional applications but remains challenging due to complex color tuning and trap engineering. In this study, we report Pr 3+ and Gd 3+ co‐doped LaAlO 3 (LAO:PG) phosphors, designed through careful modulation of multilevel traps and Pr 3+ → Gd 3+ energy transfer dynamics. These materials exhibit diverse luminescence modes, including down‐conversion luminescence (DCL), up‐conversion luminescence (UCL), persistent luminescence (PersL), optically stimulated luminescence (OSL), and thermally stimulated luminescence (TSL) across a wide spectral range. Unlike previously studied Pr 3+ ‐doped LAO, the co‐doped LAO:PG shows DCL in both UV‐visible and NIR regions and displays ultraviolet‐C UCL under visible excitation. Notably, we observe, for the first time, PersL lasting several minutes in these phosphors—an improvement over the non‐PersL behavior of Pr 3+ ‐only doped LAO. Additionally, the LAO:PG phosphors exhibit strong OSL response. TSL analysis reveals five distinct trap levels linked to these properties. Density functional theory calculations further correlate intrinsic defects to these traps, supporting a proposed mechanism for the observed multimodal luminescence. These findings highlight LAO:PG as a promising platform for developing advanced phosphors with integrated luminescence modes, paving the way for future applications in data storage, phototherapy, and anti‐counterfeiting technologies.

Chemistry↗
DOE OSTI2021

Hurricanes alter 10 Be concentrations in tropical river sediment but do not change regional erosion rate estimates

Tropical islands, including many in island arcs, are subjected to recurring disturbances from extreme storms such as tropical cyclones. To test whether such storms influence cosmogenic nuclide concentrations such that they do not reflect long-term rates of erosion, we measured meteoric and in situ 10Be in river sediment samples from Dominica, an andesitic island in the Caribbean, before and after category five Hurricane Maria (2017). Populations of before- and after-storm concentrations are statistically indistinguishable (n=7 pairs for in situ 10 Be, n=11 pairs for meteoric 10 Be). 10 Be concentrations vary from -138% to +73% within before-after sample pairs relative to the mean of the pair. These new data suggest that the effects of extreme storms on the depth and amount of near-surface erosion on Dominica vary spatially. Our data support calculations of Niemi, et al. (2005) and Yanites, et al. (2009) suggesting that basin-by-basin comparisons of erosion rates based on cosmogenic nuclides should be approached with caution in small (<~100 km 2 ) watersheds affected by mass movements and extreme storms. Erosion rates determined from in situ 10Be in Dominica (geometric mean=0.102 mm y-1, n=12) are low compared to similarly steep and wet areas globally and correlate positively with the spatial density of mass movements.

58 GEOSCIENCES↗
DOE OSTI2024

Stock-specific spatial overlap among seabird predators and Columbia River juvenile Chinook Salmon suggests a mechanism for predation during early marine residence

Abstract Objective Because predation is thought to be the primary source of natural mortality for juvenile salmon first entering the ocean, we sought to identify regions where, on average, stock-specific spatial overlap between the distribution of threatened and endangered juvenile Chinook Salmon Oncorhynchus tshawytscha and abundant fish-eating seabirds (common murres Uria aalge and sooty shearwaters Ardenna grisea) suggests the greatest potential for ocean predation risk to juvenile Chinook Salmon. Methods The relative abundance and spatial distribution of seabird predators and juvenile Chinook Salmon were quantified as part of long-term ecosystem surveys during May 2003–2012 and June 2003–2022. Genetic stock identification methods were used to assign individual fish to their respective stock groups. Stock-specific species distribution models then generated maps and indices of average annual spatial overlap between predators and prey within the survey area. Result There is unequivocal evidence for spatial overlap between common murres, sooty shearwaters, and five genetic groups of interior and lower Columbia River juvenile Chinook Salmon. We found strongly positive (≥0.70) spatial correlations between predator and prey densities in both May and June, although spatial overlap was, in general, greater during May. The region of highest spatial overlap occurred on the inner continental shelf between the Columbia River mouth (46.2°N) and Grays Harbor (47.0°N), a region at the beginning of the juvenile salmon migratory pathway that is strongly affected by freshwater outflow from the river. Conclusion Our findings support the idea that ocean avian predation during early marine residence has the potential to affect marine survival of juvenile Chinook Salmon and should be further investigated to better inform and implement ecological models and possible recovery actions for Chinook Salmon populations of the Columbia River basin.

Zamon, Jeannette E.↗
DOE OSTI2025

Enhancing low-temperature sintering in the MgO-LiF system: Mechanistic insights

In the present article, we provide compelling evidence that minor (1 wt%) additions of micron and nanometre-sized LiF particles in MgO decompose leading to free Li diffusing into MgO surfaces enhancing vacancy production during direct current sintering. The addition of nanometre-sized LiF particles leads to a > 500 °C reduction in the sintering temperature and over 99 % theoretical density of final consolidated compacts. Correlating differential scanning calorimetry with in-situ x-ray diffraction, together with Schottky calculations, post sintering electron microscopy and laser induced breakdown spectroscopy, we uncover critical insights into this impressive reduction in sintering temperature. Our quantitative analysis reveals that MgO and LiF interact at low temperatures with the diffusion of Li into the surface of MgO particles due to the intrinsic structural disorder of the LiF and MgO crystallites. Nanometre-sized LiF particles were found to react the most at low temperatures due to their increased structural disorder. Our multimodal characterization points to a Li-promoted densification and sintering mechanism. This impressive reduction in sintering temperature can be harnessed to promote low-temperature fabrication of MgO-based composites for technological applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Oxide dispersoid coherency of a ferritic-martensitic 12Cr oxide-dispersion-strengthened alloy under self-ion irradiation

Ferritic-martensitic oxide-dispersion-strengthened (ODS) alloy has shown excellent mechanical property and high radiation tolerance. However, the stability of dispersoids during displacive irradiation in the individual ferritic and martensitic phases is unclear. In this work, the correlation among dispersoid coherency, size, density, and matrix phase are studied in dual-phase 12Cr ODS after 100 peak displacements-per-atom (dpa) irradiation at 475°C, using 3.5 MeV Fe 2+ self-ions. The size and density changes of coherent and incoherent dispersoids were analyzed as a function of irradiation depth. The average dispersoid size decreased after irradiation in both phases, and the large incoherent dispersoids in the tempered martensite phase underwent a more dramatic change than those in the ferrite phase. The dispersoid density significantly increased in the ferrite phase within the irradiated region, mostly resulting in an increase in coherent dispersoid density. On the other hand, only a small change in density was observed in the tempered martensite phase for both coherent and incoherent dispersoids. This study shows that dispersoid evolution is dramatically different in ferrite and tempered martensite phases.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Defective Sn-Zn perovskites through bio-directed routes for modulating CO 2 RR

Sn and Zn oxides can catalyze electrochemical CO 2 reduction reactions (CO 2 RR) to produce useful chemicals. Herein, we report on perovskite-type tin-zinc oxides (TZO) synthesized by a biomineralization approach under benign aqueous conditions. The resulting TZO catalysts are disordered and possess a high quantity of defects, which serve as highly reactive active sites in CO 2 RR. Selectively toward formate or CO is tunable toward the relative amounts of Sn and Zn and their associated O defects. In-situ X-ray absorption spectroscopy measurements reveal electronic and structural transformations consistent with oxygen vacancy formation at Sn-O-Zn bridges as a function of cathodic potential that strongly correlates to CO2RR selectivity. Further, Density functional theory (DFT) calculations demonstrate that introducing an oxygen vacancy in TZO influences drastically influences formate/CO selectivity, consistent with our in-situ XAS findings and catalytic performance. Overall, these findings demonstrate the potential of biomineralization to fabricate perovskite-type metal oxides for energy conversion reactions.

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

Analysis of soot formation from aviation fuels in laminar counterflow flames

Combustion emissions from aviation contribute to the formation of condensation trail (contrail) that can lead to the formation of anthropogenic cirrus clouds. Ice particles that form contrails are observed to have a linear correlation with soot particle number density. Synthetic aviation fuels (SAFs) offer a promising route to mitigate the production of soot particles while also increasing energy security. Although studies have focused on combustion and spray behavior, the detailed investigation of soot formation processes for different jet fuels and their impact on models for computational fluid dynamics (CFD) applications is not well understood. Moreover, experimental measurements of soot for canonical flames using Synthetic aviation fuels (SAF) for model validation remain scarce. To address this, we use employed the Lawrence Livermore National Laboratory (LLNL) detailed soot model based on the discrete sectional method. Additionally, we develop two reduced chemical mechanisms for Jet-A and Alcohol-to-Jet (C1) that are suitable for turbulent flame simulations and couple them with the Hybrid Method of Moments (HMOM). The detailed and reduced model frameworks are validated against experimental measurements of soot volume fraction (ƒ ν ) from a counterflow burner experiment previously reported in the literature. Given the good agreement between modeling results and experimental measurements for the (1) spatial distribution of ƒ ν and (2) the non-linear variation of peak ƒ ν with strain rate, we further investigate the modeled sub-processes (nucleation, condensation, surface growth, and oxidation) using the LLNL model to analyze the assumptions in the reduced model framework. Furthermore, the results indicate a significant contribution from resonant radicals to the surface growth of soot particles, which are not accounted for in the current implementation of HMOM and could help reconcile soot predictions by the reduced model with observations.

Counterflow↗
DOE OSTI2023

Effect of Linker Structure and Functionalization on Secondary Gas Formation in Metal–Organic Frameworks

Rare-earth terephthalic acid (BDC)-based metal–organic frameworks (MOFs) are promising candidate materials for acid gas separation and adsorption from flue gas streams. However, previous simulations have shown that acid gases (H 2 O, NO 2 , and SO 2 ) react with the hydroxyl on the BDC linkers to form protonated acid gases as a potential degradation mechanism. Herein, gas-phase computational approaches were used to identify the formation energies of these secondary protonated acid gases across multiple BDC linker molecules. Formation energies for secondary protonated acid gases were evaluated using both density functional theory (DFT) and correlated wave function methods for varying BDC–gas reaction mechanisms. Upon validation of DFT to reproduce wave function calculation results, rotated conformational linkers and chemically functionalized BDC linkers with -OH, -NH 2 , and -SH were investigated. Furthermore, the calculations show that the rotational conformation affects the molecule stability. Double-functionalized BDC linkers, where two functional groups are substituted onto BDC, showed varied reaction energies depending on whether the functional groups donate or withdraw electrons from the aromatic system. Based on these results, BDC linker design must balance adsorption performance with degradation via linker dehydrogenation for the design of stable MOFs for acid gas separations.

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

Hydrolysis Reaction Pathways of Thorium Oxide Nanoclusters

Density functional theory benchmarked by correlated molecular orbital theory is used to develop a fundamental and predictive understanding of the interaction of thorium oxide nanoclusters with gas phase water to provide insight into nuclear-waste storage, production of thorium nuclear reactor fuels, and reprocessing of spent fuel. The structures of Th n O 2n (n = 3 – 6) clusters and their interactions with water have been studied at the B3LYP, MP2, and CCSD(T) levels. Hydrolysis is initiated by the formation of Lewis acid-base adducts, with relative H 2 O binding energies (physisorption) ranging from −15 kcal/mol to −22 kcal/mol. The initial H 2 O physisorption energy is ca. −21 kcal/mol regardless of the cluster size and is consistent with the experimentally obtained initial adsorption energy on a thorium dioxide surface. The physisorption enthalpies for additional water molecules can be affected by the presence of terminal groups OH groups generated by proton transfer to a Th-O near the site of adsorption. The hydrolysis products (chemisorption) form either bridging or terminal hydroxides. More exothermic hydrolysis steps were predicted for the formation of terminal hydroxides as compared to the formation of bridging hydroxides. Here, the calculated transition state barriers for transfer of protons from bound water complexes to form the chemisorption products are very low. Overall, water readily reacts with thorium oxide clusters preferring hydroxide products over hydrated complexes. First and second order fits were predicted for the combined physisorption and chemisorption energies for the hydrolysis of thorium oxide clusters. Finally, ionization energies and electron affinities were calculated as were HOMO-LUMO gaps to provide additional insights into the properties of the thorium oxide and hydroxide clusters.

Adsorption↗