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
|
Jia L, Brémond É, Zaida L, Gaüzère B, Tognetti V, Joubert L. Predicting redox potentials by graph-based machine learning methods. J Comput Chem 2024. [PMID: 38923574 DOI: 10.1002/jcc.27380] [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: 11/08/2023] [Revised: 03/25/2024] [Accepted: 04/19/2024] [Indexed: 06/28/2024]
Abstract
The evaluation of oxidation and reduction potentials is a pivotal task in various chemical fields. However, their accurate prediction by theoretical computations, which is a complementary task and sometimes the only alternative to experimental measurement, may be often resource-intensive and time-consuming. This paper addresses this challenge through the application of machine learning techniques, with a particular focus on graph-based methods (such as graph edit distances, graph kernels, and graph neural networks) that are reviewed to enlighten their deep links with theoretical chemistry. To this aim, we establish the ORedOx159 database, a comprehensive, homogeneous (with reference values stemming from density functional theory calculations), and reliable resource containing 318 one-electron reduction and oxidation reactions and featuring 159 large organic compounds. Subsequently, we provide an instructive overview of the good practice in machine learning and of commonly utilized machine learning models. We then assess their predictive performances on the ORedOx159 dataset through extensive analyses. Our simulations using descriptors that are computed in an almost instantaneous way result in a notable improvement in prediction accuracy, with mean absolute error (MAE) values equal to 5.6 kcal mol- 1 $$ {}^{-1} $$ for reduction and 7.2 kcal mol- 1 $$ {}^{-1} $$ for oxidation potentials, which paves a way toward efficient in silico design of new electrochemical systems.
Collapse
Affiliation(s)
- Linlin Jia
- The PRG Group, Institute of Computer Science, University of Bern, Bern, Switzerland
| | - Éric Brémond
- Université Paris Cité, ITODYS, CNRS, Paris, France
| | | | - Benoit Gaüzère
- LITIS, Univ Rouen Normandie, INSA Rouen Normandie, Université Le Havre Normandie, Normandie Univ, Rouen, France
| | - Vincent Tognetti
- Normandy Univ., COBRA UMR 6014 & FR 3038, Université de Rouen, INSA Rouen, CNRS, Mont St Aignan Cedex, France
| | - Laurent Joubert
- Normandy Univ., COBRA UMR 6014 & FR 3038, Université de Rouen, INSA Rouen, CNRS, Mont St Aignan Cedex, France
| |
Collapse
|
3
|
Rodríguez-Mayorga M, Besalú-Sala P, Pérez-Jiménez ÁJ, Sancho-García JC. Application to nonlinear optical properties of the RSX-QIDH double-hybrid range-separated functional. J Comput Chem 2024; 45:995-1001. [PMID: 38206899 DOI: 10.1002/jcc.27302] [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: 11/10/2023] [Revised: 12/19/2023] [Accepted: 12/21/2023] [Indexed: 01/13/2024]
Abstract
The effective calculation of static nonlinear optical properties requires a considerably high accuracy at a reasonable computational cost, to tackle challenging organic and inorganic systems acting as precursors and/or active layers of materials in (nano-)devices. That trade-off implies to obtain very accurate electronic energies in the presence of externally applied electric fields to consequently obtain static polarizabilities (α i j ) and hyper-polarizabilities (β i j k andγ i j k l ). Density functional theory is known to provide an excellent compromise between accuracy and computational cost, which is however largely impeded for these properties without introducing range-separation techniques. We thus explore here the ability of a modern (double-hybrid and range-separated) Range-Separated eXchange Quadratic Integrand Double-Hybrid exchange-correlation functional to compete in accuracy with more costly and/or tuned methods, thanks to its robust and parameter-free nature.
Collapse
Affiliation(s)
- M Rodríguez-Mayorga
- Department of Physical Chemistry, University of Alicante, Alicante, Spain
- Université Grenoble Alpes, CNRS, Inst. NÉEL, Grenoble, France
| | - P Besalú-Sala
- Department of Chemistry and Pharmaceutical Sciences, Amsterdam Institute for Molecular and Life Sciences (AIMMS), Vrije Universiteit Amsterdam, Amsterdam, The Netherlands
| | - Á J Pérez-Jiménez
- Department of Physical Chemistry, University of Alicante, Alicante, Spain
| | - J C Sancho-García
- Department of Physical Chemistry, University of Alicante, Alicante, Spain
| |
Collapse
|
4
|
Knysh I, Raimbault D, Duchemin I, Blase X, Jacquemin D. Assessing the accuracy of TD-DFT excited-state geometries through optimal tuning with GW energy levels. J Chem Phys 2024; 160:144115. [PMID: 38602292 DOI: 10.1063/5.0203818] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2024] [Accepted: 03/25/2024] [Indexed: 04/12/2024] Open
Abstract
We study the accuracy of excited state (ES) geometries using optimally tuned LC-PBE functionals with tuning based on GW quasiparticle energies. We compare the results obtained with the PBE, PBE0, non-tuned, and tuned LC-PBE functionals with available high-level CC reference values as well as experimental data. First, we compare ES geometrical parameters obtained for three different types of systems: molecules composed of a few atoms, 4-(dimethylamino)benzonitrile (DMABN), and conjugated dyes. To this end, we used wave-function results as benchmarks. Next, we evaluate the accuracy of the theoretically simulated spectra as compared to the experimental ones for five large dyes. Our results show that, besides small compact molecules for which tuning LC-PBE does not allow obtaining geometries more accurate than those computed with standard functionals, tuned range-separated functionals are clearly to be favored, not only for ES geometries but also for 0-0 energies, band shapes, and intensities for absorption and emission spectra. In particular, the results indicate that GW-tuned LC-PBE functionals provide improved matching with experimental spectra as compared to conventionally tuned functionals. It is an open question whether TD-DFT with GW-tuned functionals can qualitatively mimic the actual many-body Bethe-Salpeter (BSE/GW) formalism for which analytic ionic gradients remain to be developed.
Collapse
Affiliation(s)
- Iryna Knysh
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Denez Raimbault
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Ivan Duchemin
- Université Grenoble Alpes, CEA, IRIG-MEM-L_Sim, 38054 Grenoble, France
| | - Xavier Blase
- Université Grenoble Alpes, CNRS, Institut, Néel F-38042, Grenoble
| | - Denis Jacquemin
- Nantes Université, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
- Institut Universitaire de France, 75005 Paris, France
| |
Collapse
|
5
|
Sburlati S, Gustafson A, Kahr B. Comparative rotatory power of bent and twisted polyynes. Chirality 2023; 35:838-845. [PMID: 37226985 DOI: 10.1002/chir.23579] [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/11/2023] [Revised: 04/29/2023] [Accepted: 05/03/2023] [Indexed: 05/26/2023]
Abstract
Linear polyynes of the formula C18 H2 (symmetry D∞h ) were bent in silico by progressively introducing CCC angles less than 180°. The bent structures (symmetry C2v ) were then twisted by introducing torsion angles across the CCCC segments by as much as 60°. The gyration tensors of these 19 structures (linear, bent, and twisted) were computed by linear response methods. Bending is massively generative of optical activity in oriented structures, even achiral structures, whereas twisting in conjunction with bending, serves to linearize the molecules and diminish maximally observable optical activity. This computational exercise is intended to unbind the infelicitous linkage of optical activity and chirality, which is only meaningful in isotropic media. Although bent structures are not optically active in solution-the spatial average of the optical activity is necessarily zero-solution measurements that deliver the spatial averages are a special class of measurements, albeit the overwhelmingly most common chiroptical measurements, that prejudice our common understanding of how π-conjugated structures generate gyration. Bending is far more effective than twisting at generating optical activity along some directions for oriented structures. The respective contributions from the transition electric dipole-magnetic dipole polarizability and the transition electric dipole-electric quadrupole polarizability are compared.
Collapse
Affiliation(s)
- Sophia Sburlati
- Department of Chemistry and Molecular Design Institute, New York University, New York City, New York, USA
| | - Afton Gustafson
- Department of Chemistry and Molecular Design Institute, New York University, New York City, New York, USA
| | - Bart Kahr
- Department of Chemistry and Molecular Design Institute, New York University, New York City, New York, USA
| |
Collapse
|
6
|
Tada T. Quantum Chemical Studies on Possible Molecular Devices Based on Electric Field-Induced Intramolecular Charge Transfer. J Phys Chem A 2023; 127:7297-7308. [PMID: 37638599 DOI: 10.1021/acs.jpca.3c02195] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/29/2023]
Abstract
We report quantum chemical studies on possible molecular devices working based on electric field-induced intramolecular charge transfer (EFIMCT). In the case of donor-acceptor (DA)-type molecular systems, intramolecular charge transfer (IMCT) can be induced by applying the external electric field to molecular systems along the charge transport direction, providing a possible switching mechanism which does not depend upon the electron-phonon coupling effect and is different from the negative differential resistance mechanism observed in the well-known NO2-substituted phenylene ethynylene oligomers. When the EFIMCT proceeds, the molecular systems have strong static electron correlation effects, where the standard nonequilibrium Green's function-density functional theory (DFT) approach cannot be applied to the molecular junction. As a first step toward practical switching devices, we do quantum chemical studies on the EFIMCT in such molecular systems as an isolated molecule, instead of using the electrode-junction-electrode open quantum system model. A prototype molecule P1 is designed as a tentative candidate molecule where the EFIMCT can proceed. The complete active space self-consistent field (CASSCF) molecular orbital calculations on P1 indicate that the EFIMCT can proceed at the external electric field intensity of 0.003 au, corresponding to about 2.25 V bias voltage. This calculated result strongly suggests that the development of this type of switching devices working at practically low bias voltage is feasible if the molecular system is properly designed. Broken symmetry unrestricted Hartree-Fock and spin-polarized Kohn-Sham DFT calculations also qualitatively reproduce the CASSCF results on P1, to some extent, indicating that these approaches can be employed for rough estimations on the EFIMCT such as the first screening of a large quantity of candidate molecules for this type of molecular devices. The possibility of molecular memory devices based on the EFIMCT is also discussed by analyzing the ground and excited potential energy surface model. Remaining challenges to develop practical molecular devices are discussed.
Collapse
Affiliation(s)
- Tsukasa Tada
- Department of Chemistry, Graduate School of Science, Tokyo Metropolitan University, Minami-Ohsawa 1-1, Hachioji, Tokyo 192-0397, Japan
| |
Collapse
|
7
|
Cohen RD, Wood JS, Lam YH, Buevich AV, Sherer EC, Reibarkh M, Williamson RT, Martin GE. DELTA50: A Highly Accurate Database of Experimental 1H and 13C NMR Chemical Shifts Applied to DFT Benchmarking. Molecules 2023; 28:molecules28062449. [PMID: 36985422 PMCID: PMC10051451 DOI: 10.3390/molecules28062449] [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: 01/25/2023] [Revised: 02/23/2023] [Accepted: 02/28/2023] [Indexed: 03/30/2023] Open
Abstract
Density functional theory (DFT) benchmark studies of 1H and 13C NMR chemical shifts often yield differing conclusions, likely due to non-optimal test molecules and non-standardized data acquisition. To address this issue, we carefully selected and measured 1H and 13C NMR chemical shifts for 50 structurally diverse small organic molecules containing atoms from only the first two rows of the periodic table. Our NMR dataset, DELTA50, was used to calculate linear scaling factors and to evaluate the accuracy of 73 density functionals, 40 basis sets, 3 solvent models, and 3 gauge-referencing schemes. The best performing DFT methodologies for 1H and 13C NMR chemical shift predictions were WP04/6-311++G(2d,p) and ωB97X-D/def2-SVP, respectively, when combined with the polarizable continuum solvent model (PCM) and gauge-independent atomic orbital (GIAO) method. Geometries should be optimized at the B3LYP-D3/6-311G(d,p) level including the PCM solvent model for the best accuracy. Predictions of 20 organic compounds and natural products from a separate probe set had root-mean-square deviations (RMSD) of 0.07 to 0.19 for 1H and 0.5 to 2.9 for 13C. Maximum deviations were less than 0.5 and 6.5 ppm for 1H and 13C, respectively.
Collapse
Affiliation(s)
- Ryan D Cohen
- Analytical Research and Development, Merck & Co., Inc., Rahway, NJ 07065, USA
- Department of Chemistry and Biochemistry, Seton Hall University, South Orange, NJ 07079, USA
| | - Jared S Wood
- Analytical Research and Development, Merck & Co., Inc., Rahway, NJ 07065, USA
- Department of Chemistry and Biochemistry, University of North Carolina Wilmington, Wilmington, NC 28409, USA
| | - Yu-Hong Lam
- Department of Computational and Structural Chemistry, Merck & Co., Inc., Rahway, NJ 07065, USA
| | - Alexei V Buevich
- Analytical Research and Development, Merck & Co., Inc., Rahway, NJ 07065, USA
| | - Edward C Sherer
- Analytical Research and Development, Merck & Co., Inc., Rahway, NJ 07065, USA
| | - Mikhail Reibarkh
- Analytical Research and Development, Merck & Co., Inc., Rahway, NJ 07065, USA
| | - R Thomas Williamson
- Department of Chemistry and Biochemistry, University of North Carolina Wilmington, Wilmington, NC 28409, USA
| | - Gary E Martin
- Department of Chemistry and Biochemistry, Seton Hall University, South Orange, NJ 07079, USA
| |
Collapse
|
8
|
Sandoval-Salinas ME, Brémond E, Pérez-Jiménez AJ, Adamo C, Sancho-García JC. Excitation energies of polycylic aromatic hydrocarbons by double-hybrid functionals: Assessing the PBE0-DH and PBE-QIDH models and their range-separated versions. J Chem Phys 2023; 158:044105. [PMID: 36725511 DOI: 10.1063/5.0134946] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023] Open
Abstract
A family of non-empirical double-hybrid (DH) density functionals, such as Perdew-Burke-Ernzerhof (PBE)0-DH, PBE-QIDH, and their range-separated exchange (RSX) versions RSX-0DH and RSX-QIDH, all using Perdew-Burke-Ernzerhof(PBE) exchange and correlationfunctionals, is applied here to calculate the excitation energies for increasingly longer linear and cyclic acenes as part of their intense benchmarking for excited states of all types. The energies for the two lowest-lying singlet 1La and 1Lb states of linear oligoacenes as well as the triplet 3La and 3Lb states, are calculated and compared with experimental results. These functionals clearly outperform the results obtained from hybrid functionals and favorably compare with other double-hybrid expressions also tested here, such as B2-PLYP, B2GP-PLYP, ωB2-PLYP, and ωB2GP-PLYP. The study is complemented by the computation of adiabatic S0-T1 singlet-triplet energy difference for linear acenes as well as the extension of the study to strained cyclic oligomers, showing how the family of non-empirical expressions robustly leads to competitive results.
Collapse
Affiliation(s)
- M E Sandoval-Salinas
- Department of Physical Chemistry, University of Alicante, E-03080 Alicante, Spain
| | - E Brémond
- ITODYS, CNRS, Université Paris Cité, F-75006 Paris, France
| | - A J Pérez-Jiménez
- Department of Physical Chemistry, University of Alicante, E-03080 Alicante, Spain
| | - C Adamo
- Chimie ParisTech, CNRS, Institute of Chemistry for Life and Health Sciences (i-CLeHS), UMR8060, PSL Research University, F-75005 Paris, France
| | - J C Sancho-García
- Department of Physical Chemistry, University of Alicante, E-03080 Alicante, Spain
| |
Collapse
|
9
|
de Oliveira MT, Alves JMA, Vrech NL, Braga AAC, Barboza CA. A comprehensive benchmark investigation of quantum chemical methods for carbocations. Phys Chem Chem Phys 2023; 25:1903-1922. [PMID: 36541431 DOI: 10.1039/d2cp04603b] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The application of various density functional approximations (DFAs) and an emphasis on popular methods without any consensus have prevailed in computational studies dedicated to carbocations. More importantly, an extensive and rigorous benchmark investigation on density functionals for the class is still lacking. To close this gap, we present a comprehensive benchmark study of quantum chemical methods on a series of classical and nonclassical carbocations, the CARBO33 dataset. We evaluate a total of 107 DFT methods from all rungs giving particular attention to double hybrid density functionals as the potential of the class has been largely undermined in the context of carbocations. To support our findings, DLPNO-CCSD(T) at the complete basis set (CBS) limit and W1-F12 are used as reference methods. Our results indicate that the composite CBS-QB3 method performs poorly and should not be adopted for target energies. Oftentimes, the tested DFAs of a lower rung perform better than several DFAs in a higher rung of Perdew's "Jacob's ladder". Nonetheless, double hybrids DSD-PBEP86-NL and ωB97X-2-D3(BJ) stand out by showing the overall best performance. Among the hybrids evaluated, about half of them show mean absolute deviation (MAD) below 1.1 kcal mol-1, including the popular hybrids M06-2X and mPW1PW91. In this family, MN15-D3(BJ) performs particularly well (MAD = 0.77 kcal mol-1) displaying reliable results across various tests. Highly popular B3LYP exhibited one of the worst performances (MAD = 4.74 kcal mol-1), and we do not recommend its application to carbocations. We also assess the 24 general-purpose basis sets of single- up to quadruple-ζ quality. The best compromise between accuracy and computational cost is achieved with cc-pVTZ followed by def2-TZVP. Computations on larger structures of general interest, including terpene carbocations, are also presented for selected DFT methods confirming general trends in the results.
Collapse
Affiliation(s)
- Marcelo T de Oliveira
- Department of Chemistry and Physics, La Trobe Institute of Molecular Sciences, La Trobe University, Melbourne, Victoria 3086, Australia. .,Chemistry Institute of São Carlos, University of São Paulo, Av. Trabalhador São Carlense 400, 13566-590, São Carlos, SP, Brazil
| | - Júlia M A Alves
- Chemistry Institute of São Carlos, University of São Paulo, Av. Trabalhador São Carlense 400, 13566-590, São Carlos, SP, Brazil
| | - Natália L Vrech
- Instituto de Química, Universidade de São Paulo, Av. Prof. Lineu Prestes 748, 05508-000, São Paulo, SP, Brazil
| | - Ataualpa A C Braga
- Instituto de Química, Universidade de São Paulo, Av. Prof. Lineu Prestes 748, 05508-000, São Paulo, SP, Brazil
| | - Cristina A Barboza
- Faculty of Chemistry, Wrocław University of Science and Technology, Wrocław, Poland.,Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, Warsaw, Poland
| |
Collapse
|
10
|
Helal W. Double Hybrid Density Functionals for the Electronic Excitation Energies of Linear Cyanines. J Phys Chem A 2023; 127:131-141. [PMID: 36537875 DOI: 10.1021/acs.jpca.2c07192] [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
The lowest bright electronic excitations of seven model linear cyanines (CN3-CN15) using 28 double-hybrid (DH) density functionals are benchmarked against accurate and recent CC3 results. Some of these DH functionals are recently designed specifically for excited electronic state calculations. In addition, CIS, CIS(D), SCS-CIS(D), and SOS-CIS(D) were also tested. Four different basis sets were used for the vertical electronic excitation calculations: cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis. Augmented basis sets (e.g. aug-cc-pVDZ and aug-cc-pVTZ) are found to be required for accurate and consistent results using DH functionals. The DH functionals tested in this work are classified into four main groups: global double-hybrids (GDH), range-separated double-hybrids (RSDH), spin-component and spin-opposite scaling global double-hybrids (SCS/SOS-GDH), and spin-component and spin-opposite scaling range-separated double-hybrids (SCS/SOS-RSDH). Within these groups, the SCS/SOS-RSDH group of functionals is found to provide the lowest mean absolute error (MAE) values (in the range 0.020-0.148 eV) in comparison with the GDH group (0.195-0.441 eV), the RSDH group (0.186-0.511 eV), and the SCS/SOS-GDH group (0.079-0.461 eV). Of all the DH functionals and ab initio methods investigated in the present contribution, the following functionals are found to be the most accurate and consistent: SCS-ωB2GPPLYP (MAE = 0.036 eV), SOS-ωB2GPPLYP (MAE = 0.020 eV), SOS-ωB88PP86 (MAE = 0.035 eV), and SOS-ωPBEPP86 (MAE = 0.037 eV). In general, the ab initio methods tested here show mediocre performance as compared to many DH functionals.
Collapse
Affiliation(s)
- Wissam Helal
- Department of Chemistry, The University of Jordan, Amman11942, Jordan
| |
Collapse
|
11
|
Spectroscopic, structural, and intermolecular interactions of 4-(2‑hydroxy-3-methoxybenzylideneamino)-N-(5-methylisoxazol-3-yl)benzenesulfonamide enol-imine and keto-amine isomers. J Mol Struct 2023. [DOI: 10.1016/j.molstruc.2023.134978] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
|
12
|
Shirani J, Farraj SA, Yuan S, Bevan KH. First-principles redox energy estimates under the condition of satisfying the general form of Koopmans’ theorem: An atomistic study of aqueous iron. J Chem Phys 2022; 157:184110. [DOI: 10.1063/5.0098476] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022] Open
Abstract
In this work, we explore the relative accuracy to which a hybrid functional, in the context of density functional theory, may predict redox properties under the constraint of satisfying the general form of Koopmans’ theorem. Taking aqueous iron as our model system within the framework of first-principles molecular dynamics, direct comparison between computed single-particle energies and experimental ionization data is assessed by both (1) tuning the degree of hybrid exchange, to satisfy the general form of Koopmans’ theorem, and (2) ensuring the application of finite-size corrections. These finite-size corrections are benchmarked through classical molecular dynamics calculations, extended to large atomic ensembles, for which good convergence is obtained in the large supercell limit. Our first-principles findings indicate that while precise quantitative agreement with experimental ionization data cannot always be attained for solvated systems, when satisfying the general form of Koopmans’ theorem via hybrid functionals, theoretically robust estimates of single-particle redox energies are most often arrived at by employing a total energy difference approach. That is, when seeking to employ a value of exact exchange that does not satisfy the general form of Koopmans’ theorem, but some other physical metric, the single-particle energy estimate that would most closely align with the general form of Koopmans’ theorem is obtained from a total energy difference approach. In this respect, these findings provide important guidance for the more general comparison of redox energies computed via hybrid functionals with experimental data.
Collapse
Affiliation(s)
- Javad Shirani
- Division of Materials Engineering, Faculty of Engineering, McGill University, Montréal, Québec H3A 0C5, Canada
| | - Sinan Abi Farraj
- Division of Materials Engineering, Faculty of Engineering, McGill University, Montréal, Québec H3A 0C5, Canada
| | - Shuaishuai Yuan
- Division of Materials Engineering, Faculty of Engineering, McGill University, Montréal, Québec H3A 0C5, Canada
| | - Kirk H. Bevan
- Division of Materials Engineering, Faculty of Engineering, McGill University, Montréal, Québec H3A 0C5, Canada
- Centre for the Physics of Materials, Department of Physics, McGill University, Montréal, Québec H3A 2T8, Canada
| |
Collapse
|
13
|
Bursch M, Mewes J, Hansen A, Grimme S. Best-Practice DFT Protocols for Basic Molecular Computational Chemistry. Angew Chem Int Ed Engl 2022; 61:e202205735. [PMID: 36103607 PMCID: PMC9826355 DOI: 10.1002/anie.202205735] [Citation(s) in RCA: 128] [Impact Index Per Article: 64.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2022] [Indexed: 01/11/2023]
Abstract
Nowadays, many chemical investigations are supported by routine calculations of molecular structures, reaction energies, barrier heights, and spectroscopic properties. The lion's share of these quantum-chemical calculations applies density functional theory (DFT) evaluated in atomic-orbital basis sets. This work provides best-practice guidance on the numerous methodological and technical aspects of DFT calculations in three parts: Firstly, we set the stage and introduce a step-by-step decision tree to choose a computational protocol that models the experiment as closely as possible. Secondly, we present a recommendation matrix to guide the choice of functional and basis set depending on the task at hand. A particular focus is on achieving an optimal balance between accuracy, robustness, and efficiency through multi-level approaches. Finally, we discuss selected representative examples to illustrate the recommended protocols and the effect of methodological choices.
Collapse
Affiliation(s)
- Markus Bursch
- Max-Planck-Institut für KohlenforschungKaiser-Wilhelm-Platz 145470Mülheim an der RuhrGermany
| | - Jan‐Michael Mewes
- Mulliken Center for Theoretical ChemistryInstitut für Physikalische und Theoretische ChemieUniversität BonnBeringstraße 453115BonnGermany
| | - Andreas Hansen
- Mulliken Center for Theoretical ChemistryInstitut für Physikalische und Theoretische ChemieUniversität BonnBeringstraße 453115BonnGermany
| | - Stefan Grimme
- Mulliken Center for Theoretical ChemistryInstitut für Physikalische und Theoretische ChemieUniversität BonnBeringstraße 453115BonnGermany
| |
Collapse
|
14
|
Bursch M, Mewes J, Hansen A, Grimme S. Best‐Practice DFT Protocols for Basic Molecular Computational Chemistry**. Angew Chem Int Ed Engl 2022. [DOI: 10.1002/ange.202205735] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Markus Bursch
- Max-Planck-Institut für Kohlenforschung Kaiser-Wilhelm-Platz 1 45470 Mülheim an der Ruhr Germany
| | - Jan‐Michael Mewes
- Mulliken Center for Theoretical Chemistry Institut für Physikalische und Theoretische Chemie Universität Bonn Beringstraße 4 53115 Bonn Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry Institut für Physikalische und Theoretische Chemie Universität Bonn Beringstraße 4 53115 Bonn Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry Institut für Physikalische und Theoretische Chemie Universität Bonn Beringstraße 4 53115 Bonn Germany
| |
Collapse
|
15
|
Santra G, Calinsky R, Martin JML. Benefits of Range-Separated Hybrid and Double-Hybrid Functionals for a Large and Diverse Data Set of Reaction Energies and Barrier Heights. J Phys Chem A 2022; 126:5492-5505. [PMID: 35930677 PMCID: PMC9393870 DOI: 10.1021/acs.jpca.2c03922] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
To better understand the thermochemical kinetics and
mechanism
of a specific chemical reaction, an accurate estimation of barrier
heights (forward and reverse) and reaction energies is vital. Because
of the large size of reactants and transition state structures involved
in real-life mechanistic studies (e.g., enzymatically catalyzed reactions),
density functional theory remains the workhorse for such calculations.
In this paper, we have assessed the performance of 91 density functionals
for modeling the reaction energies and barrier heights on a large
and chemically diverse data set (BH9) composed of 449 organic chemistry
reactions. We have shown that range-separated hybrid functionals perform
better than the global hybrids for BH9 barrier heights and reaction
energies. Except for the PBE-based range-separated nonempirical double
hybrids, range separation of the exchange term helps improve the performance
for barrier heights and reaction energies. The 16-parameter Berkeley
double hybrid, ωB97M(2), performs remarkably well for both properties.
However, our minimally empirical range-separated double hybrid functionals
offer marginally better accuracy than ωB97M(2) for BH9 barrier
heights and reaction energies.
Collapse
Affiliation(s)
- Golokesh Santra
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Rivka Calinsky
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| | - Jan M L Martin
- Department of Molecular Chemistry and Materials Science, Weizmann Institute of Science, 7610001 Reḥovot, Israel
| |
Collapse
|
16
|
Van Dijk J, Casanova-Páez M, Goerigk L. Assessing Recent Time-Dependent Double-Hybrid Density Functionals on Doublet-Doublet Excitations. ACS PHYSICAL CHEMISTRY AU 2022; 2:407-416. [PMID: 36855692 PMCID: PMC9955292 DOI: 10.1021/acsphyschemau.2c00014] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
This work is the first thorough investigation of time-dependent double-hybrid density functionals (DHDFs) for the calculation of doublet-doublet excitation energies. It sheds light on the current state-of-the-art techniques in the field and clarifies if there is still room for future improvements. Overall, 29 hybrid functionals and DHDFs are investigated. We separately analyze the individual impacts of the Tamm-Dancoff approximation (TDA), range separation, and spin-component/opposite scaling (SCS/SOS) on 45 doublet-doublet excitations in 23 radicals before concluding with an overarching analysis that includes and excludes challenging excitations with double-excitation or multireference character. Our results show again that so-called "nonempirical" DHDFs are outperformed by semiempirical ones. While the best assessed functionals are DHDFs, some of the worst are also DHDFs and outperformed by all assessed hybrids. SCS/SOS is particularly beneficial for range-separated DHDFs. Spin-scaled, range-separated DHDFs paired with the TDA belong to the best tested methods here, and we particularly highlight SCS-ωB2GP-PLYP, SOS-ωB2PLYP, SOS-ωB2GP-PLYP, SOS-ωB88PP86, SOS-RSX-QIDH, and SOS-ωPBEPP86. When comparing our functional rankings with previous studies on singlet-singlet and singlet-triplet excitations, we recommend TDA-SOS-ωB88PP86 and TDA-SOS-ωPBEPP86 as robust methods for excitation energies in general until further improvements have been achieved that surpass the chemical accuracy threshold for challenging open-shell excitations without increasing the computational effort.
Collapse
Affiliation(s)
- Joshua Van Dijk
- School
of Chemistry, The University of Melbourne, Victoria 3010, Australia
| | - Marcos Casanova-Páez
- School
of Chemistry, The University of Melbourne, Victoria 3010, Australia,Max-Planck-Institut
für Kohlenforschung, Kaiser-Wilhelm Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Lars Goerigk
- School
of Chemistry, The University of Melbourne, Victoria 3010, Australia,. Phone: +61 3 834 46784
| |
Collapse
|
17
|
Kilic E, Elmazoglu Z, Almammadov T, Kepil D, Etienne T, Marion A, Gunbas G, Kolemen S. Activity-Based Photosensitizers with Optimized Triplet State Characteristics Toward Cancer Cell Selective and Image Guided Photodynamic Therapy. ACS APPLIED BIO MATERIALS 2022; 5:2754-2767. [PMID: 35537187 PMCID: PMC9214761 DOI: 10.1021/acsabm.2c00202] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
![]()
Activity-based theranostic
photosensitizers are highly attractive
in photodynamic therapy as they offer enhanced therapeutic outcome
on cancer cells with an imaging opportunity at the same time. However,
photosensitizers (PS) cores that can be easily converted to activity-based
photosensitizers (aPSs) are still quite limited in the literature.
In this study, we modified the dicyanomethylene-4H-chromene (DCM) core with a heavy iodine atom to get two different
PSs (DCMO-I, I-DCMO-Cl) that can be further
converted to aPS after simple modifications. The effect of iodine
positioning on singlet oxygen generation capacity was also evaluated
through computational studies. DCMO-I showed better performance
in solution experiments and further proved to be a promising phototheranostic
scaffold via cell culture studies. Later, a cysteine (Cys) activatable
PS based on the DCMO-I core (DCMO-I-Cys) was
developed, which induced selective photocytotoxicity along with a
fluorescence turn-on response in Cys rich cancer cells.
Collapse
Affiliation(s)
- Eda Kilic
- Department of Chemistry, Koç University, 34450 Istanbul, Turkey
| | - Zubeyir Elmazoglu
- Department of Chemistry, Middle East Technical University (METU), 06800, Ankara, Turkey
| | | | - Dilay Kepil
- Department of Chemistry, Middle East Technical University (METU), 06800, Ankara, Turkey
| | | | - Antoine Marion
- Department of Chemistry, Middle East Technical University (METU), 06800, Ankara, Turkey
| | - Gorkem Gunbas
- Department of Chemistry, Middle East Technical University (METU), 06800, Ankara, Turkey
| | - Safacan Kolemen
- Department of Chemistry, Koç University, 34450 Istanbul, Turkey.,Surface Science and Technology Center (KUYTAM), Koç University, 34450 Istanbul, Turkey.,Boron and Advanced Materials Application and Research Center, Koç University, 34450 Istanbul, Turkey.,TUPRAS Energy Center (KUTEM), Koç University, 34450 Istanbul, Turkey
| |
Collapse
|
18
|
Brémond É, Li H, Pérez-Jiménez ÁJ, Sancho-García JC, Adamo C. Tackling an accurate description of molecular reactivity with double-hybrid density functionals. J Chem Phys 2022; 156:161101. [DOI: 10.1063/5.0087586] [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
In this Communication, we assess a panel of 18 double-hybrid density functionals for the modeling of the thermochemical and kinetic properties of an extended dataset of 449 organic chemistry reactions belonging to the BH9 database. We show that most of DHs provide a statistically robust performance to model barrier height and reaction energies in reaching the “chemical accuracy.” In particular, we show that nonempirical DHs, such as PBE0-DH and PBE-QIDH, or minimally parameterized alternatives, such as ωB2PLYP and B2K-PLYP, succeed to accurately model both properties in a balanced fashion. We demonstrate, however, that parameterized approaches, such as ωB97X-2 or DSD-like DHs, are more biased to only one of both properties.
Collapse
Affiliation(s)
- Éric Brémond
- ITODYS, CNRS, Université de Paris, F-75006 Paris, France
| | - Hanwei Li
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences (i-CLeHS), F-75005 Paris, France
| | | | | | - Carlo Adamo
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences (i-CLeHS), F-75005 Paris, France
- Institut Universitaire de France, 103 Boulevard Saint Michel, F-75005 Paris, France
| |
Collapse
|
19
|
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
|
20
|
Bursch M, Neugebauer H, Ehlert S, Grimme S. Dispersion corrected r 2SCAN based global hybrid functionals: r 2SCANh, r 2SCAN0, and r 2SCAN50. J Chem Phys 2022; 156:134105. [PMID: 35395897 DOI: 10.1063/5.0086040] [Citation(s) in RCA: 25] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
The regularized and restored semilocal meta-generalized gradient approximation (meta-GGA) exchange-correlation functional r2SCAN [Furness et al., J. Phys. Chem. Lett. 11, 8208-8215 (2020)] is used to create three global hybrid functionals with varying admixtures of Hartree-Fock "exact" exchange (HFX). The resulting functionals r2SCANh (10% HFX), r2SCAN0 (25% HFX), and r2SCAN50 (50% HFX) are combined with the semi-classical D4 London dispersion correction. The new functionals are assessed for the calculation of molecular geometries, main-group, and metalorganic thermochemistry at 26 comprehensive benchmark sets. These include the extensive GMTKN55 database, ROST61, and IONPI19 sets. It is shown that a moderate admixture of HFX leads to relative improvements of the mean absolute deviations for thermochemistry of 11% (r2SCANh-D4), 16% (r2SCAN0-D4), and 1% (r2SCAN50-D4) compared to the parental semi-local meta-GGA. For organometallic reaction energies and barriers, r2SCAN0-D4 yields an even larger mean improvement of 35%. The computation of structural parameters (geometry optimization) does not systematically profit from the HFX admixture. Overall, the best variant r2SCAN0-D4 performs well for both main-group and organometallic thermochemistry and is better or on par with well-established global hybrid functionals, such as PW6B95-D4 or PBE0-D4. Regarding systems prone to self-interaction errors (SIE4x4), r2SCAN0-D4 shows reasonable performance, reaching the quality of the range-separated ωB97X-V functional. Accordingly, r2SCAN0-D4 in combination with a sufficiently converged basis set [def2-QZVP(P)] represents a robust and reliable choice for general use in the calculation of thermochemical properties of both main-group and organometallic chemistry.
Collapse
Affiliation(s)
- Markus Bursch
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Hagen Neugebauer
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, D-53115 Bonn, Germany
| | - Sebastian Ehlert
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, D-53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, D-53115 Bonn, Germany
| |
Collapse
|
21
|
EPR Spectroscopy of Cu(II) Complexes: Prediction of g-Tensors Using Double-Hybrid Density Functional Theory. MAGNETOCHEMISTRY 2022. [DOI: 10.3390/magnetochemistry8040036] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Computational electron paramagnetic resonance (EPR) spectroscopy is an important field of applied quantum chemistry that contributes greatly to connecting spectroscopic observations with the fundamental description of electronic structure for open-shell molecules. However, not all EPR parameters can be predicted accurately and reliably for all chemical systems. Among transition metal ions, Cu(II) centers in inorganic chemistry and biology, and their associated EPR properties such as hyperfine coupling and g-tensors, pose exceptional difficulties for all levels of quantum chemistry. In the present work, we approach the problem of Cu(II) g-tensor calculations using double-hybrid density functional theory (DHDFT). Using a reference set of 18 structurally and spectroscopically characterized Cu(II) complexes, we evaluate a wide range of modern double-hybrid density functionals (DHDFs) that have not been applied previously to this problem. Our results suggest that the current generation of DHDFs consistently and systematically outperform other computational approaches. The B2GP-PLYP and PBE0-DH functionals are singled out as the best DHDFs on average for the prediction of Cu(II) g-tensors. The performance of the different functionals is discussed and suggestions are made for practical applications and future methodological developments.
Collapse
|
22
|
Mester D, Kállay M. Charge-Transfer Excitations within Density Functional Theory: How Accurate Are the Most Recommended Approaches? J Chem Theory Comput 2022; 18:1646-1662. [PMID: 35200021 PMCID: PMC8908740 DOI: 10.1021/acs.jctc.1c01307] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
The performance of
the most recent density functionals is assessed
for charge-transfer (CT) excitations using comprehensive intra- and
intermolecular CT benchmark sets with high-quality reference values.
For this comparison, the state-of-the-art range-separated (RS) and
long-range-corrected (LC) double hybrid (DH) approaches are selected,
and global DH and LC hybrid functionals are also inspected. The correct
long-range behavior of the exchange–correlation (XC) energy
is extensively studied, and various CT descriptors are compared as
well. Our results show that the most robust performance is attained
by RS-PBE-P86/SOS-ADC(2), as it is suitable to describe both types
of CT excitations with outstanding accuracy. Furthermore, concerning
the intramolecular transitions, unexpectedly excellent results are
obtained for most of the global DHs, but their limitations are also
demonstrated for bimolecular complexes. Despite the outstanding performance
of the LC-DH methods for common intramolecular excitations, serious
deficiencies are pointed out for intermolecular CT transitions, and
the wrong long-range behavior of the XC energy is revealed. The application
of LC hybrids to such transitions is not recommended in any respect.
Collapse
Affiliation(s)
- Dávid Mester
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| |
Collapse
|
23
|
Brémond E, Pérez-Jiménez AJ, Adamo C, Sancho-García JC. Stability of the polyynic form of C 18, C 22, C 26, and C 30 nanorings: a challenge tackled by range-separated double-hybrid density functionals. Phys Chem Chem Phys 2022; 24:4515-4525. [PMID: 35119058 DOI: 10.1039/d1cp04996h] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We calculate the relative energy between the cumulene and polyyne structures of a set of C4k+2 (k = 4-7) rings (C18, C22, C26, and C30 prompted by the recent synthesis of the cyclo[18]carbon (or simply C18) compounds. Reference results were obtained by a costly Quantum Monte-Carlo (QMC) approach, providing thus very accurate values allowing to systematically compare the performance of a variety of wavefunction methods [(i.e., MP2, SCS-MP2, SOS-MP2, DLPNO-CCSD, and DLPNO-CCSD(T)] as well as DFT approaches, applying for the latter a diversity of density functionals covering global and range-separated hybrid and double-hybrid models. The influence of the use of a range-separation scheme for density functionals, for both hybrid and double-hybrid expressions, is discussed according to its key role. Overall, range-separated double-hybrid functionals (e.g., RSX-QIDH) behave very accurately and provide competitive results compared with DLPNO-CCSD(T), at a more reasonable computational cost.
Collapse
Affiliation(s)
- E Brémond
- Université de Paris, ITODYS, CNRS, F-75006 Paris, France
| | - A J Pérez-Jiménez
- Department of Physical Chemistry, University of Alicante, E-03080 Alicante, Spain.
| | - C 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
| | - J C Sancho-García
- Department of Physical Chemistry, University of Alicante, E-03080 Alicante, Spain.
| |
Collapse
|
24
|
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
|
25
|
Flores-Holguín N, Ortega-Castro J, Frau J, Glossman-Mitnik D. Conceptual DFT-Based Computational Peptidology, Pharmacokinetics Study and ADMET Report of the Veraguamides A–G Family of Marine Natural Drugs. Mar Drugs 2022; 20:md20020097. [PMID: 35200627 PMCID: PMC8874632 DOI: 10.3390/md20020097] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2021] [Revised: 01/15/2022] [Accepted: 01/19/2022] [Indexed: 12/17/2022] Open
Abstract
As a continuation of our research on the chemical reactivity, pharmacokinetics and ADMET properties of cyclopeptides of marine origin with potential therapeutic abilities, in this work our already presented integrated molecular modeling protocol has been used for the study of the chemical reactivity and bioactivity properties of the Veraguamides A–G family of marine natural drugs. This protocol results from the estimation of the conceptual density functional theory (CDFT) chemical reactivity descriptors together with several chemoinformatics tools commonly considered within the process of development of new therapeutic drugs. CP-CDFT is a branch of computational chemistry and molecular modeling dedicated to the study of peptides, and it is a protocol that allows the estimation with great accuracy of the CDFT-based reactivity descriptors and the associated physical and chemical properties, which can aid in determining the ability of the studied peptides to behave as potential useful drugs. Moreover, the superiority of the MN12SX density functional over other long-range corrected density functionals for the prediction of chemical and physical properties in the presence of water as the solvent is clearly demonstrated. The research was supplemented with an investigation of the bioactivity of the molecular systems and their ADMET (absorption, distribution, metabolism, excretion, and toxicity) parameters, as is customary in medicinal chemistry. Some instances of the CDFT-based chemical reactivity descriptors’ capacity to predict the pKas of peptides as well as their potential as AGE inhibitors are also shown.
Collapse
Affiliation(s)
- Norma Flores-Holguín
- Laboratorio Virtual NANOCOSMOS, Departamento de Medio Ambiente y Energía, Centro de Investigación en Materiales Avanzados, Chihuahua 31136, Mexico;
| | - Joaquín Ortega-Castro
- Departament de Química, Facultat de Ciènces, Universitat de les Illes Balears, E-07122 Palma de Mallorca, Spain; (J.O.-C.); (J.F.)
| | - Juan Frau
- Departament de Química, Facultat de Ciènces, Universitat de les Illes Balears, E-07122 Palma de Mallorca, Spain; (J.O.-C.); (J.F.)
| | - Daniel Glossman-Mitnik
- Laboratorio Virtual NANOCOSMOS, Departamento de Medio Ambiente y Energía, Centro de Investigación en Materiales Avanzados, Chihuahua 31136, Mexico;
- Correspondence: ; Tel.: +52-614-439-1151
| |
Collapse
|
26
|
Sancho-García JC, Brémond E, Ricci G, Pérez-Jiménez AJ, Olivier Y, Adamo C. Violation of Hund’s rule in molecules: Predicting the excited-state energy inversion by TD-DFT with double-hybrid methods. J Chem Phys 2022; 156:034105. [DOI: 10.1063/5.0076545] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- J. C. Sancho-García
- Department of Physical Chemistry, University of Alicante, E-03080 Alicante, Spain
| | - E. Brémond
- Université de Paris, ITODYS, CNRS, F-75006 Paris, France
| | - G. Ricci
- Laboratory for Computational Modeling of Functional Materials, Namur Institute of Structured Matter, Université de Namur, Rue de Bruxelles, B-5000 Namur, Belgium
| | - A. J. Pérez-Jiménez
- Department of Physical Chemistry, University of Alicante, E-03080 Alicante, Spain
| | - Y. Olivier
- Laboratory for Computational Modeling of Functional Materials, Namur Institute of Structured Matter, Université de Namur, Rue de Bruxelles, B-5000 Namur, Belgium
| | - C. Adamo
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences (i-CLeHS), FRE 2027, F-75005 Paris, France
- Institut Universitaire de France, 103 Boulevard Saint Michel, F-75005 Paris, France
| |
Collapse
|
27
|
Wu R, Matta M, Paulsen BD, Rivnay J. Operando Characterization of Organic Mixed Ionic/Electronic Conducting Materials. Chem Rev 2022; 122:4493-4551. [PMID: 35026108 DOI: 10.1021/acs.chemrev.1c00597] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Operando characterization plays an important role in revealing the structure-property relationships of organic mixed ionic/electronic conductors (OMIECs), enabling the direct observation of dynamic changes during device operation and thus guiding the development of new materials. This review focuses on the application of different operando characterization techniques in the study of OMIECs, highlighting the time-dependent and bias-dependent structure, composition, and morphology information extracted from these techniques. We first illustrate the needs, requirements, and challenges of operando characterization then provide an overview of relevant experimental techniques, including spectroscopy, scattering, microbalance, microprobe, and electron microscopy. We also compare different in silico methods and discuss the interplay of these computational methods with experimental techniques. Finally, we provide an outlook on the future development of operando for OMIEC-based devices and look toward multimodal operando techniques for more comprehensive and accurate description of OMIECs.
Collapse
Affiliation(s)
- Ruiheng Wu
- Department of Chemistry, Northwestern University, Evanston, Illinois 60208, United States
| | - Micaela Matta
- Department of Chemistry, University of Liverpool, Liverpool L69 7ZD, United Kingdom
| | - Bryan D Paulsen
- Department of Biomedical Engineering, Northwestern University, Evanston, Illinois 60208, United States
| | - Jonathan Rivnay
- Department of Biomedical Engineering, Northwestern University, Evanston, Illinois 60208, United States.,Simpson Querrey Institute, Northwestern University, Chicago, Illinois 60611, United States
| |
Collapse
|
28
|
Mester D, Kállay M. Accurate Spectral Properties within Double-Hybrid Density Functional Theory: A Spin-Scaled Range-Separated Second-Order Algebraic-Diagrammatic Construction-Based Approach. J Chem Theory Comput 2022; 18:865-882. [PMID: 35023739 PMCID: PMC8830052 DOI: 10.1021/acs.jctc.1c01100] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Our second-order algebraic-diagrammatic construction [ADC(2)]-based double-hybrid (DH) ansatz (J. Chem. Theory Comput. 2019, 15, 4440. DOI: 10.1021/acs.jctc.9b00391) is combined with range-separation techniques. In the present scheme, both the exchange and the correlation contributions are range-separated, while spin-scaling approaches are also applied. The new methods are thoroughly tested for the most popular benchmark sets including 250 singlet and 156 triplet excitations, as well as 80 oscillator strengths. It is demonstrated that the range separation for the correlation contributions is highly recommended for both the genuine and the ADC(2)-based DH approaches. Our results show that the latter scheme slightly but consistently outperforms the former one for single excitation dominated transitions. Furthermore, states with larger fractions of double excitations are assessed as well, and challenging charge-transfer excitations are also discussed, where the recently proposed spin-scaled long-range corrected DHs fail. The suggested iterative fourth-power scaling RS-PBE-P86/SOS-ADC(2) method, using only three adjustable parameters, provides the most robust and accurate excitation energies within the DH theory. In addition, the relative error of the oscillator strengths is reduced by 65% compared to the best genuine DH functionals.
Collapse
Affiliation(s)
- Dávid Mester
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| |
Collapse
|
29
|
Flores-Holguín N, Frau J, Glossman-Mitnik D. Computational peptidology approach to the study of the chemical reactivity and bioactivity properties of Aspergillipeptide D, a cyclopentapeptide of marine origin. Sci Rep 2022; 12:506. [PMID: 35017576 PMCID: PMC8752680 DOI: 10.1038/s41598-021-04513-z] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2021] [Accepted: 12/24/2021] [Indexed: 12/19/2022] Open
Abstract
Aspergillipeptide D is a cyclic pentapeptide isolated from the marine gorgonian Melitodes squamata-derived fungus Aspergillus sp. SCSIO 41501 that it has been shown to present moderate activity against herpes virus simplex type 1 (HSV-1). Thus, this paper presents the results of a computational study of this cyclopentapeptide's chemical reactivity and bioactivity properties using a CDFT-based computational peptidology (CDFT-CP) methodology, which is derived from combining chemical reactivity descriptors derived from Conceptual Density Functional Theory (CDFT) and some Cheminformatics tools which may be used. This results in an improvement of the virtual screening procedure by a similarity search allowing the identification and validation of the known ability of the peptide to act as a possible useful drug. This was followed by an examination of the drug's bioactivity and pharmacokinetics indices in relation to the ADMET (Absorption, Distribution, Metabolism, Excretion, and Toxicity) characteristics. The findings provide further evidence of the MN12SX density functional's superiority in proving the Janak and Ionization Energy theorems using the proposed KID approach. This has proven to be beneficial in accurately predicting CDFT reactivity characteristics, which aid in the understanding of chemical reactivity. The Computational Pharmacokinetics study revealed the potential ability of Aspergillipeptide D as a therapeutic drug through the interaction with different target receptors. The ADMET indices confirm this assertion through the absence of toxicity and good absorption and distribution properties.
Collapse
Affiliation(s)
- Norma Flores-Holguín
- Laboratorio Virtual NANOCOSMOS, Departamento de Medio Ambiente y Energía, Centro de Investigación en Materiales Avanzados, 31136, Chihuahua, CHIH, Mexico
| | - Juan Frau
- Departament de Química, Universitat de les Illes Balears, Palma de Mallorca, 07122, Spain
| | - Daniel Glossman-Mitnik
- Laboratorio Virtual NANOCOSMOS, Departamento de Medio Ambiente y Energía, Centro de Investigación en Materiales Avanzados, 31136, Chihuahua, CHIH, Mexico.
| |
Collapse
|
30
|
Thor W, Zhang Y, Wong KL, Tanner PA. Orbital transitions: insight into energy transfer through an antenna for an organo-lanthanide complex. Chem Commun (Camb) 2021; 57:10727-10730. [PMID: 34585177 DOI: 10.1039/d1cc05246b] [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
Supported by experimental work, wavefunction theory (WFT) calculations and density functional theory (DFT) calculations employing a range of functionals have been performed for two lanthanide complexes to investigate, in gas and solution phases, the representations of frontier orbitals and the orbital transitions between singlet states. The orbital transitions calculated using CASSCF/NEVTP2 served as reference. Functionals with a higher proportion of Hartree-Fock exchange gave better agreement with WFT. The choice of functional is therefore important for understanding the nature of orbital transitions and this is especially relevant in formulating antenna-metal ion energy transfer (ET) mechanisms.
Collapse
Affiliation(s)
- Waygen Thor
- Department of Chemistry, Hong Kong Baptist University, Waterloo Road, Kowloon Tong, Hong Kong S.A.R, Urumqi 830046, Xinjiang, P. R. China.
| | - Yonghong Zhang
- Department of Chemistry, Hong Kong Baptist University, Waterloo Road, Kowloon Tong, Hong Kong S.A.R, Urumqi 830046, Xinjiang, P. R. China. .,State Key Laboratory of Chemistry and Utilization of Carbon Based Energy Resources, Key Laboratory of Oil and Gas Fine Chemicals, Ministry of Education & Xinjiang Uygur Autonomous Region, Urumqi Key Laboratory of Green Catalysis and Synthesis Technology, College of Chemistry, Xinjiang University, Urumqi 830046, Xinjiang, P. R. China
| | - Ka-Leung Wong
- Department of Chemistry, Hong Kong Baptist University, Waterloo Road, Kowloon Tong, Hong Kong S.A.R, Urumqi 830046, Xinjiang, P. R. China.
| | - Peter A Tanner
- Department of Chemistry, Hong Kong Baptist University, Waterloo Road, Kowloon Tong, Hong Kong S.A.R, Urumqi 830046, Xinjiang, P. R. China.
| |
Collapse
|
31
|
Robertson CM, Winter SM, Howard JAK, Probert MR, Oakley RT. Low temperature insights into the crystal and magnetic structure of a neutral radical ferromagnet. Chem Commun (Camb) 2021; 57:10238-10241. [PMID: 34528050 DOI: 10.1039/d1cc03842g] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The crystal structure of the radical ferromagnet 1a at 2 K reveals a contraction in the unit cell c constant which, at the molecular level, translates into a decrease in slippage of the radical π-stacks and an increase in ferromagnetic exchange interactions along the stacking axis. The results of BS-DFT calculations using long-range corrected functionals are consistent with an overall ferromagnetic topology.
Collapse
Affiliation(s)
- Craig M Robertson
- Department of Chemistry, University of Liverpool, Liverpool L69 7ZD, UK
| | - Stephen M Winter
- Department of Physics, Wake Forest University, Winston-Salem, NC 27109, USA
| | | | - Michael R Probert
- Chemistry, School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne, NE1 7RU, UK
| | - Richard T Oakley
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada.
| |
Collapse
|
32
|
Smith BA, Vogiatzis KD. σ-Donation and π-Backdonation Effects in Dative Bonds of Main-Group Elements. J Phys Chem A 2021; 125:7956-7966. [PMID: 34477393 DOI: 10.1021/acs.jpca.1c05956] [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
The nature of donor-acceptor interactions is important for the understanding of dative bonding and can provide vital insights into many chemical processes. Here, we have performed a computational study to elucidate substantial differences between different types of dative interactions. For this purpose, a data set of 20 molecular complexes stabilized by dative bonds was developed (DAT20). A benchmark study that considers many popular density functionals with respect to accurate quantum chemical interaction energies and geometries revealed two different trends between the complexes of DAT20. This behavior was further explored by means of frontier molecular orbitals, extended-transition-state natural orbitals for chemical valence (ETS-NOCV), and natural energy decomposition analysis (NEDA). These methods revealed the extent of the forward and backdonation between the donor and acceptor molecules and how they influence the total interaction energies and molecular geometries. A new classification of dative bonds is suggested.
Collapse
Affiliation(s)
- Brett A Smith
- Department of Chemistry, University of Tennessee, Knoxville, Tennessee 37996, United States
| | | |
Collapse
|
33
|
A CDFT-Based Computational Peptidology (CDFT-CP) Study of the Chemical Reactivity and Bioactivity of the Marine-Derived Alternaramide Cyclopentadepsipeptide. J CHEM-NY 2021. [DOI: 10.1155/2021/2989611] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
Abstract
Alternaramide is a cyclic pentadepsipeptide isolated from marine sources that has been shown to present weak antibiotic activity against Bacillus subtilis and Staphylococcus aureus as well as inhibitory effects on inflammatory mediator expressions. Thus, this work reports the results of a computational study of the chemical reactivity and bioactivity properties of this cyclopentadepsipeptide considering a CDFT-based computational peptidology (CDFT-CP) methodology that results from the combination of the chemical reactivity descriptors that arise from conceptual density functional theory (CDFT) together with some cheminformatics tools that can be used to estimate the associated physicochemical parameters, to improve the process of virtual screening through a similarity search, and to identify the ability of the peptide to behave as a potential useful drug, complemented with an analysis of its bioactivity and pharmacokinetics indices related to the ADMET (absorption, distribution, metabolism, excretion, and toxicity) features. The results represent a new confirmation of the superiority of the MN12SX density functional in the fulfilment of the Janak and ionization energy theorems through the proposed KID procedure. This has been useful for the accurate prediction of the CDFT reactivity descriptors that help in understanding the chemical reactivity. The computational pharmacokinetics study revealed the potential ability of alternaramide as a therapeutic drug by interacting with GPCR ligands and protease inhibitors. The ADMET indices confirm this assertion through the absence of toxicity and good absorption and distribution properties.
Collapse
|
34
|
Casanova-Páez M, Goerigk L. Time-Dependent Long-Range-Corrected Double-Hybrid Density Functionals with Spin-Component and Spin-Opposite Scaling: A Comprehensive Analysis of Singlet-Singlet and Singlet-Triplet Excitation Energies. J Chem Theory Comput 2021; 17:5165-5186. [PMID: 34291643 DOI: 10.1021/acs.jctc.1c00535] [Citation(s) in RCA: 40] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Following the work on spin-component and spin-opposite scaled (SCS/SOS) global double hybrids for singlet-singlet excitations by Schwabe and Goerigk [ J. Chem. Theory Comput. 2017, 13, 4307-4323] and our own works on new long-range corrected (LC) double hybrids for singlet-singlet and singlet-triplet excitations [ J. Chem. Theory Comput. 2019, 15, 4735-4744 and J. Chem. Phys. 2020, 153, 064106], we present new LC double hybrids with SCS/SOS that demonstrate further improvement over previously published results and methods. We introduce new unscaled and scaled versions of different global and LC double hybrids based on Becke88 or PBE exchange combined with LYP, PBE, or P86 correlation. For singlet-singlet excitations, we cross-validate them on six benchmark sets that cover small to medium-sized chromophores with different excitation types (local-valence, Rydberg, and charge transfer). For singlet-triplet excitations, we perform the cross-validation on three different benchmark sets following the same analysis as in our previous work in 2020. In total, 203 excitations are analyzed. Our results confirm and extend those of Schwabe and Goerigk regarding the superior performance of SCS and SOS variants compared to their unscaled parents by decreasing mean absolute deviations, root-mean-square deviations, or error spans by more than half and bringing absolute mean deviations closer to zero. Our SCS/SOS variants are shown to be highly efficient and robust for the computation of vertical excitation energies, which even outperform specialized double hybrids that also contain an LC in their perturbative part. In particular, our new SCS/SOS-ωPBEPP86 and SCS/SOS-ωB88PP86 functionals are four of the most accurate and robust methods tested in this work, and we fully recommend them for future applications. However, if the relevant SCS and SOS algorithms are not available to the user, we suggest ωPBEPP86 as the best unscaled method in this work.
Collapse
Affiliation(s)
- Marcos Casanova-Páez
- School of Chemistry, The University of Melbourne, Parkville, Victoria 3010, Australia
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Parkville, Victoria 3010, Australia
| |
Collapse
|
35
|
Kanchanakungwankul S, Truhlar DG. Examination of How Well Long-Range-Corrected Density Functionals Satisfy the Ionization Energy Theorem. J Chem Theory Comput 2021; 17:4823-4830. [PMID: 34319716 DOI: 10.1021/acs.jctc.1c00440] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We calculated the vertical ionization energies (VIE) of 99 species in two ways to examine the accuracy of several long-range-corrected (LC) hybrid meta functionals in comparison with a gradient approximation (GA), global hybrids, and doubly hybrids. In the category of LC functionals, we examined both those with meta ingredients (i.e., that depend on the kinetic energy density) and those without them. The LC-hybrid meta functionals examined are M11, revM11, M11plus, and ωB97M-V. The reference data used to assess accuracy consist of 95 molecules and 4 atoms in the GW100 set. The two methods studied are the ΔSCF method (involving the difference of neutral and cation self-consistent field (SCF) energies) and the ionization energy theorem (involving the orbital energy of the highest occupied molecular orbital, HOMO). We calculated linear correlation coefficients (r2) and mean absolute deviations (MADs) between each approach and the reference VIE value from the CCSD(T)/def2-TZVPP level of theory. We compared the new LC-hybrid meta calculations to calculations with the 10 functionals in a previous VIE study by Brémond et al. and to the calculations with LC-BLYP (LC-Becke, Lee-Yang-Parr), CAM-B3LYP (Coulomb-attenuating-method Becke-3-parameter Lee-Yang-Parr), LC-ωHPBE, and ωB97X-D. The results show that Minnesota LC-hybrid meta functionals have the smallest mean absolute deviation of ionization energy theorem VIEs with the reference data; the LC-ωHPBE functional also does quite well in this test. This is very encouraging and indicates that LC-hybrid meta functionals would be the best starting points for the tuning strategy that has been shown to be a very good procedure for improving time-dependent density functional calculations, and it also helps explain the good success of LC-hybrid meta functionals for molecular excitation energies.
Collapse
Affiliation(s)
- Siriluk Kanchanakungwankul
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
36
|
Mester D, Kállay M. Spin-Scaled Range-Separated Double-Hybrid Density Functional Theory for Excited States. J Chem Theory Comput 2021; 17:4211-4224. [PMID: 34152771 PMCID: PMC8280718 DOI: 10.1021/acs.jctc.1c00422] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2021] [Indexed: 11/28/2022]
Abstract
Our recently presented range-separated (RS) double-hybrid (DH) time-dependent density functional approach [J. Chem. Theory Comput. 17, 927 (2021)] is combined with spin-scaling techniques. The proposed spin-component-scaled (SCS) and scaled-opposite-spin (SOS) variants are thoroughly tested for almost 500 excitations including the most challenging types. This comprehensive study provides useful information not only about the new approaches but also about the most prominent methods in the DH class. The benchmark calculations confirm the robustness of the RS-DH ansatz, while several tendencies and deficiencies are pointed out for the existing functionals. Our results show that the SCS variant consistently improves the results, while the SOS variant preserves the benefits of the original RS-DH method reducing its computational expenses. It is also demonstrated that, besides our approaches, only the nonempirical functionals provide balanced performance for general applications, while particular methods are only suggested for certain types of excitations.
Collapse
Affiliation(s)
- Dávid Mester
- Department of Physical Chemistry and
Materials Science, Budapest University of
Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and
Materials Science, Budapest University of
Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| |
Collapse
|
37
|
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
|
38
|
Santra G, Cho M, Martin JML. Exploring Avenues beyond Revised DSD Functionals: I. Range Separation, with xDSD as a Special Case. J Phys Chem A 2021; 125:4614-4627. [PMID: 34009986 PMCID: PMC8279641 DOI: 10.1021/acs.jpca.1c01294] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2021] [Revised: 05/06/2021] [Indexed: 01/16/2023]
Abstract
We have explored the use of range separation as a possible avenue for further improvement on our revDSD minimally empirical double hybrid functionals. Such ωDSD functionals encompass the XYG3 type of double hybrid (i.e., xDSD) as a special case for ω → 0. As in our previous studies, the large and chemically diverse GMTKN55 benchmark suite was used for evaluation. Especially when using the D4 rather than D3BJ dispersion model, xDSD has a slight performance advantage in WTMAD2. As in previous studies, PBEP86 is the winning combination for the semilocal parts. xDSDn-PBEP86-D4 marginally outperforms the previous "best in class" ωB97M(2) Berkeley double hybrid but without range separation and using fewer than half the number of empirical parameters. Range separation turns out to offer only marginal further improvements on GMTKN55 itself. While ωB97M(2) still yields better performance for small-molecule thermochemistry, this is compensated in WTMAD2 by the superior performance of the new functionals for conformer equilibria. Results for two external test sets with pronounced static correlation effects may indicate that range-separated double hybrids are more resilient to such effects.
Collapse
Affiliation(s)
- Golokesh Santra
- Department
of Organic Chemistry, Weizmann Institute
of Science, 7610001 Reḥovot, Israel
| | - Minsik Cho
- Department
of Organic Chemistry, Weizmann Institute
of Science, 7610001 Reḥovot, Israel
- Department
of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| | - Jan M. L. Martin
- Department
of Organic Chemistry, Weizmann Institute
of Science, 7610001 Reḥovot, Israel
| |
Collapse
|
39
|
Najibi A, Casanova-Páez M, Goerigk L. Analysis of Recent BLYP- and PBE-Based Range-Separated Double-Hybrid Density Functional Approximations for Main-Group Thermochemistry, Kinetics, and Noncovalent Interactions. J Phys Chem A 2021; 125:4026-4035. [DOI: 10.1021/acs.jpca.1c02549] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Asim Najibi
- School of Chemistry, The University of Melbourne, Parkville 3010, Australia
| | | | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Parkville 3010, Australia
| |
Collapse
|
40
|
Brémond É, Ottochian A, Pérez-Jiménez ÁJ, Ciofini I, Scalmani G, Frisch MJ, Sancho-García JC, Adamo C. Assessing challenging intra- and inter-molecular charge-transfer excitations energies with double-hybrid density functionals. J Comput Chem 2021; 42:970-981. [PMID: 33748983 DOI: 10.1002/jcc.26517] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2021] [Accepted: 02/15/2021] [Indexed: 12/24/2022]
Abstract
We investigate the performance of a set of recently introduced range-separated double-hybrid functionals, namely ωB2-PLYP, ωB2GP-PLYP, RSX-0DH, and RSX-QIDH models for hard-to-calculate excitation energies. We compare with the parent (B2-PLYP, B2GP-PLYP, PBE0-DH, and PBE-QIDH) and other (DSD-PBEP86) double-hybrid models as well as with some of the most widely employed hybrid functionals (B3LYP, PBE0, M06-2X, and ωB97X). For this purpose, we select a number of medium-sized intra- and inter-molecular charge-transfer excitations, which are known to be challenging to calculate using time-dependent density-functional theory (TD-DFT) and for which accurate reference values are available. We assess whether the high accuracy shown by the newest double-hybrid models is also confirmed for those cases too. We find that asymptotically corrected double-hybrid models yield a superior performance, especially for the inter-molecular charge-transfer excitation energies, as compared to standard double-hybrid models. Overall, the PBE-QIDH and its corresponding range-separated RSX-QIDH functional are recommended for general-purpose TD-DFT applications, depending on whether long-range effects are expected to play a significant role.
Collapse
Affiliation(s)
- Éric Brémond
- Université de Paris, ITODYS, CNRS, Paris, F-75006, France
| | - Alistar Ottochian
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences (i-CLeHS, Paris, France
| | | | - Ilaria Ciofini
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences (i-CLeHS, Paris, France
| | | | | | | | - Carlo Adamo
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences (i-CLeHS, Paris, France.,Institut Universitaire de France, Paris, France
| |
Collapse
|
41
|
Zhang IY, Xu X. Exploring the Limits of the XYG3-Type Doubly Hybrid Approximations for the Main-Group Chemistry: The xDH@B3LYP Model. J Phys Chem Lett 2021; 12:2638-2644. [PMID: 33689361 DOI: 10.1021/acs.jpclett.1c00360] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
Despite being the most accurate class of density functional approximations for the main-group chemistry, doubly hybrid approximations (DHAs) are generally considered to be incomplete in describing the medium- to long-range dispersive interactions. The existing DHAs are often supplemented with empirical long-range dispersion corrections. By using the extensive and chemically diverse GMTKN55 database, we explore the limits of the XYG3-type DHAs using the B3LYP reference orbitals, namely, xDH@B3LYP, with a gradually relaxed constraint on the mixing parameters of DHAs. Our results demonstrate that the xDH@B3LYP model can provide a balanced description of both covalent and noncovalent interactions with the accuracy and robustness comparable to or even better than the very expensive composite methods in wave function theory. Such an accuracy can be achieved without resorting to the use of any long-range correction scheme, shedding new light on the development of DHAs.
Collapse
Affiliation(s)
- Igor Ying Zhang
- Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, MOE Key Laboratory of Computational Physical Sciences, Shanghai Key Laboratory of Bioactive Small Molecules, Department of Chemistry, Fudan University, Shanghai 200433, People's Republic of China
| | - Xin Xu
- Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, MOE Key Laboratory of Computational Physical Sciences, Shanghai Key Laboratory of Bioactive Small Molecules, Department of Chemistry, Fudan University, Shanghai 200433, People's Republic of China
| |
Collapse
|
42
|
Mester D, Kállay M. A Simple Range-Separated Double-Hybrid Density Functional Theory for Excited States. J Chem Theory Comput 2021; 17:927-942. [PMID: 33400872 PMCID: PMC7884002 DOI: 10.1021/acs.jctc.0c01135] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2020] [Indexed: 01/12/2023]
Abstract
A simple and robust range-separated (RS) double-hybrid (DH) time-dependent density functional approach is presented for the accurate calculation of excitation energies of molecules within the Tamm-Dancoff approximation. The scheme can be considered as an excited-state extension of the ansatz proposed by Toulouse and co-workers [J. Chem. Phys. 2018, 148, 164105], which is based on the two-parameter decomposition of the Coulomb potential, for which both the exchange and correlation contributions are range-separated. A flexible and efficient implementation of the new scheme is also presented, which facilitates its extension to any combination of exchange and correlation functionals. The performance of the new approximation is tested for singlet excitations on several benchmark compilations and thoroughly compared to that of representative DH, RS hybrid, and RS DH functionals. The one-electron basis set dependence and computation times are also assessed. Our results show that the new approach improves on standard DHs in most cases, and it can provide a more robust and accurate alternative. In addition, on average, it noticeably surpasses the existing RS hybrid and RS DH functionals.
Collapse
Affiliation(s)
- Dávid Mester
- Department of Physical Chemistry and
Materials Science, Budapest University of
Technology and Economics, P.O. Box 91, Budapest, H-1521, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and
Materials Science, Budapest University of
Technology and Economics, P.O. Box 91, Budapest, H-1521, Hungary
| |
Collapse
|
43
|
Goerigk L, Casanova-Paéz M. The Trip to the Density Functional Theory Zoo Continues: Making a Case for Time-Dependent Double Hybrids for Excited-State Problems. Aust J Chem 2021. [DOI: 10.1071/ch20093] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
This account is written for general users of time-dependent density functional theory (TD-DFT) methods as well as chemists who are unfamiliar with the field. It includes a brief overview of conventional TD-DFT approaches and recommendations for applications to organic molecules based on our own experience. The main emphasis of this work, however, lies in providing the first in-depth review of time-dependent double-hybrid density functionals. They were first established in 2007 with very promising follow-up studies in the subsequent four years before developments or applications became scarce. The topic has regained more interest since 2017, and this account reviews those latest developments led by our group. These developments have shown unprecedented robustness for a variety of different types of electronic excitations when compared to more conventional TD-DFT methods. In particular, time-dependent double hybrids do not suffer from artificial ghost states and are able to reproduce exciton-coupled absorption spectra. Our latest methods include range separation and belong to the currently best TD-DFT methods for singlet-singlet excitations in organic molecules. While there is still room for improvement and further development in this space, we hope that this account encourages users to adjust their computational protocols to such new methods to provide more real-life testing and scenarios.
Collapse
|
44
|
Graf D, Ochsenfeld C. A range-separated generalized Kohn-Sham method including a long-range nonlocal random phase approximation correlation potential. J Chem Phys 2020; 153:244118. [PMID: 33380112 DOI: 10.1063/5.0031310] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
Based on our recently published range-separated random phase approximation (RPA) functional [Kreppel et al., "Range-separated density-functional theory in combination with the random phase approximation: An accuracy benchmark," J. Chem. Theory Comput. 16, 2985-2994 (2020)], we introduce self-consistent minimization with respect to the one-particle density matrix. In contrast to the range-separated RPA methods presented so far, the new method includes a long-range nonlocal RPA correlation potential in the orbital optimization process, making it a full-featured variational generalized Kohn-Sham (GKS) method. The new method not only improves upon all other tested RPA schemes including the standard post-GKS range-separated RPA for the investigated test cases covering general main group thermochemistry, kinetics, and noncovalent interactions but also significantly outperforms the popular G0W0 method in estimating the ionization potentials and fundamental gaps considered in this work using the eigenvalue spectra obtained from the GKS Hamiltonian.
Collapse
Affiliation(s)
- Daniel Graf
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), D-81377 Munich, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), D-81377 Munich, Germany
| |
Collapse
|
45
|
Gerrits N, Smeets EWF, Vuckovic S, Powell AD, Doblhoff-Dier K, Kroes GJ. Density Functional Theory for Molecule-Metal Surface Reactions: When Does the Generalized Gradient Approximation Get It Right, and What to Do If It Does Not. J Phys Chem Lett 2020; 11:10552-10560. [PMID: 33295770 PMCID: PMC7751010 DOI: 10.1021/acs.jpclett.0c02452] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/10/2023]
Abstract
While density functional theory (DFT) is perhaps the most used electronic structure theory in chemistry, many of its practical aspects remain poorly understood. For instance, DFT at the generalized gradient approximation (GGA) tends to fail miserably at describing gas-phase reaction barriers, while it performs surprisingly well for many molecule-metal surface reactions. GGA-DFT also fails for many systems in the latter category, and up to now it has not been clear when one may expect it to work. We show that GGA-DFT tends to work if the difference between the work function of the metal and the molecule's electron affinity is greater than ∼7 eV and to fail if this difference is smaller, with sticking of O2 on Al(111) being a spectacular example. Using dynamics calculations we show that, for this system, the DFT problem may be solved as done for gas-phase reactions, i.e., by resorting to hybrid functionals, but using screening at long-range to obtain a correct description of the metal. Our results suggest the GGA error in the O2 + Al(111) barrier height to be functional driven. Our results also suggest the possibility to compute potential energy surfaces for the difficult-to-treat systems with computationally cheap nonself-consistent calculations in which a hybrid functional is applied to a GGA density.
Collapse
Affiliation(s)
- Nick Gerrits
- Leiden
Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Egidius W. F. Smeets
- Leiden
Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Stefan Vuckovic
- Department
of Chemistry, University of California, Irvine, California 92697, United States
| | - Andrew D. Powell
- Leiden
Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Katharina Doblhoff-Dier
- Leiden
Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| | - Geert-Jan Kroes
- Leiden
Institute of Chemistry, Leiden University, Gorlaeus Laboratories, P.O. Box 9502, 2300 RA Leiden, The Netherlands
| |
Collapse
|
46
|
Jana S, Śmiga S, Constantin LA, Samal P. Generalizing Double-Hybrid Density Functionals: Impact of Higher-Order Perturbation Terms. J Chem Theory Comput 2020; 16:7413-7430. [PMID: 33205659 PMCID: PMC7735712 DOI: 10.1021/acs.jctc.0c00823] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2020] [Indexed: 11/28/2022]
Abstract
Connections between the Görling-Levy (GL) perturbation theory and the parameters of double-hybrid (DH) density functional are established via adiabatic connection formalism. Moreover, we present a more general DH density functional theory, where the higher-order perturbation terms beyond the second-order GL2 one, such as GL3 and GL4, also contribute. It is shown that a class of DH functionals including previously proposed ones can be formed using the present construction. Based on the proposed formalism, we assess the performance of higher-order DH and long-range corrected DH formed on the Perdew-Burke-Ernzerhof (PBE) semilocal functional and second-order GL2 correlation energy. The underlying construction of DH functionals based on the generalized many-body perturbation approaches is physically appealing in terms of the development of the non-local forms using more accurate and sophisticated semilocal functionals.
Collapse
Affiliation(s)
- Subrata Jana
- School
of Physical Sciences, National Institute
of Science Education and Research, HBNI, Bhubaneswar 752050, India
| | - Szymon Śmiga
- Institute
of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziadzka 5, 87-100 Toruń, Poland
| | - Lucian A. Constantin
- Consiglio
Nazionale delle Ricerche CNR-NANO, Istituto
di Nanoscienze, 41125 Modena, Italy
| | - Prasanjit Samal
- School
of Physical Sciences, National Institute
of Science Education and Research, HBNI, Bhubaneswar 752050, India
| |
Collapse
|
47
|
Ediger MD, Jensen L, Manolopoulos DE, Martinez TJ, Michaelides A, Reichman DR, Sherrill CD, Shi Q, Straub JE, Vega C, Wang LS, Brigham EC, Lian T. JCP Emerging Investigator Special Collection 2019. J Chem Phys 2020; 153:110402. [PMID: 32962387 DOI: 10.1063/5.0021946] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Affiliation(s)
- Mark D Ediger
- Department of Chemistry, University of Wisconsin-Madison, Madison, Wisconsin 53706, USA
| | - Lasse Jensen
- Department of Chemistry, Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | - David E Manolopoulos
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| | - Todd J Martinez
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Angelos Michaelides
- Thomas Young Centre and London Centre for Nanotechnology, 17-19 Gordon Street, London WC1H 0AH, United Kingdom and Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - David R Reichman
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Qiang Shi
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences, Zhongguancun, Beijing 100190, China; University of Chinese Academy of Sciences, Beijing 100049, China; and Physical Science Laboratory, Huairou National Comprehensive Science Center, Beijing 101407, China
| | - John E Straub
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, USA
| | - Carlos Vega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Lai-Sheng Wang
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| | | | - Tianquan Lian
- Department of Chemistry, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
48
|
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
|
49
|
Brémond É, Pérez-Jiménez ÁJ, Sancho-García JC, Adamo C. Range-separated hybrid and double-hybrid density functionals: A quest for the determination of the range-separation parameter. J Chem Phys 2020; 152:244124. [PMID: 32610956 DOI: 10.1063/5.0010976] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We recently derived a new and simple route to the determination of the range-separation parameter in range-separated exchange hybrid and double-hybrid density functionals by imposing an additional constraint to the exchange-correlation energy to recover the total energy of the hydrogen atom [Brémond et al., J. Chem. Phys. 15, 201102 (2019)]. Here, we thoroughly assess this choice by statistically comparing the derived values of the range-separation parameters to the ones obtained using the optimal tuning (OT) approach. We show that both approaches closely agree, thus, confirming the reliability of ours. We demonstrate that it provides very close performances in the computation of properties particularly prone to the one- and many-electron self-interaction errors (i.e., ionization potentials). Our approach arises as an alternative to the OT procedure, conserving the accuracy and efficiency of a standard Kohn-Sham approach to density-functional theory computation.
Collapse
Affiliation(s)
- Éric Brémond
- Université de Paris, ITODYS, CNRS, F-75006 Paris, France
| | | | | | - Carlo Adamo
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences (i-CLeHS), FRE 2027, F-75005 Paris, France
| |
Collapse
|
50
|
Brémond É, Pérez-Jiménez ÁJ, Adamo C, Sancho-García JC. sp-hybridized carbon allotrope molecular structures: An ongoing challenge for density-functional approximations. J Chem Phys 2019; 151:211104. [PMID: 31822105 DOI: 10.1063/1.5133639] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The recent synthesis of a C18 monocyclic ring constitutes a major breakthrough as a new all-carbon disclosed form. However, modern density functional theory approaches do not lead to the correct experimental polyynic structure and favor the cumulenic one instead. We demonstrate here that this serious drawback can be solved by recently developed range-separated nonempirical schemes, independently of which kind of functional is being applied (i.e., semilocal, hybrid, or double-hybrid).
Collapse
Affiliation(s)
- Éric Brémond
- Université de Paris, ITODYS, UMR CNRS 7086, 15 Rue J.-A. de Baïf, F-75013 Paris, France
| | | | - Carlo Adamo
- Chimie ParisTech, PSL Research University, CNRS, Institute of Chemistry for Life and Health Sciences (i-CLeHS), FRE 2027, F-75005 Paris, France
| | | |
Collapse
|