Calibrating Reaction Enthalpies: Use of Density ... - ACS Publications

12 Dec 2016 - Department of Chemistry and Computational Research on Materials Institute, The University of Memphis, Memphis, Tennessee. 38152, United ...
0 downloads 0 Views 1MB Size
Article pubs.acs.org/JPCA

Calibrating Reaction Enthalpies: Use of Density Functional Theory and the Correlation Consistent Composite Approach in the Design of Photochromic Materials Roger G. Letterman,† Nathan J. DeYonker,*,† Theodore J. Burkey,† and Charles Edwin Webster*,†,‡ †

Department of Chemistry and Computational Research on Materials Institute, The University of Memphis, Memphis, Tennessee 38152, United States ‡ Department of Chemistry and Center for Computational Sciences, Mississippi State University, Mississippi State, Mississippi 39762-9573, United States S Supporting Information *

ABSTRACT: Acquisition of highly accurate energetic data for chromium-containing molecules and various chromium carbonyl complexes is a major step toward calibrating bond energies and thermal isomerization energies from mechanisms for Cr-centered photochromic materials being developed in our laboratories. The performance of six density functionals in conjunction with seven basis sets, utilizing Gaussian-type orbitals, has been evaluated for the calculation of gas-phase enthalpies of formation and enthalpies of reaction at 298.15 K on various chromium-containing systems. Nineteen molecules were examined: Cr(CO)6, Cr(CO)5, Cr(CO)5(C2H4), Cr(CO)5(C2ClH3), Cr(CO)5(cis-(C2Cl2H2)), Cr(CO)5(gem-(C2Cl2H2)), Cr(CO)5(trans-(C2Cl2H2)), Cr(CO)5(C2Cl3H), Cr(CO)5(C2Cl4), CrO2, CrF2, CrCl2, CrCl4, CrBr2, CrBr4, CrOCl2, CrO2Cl2, CrOF2, and CrO2F2. The performance of 69 density functionals in conjunction with a single basis set utilizing Slater-type orbitals (STO) and a zeroth-order relativistic approximation was also evaluated for the same test set. Values derived from density functional theory were compared to experimental values where available, or values derived from the correlation consistent composite approach (ccCA). When all reactions were considered, the functionals that exhibited the smallest mean absolute deviations (MADs, in kcal mol−1) from ccCA-derived values were B97-1 (6.9), VS98 (9.0), and KCIS (9.4) in conjunction with quadruple-ζ STO basis sets and B97-1 (9.3) in conjunction with cc-pVTZ basis sets. When considering only the set of gas-phase reaction enthalpies (ΔrH°gas), the functional that exhibited the smallest MADs from ccCA-derived values were B97-1 in conjunction with cc-pVTZ basis sets (9.1) and PBEPBE in conjunction with polarized valence triple-ζ basis set/effective core potential combination for Cr and augmented and multiple polarized triple-ζ Pople style basis sets (9.5). Also of interest, certainly because of known cancellation of errors, PBEPBE with the least-computationally expensive basis set combination considered in the present study (valence double-ζ basis set/effective core potential combination for Cr and singly-polarized double-ζ Pople style basis sets) also provided reasonable accuracy (11.1). An increase in basis set size was found to have an improvement in accuracy for the best performing functional (B97-1).



INTRODUCTION

functional groups, steric effects, and the proximity of each functional group to the metal center preceding chelation. Also, the bistability of these organometallic photochromes requires an adequate free energy of activation (ΔG⧧ > ∼30 kcal mol−1) for rearrangement such that the complex will not thermally isomerize between chelates at room temperature. Hence, the ground-state bond dissociation enthalpy (BDE) required to cleave a metal−ligand bond is of particular interest to researchers to aid in the rational design of new organometallic photochromes.3,4,9−16 Our lab and others have used photo-

The organometallic photochromic systems under development in our laboratories utilize light-driven linkage isomerizations of tethered bifunctional chelates to selectively cleave the metal− ligand bond and induce a subsequent rearrangement to form a chelate with a different functional group.1−8 In an ideal photochromic system, the photoisomerization will induce spectral changes in the organometallic complex when the second functional group chelates to the metal. Initiation of the photolinkage isomerization is, in part, dependent on the metal−ligand bond strength of the chelated functional group. For example, near-UV irradiation tends to break the weakest bond in many of these complexes.2 The thermodynamics of the linkage isomerization are dependent on several factors, including the relative metal−ligand bond strengths of tethered © XXXX American Chemical Society

Received: September 13, 2016 Revised: November 9, 2016

A

DOI: 10.1021/acs.jpca.6b09278 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

relativistic effects, and an increased complexity of the bonding, (e.g., π back-bonding from the metal to a CO ligand). Wilson, DeYonker, and co-workers suggest that the target accuracy for computed enthalpies of formation should be within 3.0 kcal mol−1 for transition-metal complexes. Jiang, DeYonker, and Wilson applied ccCA-TM toward the calculation of the ΔfH°gas of 225 molecules containing the third-row d-block (elements Sc to Zn).29,30 They reported that transition-metal chemical accuracy (MAD less than or equal to 3.0 kcal mol−1) was achieved for a subset of 70 of these molecules with experimental uncertainties less than or equal to 2.0 kcal mol−1. It was also observed that the MADs decrease monotonically with decreasing thresholds in experimental uncertainties. It should also be noted that the experimental uncertainties in the measured ΔfH°gas for many of the species in the 225 molecule test set are greater than 3.0 kcal mol−1 (and in some cases much greater) or are unreported. The ccCA is computationally expensive and scales as N7 (where N is the number of basis functions); thus, many of the organometallic photochromes studied in our lab are simply too large to examine via ccCA-TM. DFT methods based on a SCFtype algorithm formally scale as N4, while density-fitting reduces the scaling to N3, but the reliability of a given functional for computationally derived BDEs can vary greatly depending on the chemical environment. Numerous studies have been conducted in recent years to evaluate the performance of DFT methods toward the prediction of various properties of transition-metal-containing systems.31−41 However, previous studies differ in their recommendations of the ideal DFT functional. For example, Furche and Perdew reported the computed bond energetics of 64 3d transitionmetal compounds examined as part of an investigation into the performance of semilocal and hybrid DFT functionals. They concluded, in part, that TPSS is a general purpose functional which exhibits nearly uniform accuracy for their selected systems.32 Landis and co-workers reported a systematic computational study of transition metal−ligand bond enthalpies via B3LYP and CCSD(T) for a wide assortment of ligands in saturated transition-metal hydride species covering the entire d-block (including 79 M−H and M−CH 3 compounds). The reported differences in the metal−ligand bonding within a row were significant.33 Riley and Merz examined the performance of 12 DFT functionals in the calculation of 94 ΔfH°gas and 58 ionization potentials (IP). They reported that the TPSS1KCIS functional in combination with the TZVP basis set yielded the lowest average errors in the calculation of ΔfH°gas for the transition-metal-containing species in the test set. Cr-containing systems provided the worst results, and carbonyl-bearing compounds were shown to be problematic for functionals without exact exchange.34 Wilson, Cundari, and co-workers used correlation-consistent basis sets to assess the performance of 44 DFT functionals toward the prediction of ΔfH°gas for 19 3d TM complexes with small (0.05−3.3 kcal mol−1) experimental uncertainties, including carbonyl-containing compounds. They reported significantly larger average MADs for compounds containing carbonyls (e.g., the avg MAD (from experimental values) for all meta-GGA functionals tested was greater than 80 kcal mol−1).35 Wilson and co-workers in a separate study tested the performance of 13 single- and double-hybrid density functionals in combination with correlation-consistent basis sets in the calculation of ΔfH°gas for a set of 193 3d transition-metalcontaining molecules. They reported that B97-1 provided

acoustic calorimetry (PAC) to study the BDEs of various ligands in transition-metal (TM) complexes.1−4,9,10,17−19 Several aspects of PAC techniques have led our group to search for an alternate means to determine metal−ligand bond strengths. For example, ignoring the effects of reaction volumes (defined as the difference in molar volumes between products and reactants) can cause significant error in the determination of BDEs via PAC. Establishing the contribution of reaction volume to the BDE requires numerous PAC experiments for a single complex, wherein the elastic parameters of the solvent are altered.3 Our attempts to address the reaction volume component via computational methods have been met with limited success because the errors in the experimental (ranging from 5 to 59 mL mol−1) and computational molar volumes (22.8 mL mol−1 standard error in the estimated molar volume) are larger than the experimentally determined reaction volumes (absolute value ranging from 0.9 to 14.4 mL mol−1 for that study).2 For example, at room temperature, a change in reaction volume of 3 mL mol−1 in heptane leads to an error of 1 kcal mol−1 in the BDE derived from PAC.2 One of our goals is to alleviate the need to perform these particularly expensive and time-consuming experiments. Use of an appropriate computational method could provide an alternative approach to obtaining BDEs. In the current work, we seek to ascertain a less-expensive density functional theory (DFT) method calibrated to a composite method known to provide accurate energetic properties for transition-metal-containing species. The correlation consistent composite approach (ccCA) developed by Wilson, Cundari, and co-workers is one example of a highly accurate composite method for determining atomic and molecular energies.20−25 The ccCA-derived energies can be used to calculate the gas-phase enthalpy of reaction (ΔrH°gas), including BDEs and enthalpies of formation (ΔfH°gas). The ccCA is based on complete basis set (CBS) extrapolations of MP2 energies and additive CCSD(T) energies. Geometries and vibrational frequencies are obtained at the B3LYP/cc-pVTZ level of theory. Other than the inherent parameters in the B3LYP functional, ccCA is an entirely ab initio model chemistry. The ccCA has been shown to have a mean absolute deviation (MAD) from experiment of 1.0 kcal mol−1 for 454 energetic properties in the main group G3/05 training set for atoms and molecules.20−24,26 The prediction of ΔfH°gas for transition-metal containing species has also been explored via ccCA (by means of a slightly altered methodology called ccCA-TM). Many factors contribute to potential inaccuracies in the experimental determination of enthalpies, for example, material instability and incomplete combustion in the case of enthalpy of formation or the determination of reaction volume for photocalorimetric methods in the case of reaction enthalpy. Additionally, much of the primary literature of known gas-phase inorganic enthalpies of formation, where such information is available, is of questionable quality. For example, ΔfH°gas values for only two of the nine organometallic complexes examined in this work have been determined experimentally, and the reported errors for Cr(CO)5 and Cr(CO)6 are ±3.1 and ±20 kcal mol−1, respectively.27 In some cases, thermochemical information provided in the literature lacks an adequate assessment of experimental errors or exhibits large uncertainties, thereby hindering accurate benchmarking of theoretical methods.28,29 For theory, several factors lead to significant inaccuracies in prediction of energetic properties for species which contain transition metals: larger valence electron space, stronger B

DOI: 10.1021/acs.jpca.6b09278 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A Table 1. Functionals Utilized for the Post-SCF Gradient Corrections to LDA Density Functional Energy GGA

meta-GGA

hybrid-GGA

hybrid-meta-GGA

OPerdew52,68,69 BOP49,73 BLYP48,49 BLYP-D79 OLYP48,68 OPBE53,68,86 XLYP88 BP8649,52,69 BP86-D79 HCTH/9351 HCTH/12097 HCTH/14797 HCTH/40799 PBE53,54 PBE-D79 mPBE104 revPBE105 RPBE107 mPW91 PW91110−112 FT97113 KT1114 KT2114 Becke88x+BR89c115−117

KCIS39 KCIS-modified74 mPBEKCIS76 PKZBx-KCIScor80 PKZB60,82 TPSS71,72 TPSS-D79 BmTau190 M06-L92 OLAP368,94 BLAP394 VS9898 Becke00100 τ-HCTH89

B9750,70 B97-151,70 B97-270,77 B97-D79 PBE070,83−85 OPBE053,54,68,86,87 τ-HCTH-hybrid70,89 mPW1PW70,91 mPW1K70,93 B1PW91(VWN5)70,95,96 B1LYP(VWN5)70,95,96 B3LYP-D79 B3LYP(VWN5)70,96,101 B3LYP*(VWN5)70,96,102 O3LYP(VWN5)103 X3LYP(VWN5)88 KMLYP(VWN5)70,96,106 BHandH48,70,108 BHandHLYP48,70,109

TPSSh70−72 M0575 M05-2X78 M0681 M06-2X81 mPBE0KCIS70,76 mPBE1KCIS70,76

experimental values. It is unclear how greater uncertainty in the calculation of ΔfH°gas for chromium-centered carbonyl-bearing compounds will affect the choice of DFT functional when examining the organometallic complexes of interest to our lab because calculating ΔfH°gas using the atomization scheme is a test of methodology much more stringent than calculating BDEs where chemical environments of the product and reactant are typically quite similar.45 Furthermore, we compare the DFT- and ccCA-derived enthalpies for systems where the experimental data are not available. This comparison will also be used to establish the fitness of various DFT methodologies for larger systems that are intractable for ccCA.

