Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Quantum chemistry methods”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Theoretical Modeling of Electrochemical Proton-Coupled Electron Transfer

Proton-coupled electron transfer (PCET) plays an essential role in a wide range of electrocatalytic processes. Here, a vast array of theoretical and computational methods have been developed to study electrochemical PCET. These methods can be used to calculate redox potentials and pKa values for molecular electrocatalysts, proton-coupled redox potentials and bond dissociation free energies for PCET at metal and semiconductor interfaces, and reorganization energies associated with electrochemical PCET. Periodic density functional theory can also be used to compute PCET activation energies and perform molecular dynamics simulations of electrochemical interfaces. Various approaches for maintaining a constant electrode potential in electronic structure calculations and modeling complex interactions in the electric double layer (EDL) have been developed. Theoretical formulations for both homogeneous and heterogeneous electrochemical PCET spanning the adiabatic, nonadiabatic, and solvent-controlled regimes have been developed and provide analytical expressions for the rate constants and current densities as functions of applied potential. The quantum mechanical treatment of the proton and inclusion of excited vibronic states have been shown to be critical for describing experimental data, such as Tafel slopes and potential dependent kinetic isotope effects. The calculated rate constants can be used as input to microkinetic models and voltammogram simulations to elucidate complex electrocatalytic processes.

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

Local self-interaction correction method with a simple scaling factor

A recently proposed local self-interaction correction (LSIC) method [Zope et al., J. Chem. Phys., 2019, 151, 214108] when applied to the simplest local density approximation provides a significant improvement over standard Perdew–Zunger SIC (PZSIC) for both equilibrium properties such as total or atomization energies as well as properties involving stretched bond such as barrier heights. The method uses an iso-orbital indicator to identify the single-electron regions. To demonstrate the LSIC method, Zope et al. used the ratio z σ of von Weizsäcker τ W σ and total kinetic energy densities τ σ , (z σ = τ W σ /τ σ ) as a scaling factor to scale the self-interaction correction. The present work further explores the LSIC method using a ratio of orbital and spin densities as a simpler scaling factor in place of the ratio of kinetic energy densities. We compute a wide array of both, equilibrium and non-equilibrium properties using LSIC and orbital scaling methods using this simple scaling factor and compare them with previously reported results. Our study shows that LSIC with the simple scaling factor performs better than PZSIC, with results comparable to those obtained by LSIC( z σ ) for most properties, but has slightly larger errors than LSIC( z σ ). Furthermore, we study the binding energies of small water clusters using both scaling factors. Our results show that LSIC with zσ has limitations in predicting the cluster binding energies of weakly bonded systems due to the inability of z σ to distinguish weakly bonded regions from slowly varying density regions. LSIC when used with the density ratio as a scaling factor, on the other hand, provides a good description of water cluster binding energies, thus highlighting the appropriate choice of the iso-orbital indicator.

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

Iron(III)–Oxo Cluster Chemistry with Dimethylarsinate Ligands: Structures, Magnetic Properties, and Computational Studies

A program has been initiated to develop Fe III /oxo cluster chemistry with the ‘pseudo-carboxylate’ ligand dimethylarsinate (Me 2 AsO 2 - ) for comparison with the well investigated Fe III /oxo/carboxylate cluster area. The synthesis and characterization of three polynuclear Fe III complexes are reported, [Fe 12 O 4 (O 2 C t Bu) 8 (O 2 AsMe 2 ) 17 (H 2 O) 3 ]Cl 3 (1), Na 2 [Fe 12 Na 2 O 4 (O 2 AsMe 2 ) 20 (NO 3 ) 6 (Me 2 AsO 2 H) 2 (H 2 O) 4 ](NO 3 ) 6 (2) and [Fe 3 (O 2 AsMe 2 ) 6 (Me 2 AsO 2 H) 2 (hqn) 2 ](NO 3 ) (3), where hqn is 8-hydroxyquinoline. The Fe 12 core of 1 is a type never previously encountered in Fe III carboxylate chemistry, consisting of two Fe 6 units each of which comprises two {Fe 3 (μ 3 -O 2- )} units bridged by three Me 2 AsO 2 - groups and linked into an Fe 12 loop structure by two anti-anti η 1 :η 1 :μ Me 2 AsO 2 - groups, a bridging mode extremely rare with carboxylates. 2 also consists of two Fe 6 units, differing in their ligation from those in 1, and this time linked together into a linear structure by a central {Na 2 (NO 3 ) 2 } bridging unit. 3 is a linear Fe 3 complex with no monoatomic bridges between Fe III ions, a very rare situation in Fe III chemistry with any ligands, and unprecedented in Fe carboxylate chemistry. The distinct differences observed in arsinate vs carboxylate ligation modes are rationalized largely based on the greater basicity of the former vs the latter. Variable-temperature dc and ac magnetic susceptibility data reveal all Fe 2 pairwise interactions to be antiferromagnetic. For 1 and 2, the different J ij couplings were estimated by use of a magnetostructural correlation for high nuclearity Fe III -oxo clusters and by density functional theory calculations using broken symmetry methods, allowing identification of their relative spin vector alignments and thus rationalization of their S = 0 ground states. The J ij values were then used as input values to give excellent fits of the experimental χM T vs T data. For 3, the fits of the experimental χM T vs T data to the Van Vleck equation or with PHI gave a very weak J 12 = -0.8(1) cm -1 (H = –2JŠ i ·Š j convention) between adjacent Fe III ions, and an S = 5/2 ground state. Furthermore, these initial Fe III arsinate complexes also provide structural parameters that help validate literature assignments of arsinate binding modes to iron oxide/hydroxide minerals as part of environmental concerns of using arsenic-containing herbicides in agriculture.

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

