Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DFT”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Hybrid Quantum Mechanical, Molecular Mechanical, and Machine Learning Potential for Computing Aqueous-Phase Adsorption Free Energies on Metal Surfaces

Performing reliable computer simulations of elementary processes occurring at metal–water interfaces is pivotal for novel catalyst design in sustainable energy applications. Computational catalyst design hinges on the ability to reliably and efficiently compute the potential energy surface (PES) of the system. Here, due to the large system sizes needed for studying processes at liquid water–metal interfaces, these systems can currently not be described using density functional theory (DFT). In this work, we used a hybrid quantum mechanical, molecular mechanical, and machine learning potential for studying the adsorption behavior of phenol, atomic hydrogen, 2-butanol, and 2-butanone on the (0001) facet of Ru under reducing conditions when Ru is not oxidized. Specifically, we describe the adsorbate and the surrounding metal atoms at the DFT level of theory. Here, we also considered the electrostatic field effect of the water molecules on adsorbate–metal interactions. Next, for the water–water and water–adsorbate interactions, we used established classical force fields. Finally, for the water–Ru surface interaction, for which no reliable force fields have been published, we used Behler–Parrinello high-dimensional neural network potentials (HDNNPs). Employing this setup, we used our explicit solvation for metal surface (eSMS) approach to compute the aqueous-phase effect on the low-coverage adsorption of selected molecules and atoms on the (0001) facet of Ru. In agreement with previous experimental and computational studies of oxygenated molecules over transition metal facets, we found that liquid water destabilizes the tested adsorbates on Ru(0001). Interestingly, our findings indicate that adsorbates on Ru are less affected by the presence of an aqueous phase than on other transition metals (e.g., Pt), highlighting the necessity of experimental investigations of Ru-based catalytic systems in liquid water.

Adsorption↗
DOE OSTI2025

Identifying Band Inversions in Topological Materials Using Diffusion Monte Carlo

Topological insulators are characterized by insulating bulk states and robust metallic surface states. Band inversion is a hallmark of topological insulators. At time-reversal invariant points in the Brillouin zone, spin–orbit coupling (SOC) induces a swapping of orbital character at the bulk band edges. Reliably detecting band inversion in solid-state systems with many-body methods would aid in identifying possible candidates for spintronics and quantum computing applications and improve our understanding of the physics behind topologically nontrivial systems. Density functional theory (DFT) methods are a well-established means of investigating these interesting materials due to their favorable balance of computational cost and accuracy but often struggle to accurately model the electron–electron correlations present in the many materials containing heavier elements. In this work, we develop a novel method to detect band inversion within continuum quantum Monte Carlo (QMC) methods that can accurately treat the electron correlation and spin–orbit coupling that are crucial to the physics of topological insulators. Our approach applies a momentum-space-resolved atomic population analysis throughout the first Brillouin zone utilizing the Löwdin method and the one-body reduced density matrix produced with diffusion Monte Carlo (DMC). We integrate this method into QMCPACK, an open source ab initio QMC package, so that these ground-state methods can be used to complement experimental studies and validate prior DFT work on predicting the band structures of correlated topological insulators. Here, we demonstrate this new technique on the topological insulator bismuth telluride, which displays band inversion between its Bi-p and Te-p states at the Γ-point. We show an increase in charge on the bismuth-p orbital and a decrease in charge on the tellurium-p orbital when comparing band structures with and without SOC. Additionally, we use our method to compare the degree of band inversion present in monolayer Bi 2 Te 3 , which has no interlayer van der Waals interactions, to that seen in the bilayer and bulk. The method presented here will enable future many-body studies of band inversion that can shed light on the delicate interplay between correlation and topology in correlated topological materials.

Band structure↗
DOE OSTI2026

Characterizing Defect Dynamics in Silicon Carbide Using Symmetry-Adapted Collective Variables and Machine Learning Interatomic Potentials

Silicon carbide (SiC) divacancies are attractive candidates for spin-defect qubits possessing long coherence times and optical addressability. The high activation barriers associated with SiC defect formation and motion pose challenges for their study by first-principles molecular dynamics. In this work, we develop and deploy machine learning interatomic potentials (MLIPs) to accelerate defect dynamics simulations while retaining ab initio accuracy. We employ an active learning strategy comprising symmetry-adapted collective variable discovery and enhanced sampling to compile configurationally diverse training data, calculation of energies and forces using density functional theory (DFT), and training of an E(3)-equivariant MLIP based on the Allegro model. Here, the trained MLIP reproduces DFT-level accuracy in defect transition activation free energy barriers, enables the efficient and stable simulation of multidefect 216-atom supercells, and permits an analysis of the temperature dependence of defect thermodynamic stability and formation/annihilation kinetics to propose an optimal annealing temperature to maximally stabilize VV divacancies.

