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

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.

Gas Phase Silver Thermochemistry from First Principles

Domain-based local pair natural orbital coupled cluster approach with single, double, and perturbative triple excitations, DLPNO-CCSD(T), has been applied within a framework of a reduced version of the reaction-based Feller-Peterson-Dixon (FPD) scheme to predict gas phase heats of formation and absolute entropies of silver inorganic and organometallic compounds. First, we evaluated all existing experimental data currently limited by thermodynamic functions of 10 silver substances (AgH, AgF, AgBr, AgI, Ag₂, Ag₂S, Ag₂Se, Ag₂Te, AgCN, AgPO₂). The mean average deviation between computed and experimental heats of formation was found to be 1.9 kcal/mol. Notably, all predicted heats of formation turned out to be within the error bounds of their experimental counterparts. Second, we predicted heats of formation and entropies for additional 90 silver species with no experimental data available, substantially enriching silver thermochemistry. Combination of gas phase heats of formation Δ H_f and entropies S° of AgNO₂, AgSCN, Ag₂SO₄, and Ag₂SeO₄ obtained in this work, with respective solid-state information, resulted in accurate sublimation thermochemistry of these compounds. Complementation of predicted Δ H_f with heats of formation of some neutrals and positive ions produced 33 silver bond strengths of high reliability. Obtained thermochemical data are promising for developing the concepts of silver chemistry. In addition, derived heats of formation and bond dissociation enthalpies, due to their high diversity, are found to be relevant for testing and training of computational chemistry methods.

The interplay and strength of the π⋯HF, C⋯HF, F⋯HF and F⋯HC hydrogen bonds upon the formation of multimolecular complexes based on C2H2⋯HF and C2H4⋯HF small dimers

The conception of this theoretical research was idealized aiming to unveil the intermolecular structures of complexes formed by acetylene or ethylene and hydrofluoric acid. At light of computational calculations by using the B3LYP/6-311++G(d,p) method, the geometries of the C₂H₂⋯(HF), C₂H₂⋯2(HF), C₂H₂⋯4(HF), C₂H₄⋯(HF), C₂H₄⋯2(HF) and C₂H₄⋯4(HF) hydrogen-bonded complexes were fully optimized. Moreover, the Post-Hartree-Fock calculations MP2/6-311++G(d,p), MP2/aug-cc-pVTZ, MP4(SDQ)/6-311++G(d,p) and CCSD/6-311++G(d,p) also were also used. The infrared spectra were analyzed in order to identify the new vibrational modes and frequencies of the proton donors shifted to red region. Through the modeling of charge-fluxes on the basis of the Quantum Theory of Atoms In Molecules (QTAIM) and, by contradicting the expectation of the hydrofluorination mechanisms of acetylene or ethylene, C⋯HF was recognized as a new type of hydrogen bond instead of the already well known π⋯H. The calculations of the Natural Bonding Orbital (NBO) and Charges derived from the Electrostatic Potential Grid-based (ChElPG) were also applied to interpret the shifting frequencies as well as measuring of the punctual charge-transfer after the formation of the complexes. Finally, the determination of the stabilization energy was carried out through the arguments of the Fock matrix in NBO basis and through the supermolecule approach. Also it is worthwhile to notice that some algebraic formulations were used for determining the electronic cooperative effect (CE).

Benson group additivity values of phosphines and phosphine oxides: Fast and accurate computational thermochemistry of organophosphorus species

Composite quantum chemical methods W1X-1 and CBS-QB3 are used to calculate the gas phase standard enthalpy of formation, entropy, and heat capacity of 38 phosphines and phosphine oxides for which reliable experimental thermochemical information is limited or simply nonexistent. For alkyl phosphines and phosphine oxides, the W1X-1, and CBS-QB3 results are mutually consistent and in excellent agreement with available G3X values and empirical data. In the case of aryl-substituted species, different computational methods show more variation, with G3X enthalpies being furthest from experimental values. The calculated thermochemical data are subsequently used to determine Benson group additivity contributions for 24 Benson groups and group pairs involving phosphorus, thereby allowing fast and accurate estimations of thermochemical data of many organophosphorus compounds of any complexity. Such data are indispensable, for example, in chemical process design or estimating potential hazards of new chemical compounds. © 2018 Wiley Periodicals, Inc.