optimal results for Cr-containing species. They also noted that the results for Cr-containing species displayed the largest MADs (approximately 10−20 kcal mol−1) compared to all other metal-centered species, which is a conclusion also noted by Riley and Merz.34,36 Variations in the performance of DFT functionals, especially for Cr-centered organometallic complexes, may make the best choice of DFT method unclear in the absence of both reliable experimental data and benchmarking with highly accurate ab initio techniques. The atomization energies of CrO, CrO2, and CrO3 were evaluated by Chan, Raghavachari, Radom, and Karton via composite methods with post-CCSD(T) terms up to CCSDTQ(5)Λ and CCSDTQ(5). Heats of formation from best estimated atomization energies for CrO, CrO2, and CrO3 were determined to be 47.4 ± 1, −19.4 ± 1, and −68.5 ± 5 kcal mol−1, respectively.42 Zhang, Truhlar, and Tang reported the results of 42 exchange-correlation functionals (11 different types) used to predict the average bond energies of 70 (3d) transition-metalcontaining molecules, 19 of which were single-reference, and 51 were multireference. Larger errors on average were reported for middle transition metals (V, Cr, and Mn) than for either early or late transition metals.43 Xu, Zhang, Tang, and Truhlar tested and compared the results of coupled cluster (CC) and Kohn−Sham DFT methods on 20 3d transition-metal-containing diatomic molecules against the most reliable experimental data available. On average, nearly half of the exchange-correlation functionals tested performed better than CCSD(T) when utilizing the same extended basis set for a given molecule.44 In the current study, for mostly Cr-containing species with small experimental errors, we present a comparison of gasphase enthalpies of reaction (ΔrH°gas) [BDEs, exchange reaction enthalpies (ΔexH°gas), and enthalpies of formation (ΔfH°gas)] of ccCA and various DFT model chemistries to



COMPUTATIONAL METHODS Density Functional Theory. Computations were carried out using the Gaussian0946 implementations of B3LYP (the B347 exchange and the LYP correlation48), BHandH (the Becke half and half exchange correlation),49 B97-1 (the Becke B9750 exchange and correlation modified by Handy, Tozer, and coworkers51), BVP86 (the Becke exchange (B)49 and the Perdew correlation (VP86)52), and PBEPBE (the PBE exchange and correlation53,54) functionals of density functional theory55 using the default pruned fine grids for energies (75, 302), default pruned course grids for gradients and Hessians (35, 110) (neither grid is pruned for chromium), and default SCF convergence for geometry optimizations (10−8). Single-point calculations on B3LYP/BS3 geometries were also carried out using the Amsterdam Density Functional program (ADF, versions 2008.01 and 2009.01).56−58 ADF uses Slater-type orbital (STO) basis sets rather than Gaussian type orbitals (GTO) and contains a density fitting procedure using auxiliary functions (fit functions) for the evaluation of the Coulomb potential and molecular density. All reported ADF calculations were post-SCF gradient corrections to the LDA59 (local density approximation) density functional energy, utilized a singleC

DOI: 10.1021/acs.jpca.6b09278 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

correlation-consistent basis sets of Dunning and co-workers.124−128 BS7: The basis set for chromium was derived from the completely uncontracted Hay and Wadt basis set and ECP combination (LANL08).132,133,138,139 The basis set for all other atoms were the cc-pVTZ correlation-consistent basis sets of Dunning and co-workers.124−128 BS8 (used solely in single-point ADF calculations): The basis set for all atoms was the standard ADF ZORA QZ4P basis set (an uncontracted quadruple-ζ STO basis set with four polarization functions). ADF ZORA basis sets were optimized for use within ZORA relativistic calculations. The density fitting approximation140−143 for the fitting of the Coulomb potential was used for all BVP86/GTO and PBEPBE/GTO calculations; auxiliary density-fitting basis functions were generated automatically (by the procedure implemented in Gaussian 09) for the specified AO basis set.144 Spherical harmonic d functions were used throughout; that is, there are five angular basis functions per d function, etc. All structures were fully optimized, and analytical frequency calculations were performed on all structures to ensure a zeroth-order saddle point (a local minimum) was achieved. All DFT enthalpies reported were calculated using the respective functional with unscaled harmonic vibrational frequencies. All thermal corrections were performed at 298.15 K.

point grid accuracy factor (4.0), and used the default SCF convergence (10−8). The functionals utilized in the post-SCF gradient corrections to the LDA density functional energy are listed in Table 1.39,60 Relativistic effects were treated with the zeroth-order regular approximation (ZORA)61−66 with the default sum of atoms potential approximation (SAPA) potential.67 The frozen core approximation was used in all ADF computations. Description of the Correlation Consistent Composite Approach. The ccCA is a nonparameterized MP2-based model chemistry known to provide chemically accurate results. Unlike many other model chemistries, ccCA can predict reliable ΔfH°gas and BDEs for TM-containing molecules.20−22,28,30 The ccCA was developed to employ geometries optimized with B3LYP utilizing cc-pVTZ (BS3, defined in Table 2) followed by a series of single-point energy calculations for the general ccCATM formulation.28,29 A summary of the ccCA may be found in the Supporting Information. Table 2. Basis Set Combinations basis set

transition metal

other atoms

BS1 BS2 BS3 BS4 BS5 BS6 BS7 BS8

mod-LANL2DZ mod-LANL2DZ cc-pVTZ (Cr) mod-LANL2TZ(f) sc-stuttgart nonrel (Cr) sc-stuttgart fully relativistic LANL08 ZORA QZ4P

6-31G(d′) 6-31++G(d′,p′) cc-pVTZ 6-311++G(2df,2p) cc-pVTZ cc-pVTZ cc-pVTZ ZORA QZ4P



RESULTS Standard Enthalpy of Formation. For the 32 molecules relevant to our study (the 19 chromium-containing molecules and 13 ligand “fragment” molecules, listed in Table 3), the gasphase enthalpy of formation (ΔfH°gas) was calculated in the manner described by Curtiss and co-workers145 from total electronic energies computed via ccCA and the temperaturedependent enthalpy correction computed via B3LYP/BS3. MAD is the Mean Absolute Deviation, MSD is the Mean Signed Deviation, and RMSD is the Root-Mean-Square Deviation. The ccCA-derived ΔfH°gas values are shown in Table 3 along with experimental measurements, where available. Next, the ΔfH°gas was calculated for these 32 molecules via GTO methods (B3LYP, BHandH, BVP86, PBEPBE, B97-1, and ωB97X-D functionals utilizing BS1 through BS7) and via STO methods (64 DFT functionals listed in Table 1 utilizing BS8 and the temperature-dependent enthalpy corrections computed via B3LYP/BS3). These DFT-derived ΔfH°gas values were compared to ccCA-derived values, and the complete sets of GTO and STO results may be found in the Supporting Information in Tables S1−S6 and S7−S18, respectively. Summaries of the MADs for the ΔfH°gas of each GTO and STO method from ccCA are shown in Tables 4 and 5, respectively. Standard Enthalpy of Reaction. The gas-phase standard enthalpy of reaction (ΔrH°gas) for each reaction shown in Table 6 was calculated as shown in eq 5 from ccCA-derived ΔfH°gas values.

Basis Sets. Several basis set combinations were utilized (see the following description and Table 2). In each case, where an all-electron main group cc-pVTZ basis set was used, the redundant functions were removed and then the basis set was linearly transformed, as suggested by Davidson.118 BS1: The basis set for chromium (341/341/41) was the Hay and Wadt basis set and effective core potential (ECP) combination (LANL2DZ)47 as modified by Couty and Hall, where the two outermost p functions have been replaced by a (41) split of the optimized chromium 4p function.119 The 631G(d′)120 basis sets were used for all other atoms. BS2: The metal basis set of BS1 was used for BS2. The 6-31+ +G(d′,p′)121−123 basis sets were used for all other atoms. BS3: The basis sets for all atoms were the cc-pVTZ correlation-consistent basis sets of Dunning and co-workers and Peterson and co-workers.124−130 BS4: The basis sets for chromium were the Hay and Wadt basis set and effective core potential (ECP) combination (LANL2TZ(f))131−133 that was uncontracted to [4s4p3d]. The 6-311++G(2df,2p)134,135 basis sets developed by Pople and coworkers were used for all other atoms. BS5: The basis set for chromium136 was the nonrelativistic contracted basis set developed by Peterson, Dolg, and Stoll and ECP (Stuttgart/Cologne group ECP10MHF118,137) combination. The basis sets for all other atoms were the cc-pVTZ correlation-consistent basis sets of Dunning and co-workers.124−128 BS6: The basis set for chromium136 was the fully relativistic contracted basis set developed by Peterson, Dolg, and Stoll and ECP (Stuttgart/Cologne group ECP10MDF118,137) combination. The basis sets for all other atoms were the cc-pVTZ

° = Δr Hgas

° ∑ ΔHgas

products



° ∑ ΔHgas

reactants

(5)

The ccCA-derived ΔrH°gas values are shown in Table 6 along with experimental measurements, where available. Next, the ΔrH°gas was calculated for each reaction via GTO methods (B3LYP, BHandH, BVP86, PBEPBE, B97-1, and ωB97X-D functionals utilizing BS1 through BS7) and via STO D

DOI: 10.1021/acs.jpca.6b09278 J. Phys. Chem. A XXXX, XXX, XXX−XXX

The Journal of Physical Chemistry A



Table 3. Experimental and ccCA-Derived ΔfH°gas with MSD, MAD, and RMSD Displayed for the Set of Species Presented (kcal mol−1)a species

experiment

ccCA

Cr(CO)5 Cr(CO)6 Cr(CO)5(C2H4) Cr(CO)5(C2H3Cl) Cr(CO)5(C2H2Cl2)b Cr(CO)5(C2H2Cl2)c Cr(CO)5(C2H2Cl2)d Cr(CO)5(C2HCl3) Cr(CO)5(C2Cl4) CO C2H4 C2H3Cl C2H2Cl2b C2H2Cl2c C2H2Cl2d C2HCl3 C2Cl4 CrBr2 CrBr4 CrCl2 CrCl4 CrF2 CrO2 CrO2Cl2 CrO2F2 CrOCl2 CrOF2 H2 O2 F2 Cl2 Br2

−153.9 (3.1) −240.0 (1.1)27,147 n.d. n.d. n.d. n.d. n.d. n.d. n.d. −26.4 (0.007)27 12.6 (0.036)148−150 5.3 (0.08)148−150 −0.7 (0.5)27 0.5 (0.33)27 −0.1 (0.5)27 −4.2 (0.5)27 −5.8 (1.0)27 −14.1 (4.3)151 −51.3 (4.5)151 −28.1 (0.4)152 −94.8 (3.3)152 −103.2 (3.0)152 −23.7 (1.2)153 −124.1 (1.0)152 −201.6 (5.0)152 −80.4 (5.4)152 −153.3 (3.4)152 0.0 0.0 0.0 0.0 7.4 (0.026)27 MSD MAD RMSD

−170.0 −239.5 −187.6 −192.1 −195.4 −193.6 −196.5 −198.0 −199.0 −26.3 12.7 5.3 −0.1 0.6 0.5 −3.4 −5.5 −9.1 −53.9 −32.1 −100.8 −107.4 −15.9 −124.1 −202.0 −79.7 −157.6 −0.5 0.8 0.2 0.2 8.0 0.8 2.3 4.2

27,146

Uncertainty values, when known, are given in parentheses; n.d. denotes no determined or known experimental value. bcis-Olefin. c Geminal-olefin. dtrans-Olefin.

Table 4. Summary of MADs (kcal mol−1) for ΔfH°gas from ccCA for Each GTO Method Listed B3LYP B97-1 BHandH BVP86 PBEPBE ωB97X-D

BS2

BS3

BS4

BS5

BS6

BS7

39.1 19.1 29.1 18.1 34.7 35.8

47.8 25.9 31.8 12.3 24.3 42.4

29.2 9.2 20.6 23.9 39.0 21.9

31.6 8.3 24.6 22.0 35.6 24.4

47.8 26.0 35.5 23.8 35.7 45.0

31.8 9.4 21.2 20.2 34.8 27.2

32.4 9.9 23.9 21.3 35.6 27.3

DISCUSSION

