A pseudopotential-based composite method: The relativistic pseudopotential correlation consistent composite approach for molecules containing 4d transition metals (Y–Cd)

Unraveling the electronic structure of LuH, LuN, and LuNH: building blocks of new materials

Advances in superconductor technology have been pursued for decades, moving towards room temperature models, such as a postulated nitrogen-doped lutetium hydride network. While experimental observations have been contradictory, insight into the building blocks of potential new superconductor materials can be gained theoretically, unravelling the fascinating electronic structure of these compounds at a molecular level. Here, the fundamental building blocks of lutetium materials (LuH, LuN, and LuNH) have been examined. The structures, spectroscopic constants for the ground and excited states, and the potential energy curves have been obtained for these species using complete active self-consistent field (CASSCF) and multireference configuration interaction with Davidson's correction (MRCI+Q) methods. For LuNH, the energetic properties of its isomers are determined. The bond dissociation energies of the three building blocks are calculated with the state-of-the-art f-block ab initio correlation consistent composite approach (f-ccCA) and the high accuracy extrapolated ab initio thermochemistry (HEAT) scheme. As well, an analysis of different formation pathways of LuNH has been provided.

Critical benchmarking of popular composite thermochemistry models and density functional approximations on a probabilistically pruned benchmark dataset of formation enthalpies

First-principles calculation of the standard formation enthalpy, ΔH_f° (298 K), in such a large scale as required by chemical space explorations, is amenable only with density functional approximations (DFAs) and certain composite wave function theories (cWFTs). Unfortunately, the accuracies of popular range-separated hybrid, "rung-4" DFAs, and cWFTs that offer the best accuracy-vs-cost trade-off have until now been established only for datasets predominantly comprising small molecules; their transferability to larger systems remains vague. In this study, we present an extended benchmark dataset of ΔH_f° for structurally and electronically diverse molecules. We apply quartile-ranking based on boundary-corrected kernel density estimation to filter outliers and arrive at probabilistically pruned enthalpies of 1694 compounds (PPE1694). For this dataset, we rank the prediction accuracies of G4, G4(MP2), ccCA, CBS-QB3, and 23 popular DFAs using conventional and probabilistic error metrics. We discuss systematic prediction errors and highlight the role an empirical higher-level correction plays in the G4(MP2) model. Furthermore, we comment on uncertainties associated with the reference empirical data for atoms and the systematic errors stemming from these that grow with the molecular size. We believe that these findings will aid in identifying meaningful application domains for quantum thermochemical methods.

Thermochemistry of Tungsten-3p Elements for Density Functional Theory, Caveat Lector!

There are two primary foci in this research on WE (E = Si, P, and S) bonds: prediction of their bond dissociation enthalpies (BDEs), including σ- and π-bond energy components, and assessing the uncertainty of these BDE predictions for levels of theory commonly used in the literature. The internal standards for computational accuracy include metal-element bond lengths (mean absolute error = 1.8 ± 1.2%), main group homolog BDEs versus higher levels of ab initio theory (W1U and G4 BDEs, R² = 0.98), and DLPNO-CCSD(T)/def2-QZVPP calculations for metal-ligand BDEs (R² = 0.88). The W═Si first π-bond is underreported for density functional theory (DFT)/MP2 methods versus DLPNO-CCSD(T), while the latter shows negligible strength for the W;Si second π-bond, consistent with the literature. This research highlights clear issues with the underlying assumptions required for the use of perturbation theory methods for the fragments derived from W-P homolysis. The difficulties associated with modeling the metal thermochemistry with DFT (and MP2) levels of theory are manifest in the broad standard deviations observed. However, the average BDEs found using 48 popular DFT and MP2 levels of theory are reliable, 10.8 ± 6.8% mean absolute error (with W-P removed) versus DLPNO-CCSD(T), with the caveat that the individual basis set/pseudopotential/valence basis set combination can vary wildly. Analysis of the absolute error percentages with respect to the level of theory indicates little benefit to going higher on Jacob's Ladder, as simpler methods have lower error versus high-level ab initio techniques such as G4 and DLPNO-CCSD(T).

Predicting Bond Dissociation Energies and Bond Lengths of Coordinatively Unsaturated Vanadium-Ligand Bonds