Computer simulations↗
DOE OSTI2020

A Multiple Filter Based Neural Network Approach to the Extrapolation of Adsorption Energies on Metal Surfaces for Catalysis Applications

Computational catalyst discovery involves the development of microkinetic reactor models based on estimated parameters determined from density functional theory (DFT). For complex surface chemistries, the number of reaction intermediates can be very large, and the cost of calculating the adsorption energies by DFT for all surface intermediates even for one active site model can become prohibitive. In this paper, we have identified appropriate descriptors and machine learning models that can be used to predict a significant part of these adsorption energies given data on the rest of them. Furthermore, our investigations also included the case when the species data used to train the predictive model are of different size relative to the species the model tries to predict—this is an extrapolation in the data space which is typically difficult with regular machine learning models. Due to the relative size of the available data sets, we have attempted to extrapolate from the larger species to the smaller ones in the current work. Here, we have developed a neural network based predictive model that combines an established additive atomic contribution based model with the concepts of a convolutional neural network that, when extrapolating, achieves a statistically significant improvement over the previous models.

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

Excited State Orbital Optimization via Minimizing the Square of the Gradient: General Approach and Application to Singly and Doubly Excited States via Density Functional Theory

We present a general approach to converge excited state solutions to any quantum chemistry orbital optimization process, without the risk of variational collapse. The resulting square gradient minimization (SGM) approach only requires analytic energy/Lagrangian orbital gradients and merely costs 3 times as much as ground state orbital optimization (per iteration), when implemented via a finite difference approach. SGM is applied to both single determinant ΔSCF and spin-purified restricted open-shell Kohn-Sham (ROKS) approaches to study the accuracy of orbital optimized DFT excited states. It is found that SGM can converge challenging states where the maximum overlap method (MOM) or analogues either collapse to the ground state or fail to converge. We also report that ΔSCF/ROKS predict highly accurate excitation energies for doubly excited states (which are inaccessible via TDDFT). Singly excited states obtained via ROKS are also found to be quite accurate, especially for Rydberg states that frustrate (semi)local TDDFT. Our results suggest that orbital optimized excited state DFT methods can be used to push past the limitations of TDDFT to doubly excited, charge-transfer, or Rydberg states, making them a useful tool for the practical quantum chemist's toolbox for studying excited states in large systems.

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

Automation of Active Space Selection for Multireference Methods via Machine Learning on Chemical Bond Dissociation

Predicting and understanding the chemical bond is one of the major challenges of computational quantum chemistry. Kohn–Sham density functional theory (KS-DFT) is the most common method, but approximate density functionals may not be able to describe systems where multiple electronic configurations are equally important. Multiconfigurational wave functions, on the other hand, can provide a detailed understanding of the electronic structures and chemical bonds of such systems. In the complete active space self-consistent field (CASSCF) method, one performs a full configuration interaction calculation in an active space consisting of active electrons and active orbitals. However, CASSCF and its variants require the selection of these active spaces. This choice is not black box; it requires significant experience and testing by the user, and thus active space methods are not considered particularly user-friendly and are employed only by a minority of quantum chemists. Our goal is to popularize these methods by making it easier to make good active space choices. Here, we present a machine learning protocol that performs an automated selection of active spaces for chemical bond dissociation calculations of main group diatomic molecules. The protocol shows high prediction performance for a given target system as long as a properly correlated system is chosen for training. Good active spaces are correctly predicted with a considerably better success rate than random guess (larger than 80% precision for most systems studied). Our automated machine learning protocol shows that a “black-box” mode is possible for facilitating and accelerating the large-scale calculations on multireference systems where single-reference methods such as KS-DFT cannot be applied.

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

Domain Separated Density Functional Theory for Reaction Energy Barriers and Optical Excitations

