Search NASA⌕ Search

SEARCH · Search NASA

Results for “CBS”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Geometric Energy Derivatives at the Complete Basis Set Limit: Application to the Equilibrium Structure and Molecular Force Field of Formaldehyde

Here, geometric energy derivatives which rely on core-corrected focal-point energies extrapolated to the complete basis set (CBS) limit of coupled cluster theory with iterative and noniterative quadruple excitations, CCSDTQ and CCSDT(Q), are used as elements of molecular gradients and, in the case of CCSDT(Q), expansion coefficients of an anharmonic force field. These gradients are used to determine the CCSDTQ/CBS and CCSDT(Q)/CBS equilibrium structure of the S 0 ground state of H 2 CO where excellent agreement is observed with previous work and experimentally derived results. A fourth-order expansion about this CCSDT(Q)/CBS reference geometry using the same level of theory produces an exceptional level of agreement to spectroscopically observed vibrational band origins with a MAE of 0.57 cm –1 . Second-order vibrational perturbation theory (VPT2) and variational discrete variable representation (DVR) results are contrasted and discussed. Vibration–rotation, anharmonicity, and centrifugal distortion constants from the VPT2 analysis are reported and compared to previous work. Additionally, an initial application of a sum-over-states fourth-order vibrational perturbation theory (VPT4) formalism is employed herein, utilizing quintic and sextic derivatives obtained with a recursive algorithmic approach for response theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Benchmark Electronic Structure Calculations for H 3 O + (H 2 O) n , n = 0–5, Clusters and Tests of an Existing 1,2,3-Body Potential Energy Surface with a New 4-Body Correction

We report a revision of a previously developed many-body (up to 3-body inter- actions) CCSD(T)-based potential energy surface (PES) for the hydrated proton and its extension that incorporates a simple 4-body term for the hydronium-water-water-water interactions. The new PES is parametrized from and tested against new extensive benchmark CCSD(T) Complete Basis Set (CBS) estimates of the binding energies of H3O+(H2O)n clusters, n=0-5, which include all known to date low-lying energy isomers. We report the energetics, structures and harmonic frequencies of the hydronium-water clusters using the PES and compare them against the new benchmarks as well as other best available ab initio results. This critical comparison demonstrates the high accuracy of the new many-body hydronium-water PES.

Heindel, Joseph P.↗
DOE OSTI2020

The many-body expansion for aqueous systems revisited: 1. Water-water interactions

We revisit the Many-Body Expansion (MBE) for water-water interactions by examining the effects of the basis set, including those resulting from the Basis Set Superposition Error (BSSE) correction, on the various terms for selected sizes of water clusters up to n = 21. The analysis is performed at the second order Møller-Plesset (MP2) perturbation theory with the family of augmented correlation consistent basis sets up to five zeta quality for the (H2O)n, n = 7, 10, 13, 16 and 21 clusters, for which we report either the complete MBE (n = 7, 10) or the one through the 6-body (n = 13) and the 5-body terms (n = 16, 21). Our results suggest that any sizeable contributions to the total cluster binding energy arising from the 5-body and larger terms are solely an artifact of the size of the basis set. Indeed, all terms above the 4-body converge monotonically to practically zero at the Complete Basis Set (CBS) limit. The magnitudes of the 5-body and above terms are accurately reproduced even with the smaller basis set of the series (aug-cc-pVDZ) once the BSSE correction is taken into account. The same level of theory (MP2/aug-cc-pVDZ, BSSEcorrected) also accurately reproduces the MP2/CBS values of the 3- and 4-body terms. The contribution of electron correlation to the 3- and 4-body terms is quite small, so that neglecting the correlation contribution in all terms above the 3-body results in an error of ??(0.1%). The BSSE correction to the largest 2-body term in the MBE was accurately estimated from the function ?? [1 + erf(-?? ??)], which is proportional to the common (overlapping) area between two Gaussian distributions whose centers are separated by R, with the constants ??, ?? fitted to the calculated BSSE corrections for the individual 2-body terms of the clusters with each basis set. Our results unquestionably demonstrate that the MBE for water-water interactions converges to the 4-body term, since any finite terms above the 4-body are artifacts of the size of the basis set and can be safely truncated at the 4-body term when either a very large basis set is used or BSSE corrections are taken into account even with the aug-cc-pVDZ basis set. We expect these findings to have important consequences in the pursuit of accurate many-body molecular dynamics simulations for aqueous systems.

