A Practical Guide to Reliable First Principles Computational Thermochemistry Predictions Across the Periodic Table

Pyrolysis of bioethanol and biobutanol: A thermodynamic and kinetic study

CONTEXT

Bioethanol and biobutanol are renewable oxygenated fuels derived from biomass, commonly blended with gasoline for use in gasoline engines. These alcohol-based fuels have high oxygen content, promoting more complete combustion and reducing carbon dioxide emissions compared to petroleum fuels. However, during combustion, oxygenated radicals can interact and lead to the formation of formaldehyde, a highly toxic compound. This study delves into the thermodynamic and kinetic study of biofuel pyrolysis using quantum chemical methods. Our results identify C-C bond as the weakest in the initiation step, with bond dissociation enthalpy around 86 kcal/mol. Notably, ethanol exhibits higher bond dissociation energies than butanol. While the initiation step predominantly involves C-C bond breaking, the propagation step reveals a competition between H abstraction and C-C bond cleavage. Analyzing the computed rate constants and Gibbs free energies for radical reactions in the propagation steps indicates the likelihood formation of acetaldehyde, formaldehydes, methane, and ethylene. These products indeed present significant risks to both human health and the environment. This emphasizes the importance of carefully controlling macroscopic thermodynamic variables, such as temperature and pressure, during the pyrolysis of alcohol. Proper regulation of these factors is crucial in minimizing the formation of harmful aldehydes and ensuring a safer and more sustainable process.

METHODS

The reaction mechanisms of thermal decomposition are analyzed using UωB97XD/6-311 + G(3 df,2p), G4MP2, and G4 computational methods. The latter offers highly accurate enthalpies of formation, with a deviation from experiment values approximately 1 kcal/mol, though it is computationally expensive compared to DFT. To evaluate the diradical character of certain open-shell intermediate species, CASSCF and MP2-CASSCF methods, which effectively account for static correlation effects, are employed with the cc-pVDZ basis set. Thermodynamic and kinetic analyses are carried out using both ab initio and semi-empirical approaches through Gaussian 09 and OpenSMOKE + + 0.21.0 programs.

Exploring ground-state and ionization potentials of the H2CO⋯HNO dimers

In this work, we performed an investigation on the structures and ionization potentials (IPs) of the H2CO⋯HNO dimer. The ground-state properties of six conformations were explored using the coupled-cluster with single, double, and perturbative triple excitations, CCSD(T), approach with large correlation consistent basis sets. Conformation III presented the strongest hydrogen-bonding interaction (with the NH⋯O distance being 2.016 Å) and was assigned as the most stable among the conformations. In addition, twelve lowest-lying IPs of all the H2CO⋯HNO conformations were determined using the equation-of-motion ionization potential coupled-cluster with single and double excitations method (EOMIP-CCSD) combined with correlation consistent basis sets, extrapolation to the complete basis set limit, and consideration of core correlation effects. The first IP of conformation III was determined to be 10.46 eV, while the corresponding values for conformations I and II were found as being lower than the value obtained for conformation III by 0.41 eV and 0.24 eV, respectively. These differences (that were also noticed for other low-lying IPs) may be helpful for the assignments of experimental results and, thus, it is expected that the outcomes from this work may serve as motivation for other experimental and theoretical investigations involving H2CO⋯HNO dimers (particularly studies that rely on the quantities obtained here).

Energetic and Electronic Properties of AcX and LaX (X = O and F)

The bonding and spectroscopic properties of LaX and AcX (X = O and F) diatomic molecules were studied by high-level ab initio CCSD(T) and SO-CASPT2 electronic structure calculations. Bond dissociation energies (BDEs) were calculated at the Feller-Peterson-Dixon (FPD) level. Potential energy curves and spectroscopic constants for the lowest-lying spin-orbit Ω states were obtained at the SO-CASPT2/aQ-DK level. A dense manifold of excited states was described for the monofluorides with the ground states well separated from the excited states. The spectroscopic parameters were in good agreement with those reported experimentally for LaO and LaF. For the diatomic molecules containing actinides, no experimental data of these parameters was found, but the results were consistent with other high-level calculations. The BDEs calculated at the FPD level were 791.3 (LaO), 705.2 (AcO), 650.0 (LaF), and 678.6 (AcF) kJ/mol. The NBO analysis showed that the monofluorides are essentially ionic, which explains why the BDE(AcF) is higher than BDE(LaF); for the monoxides, covalent contributions involving the d orbitals of the metal and the p orbitals of the oxygen are stronger for LaO than AcO, which explains the higher BDE for LaO. The bond orders are predicted to be 2 for LaF and AcF, 3 for AcO, and higher than 3 for LaO.

Accurate and reliable thermochemistry by data analysis of complex thermochemical networks using Active Thermochemical Tables: the case of glycine thermochemistry

Active Thermochemical Tables (ATcT) were successfully used to resolve the existing inconsistencies related to the thermochemistry of glycine, based on statistically analyzing and solving a thermochemical network that includes >3350 chemical species interconnected by nearly 35 000 thermochemically-relevant determinations from experiment and high-level theory. The current ATcT results for the 298.15 K enthalpies of formation are -394.70 ± 0.55 kJ mol-1 for gas phase glycine, -528.37 ± 0.20 kJ mol-1 for solid α-glycine, -528.05 ± 0.22 kJ mol-1 for β-glycine, -528.64 ± 0.23 kJ mol-1 for γ-glycine, -514.22 ± 0.20 kJ mol-1 for aqueous undissociated glycine, and -470.09 ± 0.20 kJ mol-1 for fully dissociated aqueous glycine at infinite dilution. In addition, a new set of thermophysical properties of gas phase glycine was obtained from a fully corrected nonrigid rotor anharmonic oscillator (NRRAO) partition function, which includes all conformers. Corresponding sets of thermophysical properties of α-, β-, and γ-glycine are also presented.

Prospects for rank-reduced CCSD(T) in the context of high-accuracy thermochemistry