Standard Enthalpy of Formation. DFT methods often exhibit large errors in the exchange and correlation energies. Furthermore, the exchange energy may be underestimated, while the correlation energy is overestimated, leading to a cancellation of errors.154,155 The prediction of reaction enthalpies where the products and reactants exhibit very similar electronic structures is a less difficult task compared to predicting enthalpies of formation from standard state elements. The computation of ΔfH°gas from atomization energies is a daunting task for computational methods because ΔfH°gas is typically calculated from gas-phase atomization energies, and the electronic structures of the product and reactants are quite different. Therefore, inaccuracies in the prediction of this property are addressed first. While data on gas phase enthalpy of formation under standard conditions (ΔfH°gas) are unavailable for the Cr(CO)5olefin complexes investigated in the course of this work, ΔfH°gas values of −218 ± 20 and −153.9 ± 3.1 kcal mol−1 have been reported for Cr(CO)6 and Cr(CO)5, respectively.27,158 The value for Cr(CO)6 published in the NIST WebBook represents the average of nine ΔfH°gas values; however, the ccCA-derived ΔfH°gas for Cr(CO)6 (−239.5 kcal mol−1) compares favorably to the most exothermic of the nine measurements (−240.4 ± 1.1 kcal mol−1). This value is adopted as the experimental ΔfH°gas in this work. The ccCA-derived ΔfH°gas for Cr(CO)5 (−170.0 kcal mol−1) varies from experiment by 16.1 kcal mol−1 and represents the greatest deviation from experiment for any ccCA-derived ΔfH°gas examined in the course of this study. This significant deviation in the ccCA-derived and experimental values for Cr(CO)5 is a statistical outlier; however, the experimental ΔfH°gas for Cr(CO)5 published in the NIST WebBook relies on the ΔfH°gas for Cr(CO)6 of −217.1 ± 0.62 kcal mol−1 and, as previously noted, the ΔfH°gas for Cr(CO)6 may not be reliable.28,29 The MAD of ccCA-derived versus experimental ΔfH°gas shown in Table 3 for the 25 molecules where experimental values are known is 2.3 kcal mol−1. If the ΔfH°gas value for Cr(CO)5 is excluded, the MAD is improved to 1.7 kcal mol−1. The ccCA-derived and experimental ΔfH°gas values for each species (Table 3) were compared, and the MSD, MAD, and RMSD were determined to be 0.8, 2.3, and 4.2 kcal mol−1, respectively, when including Cr(CO)5. When Cr(CO)5 is removed from the data set, the MSD, MAD, and RMSD improved to 0.2, 1.7, and 2.8 kcal mol−1, respectively. As shown in Figure 1, the ccCA-derived and experimental ΔfH°gas values were well-correlated, as indicated by the R2 value of 0.997, a standard error in the y-estimate of 3.8 kcal mol−1, and a standard error in the y-intercept of 0.94 kcal mol−1. ΔfH°gas via Methods Employing GTO Basis Sets. Only the B97-1 functional produced a MAD below 10 kcal mol−1 when considering the entire test set with the best overall results corresponding to BS4, BS3, BS6, and BS7 having MADs of 8.3, 9.2, 9.4, and 9.9 kcal mol−1, respectively. Interestingly, Wilson and co-workers reported that the B97-1 functional produced a MAD of 8.0 kcal mol−1 at the B97-1/aug-cc-pVQZ//B97-1/ aug-cc-pVTZ level of theory for 25 complexes with coordinated bonds or organic ligands from the ccCA-TM/11 test set.36 They also reported a MAD of 6.9 kcal mol−1 for the B97-1 functional when paired with the cc-pVTZ basis set considering 19 TM species. The MAD was improved (3.1 kcal mol−1) when the size of the basis set was increased to cc-pVQZ.35 The next

a

BS1

Article

methods (64 DFT functionals, given in Table 1, utilizing BS8 and the temperature-dependent enthalpy corrections computed via B3LYP/BS3). These DFT-derived ΔrH°gas values were compared to ccCA-derived values. The GTO and STO results are shown in Tables S19−S24 and S25−S33, respectively. Summaries of the MADs for each GTO and STO method from ccCA are shown in Tables 7 and 8, respectively. E

DOI: 10.1021/acs.jpca.6b09278 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A Table 5. Summary of MADs (kcal mol−1) for ΔfH°gas from ccCA for Each STO Method Listeda

a

GGA/BS8

MAD

meta-GGA/BS8

MAD

hybrid/BS8

MAD

revPBE RPBE HCTH/407 Becke88x+BR89c HCTH/93 OLYP OPBE XLYP HCTH/147 FT97 BLYP-D BLYP HCTH/120 BOP OPerdew KT2 mPW BP86-D BP86 mPBE PBE PBE-D PW91 KT1

8.5 8.7 9.0 10.1 11.0 11.3 12.4 13.2 13.7 13.9 14.3 14.3 14.6 15.0 18.8 28.6 28.6 32.4 32.4 35.2 46.7 46.7 48.1 56.2

KCIS BLAP3 VS98 PKZB TPSS-D TPSS BmTau1 Becke00 M06-L PKZBx-KCIScor KCIS-modified OLAP3 τ-HCTH mPBEKCIS

8.0 9.4 11.6 12.0 12.8 12.8 13.0 13.0 13.1 15.1 15.8 17.4 18.1 42.2

B97-1 B97 B97-2 B97-D B3LYP*(VWN5) M06 M05 PBE0 mPBE1KCIS mPBE0KCIS TPSSh τ-HCTH-hybrid B3LYP-D B3LYP(VWN5) O3LYP(VWN5) X3LYP(VWN5) mPW1PW B1PW91(VWN5) BHandH B1LYP(VWN5) M05-2X M06-2X KMLYP(VWN5) OPBE0 mPW1K BHandHLYP

5.5 10.1 10.3 11.3 14.5 15.4 16.6 18.5 20.0 20.8 21.0 21.2 22.8 23.2 25.9 26.4 28.3 32.2 33.8 41.5 43.2 48.1 52.0 55.2 65.8 90.8

All DFT/BS8(STO) computations utilized B3LYP/BS3 optimized geometries and temperature-dependent enthalpy corrections.

the hybrid functionals (B3LYP and BHandH) for the prediction of ΔfH°gas in the subset of small Cr-containing species, yet B97-1 provides the best performance.32 A comparison of the MSD, MAD, and RMSD for the hybrid functionals (B3LYP and BHandH) and pure functionals (BVP86 and PBEPBE) with regard to ΔfH°gas shows that the levels of theory utilizing pure functionals are more accurate. The hybrid functionals also display a bias toward underestimating ΔfH°gas, while the pure functionals display a bias for overestimating ΔfH°gas, as shown by the MSDs reported for each functional. The MADs of the computed ΔfH°gas for all DFT methods employing STO basis sets examined in this study are shown in Figure S3. The most accurate level of theory is B97-1/BS4 with a MAD of 8.3 kcal mol−1 compared to ccCA-derived values. Three other levels of theory are also competitive: B97-1/BS3, B97-1/BS6, and B97-1/BS7 with MADs of 9.2, 9.4, and 9.7 kcal mol−1, respectively. A plot of the ΔfH°gas computed via B97-1/ BS4 versus ccCA (Figure S4) demonstrates that these values are well-correlated, as shown by the R2 value of 0.994. ΔfH°gas via Methods Employing STO Basis Sets. Three GGA functionals (revPBE, RPBE, and HCTH/407), two metaGGA functionals (KCIS and BLAP3), and one hybrid functional (B97-1) produce MADs of less than 10.0 kcal mol−1 when considering the entire test set, with the lowest MAD (5.5 kcal mol−1) from the B97-1 functional. The differences between ΔfH°gas calculated via DFT methods utilizing STO basis sets and those calculated via ccCA for each species along with corresponding MSD, MAD, and RMSD statistics (Tables S7−S18) were compared. The metrics of energetic data for functionals utilizing the generalized gradient approximation are displayed in Figures S5

best levels of theory are BVP86/BS2 and BVP86/BS1 with MADs of 12.3 and 18.1 kcal mol−1, respectively. The least accurate level of theory is B3LYP/BS5 with MAD of 47.8 kcal mol−1. B3LYP results were unfortunately poor, considering this is arguably among the most widely used DFT functionals. Wilson, Cundari, and co-workers reported that hybrid functionals (which include various amounts of HF exchange) performed better in the determination of ΔfH°gas for the 19 TM species they examined.35 However, Jiang and Wilson have more recently shown that B3LYP exhibited the worst performance of the 13 functionals tested with correlation-consistent basis sets for the 34 Cr-containing species in their test set.36 Riley and Merz concluded that DFT methods used to calculate ΔfH°gas give errors for TM compounds much larger than those for organic ones, at least for the 6-31G** basis set.34 This conclusion proved true for each of the functionals tested in this work with the exception of the BHandH functional. This functional, perhaps fortuitously, provides MADs for the olefins larger than those for the organometallic complexes when paired with BS1 through BS7. The MSD, MAD, and RMSD for each level of theory shown in Tables S1−S6 are summarized as vertical bar charts; see Figure S1 for the hybrid functionals B3LYP, B97-1, and BHandH and Figure S2 for the pure functionals BVP86 and PBEPBE. For DFT methods paired with GTO-type basis sets, we find that MSDs vary from −47.8 to 37.8 kcal mol−1, MADs vary from 12.3 to 47.8 kcal mol−1, and RMSDs vary from 15.1 to 61.3 kcal mol−1. Furche and Perdew suggested that less-exact HF exchange is needed to describe systems with stronger correlation when determining ΔfH°gas. Therefore, it is interesting to note that the pure functionals (BVP86 and PBEPBE) in this study provide performance better than that of F

DOI: 10.1021/acs.jpca.6b09278 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

The MADs of the computed ΔfH°gas for all DFT methods employing STO basis sets examined in this study are shown in Figure S10. The B97-1/BS8//B3LYP/BS3 level of theory is most accurate, exhibiting a MAD of 5.5 kcal mol−1 compared to ccCA-derived values. A plot of the ΔfH°gas computed via B971/BS8//B3LYP/BS3 versus ccCA (Figure S11) demonstrates that these values are well-correlated, as indicated by the R2 value of 0.995. Standard Enthalpy of Reaction. In the course of our development of photochromic organometallic complexes, we typically study isomerizations involving the reversible formation of two chelates via the breaking and formation of single bonds in the tethered bifunctional moiety of a polydentate ligand.1,4−6,8 Therefore, a comparison of reaction enthalpies (ΔrH°gas) may be a better metric than enthalpies of formation/ atomization energies for the determination of appropriate methods to model organometallic mechanisms and bond dissociation enthalpies in our lab. When computing reaction energies or reaction enthalpies, differences in the number of bonds and the nature of the bonding in products versus reactants may offset, thereby leading to cancellation of electronic structure errors. Wheeler et al. showed that, within the hierarchy framework of hydrocarbon reaction classes, the hyperhomodesmotic and homodesmotic reactions produced the least error for the closed-shell molecules examined.45 In other words, methodological dependencies are mostly eliminated via cancellation of error when products and reactants come from an increasingly similar chemical environment. A systematic hierarchy of reaction energies has not been derived for transition-metal species. However, the underlying principles pertaining to the number and type of bond in products and reactants as applied to closed-shell hydrocarbons holds true when the comparison of reaction enthalpies and bond dissociation enthalpies of Cr-containing species using ccCA and various DFT is made. The first step in the formation of these organometallic photochromic complexes is the photodissociation of a carbonyl ligand from the tricarbonyl parent, thereby allowing the coordination of a tethered ligand and the resultant formation of the desired chelate. This process is similar to the first photodissociation of CO from Cr(CO)6 and subsequent coordination of solvent or other ligand. The gas-phase BDE of CO from Cr(CO)6 was reported to be 36.8 ± 2 kcal mol−1 using laser pyrolysis by Smith and co-workers,156 and Fletcher and Rosenfeld have used time-resolved infrared laser absorption spectroscopy to estimate a BDE of 37.1 ± 5.0 kcal mol−1.157 The ccCA-derived ΔrH°gas for this reaction is 43.3 kcal mol−1, a difference of 6.5 kcal mol−1 from the experimental value with the smallest reported uncertainty. Experimental BDEs were reported for a series of Cr(CO)5ligand complexes [ligand = 1,1-dichloroethylene (1,1-C2H2Cl2), 1,1,2-trichloroethylene (C2HCl3), and 1,1,2,2-tetrachloroethylene (C2Cl4)] by Cedeño and co-workers via laser photoacoustic calorimetry (LPAC) in n-hexane.158 BDEs for these complexes were determined to be 14.3, 13.7, and 13.5 kcal mol−1, respectively. However, the authors did not utilize the correct quantum yield in the calculation of BDEs. The authors also erroneously assumed that displacement of solvent by olefin would result in no change in reaction volume. PAC results presented by Burkey and co-workers suggest that the reaction volume for displacement of heptane solvent by an olefin in Mo(CO)5-L is not zero.2 Lastly, the BDEs reported by Cedeño and co-workers lacked experimental uncertainties; therefore,

Table 6. Experimental and ccCA-Derived Gas-Phase Enthalpy of Reaction with MSD, MAD, and RMSD Displayed for the Set of Species Presented (kcal mol−1)a reaction

experiment

ccCA