Automated discovery of a robust interatomic potential for aluminum

Machine learning, trained on quantum mechanics (QM) calculations, is a powerful tool for modeling potential energy surfaces. A critical factor is the quality and diversity of the training dataset. Here we present a highly automated approach to dataset construction and demonstrate the method by building a potential for elemental aluminum (ANI-Al). In our active learning scheme, the ML potential under development is used to drive non-equilibrium molecular dynamics simulations with time-varying applied temperatures. Whenever a configuration is reached for which the ML uncertainty is large, new QM data is collected. The ML model is periodically retrained on all available QM data. The final ANI-Al potential makes very accurate predictions of radial distribution function in melt, liquid-solid coexistence curve, and crystal properties such as defect energies and barriers. We perform a 1.3M atom shock simulation and show that ANI-Al force predictions shine in their agreement with new reference DFT calculations.

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

Non-Solvent Induced Phase Separation Enables Designer Redox Flow Battery Electrodes

Porous carbonaceous electrodes are performance-defining components in redox flow batteries (RFBs), where their properties impact the efficiency, cost, and durability of the system. The overarching challenge is to simultaneously fulfill multiple seemingly contradictory requirements—i.e., high surface area, low pressure drop, and facile mass transport—without sacrificing scalability or manufacturability. Here, non-solvent induced phase separation (NIPS) is proposed as a versatile method to synthesize tunable porous structures suitable for use as RFB electrodes. The variation of the relative concentration of scaffold-forming polyacrylonitrile to pore-forming poly(vinylpyrrolidone) is demonstrated to result in electrodes with distinct microstructure and porosity. Tomographic microscopy, porosimetry, and spectroscopy are used to characterize the 3D structure and surface chemistry. Flow cell studies with two common redox species (i.e., all-vanadium and Fe2+/3+) reveal that the novel electrodes can outperform traditional carbon fiber electrodes. It is posited that the bimodal porous structure, with interconnected large (>50 µm) macrovoids in the through-plane direction and smaller (<5 µm) pores throughout, provides a favorable balance between offsetting traits. Although nascent, the NIPS synthesis approach has the potential to serve as a technology platform for the development of porous electrodes specifically designed to enable electrochemical flow technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH) 4 H]- model employing state-of-the-art ab initio methods

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ΔE estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree–Fock included (0%–54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ΔE is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn–Sham orbitals gives results within 4–11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ~5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

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

Polariton Control of Molecular Charge Transfer in Perylene Diimide Semiconductors

We report the modulation of molecular charge transfer in a bay-substituted perylene diimide derivative embedded in a planar distributed Bragg reflector microcavity. Angle-resolved reflectance spectra confirm the formation of upper and lower polaritons with clear Rabi splitting, indicating strong coupling between the cavity mode and molecular excitons. Using broadband transient absorption spectroscopy, we compare the ultrafast dynamics of cavity and non-cavity films. While excited-state absorption and stimulated emission pathways remain largely unchanged, kinetic modeling reveals a moderate increase in the charge transfer rate and yield under strong coupling. This enhancement is attributed to a reduction in the effective driving force via the formation of the lower polariton, placing the system deeper into the Marcus inverted regime. Our results demonstrate a promising non-chemical method leveraging cavity quantum electrodynamics to modulate charge separation processes in molecular semiconductors.

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

Understanding the Electronic Behavior of f-Element Intermetallics

Due to their partially filled f-orbitals, lanthanide and actinide intermetallic materials exhibit intriguing electronic, structural, and magnetic properties. These materials provide a unique platform to explore exotic magnetic and quantum phenomena that arise from their strongly correlated f-electrons, offering insights into f-electron behavior. Incorporating light elements such as boron and carbon can significantly impact the electronic, chemical, and physical properties of these materials, with varying influence between lanthanides and actinides due to the more localized f-electrons in lanthanides. In this work, lanthanide-based intermetallic single crystals containing light atoms were synthesized using the molten metal flux growth method. In this method, one or more low-melting metals such as aluminum, gallium, tin, and bismuth are used in excess as the reaction medium. This solution-state approach allows for lower reaction temperatures compared to solid-state reactions and enables the isolation of kinetic products rather than thermodynamically stable compounds. The synthesized crystals were analyzed using an extensive array of analytical and computational characterization techniques, including single crystal X-ray diffraction, scanning electron microscopy, energy-dispersive spectroscopy, magnetization measurements, and density functional theory calculations. These analyses provided detailed insights into the structural and magnetic properties of the synthesized materials.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2023

Butterfly, Vinylidene-Like, Monobridged and Trans Structures of Si 2 H 2 + : Comparison to the Well-Characterized Neutral Si 2 H 2

This paper investigates four of the constitutional isomers of Si 2 H 2 + , namely the butterfly, vinylidene-like, monobridged, and trans structures. These isomer geometries were all studied using the CCSD(T) method with basis sets as large as cc-pV5Z. Higher level methods CCSDT and CCSDT(Q) were used for final energetics. It is found that the butterfly isomer has the lowest energy, followed by vinylidene-like and monobridged structures; the trans structure has the highest energy of the four. All structures were compared with their neutral counterparts. Partial charges are computed and it is found that all isomers share the positive charge between the silicon atoms. NBO analysis shows that the cation Si-Si bond order is reduced by 0.5 relative to the neutral versions of the same isomer in the case of the butterfly, vinylidene-like, and monobridged; and reduced by 0.6 for the trans isomer. Vibrational frequencies, infrared intensities, and dipole moments are predicted to encourage the spectroscopic identification of Si 2 H 2 + .

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

Towards complete assignment of the infrared spectrum of the protonated water cluster H + (H 2 O) 21

The spectroscopic features of protonated water species in dilute acid solutions have been long sought after for understanding the microscopic behavior of the proton in water with gas-phase water clusters H + (H 2 O) n extensively studied as bottom-up model systems. We present a new protocol for the calculation of the infrared (IR) spectra of complex systems, which combines the fragment-based Coupled Cluster method and anharmonic vibrational quasi-degenerate perturbation theory, and demonstrate its accuracy towards the complete and accurate assignment of the IR spectrum of the H + (H 2 O) 21 cluster. The site-specific IR spectral signatures reveal two distinct structures for the internal and surface four-coordinated water molecules, which are ice-like and liquid-like, respectively. The effect of inter-molecular interaction between water molecules is addressed, and the vibrational resonance is found between the O-H stretching fundamental and the bending overtone of the nearest neighboring water molecule. The revelation of the spectral signature of the excess proton offers deeper insight into the nature of charge accommodation in the extended hydrogen-bonding network underpinning this aqueous cluster.

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

Colloidal Synthesis Path to 2D Crystalline Quantum Dot Superlattices

By combining colloidal nanocrystal synthesis, self-assembly, and solution phase epitaxial growth techniques, we developed a general method for preparing single dot thick atomically attached quantum dot (QD) superlattices with high-quality translational and crystallographic orientational order along with state-of-the-art uniformity in the attachment thickness. The procedure begins with colloidal synthesis of hexagonal prism shaped core/shell QDs (e.g., CdSe/CdS), followed by liquid subphase self-assembly and immobilization of superlattices on a substrate. Solution phase epitaxial growth of additional semiconductor material fills in the voids between the particles, resulting in a QD-in-matrix structure. The photoluminescence emission spectra of the QD-in-matrix structure retains characteristic 0D electronic confinement. Importantly, annealing of the resulting structures removes inhomogeneities in the QD–QD inorganic bridges, which our atomistic electronic structure calculations demonstrate would otherwise lead to Anderson-type localization. The piecewise nature of this procedure allows one to independently tune the size and material of the QD core, shell, QD–QD distance, and the matrix material. These four choices can be tuned to control many properties (degree of quantum confinement, quantum coupling, band alignments, etc.) depending on the specific applications. Lastly, cation exchange reactions can be performed on the final QD-in-matrix, as demonstrated herein with a CdSe/CdS to HgSe/HgS conversion.

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

Evaluating fast methods for static polarizabilities on extended conjugated oligomers

Given the importance of accurate polarizability calculations to many chemical applications, coupled with the need for efficiency when calculating the properties of sets of molecules or large oligomers, we present a benchmark study examining possible calculation methods for polarizable materials. We first investigate the accuracy of the additive model used in GFN2, a highly-efficient semi-empirical tight-binding method, and the D4 dispersion model, comparing its predicted additive polarizabilities to ωB97XD results for a subset of PubChemQC and a compiled benchmark set of molecules spanning polarizabilities from approximately 3 Å 3 to 600 Å 3 , with some compounds in the range of approximately 1200–1400 Å 3 . Although we find additive GFN2 polarizabilities, and thus D4, to have large errors with polarizability calculations on large conjugated oligomers, it would appear an empirical quadratic correction can largely remedy this. We also compare the accuracy of DFT polarizability calculations run using basis sets of varying size and level of augmentation, determining that a non-augmented basis set may be used for large, highly polarizable species in conjunction with a linear correction factor to achieve accuracy extremely close to that of aug-cc-pVTZ.

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

Characterization of mApple as a Red Fluorescent Protein for Cryogenic Single-Molecule Imaging with Turn-Off and Turn-On Active Control Mechanisms

Single-molecule superresolution microscopy is a powerful tool for the study of biological structures on size scales smaller than the optical diffraction limit. Imaging samples at cryogenic temperatures (77 K) reduces the quantum yield of photobleaching for many fluorescent labels, yielding localization precisions below 10 nm. Cryogenic imaging further enables correlation with cryogenic electron tomography. A key limitation in applying methods such as PALM and STORM to samples maintained at 77 K is the limited number of fluorophores known to undergo efficient turn-on and turn-off mechanisms necessary to control the sparsity of active emitters. Here we find that mApple, a red-emitting fluorescent protein, undergoes a novel turn-off mechanism in response to simultaneous illumination with two colors of light. This turn-off mechanism enables localization of many individual molecules in initially bright samples, but the final density of localizable emitters is limited by relatively inefficient turn-on (photoactivation). Bulk excitation and emission spectroscopy shows that mApple has access to two distinct emissive states as well as dark states accessible optically or through changes in pH. The bright and stable emission of mApple enables widefield collection of single-molecule emission spectra, which highlight the complex nature and environmental sensitivity of states observed in red fluorescent proteins.

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

Plasma and trap-based techniques for science with antimatter

Positrons (i.e., antielectrons) find use in a wide variety of applications, and antiprotons are required for the formation and study of antihydrogen. Available sources of these antiparticles are relatively weak. To optimize their use, most applications require that the antiparticles be accumulated into carefully prepared plasmas. We present an overview of the techniques that have been developed to efficiently accumulate low energy antiparticles and create, in particular, tailored antiparticle plasmas. Techniques are also described to create tailored antiparticle beams. Many of these techniques are based on methods first developed by the nonneutral plasma community using electron plasmas for increased data rate. They have enabled the creation and trapping of antihydrogen, have been critical to studies of positron and positronium interactions with matter, including advanced techniques to characterize materials and material surfaces, and have led to the creation and study of the positronium molecule. Rather than attempting to be comprehensive, we focus on techniques that have proven most useful, applications where there has been significant, recent progress, and areas that hold promise for future advances. Examples of the latter include the ever more precise comparisons of the properties of antihydrogen and hydrogen, tests of gravity using antihydrogen and positronium atoms, and efforts to create and study phases of the many-electron, many-positron system.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Static Subspace Approximation for Random Phase Approximation Correlation Energies: Applications to Materials for Catalysis and Electrochemistry

Modeling complex materials using high-fidelity, ab initio methods at low cost is a fundamental goal for quantum chemical software packages. The GW approximation and random phase approximation (RPA) provide a unified description of both electronic structure and total energies using the same physics in a many-body perturbative approach that can be more accurate than generalized-gradient density functional theory (DFT) methods. However, GW/RPA implementations have historically been limited to either specific materials classes or application toward small chemical systems. Here, the static subspace approximation allows for reduced cost full-frequency GW/RPA calculations and has previously been benchmarked thoroughly for GW calculations. Here, we describe our approach to including partial occupations of electronic orbitals in full-frequency GW and RPA calculations for the study of electrocatalysts. We benchmarked RPA total energy calculations using the subspace approximation across a diverse test suite of materials for a variety of computational parameters. The benchmarking quantifies the impact of different extrapolation procedures for representing the static polarizability at infinite screened cutoff, and shows that using screened cutoffs above 20-25 Ryd result in diminishing accuracy returns for predicting RPA total energies. Additionally, for moderately sized electrocatalytic models, 2-3 times fewer computational resources are used to compute RPA total energies by representing the static polarizability with 20-30% of the static subspace basis, with an error of approximately 0.01 eV or better in RPA adsorption energy calculations. Finally, we show that for these electrochemical models RPA can shift DFT adsorption energy shifts by up to 0.5 eV and that GW can frequently shift DFT eigenvalues of surface and adsorbate states by approximately 0.5-1 eV.

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

Resolving femtosecond photoinduced energy flow: capture of nonadiabatic reaction pathway topography and wavepacket dynamics from photoexcitation through the conical intersection seam (Final Technical Report)

The dynamics that take place within just tens to hundreds of femtoseconds following the absorption of light by a molecule can play a critical role in how the absorbed energy is directed, allowing it to be used for a specific function or dissipated harmlessly. The form of chemical change that occurs rapidly in these molecules is called a “nonadiabatic electronic transition.” Such transitions are known to mediate energy flow in natural biological systems such as the ultraviolet photoprotection mechanism of DNA and the first step of the human vision response. Understanding how these mechanisms work precisely may help scientists achieve controlled manipulation of solar energy or optical control of a wide range of energy management functions in artificial systems. Experimental methods, however, have not yet allowed a precisely resolved and complete measurement of nonadiabatic electronic transitions. This constitutes a major obstacle to progress in the field. For progress to occur that would inform a wide body of research aiming to efficiently harness the energy of light for practical purposes, it is especially important to benchmark computational models of the molecules undergoing these rapid changes with experimental measurements, in order to learn which models are accurate. With Dept. of Energy funding, we have made strong progress towards establishing a new optical method for experimentally detecting the full nonadiabatic electronic transition. This requires having coordinated pulses of light covering the visible through the mid-infrared range of the electromagnetic spectrum that last only ten femtoseconds. We have developed a new, relatively simple approach for generating such pulses of laser light, and have incorporated them into a time-resolved spectrometer for measuring rapid changes in molecules. These tools can provide the greater precision and new types of data that are needed to benchmark computational models of molecular change and thus to make progress in the field. Our tools were tested on graphene, an excellent solid-state sample for verifying the capabilities and limitations of our instrumentation. The investment made in these tools by the Dept. of Energy Office of Science will allow new fundamental scientific understanding of energy dynamics in molecules in future studies.

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

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

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

36 MATERIALS SCIENCE↗
DOE OSTI2025

The Electronic Structure of Zirconium and Hafnium Monochalcogenides

High-level ab initio CCSD(T) and spin–orbit icMRCI+Q calculations were used to predict potential energy curves (PECs) for the lowest-lying states of ZrO, ZrS, HfO, and HfS. The prediction of the ground state is basis set dependent at the icMRCI+Q level for ZrO and ZrS due to the small singlet–triplet splitting between the lowest 1 Σ + and 3 Δ states. CCSD(T) with a spin orbit correction predicted the 1 Σ + ground state in agreement with experiment. New all-electron basis sets were developed for Hf to improve the results over those predicted by use of effective core potentials (ECPs) that subsume the 4f electrons into the definition of the core. The use of the new DK-4f basis sets rather than ECPs became more important for HfO and HfS where there is a lack of a good core–valence separation. icMRCI+Q, CCSD(T), and DFT calculations for the spectroscopic parameters of ZrO, ZrS, HfO, and HfS were benchmarked with available experimental data. Bond dissociation energies (BDEs) of these four systems were calculated at the Feller–Peterson–Dixon (FPD) level to be 762.1 (ZrO), 543.5 (ZrS), 803.8 (HfO), and 575.1 kJ/mol (HfS), in excellent agreement with experiment. The HfS BDE was remeasured using the R3PI method, providing an updated experimental measurement of D 0 (HfS) = 5.978 ± 0.002 eV = 576.8 ± 0.2 kJ/mol. This experimental value, combined with experimental measurements of the ionization energies of Hf and HfS, gives the cationic BDE of D 0 (Hf + -S) = 5.124 ± 0.002 eV = 494.4 ± 0.2 kJ/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