Heindel, Joseph↗
DOE OSTI2023

Coupled Cluster Study of the Heats of Formation of UF 6 and the Uranium Oxyhalides, UO 2 X 2 (X = F, Cl, Br, I, and At)

The atomization enthalpies of the U(VI) species UF6 and the uranium oxyhalides UO 2 X 2 (X=F, Cl, Br, I, At) were calculated using a composite relativistic Feller-Peterson-Dixon (FPD) approach based on scalar relativistic DKH3-CCSD(T) with extrapolations to the CBS limit. The inherent multideterminantal nature of the U atom was mitigated by utilizing the singly charged atomic cation in all calculations with correction back to the neutral asymptote via the accurate ionization energy of the U atom. The effects of SO coupling were recovered using full 4-component CCSD(T) with contributions due to the Gaunt Hamiltonian calculated using Dirac-Hartree-Fock. The final atomization enthalpy for UF 6 (752.2 kcal/mol) was within 2.5 kcal/mol of the experimental value, but unfortunately the latter carries a ±2.4 kcal/mol uncertainty that is predominantly due to the experimental uncertainty in the formation enthalpy of U atom. The analogous value for UO 2 F 2 (607.6 kcal/mol) was in nearly exact agreement with experiment, but the latter has a stated experimental uncertainty of ±4.3 kcal/mol. The FPD atomization enthalpy for UO 2 Cl 2 (540.4 kcal/mol) was within the experimental error limits of ±5.5 kcal/mol. FPD atomization energies for the non-U-containing molecules (used for reaction enthalpies) H 2 O and HX (X=F, Cl, Br, I, At) were within at most 0.3 kcal/mol of their experimental values where available. The FPD atomization enthalpies, together with FPD reaction enthalpies for two different reactions, were used to determine heats of formation for all species of this work with estimated uncertainties of ±4 kcal/mol. The calculated heat of formation for UF 6 (-511.0 kcal/mol) is within 2.5 kcal/mol of the accurately-known (±0.45 kcal/mol) experimental value.

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

Pyrolysis Pathways of the Furanic Ether 2-Methoxyfuran

Substituted furans, including furanic ethers, derived from nonedible biomass have been proposed as second-generation biofuels. In order to use these molecules as fuels, it is important to understand how they break apart thermally. In this work, a series of experiments were conducted to study the unimolecular and low-pressure bimolecular decomposition mechanisms of the smallest furanic ether, 2-methoxyfuran. Electronic structure (CBS-QB3) calculations indicate this substituted furan has an unusually weak O - CH 3 bond, approximately 190 kJ mol -1 (45 kcal mol -1 ); thus, the primary decomposition pathway is through bond scission resulting in CH 3 and 2-furanyloxy (O - C 4 H 3 O) radicals. Final products from the ring opening of the furanyloxy radical include 2 CO, HC$\equiv$CH, and H. The decomposition of methoxyfuran is studied over a range of concentrations (0.0025-0.1%) in helium or argon in a heated silicon carbide (SiC) microtubular flow reactor (0.66-1 mm i.d., 2.5-3.5 cm long) with reactor wall temperatures from 300 to 1300 K. Inlet pressures to the reactor are 150-1500 Torr, and the gas mixture emerges as a skimmed molecular beam at a pressure of approximately 10 μTorr. Products formed at early pyrolysis times (100 μs) are detected by 118.2 nm (10.487 eV) photoionization mass spectrometry (PIMS), tunable synchrotron VUV PIMS, and matrix infrared absorption spectroscopy. Secondary products resulting from H or CH 3 addition to the parent and reaction with 2-furanyloxy were also observed and include CH 2 =CH-CHO, CH 3 -CH=CH-CHO, CH 3 =CO-CH=CH 2 , and furanones; under the conditions in the reactor, we estimate these reactions contribute to at most 1-3% of total methoxyfuran decomposition. This work also includes observation and characterization of an allylic lactone radical, 2-furanyloxy (O-C 4 H 3 O), with the assignment of several intense vibrational bands in an Ar matrix, an estimate of the ionization threshold, and photoionization efficiency. A pressure-dependent kinetic mechanism is also developed to model the decomposition behavior of methoxyfuran and provide pathways for the minor bimolecular reaction channels that are observed experimentally.