Obtaining sub-chemical accuracy (1 kJ mol-1) for reaction energies of medium-sized gas-phase molecules is a longstanding challenge in the field of thermochemical modeling. The perturbative triples correction to coupled-cluster single double triple [CCSD(T)] constitutes an important component of all high-accuracy composite model chemistries that obtain this accuracy but can be a roadblock in the calculation of medium to large systems due to its O(N7) scaling, particularly in HEAT-like model chemistries that eschew separation of core and valence correlation. This study extends the work of Lesiuk [J. Chem. Phys. 156, 064103 (2022)] with new approximate methods and assesses the accuracy of five different approximations of (T) in the context of a subset of molecules selected from the W4-17 dataset. It is demonstrated that all of these approximate methods can achieve sub-0.1 kJ mol-1 accuracy with respect to canonical, density-fitted (T) contributions with a modest number of projectors. The approximation labeled Z̃T appears to offer the best trade-off between cost and accuracy and shows significant promise in an order-of-magnitude reduction in the computational cost of the CCSD(T) component of high-accuracy model chemistries.

Anion photoelectron spectroscopy and chemical bonding of ThS2- and ThSO

Anion photoelectron spectra of ThSO- and ThS2- were recorded using the third (355 nm) harmonic of an Nd-YAG laser; these provided the measured vertical detachment energies of each anion. The experiments are supported by extensive coupled cluster calculations on ThSO, ThSO-, ThS2, and ThS2-, as well as the oxygen congeners ThO2 and ThO2-. The ab initio calculations, which included complete basis set extrapolations, spin-orbit effects using four-component coupled cluster, and higher-order correlation contributions through CCSDT(Q), yielded an adiabatic electron affinity for ThO2 that was within 0.02 eV of the previously determined experimental value. The singly occupied molecular orbital (SOMO) in all three anions corresponds primarily to the 7s orbital on Th. Successive substitution of S for each O in ThO2 leads to larger electron affinities and smaller bond angles in the neutral molecules, but larger angles in the anions. As demonstrated by Franck-Condon simulations of the spectra using the CCSD(T) spectroscopic constants, substitution of O by S significantly complicates the resulting detachment spectra due to the lower vibrational frequencies in the sulfur species. Overall the calculated vertical detachment energies are in very good agreement with the experiment. In addition to the adiabatic electron affinities of each species, atomization energies and heats of formation have also been determined via the FPD approach with expected uncertainties of 1-2 kcal/mol.

Composite ab initio vibrational spectroscopy and thermochemistry of low-valency lanthanide compounds: europium dihalides EuX2 (X = F, Cl, Br, I)

A unique diversity in chemical properties of lanthanide compounds making them very useful in various important applications stems not least from the existence of less common Ln oxidation states, in particular +2. However, the molecular properties of low-valency lanthanide compounds are not yet well studied even for the most stable divalent Ln-containing species, such as europium dihalides. In this paper, highly accurate molecular structures, vibrational spectra, and atomisation energies of the europium dihalides EuX2 (X = F, Cl, Br, I) are obtained at the complete basis set CCSD(T) level. All of the EuX2 species are calculated to be non-linear (C2v) exhibiting a regular increase in X-Eu-X bond angle on passing through the halogen series (F → Cl → Br → I), from 117° in EuF2 to 141° in EuI2, which is accompanied by a rapid decrease in the barriers to linearity, h = E(D∞h) - E(C2v), from 2180 cm-1 to 166 cm-1, respectively. The Eu-X bonds in EuX2 appear to be longer but energetically stronger than those in the respective monohalides EuX, whose properties are studied in this work at the same level of theory as EuX2. The performance of the Ln 4f-in-core pseudopotential (PP) approximation for EuX2 is assessed by comparing the Eu PP-based coupled-cluster (CC) calculations with all-electron CC benchmarks. Incorporating the Eu 4f electrons into the PP is shown to cause large errors: up to 10 deg. in X-Eu-X bond angles, 80% in barriers to linearity, and 0.05 Å in the Eu-X bond lengths, which proves the necessity of an explicit treatment of 4f electrons if high accuracy is the goal. The present results set a benchmark in the field of low-valency lanthanide chemistry, in particular for the calibration of DFT functionals.

Factorized Quadruples and a Predictor of Higher-Level Correlation in Thermochemistry

Coupled cluster theory has had a momentous impact on the ab initio prediction of molecular properties, and remains a staple ingratiate in high-accuracy thermochemical model chemistries. However, these methods require inclusion of at least some connected quadruple excitations, which generally scale at best as O ( N 9 ) with the number of basis functions. It is very difficult to predict, a priori, the effect correlation of past CCSD(T) on a given reaction energy. The purpose of this work is to examine cost-effective quadruple corrections based on the factorization theorem of the many-body perturbation theory that may address these challenges. We show that the O ( N 7 ) factorized CCSD(TQf) method introduces minimal error to predicted correlation and reaction energies as compared to the O ( N 9 ) CCSD(TQ). Further, we examine the performance of Goodson's continued fraction method in the estimation of CCSDT(Q)Λ contributions to reaction energies as well as a "new" method related to %TAE[(T)] that we refer to as a scaled perturbation estimator. We find that the scaled perturbation estimator based upon CCSD(TQf)/cc-pVDZ is capable of predicting CCSDT(Q)Λ/cc-pVDZ contributions to reaction energies with an average error of 0.07 kcal mol-1 and an L2D of 0.52 kcal mol-1 when applied to a test-suite of nearly 3000 reactions. This offers a means by which to reliably "ballpark" how important post-CCSD(T) contributions are to reaction energies while incurring no more than CCSD(T) formal cost and a little mental math.

Probing Antioxidant-Related Properties for Phenolic Compounds

In this work, properties related to antioxidant-potential mechanisms (such as the bond dissociation enthalpy, BDE, for the homolytic cleavage of the O-H bond and ionization energies, IEs) were determined for phenol, pyrocatechol, and gallic acid (GA). Both the protonated and deprotonated forms of GA were investigated. The Feller-Peterson-Dixon (FPD) composite method was employed with a variety of computational approaches, i.e., density functional theory, Möller-Plesset perturbation theory, and coupled-cluster-based methods, in combination with large correlation consistent basis sets with extrapolation to the complete basis set limit and consideration of core electron correlation effects. FPD results were compared to experimental and computational data available in the literature, presenting good agreement. For example, the FPD BDE (298 K) obtained for phenol, which was based on valence-correlated MP2/CBS calculations with contributions from correlating all electrons, was determined to be 87.56 kcal/mol, a value that is 0.42 kcal/mol lower than the result obtained in the most recent experiments, 87.98 ± 0.62. Calibration against coupled-cluster calculations was also carried out for phenol. We expect that the outcomes gathered here may help in establishing a general protocol for computational chemists that are interested in determining antioxidant-related properties for phenolic compounds with considerable accuracy as well as to motivate future IE measurements (particularly for GA) to be accomplished in the near future.

