1
|
de Moura CEV, Sokolov AY. Efficient Spin-Adapted Implementation of Multireference Algebraic Diagrammatic Construction Theory. I. Core-Ionized States and X-ray Photoelectron Spectra. J Phys Chem A 2024; 128:5816-5831. [PMID: 38962857 DOI: 10.1021/acs.jpca.4c03161] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/05/2024]
Abstract
We present an efficient implementation of multireference algebraic diagrammatic construction theory (MR-ADC) for simulating core-ionized states and X-ray photoelectron spectra (XPS). Taking advantage of spin adaptation, automatic code generation, and density fitting, our implementation can perform calculations for molecules with more than 1500 molecular orbitals, incorporating static and dynamic correlation in the ground and excited electronic states. We demonstrate the capabilities of MR-ADC methods by simulating the XPS spectra of substituted ferrocene complexes and azobenzene isomers. For the ground electronic states of these molecules, the XPS spectra computed using the extended second-order MR-ADC method (MR-ADC(2)-X) are in a very good agreement with available experimental results. We further show that MR-ADC can be used as a tool for interpreting or predicting the results of time-resolved XPS measurements by simulating the core ionization spectra of azobenzene along its photoisomerization, including the XPS signatures of excited states and the minimum energy conical intersection. This work is the first in a series of publications reporting the efficient implementations of MR-ADC methods.
Collapse
Affiliation(s)
- Carlos E V de Moura
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Alexander Yu Sokolov
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
2
|
Vysotskiy VP, Filippi C, Ryde U. Scalar Relativistic All-Electron and Pseudopotential Ab Initio Study of a Minimal Nitrogenase [Fe(SH) 4H] - Model Employing Coupled-Cluster and Auxiliary-Field Quantum Monte Carlo Many-Body Methods. J Phys Chem A 2024; 128:1358-1374. [PMID: 38324717 PMCID: PMC10895656 DOI: 10.1021/acs.jpca.3c05808] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2023] [Revised: 01/03/2024] [Accepted: 01/03/2024] [Indexed: 02/09/2024]
Abstract
Nitrogenase is the only enzyme that can cleave the triple bond in N2, making nitrogen available to organisms. The detailed mechanism of this enzyme is currently not known, and computational studies are complicated by the fact that different density functional theory (DFT) methods give very different energetic results for calculations involving nitrogenase models. Recently, we designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S to mimic different possible protonation states of the nitrogenase active site. We showed that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical methods. Based on nonrelativistic single-reference coupled-cluster (CC) calculations, we estimated that the ΔE is 101 kJ/mol. In this study, we demonstrate that scalar relativistic effects play an important role and significantly affect ΔE. Our best revised single-reference CC estimates for ΔE are 85-91 kJ/mol, including energy corrections to account for contributions beyond triples, core-valence correlation, and basis-set incompleteness error. Among coupled-cluster approaches with approximate triples, the canonical CCSD(T) exhibits the largest error for this problem. Complementary to CC, we also used phaseless auxiliary-field quantum Monte Carlo calculations (ph-AFQMC). We show that with a Hartree-Fock (HF) trial wave function, ph-AFQMC reproduces the CC results within 5 ± 1 kJ/mol. With multi-Slater-determinant (MSD) trials, the results are 82-84 ± 2 kJ/mol, indicating that multireference effects may be rather modest. Among the DFT methods tested, τ-HCTH, r2SCAN with 10-13% HF exchange with and without dispersion, and O3LYP/O3LYP-D4, and B3LYP*/B3LYP*-D4 generally perform the best. The r2SCAN12 (with 12% HF exchange) functional mimics both the best reference MSD ph-AFQMC and CC ΔE results within 2 kJ/mol.
Collapse
Affiliation(s)
- Victor P. Vysotskiy
- Department
of Computational Chemistry, Lund University,
Chemical Centre, SE-221 00 Lund, Sweden
| | - Claudia Filippi
- MESA+
Institute for Nanotechnology, University
of Twente, P.O. Box 217, Enschede 7500 AE, Netherlands
| | - Ulf Ryde
- Department
of Computational Chemistry, Lund University,
Chemical Centre, SE-221 00 Lund, Sweden
| |
Collapse
|
3
|
Boychuk BTA, Wetmore SD. Assessment of Density Functional Theory Methods for the Structural Prediction of Transition and Post-Transition Metal-Nucleic Acid Complexes. J Chem Theory Comput 2023. [PMID: 37399186 DOI: 10.1021/acs.jctc.3c00127] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/05/2023]
Abstract
Understanding the structure of metal-nucleic acid systems is important for many applications such as the design of new pharmaceuticals, metal detection platforms, and nanomaterials. Herein, we explore the ability of 20 density functional theory (DFT) functionals to reproduce the crystal structure geometry of transition and post-transition metal-nucleic acid complexes identified in the Protein Data Bank and Cambridge Structural Database. The environmental extremes of the gas phase and implicit water were considered, and analysis focused on the global and inner coordination geometry, including the coordination distances. Although gas-phase calculations were unable to describe the structure of 12 out of the 53 complexes in our test set regardless of the DFT functional considered, accounting for the broader environment through implicit solvation or constraining the model truncation points to crystallographic coordinates generally afforded agreement with the experimental structure, suggesting that functional performance for these systems is likely due to the models rather than the methods. For the remaining 41 complexes, our results show that the reliability of functionals depends on the metal identity, with the magnitude of error varying across the periodic table. Furthermore, minimal changes in the geometries of these metal-nucleic acid complexes occur upon use of the Stuttgart-Dresden effective core potential and/or inclusion of an implicit water environment. The overall top three performing functionals are ωB97X-V, ωB97X-D3(BJ), and MN15, which reliably describe the structure of a broad range of metal-nucleic acid systems. Other suitable functionals include MN15-L, which is a cheaper alternative to MN15, and PBEh-3c, which is commonly used in QM/MM calculations of biomolecules. In fact, these five methods were the only functionals tested to reproduce the coordination sphere of Cu2+-containing complexes. For metal-nucleic acid systems that do not contain Cu2+, ωB97X and ωB97X-D are also suitable choices. These top-performing methods can be utilized in future investigations of diverse metal-nucleic acid complexes of relevance to biology and material science.
Collapse
Affiliation(s)
- Briana T A Boychuk
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, AB T1K 3M4, Canada
| | - Stacey D Wetmore
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, AB T1K 3M4, Canada
| |
Collapse
|
4
|
S Almeida NM, Melin TRL, North SC, Welch BK, Wilson AK. Ab initio composite strategies and multireference approaches for lanthanide sulfides and selenides. J Chem Phys 2022; 157:024105. [PMID: 35840393 DOI: 10.1063/5.0094367] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The f-block ab initio correlation consistent composite approach was used to predict the dissociation energies of lanthanide sulfides and selenides. Geometry optimizations were carried out using density functional theory and coupled cluster singles, doubles, and perturbative triples with one- and two-component Hamiltonians. For the two-component calculations, relativistic effects were accounted for by utilizing a third-order Douglas-Kroll-Hess Hamiltonian. Spin-orbit coupling was addressed with the Breit-Pauli Hamiltonian within a multireference configuration interaction approach. The state averaged complete active space self-consistent field wavefunctions obtained for the spin-orbit coupling energies were used to assign the ground states of diatomics, and several diagnostics were used to ascertain the multireference character of the molecules.
Collapse
Affiliation(s)
- Nuno M S Almeida
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48864, USA
| | - Timothé R L Melin
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48864, USA
| | - Sasha C North
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48864, USA
| | - Bradley K Welch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48864, USA
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48864, USA
| |
Collapse
|
5
|
Otlyotov AA, Minenkov Y, Zaitsau DH, Zherikova KV, Verevkin SP. "In Vitro" and "In Vivo" Diagnostic Check for the Thermochemistry of Metal-Organic Compounds. Inorg Chem 2022; 61:10743-10755. [PMID: 35797430 DOI: 10.1021/acs.inorgchem.2c00959] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Volatile metal β-diketonates are of interest from both practical and theoretical perspectives (manufacturing of film materials, catalysis, and the nature of metal-ligand bonding). Knowledge of their reliable thermochemical properties is essential for effective applications. However, there is an unacceptable scattering of the available data on the enthalpies of formation. In this work, we proposed "in vitro" and "in vivo" diagnostic tools to verify the available enthalpies of formation in both the crystalline and gaseous states for metal tris-β-diketonates. The "in vitro" procedure involved high-level quantum-chemical calculations and was applied to define a consistent data set on the enthalpies of formation for iron(III) β-diketonates. This data set has provided the basis for "in vivo" structure-property-based diagnostics to evaluate the robustness of the thermochemical data for β-diketonate tris-complexes with metals other than iron.
Collapse
Affiliation(s)
- Arseniy A Otlyotov
- N. N. Semenov Federal Research Center for Chemical Physics RAS, Moscow 119991, Russian Federation
| | - Yury Minenkov
- N. N. Semenov Federal Research Center for Chemical Physics RAS, Moscow 119991, Russian Federation.,Joint Institute for High Temperatures, Russian Academy of Sciences, Moscow 125412, Russian Federation
| | - Dzmitry H Zaitsau
- Department of Physical Chemistry and Faculty of Interdisciplinary Research, Competence Centre CALOR, University of Rostock, 18059 Rostock, Germany
| | - Kseniya V Zherikova
- Nikolaev Institute of Inorganic Chemistry of Siberian Branch of Russian Academy of Sciences, Novosibirsk 630090, Russian Federation
| | - Sergey P Verevkin
- Department of Physical Chemistry and Faculty of Interdisciplinary Research, Competence Centre CALOR, University of Rostock, 18059 Rostock, Germany.,Department of Physical Chemistry, Kazan Federal University, Kazan 420008, Russian Federation
| |
Collapse
|
6
|
North SC, Wilson AK. Ab Initio Composite Approaches for Heavy Element Energetics: Ionization Potentials for the Actinide Series of Elements. J Phys Chem A 2022; 126:3027-3042. [PMID: 35427146 DOI: 10.1021/acs.jpca.2c01007] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The first, second, and third gas-phase ionization potentials have been determined for the actinide series of elements using an ab initio composite scalar and fully relativistic approach, employing the coupled cluster with single, double, and perturbative triple excitations (CCSD(T)) and Dirac Hartree-Fock (DHF) methods, extrapolated to the complete basis set (CBS) limit. The impact of electron correlation and basis set choice within this framework are examined. Additionally, the first three ionization potentials were obtained using an ab initio heavy element correlation-consistent Composite Approach (here referred to as α-ccCA). This is the first utilization of a ccCA for actinide species.
Collapse
Affiliation(s)
- Sasha C North
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
7
|
Pham BQ, Datta D, Gordon MS. PDG: A Composite Method Based on the Resolution of the Identity. J Phys Chem A 2021; 125:9421-9429. [PMID: 34658243 DOI: 10.1021/acs.jpca.1c06186] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The Gaussian-3 (G3) composite approach for thermochemical properties is revisited in light of the enhanced computational efficiency and reduced memory costs by applying the resolution-of-the-identity (RI) approximation for two-electron repulsion integrals (ERIs) to the computationally demanding component methods in the G3 model: the energy and gradient computations via the second-order Møller-Plesset perturbation theory (MP2) and the energy computations using the coupled-cluster singles-doubles method augmented with noniterative triples corrections [CCSD(T)]. Efficient implementation of the RI-based methods is achieved by employing a hybrid distributed/shared memory model based on MPI and OpenMP. The new variant of the G3 composite approach based on the RI approximation is termed the RI-G3 scheme, or alternatively the PDG method. The accuracy of the new RI-G3/PDG scheme is compared to the "standard" G3 composite approach that employs the memory-expensive four-center ERIs in the MP2 and CCSD(T) calculations. Taking the computation of the heats of formation of the closed-shell molecules in the G3/99 test set as a test case, it is demonstrated that the RI approximation introduces negligible changes to the mean absolute errors relative to the standard G3 model (less than 0.1 kcal/mol), while the standard deviations remain unaltered. The efficiency and memory requirements for the RI-MP2 and RI-CCSD(T) methods are compared to the standard MP2 and CCSD(T) approaches, respectively. The hybrid MPI/OpenMP-based RI-MP2 energy plus gradient computation is found to attain a 7.5× speedup over the standard MP2 calculations. For the most demanding CCSD(T) calculations, the application of the RI approximation is found to nearly halve the memory demand, confer about a 4-5× speedup for the CCSD iterations, and reduce the computational time for the compute-intensive triples correction step by several hours.
Collapse
|
8
|
Maurer LR, Bursch M, Grimme S, Hansen A. Assessing Density Functional Theory for Chemically Relevant Open-Shell Transition Metal Reactions. J Chem Theory Comput 2021; 17:6134-6151. [PMID: 34546754 DOI: 10.1021/acs.jctc.1c00659] [Citation(s) in RCA: 55] [Impact Index Per Article: 18.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Due to the principle lack of systematic improvement possibilities of density functional theory, careful assessment of the performance of density functional approximations (DFAs) on well-designed benchmark sets, for example, for reaction energies and barrier heights, is crucial. While main-group chemistry is well covered by several available sets, benchmark data for transition metal chemistry is sparse. This is especially the case for larger, chemically relevant molecules. Addressing this issue, we recently introduced the MOR41 benchmark which covers chemically relevant reactions of closed-shell complexes. In this work, we extend these efforts to single-reference open-shell systems and introduce the "reactions of open-shell single-reference transition metal complexes" (ROST61) benchmark set. ROST61 includes accurate coupled-cluster reference values for 61 reaction energies with a mean reaction energy of -42.8 kcal mol-1. Complexes with 13-93 atoms covering 20 d-block elements are included, but due to the restriction to single-reference open-shell systems, important elements such as iron or platinum could not be taken into account, or only to a small extent. We assess the performance of 31 DFAs in combination with three London dispersion (LD) correction schemes. Further, DFT-based composite methods, MP2, and a few semiempirical quantum chemical methods are evaluated. Consistent with the results for the MOR41 closed-shell benchmark, we find that the ordering of DFAs according to Jacob's ladder is preserved and that adding an LD correction is crucial, clearly improving almost all tested methods. The recently introduced r2SCAN-3c composite method stands out with a remarkable mean absolute deviation (MAD) of only 2.9 kcal mol-1, which is surpassed only by hybrid DFAs with low amounts of Fock exchange (e.g., 2.3 kcal mol-1 for TPSS0-D4/def2-QZVPP) and double-hybrid (DH) DFAs but at a significantly higher computational cost. The lowest MAD of only 1.6 kcal mol-1 is obtained with the DH DFA PWPB95-D4 in the def2-QZVPP basis set approaching the estimated accuracy of the reference method. Overall, the ROST61 set adds important reference data to a sparsely sampled but practically relevant area of chemistry. At this point, it provides valuable orientation for the application and development of new DFAs and electronic structure methods in general.
Collapse
Affiliation(s)
- Leonard R Maurer
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Markus Bursch
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| |
Collapse
|
9
|
Olive LN, Dornshuld EV, Webster CE. The curious case of DMSO: A CCSD(T)/CBS(aQ56+d) benchmark and DFT study. J Chem Phys 2021; 155:114304. [PMID: 34551533 DOI: 10.1063/5.0057327] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
This work addresses the pathological behavior of the energetics of dimethyl sulfoxide and related sulfur-containing compounds by providing the computational benchmark energetics of R2E2 species, where R = H/CH3 and E = O/S, with bent and pyramidal geometries using state-of-the-art methodologies. These 22 geometries were fully characterized with coupled-cluster with single, double, and perturbative triple excitations [CCSD(T)], second-order Møller-Plesset perturbation theory (MP2), and 22 density functional theory (DFT) methods with 8, 12, and 12, respectively, correlation consistent basis sets of double-, triple-, or quadruple-ζ quality. The relative energetics were determined at the MP2 and CCSD(T) complete basis set (CBS) limits using 17 basis sets up to sextuple-ζ and include augmented, tight-d, and core-valence correlation consistent basis sets. The relative energies of oxygen-/sulfur-containing compounds exhibit exceptionally slow convergence to the CBS limit with canonical methods as well as significant basis set dependence. CCSD(T) with quadruple-ζ basis sets can give qualitatively incorrect relative energies. Explicitly correlated MP2-F12 and CCSD(T)-F12 methods dramatically accelerate the convergence of the relative energies to the CBS limit for these problematic compounds. The F12 methods with a triple-ζ quality basis set give relative energies that deviate no more than 0.41 kcal mol-1 from the benchmark CBS limit. The correlation consistent Composite Approach (ccCA), ccCA-TM (TM for transition metals), and G3B3 deviated by no more than 2 kcal mol-1 from the benchmark CBS limits. Relative energies for oxygen-/sulfur-containing systems fully characterized with DFT are quite unreliable even with triple-ζ quality basis sets, and 13 out of 45 combinations fortuitously give a relative energy that is within 1 kcal mol-1 on average from the benchmark CCSD(T) CBS limit for these systems.
Collapse
Affiliation(s)
- Laura N Olive
- Department of Chemistry, Mississippi State University, Mississippi State, Mississippi 39762-9573, USA
| | - Eric Van Dornshuld
- Department of Chemistry, Mississippi State University, Mississippi State, Mississippi 39762-9573, USA
| | - Charles Edwin Webster
- Department of Chemistry, Mississippi State University, Mississippi State, Mississippi 39762-9573, USA
| |
Collapse
|
10
|
Lam CS, Lau KC. High-Level ab Initio Predictions for the Ionization Energies, Bond Dissociation Energies, and Heats of Formation of Vanadium Methylidene, Vanadium Methyl Species, and Their Cations (VCH 2/VCH 2+, VCH 3/VCH 3+). J Phys Chem A 2021; 125:4957-4966. [PMID: 34076442 DOI: 10.1021/acs.jpca.1c01381] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The ionization energies of VCH2 and VCH3, the various 0 K bond dissociation energies (D0s) in their neutrals and cations, and their respective heats of formation at 0 and 298 K are computed by the single-reference, wave function-based CCSDTQ/CBS procedure. The core of the composite method is the approximation to the complete basis set (CBS) limit at the coupled cluster (CC) level which includes up to full quadruple excitations. The zero-point vibrational energy, core-valence correlation, spin-orbit coupling, and scalar relativistic effects have their contributions incorporated in an additive manner. For the species in the current study, this protocol requires geometry optimizations and harmonic frequency calculations practically no higher than the CCSD(T)/aug-cc-pwCVTZ and CCSD(T)/aug-cc-pVTZ levels, respectively. The present calculations successfully predict D0(V+-CH3) = 2.126 eV and D0(V+-CH2) = 3.298 eV in remarkable agreement with the data recently measured by a spin-orbit state selected V+ + CH4 collision experiment (Phys. Chem. Chem. Phys. 2021, 23, 273-286). The good accord encourages the use of CCSDTQ/CBS protocol in thermochemical predictions of various feasible product channels identified in methane activation by transition metal species.
Collapse
Affiliation(s)
- Chow-Shing Lam
- Department of Chemistry, City University of Hong Kong, Tat Chee Avenue, Kowloon, Hong Kong
| | - Kai-Chung Lau
- Department of Chemistry, City University of Hong Kong, Tat Chee Avenue, Kowloon, Hong Kong
| |
Collapse
|
11
|
Frey NC, Dornshuld EV, Webster CE. Benchmarking the Fluxional Processes of Organometallic Piano-Stool Complexes. Molecules 2021; 26:2310. [PMID: 33923446 PMCID: PMC8073612 DOI: 10.3390/molecules26082310] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2021] [Revised: 04/09/2021] [Accepted: 04/13/2021] [Indexed: 11/16/2022] Open
Abstract
The correlation consistent Composite Approach for transition metals (ccCA-TM) and density functional theory (DFT) computations have been applied to investigate the fluxional mechanisms of cyclooctatetraene tricarbonyl chromium ((COT)Cr(CO)3) and 1,3,5,7-tetramethylcyclooctatetraene tricarbonyl chromium, molybdenum, and tungsten ((TMCOT)M(CO)3 (M = Cr, Mo, and W)) complexes. The geometries of (COT)Cr(CO)3 were fully characterized with the PBEPBE, PBE0, B3LYP, and B97-1 functionals with various basis set/ECP combinations, while all investigated (TMCOT)M(CO)3 complexes were fully characterized with the PBEPBE, PBE0, and B3LYP methods. The energetics of the fluxional dynamics of (COT)Cr(CO)3 were examined using the correlation consistent Composite Approach for transition metals (ccCA-TM) to provide reliable energy benchmarks for corresponding DFT results. The PBE0/BS1 results are in semiquantitative agreement with the ccCA-TM results. Various transition states were identified for the fluxional processes of (COT)Cr(CO)3. The PBEPBE/BS1 energetics indicate that the 1,2-shift is the lowest energy fluxional process, while the B3LYP/BS1 energetics (where BS1 = H, C, O: 6-31G(d'); M: mod-LANL2DZ(f)-ECP) indicate the 1,3-shift having a lower electronic energy of activation than the 1,2-shift by 2.9 kcal mol-1. Notably, PBE0/BS1 describes the (CO)3 rotation to be the lowest energy process, followed by the 1,3-shift. Six transition states have been identified in the fluxional processes of each of the (TMCOT)M(CO)3 complexes (except for (TMCOT)W(CO)3), two of which are 1,2-shift transition states. The lowest-energy fluxional process of each (TMCOT)M(CO)3 complex (computed with the PBE0 functional) has a ΔG‡ of 12.6, 12.8, and 13.2 kcal mol-1 for Cr, Mo, and W complexes, respectively. Good agreement was observed between the experimental and computed 1H-NMR and 13C-NMR chemical shifts for (TMCOT)Cr(CO)3 and (TMCOT)Mo(CO)3 at three different temperature regimes, with coalescence of chemically equivalent groups at higher temperatures.
Collapse
Affiliation(s)
| | | | - Charles Edwin Webster
- Department of Chemistry, Mississippi State University, 310 President’s Circle, Starkville, MS 39762-9573, USA; (N.C.F.); (E.V.D.)
| |
Collapse
|
12
|
Das SK, Chakraborty S, Ramakrishnan R. Critical benchmarking of popular composite thermochemistry models and density functional approximations on a probabilistically pruned benchmark dataset of formation enthalpies. J Chem Phys 2021; 154:044113. [PMID: 33514111 DOI: 10.1063/5.0032713] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
First-principles calculation of the standard formation enthalpy, ΔHf° (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 ΔHf° 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.
Collapse
Affiliation(s)
- Sambit Kumar Das
- Tata Institute of Fundamental Research, Centre for Interdisciplinary Sciences, Hyderabad 500107, India
| | - Sabyasachi Chakraborty
- Tata Institute of Fundamental Research, Centre for Interdisciplinary Sciences, Hyderabad 500107, India
| | - Raghunathan Ramakrishnan
- Tata Institute of Fundamental Research, Centre for Interdisciplinary Sciences, Hyderabad 500107, India
| |
Collapse
|
13
|
Moulder CA, Kafle K, Zhou CX, Cundari TR. Thermochemistry of Tungsten-3p Elements for Density Functional Theory, Caveat Lector! J Phys Chem A 2021; 125:681-690. [PMID: 33405918 DOI: 10.1021/acs.jpca.0c05351] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
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, R2 = 0.98), and DLPNO-CCSD(T)/def2-QZVPP calculations for metal-ligand BDEs (R2 = 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).
Collapse
Affiliation(s)
- Catherine A Moulder
- Department of Chemistry, Center for Advanced Scientific Computing and Modeling (CASCaM), University of North Texas, 1155 Union Circle, #305070, Denton, Texas 76203-5017, United States
| | - Kristina Kafle
- Department of Chemistry, Center for Advanced Scientific Computing and Modeling (CASCaM), University of North Texas, 1155 Union Circle, #305070, Denton, Texas 76203-5017, United States
| | - Christopher X Zhou
- Department of Chemistry, Center for Advanced Scientific Computing and Modeling (CASCaM), University of North Texas, 1155 Union Circle, #305070, Denton, Texas 76203-5017, United States
| | - Thomas R Cundari
- Department of Chemistry, Center for Advanced Scientific Computing and Modeling (CASCaM), University of North Texas, 1155 Union Circle, #305070, Denton, Texas 76203-5017, United States
| |
Collapse
|
14
|
Fortenberry RC, DeYonker NJ. Rovibrational Quantum Chemical Treatment of Inorganic and Organometallic Astrochemicals. Acc Chem Res 2021; 54:271-279. [PMID: 33356121 DOI: 10.1021/acs.accounts.0c00631] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
ConspectusOur two groups have both independently and collaboratively been pushing quantum-chemical techniques to produce highly accurate predictions of anharmonic vibrational frequencies and spectroscopic constants for molecules containing atoms outside of the typical upper p block. Methodologies employ composite approaches, relying on various levels of coupled cluster theory-most often at the singles, doubles, and perturbative triples level-and quartic force field constructions of the potential portion of the intramolecular Watson Hamiltonian. Such methods are known to perform well for organic species, and we have extended this to molecules containing atoms outside of this realm.One notable atom that has received much attention in this application is magnesium. Mg is the second-most-abundant element in the Earth's mantle, and while molecules containing this element are among the confirmed astrochemicals, its further atomic abundance in the galaxy implies that many more molecules (both purely inorganic and organometallic) containing element 12 exist in astrophysical regions in chemical sizes between those of atoms and dust-sized nanocrystals. Our approach discussed herein is producing quality benchmarks and predicting novel data for magnesium-bearing molecules.The story is similar for Al and Si, which are also notably abundant in both rocky bodies and the universe at large. While Na, Sc, and Cu may not be as abundant as Mg, Al, and Si, molecules containing Na and transition metals have also previously been reported to be detected beyond the Earth. Consequently, the need to produce spectral reference data for molecules containing such atoms is growing. While several experimental groups (including, notably, the groups in Arizona, Boston, and France/Spain) have clearly led the way in detection of inorganic/organometallic molecules in space, computational support and even rational design can provide novel avenues for the detection of molecules containing atoms not typically studied in most laboratories. The application of quantum chemistry to other elements beyond carbon and its cronies at the top right of the periodic table promises a better understanding of the observable universe. It will also provide novel and fundamental chemical insights pushing the "central science" into new molecular territory.
Collapse
Affiliation(s)
- Ryan C. Fortenberry
- Department of Chemistry and Biochemistry, University of Mississippi, University, Mississippi 38677-1848, United States
| | - Nathan J. DeYonker
- Department of Chemistry, University of Memphis, Memphis, Tennessee 38152, United States
| |
Collapse
|
15
|
Probing the electronic structures and properties of neutral and charged FeSin(−1,0,+1) (n = 1–6) clusters using ccCA theory. J Mol Model 2020; 26:283. [DOI: 10.1007/s00894-020-04551-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2019] [Accepted: 09/17/2020] [Indexed: 10/23/2022]
|
16
|
Affiliation(s)
- Karl K. Irikura
- Chemical Sciences Division, National Institute of Standards and Technology, Gaithersburg, Maryland 20899-8320, United States
| |
Collapse
|
17
|
Burton MA, Cheng Q, Halfen DT, Lane JH, DeYonker NJ, Ziurys LM. The structure of ScC 2 (X̃ 2A 1): A combined Fourier transform microwave/millimeter-wave spectroscopic and computational study. J Chem Phys 2020; 153:034304. [PMID: 32716169 DOI: 10.1063/5.0008746] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023] Open
Abstract
Pure rotational spectra of Sc13C2 (X̃2A1) and Sc12C13C (X̃2A') 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 13CH4 or 13CH4/12CH4, diluted in argon. The NKa,Kc = 10,1 → 00,0, 20,2 → 10,1, 30,3 → 20,2, and 40,4 → 30,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 13C (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 Sc13C2 and Sc12C13C. In addition, a quartic force field was calculated for ScC2 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 Sc12C2, an accurate structure for ScC2 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 ScC2 contains a C2 - moiety and is largely ionic.
Collapse
Affiliation(s)
- M A Burton
- Department of Chemistry and Biochemistry, Department of Astronomy, Steward Observatory, University of Arizona, 1305 E. 4th Street, Tucson, Arizona 85719, USA
| | - Q Cheng
- Department of Chemistry, University of Memphis, Memphis, Tennessee 38152, USA
| | - D T Halfen
- Department of Chemistry and Biochemistry, Department of Astronomy, Steward Observatory, University of Arizona, 1305 E. 4th Street, Tucson, Arizona 85719, USA
| | - J H Lane
- Department of Chemistry and Biochemistry, Department of Astronomy, Steward Observatory, University of Arizona, 1305 E. 4th Street, Tucson, Arizona 85719, USA
| | - N J DeYonker
- Department of Chemistry, University of Memphis, Memphis, Tennessee 38152, USA
| | - L M Ziurys
- Department of Chemistry and Biochemistry, Department of Astronomy, Steward Observatory, University of Arizona, 1305 E. 4th Street, Tucson, Arizona 85719, USA
| |
Collapse
|
18
|
Moulder CA, Kafle K, Cundari TR. Tungsten-Ligand Bond Strengths for 2p Elements Including σ- and π-Bond Strength Components, A Density Functional Theory and ab Initio Study. J Phys Chem A 2019; 123:7940-7949. [PMID: 31240921 DOI: 10.1021/acs.jpca.9b03272] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Three WVI crystal structures with multifarious metal-ligand bond types are used to theoretically predict homolytic metal-element bond enthalpies with 11 popular DFT functionals, MP2 wave function methods, and four common valence basis set/pseudopotentials in order to evaluate the accuracy and precision of the resultant bond enthalpy data. To our knowledge, for the first time, estimates of component metal-ligand σ- and π-bond strengths are computed. The WE (E = C, N, O) bond enthalpies have the consistent trend σ > second π > first π. In contrast, the element-element BDE trend for the 2p homologues is second π > first π > σ for nitrogen and oxygen, and σ > first π > second π for carbon. These differences may underpin the differences in stability trends and thus reactivity behavior for metal-element multiple bonds as compared to the element-element multiple bonds, and metal-element triple bonds versus their corresponding double bonded counterparts. For example, Odom et al. show that MeI nucleophilically attacks at the imide (M═N) rather than the nitride (M ≡ N) ligand; the relative π-bond strengths derived herein provide a thermodynamic rationalization for this site preference. In this study, it is deduced from the calculated thermodynamics that the W-oxo ligand is more congruous with a triple bond than a double bond, consistent with the bonding model set forth in the seminal 1961 Ballhausen-Gray paper.
Collapse
Affiliation(s)
- Catherine A Moulder
- Department of Chemistry & Center for Advanced Scientific Computing and Modeling (CASCaM) , University of North Texas , 1155 Union Circle, #305070 , Denton , Texas 76203-5017 , United States
| | - Kristina Kafle
- Department of Chemistry & Center for Advanced Scientific Computing and Modeling (CASCaM) , University of North Texas , 1155 Union Circle, #305070 , Denton , Texas 76203-5017 , United States
| | - Thomas R Cundari
- Department of Chemistry & Center for Advanced Scientific Computing and Modeling (CASCaM) , University of North Texas , 1155 Union Circle, #305070 , Denton , Texas 76203-5017 , United States
| |
Collapse
|
19
|
Lam CS, Lau KC, Ng CY. High-Level Ab Initio Predictions for the Ionization Energy, Bond Dissociation Energies, and Heats of Formation of Vanadium Methylidyne Radical and Its Cation (VCH/VCH +). J Phys Chem A 2019; 123:7454-7462. [PMID: 31414807 DOI: 10.1021/acs.jpca.9b05493] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The ionization energy (IE) of VCH, the 0 K V-CH/VC-H bond dissociation energies (D0s), and the heats of formation at 0 K (ΔHf0°) and 298 K (ΔHf298°) for VCH/VCH+ are predicted by the wave function-based CCSDTQ/CBS approach. This composite-coupled cluster method includes full quadruple excitations in conjunction with the approximation to the complete basis set (CBS) limit. The contributions of zero-point vibrational energy, core-valence (CV) correlation, spin-orbit coupling, and scalar relativistic corrections are taken into account. The present calculations show that adiabatic IE(VCH) = 6.785 eV and demonstrate excellent agreement with an IE value of 6.774 7 ± 0.000 1 eV measured with two-color laser-pulsed field ionization-photoelectron spectroscopy. The CCSDT and MRCI+Q methods which include CV correlations give the best predictions of harmonic frequencies: ω2 (ω2+) (bending) = 534 (650) and 564 (641) cm-1 and the V-CH stretching ω3 (ω3+) = 835 (827) and 856 (857) cm-1 compared with the experimental values. In this work, we offer a streamlined CCSDTQ/CBS approach which shows an error limit (≤20 meV) matching with previous benchmarking efforts for reliable IE and D0 predictions for VCH/VCH+. The CCSDTQ/CBS D0(V+-CH) - D0(V-CH) = -0.012 eV and D0(VC+-H) - D0(VC-H) = 0.345 eV are in good accord with the experimentally derived values of -0.028 4 ± 0.000 1 and 0.355 9 ± 0.000 1 eV, respectively. The present study has demonstrated that the CCSDTQ/CBS protocol can be readily extended to investigate triatomic molecules containing 3d-metals.
Collapse
Affiliation(s)
- Chow-Shing Lam
- Department of Chemistry , City University of Hong Kong , Tat Chee Avenue , Kowloon , Hong Kong
| | - Kai-Chung Lau
- Department of Chemistry , City University of Hong Kong , Tat Chee Avenue , Kowloon , Hong Kong
| | - Cheuk-Yiu Ng
- Department of Chemistry , University of California, Davis , Davis , California 95616 , United States
| |
Collapse
|
20
|
Abstract
A fundamental need in chemistry is understanding the chemical bond, for which the most quantitative measure is the bond dissociation energy (BDE). While BDEs of chemical bonds formed from the lighter main group elements are generally well-known and readily calculated by modern computational chemistry, chemical bonds involving the transition metals, lanthanides, and actinides remain computationally extremely challenging. This is due to the simultaneous importance of electron correlation, spin-orbit interaction, and other relativistic effects, coupled with the large numbers of low-lying states that are accessible in systems with open d or f subshells. The development of efficient and accurate computational methods for these species is currently a major focus of the field. An obstacle to this effort has been the scarcity of highly precise benchmarks for the BDEs of M-X bonds. For most of the transition metal, lanthanide, or actinide systems, tabulated BDEs of M-X bonds have been determined by Knudsen effusion mass spectrometric measurements of high-temperature equilibria. The measured ion signals are converted to pressures and activities of the species involved in the equilibrium, and the equilibrium constants are then analyzed using a van't Hoff plot or the third-law method to extract the reaction enthalpy, which is extrapolated to 0 K to obtain the BDE. This procedure introduces errors at every step and ultimately leads to BDEs that are typically uncertain by 2-20 kcal mol-1 (0.1-1 eV). A second method in common use employs a thermochemical cycle in which the ionization energies of the MX molecule and M atom are combined with the BDE of the M+-X bond, obtained via guided ion beam mass spectrometry, to yield the BDE of the neutral, M-X. When accurate values of all three components of the cycle are available, this method yields good results-but only rarely are all three values available. We have recently implemented a new method for the precise measurement of BDEs in molecules with large densities of electronic states that is based on the rapid predissociation of these species as soon as the ground separated atom limit is exceeded. When a sharp predissociation threshold is observed, its value directly provides the BDE of the system. With this method, we are able in favorable cases to determine M-X BDEs to an accuracy of ∼0.1 kcal mol-1 (0.004 eV). The method is generally applicable to species that have a high density of states at the ground separated atom limit and has been used to measure the BDEs of more than 50 transition metal-main group MX molecules thus far. In addition, a number of metal-metal BDEs have also been measured with this method. There are good prospects for the extension of the method to polyatomic systems and to lanthanide and actinide-containing molecules. These precise BDE measurements provide chemical trends for the BDEs across the transition metal series, as well as crucial benchmarks for the development of efficient and accurate computational methods for the d- and f-block elements.
Collapse
Affiliation(s)
- Michael D. Morse
- Department of Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112, United States
| |
Collapse
|
21
|
Janesko BG, Proynov E, Scalmani G, Frisch MJ. Long-range-corrected Rung 3.5 density functional approximations. J Chem Phys 2018; 148:104112. [DOI: 10.1063/1.5017981] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Benjamin G. Janesko
- Department of Chemistry and Biochemistry, Texas Christian University, Fort Worth, Texas 76110, USA
| | - Emil Proynov
- Department of Chemistry and Biochemistry, Texas Christian University, Fort Worth, Texas 76110, USA
| | | | | |
Collapse
|
22
|
Minenkov Y, Bistoni G, Riplinger C, Auer AA, Neese F, Cavallo L. Pair natural orbital and canonical coupled cluster reaction enthalpies involving light to heavy alkali and alkaline earth metals: the importance of sub-valence correlation. Phys Chem Chem Phys 2018; 19:9374-9391. [PMID: 28327742 DOI: 10.1039/c7cp00836h] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
In this work, we tested canonical and domain based pair natural orbital coupled cluster methods (CCSD(T) and DLPNO-CCSD(T), respectively) for a set of 32 ligand exchange and association/dissociation reaction enthalpies involving ionic complexes of Li, Be, Na, Mg, Ca, Sr, Ba and Pb(ii). Two strategies were investigated: in the former, only valence electrons were included in the correlation treatment, giving rise to the computationally very efficient FC (frozen core) approach; in the latter, all non-ECP electrons were included in the correlation treatment, giving rise to the AE (all electron) approach. Apart from reactions involving Li and Be, the FC approach resulted in non-homogeneous performance. The FC approach leads to very small errors (<2 kcal mol-1) for some reactions of Na, Mg, Ca, Sr, Ba and Pb, while for a few reactions of Ca and Ba deviations up to 40 kcal mol-1 have been obtained. Large errors are both due to artificial mixing of the core (sub-valence) orbitals of metals and the valence orbitals of oxygen and halogens in the molecular orbitals treated as core, and due to neglecting core-core and core-valence correlation effects. These large errors are reduced to a few kcal mol-1 if the AE approach is used or the sub-valence orbitals of metals are included in the correlation treatment. On the technical side, the CCSD(T) and DLPNO-CCSD(T) results differ by a fraction of kcal mol-1, indicating the latter method as the perfect choice when the CPU efficiency is essential. For completely black-box applications, as requested in catalysis or thermochemical calculations, we recommend the DLPNO-CCSD(T) method with all electrons that are not covered by effective core potentials included in the correlation treatment and correlation-consistent polarized core valence basis sets of cc-pwCVQZ(-PP) quality.
Collapse
Affiliation(s)
- Yury Minenkov
- Physical Sciences and Engineering Division, King Abdullah University of Science and Technology, KAUST Catalysis Center, Thuwal 23955-6900, Saudi Arabia.
| | - Giovanni Bistoni
- Department of Molecular Theory and Spectroscopy, Max Planck Institute for Chemical Energy Conversion, Stiftstr. 34-36, D-45470 Mülheim an der Ruhr, Germany.
| | - Christoph Riplinger
- Department of Molecular Theory and Spectroscopy, Max Planck Institute for Chemical Energy Conversion, Stiftstr. 34-36, D-45470 Mülheim an der Ruhr, Germany.
| | - Alexander A Auer
- Department of Molecular Theory and Spectroscopy, Max Planck Institute for Chemical Energy Conversion, Stiftstr. 34-36, D-45470 Mülheim an der Ruhr, Germany.
| | - Frank Neese
- Department of Molecular Theory and Spectroscopy, Max Planck Institute for Chemical Energy Conversion, Stiftstr. 34-36, D-45470 Mülheim an der Ruhr, Germany.
| | - Luigi Cavallo
- Physical Sciences and Engineering Division, King Abdullah University of Science and Technology, KAUST Catalysis Center, Thuwal 23955-6900, Saudi Arabia.
| |
Collapse
|
23
|
Cheng Q, Fortenberry RC, DeYonker NJ. Towards a quantum chemical protocol for the prediction of rovibrational spectroscopic data for transition metal molecules: Exploration of CuCN, CuOH, and CuCCH. J Chem Phys 2017; 147:234303. [PMID: 29272934 DOI: 10.1063/1.5006931] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Qianyi Cheng
- Department of Chemistry, University of Memphis, Memphis, Tennessee 38152, USA
| | - Ryan C Fortenberry
- Department of Chemistry and Biochemistry, Georgia Southern University, Statesboro, Georgia 30460, USA
| | - Nathan J DeYonker
- Department of Chemistry, University of Memphis, Memphis, Tennessee 38152, USA
| |
Collapse
|
24
|
Plascencia C, Wang J, Wilson AK. Importance of the ligand basis set in ab initio thermochemical calculations of transition metal species. Chem Phys Lett 2017. [DOI: 10.1016/j.cplett.2017.08.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
25
|
Determan JJ, Poole K, Scalmani G, Frisch MJ, Janesko BG, Wilson AK. Comparative Study of Nonhybrid Density Functional Approximations for the Prediction of 3d Transition Metal Thermochemistry. J Chem Theory Comput 2017; 13:4907-4913. [PMID: 28877436 DOI: 10.1021/acs.jctc.7b00809] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The utility of several nonhybrid density functional approximations (DFAs) is considered for the prediction of gas phase enthalpies of formation for a large set of 3d transition metal-containing molecules. Nonhybrid DFAs can model thermochemical values for 3d transition metal-containing molecules with accuracy comparable to that of hybrid functionals. The GAM-generalized gradient approximation (GGA); the TPSS, M06-L, and MN15-L meta-GGAs; and the Rung 3.5 PBE+ΠLDA(s) DFAs all give root-mean-square deviations below that of the widely used B3LYP hybrid. Modern nonhybrid DFAs continue to show utility for transition metal thermochemistry.
Collapse
Affiliation(s)
- John J Determan
- Department of Chemistry, Texas Christian University , Fort Worth, Texas 76129, United States
| | - Katelyn Poole
- Department of Chemistry, Texas Christian University , Fort Worth, Texas 76129, United States.,Department of Chemistry and Center for Advanced Scientific Computing and Modeling (CASCaM), University of North Texas , Denton, Texas 76203-5017, United States
| | - Giovanni Scalmani
- Gaussian, Inc. , 340 Quinnipiac Street, Building 40, Wallingford, Connecticut 06492, United States
| | - Michael J Frisch
- Gaussian, Inc. , 340 Quinnipiac Street, Building 40, Wallingford, Connecticut 06492, United States
| | - Benjamin G Janesko
- Department of Chemistry, Texas Christian University , Fort Worth, Texas 76129, United States
| | - Angela K Wilson
- Department of Chemistry, Michigan State University , East Lansing, Michigan 48824 United States
| |
Collapse
|
26
|
A CASSCF/CASPT2 investigation on electron detachments from ScSi
n
−
(n = 4–6) clusters. J Mol Model 2017; 23:282. [DOI: 10.1007/s00894-017-3461-2] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2017] [Accepted: 09/05/2017] [Indexed: 10/18/2022]
|
27
|
Abstract
The exact exchange energy and its energy density are useful but computationally expensive ingredients in density functional approximations for Kohn-Sham density functional theory. We present detailed tests of some exact nonempirical upper bounds to exact exchange. These "Rung 3.5" upper bounds contract the Kohn-Sham one-particle density matrix with model density matrices used to construct semilocal model exchange holes and invoke the Cauchy-Schwarz inequality. The contraction automatically eliminates the computationally expensive long-range component of the exact exchange hole. Numerical tests show that the exchange upper bounds underestimate total exchange energies while predicting other properties with accuracy approaching standard hybrid approximations.
Collapse
Affiliation(s)
- Emil Proynov
- Department of Chemistry and Biochemistry, Texas Christian University , Fort Worth, Texas 76110, United States
| | - Benjamin G Janesko
- Department of Chemistry and Biochemistry, Texas Christian University , Fort Worth, Texas 76110, United States
| |
Collapse
|
28
|
Matthew DJ, Tieu E, Morse MD. Determination of the bond dissociation energies of FeX and NiX (X = C, S, Se). J Chem Phys 2017; 146:144310. [DOI: 10.1063/1.4979679] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Daniel J. Matthew
- Department of Chemistry, University of Utah, Salt Lake City, Utah UT 84112, USA
| | - Erick Tieu
- Department of Chemistry, University of Utah, Salt Lake City, Utah UT 84112, USA
| | - Michael D. Morse
- Department of Chemistry, University of Utah, Salt Lake City, Utah UT 84112, USA
| |
Collapse
|
29
|
|
30
|
Abstract
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.
Collapse
Affiliation(s)
| | - Kenneth M. Merz
- Department of Chemistry, Department of Biochemistry and Molecular Biology, and Institute of Cyber-Enabled Research, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
31
|
Fang Z, Vasiliu M, Peterson KA, Dixon DA. Prediction of Bond Dissociation Energies/Heats of Formation for Diatomic Transition Metal Compounds: CCSD(T) Works. J Chem Theory Comput 2017; 13:1057-1066. [PMID: 28080051 DOI: 10.1021/acs.jctc.6b00971] [Citation(s) in RCA: 75] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
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 T1 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.
Collapse
Affiliation(s)
- Zongtang Fang
- Department of Chemistry, The University of Alabama , Shelby Hall, Tuscaloosa, Alabama 35487-0336, United States
| | - Monica Vasiliu
- Department of Chemistry, The University of Alabama , Shelby Hall, Tuscaloosa, Alabama 35487-0336, United States
| | - Kirk A Peterson
- Department of Chemistry, Washington State University , Pullman Washington 99164-4630 United States
| | - David A Dixon
- Department of Chemistry, The University of Alabama , Shelby Hall, Tuscaloosa, Alabama 35487-0336, United States
| |
Collapse
|
32
|
Pan Y, Luo Z, Chang YC, Lau KC, Ng CY. High-Level ab Initio Predictions for the Ionization Energies, Bond Dissociation Energies, and Heats of Formation of Titanium Oxides and Their Cations (TiOn/TiOn+, n = 1 and 2). J Phys Chem A 2017; 121:669-679. [DOI: 10.1021/acs.jpca.6b09491] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yi Pan
- Department
of Biology and Chemistry, City University of Hong Kong, Tat Chee
Avenue, Kowloon, Hong Kong
| | - Zhihong Luo
- Department
of Chemistry, University of California, Davis, California 95616, United States
| | - Yih-Chung Chang
- Department
of Chemistry, University of California, Davis, California 95616, United States
| | - Kai-Chung Lau
- Department
of Biology and Chemistry, City University of Hong Kong, Tat Chee
Avenue, Kowloon, Hong Kong
| | - C. Y. Ng
- Department
of Chemistry, University of California, Davis, California 95616, United States
| |
Collapse
|
33
|
Letterman RG, DeYonker NJ, Burkey TJ, Webster CE. Calibrating Reaction Enthalpies: Use of Density Functional Theory and the Correlation Consistent Composite Approach in the Design of Photochromic Materials. J Phys Chem A 2016; 120:9982-9997. [PMID: 27936738 DOI: 10.1021/acs.jpca.6b09278] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
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).
Collapse
Affiliation(s)
- Roger G Letterman
- Department of Chemistry and Computational Research on Materials Institute, The University of Memphis , Memphis, Tennessee 38152, United States
| | - Nathan J DeYonker
- Department of Chemistry and Computational Research on Materials Institute, The University of Memphis , Memphis, Tennessee 38152, United States
| | - Theodore J Burkey
- Department of Chemistry and Computational Research on Materials Institute, The University of Memphis , Memphis, Tennessee 38152, United States
| | - 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
| |
Collapse
|
34
|
Calbo J, Viruela R, Ortí E, Aragó J. Relationship between Electron Affinity and Half-Wave Reduction Potential: A Theoretical Study on Cyclic Electron-Acceptor Compounds. Chemphyschem 2016; 17:3881-3890. [DOI: 10.1002/cphc.201600778] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2016] [Indexed: 11/10/2022]
Affiliation(s)
- Joaquín Calbo
- Instituto de Ciencia Molecular; Universidad de Valencia; Catedrático José Beltrán 2 46980 Paterna Spain
| | - Rafael Viruela
- Instituto de Ciencia Molecular; Universidad de Valencia; Catedrático José Beltrán 2 46980 Paterna Spain
| | - Enrique Ortí
- Instituto de Ciencia Molecular; Universidad de Valencia; Catedrático José Beltrán 2 46980 Paterna Spain
| | - Juan Aragó
- Instituto de Ciencia Molecular; Universidad de Valencia; Catedrático José Beltrán 2 46980 Paterna Spain
| |
Collapse
|
35
|
Spin-Adapted Formulation and Implementation of Density Cumulant Functional Theory with Density-Fitting Approximation: Application to Transition Metal Compounds. J Chem Theory Comput 2016; 12:4833-4842. [PMID: 27606799 DOI: 10.1021/acs.jctc.6b00589] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Density cumulant functional theory (DCT) has recently emerged as an attractive ab initio approach for the treatment of electron correlation. In its orbital-optimized formulation (ODC-12) [J. Chem. Phys. 139, 204110 (2013)], DCT has been shown to provide reliable results for a variety of challenging chemical systems. Among the attractive properties of DCT are its size-consistency and size-extensivity, as well as the efficient computation of the molecular properties and analytic gradients. In this work, we present a new formulation and implementation of DCT that takes advantage of spin adaptation and the density-fitting approximation (DF-ODC-12). Our new spin-adapted DF-ODC-12 implementation is more efficient than the previous ODC-12 implementation with up to a ∼12-fold speed-up. We demonstrate the capabilities of DF-ODC-12 with a study of transition metal compounds, which require high levels of electron correlation treatment. For transition metal carbonyl complexes [Fe(CO)5, Cr(CO)6] and the ferrocene molecule [Fe(Cp)2], the DF-ODC-12 equilibrium parameters and bond dissociation energies extrapolated to the complete basis set limit are in very good agreement with reference data derived from experiment.
Collapse
|
36
|
Fang Z, Both J, Li S, Yue S, Aprà E, Keçeli M, Wagner AF, Dixon DA. Benchmark Calculations of Energetic Properties of Groups 4 and 6 Transition Metal Oxide Nanoclusters Including Comparison to Density Functional Theory. J Chem Theory Comput 2016; 12:3689-710. [DOI: 10.1021/acs.jctc.6b00464] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Zongtang Fang
- Department
of Chemistry, The University of Alabama, Shelby Hall, Box
870336, Tuscaloosa, Alabama 35487-0336, United States
| | - Johan Both
- Department
of Chemistry, The University of Alabama, Shelby Hall, Box
870336, Tuscaloosa, Alabama 35487-0336, United States
| | - Shenggang Li
- Department
of Chemistry, The University of Alabama, Shelby Hall, Box
870336, Tuscaloosa, Alabama 35487-0336, United States
| | - Shuwen Yue
- Department
of Chemistry, The University of Alabama, Shelby Hall, Box
870336, Tuscaloosa, Alabama 35487-0336, United States
| | - Edoardo Aprà
- William
R. Wiley Environmental Molecular Sciences Laboratory, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Murat Keçeli
- Chemical
Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Albert F. Wagner
- Chemical
Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - David A. Dixon
- Department
of Chemistry, The University of Alabama, Shelby Hall, Box
870336, Tuscaloosa, Alabama 35487-0336, United States
| |
Collapse
|
37
|
Johnson EL, Davis QC, Morse MD. Predissociation measurements of bond dissociation energies: VC, VN, and VS. J Chem Phys 2016; 144:234306. [DOI: 10.1063/1.4953782] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Eric L. Johnson
- Department of Chemistry, University of Utah, Salt Lake City, Utah 84112, USA
| | - Quincy C. Davis
- Department of Chemistry, University of Utah, Salt Lake City, Utah 84112, USA
| | - Michael D. Morse
- Department of Chemistry, University of Utah, Salt Lake City, Utah 84112, USA
| |
Collapse
|
38
|
Tran QT, Tran VT. Quantum chemical study of the geometrical and electronic structures of ScSi3−/0 clusters and assignment of the anion photoelectron spectra. J Chem Phys 2016; 144:214305. [DOI: 10.1063/1.4953082] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Quoc Tri Tran
- Theoretical and Physical Chemistry Division, Dong Thap University, 783-Pham Huu Lau, Cao Lanh City, Ward 6, Dong Thap, VietNam
| | - Van Tan Tran
- Theoretical and Physical Chemistry Division, Dong Thap University, 783-Pham Huu Lau, Cao Lanh City, Ward 6, Dong Thap, VietNam
| |
Collapse
|
39
|
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-60. [DOI: 10.1021/acs.jctc.5b01163] [Citation(s) in RCA: 39] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yury Minenkov
- King Abdullah University of Science and Technology (KAUST), Physical
Science and Engineering Division (PSE), KAUST Catalysis Center (KCC), Thuwal 23955-6900, Saudi Arabia
| | - Edrisse Chermak
- King Abdullah University of Science and Technology (KAUST), Physical
Science and Engineering Division (PSE), KAUST Catalysis Center (KCC), Thuwal 23955-6900, Saudi Arabia
| | - Luigi Cavallo
- King Abdullah University of Science and Technology (KAUST), Physical
Science and Engineering Division (PSE), KAUST Catalysis Center (KCC), Thuwal 23955-6900, Saudi Arabia
| |
Collapse
|
40
|
Peterson C, Penchoff D, Wilson A. Prediction of Thermochemical Properties Across the Periodic Table. ANNUAL REPORTS IN COMPUTATIONAL CHEMISTRY 2016. [DOI: 10.1016/bs.arcc.2016.04.001] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
41
|
|
42
|
Raghavachari K, Saha A. Accurate Composite and Fragment-Based Quantum Chemical Models for Large Molecules. Chem Rev 2015; 115:5643-77. [PMID: 25849163 DOI: 10.1021/cr500606e] [Citation(s) in RCA: 183] [Impact Index Per Article: 20.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Krishnan Raghavachari
- Department of Chemistry, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Arjun Saha
- Department of Chemistry, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
43
|
DeYonker NJ, Halfen DT, Allen WD, Ziurys LM. The electronic structure of vanadium monochloride cation (VCl+): Tackling the complexities of transition metal species. J Chem Phys 2014; 141:204302. [DOI: 10.1063/1.4901239] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Affiliation(s)
- Nathan J. DeYonker
- Department of Chemistry, The University of Memphis, Memphis, Tennessee 38152, USA
| | - DeWayne T. Halfen
- Department of Chemistry, Department of Astronomy, Arizona Radio Observatory, and Steward Observatory, University of Arizona, Tucson, Arizona 85721, USA
| | - Wesley D. Allen
- Department of Chemistry and Center for Computational Chemistry, University of Georgia, Athens, Georgia 30602, USA
| | - Lucy M. Ziurys
- Department of Chemistry, Department of Astronomy, Arizona Radio Observatory, and Steward Observatory, University of Arizona, Tucson, Arizona 85721, USA
| |
Collapse
|
44
|
Otten BM, Figg TM, Cundari TR. The curious case of mesityl azide and its reactivity with bpyNiEt2. Inorg Chem 2014; 53:11633-9. [PMID: 25325403 DOI: 10.1021/ic501766k] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
A DFT analysis of the reaction of bpyNiEt2 with ArN3 was performed for para-tolyl-azide (Ar = pTol), 3,5-dimethyl-phenyl-azide (Ar = mXy) and ortho-tolyl-azide (Ar = oTol), and mesityl-azide (MesN3). Of particular interest were the different products obtained for the latter (ethylene, butane, azomesitylene, mesityl-ethylamine, etc.) versus the other reagents, i.e., (bpyNi(N(Ar)Et)(Et)). Calculated thermodynamics and kinetics for metal-free reactions did not differentiate MesN3 from the other aryl azides. Once (2)bpyNiEt(•) was generated via bond homolysis, formation of ethylene by β-H elimination was facile, as was formation of nickel-imidyl (NR(-•)) intermediates by reaction of ArN3 with bpyNiEtx (x = 0-2). On the basis of computed energetics, three reactions of bpyNiEt2 were proposed to compete: Ni-C bond homolysis, reductive elimination of butane, and nucleophilic attack (NA) by ArN3. Inspection of their temperature dependence suggested that NA and Ni-Et bond homolysis dominated at lower and higher temperatures, respectively. Calculated Ni-N and Ni-C bond dissociation free energies (BDFEs) suggested the role of radical pathways in discriminating bpyNiEt2/ArN3 reactions, and implied that the concentration of radicals such as aminyl (ArN(•)(Et)), (2)bpyNiEt(•), and Et(•) will be greatest for MesN3.
Collapse
Affiliation(s)
- Brooke M Otten
- Department of Chemistry, Center for Advanced Scientific Computing and Modeling, University of North Texas , 1508 Mulberry, Denton, Texas 76203, United States
| | | | | |
Collapse
|
45
|
MR-ccCA: A route for accurate ground and excited state potential energy curves and spectroscopic properties for third-row diatomic molecules. COMPUT THEOR CHEM 2014. [DOI: 10.1016/j.comptc.2014.04.013] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
46
|
DeYonker NJ, Shah SA. The role of core–valence electron correlation in gallium halides: a comparison of composite methods. Theor Chem Acc 2014. [DOI: 10.1007/s00214-014-1518-1] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
47
|
Marenich AV, Ho J, Coote ML, Cramer CJ, Truhlar DG. Computational electrochemistry: prediction of liquid-phase reduction potentials. Phys Chem Chem Phys 2014; 16:15068-106. [PMID: 24958074 DOI: 10.1039/c4cp01572j] [Citation(s) in RCA: 314] [Impact Index Per Article: 31.4] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
This article reviews recent developments and applications in the area of computational electrochemistry. Our focus is on predicting the reduction potentials of electron transfer and other electrochemical reactions and half-reactions in both aqueous and nonaqueous solutions. Topics covered include various computational protocols that combine quantum mechanical electronic structure methods (such as density functional theory) with implicit-solvent models, explicit-solvent protocols that employ Monte Carlo or molecular dynamics simulations (for example, Car-Parrinello molecular dynamics using the grand canonical ensemble formalism), and the Marcus theory of electronic charge transfer. We also review computational approaches based on empirical relationships between molecular and electronic structure and electron transfer reactivity. The scope of the implicit-solvent protocols is emphasized, and the present status of the theory and future directions are outlined.
Collapse
Affiliation(s)
- Aleksandr V Marenich
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant Street S.E., Minneapolis, MN 55455-0431, USA.
| | | | | | | | | |
Collapse
|
48
|
Riojas AG, Wilson AK. Solv-ccCA: Implicit Solvation and the Correlation Consistent Composite Approach for the Determination of pKa. J Chem Theory Comput 2014; 10:1500-10. [PMID: 26580366 DOI: 10.1021/ct400908z] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
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.
Collapse
Affiliation(s)
- Amanda G Riojas
- Department of Chemistry and Center for Advanced Scientific Computing and Modeling (CASCaM), University of North Texas , Denton, Texas 76203-5017, United States
| | - Angela K Wilson
- Department of Chemistry and Center for Advanced Scientific Computing and Modeling (CASCaM), University of North Texas , Denton, Texas 76203-5017, United States
| |
Collapse
|
49
|
Lu J, Yang J, Kang Y, Ning H. Probing the electronic structures and properties of neutral and anionic ScSi n (0,−1) (n = 1–6) clusters using ccCA-TM and G4 theory. J Mol Model 2014; 20:2114. [DOI: 10.1007/s00894-014-2114-y] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2013] [Accepted: 12/12/2013] [Indexed: 12/01/2022]
|
50
|
Pudasaini B, Janesko BG. Agostic Interactions in Nickel(II) Complexes: Trans Influence of Ancillary Ligands on the Strength of the Bond. Organometallics 2013. [DOI: 10.1021/om400731j] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Bimal Pudasaini
- Department
of Chemistry, Texas Christian University, Box 298860, Fort
Worth, Texas 76129, United States
| | - Benjamin G. Janesko
- Department
of Chemistry, Texas Christian University, Box 298860, Fort
Worth, Texas 76129, United States
| |
Collapse
|