09 BIOMASS FUELS↗
DOE OSTI2025

The Remarkable I 2 O 3 Molecule: A New View from Theory

Atmospheric iodine chemistry has garnered increasing attention as a result of increased iodine emissions. A key subset of this chemistry involves iodine oxides (I 2 O 2–5 ), which serve as precursors to particle formation. Among these, I 2 O 3 is the simplest iodine oxide involved in particle formation, but it has remained undetected in the atmosphere. Previous theoretical studies have characterized this peculiar molecule, primarily using energies to refine geometries obtained at low levels of theory. Due to the reemerging interest in I 2 O 3 , this study presents geometries optimized at the CCSD(T)/aug-cc-pwCVTZ-PP level of theory─marking the first instance, to the best of our knowledge, where this system has been studied exclusively with CCSD(T). Harmonic vibrational frequencies were computed at the same level of theory. Final energetics were obtained using the very high level CCSDT(Q) method with basis sets up to quintuple-zeta cardinality (aug-cc-pwCV5Z-PP) and extrapolated to the CBS limit to yield CCSDT(Q)/CBS//CCSD(T)/aug-cc-pwCVTZ-PP energies. These energies include harmonic zero-point vibrational energy corrections and scalar relativistic energy corrections. Additionally, this study discovers new isomers along the I 2 O 3 potential energy surface, a novel contribution to the field. The performance of different computational methods and DFT functionals commonly used in atmospheric chemistry is also assessed relative to high-level theoretical methods.

basis sets↗
DOE OSTI2025

Ultrafast Correlation Energy Estimator

A virtually no-cost method is proposed that can compute the correlation energies of general, covalently bonded, organic, and inorganic molecules (including conjugated π-electron systems) with a well-defined dominant Lewis structure at the accuracy of 99.5% of the near-exact values determined by the coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] in the complete-basis-set (CBS) limit. This Correlation Energy Per Bond (CEPB) method assigns a partial correlation energy to each bond type (characterized by the identities of the two atoms forming the bond and its integer bond order) and to a lone pair, regardless of the bond length, bond angle, sp-hybridization, π-electron conjugation, ionicity, noncovalent interactions, etc. At its current stage, the method is mainly suitable for near-equilibrium geometries. The correlation energies per bond are determined by a fit to the CCSD(T)/CBS benchmarks. It can neither improve the equilibrium structures nor discern conformers or positional isomers, yet its accuracy for reaction energies rivals that of the second-order Møller–Plesset perturbation theory, which is far more expensive. Its promising performance underscores the possibility that surprisingly compact, chemically intuitive molecular fragments exist into which correlation energies can be partitioned, leading to various ultrafast correlation-energy estimators tailored to different purposes.

Binding energy↗
DOE OSTI2016

Low Temperature Chlorine-Initiated Oxidation of Small-Chain Methyl Esters: Quantification of Chain-Terminating HO 2 -Elimination Channels