W4Λ: Leveraging Λ Coupled-Cluster for Accurate Computational Thermochemistry Approaches

High-accuracy composite wave function methods like Weizmann-4 (W4) theory, high-accuracy extrapolated ab initio thermochemistry (HEAT), and the Feller-Peterson-Dixon (FPD) approach enable sub-kJ/mol accuracy in gas-phase thermochemical properties. Their biggest computational bottleneck is the evaluation of the valence post-CCSD(T) correction term. We demonstrate here, for the W4-17 thermochemistry benchmark and subsets thereof, that the Λ coupled-cluster expansion converges more rapidly and smoothly than the regular coupled-cluster series. By means of CCSDT(Q)Λ and CCSDTQ(5)Λ, we can considerably (up to an order of magnitude) accelerate W4- and W4.3-type calculations without loss in accuracy, leading to the W4Λ and W4.3Λ computational thermochemistry protocols.

Supersonically expanded sodium metal-dilute halogen gas interactions. The importance of reaction populated and energy storing reservoir states and population inversion created amplification in Na2

The reactions of Cl, Br, and I with Nan=2,3 produced in a supersonic expansion form Na2* and Na* excited states extending across the visible and ultraviolet regions. Emission in the region extending from 410 to 600 nm indicates selectively formed excited state Na2 emission features. Experimental evidence suggests that this emission is associated with Na3 + X reactions. Broadband (0.5 cm-1) laser measurements demonstrate gain (population inversion) for select features at∼524-528(1%), ∼492(0.3%), and ∼458.7-461(0.8%) nm. Single mode (0.007 cm-1) measurements extending from 528.03 to 527.63 nm demonstrate amplification involving five to six individual rovibronic levels with a maximum gain close to 3% recorded at 527.9 nm. The observed gain is associated with select transitions from levels of the Na2 11Πu state populated, via identified curve crossings, through collision induced transfer from long-lived Na2 21Σg+ and 11Πg reservoir states. Collision induced population buildup in the lowest vibrational levels of these reservoir states and collision induced transfer to the Na2 11Πu state create a population inversion in transitions to the X 1Σg+ state of Na2. The observed amplification is aided by rapid vibrational and rotational relaxation in both the Na2 ground and excited reservoir states producing amplifiers in the visible region like the HF amplifier in the infrared. This study suggests the importance of reaction populated and energy storing long-lived reservoir states in small sodium molecule combustion processes and indicates the potential for providing new short wavelength visible and ultraviolet amplifiers for future laser-based chemical propulsion concepts.

Prediction of Aqueous Reduction Potentials of X•, ChH•, and XO• Radicals with X = Halogen and Ch = Chalcogen

The aqueous electron affinity and aqueous reduction potentials for F•, Cl•, Br•, I•, OH•, SH•, SeH•, TeH•, ClO•, BrO•, and IO• were calculated using electronic structure methods for explicit cluster models coupled with a self-consistent reaction field (SMD) to treat the aqueous solvent. Calculations were conducted using MP2 and correlated molecular orbital theory up to the CCSD(T)-F12b level for water tetramer clusters and MP2 for octamer cluster. Inclusion of explicit waters was found to be important for accurately predicting the redox potentials in a number of cases. The calculated reduction potentials for X• and ChH• were predicted to within ∼0.1 V of the reported literature values. Fluorine is anomalous due to abstraction of a hydrogen from one of the surrounding water molecules to form a hydroxyl radical and hydrogen fluoride, so its redox potential was calculated using only an implicit model. Larger deviations from experiment were predicted for ClO• and BrO•. These deviations are due to the free energy of solvation of the anion being too negative, as found in the pKa calculations, and that for the neutral being too positive with the current approach.

A larger basis set describes atomization energy core-valence correction better than a higher-order coupled-cluster method

The accuracy of coupled-cluster methods for the computation of core-valence correction to atomization energy was assessed. Truncation levels up to CCSDTQP were considered together with (aug-)cc-pwCVnZ (n = D, T, Q, 5) basis sets and three different extrapolation techniques (canonical and flexible Helgaker formula and Riemann zeta function extrapolation). With the exception of CCSD, a more accurate correction can be obtained from a larger basis set using a lower-level coupled-cluster method, and not vice versa. For the CCSD(T) level, it also implies faster computations with modern codes. We also discussed the importance of moving to higher-order or all-electron methods for geometry optimizations. The present study provides the general knowledge needed for the most accurate state-of-the-art computations.

The Crucial Role of Sb2 F10 in the Chemical Synthesis of F2

The purely chemical synthesis of fluorine is a spectacular reaction which for more than a century had been believed to be impossible. In 1986, it was finally experimentally achieved, but since then this important reaction has not been further studied and its detailed mechanism had been a mystery. The known thermal stability of MnF4 casts serious doubts on the originally proposed hypothesis that MnF4 is thermodynamically unstable and decomposes spontaneously to a lower manganese fluoride and F2 . This apparent discrepancy has now been resolved experimentally and by electronic structure calculations. It is shown that the reductive elimination of F2 requires a large excess of SbF5 and occurs in the last reaction step when in the intermediate [SbF6 ][MnF2 ][Sb2 F11 ] the addition of one more SbF5 molecule to the [SbF6 ]- anion generates a second tridentate [Sb2 F11 ]- anion. The two tridentate [Sb2 F11 ]- anions then provide six fluorine bridges to the Mn atom thereby facilitating the reductive elimination of the two fluorine ligands as F2 .

Coupled Cluster Study of the Heats of Formation of UF6 and the Uranium Oxyhalides, UO2X2 (X = F, Cl, Br, I, and At)