Cr(CO)6 → Cr(CO)5 + CO Cr(CO)5(C2H4) → Cr(CO)5 + C2H4 Cr(CO)5(C2H3Cl) → Cr(CO)5 + C2H3Cl Cr(CO)5(C2H2Cl2) → Cr(CO)5 + C2H2Cl2b Cr(CO)5(C2H2Cl2) → Cr(CO)5 + C2H2Cl2c Cr(CO)5(C2H2Cl2) → Cr(CO)5 + C2H2Cl2d Cr(CO)5(C2HCl3) → Cr(CO)5 + C2HCl3 Cr(CO)5(C2Cl4) → Cr(CO)5 + C2Cl4 CrO2 → Cr + O2 CrF2 → Cr + F2 CrCl2 → Cr + Cl2 CrBr2 → Cr + Br2 CrO2Cl2 → Cr + O2 + Cl2 CrOCl2 → Cr + 1/2O2 + Cl2 CrOCl2 + F2 → CrOF2 + Cl2 CrO2F2 + Cl2 → CrO2Cl2 + F2 CrO2 + Br2 → CrBr2 + O2 CrO2 + F2 → CrF2 + O2 CrO2 + Cl2 → CrCl2 + O2 CrF2 + Cl2 → CrCl2 + F2 CrCl2 + Br2 → CrBr2 + Cl2 CrF2 + Br2 → CrBr2 + F2 CrCl2 + O2 + F2 → CrO2F2 + Cl2 CrF2 + O2 + Cl2 → CrO2Cl2 + F2

36.8 (3.0) 25.7 (1.0) n.d. n.d. 14.3 (1.0) n.d. 13.7 (0.5) 12.8 (1.6) 119.2 (2.2) 198.7 (4.0) 123.6 (1.4) 116.9 (5.3) 219.6 (2.0) 175.9 (4.6) −72.9 (8.8) 77.5 (6.0) 2.2 (5.5) −79.5 (4.2) −4.4 (1.6) 75.1 (3.4) 6.6 (4.7) 81.7 (7.3) −173.5 (5.4) −20.9 (4.0) MSD MAD RMSD

43.3 30.3 27.4 25.3 24.2 26.9 24.7 23.5 111.8 202.7 127.3 112.1 220.1 175.2 −78.0 78.0 −0.3 −90.9 −15.5 75.4 15.2 90.6 −170.8 −17.5 −0.6 4.9 6.3

a

Uncertainty values, when known, are given in parentheses; n.d. denotes no determined or known experimental value. bcis-Olefin. c Geminal-olefin. dtrans-Olefin.

Table 7. Summary of MADs (kcal mol−1) for ΔrH°gas from ccCA for Each GTO Method Listed B3LYP B97-1 BHandH BVP86 PBEPBE ωB97X-D

BS1

BS2

BS3

BS4

BS5

BS6

BS7

20.8 17.8 33.2 12.4 11.1 22.5

20.4 16.9 31.1 11.9 10.8 21.0

13.3 9.3 21.1 12.2 11.4 12.6

16.2 11.9 24.2 10.4 9.5 15.6

28.9 25.4 41.1 18.0 17.1 29.7

15.0 10.3 21.7 11.2 10.6 15.4

16.5 12.6 24.1 11.0 10.1 16.4

and S6. For DFT methods utilizing GGA functionals and STO basis sets, we find that the MSDs varied from −4.8 to 56.1 kcal mol−1, MADs varied from 8.5 to 56.2 kcal mol−1, and RMSDs varied from 11.8 to 69.2 kcal mol−1, as shown in Tables S8−S11. A similar comparison for functionals utilizing metageneralized gradient approximations and hybrid approximations was also performed. For DFT methods utilizing meta-GGA functionals and STO basis sets, we find that the MSDs vary from −13.2 to 42.0 kcal mol−1, MADs vary from 8.0 to 42.2 kcal mol−1, and RMSDs vary from 11.4 to 55.4 kcal mol−1. The metrics for meta-GGA functionals are summarized for each method in Figure S7. DFT methods utilizing hybrid functionals and STO basis sets exhibit MSDs, MADs, and RMSDs from −83.8 to 18.0, 5.5 to 90.8, and 7.4 to 124.1 kcal mol−1, respectively. The statistics for hybrid functionals are summarized in Figures S8 and S9. G

DOI: 10.1021/acs.jpca.6b09278 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A Table 8. Summary of MADs (kcal mol−1) for ΔrH°gas from ccCA for Each STO Method Listed GGA/BS8

MAD

meta-GGA/BS8

MAD

hybrid/BS8

MAD

KT2 BP86-D PBE-D KT1 PW91 PBE HCTH/120 BLYP-D mPBE HCTH/147 BP86 mPW OPerdew RPBE revPBE Becke88x+BR89c HCTH/407 OPBE OLYP HCTH/93 BLYP XLYP BOP FT97

9.1 9.9 10.2 12.1 12.4 12.4 12.4 12.6 12.7 12.9 13.0 13.4 13.5 13.7 13.8 14.3 14.5 15.1 15.2 15.2 16.2 16.2 17.7 18.2

VS98 TPSS-D PKZB KCIS TPSS PKZBx-KCIScor Becke00 KCIS-modified M06-L τ-HCTH mPBEKCIS BLAP3 BmTau1 OLAP3

5.6 8.5 11.1 11.3 11.4 11.8 12.1 12.3 12.7 12.9 14.4 17.7 17.9 20.7

TPSSh B97-1 B3LYP-D M06 B97-D B97 B3LYP*(VWN5) M05 B97-2 mPBE1KCIS B3LYP(VWN5) X3LYP(VWN5) PBE0 O3LYP(VWN5) mPBE0KCIS mPW1PW B1LYP(VWN5) B1PW91(VWN5) BHandH OPBE0 M06-2X M05-2X mPW1K KMLYP(VWN5) BHandHLYP τ-HCTH-hybrid

7.8 8.9 9.5 9.5 10.1 10.4 10.6 10.8 11.8 11.9 13.1 13.2 14.8 15.9 16.4 16.6 16.7 17.6 22.1 24.0 26.6 27.4 27.9 29.8 32.3 220.8

respectively. It is not surprising that the gas-phase computed and solvated experimental results differ by 10 kcal mol−1 or more for aforementioned reasons. Moreover, our group recently performed a comparison of computed and experimental substitution enthalpies for a series of chromium arene tricarbonyl species in which a CO was substituted for a pyridyl ligand.1 The computed ΔrH°gas values were consistently greater (12 to 21 kcal mol−1) than the LPAC results in heptane and, in these experiments, the contributions of reaction volume, while recognized, were not determined. While the trends in the LPAC results were consistent with the reaction enthalpy trends, it was clear the absolute results were not an accurate representation of ΔrH°gas. It is necessary to consider the change in volume between reactants and products in the accurate determination of substitution enthalpies via LPAC. We concluded that the coordination of solvent plays an important role in the determination of the substitution enthalpy; therefore, ΔrH°gas cannot be compared directly with the LPAC experiments in the condensed phase. The aforementioned experimentally derived BDE of 13.5 kcal mol−1 for Cr(CO)5-(C2Cl4) determined directly from LPAC in n-hexane was in close agreement with a value of 12.8 ± 1.6 kcal mol−1 determined for gas phase BDE for this bond via transient infrared spectroscopy; however, this may be serendipitous.159 Gas-phase BDEs of 25.7 ± 1 and 24.7 ± 2.4 kcal mol−1 were reported for the Cr(CO)5-(C2H4) bond as identified via timeresolved infrared laser absorption spectroscopy by Grant and co-workers and Weitz and co-workers, respectively.160,161 We adopted the gas-phase values with the least experimental uncertainty, where available, for this study. A series of gas-phase decomposition and exchange reaction enthalpies were also calculated utilizing eight of the small Crcontaining species in the ΔfH°gas test set. These reactions are

Figure 1. Experimental ΔfH°gas versus ccCA-derived values.

Figure 2. Calculated values of ΔrH°gas derived via ccCA-derived versus experimentally determined values.

these values are used herein with reservation due to the paucity of experimentally determined BDEs for the organometallic complexes examined in this work. The ccCA-derived ΔrH°gas values for these reactions are 25.5, 25.8, and 24.3 kcal mol−1, H

DOI: 10.1021/acs.jpca.6b09278 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

Figure 3. MADs from ccCA values (kcal mol−1) for ΔH°gas when all reactions are considered, computed via B3LYP, B97-1, BHandH, BVP86, PBEPBE, and ωB97X-D utilizing BS1−BS7.

Table 9. Summary of MADs (kcal mol−1) for ΔH°gas (All Reactions) from ccCA for Each GTO Method Listed B3LYP B97-1 BHandH BVP86 PBEPBE ωB97X-D

BS1

BS2

BS3

BS4

BS5

BS6

BS7

31.3 18.5 30.8 15.7 24.6 30.1

36.0 22.0 31.5 12.1 18.5 33.2

22.4 9.3 20.8 18.9 27.2 17.9

25.0 9.8 24.4 17.0 24.4 20.7

39.7 25.7 37.9 21.3 27.8 38.4

24.6 9.8 21.4 16.4 24.4 22.2

25.6 11.0 24.0 16.9 24.6 22.6

uncertainties in the ΔfH°gas were propagated to determine the total uncertainty for the ΔrH°gas. The average uncertainty for the test set of 16 reactions involving small Cr-containing species was 4.6 kcal mol−1 with the largest uncertainty (8.8 kcal mol−1) obtained from the substitution of Cl2 by F2 in CrOCl2. The ccCA-derived values of ΔrH°gas compare favorably with available experimental values for the reactions shown in Table 6, and the MSD, MAD, and RMSD for this set of reactions were −1.7, 5.0, and 6.6 kcal mol−1, respectively. Again, in the same manner as in Figure 1, the ccCA-derived versus experimental ΔrH°gas values shown in Figure 2 were wellcorrelated, as demonstrated by the R2 value of 0.995, a standard error in the y-estimate of 6.9 kcal mol−1, and a standard error in the y-intercept of 1.7 kcal mol−1

listed in Table 6 with the ccCA-derived and experimental ΔrH°gas. The experimental values for these reactions involving small Cr-containing species were calculated as the difference between the Δ f H° gas of products and reactants. The

Table 10. Summary of MADs (kcal mol−1) for ΔH°gas (All Reactions) from ccCA for Each STO Method Listed GGA/BS8

MAD

meta-GGA/BS8

MAD

hybrid/BS8

MAD

revPBE RPBE HCTH/407 Becke88x+BR89c HCTH/93 OLYP HCTH/147 OPBE BLYP-D HCTH/120 XLYP BLYP FT97 BOP OPerdew KT2 mPW BP86-D BP mPBE PBE-D PBE PW91 KT1

10.8 10.8 11.4 11.9 12.8 13.0 13.3 13.5 13.6 13.7 14.5 15.1 15.7 16.2 16.5 20.2 22.1 22.7 24.1 25.6 31.0 32.0 32.8 37.3

VS98 KCIS-original TPSS-D PKZB TPSS Becke00 BLAP3 M06-L PKZBx-KCIScor KCIS-modified BmTau1 TPSSh τ-HCTH mPBE1KCIS OLAP3 mPBEKCIS

9.0 9.4 10.9 11.6 12.2 12.6 12.9 12.9 13.7 14.3 15.1 15.3 15.9 16.5 18.8 30.3

B97-1 B97 B97-D B97-2 B3LYP*(VWN5) M06 M05 PBE0 B3LYP-D B3LYP(VWN5) mPBE0KCIS X3LYP(VWN5) O3LYP(VWN5) mPW1PW B1PW91(VWN5) BHandH B1LYP(VWN5) M05-2X M06-2X OPBE0 KMLYP(VWN5) mPW1K BHandHLYP τ-HCTH-hybrid

6.9 10.2 10.8 10.9 12.8 12.8 14.1 16.9 17.1 18.8 18.9 20.7 21.6 23.3 25.9 28.8 30.9 36.5 38.9 41.8 42.5 49.6 65.7 106.7

I

DOI: 10.1021/acs.jpca.6b09278 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