In this paper, Cl-initiated oxidation reactions of three small-chain methyl esters, methyl propanoate (CH 3 CH 2 COOCH 3 ; MP), methyl butanoate (CH 3 CH 2 CH 2 COOCH 3 ; MB), and methyl valerate (CH 3 CH 2 CH 2 CH 2 COOCH 3 ; MV), are studied at 1 or 8 Torr and 550 and 650 K. Products are monitored as a function of mass, time, and photoionization energy using multiplexed photoionization mass spectrometry coupled to tunable synchrotron photoionization radiation. Pulsed photolysis of molecular chlorine is the source of Cl radicals, which remove an H atom from the ester, forming a free radical. In each case, after addition of O 2 to the initial radicals, chain-terminating HO 2 -elimination reactions are observed to be important. Branching ratios among competing HO 2 -elimination channels are determined via absolute photoionization spectra of the unsaturated methyl ester coproducts. At 550 K, HO 2 -elimination is observed to be selective, resulting in nearly exclusive production of the conjugated methyl ester coproducts, methyl propenoate, methyl-2-butenoate, and methyl-2-pentenoate, respectively. However, in MV, upon raising the temperature to 650 K, other HO 2 -elimination pathways are observed that yield methyl-3-pentenoate and methyl-4-pentenoate. In each methyl ester oxidation reaction, a peak is observed at a mass consistent with cyclic ether formation, indicating chain-propagating OH loss/ring formation pathways via QOOH intermediates. Evidence is observed for the participation of resonance-stabilized QOOH in the most prominent cyclic ether pathways. Finally, stationary point energies for HO 2 -elimination pathways and select cyclic ether formation channels are calculated at the CBS-QB3 level of theory and assist in the assignment of reaction pathways and final products.

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↗
DOE OSTI2019

Reaction Pathways, Thermodynamics, and Kinetics of Cyclopentanone Oxidation Intermediates: A Theoretical Approach

Despite the promising role of cyclopentanone as a bio-derived fuel, thermodynamic and kinetic data are lacking for low temperature oxidation regimes. In this study, ab initio calculations at the CBS-QB3 level explore the subsequent reactivity that results from O 2 -addition to 2- and 3-oxo cyclopentyl radicals, including expected reaction classes such as intra-H migration, HO 2 -elimination, cyclic ether formation and β-scission along with their thermodynamic parameters. Some of the rates are similar to the analogous reactions of cyclopentane, but some other reactions of cyclopentanone are very different. The carbonyl group hinders H-migration from the α’ position but promotes HO 2 -elimination. Enol peroxy formation from some hydroperoxy alkyl radicals of cyclopentanone is unexpectedly important, and so is HO 2 -elimination by β-scission. Our calculations also indicated that at engine relevant conditions the α-RO 2 prefers to go back to reactants 2-oxo cyclopentyl radical and O 2 . Therefore, the reactions resulting from HO2-addition to 2-oxo cyclopentyl is also provided. The lowest barrier channel identified on the singlet surface corresponds to an unexpected intra OH-migration path concerted with ring opening. In conclusion, this valuable information will advance the construction of improved kinetic models for the oxidation of cylopentanone.

09 BIOMASS FUELS↗
DOE OSTI2023

Barrier Heights for Diels–Alder Transition States Leading to Pentacyclic Adducts: A Benchmark Study of Crowded, Strained Transition States of Large Molecules

Theoretical characterization of reactions of complex molecules depends on providing consistent accuracy for the relative energies of intermediates and transition states. Here we employ the DLPNO-CCSD(T) method with core–valence correlation, large basis sets, and extrapolation to the CBS limit to provide benchmark values for Diels–Alder transition states leading to competitive strained pentacyclic adducts. We then used those benchmarks to test a diverse set of wave function and density functional methods for the absolute and relative barrier heights of these transition states. Our results show that only a few of the tested density functionals can predict the absolute barrier heights satisfactorily, although relative barrier heights are more accurate. The most accurate functionals tested are ωB97M-V, M11plus, ωB97X-V, PBE-D3(0), M11, and MN15 with MUDs from best estimates less than 3.0 kcal. Furthermore, these findings can guide selection of density functionals for future studies of crowded, strained transition states of large molecules.

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