The atomization enthalpies of the U(VI) species UF6 and the uranium oxyhalides UO2X2 (X = F, Cl, Br, I, and 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 multideterminant 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 UF6 (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 the U atom. The analogous value for UO2F2 (607.6 kcal/mol) was in nearly exact agreement with the experiment, but the latter has a stated experimental uncertainty of ±4.3 kcal/mol. The FPD atomization enthalpy for UO2Cl2 (540.4 kcal/mol) was within the experimental error limit of ±5.5 kcal/mol. FPD atomization energies for the non-U-containing molecules (used for reaction enthalpies) H2O and HX (X = F, Cl, Br, I, and 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 UF6 (-511.0 kcal/mol) is within 2.5 kcal/mol of the accurately known (±0.45 kcal/mol) experimental value.

Ab Initio Rovibrational Spectroscopy of the Acetylide Anion

In this work the rovibrational spectrum of the acetylide anion HCC- is investigated using high-level electronic structure methods and variational rovibrational calculations. Using a composite approach the potential energy surface and dipole surface is constructed from explicitly correlated coupled-cluster accounting for corrections due to core-valence correlation, scalar relativistic effects and higher-order excitation effects. Previous approaches for approximating the latter are critically evaluated. Employing the composite potential, accurate spectroscopic parameters determined from variational calculations are presented. In comparison to the few available reference data the present results show excellent agreement with ground state rotational constants within 0.005% of the experimental value. Intensities determined from the variational calculations suggest the bending fundamental transition ν2 around 510 cm-1 to be the best target for detection. The rather weak CD stretching fundamental ν1 in deuterated isotopologues show a second-order resonance with the (0,20,1) state and the consequences are discussed in some detail. The spectroscopic parameters and band intensities provided for a number of vibrational bands in isotopologues of the acetylide anion should facilitate future spectroscopic investigations.

Influence of the complete basis set approximation, tight weighted-core, and diffuse functions on the DLPNO-CCSD(T1) atomization energies of neutral H,C,O-compounds

The impact of complete basis set extrapolation schemes (CBS), diffuse functions, and tight weighted-core functions on enthalpies of formation predicted via the DLPNO-CCSD(T1) reduced Feller-Peterson-Dixon approach has been examined for neutral H,C,O-compounds. All tested three-point (TZ/QZ/5Z) extrapolation schemes result in mean unsigned deviation (MUD) below 2 kJ mol^-1 relative to the experiment. The two-point QZ/5Z and TZ/QZ CBS 1 / l max 3 extrapolation schemes are inferior to their inverse power counterpart ( 1 / l max + 1 / 2 4 ) by 1.3 and 4.3 kJ mol^-1 . The CBS extrapolated frozen core atomization energies are insensitive (within 1 kJ mol^-1 ) to augmentation of the basis set with tight weighted core functions. The core-valence correlation effects converge already at triple-ζ, although double-ζ/triple-ζ CBS extrapolation performs better and is recommended. The effect of diffuse function augmentation converges slowly, and cannot be reproduced with double- ζ or triple- ζ calculations as these are plagued with basis set superposition and incompleteness errors.

Active Thermochemical Tables: Enthalpies of Formation of Bromo- and Iodo-Methanes, Ethenes and Ethynes

The thermochemistry of halocarbon species containing iodine and bromine is examined through an extensive interplay between new Feller-Peterson-Dixon (FPD) style composite methods and a detailed analysis of all available experimental and theoretical determinations using the thermochemical network that underlies the Active Thermochemical Tables (ATcT). From the computational viewpoint, a slower convergence of the components of composite thermochemistry methods is observed relative to species that solely contain first row elements, leading to a higher computational expense for achieving comparable levels of accuracy. Potential systematic sources of computational uncertainty are investigated, and, not surprisingly, spin-orbit coupling is found to be a critical component, particularly for iodine containing molecular species. The ATcT analysis of available experimental and theoretical determinations indicates that prior theoretical determinations have significantly larger uncertainties than originally reported, particularly in cases where molecular spin-orbit effects were ignored. Accurate and reliable heats of formation are reported for 38 halogen containing systems, based on combining the current computations with previous experimental and theoretical work via the ATcT approach.

A high level theory investigation on the lowest-lying ionization potentials of glycine (NH2CH2COOH)

In this work, an investigation on the ionization potentials (IPs) of the glycine molecule (NH₂CH₂COOH) is presented. IPs ranging up to ∼20 eV were probed for each of the six conformations considered, with the referred threshold being chosen based on both: (i) the observations by recent photoelectron-photoion coincidence (PEPICO) experiments and (ii) the energy range of relevance to the modeling of other photo-induced processes (e.g., photoionization). For computing the IPs, the equation-of-motion ionization potential coupled-cluster with single and double excitations method (EOMIP-CCSD) was employed with large correlation consistent aug-cc-pVXZ and aug-cc-pCVXZ (X = D, T, and Q) basis sets. Extrapolation to the complete basis set limit and consideration of core electron correlation effects were also taken into account. Subsequently, the Feller-Peterson-Dixon (FPD) approach was used for considering all the contributions and to obtain accurate IPs. In addition, coupled-cluster with single and double excitations as well as perturbative triples, CCSD(T), was also used with the aug-cc-pVTZ basis set. When compared to each other, results obtained through the use of these approaches yielded excellent agreement. In general, the outcomes from the present work provide additional information to the insights gathered from the recent PEPICO experiments as well as accurate IPs for all 6 conformations of glycine using an approach based on high levels of theory. Hence, it is expected that other investigations focusing on photo-induced processes originating from NH₂CH₂COOH (for instance, the computational modeling of its photoionization) will be motivated for study in the future.

Ab initio composite strategies and multireference approaches for lanthanide sulfides and selenides

The f-block ab initio correlation consistent composite approach was used to predict the dissociation energies of lanthanide sulfides and selenides. Geometry optimizations were carried out using density functional theory and coupled cluster singles, doubles, and perturbative triples with one- and two-component Hamiltonians. For the two-component calculations, relativistic effects were accounted for by utilizing a third-order Douglas-Kroll-Hess Hamiltonian. Spin-orbit coupling was addressed with the Breit-Pauli Hamiltonian within a multireference configuration interaction approach. The state averaged complete active space self-consistent field wavefunctions obtained for the spin-orbit coupling energies were used to assign the ground states of diatomics, and several diagnostics were used to ascertain the multireference character of the molecules.

MP2-F12 Basis Set Convergence near the Complete Basis Set Limit: Are h Functions Sufficient?

We have investigated the title question for the W4-08 thermochemical benchmark using l-saturated truncations of a large reference (REF) basis set, as well as for standard F12-optimized basis sets. With the REF basis set, the root-mean-square (RMS) contribution of i functions to the MP2-F12 total atomization energies (TAEs) is about 0.01 kcal/mol, the largest individual contributions being 0.04 kcal/mol for P₂ and P₄. However, even for these cases, basis set extrapolation from {g,h} basis sets adequately addresses the problem. Using basis sets insufficiently saturated in the spdfgh angular momenta may lead to exaggerated i function contributions. For extrapolation from spdfg and spdfgh basis sets, basis set convergence appears to be quite close to the theoretical asymptotic ∝ L^–7 behavior. We hence conclude that h functions are sufficient even for highly demanding F12 applications. With one-parameter extrapolation, spdf and spdfg basis sets are adequate, aug-cc-pV{T,Q}Z-F12 yielding a RMSD = 0.03 kcal/mol. A limited exploration of CCSD(F12*) and CCSD-F12b suggests our conclusions are applicable to higher-level F12 methods as well.

Ab Initio Composite Approaches for Heavy Element Energetics: Ionization Potentials for the Actinide Series of Elements

The first, second, and third gas-phase ionization potentials have been determined for the actinide series of elements using an ab initio composite scalar and fully relativistic approach, employing the coupled cluster with single, double, and perturbative triple excitations (CCSD(T)) and Dirac Hartree-Fock (DHF) methods, extrapolated to the complete basis set (CBS) limit. The impact of electron correlation and basis set choice within this framework are examined. Additionally, the first three ionization potentials were obtained using an ab initio heavy element correlation-consistent Composite Approach (here referred to as α-ccCA). This is the first utilization of a ccCA for actinide species.

Tackling an accurate description of molecular reactivity with double-hybrid density functionals

In this Communication, we assess a panel of 18 double-hybrid density functionals for the modeling of the thermochemical and kinetic properties of an extended dataset of 449 organic chemistry reactions belonging to the BH9 database. We show that most of DHs provide a statistically robust performance to model barrier height and reaction energies in reaching the “chemical accuracy.” In particular, we show that nonempirical DHs, such as PBE0-DH and PBE-QIDH, or minimally parameterized alternatives, such as ωB2PLYP and B2K-PLYP, succeed to accurately model both properties in a balanced fashion. We demonstrate, however, that parameterized approaches, such as ωB97X-2 or DSD-like DHs, are more biased to only one of both properties.

Electronically Excited Complex Formation in Magnesium Cluster-Halogen Atom Reactions

A near ultraviolet transition of Mg₂F has been observed in emission from the reaction between magnesium clusters, most likely Mg₃, and fluorine atoms. Because there is little evidence for upper-state internal excitation, the spectrum is assigned assuming that the upper state is quenched to its lowest vibrational levels. Two of possibly three ground-state vibrational frequencies, υ₁ = 516 ± 10 cm^-1 and υ₂ = 104 ± 10 cm^-1, have been established. Dispersed laser-induced fluorescence studies extrapolating on the observed chemiluminescence indicate an excited-state symmetric stretch frequency of order 370 ± 30 cm^-1. Electronic structure calculations at the CCSD(T)/CBS level predict that the ground state of Mg₂F has C_2v symmetry and can be described as an Mg₂⁺F^- ion pair with two Mg-F bonds. Like the MgF A-X transition that is largely a transition between Mg orbitals, the observed transition in Mg₂F is largely between orbitals on the magnesium dimer ion. The asymmetric C_∞v Mg₂⁺F^- complex is also a minimum and is predicted to be 6.7 kcal/mol higher in energy. Calculated structures for the Mg₂Cl isomers are also presented and used to further interpret the experimental results for the reaction of Mg clusters with Cl atoms. In contrast to Mg₂F, the ground state of Mg₂Cl is a linear C_∞v MgMgCl structure with the C_2v and D_∞h isomers of the MgClMg structure slightly higher in energy.

Predicting the Mechanism and Products of CO2 Capture by Amines in the Presence of H2O

An extensive correlated molecular orbital theory study of the reactions of CO₂ with a range of substituted amines and H₂O in the gas phase and aqueous solution was performed at the G3(MP2) level with a self-consistent reaction field approach. The G3(MP2) calculations were benchmarked at the CCSD(T)/CBS level for NH₃ reactions. A catalytic NH₃ reduces the energy barrier more than a catalytic H₂O for the formation of H₂NCOOH and H₂CO₃. In aqueous solution, the barriers to form both H₂NCOOH and H₂CO₃ are reduced, with HCO₃^- formation possible with one amine present and H₂NCOO^- formation possible only with two amines. Further reactions of H₂NCOOH to form HNCO and urea via the Bazarov reaction have high barriers and are unlikely in both the gas phase and aqueous solution. Reaction coordinates for CH₃NH₂, CH₃CH₂NH₂, (CH₃)₂NH, CH₃CH₂CH₂NH₂, (CH₃)₃N, and DMAP were also calculated. The barrier for proton transfer correlates with amine basicity for alkylammonium carbamate (ΔG^‡_aq < 15 kcal/mol) and alkylammonium bicarbonate (ΔG^‡_aq < 30 kcal/mol) formation. In aqueous solution, carbamic acids, carbamates, and bicarbonates can all form in small amounts with ammonium carbamates dominating for primary and secondary alkylamines. These results have implications for CO₂ capture by amines in both the gas phase and aqueous solution as well as in the solid state, if enough water is present.

Ionization potentials for the H2CO trimer

In this work, a computational study on the ionization potentials (IPs) of the formaldehyde trimer, (H₂CO)₃, is presented. Twelve lowest-lying vertical IPs were determined through the use of the coupled-cluster level of theory using correlation consistent basis sets with extrapolation to the complete basis set limit and consideration of core electron correlation effects. Specifically, the equation-of-motion ionization potential coupled-cluster with single and double excitations method with the aug-cc-pVnZ and aug-cc-pCVnZ (n = D and T) basis sets was used. The Feller-Peterson-Dixon (FPD) composite approach was employed to provide accurate IPs, and eight conformations of (H₂CO)₃ were considered. The FPD IPs determined for (H₂CO)₃ were found to be systematically lower than those computed for the dimer and monomer of H₂CO in the pattern IP_(monomer) > IP_(dimer) > IP_(trimer) for a given IP. In addition, the IPs calculated when considering only the more stable conformation (C₀) are in good agreement with those obtained using the eight conformations of the H₂CO trimer, and thus, the actual conformation played only a minor role in determining such properties in the present case. By providing first accurate IP results for the H₂CO trimer, we hope to motivate future experimental and computational investigations (e.g., studies involving photoionization) that rely on such quantities.

The electron affinity of the uranium atom

The results of a combined experimental and computational study of the uranium atom are presented with the aim of determining its electron affinity. Experimentally, the electron affinity of uranium was measured via negative ion photoelectron spectroscopy of the uranium atomic anion, U^-. Computationally, the electron affinities of both thorium and uranium were calculated by conducting relativistic coupled-cluster and multi-reference configuration interaction calculations. The experimentally determined value of the electron affinity of the uranium atom was determined to be 0.309 ± 0.025 eV. The computationally predicted electron affinity of uranium based on composite coupled cluster calculations and full four-component spin-orbit coupling was found to be 0.232 eV. Predominately due to a better convergence of the coupled cluster sequence for Th and Th^-, the final calculated electron affinity of Th, 0.565 eV, was in much better agreement with the accurate experimental value of 0.608 eV. In both cases, the ground state of the anion corresponds to electron attachment to the 6d orbital.

Coupled Cluster Studies of Platinum-Actinide Interactions. Thermochemistry of PtAnOn+ (n = 0-2 and An = U, Np, Pu)

Accurate Pt-An bond dissociation enthalpies (BDEs) for PtAnOⁿ⁺ (An = U, Np, Pu and n = 0-2) and the corresponding enthalpies for the Pt + OAnOⁿ⁺ substitution reactions have been studied for the first time using an accurate composite coupled cluster approach. Analogous O-AnOⁿ⁺ bond dissociation enthalpies are also presented. To make the study possible, new correlation consistent basis sets optimized using the all-electron third-order Douglas-Kroll-Hess (DKH3) scalar relativistic Hamiltonian are developed and reported for Pt and Au, with accompanying benchmark calculations of their atomic ionization potentials to demonstrate the effectiveness of the new basis sets. For the charged PtAnOⁿ⁺ species (n = 1, 2), a low-spin state (LSS) for which the Pt-An σ bond is doubly occupied is studied together with a high-spin state (HSS) obtained by unpairing the σ bond orbital and placing one electron into the An 5f shell. The relative energies of the two spin states have been compared and qualitatively assessed via natural population and natural bond analyses. The enthalpies for the Pt substitution reactions, i.e., Pt + OAnOⁿ⁺ → PtAnOⁿ⁺ + O, are calculated to range from about 14-62 kcal/mol, and the Pt-AnOⁿ⁺ bond dissociation enthalpies range from about 78-149 kcal/mol for the ground electronic states. For the PtAnO⁺ species, the LSSs were all predicted to be the ground state, whereas the PtAnO²⁺ molecules all favored the HSSs. The prediction for PtUO²⁺ is consistent with previous theoretical findings. The natural bond orbital analyses indicate a triple bond between An and O, with a double to quadruple bond between the An and Pt.

Beyond Chemical Accuracy for Alkane Thermochemistry: The DHthermo Approach

The so-called protobranching phenomenon, that is the greater stability of branched alkanes with respect to their linear isomers, represents an interesting challenge for approaches based on density functional theory (DFT), since it requires a balanced description of several electronic effects, including (intramolecular) dispersion forces. Here, we investigate this problem using a protocol recently developed based on double-hybrid functionals and a small basis set, DH-SVPD, suited for noncovalent interactions. The energies of bond separation reactions (BSR), defined on the basis of an isodesmic principle, are taken as reference properties for the evaluation of 15 DFT approaches. The obtained results show that error lower than the so-called "chemical accuracy" (<1.0 kcal/mol) can be obtained by the proposed protocol on both relative reaction energies and enthalpies. These results are then verified on the standard BSR36 data set and support the proposition of our computational protocol, named DHthermo, where any DH functional, such as PBE-QIDH or B2PLYP, provides accurate results when coupled to an empirical dispersion correction and the DH-SVPD basis set. This protocol not only gives subchemical accuracy on the thermochemistry of alkanes but it is extremely easy to use with common quantum-chemistry codes.

Looking for the bricks of the life in the interstellar medium: The fascinating world of astrochemistry

The discovery in the interstellar medium of molecules showing a certain degree of complexity, and in particular those with a prebiotic character, has attracted great interest. A complex chemistry takes place in space, but the processes that lead to the production of molecular species are a matter of intense discussion, the knowledge still being at a rather primitive stage. Debate on the origins of interstellar molecules has been further stimulated by the identification of biomolecular building blocks, such as nucleobases and amino acids, in meteorites and comets. Since many of the molecules found in space play a role in the chemistry of life, the issue of their molecular genesis and evolution might be related to the profound question of the origin of life itself. Understanding the underlying chemical processes, including the production, reactions and destruction of compounds, requires the concomitant study of spectroscopy, gas-phase reactivity, and heterogeneous processes on dust-grains. The aim of this contribution is to provide a general view of a complex and multifaceted challenge, while focusing on the role played by molecular spectroscopy and quantum-chemical computations. In particular, the derivation of the molecular spectroscopic features and the investigation of gas-phase formation routes of prebiotic species in the interstellar medium are addressed from a computational point of view.

Gas-Phase Thermochemistry of MX3 and M2X6 (M = Sc, Y; X = F, Cl, Br, I) from a Composite Reaction-Based Approach: Homolytic versus Heterolytic Cleavage

A domain-based local-pair natural-orbital coupled-cluster approach with single, double, and improved linear-scaling perturbative triple correction via an iterative algorithm, DLPNO-CCSD(T₁), was applied within the framework of the Feller-Peterson-Dixon approach to derive gas-phase heats of formation of scandium and yttrium trihalides and their dimers via a set of homolytic and heterolytic dissociation reactions. All predicted heats of formation moderately depend on the reaction type with the most and least negative values obtained for homolytic and heterolytic dissociation, respectively. The basis set size dependence, as well as the influence of static correlation effects not covered by the standard (DLPNO-)CCSD(T) approach, suggests that exploitation of the heterolytic dissociation reactions with the formation of M³⁺ and X^- ions leads to the most robust heats of formation. The gas-phase formation enthalpies ΔH_f°(0 K)/ΔH_f°(298.15 K) and absolute entropies S°(298.15 K) were obtained for the first time for the Sc₂F₆, Sc₂Br₆, and Sc₂I₆ species. For ScBr₃, ScI₃, Sc₂Cl₆, and Y₂Cl₆, we suggest a reexamination of the experimental heats of formation available in the literature. For other compounds, the predicted values were found to be in good agreement with the experimental estimates. Extracted MX₃ (M = Sc, Y; X = F, Cl, Br, and I) 0 K atomization enthalpies indicate weaker bonding when moving from fluorine to iodine and from yttrium to scandium. Likewise, the stability of yttrium trihalide dimers degrades when going from fluorine to iodine. Respective scandium trihalide dimers are less stable, with 0 K dimer dissociation energy decreasing in the row fluorine - chlorine - bromine ≈ iodine. Correlation of the (n - 1)s²p⁶ electrons on bromine and iodine, inclusion of zero-point energy, relativistic effects, and the effective-core-potential correction as well as amelioration of the DLPNO localization inaccuracy are shown to be of similar magnitude, which is critical if accurate heats of formation are a goal.

Structural investigations, quantum mechanical studies on proton and metal affinity and biological activity predictions of selpercatinib

Cancer of the lungs and thyroid is particularly difficult to manage and treat. Notably, selpercatinib has recently been suggested as an effective drug to combat these diseases. The entire world is currently tackling the pandemic caused by the SARS-CoV-19 virus. Numerous pharmaceuticals have been evaluated for the management of the disease caused by SARS-CoV-19 (i.e., COVID-19). In this study, selpercatinib was proposed as a potential inhibitor of different SARS-CoV-19 proteins. Several intriguing effects of the molecule were found during the conducted computational investigations. Selpercatinib could effectively act as a proton sponge and exhibited high proton affinity in solution. Moreover, it was able to form complexes with metal ions in aqueous solutions. Specifically, the compound displayed high affinity towards zinc ions, which are important for the prevention of virus multiplication inside human cells. However, due to their charge, zinc ions are not able to pass the lipid bilayer and enter the cell. Thus, it was determined that selpercatinib could act as an ionophore, effectively transporting active zinc ions into cells. Furthermore, various quantum mechanical analyses, including energy studies, evaluation of the reactivity parameters, examination of the electron localisation and delocalisation properties, as well as assessment of the nonlinear optical (NLO) properties and information entropy, were conducted herein. The performed docking studies (docking scores -9.3169, -9.1002, -8.1853 and -8.1222 kcal mol^-1) demonstrated that selpercatinib strongly bound with four isolated SARS-CoV-2 proteins.

Effective Molecular Descriptors for Chemical Accuracy at DFT Cost: Fragmentation, Error-Cancellation, and Machine Learning

Recent advances in theoretical thermochemistry have allowed the study of small organic and bio-organic molecules with high accuracy. However, applications to larger molecules are still impeded by the steep scaling problem of highly accurate quantum mechanical (QM) methods, forcing the use of approximate, more cost-effective methods at a greatly reduced accuracy. One of the most successful strategies to mitigate this error is the use of systematic error-cancellation schemes, in which highly accurate QM calculations can be performed on small portions of the molecule to construct corrections to an approximate method. Herein, we build on ideas from fragmentation and error-cancellation to introduce a new family of molecular descriptors for machine learning modeled after the Connectivity-Based Hierarchy (CBH) of generalized isodesmic reaction schemes. The best performing descriptor ML(CBH-2) is constructed from fragments preserving only the immediate connectivity of all heavy (non-H) atoms of a molecule along with overlapping regions of fragments in accordance with the inclusion-exclusion principle. Our proposed approach offers a simple, chemically intuitive grouping of atoms, tuned with an optimal amount of error-cancellation, and outperforms previous structure-based descriptors using a much smaller input vector length. For a wide variety of density functionals, DFT+ΔML(CBH-2) models, trained on a set of small- to medium-sized organic HCNOSCl-containing molecules, achieved an out-of-sample MAE within 0.5 kcal/mol and 2σ (95%) confidence interval of <1.5 kcal/mol compared to accurate G4 reference values at DFT cost.

Canonical and DLPNO-Based G4(MP2)XK-Inspired Composite Wave Function Methods Parametrized against Large and Chemically Diverse Training Sets: Are They More Accurate and/or Robust than Double-Hybrid DFT?

The large and chemically diverse GMTKN55 benchmark was used as a training set for parametrizing composite wave function thermochemistry protocols akin to G4(MP2)XK theory (Chan, B.; Karton, A.; Raghavachari, K. J. Chem. Theory Comput. 2019, 15, 4478–4484). On account of their availability for elements H through Rn, Karlsruhe def2 basis sets were employed. Even after reparametrization, the GMTKN55 WTMAD2 (weighted mean absolute deviation, type 2) for G4(MP2)-XK is actually inferior to that of the best rung-4 DFT functional, ωB97M-V. By increasing the basis set for the MP2 part to def2-QZVPPD, we were able to substantially improve performance at modest cost (if an RI-MP2 approximation is made), with WTMAD2 for this G4(MP2)-XK-D method now comparable to the better rung-5 functionals (albeit at greater cost). A three-tier approach with a scaled MP3/def2-TZVPP intermediate step, however, leads to a G4(MP3)-D method that is markedly superior to even the best double hybrids ωB97M(2) and revDSD-PBEP86-D4. Evaluating the CCSD(T) component with a triple-ζ, rather than split-valence, basis set yields only a modest further improvement that is incommensurate with the drastic increase in computational cost. G4(MP3)-D and G4(MP2)-XK-D have about 40% better WTMAD2, at similar or lower computational cost, than their counterparts G4 and G4(MP2), respectively: detailed comparison reveals that the difference lies in larger molecules due to basis set incompleteness error. An E2/{T,Q} extrapolation and a CCSD(T)/def2-TZVP step provided the G4-T method of high accuracy and with just three fitted parameters. Using KS orbitals in MP2 leads to the G4(MP3|KS)-D method, which entirely eliminates the CCSD(T) step and has no steps costlier than scaled MP3; this shows a path forward to further improvements in double-hybrid density functional methods. None of our final selections require an empirical HLC correction; this cuts the number of empirical parameters in half and avoids discontinuities on potential energy surfaces. G4-T-DLPNO, a variant in which post-MP2 corrections are evaluated at the DLPNO-CCSD(T) level, achieves nearly the accuracy of G4-T but is applicable to much larger systems.

Probing the ionization potentials of the formaldehyde dimer

In this work, we present a computational investigation on the ionization potentials (IPs) of the formaldehyde dimer, (H₂CO)₂. Twelve lowest lying IPs (corresponding to the entire valence orbitals) for both C_2h and C_s symmetry conformers have been computed at the coupled cluster level of theory using large correlation consistent basis sets with extrapolation to the complete basis set limit and consideration of core electron correlation effects. Specifically, the equation-of-motion ionization potential coupled-cluster with single and double (EOMIP-CCSD) excitations method with the aug-cc-pVXZ and aug-cc-pCVXZ (X = T, Q, and 5) basis sets combined with the Feller-Peterson-Dixon approach was employed, as well as CCSD with perturbative triples [CCSD(T)] with the aug-cc-pVTZ basis sets. In general, excellent agreement was observed from the comparison between the results obtained through the use of these approaches. In addition, the IPs for the formaldehyde monomer were also obtained using such methodologies and the results compared with existing experimental data; excellent agreement was also observed in this case. To the best of our knowledge, this work represents the first of its kind to determine the IPs for all these systems using a high level theory approach and is presented to motivate experimental investigations, e.g., studies involving photoionization, particularly for the formaldehyde dimer. The equilibrium binding energy of the C_2h dimer is calculated in this work at the CCSD(T)/aug-cc-pVTZ level of theory to be -4.71 kcal/mol. At this same level of theory, the equilibrium isomerization energy between C_2h and C_s conformers is 0.76 kcal/mol (C_s conformer being more stable).

Spectroscopic and theoretical studies of UN and UN. J Chem Phys 2020;152:094302. [PMID: 33480743 DOI: 10.1063/1.5144299]

The low-energy electronic states of UN and UN⁺ have been examined using high-level electronic structure calculations and two-color photoionization techniques. The experimental measurements provided an accurate ionization energy for UN (IE = 50 802 ± 5 cm^-1). Spectra for UN⁺ yielded ro-vibrational constants and established that the ground state has the electronic angular momentum projection Ω = 4. Ab initio calculations were carried out using the spin-orbit state interacting approach with the complete active space second-order perturbation theory method. A series of correlation consistent basis sets were used in conjunction with small-core relativistic pseudopotentials on U to extrapolate to the complete basis set limits. The results for UN correctly obtained an Ω = 3.5 ground state and demonstrated a high density of configurationally related excited states with closely similar ro-vibrational constants. Similar results were obtained for UN⁺, with reduced complexity owing to the smaller number of outer-shell electrons. The calculated IE for UN was in excellent agreement with the measured value. Improved values for the dissociation energies of UN and UN⁺, as well as their heats of formation, were obtained using the Feller-Peterson-Dixon composite thermochemistry method, including corrections up through coupled cluster singles, doubles, triples and quadruples. An analysis of the ab initio results from the perspective of the ligand field theory shows that the patterns of electronic states for both UN and UN⁺ can be understood in terms of the underlying energy level structure of the atomic metal ion.

The challenging playground of astrochemistry: an integrated rotational spectroscopy - quantum chemistry strategy

While it is now well demonstrated that the interstellar medium (ISM) is characterized by a diverse and complex chemistry, a significant number of features in radioastronomical spectra are still unassigned and call for new laboratory efforts, which are increasingly based on integrated experimental and computational strategies. In parallel, the identification of an increasing number of molecules containing more than five atoms and at least one carbon atom (the so-called "interstellar" complex organic molecules), which can play a relevant role in the chemistry of life, raises the additional issue of how these species can be produced in the typical harsh conditions of the ISM. On these grounds, this perspective aims to present an integrated rotational spectroscopy - quantum chemistry approach for supporting radioastronomical observations and a computational strategy for contributing to the elucidation of chemical reactivity in the interstellar space.

Domain-based local pair natural orbital methods within the correlation consistent composite approach

Ab initio composite approaches have been utilized to model and predict main group thermochemistry within 1 kcal mol^-1 , on average, from well-established reliable experiments, primarily for molecules with less than 30 atoms. For molecules of increasing size and complexity, such as biomolecular complexes, composite methodologies have been limited in their application. Therefore, the domain-based local pair natural orbital (DLPNO) methods have been implemented within the correlation consistent composite approach (ccCA) framework, namely DLPNO-ccCA, to reduce the computational cost (disk space, CPU (central processing unit) time, memory) and predict energetic properties such as enthalpies of formation, noncovalent interactions, and conformation energies for organic biomolecular complexes including one of the largest molecules examined via composite strategies, within 1 kcal mol^-1 , after calibration with 119 molecules and a set of linear alkanes. © 2019 Wiley Periodicals, Inc.

Electronic Structure Predictions of the Energetic Properties of Tellurium Fluorides

The heats of formation, bond dissociation energies (BDEs), fluoride affinities (FA), fluorocation affinities (FCA), electron affinities (EA), and ionization energies (IP) of TeF _n ( n = 1-6) have been predicted using the Feller-Peterson-Dixon (FPD) approach. To benchmark the approach, the bond dissociation energies of Te₂ and TeO, the heats of formation of Te₂, TeH₂, TeO, and TeO₂, and the electron affinity for TeO and TeO₂ were calculated as there are experimental thermodynamic data available for these tellurium compounds, which allow confirmation of the heat of formation of Te gas as Δ H_f,0K(Te) = 50.7 ± 0.6 kcal/mol. Spin-orbit corrections are required for good results and cannot be ignored. A comparison among fluoride affinities, fluorocation affinities, electron affinities, and ionization energies of TeF _n and SeF _n is reported.