Understanding the electronic structure of coordinatively unsaturated transition-metal compounds and predicting their physical properties are of great importance for catalyst design. Bond dissociation energy D_e and bond length r_e are two of the fundamental quantities for which good predictions are important for a successful design strategy. In the present work, recent experimentally measured bond energies and bond lengths of VX diatomic molecules (X = C, N, S) are used as a gauge to consider the utility of a number of electronic structure methods. Single-reference methods are one focus because of their efficiency and utility in practical calculations, and multireference configuration interaction (MRCISD) methods and a composite coupled cluster (CCC) method are a second focus because of their potential high accuracy. The comparison is especially challenging because of the large multireference M diagnostics of these molecules, in the range 0.15-0.19. For the single-reference methods, Kohn-Sham density functional theory (KS-DFT) has been tested with a variety of approximate exchange-correlation functionals. Of these, MOHLYP provides the bond dissociation energies in best agreement with experiments, and BLYP provides the bond lengths that are in best agreement with experiments; but by requiring good performance for both the D_e and r_e of the vanadium compounds, MOHLYP, MN12-L, MGGA_MS1, MGGA_MS0, O3LYP, and M06-L are the most highly recommended functionals. The CCC calculations include up to connected pentuple excitations for the valence electrons and up to connected quadruple excitations for the core-valence terms; this results in highly accurate dissociation energies and good bond lengths. Averaged over the three molecules, the mean unsigned deviation of CCC bond energies from experimental ones is only 0.4 kcal/mol, demonstrating excellent convergence of theory and experiments.

The structure of ScC2 (X̃2A1): A combined Fourier transform microwave/millimeter-wave spectroscopic and computational study