Low temperature (550-700 K) oxidation pathways of cyclic ketones: Dominance of HO 2 -elimination channels yielding conjugated cyclic coproducts

The low-temperature oxidation of three cyclic ketones, cyclopentanone (CPO; C 5 H 8 O), cyclohexanone (CHO; C 6 H 10 O), and 2-methyl-cyclopentanone (2-Me-CPO; CH 3 –C 5 H7 O), is studied between 550 and 700 K and at 4 or 8 Torr total pressure. Initial fuel radicals R are formedvia fast H-abstraction from the ketones by laser-photolytically generated chlorine atoms. Intermediates and products from the subsequent reactions of these radicals in the presence of excess O 2 are probed with time and isomeric resolution using multiplexed photoionization mass spectrometry with tunable synchrotron ionizing radiation. For CPO and CHO the dominant product channel in the R + O 2 reactions is chain-terminating HO 2 -elimination yielding the conjugated cyclic coproducts 2-cyclopentenone and 2-cyclohexenone, respectively. Results on oxidation of 2-Me-CPO also show a dominant contribution from HO 2 -elimination. Moreover, the photoionization spectrum of the co-product suggests formation of 2-methyl-2-cyclopentenone and/or 2-cyclohexenone, resulting from a rapid Dowd–Beckwith rearrangement, preceding addition to O 2 , of the initial (2-oxocyclopentyl)methyl radical to 3-oxocyclohexyl. Cyclic ethers, markers for hydroperoxyalkyl radicals (QOOH), key intermediates in chain-propagating and chain-branching low-temperature combustion pathways, are only minor products. The interpretation of the experimental results is supported by stationary point calculations on the potential energy surfaces of the associated R + O 2 reactions at the CBS-QB3 level. Furthermore, the calculations indicate that HO 2 -elimination channels are energetically favored and product formation via QOOH is disfavored. Lastly, the prominence of chain-terminating pathways linked with HO 2 formation in low-temperature oxidation of cyclic ketones suggests little low-temperature reactivity of these species as fuels in internal combustion engines.

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

Experimental and computational investigation of the bond energy of thorium dicarbonyl cation and theoretical elucidation of its isomerization mechanism to the thermodynamically most stable isomer, thorium oxide ketenylidene cation, OTh + CCO

Collision-induced dissociation (CID) of [Th,2C,2O] + with Xe is performed using a guided ion beam tandem mass spectrometer (GIBMS). The only products observed are ThCO + and Th + by sequential loss of CO ligands. The experimental findings and theoretical calculations support that the structure of [Th,2C,2O] + is the bent homoleptic thorium dicarbonyl cation, Th + (CO) 2 , having quartet spin, which is both thermodynamically and kinetically stable enough in the gas phase to be observed in our GIBMS instrument. Analysis of the kinetic energy-dependent cross sections for this CID reaction yields the first experimental determination of the bond dissociation energy (BDE) of (CO)Th + –CO at 0 K as 1.05 ± 0.09 eV. A theoretical BDE calculated at the CCSD(T) level with cc-pVXZ (X = T and Q) basis sets and a complete basis set (CBS) extrapolation is in very good agreement with the experimental result. Although the doublet spin bent thorium oxide ketenylidene cation, OTh + CCO, is calculated to be the most thermodynamically stable structure, it is not observed in our experiment where [Th,2C,2O] + is formed by association of Th + and CO in a direct current discharge flow tube (DC/FT) ion source. Potential energy profiles of both quartet and doublet spin are constructed to elucidate the isomerization mechanism of Th + (CO) 2 to OTh + CCO. The failure to observe OTh + CCO is attributed to a barrier associated with C–C bond formation, which makes OTh + CCO kinetically inaccessible under our experimental conditions. Furthermore, chemical bonding patterns in low-lying states of linear and bent Th + (CO) 2 and OTh + CCO isomers are also investigated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Theoretical study of the CO 2 –O 2 van der Waals complex: potential energy surface and applications