We recently proposed Domain Separated Density Functional Theory (DS-DFT), a framework that allows for the combination of different levels of theory for the computation of the electronic structure of molecules. Here, this work discusses the application of DS-DFT to the computation of transition-state energy barriers and optical absorption spectra. We considered several hydrogen abstraction reactions and optical spectra of molecule/metal cluster systems, including the absorption of individual species such as carbon monoxide, methane, and molecular hydrogen to a Li6 cluster. We present and discuss two domain-separated methods: i), the screened-density approximation (SDA), and, ii), linearly weighted exchange (LWE). We find that SDA, which is applied as a hybridization based on atomic domains, could be useful to computing energy barriers, whereas LWE is suited for the analysis of electronic properties such as ground-state gaps, excitation energies, and oscillator strengths.

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

Analysis of Intermediates and Products from the Dehydrogenation of Mg(BH 4 ) 2

The thermodynamic properties of key compounds, Mg(B 3 H 8 ) 2 , MgB 2 H 6 , MgB 10 H 10 , Mg(B 11 H 14 ) 2 , Mg 3 (B 3 H 6 ) 2 , and MgB 12 H 12 , proposed to be formed in the release of hydrogen from magnesium borohydride Mg(BH 4 ) 2 and uptake of hydrogen by MgB 2 , have been investigated using solid–state density functional theory (DFT) calculations. Further, more accurate treatment of cell–size effects to the entropies were also investigated, in order to improve the accuracy of the thermodynamic properties of complex borohydrides. We find that the zero–point energy corrections can lower the electronic energies of reaction by 20–30 kJ/(mol H 2 ) for these intermediates, while adding the thermal and entropy contributions results in the total decrease up to ~50 kJ/(mol H 2 ). Although our treatment lowers the calculated formation energy of Mg(B 3 H 8 ) 2 , it is still too high to explain the experimental observation of B 3 H 8 – . We discuss possible reasons for this disparity and propose that the formation of B 3 H 8 – and H – in a disordered amorphous phase has a large energy difference compared to the phase–separated Mg(B 3 H 8 ) 2 and MgH 2 considered in calculations. Comparison of the experimental and NMR chemical shifts calculated within a DFT approach for known species Mg(BH 4 ) 2 , Mg(B 3 H 8 ) 2 , Mg(B 11 H 14 ) 2 , MgB 10 H 10 and MgB 12 H 12 provides validation for predicting the chemical shifts of the other compounds which are yet to be confirmed experimentally. These include MgB 2 H 6 and the proposed tri–anion species Mg 3 (B 3 H 6 ) 2 that both have favorable thermodynamics for reversible hydrogen storage in Mg(BH 4 ) 2 without the formation of MgH 2 as a co–product which could phase–separate and inhibit rehydrogenation.

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

Computational Predictions of the Hydrolysis of 2,4,6-Trinitrotoluene (TNT) and 2,4-Dinitroanisole (DNAN)

Hydrolysis is a common transformation reaction that can affect the environmental fate of many organic compounds. In this study, three proposed mechanisms of alkaline hydrolysis of 2,4,6-trinitrotoluene (TNT) and 2,4-dinitroaniline (DNAN) were investigated with plane-wave density functional theory (DFT) combined with ab initio and classical molecular dynamics (AIMD/MM) free energy simulations, Gaussian basis set DFT calculations, and correlated molecular orbital theory calculations. Most of the computations in this study were carried out using the Arrows web based tools. For each mechanism: Meisenheimer complex formation, nucleophilic aromatic substitution, and proton abstraction reaction energies and activation barriers were calculated for the reaction at each relevant site. For TNT, it was found that the most kinetically favorable first hydrolysis steps involve Meisenheimer complex formation by attachment of OH$^{-}$ at the C1 and C3 arene carbons and proton abstraction from the methyl group. The nucleophilic aromatic substitution reactions at the C2 and C4 arene carbons were found to be thermodynamically favorable. However, the calculated activation barriers were slightly lower than previous studies, but still found to be $\Delta$G$^{\ddagger} \approx 18$ kcal/mol using PBEo AIMD/MM free energy simulations, suggesting the reactions are not kinetically significant. For DNAN, the barriers of nucleophilic aromatic substitution were even greater ($\Delta$G$^{\ddagger} > 29$ kcal/mol PBEo AIMD/MM). Further, the most favorable hydrolysis reaction for DNAN was found to be a two-step process in which the hydroxyl first attacks the C1 carbon to form a Meisenheimer complex at the C1 arene carbon C1-(OCH$_3$)OH$^{-}$, subsequently, the methoxy anion (-OCH$_3$) at the C1 arene carbon dissociates and the proton shuttles from the C1-OH to the dissociated methoxy group, resulting in methanol and an aryloxy anion.

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

