1
|
Aðalsteinsson HM, Bjornsson R. Ionization energies of metallocenes: a coupled cluster study of cobaltocene. Phys Chem Chem Phys 2023; 25:4570-4587. [PMID: 36723003 DOI: 10.1039/d2cp04715b] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Open-shell transition metal chemistry presents challenges to contemporary electronic structure methods, based on either density functional or wavefunction theory. While CCSD(T) is the well-trusted gold standard for maingroup thermochemistry, the accuracy and robustness of the method is less clear for open-shell transition metal chemistry, requiring benchmarking of CCSD(T)-based protocols against either higher-level theory or experiment. Ionization energies (IEs) of metallocenes provide an interesting test case with metallocenes being common redox reagents as well as playing roles as redox mediators and cocatalysts in redox catalysis. Using highly accurate ZEKE-MATI experimental measurements of gas phase adiabatic (5.3275 ± 0.0006 eV) and vertical (5.4424 ± 0.0006 eV) ionization energies of cobaltocene, we systematically assessed the accuracy of the local coupled-cluster method DLPNO-CCSD(T) with respect to geometry, reference determinant, basis set size and extrapolation schemes, PNO cut-off and extrapolation, local triples approximation, relativistic effects and core-valence correlation. We show that PNO errors are controllable via the recently introduced PNO extrapolation schemes and that the expensive iterative triples (T1) contribution can be made more manageable by calculating it as a smaller-basis/smaller PNO-cutoff correction. The reference determinant turns out to be a critical aspect in these calculations with the HF determinant resulting in large DLPNO-CCSD(T) errors, likely due to the qualitatively flawed molecular orbital spectrum. The BP86 functional on the other hand was found to provide reference orbitals giving small DLPNO-CCSD(T) errors, likely due to more realistic orbitals as suggested by the more consistent MO spectrum compared to HF. A protocol including complete basis set extrapolations with correlation-consistent basis sets, complete PNO space extrapolations, iterative triples- and core-valence correlation corrections was found to give errors of -0.07 eV and -0.03 eV for adiabatic- and vertical-IE of cobaltocene, respectively, giving close to chemical accuracy for both properties. A computationally efficient DLPNO-CCSD(T) protocol was devised and tested against adiabatic ionization energies of 6 different metallocenes (V, Cr, Mn, Fe, Co, Ni). For the other metallocenes, the iterative triples (T1) and PNO extrapolation contributions turn out to be even more important. The results give errors close to the experimental uncertainty, similar to recent auxiliary-field quantum Monte Carlo results. The quality of the reference determinant orbitals is identified as the main source of uncertainty in CCSD(T) calculations of metallocenes.
Collapse
Affiliation(s)
| | - Ragnar Bjornsson
- Science Institute, University of Iceland, 107 Reykjavik, Iceland.,Univ Grenoble Alpes, CNRS, CEA, IRIG, Laboratoire de Chimie et Biologie des Métaux, 17 Rue des Martyrs, F-38054 Grenoble Cedex, France.
| |
Collapse
|
2
|
Accuracy of electronic density calculated using an optimally tuned range-separated hybrid functional. Theor Chem Acc 2023. [DOI: 10.1007/s00214-022-02952-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
|
3
|
Benedek Z, Timár P, Szilvási T, Barcza G. Sensitivity of coupled cluster electronic properties on the reference determinant: Can Kohn-Sham orbitals be more beneficial than Hartree-Fock orbitals? J Comput Chem 2022; 43:2103-2120. [PMID: 36200186 DOI: 10.1002/jcc.26996] [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: 03/01/2022] [Revised: 06/28/2022] [Accepted: 06/30/2022] [Indexed: 11/09/2022]
Abstract
Coupled cluster calculations are traditionally performed over Hartree-Fock reference orbitals (HF-CC methodology). However, in the literature it has been repeatedly raised whether the use of a Kohn-Sham reference (KS-CC methodology) might result in improved performance relative to HF-CC. In the present study, we re-examine the relation of HF-CC and KS-CC methods by comparing the results of widely applied truncated CC calculations (CCSD, CCSD(T), CCSDT) to the limit of full configuration interaction (FCI), which serves as an undebatable reference point of accuracy. Based on a series of CC calculations on diatoms and transition metal complexes, we demonstrate that no systematic improvement of coupled cluster electronic energies, densities and chemical reaction energies is expected when changing from HF to a KS reference. Nevertheless, fortuitous error cancellations might occasionally result in illusory improvement compared to HF-CC. Altogether, the application of KS-CC is not advantageous over HF-CC, but it is also not unreasonable as the choice of reference has negligible influence on the results at sufficiently high CC levels. In addition, KS-CC can be a particularly useful alternative if difficulties are encountered in HF or HF-CC convergence. It is also notable that KS-CC results are found to be practically independent of the chosen density functional, which implies that almost any KS-CC method can be used in place of HF-CC.
Collapse
Affiliation(s)
- Zsolt Benedek
- Wigner Research Centre for Physics, Budapest, Hungary.,Department of Inorganic and Analytical Chemistry, Budapest University of Technology and Economics, Budapest, Hungary.,Department of Chemical and Biological Engineering, The University of Alabama, Tuscaloosa, Alabama, USA
| | - Paula Timár
- Wigner Research Centre for Physics, Budapest, Hungary.,Department of Inorganic and Analytical Chemistry, Budapest University of Technology and Economics, Budapest, Hungary
| | - Tibor Szilvási
- Department of Chemical and Biological Engineering, The University of Alabama, Tuscaloosa, Alabama, USA
| | - Gergely Barcza
- Wigner Research Centre for Physics, Budapest, Hungary.,Department of Chemical and Biological Engineering, The University of Alabama, Tuscaloosa, Alabama, USA.,Department of Physics of Complex Systems, ELTE Eötvös Loránd University, Budapest, Hungary
| |
Collapse
|
4
|
Windom ZW, Perera A, Bartlett RJ. Benchmarking isotropic hyperfine coupling constants using (QTP) DFT functionals and coupled cluster theory. J Chem Phys 2022; 156:094107. [DOI: 10.1063/5.0069928] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Significant effort has been devoted to benchmarking isotropic hyperfine coupling constants for both wavefunction and density-based approaches in recent years, as accurate theoretical predictions aid the fitting of experimental model Hamiltonians. However, literature examining the predictive quality of a Density Functional Theory (DFT) functional abiding by the Bartlett IP condition is absent. In an attempt to rectify this, we report isotropic hyperfine coupling constant predictions of 24 commonly used DFT functionals on a total of 56 radicals, with the intent of exploring the successes and failures of the Quantum Theory Project (QTP) line of DFT functionals (i.e., CAM-QTP00, CAM-QTP01, CAM-QTP02, and QTP17) for this property. Included in this benchmark study are both small and large organic radicals as well as transition metal complexes, all of which have been studied to some extent in prior work. Subsequent coupled-cluster singles and doubles (CCSD) and CCSD withperturbative triples [CCSD(T)] calculations on small and large organic radicals show modest improvement as compared to prior work and offer an additional avenue for evaluation of DFT functional performance. We find that the QTP17 and CAM-QTP00 functionals consistently underperform, despite being parameterized to satisfy an IP eigenvalue condition primarily focused on inner shell electrons. On the other hand, the CAM-QTP01 functional is the most accurate functional in both organic radical datasets. Furthermore, both CAM-QTP01 and CAM-QTP02 are the most accurate functionals tested on the transition metal dataset. A significant portion of functionals were found to have comparable errors (within 5–15 MHz), but the hybrid class of DFT functionals maintains a consistently optimal balance between accuracy and precision across all datasets.
Collapse
Affiliation(s)
- Zachary W. Windom
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| | - Ajith Perera
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| | - Rodney J. Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| |
Collapse
|
5
|
Landeros-Rivera B, Gallegos M, Munarriz J, Laplaza R, Contreras García J. New venues in electron density analysis. Phys Chem Chem Phys 2022; 24:21538-21548. [DOI: 10.1039/d2cp01517j] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We provide a comprehensive overview of the chemical information within the electron density: how to extract information, but also how to obtain and how to assess the quality of the...
Collapse
|
6
|
Brémond É, Tognetti V, Chermette H, Sancho-García JC, Joubert L, Adamo C. Electronic Energy and Local Property Errors at QTAIM Critical Points while Climbing Perdew's Ladder of Density-Functional Approximations. J Chem Theory Comput 2021; 18:293-308. [PMID: 34958205 DOI: 10.1021/acs.jctc.1c00981] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We investigate the relationships between electron-density and electronic-energy errors produced by modern exchange-correlation density-functional approximations belonging to all of the rungs of Perdew's ladder. To this aim, a panel of relevant (semi)local properties evaluated at critical points of the electron-density field (as defined within the framework of Bader's atoms-in-molecules theory) are computed on a large selection of molecular systems involved in thermodynamic, kinetic, and noncovalent interaction chemical databases using density functionals developed in a nonempirical and minimally and highly parametrized fashion. The comparison of their density- and energy-based performance, also discussed in terms of density-driven errors, casts light on the strengths and weaknesses of the most recent and efficient density-functional approximations.
Collapse
Affiliation(s)
- Éric Brémond
- Université de Paris, ITODYS, CNRS, F-75006 Paris, France
| | - Vincent Tognetti
- Normandy University, COBRA UMR 6014 and FR 3038, Université de Rouen INSA Rouen, CNRS, F-76821 Mont St Aignan, France
| | - Henry Chermette
- Université de Lyon, Institut des Sciences Analytiques, UMR 5280, CNRS, Université Lyon 1, 5 rue de la Doua, F-69100 Villeurbanne, France
| | | | - Laurent Joubert
- Normandy University, COBRA UMR 6014 and FR 3038, Université de Rouen INSA Rouen, CNRS, F-76821 Mont St Aignan, France
| | - Carlo Adamo
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences (i-CLeHS), UMR 8060, F-75005 Paris, France.,Institut Universitaire de France, 103 Boulevard Saint Michel, F-75005 Paris, France
| |
Collapse
|
7
|
Martín-Fernández C, Harvey JN. On the Use of Normalized Metrics for Density Sensitivity Analysis in DFT. J Phys Chem A 2021; 125:4639-4652. [PMID: 34018759 DOI: 10.1021/acs.jpca.1c01290] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
In the past years, there has been a discussion about how the errors in density functional theory might be related to errors in the self-consistent densities obtained from different density functional approximations. This, in turn, brings up the discussion about the different ways in which we can measure such errors and develop metrics that assess the sensitivity of calculated energies to changes in the density. It is important to realize that there cannot be a unique metric in order to look at this density sensitivity, simultaneously needing size-extensive and size-intensive metrics. In this study, we report two metrics that are widely applicable to any density functional approximation. We also show how they can be used to classify different chemical systems of interest with respect to their sensitivity to small variations in the density.
Collapse
Affiliation(s)
| | - Jeremy N Harvey
- Department of Chemistry, KU Leuven, Celestijnenlaan, 200F 3001 Leuven, Belgium
| |
Collapse
|
8
|
Hait D, Liang YH, Head-Gordon M. Too big, too small, or just right? A benchmark assessment of density functional theory for predicting the spatial extent of the electron density of small chemical systems. J Chem Phys 2021; 154:074109. [DOI: 10.1063/5.0038694] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Affiliation(s)
- Diptarka Hait
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Yu Hsuan Liang
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
9
|
Lindquist-Kleissler B, Wenger JS, Johnstone TC. Analysis of Oxygen–Pnictogen Bonding with Full Bond Path Topological Analysis of the Electron Density. Inorg Chem 2021; 60:1846-1856. [DOI: 10.1021/acs.inorgchem.0c03308] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Affiliation(s)
- Brent Lindquist-Kleissler
- Department of Chemistry and Biochemistry, University of California, Santa Cruz, Santa Cruz, California 95064, United States
| | - John S. Wenger
- Department of Chemistry and Biochemistry, University of California, Santa Cruz, Santa Cruz, California 95064, United States
| | - Timothy C. Johnstone
- Department of Chemistry and Biochemistry, University of California, Santa Cruz, Santa Cruz, California 95064, United States
| |
Collapse
|
10
|
Jakobsen P, Jensen F. Representing Exact Electron Densities by a Single Slater Determinant in Finite Basis Sets. J Chem Theory Comput 2020; 17:269-276. [PMID: 33287541 DOI: 10.1021/acs.jctc.0c01029] [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/28/2022]
Abstract
The premise for Kohn-Sham density functional theory is that the exact electron density can be generated by a set of orbitals in a single Slater determinant. While this is ensured in a complete basis set, it has been shown that it cannot hold in small basis sets. The present work probes how accurately a reference electron density of the full-CI type can be reproduced by a set of orbitals in a single Slater determinant, as a function of the basis set used for the fitting electron density. The key finding is that the fitting error may be significant for basis sets of double- or triple-ζ quality. It is also shown that it is important that the fitting basis set includes the same basis functions as used for generating the reference electron density. The main limitation in a given basis set is the lack of higher order polarization functions. The error for practical purposes becomes insignificant for basis sets of quadruple-ζ or better quality, and this should be the choice when assessing the accuracy of exchange-correlation functionals by comparing electron densities to accurate reference results generated by wave function methods. The methodology in the present work can be used to transform an electron density from a multideterminant wave function into a set of orbitals in a single Slater determinant, and this may be useful for developing and testing new exchange-correlation functionals.
Collapse
Affiliation(s)
- Philip Jakobsen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus, Denmark
| | - Frank Jensen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus, Denmark
| |
Collapse
|
11
|
Witwicki M, Walencik PK, Jezierska J. How accurate is density functional theory in predicting spin density? An insight from the prediction of hyperfine coupling constants. J Mol Model 2019; 26:10. [PMID: 31834497 DOI: 10.1007/s00894-019-4268-0] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2019] [Accepted: 11/25/2019] [Indexed: 01/30/2023]
Abstract
Electron paramagnetic resonance (EPR) spectroscopy has been proven to be an important technique for studying paramagnetic systems. Probably, the most accessible EPR parameter and the one that provides a significant amount of information about molecular structure and spin density is the hyperfine coupling constant (HFCC). Hence, accurate quantum-chemical modeling of HFCCs is frequently essential to the adequate interpretation of EPR spectra. It requires the precise spin density, which is the difference between the densities of α- and β-electrons, and thus, its quality is expected to reflect the quality of the total electron density. The question of which approximate exchange-correlation density functional yields sufficiently accurate HFCCs, and thus, the spin density remains open. To assess the performance of well-established density functionals for calculating HFCCs, we used a series of 26 small paramagnetic species and compared the obtained results to the CCSD reference values. The performance of DFT was also tested on EPR-studied o-semiquinone radical interacting with water molecules and Mg2+ cation. The HFCCs were additionally calculated by the DLPNO-CCSD method, and this wave function-based technique was found superior to all functionals we tested. Although some functionals were found, on average, to be fairly efficient, we found that the most accurate functional is system-dependent, and therefore, the DLPNO-CCSD method should be preferred for theoretical investigations of the HFCCs and spin density.
Collapse
Affiliation(s)
- Maciej Witwicki
- Faculty of Chemistry, Wrocław University, ul. F. Joliot-Curie 14, 50-383, Wrocław, Poland.
| | - Paulina K Walencik
- Faculty of Chemistry, Wrocław University, ul. F. Joliot-Curie 14, 50-383, Wrocław, Poland
| | - Julia Jezierska
- Faculty of Chemistry, Wrocław University, ul. F. Joliot-Curie 14, 50-383, Wrocław, Poland
| |
Collapse
|
12
|
Abstract
The attraction density functional theory (DFT) has for electronic structure theory is that it is easier to do computationally than ab initio, correlated wavefunction methods, due to its effective one-particle structure. On the contrary, ab initio theorists insist on the ability to converge to the right answer in appropriate limits, but this requires a treatment of the reduced two-particle density matrix. DFT avoids that by appealing to an "existence" theorem (not a constructive one) that all its effects are subsummed into a DFT functional of the one-particle density. However, the existence of thousands of DFT functionals emphasizes that there is no satisfactory way to systematically improve the Kohn-Sham (KS) version as most changes in parameterization or formulation seldom lead to a new functional that is genuinely better than others. Some researchers in the DFT community try to address this issue by imposing conditions rigorously derived from exact DFT considerations, but to date, no one has shown how this route will ever lead to converged results even for the ground state, much less for all the other electronic states obtained from time-dependent DFT that are critically important for chemistry. On the contrary, coupled-cluster (CC) theory and its equation-of-motion extensions provide rigorous results for both that KS-DFT methods are attempting to emulate. How to use them and their exact formal properties to tie CC theory to an effective one-particle form is the target of this perspective. This route addresses the devil's triangle of KS-DFT problems: the one-particle spectrum, self-interaction, and the integer discontinuity.
Collapse
Affiliation(s)
- Rodney J Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611-8435, USA
| |
Collapse
|
13
|
Laplaza R, Polo V, Contreras-García J. Localizing electron density errors in density functional theory. Phys Chem Chem Phys 2019; 21:20927-20938. [PMID: 31517339 DOI: 10.1039/c9cp02806d] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The accuracy of different density functional approximations is assessed through the use of quantum chemical topology on molecular electron densities. In particular, three simple yet ever-important systems are studied: N2, CO and ethane. Our results exemplify how real-space descriptors can help understand the sources of errors in density functional theory, avoiding unwanted error compensation present in simplified statistical metrics. Errors in "well-built" functionals are shown to be concentrated in chemically meaningful regions of space, and hence they are predictable. Conversely, strongly parametrized functionals show isotropic errors that cannot be traced back to chemically transferable units. Moreover, we will show that energetic corrections are mapped back into improvements in the density in chemically meaningful regions. These results point at the relevance of real-space perspectives when parametrizing or relating energy and density errors.
Collapse
Affiliation(s)
- Rubén Laplaza
- Sorbonne Université, CNRS, Laboratoire de Chimie Thèorique, LCT, F. 75005 Paris, France and Departamento de Química Física and Instituto de Biocomputación y Física de Sistemas Complejos (BIFI), Universidad de Zaragoza, E. 50009, Zaragoza, Spain
| | - Victor Polo
- Departamento de Química Física and Instituto de Biocomputación y Física de Sistemas Complejos (BIFI), Universidad de Zaragoza, E. 50009, Zaragoza, Spain
| | - Julia Contreras-García
- Sorbonne Université, CNRS, Laboratoire de Chimie Thèorique, LCT, F. 75005 Paris, France.
| |
Collapse
|
14
|
Nikolaienko TY. The maximum occupancy condition for the localized property-optimized orbitals. Phys Chem Chem Phys 2019; 21:5285-5294. [PMID: 30778429 DOI: 10.1039/c8cp07276k] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
It is shown analytically that the Chemist's Localized Property-optimized Orbitals (CLPOs), which are the localized orbitals obtainable from the results of ab initio calculations by the open-source program JANPA (http://janpa.sourceforge.net/) according to the recently proposed optimal property partitioning condition, form the Lewis structure with nearly maximum possible total electron occupancy. The conditions required for this additional optimality to hold are discussed. In particular, when a single-determinant wavefunction is used to describe the molecular system without a noticeable electron delocalization, CLPOs derived from this wavefunction approximately optimize the same target quantity as the Natural Bond Orbitals (NBOs), establishing in this way the link between the two sets of localized orbitals. The performance of CLPO and NBO methods is compared by using a dataset containing 7101 small molecules, and the relevant methodological features of both methods are discussed.
Collapse
Affiliation(s)
- Tymofii Yu Nikolaienko
- Faculty of Physics of Taras Shevchenko National University of Kyiv, 64/13 Volodymyrska Str., Kyiv 01601, Ukraine.
| |
Collapse
|
15
|
Mostafanejad M, Haney J, DePrince AE. Kinetic-energy-based error quantification in Kohn–Sham density functional theory. Phys Chem Chem Phys 2019; 21:26492-26501. [DOI: 10.1039/c9cp04595c] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We present a basis-independent metric to assess the quality of the electron density obtained from Kohn–Sham (KS) density functional theory (DFT).
Collapse
Affiliation(s)
| | - Jessica Haney
- Department of Chemistry and Biochemistry
- Florida State University
- Tallahassee
- USA
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry
- Florida State University
- Tallahassee
- USA
| |
Collapse
|
16
|
Bouchafra Y, Shee A, Réal F, Vallet V, Severo Pereira Gomes A. Predictive Simulations of Ionization Energies of Solvated Halide Ions with Relativistic Embedded Equation of Motion Coupled Cluster Theory. PHYSICAL REVIEW LETTERS 2018; 121:266001. [PMID: 30636145 DOI: 10.1103/physrevlett.121.266001] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/06/2018] [Indexed: 06/09/2023]
Abstract
A subsystem approach for obtaining electron binding energies in the valence region is presented and applied to the case of halide ions (X^{-},X=F-At) in water. This approach is based on electronic structure calculations combining the relativistic equation-of-motion coupled cluster method for electron detachment and density functional theory via the frozen density embedding approach, using structures from classical molecular dynamics with polarizable force fields for discrete systems (in our study, droplets containing the anion and 50 water molecules). Our results indicate that one can accurately capture both the large solvent effect observed for the halides and the splitting of their ionization signals due to the increasingly large spin-orbit coupling of the p_{3/2}-p_{1/2} manifold across the series, at an affordable computational cost. Furthermore, owing to the quantum mechanical treatment of both solute and solvent electron binding energies of semiquantitative quality are also obtained for (bulk) water as by-products of the calculations for the halogens (in droplets).
Collapse
Affiliation(s)
- Yassine Bouchafra
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, F-59000 Lille, France
| | - Avijit Shee
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, F-59000 Lille, France
| | - Florent Réal
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, F-59000 Lille, France
| | - Valérie Vallet
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, F-59000 Lille, France
| | - André Severo Pereira Gomes
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, F-59000 Lille, France
| |
Collapse
|
17
|
Lutz JJ, Duan XF, Ranasinghe DS, Jin Y, Margraf JT, Perera A, Burggraf LW, Bartlett RJ. Valence and charge-transfer optical properties for some Si nC m( m, n≤ 12) clusters: Comparing TD-DFT, complete-basis-limit EOMCC, and benchmarks from spectroscopy. J Chem Phys 2018; 148:174309. [DOI: 10.1063/1.5022701] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Jesse J. Lutz
- Air Force Institute of Technology, Wright-Patterson Air Force Base, Ohio 45433, USA
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Xiaofeng F. Duan
- Air Force Institute of Technology, Wright-Patterson Air Force Base, Ohio 45433, USA
- Air Force Research Laboratory DoD Supercomputing Resource Center, Wright-Patterson Air Force Base, Ohio 45433, USA
| | | | - Yifan Jin
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Johannes T. Margraf
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Ajith Perera
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Larry W. Burggraf
- Air Force Institute of Technology, Wright-Patterson Air Force Base, Ohio 45433, USA
| | - Rodney J. Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
18
|
Interplay between test sets and statistical procedures in ranking DFT methods: The case of electron density studies. MENDELEEV COMMUNICATIONS 2018. [DOI: 10.1016/j.mencom.2018.05.001] [Citation(s) in RCA: 33] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
19
|
Hait D, Head-Gordon M. How Accurate Is Density Functional Theory at Predicting Dipole Moments? An Assessment Using a New Database of 200 Benchmark Values. J Chem Theory Comput 2018; 14:1969-1981. [DOI: 10.1021/acs.jctc.7b01252] [Citation(s) in RCA: 126] [Impact Index Per Article: 21.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Diptarka Hait
- 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
|
20
|
Costa R, Reta D, Moreira IDPR, Illas F. Post-B3LYP Functionals Do Not Improve the Description of Magnetic Coupling in Cu(II) Dinuclear Complexes. J Phys Chem A 2018; 122:3423-3432. [DOI: 10.1021/acs.jpca.7b12663] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
| | - Daniel Reta
- School of Chemistry, The University of Manchester, Oxford Road, Manchester M13 9PL, United Kingdom
| | | | | |
Collapse
|