Pure rotational spectra of Sc¹³C₂ (X̃²A₁) and Sc¹²C¹³C (X̃²A') have been measured using Fourier transform microwave/millimeter-wave methods. These molecules were synthesized in a DC discharge from the reaction of scandium vapor, produced via laser ablation, with ¹³CH₄ or ¹³CH₄/¹²CH₄, diluted in argon. The N_Ka,Kc = 1_0,1 → 0_0,0, 2_0,2 → 1_0,1, 3_0,3 → 2_0,2, and 4_0,4 → 3_0,3 transitions in the frequency range of 14 GHz-61 GHz were observed for both species, each exhibiting hyperfine splittings due to the nuclear spins of ¹³C (I = 1/2) and/or Sc (I = 7/2). These data have been analyzed with an asymmetric top Hamiltonian, and rotational, spin-rotation, and hyperfine parameters have been determined for Sc¹³C₂ and Sc¹²C¹³C. In addition, a quartic force field was calculated for ScC₂ and its isotopologues using a highly accurate coupled cluster-based composite method, incorporating complete basis set extrapolation, scalar relativistic corrections, outer core and inner core electron correlation, and higher-order valence correlation effects. The agreement between experimental and computed rotational constants, including the effective constant (B + C), is ∼0.5% for all three isotopologues. This remarkable agreement suggests promise in predicting rotational spectra of new transition metal-carbon bearing molecules. In combination with previous work on Sc¹²C₂, an accurate structure for ScC₂ has been established using combined experimental (B, C) and theoretical (A) rotational constants. The radical is cyclic (or T-shaped) with r_(Sc-C) = 2.048(2) Å, r_(C-C) = 1.272(2) Å, and ∠_(C-Sc-C) = 36.2(1)°. The experimental and theoretical results also suggest that ScC₂ contains a C₂ ^- moiety and is largely ionic.

SAMPL6 logP challenge: machine learning and quantum mechanical approaches

Two different types of approaches: (a) approaches that combine quantitative structure activity relationships, quantum mechanical electronic structure methods, and machine-learning and, (b) electronic structure vertical solvation approaches, were used to predict the logP coefficients of 11 molecules as part of the SAMPL6 logP blind prediction challenge. Using electronic structures optimized with density functional theory (DFT), several molecular descriptors were calculated for each molecule, including van der Waals areas and volumes, HOMO/LUMO energies, dipole moments, polarizabilities, and electrophilic and nucleophilic superdelocalizabilities. A multilinear regression model and a partial least squares model were used to train a set of 97 molecules. As well, descriptors were generated using the molecular operating environment and used to create additional machine learning models. Electronic structure vertical solvation approaches considered include DFT and the domain-based local pair natural orbital methods combined with the solvated variant of the correlation consistent composite approach.

Prediction of pKa s of Late Transition-Metal Hydrides via a QM/QM Approach

Three implicit solvation models, the conductor-like polarizable continuum model (C-PCM), the conductor-like screening model (COSMO), and universal implicit solvent model (SMD), combined with a hybrid two layer QM/QM approach (ONIOM), were utilized to calculate the pK_a values, using a direct thermodynamic scheme, of a set of Group 10 transition metal (TM) hydrides in acetonitrile. To obtain the optimal combination of quantum methods for ONIOM calculations with implicit solvation models, the influence of factors, such as the choice of density functional and basis set, the atomic radii used to build a cavity in the solvent, and the size of the model system in an ONIOM scheme, was examined. Additionally, the impact of Grimme's empirical dispersion correction and exact exchange was also investigated. The results were calibrated by experimental data. This investigation provides insight about effective models for the prediction of thermodynamic properties of TM-containing complexes with bulky ligands. © 2019 Wiley Periodicals, Inc.

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.

Prediction of CB[8] host-guest binding free energies in SAMPL6 using the double-decoupling method

This study reports the results of binding free energy calculations for CB[8] host-guest systems in the SAMPL6 blind challenge (receipt ID 3z83m). Force-field parameters were developed specific for each of host and guest molecules to improve configurational sampling. We used quantum mechanical (QM) implicit solvent calculations and QM force matching to determine non-bonded (partial atomic charges) and bonded terms, respectively. Free energy calculations were carried out using the double-decoupling method (DDM) combined with Hamiltonian replica exchange method (HREM) and Bennett acceptance ratio (BAR) method. The root mean square error (RMSE) of the predicted values using DDM with respect to the experimental results was 4.32 kcal/mol. The coefficient of determination (R²) and Kendall rank coefficient (τ) were 0.49 and 0.52, respectively, highest of all submissions. In addition, these were compared to the results obtained by umbrella sampling (US) and weighted histogram analysis method (WHAM). Overall, DDM achieved a higher prediction accuracy than the US method. Results are discussed in terms of parameterization and free energy simulations.

Towards a quantum chemical protocol for the prediction of rovibrational spectroscopic data for transition metal molecules: Exploration of CuCN, CuOH, and CuCCH

High accuracy electronic structure computations for small transition metal-containing molecules have been a long term challenge. Due to coupling between electronic and nuclear wave functions, even experimental/theoretical identification of the ground electronic state requires tremendous efforts. Quartic force fields (QFFs) are effective ab initio tools for obtaining reliable anharmonic spectroscopic properties. However, the method that employs complete basis set limit extrapolation ("C"), consideration of core electron correlation ("cC"), and inclusion of scalar relativity ("R") to produce the energy points on the QFF, the composite CcCR methodology, has not yet been utilized to study inorganic spectroscopy. This work takes the CcCR methodology and adapts it to test whether such an approach is conducive for the closed-shell, copper-containing molecules CuCN, CuOH, and CuCCH. Gas phase rovibrational data are provided for all three species in their ground electronic states. Equilibrium geometries and many higher-order rovibrational properties show good agreement with earlier studies. However, there are notable differences, especially in computation of fundamental vibrational frequencies. Even with further additive corrections for the inner core electron correlation and coupled cluster with full single, double, and triple substitutions (CCSDT), the differences are still larger than expected indicating that more work should follow for predicting rovibrational properties of transition metal molecules.

Metal Ion Modeling Using Classical Mechanics

Metal ions play significant roles in numerous fields including chemistry, geochemistry, biochemistry, and materials science. With computational tools increasingly becoming important in chemical research, methods have emerged to effectively face the challenge of modeling metal ions in the gas, aqueous, and solid phases. Herein, we review both quantum and classical modeling strategies for metal ion-containing systems that have been developed over the past few decades. This Review focuses on classical metal ion modeling based on unpolarized models (including the nonbonded, bonded, cationic dummy atom, and combined models), polarizable models (e.g., the fluctuating charge, Drude oscillator, and the induced dipole models), the angular overlap model, and valence bond-based models. Quantum mechanical studies of metal ion-containing systems at the semiempirical, ab initio, and density functional levels of theory are reviewed as well with a particular focus on how these methods inform classical modeling efforts. Finally, conclusions and future prospects and directions are offered that will further enhance the classical modeling of metal ion-containing systems.

Prediction of Bond Dissociation Energies/Heats of Formation for Diatomic Transition Metal Compounds: CCSD(T) Works

It was recently reported ( J. Chem. Theory Comput. 2015 , 11 , 2036 - 2052 ) that the coupled cluster singles and doubles with perturbative triples method, CCSD(T), should not be used as a benchmark tool for the prediction of dissociation energies (heats of formation) for the first row transition metal diatomics based on a comparison with the experimental thermodynamic values for a set of 20 diatomics. In the present work the bond dissociation energies as well as the heats of formation for those diatomics have been calculated by the Feller-Peterson-Dixon approach at the CCSD(T)/complete basis set (CBS) level of theory including scalar relativistic corrections and correlation of the outer shell of core electrons in addition to the valence electrons. Revised experimental values for the hydrides are presented that are based on new heterolytic R-H bond dissociation energies, which are needed for analysis of the mass spectrometry experiments. The agreement between the calculated bond dissociation energies and the revised experimental values of the hydrides is good. Good agreement of the calculated bond dissociation energies/heats of formation is also found for most of the chlorides, oxides, and sulfides given the experimental error bars from experiment and those of the transition metal atoms in the gas phase. Thus, reliable results can be achieved by the CCSD(T) method at the CBS limit. The use of PW91 orbitals for the CCSD(T) calculations improves the predictions for some compounds with large T₁ diagnostics at the HF-CCSD(T) level. The optimized bond distances and calculated vibrational frequencies for the diatomics also agree well with the available experimental values.

Solv-ccCA: Implicit Solvation and the Correlation Consistent Composite Approach for the Determination of pKa

Direct theoretical methods are advantageous for the prediction of pKa, as relative methods rely upon the experimental values of reference acid molecules that can limit application of the method to well-characterized systems. Here, a direct route is introduced, which incorporates the SMD universal solvation model1 within the correlation consistent Composite Approach (ccCA). This Solv-ccCA methodology has been used for the prediction of theoretical pKa values for nitrogen-containing species to within a mean absolute deviation (MAD) of 1.0 pKa unit from experimental values by utilizing a thermodynamic cycle that combines gas-phase and solution-phase calculations. Several density functionals, including B3LYP, B97-1, B97-2, B98, BMK, M06, and M06-2X, were also evaluated for use with SMD and for comparison to Solv-ccCA.

Is near-"spectroscopic accuracy" possible for heavy atoms and coupled cluster theory? An investigation of the first ionization potentials of the atoms Ga-Kr

Recent developments in ab initio coupled cluster (CC) theory and correlation consistent basis sets have ushered in an era of unprecedented accuracy when studying the spectroscopy and thermodynamics of molecules containing main group elements. These same developments have recently seen application to heavier inorganic or transition metal-containing species. The present work benchmarks conventional single reference coupled cluster theory (up to full configuration interaction for valence electron correlation and coupled cluster with up to full pentuple excitations (CCSDTQP) for core-valence correlation) and explicitly correlated coupled cluster methods [CC with single, double, and perturbative triple substitutions (CCSD(T)-F12)] for the atomic ionization potentials of the six 4p elements (Ga-Kr), a property with experimental error bars no greater than a few cm(-1). When second-order spin orbit coupling effects are included, a composite methodology based on CCSD(T) calculations yielded a mean signed error of just -0.039 kcal mol(-1) and a mean unsigned error of 0.043 kcal mol(-1). Inclusion of post-CCSD(T) correlation corrections reduced both of these values to -0.008 kcal mol(-1) and 0.025 kcal mol(-1), respectively, with the latter corresponding to an average error of just 9 cm(-1). The maximum signed error in the latter scheme was just -0.043 kcal mol(-1) (15 cm(-1)).

Explicitly Correlated Methods within the ccCA Methodology

The prediction of energetic properties within "chemical accuracy" (1 kcal mol(-1) from well-established experiment) can be a major challenge in computational quantum chemistry due to the computational requirements (computer time, memory, and disk space) needed to achieve this level of accuracy. Methodologies such as coupled cluster with single, double, and perturbative triple excitations (CCSD(T)) combined with very large basis sets are often required to reach this level of accuracy. Unfortunately, such calculations quickly become cost prohibitive as system size increases. Our group has developed an ab initio composite method, the correlation consistent Composite Approach (ccCA), which enables such accuracy to be possible, on average, but at reduced computational cost as compared with CCSD(T) in combination with a large basis set. While ccCA has proven quite useful, computational bottlenecks still occur. In this study, the means to reduce the computational cost of ccCA without compromising accuracy by utilizing explicitly correlated methods within ccCA have been considered, and an alternative formulation is described.