Density Matrix Implementation of the Fermi–Löwdin Orbital Self-Interaction Correction Method

The Fermi–Löwdin orbital self-interaction correction (FLOSIC) method effectively provides a transformation from canonical orbitals to localized Fermi–Löwdin orbitals which are used to remove the self-interaction error in the Perdew–Zunger (PZ) framework. This transformation is solely determined by a set of points in space, called Fermi–Löwdin descriptors (FODs), and the occupied canonical orbitals or the density matrix. In this work, we provide a detailed workflow for the implementation of the FLOSIC method for removal of self-interaction error in DFT calculations in an orbital-by-orbital basis that takes advantage of the unitary invariant nature of the FLOSIC method. In this way, it is possible to cast the self-consistent energy minimization at fixed FODs in the same manner than standard Kohn–Sham with one additional term in the Kohn–Sham Hamiltonian that introduces the PZ self-interaction correction. Each energy minimization iteration is divided in two substeps, one for the density matrix and one for the FODs. Expressions for the effective Kohn–Sham matrix and FOD gradients are provided such that its implementation is suitable for most electronic structure codes. Here, we analyze the convergence characteristics of the algorithm and present applications for the evaluation of NMR shielding constants and real-time time-dependent DFT simulations based on the Liouville–von Neumann equation to calculate excitation energies.

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

Importance of Dispersion in the Molecular Geometries of Mn(III) Spin-Crossover Complexes

The computational investigation of the molecular geometries of a pair of manganese(III) spin-crossover complexes is reported. For the geometry of the quintet high-spin state, density functionals significantly overestimate Mn–Namine bond distances, although the geometry for the triplet intermediate-spin state is well described. Here, comparisons with several wave function-based methods demonstrate that this error is due to the limited ability of commonly used density functionals to recover dispersion beyond a certain extent. Among the methods employed for geometry optimization, restricted open-shell Møller–Plesset perturbation theory (MP2) appropriately describes the high-spin geometry but results in a slightly shorter Mn–O distance in both spin states. On the other hand, extended multistate complete active space second-order perturbation theory (XMS-CASPT2) provides a good description of the geometry for the intermediate-spin state but also sufficiently recovers dispersion, performing well for the high-spin state. Despite the fact that the electronic structure of both spin states is dominated by one-electron configuration, XMS-CASPT2 offers a balanced approach, leading to molecular geometries with much better agreement with experiment than MP2 and DFT. A scan along the Mn–N amine bond demonstrates that for these complexes coupled cluster methods (i.e., DLPNO-CCSD(T)) also yield bond distances in agreement with experiment while multiconfiguration pair density functional theory (MC-PDFT) is unable to recover dispersion well enough, analogous to single-reference DFT.

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

Energetic and Spectroscopic Properties of Astrophysically Relevant MgC 4 H Radicals Using High-Level Ab Initio Calculations

Considering the importance of magnesium-bearing hydrocarbon molecules (MgC n H; n = 2, 4, and 6) in the carbon-rich circumstellar envelopes (e.g., IRC+10216), a total of 28 constitutional isomers of MgC 4 H have been theoretically investigated using density functional theory (DFT) and coupled-cluster methods. The zero-point vibrational energy corrected relative energies at the ROCCSD(T)/cc-pCVTZ level of theory reveal that the linear isomer, 1-magnesapent-2,4-diyn-1-yl (1, 2 Σ + ), is the global minimum geometry on the MgC 4 H potential energy surface. The latter has been detected both in the laboratory and in the evolved carbon star, IRC+10216. The calculated spectroscopic data for 1 match well with the experimental observations (error ~ 0.78%) which validates our theoretical methodology. Plausible isomerization processes happening among different isomers are examined using DFT and coupled-cluster methods. CASPT2 calculations have been performed for a few isomers exhibiting multireference characteristics. The second most stable isomer, 1-ethynyl-1λ 3 -magnesacycloprop-2-ene-2,3-diyl (2, 2 A 1 , μ = 2.54 D), is 146 kJ mol –1 higher in energy than 1 and possibly the next promising candidate to be detected in the laboratory or in the interstellar medium in future.

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