Figure 4. MADs from ccCA values (kcal mol−1) for ΔH°gas of all reactions computed via methods employing the STO basis set (DFT/BS8// B3LYP/BS3). BHandHLYP and τ-HCTH-hybrid functionals were excluded due to poor performance.

functionals BVP86 and PBEPBE (Figure S13). The average MAD, calculated as the mean of the MADs for BS1−BS7 for each functional, shows that the pure functionals are more accurate on average than the hybrid functionals examined in this study. For example, the average MADs for the BVP86 and PBEPBE model chemistries are 11.5 and 12.4 kcal mol−1, respectively, compared to 18.7, 14.9, and 28.1 kcal mol−1 for the B3LYP, B97-1, and BHandH model chemistries, respectively. All model chemistries for DFT methods utilizing GTO basis sets display a bias toward overestimating the ΔrH°gas, as shown by a positive MSD. The MADs of all methods utilizing GTO basis sets employed to calculate ΔrH°gas are displayed in the vertical bar chart shown in Figure S14. B97-1/BS3 is the most accurate model chemistry that utilizes GTO basis sets and exhibits a MAD of 9.3 kcal mol−1 compared to ccCA-derived values. The accuracy of PBEPBE/BS4 with a MAD of 9.5 kcal mol−1 is also competitive. The MADs for all other model chemistries utilizing GTO basis sets range from 10.1 to 41.1 kcal mol−1. The PBEPBE/BS1 model chemistry is the most accurate of the model chemistries employing BS1 (containing the fewest number of basis functions) with a MAD of 11.1 kcal mol−1. For larger systems or in instances of limited computational resources, this model chemistry provides the best balance of efficiency and accuracy. The correlation between values of ΔrH°gas calculated via the B97-1/BS3 model chemistry versus ccCA values for the reactions displayed in Table 4 is shown in Figure S15. The R2 value of 0.988 shows that these data sets are well-correlated. ΔrH°gas via Methods Employing STO Basis Sets. ΔrH°gas values calculated via DFT methods utilizing STO basis sets were compared to those calculated via ccCA for each species, and the differences are shown along with statistical deviations in Tables S25−S33. We find that the MSDs vary from −2.9 to 11.4 kcal mol−1, MADs vary from 9.1 to 18.2 kcal mol−1, and RMSD vary from 11.1 to 21.3 kcal mol−1. The statistics for functionals utilizing the generalized gradient approximation are displayed in Figures S16 and S17.

Table 11. Best Performing Methods for the Calculation of ΔfH°gas, ΔrH°gas, and ΔH°gas (All Reactions) STOa

GTO ΔfH°gas B97-1/BS4 B97-1/BS3 B97-1/BS6 B97-1/BS7

ΔrH°gas B97-1/BS3 PBEPBE/BS4

all reactions B97-1/BS3 B97-1/BS4 B97-1/BS6

MAD (kcal mol−1) 8.3 9.2 9.4 9.7

MAD (kcal mol−1) 9.3 9.5

MAD (kcal mol−1) 9.1 9.8 9.8

ΔfH°gas B97-1/BS8 KCIS/BS8 revPBE/BS8 HCTH(407)/ BS8 BLAP3/BS8 ΔrH°gas VS98/BS8 TPSSh/BS8 TPSS-D/BS8 B97-1/BS8 KT2/BS8 B3LYP-D/BS8 M06/BS8 BP86-D/BS8 all reactions B97-1/BS8 VS98/BS8 KCIS/BS8

MAD (kcal mol−1) 5.5 8.0 8.5 9.0 9.4 MAD (kcal mol−1) 5.6 7.8 8.5 8.9 9.1 9.5 9.5 9.9 MAD (kcal mol−1) 6.9 9.0 9.4

a

All DFT/BS8(STO) computations utilized B3LYP/BS3 optimized geometries and temperature-dependent enthalpy corrections.

ΔrH°gas via Methods Employing GTO Basis Sets. ΔrH°gas was calculated via DFT methods utilizing GTO basis sets and ccCA for each reaction shown in Table 6. The differences between DFT- and ccCA-derived values are shown in Tables S19−S24. We find that the MSDs vary from 4.4 to 18.9 kcal mol−1, MADs vary from 9.3 to 41.1 kcal mol−1, and RMSDs vary from 11.7 to 61.3 kcal mol−1. These statistics are summarized as vertical bar charts for the hybrid functionals B3LYP, B97-1, and BHandH (Figure S12) and for the pure J

DOI: 10.1021/acs.jpca.6b09278 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

The best performing STO methods (MAD < 10 kcal mol−1) are summarized in Table 10 when all reactions are considered (ΔfH°gas and ΔrH°gas). The MADs of the GGA, meta-GGA, and hybrid functionals utilizing STO basis sets examined in this study are shown in Figure 4; BHandHLYP and tau-HCTHhybrid were excluded from these figures due to poor performance (MADs of 65.7 and 106.7 kcal mol−1). The B97-1/BS8//B3LYP/BS3 level of theory exhibited the lowest MAD (6.9 kcal mol−1) of all STO methods investigated herein, and interestingly, the other functionals in the B97 family exhibited MADs from 10.2 to 10.9 kcal mol−1. The VS98 and KCIS functionals, similar to that of B97-1, exhibited MADs less than 10 kcal mol−1 but larger than 6.9 kcal mol−1. Conclusions and Summary. The performances of 6 functionals utilizing 7 GTO basis sets and 69 functionals utilizing an STO basis set were evaluated based on calculation of the standard gas-phase reaction enthalpy for a variety of reactions at 298.15 K via comparison to ccCA-derived values. These reactions include formation, exchange, and ligand dissociation reactions of Cr-containing species with organic ligands. A summary of the best performing methods for the calculation of ΔfH°gas, ΔrH°gas, and overall thermochemistry is shown in Table 11. When considering GTO basis sets, the B97-1/BS4 level of theory (MAD = 8.3 kcal mol−1) performs best for the calculation of ΔfH°gas. However, BS4 employs the largest number of basis functions of the GTO basis sets examined. The B97-1 functional with BS3, BS6, and BS7 were also competitive, and the accuracy of the B97-1 functional increases directly as the number of basis functions increases. The STO basis set (BS8) is the largest basis set employed in this study, and the B97-1/BS8 level of theory (MAD = 5.5 kcal mol−1) is also quite reliable for determining ΔfH°gas of Cr-containing molecules. Several other functionals utilizing BS8 are competitive for the calculation of ΔfH°gas, including KCIS, revPBE, HCTH(407), and BLAP3. For calculation of ΔrH°gas, the B97-1/BS3 level of theory (MAD = 9.1 kcal mol−1) is promising, while PBEPBE/BS4 obtained similar accuracy. No other functional/basis set combination employing a GTO basis set for the determination of ΔrH°gas achieved a MAD of less than 10 kcal mol−1. It is interesting to note that PBEPBE obtained reasonable accuracy (MAD of 11.1 kcal mol−1) when paired with the smallest basis set (BS1); therefore, this level of theory may be acceptable if computational cost is a concern. In contrast, B97-1/BS1 is much less accurate (MAD of 18.1 kcal mol−1). Of the methods utilizing an STO basis set (BS8), the VS98 functional performed best (MAD of 5.6 kcal mol−1) for the calculation of ΔrH°gas. Several other functionals employing BS8 are competitive with MADs of less than 10 kcal mol−1. These functionals include TPSSh, TPSS-D, B97-1, KT2, B3LYP-D, M06, and BP86-D. When all reactions of Cr-containing compounds are considered, the B97-1/BS3 level of theory is optimal (MAD = 9.1 kcal mol−1) for GTO methods with a similar accuracy achieved by B97-1/BS4 and B97-1/BS6. No other level of theory utilizing a GTO basis set achieved a MAD of less than 10 kcal mol−1 overall. For STO methods employing BS8 (and B3LYP/BS3 optimized geometries), the B97-1 functional has the best accuracy (MAD = 6.9 kcal mol−1). The VS98 and KCIS functionals are competitive with MADs of 9.0 and 9.4 kcal mol−1, respectively. Our laboratories will be using the most accurate levels of theory shown here to explore ligand design

A similar comparison for functionals employing metageneralized gradient approximations and hybrid approximations was also performed. For DFT methods utilizing meta-GGA functionals and STO basis sets, we find that the MSDs vary from −0.2 to 17.0 kcal mol−1, MADs vary from 5.6 to 20.7 kcal mol−1, and RMSDs vary from 7.1 to 25.8 kcal mol−1. The statistics for meta-GGA functionals are summarized for each method as vertical bar charts in Figure S18. DFT methods utilizing hybrid functionals and STO basis sets, excluding τ-HCTH-hybrid/BS8//B3LYP/BS3, exhibited MSDs that vary from 2.0 to 16.6 kcal mol−1, MADs that vary from 5.5 to 32.3 kcal mol−1, and RMSDs that vary from 7.8 to 43.8 kcal mol−1. The statistics for hybrid functionals are summarized in Figures S19 and S20. A chart of the MADs for all DFT methods utilizing BS8 is provided for ease of comparison (Figure S21). ΔrH°gas computed via VS98/BS8//B3LYP/BS3 was the most accurate model chemistry of the DFT methods employing STO basis sets examined in this study and exhibited a MAD of 5.6 kcal mol−1 compared to ccCA-derived values. Seven other competitive model chemistries were identified with MSDs of less than 10.0 kcal mol−1. TPSSh/BS8//B3LYP/BS3, TPSS-D/ BS8//B3LYP/BS3, B97-1/BS8//B3LYP/BS3, KT2/BS8// B3LYP/BS3, B3LYP-D/BS8//B3LYP/BS3, M06/BS8// B3LYP/BS3, and BP86-D/BS8//B3LYP/BS3 displayed MADs of 7.8, 8.5, 8.9, 9.1, 9.5, 9.5, and 9.9 kcal mol−1, respectively. A plot of the ccCA-derived ΔrH°gas values versus VS98/BS8//B3LYP/BS3 is shown in Figure S22. Again, the R2 value of 0.994 indicates that these data sets are well-correlated. Overall Method Performance When All Reactions Are Considered. As discussed previously, the enthalpy of formation is a particularly difficult property to calculate accurately. However, ΔfH°gas is fundamentally a bond forming reaction. Therefore, we evaluate the performance of each method employed in this study by considering all reactions, including ΔfH°gas and ΔrH°gas, in a single data set. The best performing GTO methods (MAD < 10 kcal mol−1) are summarized in Table 5 when all reactions are considered. The MADs of the six functionals utilizing GTO basis sets examined in this study are shown in Figure 3. The B97-1/BS3 level of theory exhibited the lowest MAD (9.1 kcal mol−1) of all GTO methods investigated. The B97-1/BS4 and B97-1/BS6 levels of theory also exhibited similar MADs of 9.8 and 9.8 kcal mol−1, respectively. When averaging MADs of all basis sets from a specific functional, the average MADs of B97-1/BS1−BS7 and BVP86/BS1−BS7 are 15.2 and 16.9 kcal mol−1, respectively. The basis set type-averaged MADs of the other functionals examined with the suite of GTOs are all greater than 20 kcal mol−1. BS3 consistently achieved the best results for the hybrid functionals B3LYP, B97-1, and BHandH, while BS2 provided the best results for the pure functionals BVP86 and PBEPBE. By way of comparison, the MADs for B97-1/BS3 and BVP86/ BS2 are 9.8 and 12.1 kcal mol−1. While no direct correlation exists between basis set size and accuracy for a given functional, in general, hybrid functionals perform better with larger basis sets, and the pure functionals perform better with smaller basis sets. For example, a qualitative basis set dependency is shown for the best performing functional (B97-1) in Table 9, where BS3, BS4, and BS6 are more accurate than BS1 and BS2. As expected, the basis set employing a nonrelativistic effective core potential (BS5) consistently provided the worst results for any given functional, as shown in Figure 3. K

DOI: 10.1021/acs.jpca.6b09278 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

(7) To, T. T.; Heilweil, E. J.; Burkey, T. J. Time-Resolved Infrared Absorption Study of Cyclopentadienyl Manganese Tricarbonyl Derivatives: Chelation of Pendant Sulfides in Acetonitrile. J. Phys. Chem. A 2006, 110, 10669−10673. (8) To, T. T.; Heilweil, E. J.; Duke, C. B.; Ruddick, K. R.; Webster, C. E.; Burkey, T. J. Development of Ultrafast Photochromic Organometallics and Photoinduced Linkage Isomerization of Arene Chromium Carbonyl Derivatives. J. Phys. Chem. A 2009, 113, 2666− 2676. (9) Farrell, G. J.; Burkey, T. J. High-Pressure Photoacoustic Calorimetry of Chromium Hexacarbonyl: Volumes of Heptane Displacement from Chromium Pentacarbonyl Heptane. J. Photochem. Photobiol. A 2000, 137, 135−139. (10) Nayak, S. K.; Burkey, T. J. Photosubstitution of Iron Carbonyl Phosphine Complexes - Quantum Yield, Kinetic, and Thermochemical Studies. J. Am. Chem. Soc. 1993, 115, 6391−6397. (11) Garg, K.; Engle, J. T.; Ziegler, C. J.; Rack, J. J. Tuning Excited State Isomerization Dynamics through Ground State Structural Changes in Analogous Ruthenium and Osmium Sulfoxide Complexes. Chem. - Eur. J. 2013, 19, 11686−11695. (12) Garg, K.; King, A. W.; Rack, J. J. One Photon Yields Two Isomerizations: Large Atomic Displacements During Electronic Excited-State Dynamics in Ruthenium Sulfoxide Complexes. J. Am. Chem. Soc. 2014, 136, 1856−1863. (13) Garg, K.; Paris, S. I. M.; Rack, J. J. A Flexible Chelate Leads to Phototriggered Isomerization in an Osmium Sulfoxide Complex. Eur. J. Inorg. Chem. 2013, 2013, 1142−1148. (14) Jin, Y. H.; Rack, J. J. Molecules in Motion: From Sub-Nanoscale to Macroscale. Isr. J. Chem. 2013, 53, 280−287. (15) McClure, B. A.; Rack, J. J. Ultrafast Spectroscopy of a Photochromic Ruthenium Sulfoxide Complex. Inorg. Chem. 2011, 50, 7586−7590. (16) Porter, B. L.; McClure, B. A.; Abrams, E. R.; Engle, J. T.; Ziegler, C. J.; Rack, J. J. Photoisomerization in an Analogous Set of Ruthenium Sulfoxide Complexes. J. Photochem. Photobiol., A 2011, 217, 341−346. (17) Daffron, J. A.; Farrell, G. J.; Burkey, T. J. High-Pressure Photoacoustic Calorimetry. Rev. Sci. Instrum. 2000, 71, 3882−3885. (18) Peters, K. S.; Snyder, G. J. Time-Resolved Photoacoustic Calorimetry - Probing the Energetics and Dynamics of Fast Chemical and Biochemical Reactions. Science 1988, 241, 1053−1057. (19) Westrick, J. A.; Goodman, J. L.; Peters, K. S. A Time-Resolved Photoacoustic Calorimetry Study of the Dynamics of Enthalpy and Volume Changes Produced in the Photodissociation of CarbonMonoxide from Sperm Whale Carboxymyoglobin. Biochemistry 1987, 26, 8313−8318. (20) DeYonker, N. J.; Cundari, T. R.; Wilson, A. K. The Correlation Consistent Composite Approach (ccCA): An Alternative to the Gaussian-N Methods. J. Chem. Phys. 2006, 124, 114104. (21) DeYonker, N. J.; Grimes, T.; Yockel, S.; Dinescu, A.; Mintz, B.; Cundari, T. R.; Wilson, A. K. The Correlation-Consistent Composite Approach: Application to the G3/99 Test Set. J. Chem. Phys. 2006, 125, 104111. (22) DeYonker, N. J.; Ho, D. S.; Wilson, A. K.; Cundari, T. R. Computational s-Block Thermochemistry with the Correlation Consistent Composite Approach. J. Phys. Chem. A 2007, 111, 10776−10780. (23) DeYonker, N. J.; Wilson, B. R.; Pierpont, A. W.; Cundari, T. R.; Wilson, A. K. Towards the Intrinsic Error of the Correlation Consistent Composite Approach (ccCA). Mol. Phys. 2009, 107, 1107−1121. (24) Ho, D. S.; DeYonker, N. J.; Wilson, A. K.; Cundari, T. R. Accurate Enthalpies of Formation of Alkali and Alkaline Earth Metal Oxides and Hydroxides: Assessment of the Correlation Consistent Composite Approach (ccCA). J. Phys. Chem. A 2006, 110, 9767−9770. (25) Wilson, A. K.; DeYonker, N. J.; Cundari, T. R. The Correlation Consistent Composite Approach (ccCA): Efficient and Pan-Periodic Kinetics and Thermodynamics. Prog. Theor. Chem. Phys. 2009, 19, 197−224.

and catalysis of Cr-containing compounds. When these levels of theory are used to elucidate proposed mechanisms where the chemical environment is not changing drastically, cancellation of errors will improve the accuracy tremendously. The calibration of methodologies increases confidence in the results from investigations of novel transition-metal chemistry.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications Website at DOI: 10.1021/acs.jpca.6b09278. Supporting tabulated data and illustrative figures are provided in Tables S1−S35 and Figures S1−S22 (PDF)



AUTHOR INFORMATION

Corresponding Authors

*E-mail: [email protected]. *E-mail: [email protected]. ORCID

Charles Edwin Webster: 0000-0002-6917-2957 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work was partially supported by the National Science Foundation under grants CHE 0911528 and OIA 1539035 and the Mississippi State University Office of Research and Economic Development. Computational work was performed using resources at the University of Memphis High-Performance Computing Facility and the Computational Research on Materials Institute at the University of Memphis (CROMIUM).



REFERENCES

(1) Duke, C. B.; Letterman, R. G.; Johnson, J. O.; Barr, J. W.; Hu, S.; Ross, C. R.; Webster, C. E.; Burkey, T. J. Photochemistry of Chromium Arene Tricarbonyl Complexes with Tethered Pyridinyl and Propenyl Groups: Investigations of the Effect of Ring Size on Chelate Formation, Structure, and Linkage Isomerization. Organometallics 2014, 33, 485−497. (2) Gittermann, S. M.; Letterman, R. G.; Jiao, T.; Leu, G.-L.; DeYonker, N. J.; Webster, C. E.; Burkey, T. J. Bond Energies, Reaction Volumes, and Kinetics for σ- and π-Complexes of Mo(CO)5L. J. Phys. Chem. A 2011, 115, 9004−9013. (3) Jiao, T. J.; Leu, G. L.; Farrell, G. J.; Burkey, T. J. Photoacoustic Calorimetry and Quantum Yields of Mo(CO)6 Ligand Exchange in Linear Alkanes: Determination of Volume of Reaction, Energetics, and Kinetics of Nucleophile Displacement of Alkane from Mo(CO)5(Alkane). J. Am. Chem. Soc. 2001, 123, 4960−4965. (4) Jiao, T. J.; Pang, Z.; Burkey, T. J.; Johnston, R. F.; Heimer, T. A.; Kleiman, V. D.; Heilweil, E. J. Ultrafast Ring Closure Energetics and Dynamics of Cyclopentadienyl Manganese Tricarbonyl Derivatives. J. Am. Chem. Soc. 1999, 121, 4618−4624. (5) To, T. T.; Barnes, C. E.; Burkey, T. J. Bistable Photochromic Organometallics Based on Linkage Isomerization: Photochemistry of Dicarbonyl(η5-Methylcyclopentadienyl)Manganese(I) Derivatives with a Bifunctional, Nonchelating Ligand. Organometallics 2004, 23, 2708−2714. (6) To, T. T.; Duke, C. B.; Junker, C. S.; O’Brien, C. M.; Ross, C. R.; Barnes, C. E.; Webster, C. E.; Burkey, T. J. Linkage Isomerization as a Mechanism for Photochromic Materials: Cyclopentadienylmanganese Tricarbonyl Derivatives with Chelatable Functional Groups. Organometallics 2008, 27, 289−296. L

DOI: 10.1021/acs.jpca.6b09278 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A (26) DeYonker, N. J.; Mintz, B.; Cundari, T. R.; Wilson, A. K. Application of the Correlation Consistent Composite Approach (ccCA) to Third-Row (Ga-Kr) Molecules. J. Chem. Theory Comput. 2008, 4, 328−334. (27) Burgess, D. R. Thermochemical Data. In NIST Chemistry WebBook, NIST Standard Reference Database Number 69; Mallard, W. G., Linstrom, P. J., Eds.; NIST: Gaitherburg, MD, 2000. (28) DeYonker, N. J.; Williams, T. G.; Imel, A. E.; Cundari, T. R.; Wilson, A. K. Accurate Thermochemistry for Transition Metal Complexes from First-Principles Calculations. J. Chem. Phys. 2009, 131, 024106. (29) Jiang, W.; DeYonker, N. J.; Determan, J. J.; Wilson, A. K. Toward Accurate Theoretical Thermochemistry of First Row Transition Metal Complexes. J. Phys. Chem. A 2012, 116, 870−885. (30) DeYonker, N. J.; Peterson, K. A.; Steyl, G.; Wilson, A. K.; Cundari, T. R. Quantitative Computational Thermochemistry of Transition Metal Species. J. Phys. Chem. A 2007, 111, 11269−11277. (31) Li, S.; Dixon, D. A. Molecular Structures and Energetics of the (TiO2)n (n = 1−4) Clusters and Their Anions. J. Phys. Chem. A 2008, 112, 6646−6666. (32) Furche, F.; Perdew, J. P. The Performance of Semilocal and Hybrid Density Functionals in 3d Transition-Metal Chemistry. J. Chem. Phys. 2006, 124, 044103. (33) Uddin, J.; Morales, C. M.; Maynard, J. H.; Landis, C. R. Computational Studies of Metal−Ligand Bond Enthalpies across the Transition-Metal Series. Organometallics 2006, 25, 5566−5581. (34) Riley, K. E.; Merz, K. M., Jr. Assessment of Density Functional Theory Methods for the Computation of Heats of Formation and Ionization Potentials of Systems Containing Third Row Transition Metals. J. Phys. Chem. A 2007, 111, 6044−6053. (35) Tekarli, S. M.; Drummond, M. L.; Williams, T. G.; Cundari, T. R.; Wilson, A. K. Performance of Density Functional Theory for 3d Transition Metal-Containing Complexes: Utilization of the Correlation Consistent Basis Sets. J. Phys. Chem. A 2009, 113, 8607−8614. (36) Jiang, W. Y.; Laury, M. L.; Powell, M.; Wilson, A. K. Comparative Study of Single and Double Hybrid Density Functionals for the Prediction of 3d Transition Metal Thermochemistry. J. Chem. Theory Comput. 2012, 8, 4102−4111. (37) Harrison, J. F. Electronic Structure of Diatomic Molecules Composed of a First-Row Transition Metal and Main-Group Element (H-F). Chem. Rev. 2000, 100, 679−716. (38) Weymuth, T.; Couzijn, E. P. A.; Chen, P.; Reiher, M. New Benchmark Set of Transition-Metal Coordination Reactions for the Assessment of Density Functionals. J. Chem. Theory Comput. 2014, 10, 3092−3103. (39) Krieger, J. B.; Chen, J.; Iafrate, G. J.; Kurth, S. In Electron Correlations and Materials Properties; Gonis, A., Kioussis, N., Ciftan, M., Eds.; Kluwer Academic: New York, 1999; pp 463−477. (40) Qu, Z. W.; Hansen, A.; Grimme, S. Co-C Bond Dissociation Energies in Cobalamin Derivatives and Dispersion Effects: Anomaly or Just Challenging? J. Chem. Theory Comput. 2015, 11, 1037−1045. (41) Minenkov, Y.; Chermak, E.; Cavallo, L. Troubles in the Systematic Prediction of Transition Metal Thermochemistry with Contemporary Out-of-the-Box Methods. J. Chem. Theory Comput. 2016, 12, 1542−1560. (42) Chan, B.; Karton, A.; Raghavachari, K.; Radom, L. Heats of Formation for CrO, CrO2, and CrO3: An Extreme Challenge for Black-Box Composite Procedures. J. Chem. Theory Comput. 2012, 8, 3159−3166. (43) Zhang, W. J.; Truhlar, D. G.; Tang, M. S. Tests of ExchangeCorrelation Functional Approximations against Reliable Experimental Data for Average Bond Energies of 3d Transition Metal Compounds. J. Chem. Theory Comput. 2013, 9, 3965−3977. (44) Xu, X. F.; Zhang, W. J.; Tang, M. S.; Truhlar, D. G. Do Practical Standard Coupled Cluster Calculations Agree Better Than KohnSham Calculations with Currently Available Functionals When Compared to the Best Available Experimental Data for Dissociation Energies of Bonds to 3d Transition Metals? J. Chem. Theory Comput. 2015, 11, 2036−2052.

(45) Wheeler, S. E.; Houk, K. N.; Schleyer, P. V. R.; Allen, W. D. A Hierarchy of Homodesmotic Reactions for Thermochemistry. J. Am. Chem. Soc. 2009, 131, 2547−2560. (46) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; et al. Gaussian 09, revision B.01; Gaussian, Inc.: Wallingford, CT, 2009. (47) Becke, A. D. Density-Functional Thermochemistry 3. The Role of Exact Exchange. J. Chem. Phys. 1993, 98, 5648−5652. (48) Lee, C. T.; Yang, W. T.; Parr, R. G. Development of the ColleSalvetti Correlation-Energy Formula into a Functional of the ElectronDensity. Phys. Rev. B: Condens. Matter Mater. Phys. 1988, 37, 785−789. (49) Becke, A. D. Density-Functional Exchange-Energy Approximation with Correct Asymptotic-Behavior. Phys. Rev. A: At., Mol., Opt. Phys. 1988, 38, 3098−3100. (50) Becke, A. D. Density-Functional Thermochemistry. 5. Systematic Optimization of Exchange-Correlation Functionals. J. Chem. Phys. 1997, 107, 8554−8560. (51) Hamprecht, F. A.; Cohen, A. J.; Tozer, D. J.; Handy, N. C. Development and Assessment of New Exchange-Correlation Functionals. J. Chem. Phys. 1998, 109, 6264−6271. (52) Perdew, J. P. Density-Functional Approximation for the Correlation-Energy of the Inhomogeneous Electron-Gas. Phys. Rev. B: Condens. Matter Mater. Phys. 1986, 33, 8822−8824. (53) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (54) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1997, 78, 1396−1396. (55) Parr, R. G.; Yang, W. Density Functional Theory of Atoms and Molecules; Oxford University Press: New York, 1989. (56) Guerra, C. F.; Snijders, J. G.; te Velde, G.; Baerends, E. J. Towards an Order-N DFT Method. Theor. Chem. Acc. 1998, 99, 391− 403. (57) Bickelhaupt, F. M.; Velde, G. T.; Baerends, E. J.; Guerra, C. F.; Van Gisbergen, S. J. A.; Snijders, J. G.; Ziegler, T. Chemistry with ADF. J. Comput. Chem. 2001, 22, 931−967. (58) Baerends, E. J.; Ziegler, T.; Atkins, A. J.; Autschbach, J.; Bashford, D.; Bérces, A.; Bickelhaupt, F. M.; Bo, C.; Boerrigter, P. M.; Cavallo, L.; et al. ADF2008.01 and ADF2009.01; SCM, Theoretical Chemistry: Vrije Universiteit, Amsterdam, The Netherlands. (59) Vosko, S. H.; Wilk, L.; Nusair, M. Accurate Spin-Dependent Electron Liquid Correlation Energies for Local Spin-Density Calculations - a Critical Analysis. Can. J. Phys. 1980, 58, 1200−1211. (60) Perdew, J. P.; Kurth, S.; Zupan, A.; Blaha, P. Accurate Density Functional with Correct Formal Properties: A Step Beyond the Generalized Gradient Approximation. Phys. Rev. Lett. 1999, 82, 2544− 2547. (61) Chang, C.; Pelissier, M.; Durand, P. Regular 2-Component Pauli-Like Effective-Hamiltonians in Dirac Theory. Phys. Scr. 1986, 34, 394−404. (62) van Lenthe, E.; Baerends, E. J.; Snijders, J. G. Relativistic Regular 2-Component Hamiltonians. J. Chem. Phys. 1993, 99, 4597− 4610. (63) van Lenthe, E.; Baerends, E. J.; Snijders, J. G. Relativistic TotalEnergy Using Regular Approximations. J. Chem. Phys. 1994, 101, 9783−9792. (64) van Lenthe, E.; van Leeuwen, R.; Baerends, E. J.; Snijders, J. G. Relativistic Regular Two-Component Hamiltonians. Int. J. Quantum Chem. 1996, 57, 281−293. (65) van Lenthe, E.; Ehlers, A.; Baerends, E. J. Geometry Optimizations in the Zero Order Regular Approximation for Relativistic Effects. J. Chem. Phys. 1999, 110, 8943−8953. (66) van Lenthe, E.; Snijders, J. G.; Baerends, E. J. The Zero-Order Regular Approximation for Relativistic Effects: The Effect of SpinOrbit Coupling in Closed Shell Molecules. J. Chem. Phys. 1996, 105, 6505−6516. (67) Philipsen, P. H. T.; van Lenthe, E.; Snijders, J. G.; Baerends, E. J. Relativistic Calculations on the Adsorption of CO on the (111) Surfaces of Ni, Pd, and Pt within the Zeroth-Order Regular M

DOI: 10.1021/acs.jpca.6b09278 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A Approximation. Phys. Rev. B: Condens. Matter Mater. Phys. 1997, 56, 13556−13562. (68) Handy, N. C.; Cohen, A. J. Left-Right Correlation Energy. Mol. Phys. 2001, 99, 403−412. (69) Perdew, J. P. Correction. Phys. Rev. B: Condens. Matter Mater. Phys. 1986, 34, 7406−7406. (70) Watson, M. A.; Handy, N. C.; Cohen, A. J. Density Functional Calculations, Using Slater Basis Sets, with Exact Exchange. J. Chem. Phys. 2003, 119, 6475−6481. (71) Tao, J. M.; Perdew, J. P.; Staroverov, V. N.; Scuseria, G. E. Climbing the Density Functional Ladder: Nonempirical MetaGeneralized Gradient Approximation Designed for Molecules and Solids. Phys. Rev. Lett. 2003, 91, 1. (72) Staroverov, V. N.; Scuseria, G. E.; Tao, J. M.; Perdew, J. P. Comparative Assessment of a New Nonempirical Density Functional: Molecules and Hydrogen-Bonded Complexes. J. Chem. Phys. 2003, 119, 12129−12137. (73) Tsuneda, T.; Suzumura, T.; Hirao, K. A New One-Parameter Progressive Colle-Salvetti-Type Correlation Functional. J. Chem. Phys. 1999, 110, 10664−10678. (74) Krieger, J. B.; Chen, J.; Iafrate, G. J.; Kurth, S. Modified KriegerChen-Iafrate-Savin (KCIS) Functional. Unpublished work. (75) Zhao, Y.; Schultz, N. E.; Truhlar, D. G. Exchange-Correlation Functional with Broad Accuracy for Metallic and Nonmetallic Compounds, Kinetics, and Noncovalent Interactions. J. Chem. Phys. 2005, 123, 161103. (76) Toulouse, J.; Adamo, C. A New Hybrid Functional Including a Meta-GGA Approach. Chem. Phys. Lett. 2002, 362, 72−78. (77) Wilson, P. J.; Bradley, T. J.; Tozer, D. J. Hybrid ExchangeCorrelation Functional Determined from Thermochemical Data and Ab Initio Potentials. J. Chem. Phys. 2001, 115, 9233−9242. (78) Zhao, Y.; Schultz, N. E.; Truhlar, D. G. Design of Density Functionals by Combining the Method of Constraint Satisfaction with Parametrization for Thermochemistry, Thermochemical Kinetics, and Noncovalent Interactions. J. Chem. Theory Comput. 2006, 2, 364−382. (79) Grimme, S. Semiempirical GGA-Type Density Functional Constructed with a Long-Range Dispersion Correction. J. Comput. Chem. 2006, 27, 1787−1799. (80) Kurth, S.; Perdew, J. P.; Blaha, P. Molecular and Solid-State Tests of Density Functional Approximations: LSD, GGAs, and MetaGGAs. Int. J. Quantum Chem. 1999, 75, 889−909. (81) Zhao, Y.; Truhlar, D. G. The M06 Suite of Density Functionals for Main Group Thermochemistry, Thermochemical Kinetics, Noncovalent Interactions, Excited States, and Transition Elements: Two New Functionals and Systematic Testing of Four M06-Class Functionals and 12 Other Functionals. Theor. Chem. Acc. 2008, 120, 215−241. (82) Perdew, J. P.; Kurth, S.; Zupan, A.; Blaha, P. Erratum: Accurate Density Functional with Correct Formal Properties: A Step Beyond the Generalized Gradient Approximation [Phys. Rev. Lett. 82, 2544 (1999)]. Phys. Rev. Lett. 1999, 82, 5179−5179. (83) Ernzerhof, M.; Scuseria, G. E. Assessment of the Perdew-BurkeErnzerhof Exchange-Correlation Functional. J. Chem. Phys. 1999, 110, 5029−5036. (84) Adamo, C.; Barone, V. Toward Reliable Density Functional Methods without Adjustable Parameters: The PBE0 Model. J. Chem. Phys. 1999, 110, 6158−6170. (85) Perdew, J. P.; Ernzerhof, M.; Burke, K. Rationale for Mixing Exact Exchange with Density Functional Approximations. J. Chem. Phys. 1996, 105, 9982−9985. (86) Swart, M.; Ehlers, A. W.; Lammertsma, K. Performance of the OPBE Exchange-Correlation Functional. Mol. Phys. 2004, 102, 2467− 2474. (87) OPBE mixed with 25% exact exchange. (88) Xu, X.; Goddard, W. A. The X3LYP Extended Density Functional for Accurate Descriptions of Nonbond Interactions, Spin States, and Thermochemical Properties. Proc. Natl. Acad. Sci. U. S. A. 2004, 101, 2673−2677.

(89) Boese, A. D.; Handy, N. C. New Exchange-Correlation Density Functionals: The Role of the Kinetic-Energy Density. J. Chem. Phys. 2002, 116, 9559−9569. (90) Proynov, E.; Chermette, H.; Salahub, D. R. New TauDependent Correlation Functional Combined with a Modified Becke Exchange. J. Chem. Phys. 2000, 113, 10013−10027. (91) Adamo, C.; Barone, V. Exchange Functionals with Improved Long-Range Behavior and Adiabatic Connection Methods without Adjustable Parameters: The MPW and MPW1PW Models. J. Chem. Phys. 1998, 108, 664−675. (92) Zhao, Y.; Truhlar, D. G. A New Local Density Functional for Main-Group Thermochemistry, Transition Metal Bonding, Thermochemical Kinetics, and Noncovalent Interactions. J. Chem. Phys. 2006, 125, 194101. (93) Lynch, B. J.; Fast, P. L.; Harris, M.; Truhlar, D. G. Adiabatic Connection for Kinetics. J. Phys. Chem. A 2000, 104, 4811−4815. (94) Proynov, E. I.; Sirois, S.; Salahub, D. R. Extension of the LAP Functional to Include Parallel Spin Correlation. Int. J. Quantum Chem. 1997, 64, 427−446. (95) Adamo, C.; Barone, V. Toward Reliable Adiabatic Connection Models Free from Adjustable Parameters. Chem. Phys. Lett. 1997, 274, 242−250. (96) Hertwig, R. H.; Koch, W. On the Parameterization of the Local Correlation Functional. What Is Becke-3-LYP? Chem. Phys. Lett. 1997, 268, 345−351. (97) Boese, A. D.; Doltsinis, N. L.; Handy, N. C.; Sprik, M. New Generalized Gradient Approximation Functionals. J. Chem. Phys. 2000, 112, 1670−1678. (98) Van Voorhis, T.; Scuseria, G. E. A Novel Form for the Exchange-Correlation Energy Functional. J. Chem. Phys. 1998, 109, 400−410. (99) Boese, A. D.; Handy, N. C. A New Parametrization of ExchangeCorrelation Generalized Gradient Approximation Functionals. J. Chem. Phys. 2001, 114, 5497−5503. (100) Becke, A. D. Simulation of Delocalized Exchange by Local Density Functionals. J. Chem. Phys. 2000, 112, 4020−4026. (101) Stephens, P. J.; Devlin, F. J.; Chabalowski, C. F.; Frisch, M. J. Ab initio Calculation of Vibrational Absorption and CircularDichroism Spectra Using Density-Functional Force-Fields. J. Phys. Chem. 1994, 98, 11623−11627. (102) Reiher, M.; Salomon, O.; Hess, B. A. Reparameterization of Hybrid Functionals Based on Energy Differences of States of Different Multiplicity. Theor. Chem. Acc. 2001, 107, 48−55. (103) Cohen, A. J.; Handy, N. C. Dynamic Correlation. Mol. Phys. 2001, 99, 607−615. (104) Adamo, C.; Barone, V. Physically Motivated Density Functionals with Improved Performances: The Modified PerdewBurke-Ernzerhof Model. J. Chem. Phys. 2002, 116, 5933−5940. (105) Zhang, Y. K.; Yang, W. T. Comment on ″Generalized Gradient Approximation Made Simple″. Phys. Rev. Lett. 1998, 80, 890−890. (106) Kang, J. K.; Musgrave, C. B. Prediction of Transition State Barriers and Enthalpies of Reaction by a New Hybrid DensityFunctional Approximation. J. Chem. Phys. 2001, 115, 11040−11051. (107) Hammer, B.; Hansen, L. B.; Norskov, J. K. Improved Adsorption Energetics within Density-Functional Theory Using Revised Perdew-Burke-Ernzerhof Functionals. Phys. Rev. B: Condens. Matter Mater. Phys. 1999, 59, 7413−7421. (108) 50% exact exchange and 50% LDA exchange plus LYP correlation (109) 50% exact exchange and 50% Becke88 exchange plus LYP correlation (110) Perdew, J. P. Electronic Structure of Solids; Akademie Verlag: Berlin, 1991; p 11−20. (111) Perdew, J. P.; Chevary, J. A.; Vosko, S. H.; Jackson, K. A.; Pederson, M. R.; Singh, D. J.; Fiolhais, C. Atoms, Molecules, Solids, and Surfaces - Applications of the Generalized Gradient Approximation for Exchange and Correlation. Phys. Rev. B: Condens. Matter Mater. Phys. 1992, 46, 6671−6687. N

DOI: 10.1021/acs.jpca.6b09278 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A

(132) Roy, L. E.; Hay, P. J.; Martin, R. L. Revised Basis Sets for the LANL Effective Core Potentials. J. Chem. Theory Comput. 2008, 4, 1029−1031. (133) Hay, P. J.; Wadt, W. R. Ab Initio Effective Core Potentials for Molecular Calculations - Potentials for K to Au Including the Outermost Core Orbitals. J. Chem. Phys. 1985, 82, 299−310. (134) Krishnan, R.; Binkley, J. S.; Seeger, R.; Pople, J. A. SelfConsistent Molecular-Orbital Methods. 20. Basis Set for Correlated Wave-Functions. J. Chem. Phys. 1980, 72, 650−654. (135) McLean, A. D.; Chandler, G. S. Contracted Gaussian-Basis Sets for Molecular Calculations. 1. 2nd Row Atoms, Z = 11−18. J. Chem. Phys. 1980, 72, 5639−5648. (136) Dolg, M.; Wedig, U.; Stoll, H.; Preuss, H. Energy-Adjusted Abinitio Pseudopotentials for the 1st-Row Transition-Elements. J. Chem. Phys. 1987, 86, 866−872. (137) Martin, J. M. L.; Sundermann, A. Correlation Consistent Valence Basis Sets for Use with the Stuttgart-Dresden-Bonn Relativistic Effective Core Potentials: The Atoms Ga-Kr and In-Xe. J. Chem. Phys. 2001, 114, 3408−3420. (138) Hay, P. J.; Wadt, W. R. Ab Initio Effective Core Potentials for Molecular Calculations - Potentials for the Transition-Metal Atoms Sc to Hg. J. Chem. Phys. 1985, 82, 270−283. (139) Wadt, W. R.; Hay, P. J. Ab Initio Effective Core Potentials for Molecular Calculations - Potentials for Main Group Elements Na to Bi. J. Chem. Phys. 1985, 82, 284−298. (140) Dunlap, B. I. Fitting the Coulomb Potential Variationally in XAlpha Molecular Calculations. J. Chem. Phys. 1983, 78, 3140−3142. (141) Dunlap, B. I. Robust and Variational Fitting: Removing the Four-Center Integrals from Center Stage in Quantum Chemistry. J. Mol. Struct.: THEOCHEM 2000, 529, 37−40. (142) Dunlap, B. I.; Connolly, J. W. D.; Sabin, J. R. Some Approximations in Applications of X-Alpha Theory. J. Chem. Phys. 1979, 71, 3396−3402. (143) Dunlap, B. I.; Connolly, J. W. D.; Sabin, J. R. 1st-Row Diatomic-Molecules and Local Density Models. J. Chem. Phys. 1979, 71, 4993−4999. (144) Truncation of the density-fitting set was carried out at the default value for angular momentum, which is defined in Gaussian 09 as Max(MaxTyp+1,2*MaxVal), where MaxTyp is the highest angular momentum in the AO basis set and MaxVal is the highest valence angular momentum. (145) Curtiss, L. A.; Raghavachari, K.; Redfern, P. C.; Pople, J. A. Assessment of Gaussian-2 and Density Functional Theories for the Computation of Enthalpies of Formation. J. Chem. Phys. 1997, 106, 1063−1079. (146) NIST. NIST Chemistry WebBook; National Institute of Standards and Technology: Washington, D.C. (accessed May 14, 2016). (147) The experimental standard enthalpy of formation for Cr(CO)6 maintained by the NIST WebBook (−218 ± 20 kcal/mol) results from the average of nine values. The ccCA-derived value in this work compares most favorably to the most stable of the nine experimental values. (148) Ruscic, B.; Pinzon, R. E.; von Laszewski, G.; Kodeboyina, D.; Burcat, A.; Leahy, D.; Montoy, D.; Wagner, A. F. Active Thermochemical Tables: Thermochemistry for the 21st Century. J. Phys.: Conf. Ser. 2005, 16, 561−570. (149) Ruscic, B. Active Thermochemical Tables (ATcT); Argonne National Laboratory: Lemont, IL, 2015. (150) Ruscic, B.; Pinzon, R. E.; Morton, M. L.; von Laszevski, G.; Bittner, S. J.; Nijsure, S. G.; Amin, K. A.; Minkoff, M.; Wagner, A. F. Introduction to Active Thermochemical Tables: Several ″Key″ Enthalpies of Formation Revisited. J. Phys. Chem. A 2004, 108, 9979−9997. (151) Yungman, V. S. Thermal Constants of Substances; Wiley, Begell House: New York, 1999. (152) Ebbinghaus, B. B. Thermodynamics of Gas-Phase Chromium Species - the Chromium Chlorides, Oxychlorides, Fluorides, Oxyfluorides, Hydroxides, Oxyhydroxides, Mixed Oxyfluorochlorohydr-

(112) Perdew, J. P.; Chevary, J. A.; Vosko, S. H.; Jackson, K. A.; Pederson, M. R.; Singh, D. J.; Fiolhais, C. Erratum: Atoms, Molecules, Solids, and Surfaces - Applications of the Generalized Gradient Approximation for Exchange and Correlation [Phys. Rev. B 46, 6671 (1992)]. Phys. Rev. B: Condens. Matter Mater. Phys. 1993, 48, 4978. (113) Filatov, M.; Thiel, W. A New Gradient-Corrected ExchangeCorrelation Density Functional. Mol. Phys. 1997, 91, 847−859. (114) Keal, T. W.; Tozer, D. J. The Exchange-Correlation Potential in Kohn-Sham Nuclear Magnetic Resonance Shielding Calculations. J. Chem. Phys. 2003, 119, 3015−3024. (115) Becke, A. D.; Roussel, M. R. Exchange Holes in Inhomogeneous Systems - a Coordinate-Space Model. Phys. Rev. A: At., Mol., Opt. Phys. 1989, 39, 3761−3767. (116) Becke, A. D. Correlation-Energy of an Inhomogeneous Electron-Gas - a Coordinate-Space Model. J. Chem. Phys. 1988, 88, 1053−1062. (117) Becke, A. D. Thermochemical Tests of a Kinetic-Energy Dependent Exchange-Correlation Approximation. Int. J. Quantum Chem. 1994, 52, 625−632. (118) Davidson, E. R. Comment on Dunning’s CorrelationConsistent Basis Sets. Chem. Phys. Lett. 1996, 260, 514−518. (119) Couty, M.; Hall, M. B. Basis Sets for Transition Metals: Optimized Outer P Functions. J. Comput. Chem. 1996, 17, 1359−1370. (120) The 6-31G(d') basis sets for carbon, nitrogen, oxygen, and fluorine have the d polarization functions taken from the 6-311G(d) basis sets. Foresman, J. B.; Frisch, Æ. Exploring Chemistry with Electronic Structure Methods, 2nd ed.; Gaussian, Inc.: Pittsburgh, PA, p 110. For other elements, the 6-31G(d') keyword in Gaussian 09 calls the 6-31G(d) basis sets. (121) Petersson, G. A.; Allaham, M. A. A Complete Basis Set Model Chemistry. Open-Shell Systems and the Total Energies of the 1st-Row Atoms. J. Chem. Phys. 1991, 94, 6081−6090. (122) Petersson, G. A.; Bennett, A.; Tensfeldt, T. G.; Allaham, M. A.; Shirley, W. A.; Mantzaris, J. A Complete Basis Set Model Chemistry. 1. The Total Energies of Closed-Shell Atoms and Hydrides of the 1stRow Elements. J. Chem. Phys. 1988, 89, 2193−2218. (123) Clark, T.; Chandrasekhar, J.; Spitznagel, G. W.; Schleyer, P. V. Efficient Diffuse Function-Augmented Basis-Sets for Anion Calculations 3. The 3-21+G Basis-Set for 1st-Row Elements, Li-F. J. Comput. Chem. 1983, 4, 294−301. (124) Dunning, T. H. Gaussian-Basis Sets for Use in Correlated Molecular Calculations. 1. The Atoms Boron through Neon and Hydrogen. J. Chem. Phys. 1989, 90, 1007−1023. (125) Woon, D. E.; Dunning, T. H. Gaussian-Basis Sets for Use in Correlated Molecular Calculations. 4. Calculation of Static Electrical Response Properties. J. Chem. Phys. 1994, 100, 2975−2988. (126) Woon, D. E.; Dunning, T. H. Gaussian-Basis Sets for Use in Correlated Molecular Calculations. 5. Core-Valence Basis-Sets for Boron through Neon. J. Chem. Phys. 1995, 103, 4572−4585. (127) Woon, D. E.; Dunning, T. H. Gaussian-Basis Sets for Use in Correlated Molecular Calculations. 3. The Atoms Aluminum through Argon. J. Chem. Phys. 1993, 98, 1358−1371. (128) Wilson, A. K.; Woon, D. E.; Peterson, K. A.; Dunning, T. H. Gaussian Basis Sets for Use in Correlated Molecular Calculations. IX. The Atoms Gallium through Krypton. J. Chem. Phys. 1999, 110, 7667− 7676. (129) Balabanov, N. B.; Peterson, K. A. Systematically Convergent Basis Sets for Transition Metals. I. All-Electron Correlation Consistent Basis Sets for the 3d Elements Sc-Zn. J. Chem. Phys. 2005, 123, 064107. (130) Balabanov, N. B.; Peterson, K. A. Basis Set Limit Electronic Excitation Energies, Ionization Potentials, and Electron Affinities for the 3d Transition Metal Atoms: Coupled Cluster and Multireference Methods. J. Chem. Phys. 2006, 125, 074110. (131) Ehlers, A. W.; Böhme, M.; Dapprich, S.; Gobbi, A.; Höllwarth, A.; Jonas, V.; Köhler, K. F.; Stegmann, R.; Veldkamp, A.; Frenking, G. Set of f-Polarization Functions for Pseudo-Potential Basis-Sets of the Transition-Metals Sc-Cu, Y-Ag and La-Au. Chem. Phys. Lett. 1993, 208, 111−114. O

DOI: 10.1021/acs.jpca.6b09278 J. Phys. Chem. A XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry A oxides, and Volatility Calculations in Waste Incineration Processes. Combust. Flame 1995, 101, 311−338. (153) Ebbinghaus, B. B. Thermodynamics of Gas-Phase Chromium Species - the Chromium Oxides, the Chromium Oxyhydroxides, and Volatility Calculations in Waste Incineration Processes. Combust. Flame 1993, 93, 119−137. (154) Gritsenko, O. V.; Schipper, P. R. T.; Baerends, E. J. Exchange and Correlation Energy in Density Functional Theory: Comparison of Accurate Density Functional Theory Quantities with Traditional Hartree-Fock Based Ones and Generalized Gradient Approximations for the Molecules Li2, N2, F2. J. Chem. Phys. 1997, 107, 5007−5015. (155) Fouqueau, A.; Casida, M. E.; Daku, L. M. L.; Hauser, A.; Neese, F. Comparison of Density Functionals for Energy and Structural Differences between the High- [5T2g:(t2g)4(eg)2] and Low[1A1g:(t2g)6(eg)0] Spin States of Iron(II) Coordination Compounds. II. More Functionals and the Hexaminoferrous Cation, [Fe(NH3)6]2+. J. Chem. Phys. 2005, 122, 044110. (156) Lewis, K. E.; Golden, D. M.; Smith, G. P. Organometallic Bond-Dissociation Energies - Laser Pyrolysis of Fe(CO)5, Cr(CO)6, Mo(CO)6, and W(CO)6. J. Am. Chem. Soc. 1984, 106, 3905−3912. (157) Fletcher, T. R.; Rosenfeld, R. N. Recombination of Cr(CO)n with CO - Kinetics and Bond-Dissociation Energies. J. Am. Chem. Soc. 1988, 110, 2097−2101. (158) Schlappi, D. N.; Cedeno, D. L. Metal-Olefin Bond Energies in M(CO)5(C2H(4-n)Cln) M = Cr, Mo, W; n = 0−4: ElectronWithdrawing Olefins Do Not Increase the Bond Strength. J. Phys. Chem. A 2009, 113, 9692−9699. (159) Cedeno, D. L.; Weitz, E. Experimental Determination of the Cr-C2Cl4 Bond Dissociation Enthalpy in Cr(CO)5(C2Cl4): Quantifying Metal-Olefin Bonding Interactions. J. Am. Chem. Soc. 2001, 123, 12857−12865. (160) McNamara, B.; Becher, D. M.; Towns, M. H.; Grant, E. R. Elementary Arrhenius Parameters in the CO-for-Ethylene Dissociative Substitution of Cr(CO)5(C2H4). J. Phys. Chem. 1994, 98, 4622−4626. (161) Wells, J. R.; House, P. G.; Weitz, E. Interaction of H2 and Prototypical Solvent Molecules with Cr(CO)5 in the Gas-Phase. J. Phys. Chem. 1994, 98, 8343−8351.

P

DOI: 10.1021/acs.jpca.6b09278 J. Phys. Chem. A XXXX, XXX, XXX−XXX