1
|
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
|
2
|
Feller D, Stanton JF, Davidson ER. Atomic isotropic hyperfine properties for first row elements (B–F) revisited. J Chem Phys 2022; 156:034304. [DOI: 10.1063/5.0080155] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- David Feller
- Department of Chemistry, Washington State University, Pullman, Washington 99164-4630, USA and University of Alabama, Tuscaloosa, Alabama 35487-0336, USA
| | - John F. Stanton
- Quantum Theory Project, Department of Chemistry, University of Florida, Gainesville, Florida 32611, USA
| | - Ernest R. Davidson
- Department of Chemistry, University of Washington, Seattle, Washington 98195-1700, USA
| |
Collapse
|
3
|
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
|
4
|
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
|
5
|
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: 6] [Impact Index Per Article: 1.2] [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
|
6
|
Feller D. Estimating the intrinsic limit of the Feller-Peterson-Dixon composite approach when applied to adiabatic ionization potentials in atoms and small molecules. J Chem Phys 2017; 147:034103. [DOI: 10.1063/1.4993625] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
7
|
Spin-Adapted Formulation and Implementation of Density Cumulant Functional Theory with Density-Fitting Approximation: Application to Transition Metal Compounds. J Chem Theory Comput 2016; 12:4833-4842. [PMID: 27606799 DOI: 10.1021/acs.jctc.6b00589] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Density cumulant functional theory (DCT) has recently emerged as an attractive ab initio approach for the treatment of electron correlation. In its orbital-optimized formulation (ODC-12) [J. Chem. Phys. 139, 204110 (2013)], DCT has been shown to provide reliable results for a variety of challenging chemical systems. Among the attractive properties of DCT are its size-consistency and size-extensivity, as well as the efficient computation of the molecular properties and analytic gradients. In this work, we present a new formulation and implementation of DCT that takes advantage of spin adaptation and the density-fitting approximation (DF-ODC-12). Our new spin-adapted DF-ODC-12 implementation is more efficient than the previous ODC-12 implementation with up to a ∼12-fold speed-up. We demonstrate the capabilities of DF-ODC-12 with a study of transition metal compounds, which require high levels of electron correlation treatment. For transition metal carbonyl complexes [Fe(CO)5, Cr(CO)6] and the ferrocene molecule [Fe(Cp)2], the DF-ODC-12 equilibrium parameters and bond dissociation energies extrapolated to the complete basis set limit are in very good agreement with reference data derived from experiment.
Collapse
|
8
|
McKemmish LK, Yurchenko SN, Tennyson J. Ab initio calculations to support accurate modelling of the rovibronic spectroscopy calculations of vanadium monoxide (VO). Mol Phys 2016. [DOI: 10.1080/00268976.2016.1225994] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
- Laura K. McKemmish
- Department of Physics and Astronomy, University College London, London, UK
| | | | - Jonathan Tennyson
- Department of Physics and Astronomy, University College London, London, UK
| |
Collapse
|
9
|
The Impact of Larger Basis Sets and Explicitly Correlated Coupled Cluster Theory on the Feller–Peterson–Dixon Composite Method. ANNUAL REPORTS IN COMPUTATIONAL CHEMISTRY 2016. [DOI: 10.1016/bs.arcc.2016.02.001] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
|
10
|
DeYonker NJ. What a Difference a Decade Has Not Made: The Murky Electronic Structure of Iron Monocyanide (FeCN) and Iron Monoisocyanide (FeNC). J Phys Chem A 2014; 119:215-23. [DOI: 10.1021/jp5110906] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Nathan J. DeYonker
- Department of Chemistry, The University of Memphis, Memphis, Tennessee 38152, United States
| |
Collapse
|