Evaluating the Importance of Conformers for Understanding the Vacuum-Ultraviolet Spectra of Oxiranes: Experiment and Theory

Vacuum-ultraviolet (VUV) absorption spectroscopy enables electronic transitions that offer the unambiguous identification of molecules. As target molecules become more complex, multifunctional species present a great challenge to both experimental and computational spectroscopy. This research reports both experimental and theoretical studies of oxiranes. Computationally, the nuclear ensemble approach has been used to accurately predict experimental spectra for a variety of molecules. However, this approach incurs great computational cost, as ensembles generally consist of thousands of geometries. The present study aims to drastically reduce the ensemble by evaluating the significance of the conformers to the predicted spectra. This approach was applied to 11 substituted oxiranes using the Conformer Rotamer Ensemble Sampling Tool (CREST) of Grimme to generate an ensemble of unique conformers determined by their Boltzmann populations. Five TD-DFT functionals (BMK, CAM-B3LYP, M06-2X, MN15, ωB97X-D) and EOM-CCSD were used to simulate the spectrum of each substituted oxirane ensemble. Computed spectra were then compared to the experiment using both qualitative and quantitative metrics. Based on these metrics, it was observed that certain conformers may not be necessary to characterize this set of oxiranes despite the temperature (323 K) of the experiment. A single conformer can then be used with TD-DFT and EOM-CCSD to replicate the experimental spectra of these medium-sized combustion species.

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

Prediction of Redox Potentials for Different Oxidation States of U, Np, Pu, and Am in Alkaline Aqueous Solution

The redox potentials for U, Np, Pu, and Am for oxidation states +III up to +VIII in alkaline aqueous solutions were predicted using density functional theory (DFT) and small-core pseudopotentials and their basis sets, with a hybrid explicit/implicit solvent model using SHE = 4.28 V. For each oxidation state, various oxo/hydroxo complexes were evaluated resulting in a variety of one-electron redox pathways. For An(VIII/VII) couples, the predicted redox potentials for the [An(VIII)O 5 (OH)] -3 /[An(VII)O 4 (OH) 2 ] -3 , or [An(VIII)O 4 (OH) 2 ] -2 /[An(VII)O 4 (OH) 2 ] -3 couples are in good agreement with existing estimates. For An(VII/VI) redox couples, all couples, particularly [An(VII)O 4 (OH) 2 ] -3 /[An(VI)O 2 (OH) 4 ] -2 , were in agreement with experimental values for U, Np, and Pu, but the results for Am showed larger differences from the estimated potentials. The An(VI/V) couples were consistent with experiment for dioxo/tetrahydroxo couples, and the An(V/IV) couples showed acceptable agreement based on actinide-specific couples, with neutral hydroxides often favored in the +IV state. The An(IV/III) couples were consistent with literature values when modeled as soluble neutral hydroxides. The use of our approach yielded calculated redox potentials that were within ±0.2 V from experimental or estimated values consistent with our prior calculations on redox potentials of actinides from Ac to Am in acidic aqueous solutions. Finally, as a result, this supports the robustness of our DFT-based methodology for predicting actinide redox potentials, offering valuable insights for actinide chemistry in aqueous solutions.

Actinides↗
DOE OSTI2025

Photochemistry of Hypervalent Iodoazide Derivatives

The photochemical properties and reaction mechanisms of a series of hypervalent iodoazide compounds (R–IN3) were investigated, with substituents (−CH 3 , −H, −CF 3 ) tuning the electronic density of the phenyl ring. With the help of UV irradiation, ultrafast time-resolved spectroscopy, and density functional theory (DFT) calculations, we elucidated the mechanisms of azide radical (N 3 • ) release and its subsequent reactivity. UV–vis spectroscopy reveals that the photoconversion rates follow the trend CH3 > H > CF3, aligning with DFT-calculated ΔG values for ring-opening transitions. Homolytic cleavage of the I–N bond is identified as the dominant pathway for N3• generation upon UV irradiation, occurring within 400 fs. The released azide radicals react with the solvent molecules, while the iodo radical R–I C • fragments undergo a thermodynamically uphill lactone ring-opening (RO) reaction and subsequent hydrogen atom abstraction from the solvent to form carboxylic acids R–I–COOH, as validated by NMR and IR spectroscopy. The study also highlights the role of substituents in influencing reaction kinetics and intermediate stability, with electron-donating groups accelerating the release of the N 3 • species. This work bridges experimental observations with computational predictions, offering a foundation for future advancements in azide-based reactions and materials.

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