A four-dimensional-potential energy surface (4D-PES) of the atmospherically relevant carbon dioxide–oxygen molecule (CO 2 –O 2 ) van der Waals complex is mapped using the ab initio explicitly correlated coupled cluster method with single, double, and perturbative triple excitations (UCCSD(T)-F12b), and extrapolation to the complete basis set (CBS) limit using the cc-pVTZ-F12/cc-pVQZ-F12 bases and the l –3 formula. An analytic representation of the 4D-PES was fitted using the method of interpolating moving least squares (IMLS). These calculations predict that the most stable configuration of CO 2 –O 2 complex corresponds to a planar slipped-parallel structure with a binding energy of V ~ –243 cm –1 . Another isomer is found on the PES, corresponding to a non-planar cross-shaped structure, with V ~ –218 cm –1 . The transition structure connecting the two minima is found at V ~ –211 cm –1 . We also performed comparisons with some CO 2 –X van der Waals complexes. Moreover, we provide a SAPT analysis of this molecular system. Then, we discuss the complexation induced shifts of CO 2 and O 2 . Afterwards, this new 4D-PES is employed to compute the second virial coefficient including temperature dependence. A comparison between quantities obtained in our calculations and those from experiments found close agreement attesting to the high quality of the PES and to the importance of considering a full description of the anisotropic potential for the derivation of thermophysical properties of CO 2 –O 2 mixtures.

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

Prediction of the structures and heats of formation of MO 2 , MO 3 , and M 2 O 5 for M = V, Nb, Ta, Pa

Structures for the mono-, di-, and tri-bridge isomers of M 2 O 5 as well as those for the MO 2 and MO 3 fragments for M = V, Nb, Ta, and Pa were optimized at the density functional theory (DFT) level. Single point CCSD(T) calculations extrapolated to the complete basis set (CBS) limit at the DFT geometries were used to predict the energetics. Here, the lowest energy dimer isomer was the di-bridge for M = V and Nb and the tri-bridge for M = Ta and Pa. The di-bridge isomers were predicted to be composed of MO 2 + and MO 3 - fragments, whereas the mono- and tri-bridge are two MO 2 + fragments linked by an O 2 - . The heats of formation of M 2 O 5 dimers, as well as MO 2 and MO 3 neutral and ionic species were predicted using the Feller–Peterson–Dixon (FPD) approach. The heats of formation of the MF 5 species were calculated to provide additional benchmarks. Dimerization energies to form the M 2 O 5 dimers are predicted to become more negative going down group 5 and range from -29 to -45 kcal mol -1 . The ionization energies (IEs) for VO 2 and TaO 2 are essentially the same at 8.75 eV whereas the IEs for NbO 2 and PaO 2 are 8.10 and 6.25 eV, respectively. The predicted adiabatic electron affinities (AEAs) range from 3.75 eV to 4.45 eV for the MO 3 species and vertical detachment energies from 4.21 to 4.59 eV for MO 3 - . The calculated M[double bond, length as m-dash]O bond dissociation energies increase from 143 kcal mol -1 for M = V to ~170 kcal mol -1 for M = Nb and Ta to ~200 kcal mol -1 for M = Pa. The M–O bond dissociation energies are all similar ranging from 97 to 107 kcal mol -1 . Natural bond analysis provided insights into the types of chemical bonds in terms of their ionic character. Pa 2 O 5 is predicted to behave like an actinyl species dominated by the interactions of approximately linear PaO 2 + groups.

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

