1
|
Chamkin AA, Chamkina ES. Assessment of the applicability of DFT methods to [Cp*Rh]-catalyzed hydrogen evolution processes. J Comput Chem 2024. [PMID: 39052232 DOI: 10.1002/jcc.27468] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2024] [Revised: 07/04/2024] [Accepted: 07/10/2024] [Indexed: 07/27/2024]
Abstract
The present computational study provides a benchmark of density functional theory (DFT) methods in describing hydrogen evolution processes catalyzed by [Cp*Rh]-containing organometallic complexes. A test set was composed of 26 elementary reactions featuring chemical transformations and bonding situations essential for the field, including the emerging concept of non-innocent Cp* behavior. Reference values were obtained from a highly accurate 3/4 complete basis set and 6/7 complete PNO space extrapolated DLPNO-CCSD(T) energies. The performance of lower-level extrapolation procedures was also assessed. We considered 84 density functionals (DF) (including 13 generalized gradient approximations (GGA), nine meta-GGAs, 33 hybrids, and 29 double-hybrids) and three composite methods (HF-3c, PBEh-3c, and r2SCAN-3c), combined with different types of dispersion corrections (D3(0), D3BJ, D4, and VV10). The most accurate approach is the PBE0-DH-D3BJ (MAD of 1.36 kcal mol-1) followed by TPSS0-D3BJ (MAD of 1.60 kcal mol-1). Low-cost r2SCAN-3c composite provides a less accurate but much faster alternative (MAD of 2.39 kcal mol-1). The widely used Minnesota-family M06-L, M06, and M06-2X DFs should be avoided (MADs of 3.70, 3.94, and 4.01 kcal mol-1, respectively).
Collapse
Affiliation(s)
- Aleksandr A Chamkin
- A.N.Nesmeyanov Institute of Organoelement Compounds of Russian Academy of Sciences, Moscow, Russia
| | - Elena S Chamkina
- A.N.Nesmeyanov Institute of Organoelement Compounds of Russian Academy of Sciences, Moscow, Russia
| |
Collapse
|
2
|
Abraham V, Harsha G, Zgid D. Relativistic Fully Self-Consistent GW for Molecules: Total Energies and Ionization Potentials. J Chem Theory Comput 2024; 20:4579-4590. [PMID: 38778459 DOI: 10.1021/acs.jctc.4c00075] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/25/2024]
Abstract
The fully self-consistent GW (scGW) method with an iterative solution of the Dyson equation provides a consistent approach for describing the ground and excited states without any dependence on the mean-field reference. In this work, we present a relativistic version of scGW for molecules containing heavy elements using the exact two-component (X2C) Coulomb approximation. We benchmark SOC-81 data set containing closed shell heavy elements for the first ionization potential using the fully self-consistent GW as well as one-shot GW. The self-consistent GW provides superior results compared to G0W0 with PBE reference and comparable results to G0W0 with PBE0 while also removing the starting point dependence. The photoelectron spectra obtained at the X2C level demonstrate very good agreement with the experimental spectra. We also observe that scGW provides very good estimation of ionization potential for the inner d-shell orbitals. Additionally, using the well-conserved total energy, we investigate the equilibrium bond length and harmonic frequencies of a few halogen dimers using scGW. Overall, our findings demonstrate the applicability of the fully self-consistent GW method for accurate ionization potential, photoelectron spectra, and total energies in finite systems with heavy elements with a reasonable computational scaling.
Collapse
Affiliation(s)
- Vibin Abraham
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Gaurav Harsha
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Dominika Zgid
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
- Department of Physics and Astronomy, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
3
|
Drabik G, Radoń M. Approaching the Complete Basis Set Limit for Spin-State Energetics of Mononuclear First-Row Transition Metal Complexes. J Chem Theory Comput 2024; 20:3199-3217. [PMID: 38574194 PMCID: PMC11044276 DOI: 10.1021/acs.jctc.4c00092] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2024] [Revised: 03/19/2024] [Accepted: 03/20/2024] [Indexed: 04/06/2024]
Abstract
Convergence to the complete basis set (CBS) limit is analyzed for the problem of spin-state energetics in mononuclear first-row transition metal (TM) complexes by taking under scrutiny a benchmark set of 18 energy differences between spin states for 13 chemically diverse TM complexes. The performance of conventional CCSD(T) and explicitly correlated CCSD(T)-F12a/b calculations in approaching the CCSD(T)/CBS limits is systematically studied. An economic computational protocol is developed based on the CCSD-F12a approximation and (here proposed) modified scaling of the perturbative triples term (T#). This computational protocol recovers the relative spin-state energetics of the benchmark set in excellent agreement with the reference CCSD(T)/CBS limits (mean absolute deviation of 0.4, mean signed deviation of 0.2, and maximum deviation of 0.8 kcal/mol) and enables performing canonical CCSD(T) calculations for mononuclear TM complexes sized up to ca. 50 atoms, which is illustrated by application to heme-related metalloporphyrins. Furthermore, a good transferability of the basis set incompleteness error (BSIE) is demonstrated for spin-state energetics computed using CCSD(T) and other wave function methods (MP2, CASPT2, CASPT2/CC, NEVPT2, and MRCI + Q), which justifies efficient focal-point approximations and simplifies the construction of multimethod benchmark studies.
Collapse
Affiliation(s)
- Gabriela Drabik
- Jagiellonian
University, Doctoral School
of Exact and Natural Sciences, Łojasiewicza 11, 30-348 Kraków, Poland
- Jagiellonian
University, Faculty of Chemistry, Gronostajowa 2, 30-387, Kraków Poland
| | - Mariusz Radoń
- Jagiellonian
University, Faculty of Chemistry, Gronostajowa 2, 30-387, Kraków Poland
| |
Collapse
|
4
|
Chan B. DAPD Set of Pd-Containing Diatomic Molecules: Accurate Molecular Properties and the Great Lengths to Obtain Them. J Chem Theory Comput 2023; 19:9260-9268. [PMID: 38096563 DOI: 10.1021/acs.jctc.3c01060] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2023]
Abstract
In the present study, we obtained reliable bond energy, bond length, and zero-point vibrational frequency for a set of diatomic Pd species (the DAPD set). It includes PdH, Pd2, and PdX (X = B, C, N, O, F, Al, Si, P, S, and Cl). Our highest-level protocol (W4X-L) represents scalar and spin-orbit relativistic, valence- and inner-valence correlated, extrapolated CCSDTQ(5) energy. The DAPD set of molecules is challenging for computational chemistry methods in different manners; for Pd2, the spin-orbit contribution to the bond energy is fairly large, whereas for PdC and PdSi, the post-CCSD(T) correlation components are considerable. The diverse range of requirements represents a significant challenge for lower-level methods. While density functional theory (DFT) methods generally yield good agreements for bond lengths and vibrational frequencies, large deviations are found for bond energies. In general, hybrid DFT methods are more accurate than nonhybrid functionals, but the agreement in individual cases varies. This illustrates the critical role that new high-quality reference data would play in the continual development of lower-cost methods.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo 1-14, Nagasaki 852-8521, Japan
| |
Collapse
|
5
|
Radoń M. Benchmarks for transition metal spin-state energetics: why and how to employ experimental reference data? Phys Chem Chem Phys 2023; 25:30800-30820. [PMID: 37938035 DOI: 10.1039/d3cp03537a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2023]
Abstract
Accurate prediction of energy differences between alternative spin states of transition metal complexes is essential in computational (bio)inorganic chemistry-for example, in characterization of spin crossover materials and in the theoretical modeling of open-shell reaction mechanisms-but it remains one of the most compelling problems for quantum chemistry methods. A part of this challenge is to obtain reliable reference data for benchmark studies, as even the highest-level applicable methods are known to give divergent results. This Perspective discusses two possible approaches to method benchmarking for spin-state energetics: using either theoretically computed or experiment-derived reference data. With the focus on the latter approach, an extensive general review is provided for the available experimental data of spin-state energetics and their interpretations in the context of benchmark studies, targeting the possibility of back-correcting the vibrational effects and the influence of solvents or crystalline environments. With a growing amount of experience, these effects can be now not only qualitatively understood, but also quantitatively modeled, providing the way to derive nearly chemically accurate estimates of the electronic spin-state gaps to be used as benchmarks and advancing our understanding of the phenomena related to spin states in condensed phases.
Collapse
Affiliation(s)
- Mariusz Radoń
- Faculty of Chemistry, Jagiellonian University, Gronostajowa 2, 30-387 Krakow, Krakow, Poland.
| |
Collapse
|
6
|
Neugebauer H, Vuong HT, Weber JL, Friesner RA, Shee J, Hansen A. Toward Benchmark-Quality Ab Initio Predictions for 3d Transition Metal Electrocatalysts: A Comparison of CCSD(T) and ph-AFQMC. J Chem Theory Comput 2023; 19:6208-6225. [PMID: 37655473 DOI: 10.1021/acs.jctc.3c00617] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/02/2023]
Abstract
Generating accurate ab initio ionization energies for transition metal complexes is an important step toward the accurate computational description of their electrocatalytic reactions. Benchmark-quality data is required for testing existing theoretical methods and developing new ones but is complicated to obtain for many transition metal compounds due to the potential presence of both strong dynamical and static electron correlation. In this regime, it is questionable whether the so-called gold standard, coupled cluster with singles, doubles, and perturbative triples (CCSD(T)), provides the desired level of accuracy─roughly 1-3 kcal/mol. In this work, we compiled a test set of 28 3d metal-containing molecules relevant to homogeneous electrocatalysis (termed 3dTMV) and computed their vertical ionization energies (ionization potentials) with CCSD(T) and phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) in the def2-SVP basis set. A substantial effort has been made to converge away the phaseless bias in the ph-AFQMC reference values. We assess a wide variety of multireference diagnostics and find that spin-symmetry breaking of the CCSD wave function and the PBE0 density functional correlate well with our analysis of multiconfigurational wave functions. We propose quantitative criteria based on symmetry breaking to delineate correlation regimes inside of which appropriately performed CCSD(T) can produce mean absolute deviations from the ph-AFQMC reference values of roughly 2 kcal/mol or less and outside of which CCSD(T) is expected to fail. We also present a preliminary assessment of density functional theory (DFT) functionals on the 3dTMV set.
Collapse
Affiliation(s)
- Hagen Neugebauer
- Mulliken Center for Theoretical Chemistry, Clausius Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, D-53115 Bonn, Germany
| | - Hung T Vuong
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - John L Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - James Shee
- Department of Chemistry, Rice University, Houston, Texas 77005, United States
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Clausius Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, D-53115 Bonn, Germany
| |
Collapse
|
7
|
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
|
8
|
Rossomme E, Cunha LA, Li W, Chen K, McIsaac AR, Head-Gordon T, Head-Gordon M. The Good, the Bad, and the Ugly: Pseudopotential Inconsistency Errors in Molecular Applications of Density Functional Theory. J Chem Theory Comput 2023; 19:2827-2841. [PMID: 37156013 DOI: 10.1021/acs.jctc.3c00089] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
The pseudopotential (PP) approximation is one of the most common techniques in computational chemistry. Despite its long history, the development of custom PPs has not tracked with the explosion of different density functional approximations (DFAs). As a result, the use of PPs with exchange/correlation models for which they were not developed is widespread, although this practice is known to be theoretically unsound. The extent of PP inconsistency errors (PPIEs) associated with this practice has not been systematically explored across the types of energy differences commonly evaluated in chemical applications. We evaluate PPIEs for a number of PPs and DFAs across 196 chemically relevant systems of both transition-metal and main-group elements, as represented by the W4-11, TMC34, and S22 data sets. Near the complete basis set limit, these PPs are found to cleanly approach all-electron (AE) results for noncovalent interactions but introduce root-mean-squared errors (RMSEs) upwards of 15 kcal mol-1 into predictions of covalent bond energies for a number of popular DFAs. We achieve significant improvements through the use of empirical atom- and DFA-specific PP corrections, indicating considerable systematicity of the PPIEs. The results of this work have implications for chemical modeling in both molecular contexts and for DFA design, which we discuss.
Collapse
Affiliation(s)
- Elliot Rossomme
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Leonardo A Cunha
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Wanlu Li
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Kaixuan Chen
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Alexandra R McIsaac
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Teresa Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
- Department of Bioengineering, University of California, Berkeley, California 94720, United States
- Department of Chemical and Biomolecular Engineering, University of California, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
9
|
Sweeny BC, Long BA, Maffucci D, Zuo J, Guo H, Viggiano AA, Ard SG, Shuman NS. Activation of Methane by Zr +: A Deep-Dive into the Potential Surface via Pressure- and Temperature-Dependent Kinetics with Statistical Modeling. J Phys Chem A 2023; 127:1818-1830. [PMID: 36802591 DOI: 10.1021/acs.jpca.2c07584] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/22/2023]
Abstract
The kinetics of Zr+ + CH4 are measured using a selected-ion flow tube apparatus over the temperature range 300-600 K and the pressure range 0.25-0.60 Torr. Measured rate constants are small, never exceeding 5% of the Langevin capture value. Both collisionally stabilized ZrCH4+ and bimolecular ZrCH2+ products are observed. A stochastic statistical modeling of the calculated reaction coordinate is used to fit the experimental results. The modeling indicates that an intersystem crossing from the entrance well, necessary for the bimolecular product to be formed, occurs faster than competing isomerization and dissociation processes. That sets an upper limit on the lifetime of the entrance complex to crossing of 10-11 s. The endothermicity of the bimolecular reaction is derived to be 0.09 ± 0.05 eV, in agreement with a literature value. The observed ZrCH4+ association product is determined to be primarily HZrCH3+ not Zr+(CH4), indicating that bond activation has occurred at thermal energies. The energy of HZrCH3+ relative to separated reactants is determined to be -0.80 ± 0.25 eV. Inspection of the statistical modeling results under best-fit conditions reveals reaction dependences on impact parameter, translation energy, internal energy, and angular momentum. Reaction outcomes are heavily affected by angular momentum conservation. Additionally, product energy distributions are predicted.
Collapse
Affiliation(s)
- Brendan C Sweeny
- Boston College Institute for Scientific Research, Boston, Massachusetts 02549, United States
| | - Bryan A Long
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland Air Force Base, New Mexico 87117, United States
| | - Dominique Maffucci
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland Air Force Base, New Mexico 87117, United States
| | - Junxiang Zuo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Hua Guo
- Department of Chemistry and Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| | - Albert A Viggiano
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland Air Force Base, New Mexico 87117, United States
| | - Shaun G Ard
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland Air Force Base, New Mexico 87117, United States
| | - Nicholas S Shuman
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland Air Force Base, New Mexico 87117, United States
| |
Collapse
|
10
|
Liu X, McKemmish L, Pérez-Ríos J. The performance of CCSD(T) for the calculation of dipole moments in diatomics. Phys Chem Chem Phys 2023; 25:4093-4104. [PMID: 36651174 DOI: 10.1039/d2cp05060a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
This work analyzes the accuracy of the coupled cluster with single, double, and perturbative triple excitation [CCSD(T)] method for predicting dipole moments. In particular, we benchmark CCSD(T) predictions for the equilibrium bond length, vibrational frequency, and dipole moment versus accurate experimental data. As a result, we find that CCSD(T) leads to accurate dipole moments. However, in some cases, it disagrees with the experimental values, and the disagreement can not be satisfactorily explained via relativistic or multi-reference effects. Therefore, our results indicate that benchmark studies for energy and geometry properties do not accurately describe other electron density magnitudes.
Collapse
Affiliation(s)
- Xiangyue Liu
- Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany
| | - Laura McKemmish
- School of Chemistry, UNSW Sydney, Sydney, NSW 2052, Australia
| | - Jesús Pérez-Ríos
- Department of Physics and Astronomy, Stony Brook University, Stony Brook 11794, New York, USA. .,Institute for Advanced Computational Science, Stony Brook University, Stony Brook, NY 11794-3800, USA
| |
Collapse
|
11
|
Zhang C, Cheng L. Route to Chemical Accuracy for Computational Uranium Thermochemistry. J Chem Theory Comput 2022; 18:6732-6741. [PMID: 36206308 DOI: 10.1021/acs.jctc.2c00812] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Benchmark spinor-based relativistic coupled-cluster calculations for the ionization energies of the uranium atom, the uranium monoxide molecule (UO), and the uranium dioxide molecule (UO2) and for the bond dissociation energies of UO and UO2 are reported. The accuracy of these calculations in the treatments of relativistic, electron-correlation, and basis-set effects is analyzed. The intrinsic convergence of the computed results and the favorable comparison with the experimental values demonstrate the unique applicability of the spinor representation of quantum-chemical methods to open-shell uranium-containing atomic and molecular species with uranium oxidation states ranging from U(0) to U(V).
Collapse
Affiliation(s)
- Chaoqun Zhang
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Lan Cheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
12
|
Ezarfi N, Benjelloun AT, Benzakour M, Mcharfi M. Electronic structure and thermochemistry for monocarbides MC, MC+ and MC− (M=Zn, Cd, Hg): CCSD(T) and DFT works. J INDIAN CHEM SOC 2022. [DOI: 10.1016/j.jics.2022.100600] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
13
|
Varberg TD. First detection and analysis of an electronic spectrum of vanadium hydride: The D 5Π-X 5Δ (0,0) band. J Chem Phys 2022; 157:074311. [DOI: 10.1063/5.0105844] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The D5Π-X5Δ (0,0) band of vanadium hydride at 654 nm has been recorded by laser excitation spectroscopy and represents the first analyzed spectrum of VH in the gas phase. The molecules were generated using a hollow cathode discharge source, with laser-induced fluorescence detected via the D5Π-A5Π (0,0) transition. All five main (ΔΩ = ΔΛ) subbands were observed as well as with several satellite ones, which together create a rather complex and overlapped spectrum covering the region 15180-15500 cm-1. The D5Π state displays the effects of three strong local perturbations, which are likely caused by interactions with high vibrational levels of the B5Σ- and c3Σ- states, identified in a previous multiconfigurational self-consistent field study by Koseki et al. [J. Phys. Chem. A 108, 4707 (2004)]. Molecular constants describing the X5Δ, A5Π, and D5Π states were determined in three separate least-squares fits using effective Hamiltonians written in a Hund's case (a) basis. The fine structure of the ground state is found to be consistent with its assignment as a σπ2δ, 5Δ electronic state. The fitted values of its first-order spin-orbit and rotational constants in the ground state are A = 36.5378 cm-1 and B = 5.7579 cm-1, the latter of which yields a bond length of R0 = 1.7212(2) Å. This experimental value is in good agreement with previous computational studies of the molecule and fits well within the overall trend of decreasing bond length across the series of 3d transition metal monohydrides.
Collapse
Affiliation(s)
- Thomas D. Varberg
- Department of Chemistry, Macalester College, United States of America
| |
Collapse
|
14
|
Rudshteyn B, Weber JL, Coskun D, Devlaminck PA, Zhang S, Reichman DR, Shee J, Friesner RA. Calculation of Metallocene Ionization Potentials via Auxiliary Field Quantum Monte Carlo: Toward Benchmark Quantum Chemistry for Transition Metals. J Chem Theory Comput 2022; 18:2845-2862. [PMID: 35377642 PMCID: PMC9123894 DOI: 10.1021/acs.jctc.1c01071] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
The accurate ab initio prediction of ionization energies is essential to understanding the electrochemistry of transition metal complexes in both materials science and biological applications. However, such predictions have been complicated by the scarcity of gas phase experimental data, the relatively large size of the relevant molecules, and the presence of strong electron correlation effects. In this work, we apply all-electron phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) utilizing multideterminant trial wave functions to six metallocene complexes to compare the computed adiabatic and vertical ionization energies with experimental results. We find that ph-AFQMC yields mean absolute errors (MAEs) of 1.69 ± 1.02 kcal/mol for the adiabatic energies and 2.85 ± 1.13 kcal/mol for the vertical energies. We also carry out density functional theory (DFT) calculations using a variety of functionals, which yields MAEs of 3.62-6.98 kcal/mol and 3.31-9.88 kcal/mol, as well as one variant of localized coupled cluster calculations (DLPNO-CCSD(T0) with moderate PNO cutoffs), which has MAEs of 4.96 and 6.08 kcal/mol, respectively. We also test the reliability of DLPNO-CCSD(T0) and DFT on acetylacetonate (acac) complexes for adiabatic energies measured in the same manner experimentally, and we find higher MAEs, ranging from 4.56 to 10.99 kcal/mol (with a different ordering) for DFT and 6.97 kcal/mol for DLPNO-CCSD(T0). Finally, by utilizing experimental solvation energies, we show that accurate reduction potentials in solution for the metallocene series can be obtained from the AFQMC gas phase results.
Collapse
Affiliation(s)
- Benjamin Rudshteyn
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - John L Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Dilek Coskun
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Pierre A Devlaminck
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, 162 Fifth Avenue, New York, New York 10010, United States.,Department of Physics, College of William and Mary, Williamsburg, Virginia 23187, United States
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - James Shee
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
15
|
Sweeny BC, Long BA, Viggiano AA, Ard SG, Shuman NS. Effect of Intersystem Crossings on the Kinetics of Thermal Ion-Molecule Reactions: Ti + + O 2, CO 2, and N 2O. J Phys Chem A 2022; 126:859-869. [PMID: 35107288 DOI: 10.1021/acs.jpca.1c10196] [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
A selected-ion flow tube apparatus has been used to measure rate constants and product branching fractions of 2Ti+ reacting with O2, CO2, and N2O over the range of 200-600 K. Ti+ + O2 proceeds at near the Langevin capture rate constant of 6-7 × 10-10 cm3 s-1 at all temperatures to yield 4TiO+ + O. Reactions initiated on doublet or quartet surfaces are formally spin-allowed; however, the 50% of reactions initiated on sextet surfaces must undergo an intersystem crossing (ISC). Statistical theory is used to calculate the energy and angular momentum dependences of the specific rate constants for the competing isomerization and dissociation channels. This acts as an internal clock on the lifetime to ISC, setting an upper limit on the order of τISC < 1e-11 s. 2Ti+ + CO2 produces 4TiO+ + CO less efficiently, with a rate constant fit as 5.5 ± 1.3 × 10-11 (T/300)-1.1 ± 0.2 cm3 s-1. The reaction is formally spin-prohibited, and statistical modeling shows that ISC, not a submerged transition state, is rate-limiting, occurring with a lifetime on the order of 10-7 s. Ti+ + N2O proceeds at near the capture rate constant. In this case, both Ti+ON2 and Ti+N2O entrance channel complexes are formed and can interconvert over a barrier. The main product is >90% TiO+ + N2, and the remainder is TiN+ + NO. Both channels need to undergo ISC to form ground-state products but TiO+ can be formed in an excited state exothermically. Therefore, kinetic information is obtained only for the TiN+ channel, where ISC occurs with a lifetime on the order of 10-9 s. Statistical modeling indicates that the dipole-preferred Ti+ON2 complex is formed in ∼80% of collisions, and this value is reproduced using a capture model based on the generic ion-dipole + quadrupole long-range potential.
Collapse
Affiliation(s)
- Brendan C Sweeny
- Institute for Scientific Research, Boston College, Boston, Massachusetts 02467, United States
| | - Bryan A Long
- NRC Postdoc at Air Force Research Laboratory, Space Vehicles Directorate, Kirtland Air Force Base, New Mexico 87117, United States
| | - Albert A Viggiano
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland Air Force Base, New Mexico 87117, United States
| | - Shaun G Ard
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland Air Force Base, New Mexico 87117, United States
| | - Nicholas S Shuman
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland Air Force Base, New Mexico 87117, United States
| |
Collapse
|
16
|
Minenkova I, Otlyotov AA, Cavallo L, Minenkov Y. Gas-phase thermochemistry of polycyclic aromatic hydrocarbons: an approach integrating the quantum chemistry composite scheme and reaction generator. Phys Chem Chem Phys 2022; 24:3163-3181. [PMID: 35040851 DOI: 10.1039/d1cp03702a] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We introduce a protocol aimed at predicting the accurate gas-phase enthalpies of formation of polycyclic aromatic hydrocarbons (PAHs). Automatic generation of a dataset of equilibrated chemical reactions preserving the number of carbon atoms in each hybridization state on each side of equations is at the core of our scheme. The performed tests suggest the recommended enthalpy of formation to be derived via a two-step scheme. First, we consider the reactions with a minimal sum of the total number of particles involved, N, and the absolute difference between the total number of products and reactants, |ΔN|. Second, among these reactions, we identify the one with the smallest absolute reaction enthalpy change, . This approach has been applied to predict the gas-phase enthalpies of formation of 113 PAHs via the Feller-Peterson-Dixon approach. Our calculated values provide the mean absolute deviations of 1.7, 1.9, 4.2, 8.1, and 18.5 kJ mol-1 with respect to the literature group-based error corrected (GBEC) G3MP2B3, ATOMIC (HC), group equivalent M06-2X, GBEC B3LYP, and G4MP2 values. Our predicted values give the mean signed and mean absolute errors of -7.5 and 12.9 kJ mol-1 with respect to the experimental enthalpies of formation. The combination of our predicted and the experimental values provide the solid-state enthalpies of formation, , which are not available for a few species. Approaching these values as well as , producing large discrepancies from the experimental side, would be indispensable for testing and further tuning of computational chemistry approaches.
Collapse
Affiliation(s)
- Irina Minenkova
- Frumkin Institute of Physical Chemistry and Electrochemistry, Russian Academy of Sciences, Leninskii pr. 31, Moscow 119071, Russian Federation
| | - Arseniy A Otlyotov
- N. N. Semenov Federal Research Center for Chemical Physics RAS, Kosygina Street 4, 119991 Moscow, Russian Federation.
| | - Luigi Cavallo
- KAUST Catalysis Center (KCC), King Abdullah University of Science and Technology, Thuwal-23955-6900, Saudi Arabia.
| | - Yury Minenkov
- N. N. Semenov Federal Research Center for Chemical Physics RAS, Kosygina Street 4, 119991 Moscow, Russian Federation. .,Joint Institute for High Temperatures, Russian Academy of Sciences, 13-2 Izhorskaya Street, Moscow 125412, Russian Federation
| |
Collapse
|
17
|
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
|
18
|
Nandy A, Duan C, Taylor MG, Liu F, Steeves AH, Kulik HJ. Computational Discovery of Transition-metal Complexes: From High-throughput Screening to Machine Learning. Chem Rev 2021; 121:9927-10000. [PMID: 34260198 DOI: 10.1021/acs.chemrev.1c00347] [Citation(s) in RCA: 70] [Impact Index Per Article: 23.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Transition-metal complexes are attractive targets for the design of catalysts and functional materials. The behavior of the metal-organic bond, while very tunable for achieving target properties, is challenging to predict and necessitates searching a wide and complex space to identify needles in haystacks for target applications. This review will focus on the techniques that make high-throughput search of transition-metal chemical space feasible for the discovery of complexes with desirable properties. The review will cover the development, promise, and limitations of "traditional" computational chemistry (i.e., force field, semiempirical, and density functional theory methods) as it pertains to data generation for inorganic molecular discovery. The review will also discuss the opportunities and limitations in leveraging experimental data sources. We will focus on how advances in statistical modeling, artificial intelligence, multiobjective optimization, and automation accelerate discovery of lead compounds and design rules. The overall objective of this review is to showcase how bringing together advances from diverse areas of computational chemistry and computer science have enabled the rapid uncovering of structure-property relationships in transition-metal chemistry. We aim to highlight how unique considerations in motifs of metal-organic bonding (e.g., variable spin and oxidation state, and bonding strength/nature) set them and their discovery apart from more commonly considered organic molecules. We will also highlight how uncertainty and relative data scarcity in transition-metal chemistry motivate specific developments in machine learning representations, model training, and in computational chemistry. Finally, we will conclude with an outlook of areas of opportunity for the accelerated discovery of transition-metal complexes.
Collapse
Affiliation(s)
- Aditya Nandy
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Chenru Duan
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Michael G Taylor
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Fang Liu
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Adam H Steeves
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
19
|
Welch BK, Almeida NMS, Wilson AK. Super ccCA (s-ccCA): an approach for accurate transition metal thermochemistry. Mol Phys 2021. [DOI: 10.1080/00268976.2021.1963001] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Bradley K. Welch
- Department of Chemistry, Michigan State University, East Lansing, MI, USA
| | - Nuno M. S. Almeida
- Department of Chemistry, Michigan State University, East Lansing, MI, USA
| | - Angela K. Wilson
- Department of Chemistry, Michigan State University, East Lansing, MI, USA
| |
Collapse
|
20
|
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
|
21
|
Ard SG, Viggiano AA, Shuman NS. Old School Techniques with Modern Capabilities: Kinetics Determination of Dynamical Information Such as Barriers, Multiple Entrance Channel Complexes, Product States, Spin Crossings, and Size Effects in Metallic Ion–Molecule Reactions. J Phys Chem A 2021; 125:3503-3527. [DOI: 10.1021/acs.jpca.0c11395] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Affiliation(s)
- Shaun G. Ard
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland Air Force Base, Albuquerque, New Mexico 87117, United States
| | - Albert A. Viggiano
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland Air Force Base, Albuquerque, New Mexico 87117, United States
| | - Nicholas S. Shuman
- Air Force Research Laboratory, Space Vehicles Directorate, Kirtland Air Force Base, Albuquerque, New Mexico 87117, United States
| |
Collapse
|
22
|
Franke PR, Stanton JF, Douberly GE. How to VPT2: Accurate and Intuitive Simulations of CH Stretching Infrared Spectra Using VPT2+K with Large Effective Hamiltonian Resonance Treatments. J Phys Chem A 2021; 125:1301-1324. [PMID: 33506678 DOI: 10.1021/acs.jpca.0c09526] [Citation(s) in RCA: 65] [Impact Index Per Article: 21.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
This article primarily discusses the utility of vibrational perturbation theory for the prediction of X-H stretching vibrations with particular focus on the specific variant, second-order vibrational perturbation theory with resonances (VPT2+K). It is written as a tutorial, reprinting most important formulas and providing numerous simple examples. It discusses the philosophy and practical considerations behind vibrational simulations with VPT2+K, including but not limited to computational method selection, cost-saving approximations, approaches to evaluating intensity, resonance identification, and effective Hamiltonian structure. Particular attention is given to resonance treatments, beginning with simple Fermi dyads and gradually progressing to arbitrarily large polyads that describe both Fermi and Darling-Dennison resonances. VPT2+K combined with large effective Hamiltonians is shown to be a reliable framework for modeling the complicated CH stretching spectra of alkenes. An error is also corrected in the published analytic formula for the VPT2 transition moment between the vibrational ground state and triply excited states.
Collapse
Affiliation(s)
- Peter R Franke
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States.,Department of Chemistry, University of Georgia, Athens, Georgia 30602, United States
| | - John F Stanton
- Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| | - Gary E Douberly
- Department of Chemistry, University of Georgia, Athens, Georgia 30602, United States
| |
Collapse
|
23
|
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
|
24
|
Dorofeeva OV, Ryzhova ON. Accurate estimation of enthalpies of formation for C-, H-, O-, and N-containing compounds using DLPNO-CCSD(T1)/CBS method. Struct Chem 2021. [DOI: 10.1007/s11224-020-01681-y] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
|
25
|
Minenkova I, Osina EL, Cavallo L, Minenkov Y. Gas-Phase Thermochemistry of MX 3 and M 2X 6 (M = Sc, Y; X = F, Cl, Br, I) from a Composite Reaction-Based Approach: Homolytic versus Heterolytic Cleavage. Inorg Chem 2020; 59:17084-17095. [PMID: 33210914 DOI: 10.1021/acs.inorgchem.0c02292] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A domain-based local-pair natural-orbital coupled-cluster approach with single, double, and improved linear-scaling perturbative triple correction via an iterative algorithm, DLPNO-CCSD(T1), was applied within the framework of the Feller-Peterson-Dixon approach to derive gas-phase heats of formation of scandium and yttrium trihalides and their dimers via a set of homolytic and heterolytic dissociation reactions. All predicted heats of formation moderately depend on the reaction type with the most and least negative values obtained for homolytic and heterolytic dissociation, respectively. The basis set size dependence, as well as the influence of static correlation effects not covered by the standard (DLPNO-)CCSD(T) approach, suggests that exploitation of the heterolytic dissociation reactions with the formation of M3+ and X- ions leads to the most robust heats of formation. The gas-phase formation enthalpies ΔHf°(0 K)/ΔHf°(298.15 K) and absolute entropies S°(298.15 K) were obtained for the first time for the Sc2F6, Sc2Br6, and Sc2I6 species. For ScBr3, ScI3, Sc2Cl6, and Y2Cl6, we suggest a reexamination of the experimental heats of formation available in the literature. For other compounds, the predicted values were found to be in good agreement with the experimental estimates. Extracted MX3 (M = Sc, Y; X = F, Cl, Br, and I) 0 K atomization enthalpies indicate weaker bonding when moving from fluorine to iodine and from yttrium to scandium. Likewise, the stability of yttrium trihalide dimers degrades when going from fluorine to iodine. Respective scandium trihalide dimers are less stable, with 0 K dimer dissociation energy decreasing in the row fluorine - chlorine - bromine ≈ iodine. Correlation of the (n - 1)s2p6 electrons on bromine and iodine, inclusion of zero-point energy, relativistic effects, and the effective-core-potential correction as well as amelioration of the DLPNO localization inaccuracy are shown to be of similar magnitude, which is critical if accurate heats of formation are a goal.
Collapse
Affiliation(s)
- Irina Minenkova
- Frumkin Institute of Physical Chemistry and Electrochemistry, Russian Academy of Sciences, Leninskii pr. 31, Moscow 119071, Russia
| | - Evgeniya L Osina
- Joint Institute for High Temperatures, Russian Academy of Sciences, 13-2 Izhorskaya Street, Moscow 125412, Russia
| | - Luigi Cavallo
- Catalysis Center, King Abdullah University of Science and Technology (KAUST), Thuwal 23955-6900, Saudi Arabia
| | - Yury Minenkov
- Joint Institute for High Temperatures, Russian Academy of Sciences, 13-2 Izhorskaya Street, Moscow 125412, Russia.,N.N. Semenov Federal Research Center for Chemical Physics, Russian Academy of Sciences, Kosygina Street 4, Moscow 119991, Russia
| |
Collapse
|
26
|
Bao JL, Welch BK, Ulusoy IS, Zhang X, Xu X, Wilson AK, Truhlar DG. Predicting Bond Dissociation Energies and Bond Lengths of Coordinatively Unsaturated Vanadium-Ligand Bonds. J Phys Chem A 2020; 124:9757-9770. [PMID: 33180508 DOI: 10.1021/acs.jpca.0c06519] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Understanding the electronic structure of coordinatively unsaturated transition-metal compounds and predicting their physical properties are of great importance for catalyst design. Bond dissociation energy De and bond length re are two of the fundamental quantities for which good predictions are important for a successful design strategy. In the present work, recent experimentally measured bond energies and bond lengths of VX diatomic molecules (X = C, N, S) are used as a gauge to consider the utility of a number of electronic structure methods. Single-reference methods are one focus because of their efficiency and utility in practical calculations, and multireference configuration interaction (MRCISD) methods and a composite coupled cluster (CCC) method are a second focus because of their potential high accuracy. The comparison is especially challenging because of the large multireference M diagnostics of these molecules, in the range 0.15-0.19. For the single-reference methods, Kohn-Sham density functional theory (KS-DFT) has been tested with a variety of approximate exchange-correlation functionals. Of these, MOHLYP provides the bond dissociation energies in best agreement with experiments, and BLYP provides the bond lengths that are in best agreement with experiments; but by requiring good performance for both the De and re of the vanadium compounds, MOHLYP, MN12-L, MGGA_MS1, MGGA_MS0, O3LYP, and M06-L are the most highly recommended functionals. The CCC calculations include up to connected pentuple excitations for the valence electrons and up to connected quadruple excitations for the core-valence terms; this results in highly accurate dissociation energies and good bond lengths. Averaged over the three molecules, the mean unsigned deviation of CCC bond energies from experimental ones is only 0.4 kcal/mol, demonstrating excellent convergence of theory and experiments.
Collapse
Affiliation(s)
- Junwei Lucas Bao
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Bradley K Welch
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, United States
| | - Inga S Ulusoy
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, United States.,Theoretical Chemistry, Institute of Physical Chemistry, Heidelberg University, Im Neuenheimer Feld 229, Heidelberg 69120, Germany
| | - Xin Zhang
- State Key Laboratory of Chemical Resource Engineering, Beijing Advanced Innovation Center for Soft Matter Science and Engineering, Beijing University of Chemical Technology, Beijing 100029, People's Republic of China.,Department of Chemistry, Chemical Theory Center, Inorganometallic Catalyst Design Center, and Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| | - Xuefei Xu
- Center for Combustion Energy, Department of Energy and Power Engineering, and Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Tsinghua University, Beijing 100084, China
| | - Angela K Wilson
- Department of Chemistry, Michigan State University, 578 S Shaw Lane, East Lansing, Michigan 48824-1322, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, Inorganometallic Catalyst Design Center, and Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
27
|
Chamkin AA, Serkova ES. DFT, DLPNO-CCSD(T), and NEVPT2 benchmark study of the reaction between ferrocenium and trimethylphosphine. J Comput Chem 2020; 41:2388-2397. [PMID: 32812657 DOI: 10.1002/jcc.26398] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2020] [Accepted: 07/25/2020] [Indexed: 11/10/2022]
Abstract
The reaction between ferrocenium and trimethylphosphine was studied using density functional theory (DFT), domain-based local pair natural orbital coupled cluster theory with single-, double-, and perturbative triple excitations (DLPNO-CCSD(T)), and N-electron valence state perturbation theory (NEVPT2). The accuracy of the DFT functionals decreases compared to the DLPNO-CCSD(T) level in the following order: M06-L > TPSS > M06, BLYP > PBE, PBE0, B3LYP > > PWPB95 > > DSD-BLYP. The roles of thermochemical, continuum solvation (SMD), and counterpoise corrections were evaluated. Grimme's D3 empirical dispersion correction is essential for all functionals studied except M06 and M06-L. The reliability of the frequency calculations performed directly within the SMD was confirmed. The systems showed no significant multireference character according to T1 and T2 diagnostics and the fractional occupation number (FOD) weighted electron density analysis. The multireference NEVPT2 calculations gave qualitatively valid conclusions about the reaction mechanism. However, a multireference approach is generally not recommended because it requires arbitrary chosen active spaces.
Collapse
Affiliation(s)
- Aleksandr A Chamkin
- A.N. Nesmeyanov Institute of Organoelement Compounds of Russian Academy of Sciences, Moscow, Russia
| | - Elena S Serkova
- A.N. Nesmeyanov Institute of Organoelement Compounds of Russian Academy of Sciences, Moscow, Russia
| |
Collapse
|
28
|
Moltved KA, Kepp KP. Dioxygen Binding to all 3d, 4d, and 5d Transition Metals from Coupled-Cluster Theory. Chemphyschem 2020; 21:2173-2186. [PMID: 32757346 DOI: 10.1002/cphc.202000529] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2020] [Revised: 08/04/2020] [Indexed: 11/11/2022]
Abstract
Understanding how transition metals bind and activate dioxygen (O2 ) is limited by experimental and theoretical uncertainties, making accurate quantum mechanical descriptors of interest. Here we report coupled-cluster CCSD(T) energies with large basis sets and vibrational and relativistic corrections for 160 3d, 4d, and 5d metal-O2 systems. We define four reaction energies (120 in total for the 30 metals) that quantify O-O activation and reveal linear relationships between metal-oxygen and O-O binding energies. The CCSD(T) data can be combined with thermochemical cycles to estimate chemisorption and physisorption energies for each metal from metal oxide embedding energies, in good correlation with atomization enthalpies (R2 =0.75). Spin-geometry variations can break the linearities, of interest to circumventing the Sabatier principle. Pt, Pd, Co, and Fe form a distinct group with the weakest O2 binding. R2 up to 0.84 between surface adsorption energies and our energies for MO2 systems indicate relevance also to real catalytic systems.
Collapse
Affiliation(s)
- Klaus A Moltved
- Technical University of Denmark DTU Chemistry, Building 206, 2800, Kgs. Lyngby, Denmark
| | - Kasper P Kepp
- Technical University of Denmark DTU Chemistry, Building 206, 2800, Kgs. Lyngby, Denmark
| |
Collapse
|
29
|
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
|
30
|
Zhang IY, Xu X. On the top rung of Jacob's ladder of density functional theory: Toward resolving the dilemma of
SIE
and
NCE. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1490] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Igor Ying Zhang
- Shanghai Key Laboratory of Molecular Catalysis and Innovation Materials, Collaborative Innovation Centre of Chemistry for Energy Materials, MOE Laboratory for Computational Physical Science, Department of Chemistry Fudan University Shanghai China
| | - Xin Xu
- Shanghai Key Laboratory of Molecular Catalysis and Innovation Materials, Collaborative Innovation Centre of Chemistry for Energy Materials, MOE Laboratory for Computational Physical Science, Department of Chemistry Fudan University Shanghai China
| |
Collapse
|
31
|
Armentrout PB, Chang YC, Ng CY. What is the Bond Dissociation Energy of the Vanadium Hydride Cation? J Phys Chem A 2020; 124:5306-5313. [DOI: 10.1021/acs.jpca.0c04517] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- P. B. Armentrout
- Department of Chemistry, University of Utah, Salt Lake City, Utah 84112, United States
| | - Yih-Chung Chang
- Department of Chemistry, University of California, Davis, California 95616, United States
| | - Cheuk-Yiu Ng
- Department of Chemistry, University of California, Davis, California 95616, United States
| |
Collapse
|
32
|
Rudshteyn B, Coskun D, Weber JL, Arthur EJ, Zhang S, Reichman DR, Friesner RA, Shee J. Predicting Ligand-Dissociation Energies of 3d Coordination Complexes with Auxiliary-Field Quantum Monte Carlo. J Chem Theory Comput 2020; 16:3041-3054. [DOI: 10.1021/acs.jctc.0c00070] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Affiliation(s)
- Benjamin Rudshteyn
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Dilek Coskun
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - John L. Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Evan J. Arthur
- Schrodinger Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, 162 5th Avenue, New York, New York 10010, United States
- Department of Physics, College of William and Mary, Williamsburg, Virginia 23187, United States
| | - David R. Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Richard A. Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - James Shee
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
33
|
Reference spaces for multireference coupled-cluster theory: the challenge of the CoH molecule. Theor Chem Acc 2020. [DOI: 10.1007/s00214-020-2584-1] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
34
|
Süß D, Huber SE, Mauracher A. On the impact of multi-reference character of small transition metal compounds on their bond dissociation energies. J Chem Phys 2020; 152:114104. [DOI: 10.1063/1.5143495] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Daniel Süß
- Institute of Ion Physics and Applied Physics, University of Innsbruck, Technikerstraße 25, 6020 Innsbruck, Austria
| | - Stefan E. Huber
- Institute of Ion Physics and Applied Physics, University of Innsbruck, Technikerstraße 25, 6020 Innsbruck, Austria
| | - Andreas Mauracher
- Institute of Ion Physics and Applied Physics, University of Innsbruck, Technikerstraße 25, 6020 Innsbruck, Austria
| |
Collapse
|
35
|
Moltved KA, Kepp KP. Performance of Density Functional Theory for Transition Metal Oxygen Bonds. Chemphyschem 2019; 20:3210-3220. [DOI: 10.1002/cphc.201900862] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2019] [Revised: 10/01/2019] [Indexed: 11/09/2022]
Affiliation(s)
- Klaus A. Moltved
- Technical University of DenmarkDTU Chemistry, Building 206, 2800 Kgs. Lyngby DK – Denmark
| | - Kasper P. Kepp
- Technical University of DenmarkDTU Chemistry, Building 206, 2800 Kgs. Lyngby DK – Denmark
| |
Collapse
|
36
|
Hait D, Tubman NM, Levine DS, Whaley KB, Head-Gordon M. What Levels of Coupled Cluster Theory Are Appropriate for Transition Metal Systems? A Study Using Near-Exact Quantum Chemical Values for 3d Transition Metal Binary Compounds. J Chem Theory Comput 2019; 15:5370-5385. [PMID: 31465217 DOI: 10.1021/acs.jctc.9b00674] [Citation(s) in RCA: 34] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Transition metal compounds are traditionally considered to be challenging for standard quantum chemistry approximations like coupled cluster (CC) theory, which are usually employed to validate lower level methods like density functional theory (DFT). To explore this issue, we present a database of bond dissociation energies (BDEs) for 74 spin states of 69 diatomic species containing a 3d transition metal atom and a main group element, in the moderately sized def2-SVP basis. The presented BDEs appear to have an (estimated) 3σ error less than 1 kJ/mol relative to the exact solutions to the nonrelativistic Born-Oppenheimer Hamiltonian. These benchmark values were used to assess the performance of a wide range of standard single reference CC models, as the results should be beneficial for understanding the limitations of these models for transition metal systems. We find that interactions between metals and monovalent ligands like hydride and fluoride are well described by CCSDT. Similarly, CCSDTQ appears to be adequate for bonds between metals and nominally divalent ligands like oxide and sulfide. However, interactions with polyvalent ligands like nitride and carbide are more challenging, with even CCSDTQ(P)Λ yielding errors on the scale of a few kJ/mol. We also find that many perturbative and iterative approximations to higher order terms either yield disappointing results or actually worsen the performance relative to the baseline low level CC method, indicating that complexity does not always guarantee accuracy.
Collapse
Affiliation(s)
- Diptarka Hait
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry , University of California , Berkeley , California 94720 , United States.,Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
| | - Norman M Tubman
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry , University of California , Berkeley , California 94720 , United States.,Quantum Artificial Intelligence Lab. (QuAIL), Exploration Technology Directorate , NASA Ames Research Center , Moffett Field , California 94035 , United States
| | - Daniel S Levine
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry , University of California , Berkeley , California 94720 , United States
| | - K Birgitta Whaley
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry , University of California , Berkeley , California 94720 , United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry , University of California , Berkeley , California 94720 , United States.,Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
| |
Collapse
|
37
|
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
|
38
|
Ludovicy J, Mood KH, Lüchow A. Full Wave Function Optimization with Quantum Monte Carlo—A Study of the Dissociation Energies of ZnO, FeO, FeH, and CrS. J Chem Theory Comput 2019; 15:5221-5229. [DOI: 10.1021/acs.jctc.9b00241] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jil Ludovicy
- Institute of Physical Chemistry, RWTH Aachen University, Landoltweg 2, 52062 Aachen, Germany
| | - Kaveh Haghighi Mood
- Institute of Physical Chemistry, RWTH Aachen University, Landoltweg 2, 52062 Aachen, Germany
| | - Arne Lüchow
- Institute of Physical Chemistry, RWTH Aachen University, Landoltweg 2, 52062 Aachen, Germany
| |
Collapse
|
39
|
Cox RM, Kafle A, Armentrout PB, Peterson KA. Bond energy of ThN+: A guided ion beam and quantum chemical investigation of the reactions of thorium cation with N2 and NO. J Chem Phys 2019; 151:034304. [DOI: 10.1063/1.5111534] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Richard M. Cox
- Department of Chemistry, University of Utah, Salt Lake City, Utah 84112-0850, USA
| | - Arjun Kafle
- Department of Chemistry, University of Utah, Salt Lake City, Utah 84112-0850, USA
| | - P. B. Armentrout
- Department of Chemistry, University of Utah, Salt Lake City, Utah 84112-0850, USA
| | - Kirk A. Peterson
- Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, USA
| |
Collapse
|
40
|
Chan B. The CUAGAU Set of Coupled-Cluster Reference Data for Small Copper, Silver, and Gold Compounds and Assessment of DFT Methods. J Phys Chem A 2019; 123:5781-5788. [PMID: 31241947 DOI: 10.1021/acs.jpca.9b03976] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We have obtained benchmark data for a set of small molecular systems of Cu, Ag, and Au using coupled-cluster methods. Using this collection of reference data (that we termed the CUAGAU set) for assessing DFT-type methods, we find the MN15-L nonhybrid DFT to be cost-effective for geometry optimization [mean absolute deviation (MAD) in bond length = 0.20 Å], with an accuracy that is comparable to that for the double-hybrid (DH) DFT method DSD-PBEP86 (MAD = 0.19 Å). For the computation of thermochemical properties, among "conventional" (non-MP2-based) DFT methods, the best performance is found for the global-hybrid meta-GGA functional MN15, with an MAD of 11.4 kJ mol-1. We also find the nonhybrid method B97M-rV to have a reasonable performance (MAD = 14.4 kJ mol-1), and it may serve as a cost-effective means for qualitative study. If we look beyond conventional functionals, we find DSD-PBEP86 (MAD = 7.3 kJ mol-1) to be more accurate than even MN15. Nonetheless, this level of accuracy is still not sufficient for quantitative studies. In this regard, high-level wave function methods such as composite procedures that are based on coupled cluster are still indispensable for obtaining reliable reference data for transition-metal species.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering , Nagasaki University , Bunkyo 1-14 , Nagasaki 852-8521 , Japan
| |
Collapse
|
41
|
Xu Y, Chang YC, Ng CY. Chemical Activation of a Deuterium Molecule by Collision with a Quantum Electronic State-Selected Vanadium Cation. J Phys Chem A 2019; 123:5937-5944. [DOI: 10.1021/acs.jpca.9b04644] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yuntao Xu
- 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
| | - Cheuk-Yiu Ng
- Department of Chemistry, University of California, Davis, California 95616, United States
| |
Collapse
|
42
|
Cheng Q, DeYonker NJ. Theoretical study of the low-lying electronic states of iron hydride cation. J Chem Phys 2019; 150:234304. [PMID: 31228893 DOI: 10.1063/1.5096519] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Both FeH and FeH+ are predicted to be abundant in cool stellar atmospheres and proposed to be molecular components of the gas phase interstellar medium (ISM). However, experimental and simulated data for both species are lacking, which have hindered astronomical detection. There are no published laboratory data for the spectroscopy of FeH+ in any frequency regime. It is also not established if FeH+ possesses salient multireference character, which would pose significant challenges for ab initio modeling of geometric and spectroscopic properties. With a set of high-level coupled cluster and multireference configuration interaction computations, a study of the electronic structure of the ground state and seven excited states of FeH+ was carried out. An X5Δi electronic ground state of FeH+ is found, in agreement with previous theoretical studies. Including corrections for spin-orbit coupling and anharmonic vibrational effects, the Ω = 3, ν = 0 spin ladder of the A5Πi electronic state lies 872 cm-1 higher in energy than the Ω = 4, ν = 0 spin ladder of the ground state. Combined with previous work in our laboratory, the ionization energy of FeH is computed to be 7.4851 eV. With modern multireference configuration interaction and coupled cluster methods, spectroscopic constants (re, Be, ωe, ωexe, αe, and D¯e) for several bound excited states (A5Πi, B 5Σi +, a 3Σr -, b3Φi, c3Πi, d3Δr, and 7Σ+) were characterized. This study will lead efforts to identify FeH+ in the ISM and help solve important remaining questions in quantifying metal-hydride bonding.
Collapse
Affiliation(s)
- Qianyi Cheng
- Department of Chemistry, University of Memphis, Memphis, Tennessee 38152, USA
| | - Nathan J DeYonker
- Department of Chemistry, University of Memphis, Memphis, Tennessee 38152, USA
| |
Collapse
|
43
|
Minenkova I, Sliznev VV, Cavallo L, Minenkov Y. Gas Phase Silver Thermochemistry from First Principles. Inorg Chem 2019; 58:7873-7885. [PMID: 31185536 DOI: 10.1021/acs.inorgchem.9b00556] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Domain-based local pair natural orbital coupled cluster approach with single, double, and perturbative triple excitations, DLPNO-CCSD(T), has been applied within a framework of a reduced version of the reaction-based Feller-Peterson-Dixon (FPD) scheme to predict gas phase heats of formation and absolute entropies of silver inorganic and organometallic compounds. First, we evaluated all existing experimental data currently limited by thermodynamic functions of 10 silver substances (AgH, AgF, AgBr, AgI, Ag2, Ag2S, Ag2Se, Ag2Te, AgCN, AgPO2). The mean average deviation between computed and experimental heats of formation was found to be 1.9 kcal/mol. Notably, all predicted heats of formation turned out to be within the error bounds of their experimental counterparts. Second, we predicted heats of formation and entropies for additional 90 silver species with no experimental data available, substantially enriching silver thermochemistry. Combination of gas phase heats of formation Δ Hf and entropies S° of AgNO2, AgSCN, Ag2SO4, and Ag2SeO4 obtained in this work, with respective solid-state information, resulted in accurate sublimation thermochemistry of these compounds. Complementation of predicted Δ Hf with heats of formation of some neutrals and positive ions produced 33 silver bond strengths of high reliability. Obtained thermochemical data are promising for developing the concepts of silver chemistry. In addition, derived heats of formation and bond dissociation enthalpies, due to their high diversity, are found to be relevant for testing and training of computational chemistry methods.
Collapse
Affiliation(s)
- Irina Minenkova
- Moscow Institute of Physics and Technology , Institutskiy Pereulok 9 , Dolgoprudny , Moscow Region 141700 , Russian Federation.,A.N. Frumkin Institute of Physical Chemistry and Electrochemistry , Russian Academy of Sciences , 31 Leninsky Prospect , Moscow , GSP-1, 119071 , Russian Federation
| | - Valery V Sliznev
- Research Institute for Thermodynamics and Kinetics of Chemical Processes , Ivanovo State University of Chemistry and Technology , 153460 Ivanovo , Russian Federation
| | - Luigi Cavallo
- KAUST Catalysis Center (KCC) , King Abdullah University of Science and Technology , Thuwal 23955-6900 , Saudi Arabia
| | - Yury Minenkov
- Moscow Institute of Physics and Technology , Institutskiy Pereulok 9 , Dolgoprudny , Moscow Region 141700 , Russian Federation
| |
Collapse
|
44
|
Chan B, Gill PMW, Kimura M. Assessment of DFT Methods for Transition Metals with the TMC151 Compilation of Data Sets and Comparison with Accuracies for Main-Group Chemistry. J Chem Theory Comput 2019; 15:3610-3622. [PMID: 31150234 DOI: 10.1021/acs.jctc.9b00239] [Citation(s) in RCA: 78] [Impact Index Per Article: 15.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
In the present study, we have gathered a collection (that we term TMC151) of accurate reference data for transition-metal reactions for the assessment of quantum chemistry methods. It comprises diatomic dissociation energies and reaction energies and barriers for prototypical transition-metal reactions. Our assessment of a diverse range of different types of DFT methods shows that the most accurate functionals include ωB97M-V, ωB97X-V, MN15, and B97M-rV. Notably, they have also been previously validated to be highly robust for main-group chemistry. Nevertheless, even these methods show substantially worse accuracies for transition metals than for main-group chemistry. For less accurate methods, there is not a good correlation between their accuracies for main-group and transition-metal chemistries. Thus, in the development of new DFT, it is important to assess the accuracies for both types of data. In this regard, we have formulated the TMC34 model for efficient assessment of the performance for transition metals, which complements our previously developed MG8 model for main-group chemistry. Together, they provide a cost-effective means for initial assessment of new methodologies.
Collapse
Affiliation(s)
- Bun Chan
- Graduate School of Engineering , Nagasaki University , Bunkyo 1-14 , Nagasaki 852-8521 , Japan
| | - Peter M W Gill
- Research School of Chemistry , Australian National University , Canberra , Australian Capital Territory 2601 , Australia
| | - Masanari Kimura
- Graduate School of Engineering , Nagasaki University , Bunkyo 1-14 , Nagasaki 852-8521 , Japan
| |
Collapse
|
45
|
Bross DH, Yu HG, Harding LB, Ruscic B. Active Thermochemical Tables: The Partition Function of Hydroxymethyl (CH2OH) Revisited. J Phys Chem A 2019; 123:4212-4231. [DOI: 10.1021/acs.jpca.9b02295] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- David H. Bross
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Hua-Gen Yu
- Division of Chemistry, Department of Energy and Photon Sciences, Brookhaven National Laboratory, Upton, New York 11973, United States
| | - Lawrence B. Harding
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Branko Ruscic
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
- Consortium for Advanced Science and Engineering, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
46
|
Feller D, Bross DH, Ruscic B. Enthalpy of Formation of C2H2O4 (Oxalic Acid) from High-Level Calculations and the Active Thermochemical Tables Approach. J Phys Chem A 2019; 123:3481-3496. [DOI: 10.1021/acs.jpca.8b12329] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- David Feller
- Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, United States
| | - David H. Bross
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
| | - Branko Ruscic
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Argonne, Illinois 60439, United States
- Consortium for Advanced Science and Engineering, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
47
|
Moltved KA, Kepp KP. The Metal Hydride Problem of Computational Chemistry: Origins and Consequences. J Phys Chem A 2019; 123:2888-2900. [PMID: 30884233 DOI: 10.1021/acs.jpca.9b02367] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Formation and breaking of metal-hydrogen bonds are central to many important catalytic processes such as transition-metal catalyzed ammonia synthesis, hydrogenation reactions, and water splitting, and thus, they require an adequate theoretical description. We studied a data set of all 30 M-H and 30 M+-H bonds of the 3d, 4d, and 5d transition series; 50 of these systems have experimentally known bond dissociation enthalpies (BDE). To probe both the limit of low and high coordination number, we also studied a data set of 19 ML nH complexes. The BDEs were computed using Hartree-Fock (HF), MP2, CCSD, CCSD(T), and 10 diverse density functionals including local, GGA, hybrid GGA, meta hybrid, range-separated, and double hybrids. Our ten most important findings are as follows: (1) HF fails completely to describe the metal hydrogen bond due to its lack of static correlation; (2) this makes post-HF methods such as MP2 and even CCSD(T) perform worse than many density functionals; (3) DFT requires much more HF exchange (∼35% on average) to describe the pure M-H bonds than to describe other metal ligand bonds (0-20%); (4) we design a test to determine if self-interaction error (SIE) is important by correlating DFT errors against a one-electron SIE metric; (5) we show that SIE correlates directly with the DFT errors and thus causes most of the problem; (6) HF-DFT cannot handle these systems because the HF method is too pathological already at the density level; (7) instead, we define and apply a simple metric of electronic abnormality as the difference in PBE energy computed at the self-consistent PBE0 and SVWN densities, and this metric gives appropriate spread and effectively captures density-derived errors; (8) the low electronegativity of the metal enforces a diffuse hydride-like electron density, which make the metal hydrides primary examples of many-electron systems exhibiting SIE already at equilibrium geometries; (9) in the coordinatively saturated ML nH systems, much less HF exchange is required; i.e., the HF exchange requirements vary drastically with coordination number. Accordingly, DFT is unbalanced for any catalytic process involving both M-H and M-L bonds and changing coordination numbers; (10) importantly, the range-separated and double-hybrid functionals CAM-B3LYP and B2PLYP alone perform well for both M-H and M-L systems and in both limits of low and high coordination number, and at least as well as CCSD(T). This lends hope to a balanced treatment of computational chemistry for all types of M-L bonds at variable coordination number, as required for real catalytic reactions.
Collapse
Affiliation(s)
- Klaus A Moltved
- DTU Chemistry , Technical University of Denmark , Building 206 , 2800 Kgs . Lyngby , DK , Denmark
| | - Kasper P Kepp
- DTU Chemistry , Technical University of Denmark , Building 206 , 2800 Kgs . Lyngby , DK , Denmark
| |
Collapse
|
48
|
On Achieving High Accuracy in Quantum Chemical Calculations of 3d Transition Metal-Containing Systems: A Comparison of Auxiliary-Field Quantum Monte Carlo with Coupled Cluster, Density Functional Theory, and Experiment for Diatomic Molecules. J Chem Theory Comput 2019; 15:2346-2358. [DOI: 10.1021/acs.jctc.9b00083] [Citation(s) in RCA: 47] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
|
49
|
Pollice R, Chen P. Origin of the Immiscibility of Alkanes and Perfluoroalkanes. J Am Chem Soc 2019; 141:3489-3506. [DOI: 10.1021/jacs.8b10745] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
|
50
|
Ma P, Chen H. Ligand-Dependent Multi-State Reactivity in Cobalt(III)-Catalyzed C–H Activations. ACS Catal 2019. [DOI: 10.1021/acscatal.8b04532] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Pengchen Ma
- Beijing National Laboratory for Molecular Sciences (BNLMS), CAS Key Laboratory of Photochemistry, CAS Research/Education Center for Excellence in Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, P. R. China
- University of Chinese Academy of Sciences, Beijing 100049, P. R. China
| | - Hui Chen
- Beijing National Laboratory for Molecular Sciences (BNLMS), CAS Key Laboratory of Photochemistry, CAS Research/Education Center for Excellence in Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, P. R. China
| |
Collapse
|