1
|
Romero S, Baruah T, Zope RR. Spin-state gaps and self-interaction-corrected density functional approximations: Octahedral Fe(II) complexes as case study. J Chem Phys 2023; 158:054305. [PMID: 36754787 DOI: 10.1063/5.0133999] [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/18/2023] Open
Abstract
Accurate prediction of a spin-state energy difference is crucial for understanding the spin crossover phenomena and is very challenging for density functional approximations, especially for local and semi-local approximations due to delocalization errors. Here, we investigate the effect of the self-interaction error removal from the local spin density approximation (LSDA) and Perdew-Burke-Ernzerhof generalized gradient approximation on the spin-state gaps of Fe(II) complexes with various ligands using recently developed locally scaled self-interaction correction (LSIC) by Zope et al. [J. Chem. Phys. 151, 214108 (2019)]. The LSIC method is exact for one-electron density, recovers the uniform electron gas limit of the underlying functional, and approaches the well-known Perdew-Zunger self-interaction correction (PZSIC) as a particular case when the scaling factor is set to unity. Our results, when compared with reference diffusion Monte Carlo results, show that the PZSIC method significantly overestimates spin-state gaps favoring low spin states for all ligands and does not improve upon density functional approximations. The perturbative LSIC-LSDA using PZSIC densities significantly improves the gaps with a mean absolute error of 0.51 eV but slightly overcorrects for the stronger CO ligands. The quasi-self-consistent LSIC-LSDA, such as coupled-cluster single double and perturbative triple [CCSD(T)], gives a correct sign of spin-state gaps for all ligands with a mean absolute error of 0.56 eV, comparable to that of CCSD(T) (0.49 eV).
Collapse
Affiliation(s)
- Selim Romero
- Computational Science Program, The University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Tunna Baruah
- Department of Physics, University of Texas at El Paso, El Paso, Texas 79968, USA
| | - Rajendra R Zope
- Department of Physics, University of Texas at El Paso, El Paso, Texas 79968, USA
| |
Collapse
|
2
|
Louis H, Ekereke EE, Isang BB, Ikeuba AI, Amodu IO, Gber TE, Owen AE, Adeyinka AS, Agwamba EC. Assessing the Performance of Al 12N 12 and Al 12P 12 Nanostructured Materials for Alkali Metal Ion (Li, Na, K) Batteries. ACS OMEGA 2022; 7:46183-46202. [PMID: 36570229 PMCID: PMC9773795 DOI: 10.1021/acsomega.2c04319] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/08/2022] [Accepted: 09/26/2022] [Indexed: 06/17/2023]
Abstract
This study focused on the potential of aluminum nitride (Al12N12) and aluminum phosphide (Al12P12) nanomaterials as anode electrodes of lithium-ion (Li-ion), sodium-ion (Na-ion), and potassium-ion (K-ion) batteries as investigated via density functional theory (DFT) calculations at PBE0-D3, M062X-D3, and DSDPBEP86 as the reference method. The results show that the Li-ion battery has a higher cell voltage with a binding energy of -1.210 eV and higher reduction potential of -6.791 kcal/mol compared to the sodium and potassium ion batteries with binding energies of -0.749 and -0.935 eV and reduction potentials of -6.414 and -6.513 kcal/mol, respectively, using Al12N12 material. However, in Al12P12, increases in the binding energy and reduction potential were observed in the K-ion battery with values -1.485 eV and -7.535 kcal/mol higher than the Li and Na ion batteries with binding energy and reduction potential -1.483, -1.311 eV and -7.071, -7.184 eV, respectively. Finally, Al12N12 and Al12P12 were both proposed as novel anode electrodes in Li-ion and K-ion batteries with the highest performances.
Collapse
Affiliation(s)
- Hitler Louis
- Computational
and Bio-Simulation Research Group, University
of Calabar, Calabar540221, Nigeria
- Department
of Pure and Applied Chemistry, Faculty of Physical Sciences, University of Calabar, Calabar540221, Nigeria
| | - Ernest E. Ekereke
- Computational
and Bio-Simulation Research Group, University
of Calabar, Calabar540221, Nigeria
- Department
of Mathematics, Faculty of Physical Sciences, University of Calabar, Calabar540221, Nigeria
| | - Bartholomew B. Isang
- Computational
and Bio-Simulation Research Group, University
of Calabar, Calabar540221, Nigeria
- Department
of Mathematics, Faculty of Physical Sciences, University of Calabar, Calabar540221, Nigeria
| | - Alexander I. Ikeuba
- Department
of Pure and Applied Chemistry, Faculty of Physical Sciences, University of Calabar, Calabar540221, Nigeria
| | - Ismail O. Amodu
- Computational
and Bio-Simulation Research Group, University
of Calabar, Calabar540221, Nigeria
- Department
of Mathematics, Faculty of Physical Sciences, University of Calabar, Calabar540221, Nigeria
| | - Terkumbur E. Gber
- Computational
and Bio-Simulation Research Group, University
of Calabar, Calabar540221, Nigeria
- Department
of Pure and Applied Chemistry, Faculty of Physical Sciences, University of Calabar, Calabar540221, Nigeria
| | - Aniekan E. Owen
- Computational
and Bio-Simulation Research Group, University
of Calabar, Calabar540221, Nigeria
- School
of Chemistry, University of St Andrews, St AndrewsKY16 9ST, United Kingdom
| | - Adedapo S. Adeyinka
- Department
of Chemical Sciences, University of Johannesburg, Auckland Park2006South Africa
| | - Ernest C. Agwamba
- Computational
and Bio-Simulation Research Group, University
of Calabar, Calabar540221, Nigeria
- Department
of Chemistry, Covenant University, Ota50001, Nigeria
| |
Collapse
|
3
|
Reimann M, Kaupp M. Spin-State Splittings in 3d Transition-Metal Complexes Revisited: Benchmarking Approximate Methods for Adiabatic Spin-State Energy Differences in Fe(II) Complexes. J Chem Theory Comput 2022; 18:7442-7456. [PMID: 36417564 DOI: 10.1021/acs.jctc.2c00924] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
The CASPT2+δMRCI composite approach reported in a companion paper has been extended and used to provide high-quality reference data for a series of adiabatic spin gaps (defined as ΔE = Equintet - Esinglet) of [FeIIL6]2+ complexes (L = CNH, CO, NCH, NH3, H2O), either at nonrelativistic level or including scalar relativistic effects. These highly accurate data have been used to evaluate the performance of various more approximate methods. Coupled-cluster theory with singles, doubles, and perturbative triples, CCSD(T), is found to agree well with the new reference data for Werner-type complexes but exhibits larger underestimates by up to 70 kJ/mol for the π-acceptor ligands, due to appreciable static correlation in the low-spin states of these systems. Widely used domain-based local CCSD(T) calculations, DLPNO-CCSD(T), are shown to depend very sensitively on the cutoff values used to construct the localized domains, and standard values are not sufficient. A large number of density functional approximations have been evaluated against the new reference data. The B2PLYP double hybrid gives the smallest deviations, but several functionals from different rungs of the usual ladder hierarchy give mean absolute deviations below 20 kJ/mol. This includes the B97-D semilocal functional, the PBE0* global hybrid with 15% exact-exchange admixture, as well as the local hybrids LH07s-SVWN and LH07t-SVWN. Several further functionals achieve mean absolute errors below 30 kJ/mol (M06L-D4, SSB-D, B97-1-D4, LC-ωPBE-D4, LH12ct-SsirPW92-D4, LH12ct-SsifPW92-D4, LH14t-calPBE-D4, LHJ-HFcal-D4, and several further double hybrids) and thereby also still overall outperform CCSD(T) or uncorrected CASPT2. While exact-exchange admixture is a crucial factor in favoring high-spin states, the present evaluations confirm that other aspects can be important as well. A number of the better-performing functionals underestimate the spin gaps for the π-acceptor ligands but overestimate them for L = NH3, H2O. In contrast to a previous suggestion, non-self-consistent density functional theory (DFT) computations on top of Hartree-Fock orbitals are not a promising path to produce accurate spin gaps in such complexes.
Collapse
Affiliation(s)
- Marc Reimann
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| | - Martin Kaupp
- Technische Universität Berlin, Institut für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, D-10623 Berlin, Germany
| |
Collapse
|
4
|
Prokopiou G, Hartstein M, Govind N, Kronik L. Optimal Tuning Perspective of Range-Separated Double Hybrid Functionals. J Chem Theory Comput 2022; 18:2331-2340. [PMID: 35369687 PMCID: PMC9009176 DOI: 10.1021/acs.jctc.2c00082] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2022] [Indexed: 11/29/2022]
Abstract
We study the optimal tuning of the free parameters in range-separated double hybrid functionals, based on enforcing the exact conditions of piecewise linearity and spin constancy. We find that introducing the range separation in both the exchange and the correlation terms allows for the minimization of both fractional charge and fractional spin errors for singlet atoms. The optimal set of parameters is system specific, underlining the importance of the tuning procedure. We test the performance of the resulting optimally tuned functionals for the dissociation curves of diatomic molecules. We find that they recover the correct dissociation curve for the one-electron system, H2+, and improve the dissociation curves of many-electron molecules such as H2 and Li2, but they also yield a nonphysical maximum and only converge to the correct dissociation limit at very large distances.
Collapse
Affiliation(s)
- Georgia Prokopiou
- Department
of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth 76100, Israel
| | - Michal Hartstein
- Department
of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth 76100, Israel
| | - Niranjan Govind
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Leeor Kronik
- Department
of Molecular Chemistry and Materials Science, Weizmann Institute of Science, Rehovoth 76100, Israel
| |
Collapse
|
5
|
Alipour M, Izadkhast T. Do any types of double-hybrid models render the correct order of excited state energies in inverted singlet–triplet emitters? J Chem Phys 2022; 156:064302. [DOI: 10.1063/5.0077722] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Mojtaba Alipour
- Department of Chemistry, School of Science, Shiraz University, Shiraz 71946-84795, Iran
| | - Tahereh Izadkhast
- Department of Chemistry, School of Science, Shiraz University, Shiraz 71946-84795, Iran
| |
Collapse
|
6
|
Alipour M, Karimi N. Spin-Opposite-Scaled Range-Separated Exchange Double-Hybrid Models (SOS-RSX-DHs): Marriage Between DH and RSX/SOS-RSX Is Not Always a Happy Match. J Chem Theory Comput 2021; 17:4077-4091. [PMID: 34085815 DOI: 10.1021/acs.jctc.1c00271] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The range-separated version of double-hybrid density functional theory (DH-DFT) with a remarkable efficiency for both ground-state and excited-state characteristics has recently come into spotlight. In this work, based on theoretical arguments, several variants of spin-opposite-scaled range-separated exchange double-hybrid models (SOS-RSX-DHs) have been proposed and validated. More specifically, we first extend the RSX-DHs to design some other related models. Next, the SOS version of the resulting approximations is constructed and thoroughly evaluated using standard benchmark compilations of various properties. It is shown that although there are properties for which the RSX-DH and SOS-RSX-DH frameworks are rival, there are still some problems particularly prone to the self-interaction error issues where our proposed models seem to be beneficial. Furthermore, some of the presented models devoid of any additional corrections can also surpass the recently proposed approximations from different rungs of "Jacob's Ladder". Nonetheless, perusing the results of different methods and detailed comparisons with the predecessors discloses that all things may not necessarily be well with the RSX and SOS-RSX schemes, where the parent DHs as well as their SOS counterparts can still come into play.
Collapse
Affiliation(s)
- Mojtaba Alipour
- Department of Chemistry, School of Science, Shiraz University, Shiraz 71946-84795, Iran
| | - Niloofar Karimi
- Department of Chemistry, School of Science, Shiraz University, Shiraz 71946-84795, Iran
| |
Collapse
|
7
|
Mariano LA, Vlaisavljevich B, Poloni R. Improved Spin-State Energy Differences of Fe(II) Molecular and Crystalline Complexes via the Hubbard U-Corrected Density. J Chem Theory Comput 2021; 17:2807-2816. [PMID: 33831303 DOI: 10.1021/acs.jctc.1c00034] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We recently showed that the DFT+U approach with a linear-response U yields adiabatic energy differences biased toward high spin [Mariano et al. J. Chem. Theory Comput. 2020, 16, 6755-6762]. Such bias is removed here by employing a density-corrected DFT approach where the PBE functional is evaluated on the Hubbard U-corrected density. The adiabatic energy differences of six Fe(II) molecular complexes computed using this approach, named PBE[U] here, are in excellent agreement with coupled cluster-corrected CASPT2 values for both weak- and strong-field ligands resulting in a mean absolute error (MAE) of 0.44 eV, smaller than that of the recently proposed Hartree-Fock density-corrected DFT (1.22 eV) and any other tested functional, including the best performer TPSSh (0.49 eV). We take advantage of the computational efficiency of this approach and compute the adiabatic energy differences of five molecular crystals using PBE[U] with periodic boundary conditions. The results show, again, an excellent agreement (MAE = 0.07 eV) with experimentally extracted values and a superior performance compared with the best performers M06-L (MAE = 0.08 eV) and TPSSh (MAE = 0.31 eV) computed on molecular fragments.
Collapse
Affiliation(s)
- Lorenzo A Mariano
- University Grenoble Alpes, CNRS, Grenoble-INP, SIMaP, F-38042 Grenoble, France
| | - Bess Vlaisavljevich
- Department of Chemistry, University of South Dakota, Vermillion, South Dakota 57069, United States
| | - Roberta Poloni
- University Grenoble Alpes, CNRS, Grenoble-INP, SIMaP, F-38042 Grenoble, France
| |
Collapse
|
8
|
Farcaş AA, Bende A. Theoretical modeling of the singlet-triplet spin transition in different Ni(II)-diketo-pyrphyrin-based metal-ligand octahedral complexes. Phys Chem Chem Phys 2021; 23:4784-4795. [PMID: 33599640 DOI: 10.1039/d0cp05366j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The structural stability, charge transfer effects and strength of the spin-orbit couplings in different Ni(ii)-ligand complexes have been studied at the DFT (B3LYP and CAM-B3LYP) and coupled cluster (DLPNO-CCSD(T)) levels of theory. Accordingly, two different, porphyrin- and diketo-pyrphyrin-based four-coordination macrocycles as planar ligands as well as pyridine (or pyrrole) and mesylate anion molecular groups as vertical ligands were considered in order to build metal-organic complexes with octahedral coordination configurations. For each molecular system, the identification of equilibrium geometries and the intersystem crossing (the minimum energy crossing) points between the potential energy surfaces of the singlet and triplet spin states is followed by computing the spin-orbit couplings between the two spin states. Structures, based on the diketo-pyrphyrin macrocycle as the planar ligand, show stronger six-coordination metal-organic complexes due to the extra electrostatic interaction between the positively charged central metal cation and the negatively charged vertical ligands. The results also show that the magnitude of the spin-orbit coupling is influenced by the atomic positions of deprotonations of the ligands, and implicitly the direction of the charge transfer between the ligand and the central metal ion.
Collapse
Affiliation(s)
- Alex-Adrian Farcaş
- Faculty of Physics, "Babeş-Bolyai" University, Mihail Kogalniceanu Street No. 1, Ro-400084 Cluj-Napoca, Romania
| | - Attila Bende
- Molecular and Biomolecular Physics Department, National Institute for Research and Development of Isotopic and Molecular Technologies, Donat Street, No. 67-103, Ro-400293 Cluj-Napoca, Romania.
| |
Collapse
|
9
|
da Silva AR, de Almeida JS, Rivelino R. A Theoretical Assessment of Spin and Charge States in Binuclear Cobalt-Ruthenium Complexes: Implications for a Creutz-Taube Model Ion Separated by a C 60-Derivative Bridging Ligand. J Phys Chem A 2020; 124:10826-10837. [PMID: 33296201 DOI: 10.1021/acs.jpca.0c09194] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We investigate the spin-state energetics and the role of ionic charges in the electronic configuration of binuclear complexes of the form [(NH3)5Co(py)-X-(py)Ru(NH3)5]q+. In these compounds with q = 4-6, py = pyridine, and X = C≡C and C60, the Co-Ru distance varies from ∼1.4 to ∼2.1 nm. We carry out a systematic electronic structure calculation using different exchange-correlation (xc) approaches within spin-density functional theory, which are largely employed to investigate the properties of a variety of coordination complexes. To evaluate the effects of spin states and type of spacer in the bridging ligand on the valence tautomerism between Co2+/3+ and Ru2+/3+, we examine in more detail the case of Creutz-Taube-type ions [(NH3)5Co(py)-X-(py)Ru(NH3)5]5+. Our analysis shows that the stabilization of low- and high-spin states critically depends on the total charge of the complex, type of X-bridged ligand, and employed xc approach to calculate the electron spin density. Importantly, the C60-bridged group may result in a blockage of the valence tautomerism of the Creutz-Taube complex, inducing bistable charge configurations. Overall, our results also show that an adiabatic description in terms of the frontier molecular spin-orbitals for analyzing the distinct spin-charge states of these complexes may dramatically depend on the density-functional description.
Collapse
Affiliation(s)
- Alexsandro R da Silva
- Instituto de Física, Universidade Federal da Bahia, 40210-340 Salvador, Bahia, Brazil.,Instituto Federal do Maranhão, Campus São João dos Patos, 65665-000 São João dos Patos, Maranhão, Brazil
| | | | - Roberto Rivelino
- Instituto de Física, Universidade Federal da Bahia, 40210-340 Salvador, Bahia, Brazil
| |
Collapse
|
10
|
Casanova-Páez M, Goerigk L. Assessing the Tamm–Dancoff approximation, singlet–singlet, and singlet–triplet excitations with the latest long-range corrected double-hybrid density functionals. J Chem Phys 2020; 153:064106. [DOI: 10.1063/5.0018354] [Citation(s) in RCA: 36] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Affiliation(s)
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Victoria 3010, Australia
| |
Collapse
|