1
|
Vrška D, Urban M, Neogrády P, Pittner J, Blaško M, Pitoňák M. DFT Modeling of Polyethylene Chains Cross-linked by Selected ns 1 and ns 2 Metal Atoms. J Phys Chem A 2024; 128:7634-7647. [PMID: 39219502 DOI: 10.1021/acs.jpca.4c04755] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/04/2024]
Abstract
We analyze the structures, stabilities, and thermochemical properties of polyethylene (PE) oligomer chains cross-linked by metal (M) atoms through C-M-C bonds. Representative PEn-Mm-PEn complexes contain between 7 and 15 carbon atoms in each oligomer and one to three Li, Be, Mg, Zn, Ag, or Au cross-linking metal elements. PEn-Mm-PEn complexes are quasiplanar with nearly parallel PE chains. Their stability is determined by covalent C-M-C bonds accompanied by noncovalent dispersion interactions between PEn chains. Using the CAM-B3LYP+D3BJ+ABC functional, the binding energies of PE15-M-PE15 with respect to two PE15 radicals and metal fragments are -225, -230, -322, -551, -289, and -303 kJ/mol for Li, Ag, Au, Be, Mg, and Zn atoms, respectively. Entropy contributions (109 to 121 kJ/mol at 298.15 K) destabilize all complexes significantly. With two cross-linking metal elements in PE15-M2-PE15 complexes, binding energies are about double. Complexes with several open-shell Li, Ag, or Au doublet atoms have spins located on separated C-M-C bonds. High-spin PE15-Mm1-PE15-Mm2-PE15 complexes of three PE oligomers cross-linked by up to five doublet metal atoms create parallel PE tubes, which are suggested as elementary cells for modeling magnetic polymer tubes.
Collapse
Affiliation(s)
- Dávid Vrška
- Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University in Bratislava, Ilkovičova 6, Bratislava SK-84215, Slovakia
| | - Miroslav Urban
- Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University in Bratislava, Ilkovičova 6, Bratislava SK-84215, Slovakia
| | - Pavel Neogrády
- Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University in Bratislava, Ilkovičova 6, Bratislava SK-84215, Slovakia
| | - Jiří Pittner
- Jaroslav Heyrovský Institute of Physical Chemistry of the Czech Academy of Sciences, Dolejškova 3, Prague 8 CZ-18223, Czech Republic
| | - Martin Blaško
- Centre of Operations Computing Centre of SAS, Slovak Academy of Sciences, Dúbravská cesta 9, Bratislava SK-84535, Slovakia
| | - Michal Pitoňák
- Department of Physical and Theoretical Chemistry, Faculty of Natural Sciences, Comenius University in Bratislava, Ilkovičova 6, Bratislava SK-84215, Slovakia
- Slovak National Supercomputing Centre, Dúbravská cesta 9, Bratislava SK-84104, Slovakia
| |
Collapse
|
2
|
Müller M, Hansen A, Grimme S. ωB97X-3c: A composite range-separated hybrid DFT method with a molecule-optimized polarized valence double-ζ basis set. J Chem Phys 2023; 158:014103. [PMID: 36610980 DOI: 10.1063/5.0133026] [Citation(s) in RCA: 20] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
A new composite density functional theory (DFT) method is presented. It is based on ωB97X-V as one of the best-performing density functionals for the GMTKN55 thermochemistry database and completes the family of "3c" methods toward range-separated hybrid DFT. This method is consistently available for all elements up to Rn (Z = 1-86). Its further key ingredients are a polarized valence double-ζ (vDZP) Gaussian basis set, which was fully optimized in molecular DFT calculations, in combination with large-core effective core potentials and a specially adapted D4 dispersion correction. Unlike most existing double-ζ atomic orbital sets, vDZP shows only small basis set superposition errors (BSSEs) and can compete with standard sets of triple-ζ quality. Small residual BSSE effects are efficiently absorbed by the D4 damping scheme, which overall eliminates the need for an explicit treatment or empirical corrections for BSSE. Thorough tests on a variety of thermochemistry benchmark sets show that the new composite method, dubbed ωB97X-3c, is on par with or even outperforms standard hybrid DFT methods in a quadruple-zeta basis set at a small fraction of the computational cost. Particular strengths of this method are the description of non-covalent interactions and barrier heights, for which it is among the best-performing density functionals overall.
Collapse
Affiliation(s)
- Marcel Müller
- Mulliken Center for Theoretical Chemistry, Clausius-Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Clausius-Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Clausius-Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| |
Collapse
|
3
|
Villot C, Ballesteros F, Wang D, Lao KU. Coupled Cluster Benchmarking of Large Noncovalent Complexes in L7 and S12L as Well as the C 60 Dimer, DNA-Ellipticine, and HIV-Indinavir. J Phys Chem A 2022; 126:4326-4341. [PMID: 35766331 DOI: 10.1021/acs.jpca.2c01421] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In this work, we report the benchmark binding energies of the seven complexes within the L7 data set, six host-guest complexes from the S12L data set, a C60 dimer, the DNA-ellipticine intercalation complex, and the largest system of the study, the HIV-indinavir system, which contained 343 atoms or 139 heavy atoms. The high-quality values reported were obtained via a focal point method that relies on the canonical form of second-order Møller-Plesset theory and the domain-based local pair natural orbital scheme for the coupled cluster with single double and perturbative triple excitations [DLPNO-CCSD(T)] extrapolated to the complete basis set (CBS) limit. The results in this work not only corroborate but also improve upon some previous benchmark values for large noncovalent complexes albeit at a relatively steep cost. Although local CCSD(T) and the largely successful fixed-node diffusion Monte Carlo (FN-DMC) have been shown to generally agree for small- to medium-size systems, a discrepancy in their reported binding energy values arises for large complexes, where the magnitude of the disagreement is a definite cause for concern. For example, the largest deviation in the L7 data set was 2.8 kcal/mol (∼10%) on the low end in C3GC. Such a deviation only grows worse in the S12L set, which showed a difference of up to 10.4 kcal/mol (∼25%) by a conservative estimation in buckycatcher-C60. The DNA-ellipticine complex also generated a disagreement of 4.4 kcal/mol (∼10%) between both state-of-the-art methods. The disagreement between local CCSD(T) and FN-DMC in large noncovalent complexes shows that it is urgently needed to have the canonical CCSD(T), the Monte Carlo CCSD(T), or the full configuration interaction quantum Monte Carlo approaches available to large systems on the hundred-atom scale to solve this dilemma. In addition, the performances of cheaper popular computational methods were assessed for the studied complexes with respect to DLPNO-CCSD(T)/CBS. r2SCAN-3c, B97M-V, and PBE0+D4 work well in large noncovalent complexes in this work, and GFN2-xTB performs well in π-π stacking complexes. B97M-V is the most reliable computationally efficient approach to predicting noncovalent interactions for large complexes, being the only one to have binding errors within the so-called 1 kcal/mol "chemical accuracy". The benchmark interaction energies of these host-guest complexes, molecular materials, and biological systems with electronic and medicinal implications provide crucial reference data for the improvement of current and future lower-cost methods.
Collapse
Affiliation(s)
- Corentin Villot
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| | - Francisco Ballesteros
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| | - Danyang Wang
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| | - Ka Un Lao
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| |
Collapse
|
4
|
Gasevic T, Stückrath JB, Grimme S, Bursch M. Optimization of the r 2SCAN-3c Composite Electronic-Structure Method for Use with Slater-Type Orbital Basis Sets. J Phys Chem A 2022; 126:3826-3838. [PMID: 35654439 PMCID: PMC9255700 DOI: 10.1021/acs.jpca.2c02951] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The "Swiss army knife" composite density functional electronic-structure method r2SCAN-3c (J. Chem. Phys. 2021, 154, 064103) is extended and optimized for the use with Slater-type orbital basis sets. The meta generalized-gradient approximation (meta-GGA) functional r2SCAN by Furness et al. is combined with a tailor-made polarized triple-ζ Slater-type atomic orbital (STO) basis set (mTZ2P), the semiclassical London dispersion correction (D4), and a geometrical counterpoise (gCP) correction. Relativistic effects are treated explicitly with the scalar-relativistic zeroth-order regular approximation (SR-ZORA). The performance of the new implementation is assessed on eight geometry and 74 energy benchmark sets, including the extensive GMTKN55 database as well as recent sets such as ROST61 and IONPI19. In geometry optimizations, the STO-based r2SCAN-3c is either on par with or more accurate than the hybrid density functional approximation M06-2X-D3(0)/TZP. In energy calculations, the overall accuracy is similar to the original implementation of r2SCAN-3c with Gaussian-type atomic orbitals (GTO), but basic properties, intermolecular noncovalent interactions, and barrier heights are better described with the STO approach, resulting in a lower weighted mean absolute deviation (WTMAD-2(STO) = 7.15 vs 7.50 kcal mol-1 with the original method) for the GMTKN55 database. The STO-optimized r2SCAN-3c outperforms many conventional hybrid/QZ approaches in most common applications at a fraction of their cost. The reliable, robust, and accurate r2SCAN-3c implementation with STOs is a promising alternative to the original implementation with GTOs and can be generally used for a broad field of quantum chemical problems.
Collapse
Affiliation(s)
- Thomas Gasevic
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, D-53115 Bonn, Germany
| | - Julius B Stückrath
- 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
| | - Markus Bursch
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
5
|
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
|