Activation, Dehydrogenation, and Carbon–Carbon Coupling of Methane by Iridium Cations Studied by Infrared Multiple Photon Dissociation Spectroscopy and Density Functional Theory

Products resulting from the sequential activation of one, two, three, and four methane molecules by atomic iridium cations were characterized by gas-phase infrared multiple photon dissociation spectroscopy and density functional theory (DFT) calculations. Iridium cations were generated using a laser ablation source and reacted with methane in a linear radiofrequency ion trap before mass analysis and spectroscopic interrogation in a Fourier transform ion cyclotron resonance mass spectrometer coupled to the free-electron laser for intracavity experiments (FELICE) beamline. Product ions were irradiated using infrared light over the 250–1500 cm –1 range. Comparisons between the experimental and DFT-calculated spectra enabled structural determination of the products formed. The observed products include HIrCH + , Ir(CH 2 ) 2 + , H s Ir(C 3 H 5 ) + , and Ir(CH 3 ) s (C 3 H 5 ) + , where the subscript s denotes a syn orientation of the two ligands. Furthermore, formation of the latter two products provides evidence for efficient C–H bond activation and subsequent C–C coupling on the atomic iridium cation.

Chemical reactions↗
DOE OSTI2026

Accurate Prediction of pKb in Amines: Validation of the CAM-B3LYP/6-311+G(d,p)/SMD Model

Amines play several key roles in chemistry and biology and are involved in numerous industrial processes, often with significant economic impacts. Recently, amines are also garnering interest as catalysts for polymer synthesis and for CO 2 fixation, incentivizing the need to rapidly design and screen new amino compounds. Hence, developing reliable methods to predict their physicochemical properties, e.g., the base dissociation constant (pKb), is pivotal. Here, a density functional theory (DFT)-based approach was employed to compute the pKb of substituted amines, exploring the impact of several key parameters, including (i) the number of explicit water molecules at the reaction center, (ii) the van der Waals (vdW) surface, and (iii) solvent polarizability. In previous work, it was determined that including two explicit water molecules at the reaction center resulted in highly accurate pKb estimates for primary amines. Here, we find that including a third water molecule at the reaction center is essential for accurate pKb for secondary and tertiary amines. The revised methodology was then applied to a wider selection of amines, obtaining a minimum average error (MAE) < 0.4. In conclusion, this result represents an extension of our “easy-to-use method,” a simple and direct DFT approach exploiting CAM-B3LYP/SMD/6-311G+(d,p) to compute pKb without post facto modifications.

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

Thorium Monosilicide, ThSi: An Experimental and Theoretical Study

The present theoretical and experimental combination study investigates the ThSi molecule in detail. Computationally, we utilized high-level multireference and coupled-cluster levels of theory conjoined with large correlation consistent basis sets to study a series of electronic and spin–orbit states of ThSi. Here, we report potential energy curves (PECs), electron configurations at equilibrium distances, spectroscopic constants, energetics, and spin–orbit coupling effects for 16 electronic states of ThSi. The studied 16 electronic states are arranged tightly within 0.9 eV, highlighting the complexity of the electronic spectrum of ThSi. The ground electronic state of ThSi is a single-reference 1 1 Σ + state that derives from the 1σ 2 2σ 2 1π 4 electronic configuration. The Ω = 0 + spin–orbit ground state of ThSi is composed of 1 1 Σ + (47%) and 13Π (44%) electronic states. Our measured bond energy (D0) of ThSi, obtained using resonant two-photon ionization (R2PI) spectroscopy is 3.146(4) eV, where the assigned error limit is given in parentheses in units of the last quoted digits. The computed D0 of ThSi (Ω = 0 + ) at the CBS-C-CCSD(T)-δT(Q)-δDK-δSO level (3.181 eV) is in good agreement with the experimental value. Our derived enthalpy of formation for ThSi, Δ f H 0K o (ThSi(g)), is 971.8(6.0) kJ/mol. Finally, we have performed density functional theory (DFT) calculations for ThSi(1 1 Σ + ) using 16 exchange correlation functionals that span multiple rungs of “Jacob’s ladder” of density functional approximation (DFA) to assess the DFT errors on D 0 , r e , and ω e of ThSi with respect to experimental and ab initio coupled-cluster values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