GaSb doping facilitates conduction band convergence and improves thermoelectric performance in n-type PbS

P-type lead chalcogenides have superior thermoelectric performance because they exhibit the energy convergence of several valence bands. However, despite the existence of two conduction bands, there has been no report about conduction band (CB) convergence for n-type counterparts because of the large energy difference between them. Therefore, new strategies are required to manipulate the CBs if enhancing the electrical transport performance of n-type lead chalcogenides is to be achieved. PbS is a highly attractive member of the lead chalcogenides because of its high earth-abundance and low cost. Here, we report that the introduction of GaSb can successfully dope the PbS matrix with Ga and Sb atoms occupying the Pb site in its rock salt structure. GaSb doping leads to conduction band convergence and enlarged effective density of state mass for n-type PbS. This effect results in superior power factor and decreased lattice thermal conductivity caused by the soft phonon modes and point defect scattering of phonons. Consequently, a record-high average power factor PF avg of ~20.4 μW cm –1 K –2 and figure of merit ZT avg of ~0.84 in the temperature range of 400 K to 923 K were obtained, higher than any n- and p-type PbS-based thermoelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Ground and excited state properties of ThB − and ThB: a theoretical study

In the present work, we studied a series of electronic and spin–orbit states of ThB − and ThB using high-level multireference and coupled-cluster theories. We report the potential energy curves (PECs), equilibrium electron configurations, spectroscopic constants, energetics, and spin–orbit coupling effects of 17 and 19 electronic states of ThB − and ThB, respectively. The ground state of ThB − is a single-reference 1 3 Π with a 1σ 2 2σ 2 3σ 1 1π 3 electron configuration. Detachment of an electron from the doubly occupied 1π orbital of ThB − (1 3 Π) produces the single-reference ground electronic state of ThB (1 4 Σ − ). The ground spin–orbit states of ThB − and of ThB are 1 3 Π + 0 and 1 4 Σ − 3/2 , respectively. The vertical electron detachment energy (VDE) of ThB − and the adiabatic electron attachment energy (AEA) of ThB at our largest CBS-C-CCSD(T)+δT(Q)+δSO (complete basis set effect, spin–orbit effect, and triple and perturbative quadruple electron correlation effect added coupled-cluster theory with single, double, and perturbative triple excitations) level are 1.473 eV and 1.459, respectively. The reaction of Th( 3 F) + B( 2 P o ) produces the ground state of ThB with a bond energy of 2.843 eV. Finally, we estimated a heat of formation, Δ H 0 f (298 K), of 891.01 kJ mol −1 for the ThB molecule. The high-level findings of this work are expected to aid and motivate future experimental spectroscopic investigations of ThB and ThB − species.

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

An accurate potential energy surface for the F + H{sub 2} → HF + H reaction by the coupled-cluster method

A three dimensional potential energy surface for the F + H{sub 2} → HF + H reaction has been computed by the spin unrestricted coupled cluster method with singles, doubles, triples, and perturbative quadruples [UCCSDT(2){sub Q}] using the augmented correlation-consistent polarised valence quadruple zeta basis set for the fluorine atom and the correlation-consistent polarised valence quadruple zeta basis set for the hydrogen atom. All the calculations are based on the restricted open-shell Hartree-Fock orbitals, together with the frozen core approximations, and the UCCSD(T)/complete basis set (CBS) correction term was included. The global potential energy surface was calculated by fitting the sampled ab initio points without any scaling factor for the correlation energy part using a neutral network function method. Extensive dynamics calculations have been carried out on the potential energy surface. The reaction rate constants, integral cross sections, product rotational states distribution, and forward and backward scattering as a function of collision energy of the F + HD → HF + D, F + HD → DF + H, and F + H{sub 2} reaction, were calculated by the time-independent quantum dynamics scattering theory using the new surface. The satisfactory agreement with the reported experimental observations previously demonstrates the accuracy of the new potential energy surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