1
|
Zhu A, Tew DP. Wannier Function Localization Using Bloch Intrinsic Atomic Orbitals. J Phys Chem A 2024; 128:8570-8579. [PMID: 39295205 PMCID: PMC11457216 DOI: 10.1021/acs.jpca.4c04555] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2024] [Revised: 09/09/2024] [Accepted: 09/11/2024] [Indexed: 09/21/2024]
Abstract
We extend the intrinsic atomic orbital (IAO) method for the localization of molecular orbitals to calculate well-localized generalized Wannier functions in crystals in the spirit of the Pipek-Mezey method. We furthermore present a one-shot diabatic Wannierization procedure that aligns the phases of the Bloch functions, providing immediate Wannier localization, which serves as an excellent initial guess for optimization. We test our Wannier localization implementation on a number of solid-state systems, highlighting the effectiveness of the diabatic preparation, especially for localizing core bands. Partial charges of Wannier functions generated using Bloch IAOs align well with chemical intuition, which we demonstrate through the example of the adsorption of CO on a MgO surface.
Collapse
Affiliation(s)
- Andrew Zhu
- Physical & Theoretical
Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, U.K.
| | - David P. Tew
- Physical & Theoretical
Chemistry Laboratory, University of Oxford, Oxford OX1 3QZ, U.K.
| |
Collapse
|
2
|
Roth E, Listyarini RV, Hofer TS, Cziferszky M. Host-Guest Interactions of Ruthenium(II) Arene Complexes with Cucurbit[7/8]uril. Inorg Chem 2024; 63:14021-14031. [PMID: 39016439 PMCID: PMC11289748 DOI: 10.1021/acs.inorgchem.4c01755] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2024] [Revised: 07/10/2024] [Accepted: 07/11/2024] [Indexed: 07/18/2024]
Abstract
Cucurbit[n]urils (CB[n]s) have been recognized for their chemical and thermal stability, and their ability to bind many neutral and cationic guest molecules makes them excellent hosts in a range of supramolecular applications. In drug delivery, CB[n]s can enhance drug solubility, improve chemical and physical drug stability, and allow for triggered and controlled release. This study aimed to investigate the ability of CB[7] and CB[8] as molecular hosts to bind ruthenium(II) arene complexes that are current anticancer lead structures in the area of metallodrugs. Both, experimental and computational methods, led to insights into the binding preferences and geometries of [RuII(cym)Cl2]2 (1; cym = η6-p-cymene), [RuII(cym)(dmb)Cl2]) (2; cym = η6-p-cymene; dmb = 1,3-dimethylbenzimidazol-2-ylidene), and [RuII(cym)(pta)Cl2] (3, RAPTA-C; cym = η6-p-cymene; pta = 1,3,5-triaza-7-phospha-adamantane) with CB[7] and CB[8]. Competition experiments by mass spectrometry revealed clear preferences of 2 for CB[8] and 3 for CB[7]. Based on a comparison of the associated interaction energies from quantum chemical calculations as well as experimental data, 3@CB[7] clearly prefers a binding mode, where the pta ligand is located inside the cavity of the host, and the metal ion interacts with two of the carbonyl groups on the rim of CB[7]. In contrast, complex 2 binds in two different orientations with interaction energies similar to those of both CB[n]s, with the cym ligand being either inside or outside of the cavity. These findings suggest that ruthenium(II) arene complexes are able to form stable host-guest interactions with CB[n]s, which can be exploited as drug delivery vehicles in further metallodrug development to improve their chemical stability.
Collapse
Affiliation(s)
- Elisa Roth
- Institute
for Pharmacy, Pharmaceutical Chemistry, Department of Chemistry and
Pharmacy, University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| | - Risnita Vicky Listyarini
- Institute
of General, Inorganic and Theoretical Chemistry, Center for Chemistry
and Biomedicine, University of Innsbruck, Innrain 80-82, A-6020 Innsbruck, Austria
- Chemistry
Education Study Program, Sanata Dharma University, Yogyakarta 55282, Indonesia
| | - Thomas S. Hofer
- Institute
of General, Inorganic and Theoretical Chemistry, Center for Chemistry
and Biomedicine, University of Innsbruck, Innrain 80-82, A-6020 Innsbruck, Austria
| | - Monika Cziferszky
- Institute
for Pharmacy, Pharmaceutical Chemistry, Department of Chemistry and
Pharmacy, University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| |
Collapse
|
3
|
Wodyński A, Lauw B, Reimann M, Kaupp M. Spin-Symmetry Breaking and Hyperfine Couplings in Transition-Metal Complexes Revisited Using Density Functionals Based on the Exact-Exchange Energy Density. J Chem Theory Comput 2024; 20:2033-2048. [PMID: 38411554 PMCID: PMC10938646 DOI: 10.1021/acs.jctc.3c01422] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2023] [Revised: 02/02/2024] [Accepted: 02/07/2024] [Indexed: 02/28/2024]
Abstract
A small set of mononuclear manganese complexes evaluated previously for their Mn hyperfine couplings (HFCs) has been analyzed using density functionals based on the exact-exchange energy density─in particular, the spin symmetry breaking (SSB) found previously when using hybrid functionals. Employing various strong-correlation corrected local hybrids (scLHs) and strong-correlation corrected range-separated local hybrids (scRSLHs) with or without additional corrections to their local mixing functions (LMFs) to mitigate delocalization errors (DE), the SSB and the associated dipolar HFCs of [Mn(CN)4]2-, MnO3, [Mn(CN)4N]-, and [Mn(CN)5NO]2- (the latter with cluster embedding) have been examined. Both strong-correlation (sc)-correction and DE-correction terms help to diminish SSB and correct the dipolar HFCs. The DE corrections are more effective, and the effects of the sc corrections depend on their damping factors. Interestingly, the DE-corrections reduce valence-shell spin polarization (VSSP) and thus SSB by locally enhancing exact-exchange (EXX) admixture near the metal center and thereby diminishing spin-density delocalization onto the ligand atoms. In contrast, sc corrections diminish EXX admixture locally, mostly on specific ligand atoms. This then reduces VSSP and SSB as well. The performance of scLHs and scRSLHs for the isotropic Mn HFCs has also been analyzed, with particular attention to core-shell spin-polarization contributions. Further sc-corrected functionals, such as the KP16/B13 construction and the DM21 deep-neural-network functional, have been examined.
Collapse
Affiliation(s)
- Artur Wodyński
- Technische Universität
Berlin, Institut für Chemie, Theoretische
Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, Berlin, D-10623, Germany
| | - Bryan Lauw
- Technische Universität
Berlin, Institut für Chemie, Theoretische
Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, Berlin, D-10623, Germany
| | - Marc Reimann
- Technische Universität
Berlin, Institut für Chemie, Theoretische
Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, Berlin, D-10623, Germany
| | - Martin Kaupp
- Technische Universität
Berlin, Institut für Chemie, Theoretische
Chemie/Quantenchemie, Sekr. C7, Straße des 17. Juni 135, Berlin, D-10623, Germany
| |
Collapse
|
4
|
Franzke Y, Holzer C, Andersen JH, Begušić T, Bruder F, Coriani S, Della Sala F, Fabiano E, Fedotov DA, Fürst S, Gillhuber S, Grotjahn R, Kaupp M, Kehry M, Krstić M, Mack F, Majumdar S, Nguyen BD, Parker SM, Pauly F, Pausch A, Perlt E, Phun GS, Rajabi A, Rappoport D, Samal B, Schrader T, Sharma M, Tapavicza E, Treß RS, Voora V, Wodyński A, Yu JM, Zerulla B, Furche F, Hättig C, Sierka M, Tew DP, Weigend F. TURBOMOLE: Today and Tomorrow. J Chem Theory Comput 2023; 19:6859-6890. [PMID: 37382508 PMCID: PMC10601488 DOI: 10.1021/acs.jctc.3c00347] [Citation(s) in RCA: 50] [Impact Index Per Article: 50.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2023] [Indexed: 06/30/2023]
Abstract
TURBOMOLE is a highly optimized software suite for large-scale quantum-chemical and materials science simulations of molecules, clusters, extended systems, and periodic solids. TURBOMOLE uses Gaussian basis sets and has been designed with robust and fast quantum-chemical applications in mind, ranging from homogeneous and heterogeneous catalysis to inorganic and organic chemistry and various types of spectroscopy, light-matter interactions, and biochemistry. This Perspective briefly surveys TURBOMOLE's functionality and highlights recent developments that have taken place between 2020 and 2023, comprising new electronic structure methods for molecules and solids, previously unavailable molecular properties, embedding, and molecular dynamics approaches. Select features under development are reviewed to illustrate the continuous growth of the program suite, including nuclear electronic orbital methods, Hartree-Fock-based adiabatic connection models, simplified time-dependent density functional theory, relativistic effects and magnetic properties, and multiscale modeling of optical properties.
Collapse
Affiliation(s)
- Yannick
J. Franzke
- Fachbereich
Chemie, Philipps-Universität Marburg, Hans-Meerwein-Str. 4, 35032 Marburg, Germany
| | - Christof Holzer
- Institute
of Theoretical Solid State Physics, Karlsruhe
Institute of Technology (KIT), Wolfgang-Gaede-Str. 1, 76131 Karlsruhe, Germany
| | - Josefine H. Andersen
- DTU
Chemistry, Department of Chemistry, Technical
University of Denmark, Kemitorvet Building 207, DK-2800 Kongens Lyngby, Denmark
| | - Tomislav Begušić
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Florian Bruder
- Fachbereich
Chemie, Philipps-Universität Marburg, Hans-Meerwein-Str. 4, 35032 Marburg, Germany
| | - Sonia Coriani
- DTU
Chemistry, Department of Chemistry, Technical
University of Denmark, Kemitorvet Building 207, DK-2800 Kongens Lyngby, Denmark
| | - Fabio Della Sala
- Institute
for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, 73100 Lecce, Italy
- Center for
Biomolecular Nanotechnologies @UNILE, Istituto
Italiano di Tecnologia, Via Barsanti, 73010 Arnesano, Italy
| | - Eduardo Fabiano
- Institute
for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, 73100 Lecce, Italy
- Center for
Biomolecular Nanotechnologies @UNILE, Istituto
Italiano di Tecnologia, Via Barsanti, 73010 Arnesano, Italy
| | - Daniil A. Fedotov
- DTU
Chemistry, Department of Chemistry, Technical
University of Denmark, Kemitorvet Building 207, DK-2800 Kongens Lyngby, Denmark
- Institute
of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Susanne Fürst
- Institut
für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Technische Universität Berlin, Straße des 17 Juni 135, 10623, Berlin, Germany
| | - Sebastian Gillhuber
- Institute
of Inorganic Chemistry, Karlsruhe Institute
of Technology (KIT), Engesserstr. 15, 76131 Karlsruhe, Germany
| | - Robin Grotjahn
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Martin Kaupp
- Institut
für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Technische Universität Berlin, Straße des 17 Juni 135, 10623, Berlin, Germany
| | - Max Kehry
- Institute
of Physical Chemistry, Karlsruhe Institute
of Technology (KIT), Fritz-Haber-Weg 2, 76131 Karlsruhe, Germany
| | - Marjan Krstić
- Institute
of Theoretical Solid State Physics, Karlsruhe
Institute of Technology (KIT), Wolfgang-Gaede-Str. 1, 76131 Karlsruhe, Germany
| | - Fabian Mack
- Institute
of Physical Chemistry, Karlsruhe Institute
of Technology (KIT), Fritz-Haber-Weg 2, 76131 Karlsruhe, Germany
| | - Sourav Majumdar
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Brian D. Nguyen
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Shane M. Parker
- Department
of Chemistry, Case Western Reserve University, 10900 Euclid Ave, Cleveland, Ohio 44106 United States
| | - Fabian Pauly
- Institute
of Physics, University of Augsburg, Universitätsstr. 1, 86159 Augsburg, Germany
| | - Ansgar Pausch
- Institute
of Physical Chemistry, Karlsruhe Institute
of Technology (KIT), Fritz-Haber-Weg 2, 76131 Karlsruhe, Germany
| | - Eva Perlt
- Otto-Schott-Institut
für Materialforschung, Friedrich-Schiller-Universität
Jena, Löbdergraben
32, 07743 Jena, Germany
| | - Gabriel S. Phun
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Ahmadreza Rajabi
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Dmitrij Rappoport
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Bibek Samal
- Department
of Chemical Sciences, Tata Institute of
Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| | - Tim Schrader
- Otto-Schott-Institut
für Materialforschung, Friedrich-Schiller-Universität
Jena, Löbdergraben
32, 07743 Jena, Germany
| | - Manas Sharma
- Otto-Schott-Institut
für Materialforschung, Friedrich-Schiller-Universität
Jena, Löbdergraben
32, 07743 Jena, Germany
| | - Enrico Tapavicza
- Department
of Chemistry and Biochemistry, California
State University, Long Beach, 1250 Bellflower Boulevard, Long
Beach, California 90840-9507, United States
| | - Robert S. Treß
- Lehrstuhl
für Theoretische Chemie, Ruhr-Universität
Bochum, 44801 Bochum, Germany
| | - Vamsee Voora
- Department
of Chemical Sciences, Tata Institute of
Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| | - Artur Wodyński
- Institut
für Chemie, Theoretische Chemie/Quantenchemie, Sekr. C7, Technische Universität Berlin, Straße des 17 Juni 135, 10623, Berlin, Germany
| | - Jason M. Yu
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Benedikt Zerulla
- Institute
of Nanotechnology, Karlsruhe Institute of
Technology (KIT), Hermann-von-Helmholtz-Platz
1, 76344 Eggenstein-Leopoldshafen Germany
| | - Filipp Furche
- Department
of Chemistry, University of California,
Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, United States
| | - Christof Hättig
- Lehrstuhl
für Theoretische Chemie, Ruhr-Universität
Bochum, 44801 Bochum, Germany
| | - Marek Sierka
- Otto-Schott-Institut
für Materialforschung, Friedrich-Schiller-Universität
Jena, Löbdergraben
32, 07743 Jena, Germany
| | - David P. Tew
- Physical
and Theoretical Chemistry Laboratory, University
of Oxford, South Parks
Road, Oxford OX1 3QZ, United Kingdom
| | - Florian Weigend
- Fachbereich
Chemie, Philipps-Universität Marburg, Hans-Meerwein-Str. 4, 35032 Marburg, Germany
| |
Collapse
|
5
|
Schrader T, Perlt E, Fritz T, Sierka M. Performance of Common Density Functionals for Excited States of Tetraphenyldibenzoperiflanthene. J Phys Chem A 2023; 127:3265-3273. [PMID: 37037005 DOI: 10.1021/acs.jpca.2c06715] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/12/2023]
Abstract
Time-dependent density functional theory is the method of choice to efficiently calculate excitation spectra with the functional and basis set choice allowing one to compromise between accuracy and computational cost. In this work, the performance of different functionals as well as the second-order approximate coupled cluster singles and doubles model CC2 is evaluated by comparing the results to experimental results of the example molecule tetraphenyldibenzoperiflanthene (DBP). The choice of the functional has a significant impact on the calculated spectrum of DBP. The performance of a number of different functionals was evaluated, quantified, and, where possible, discussed. The best functional, tuned-CAM-B3LYP, is used to investigate DBP on a surface of hexagonal boron nitride (h-BN). The resulting spectrum shows excellent agreement with experimental results for a monolayer of DBP on h-BN.
Collapse
Affiliation(s)
- Tim Schrader
- Otto Schott Institute of Materials Research, Friedrich Schiller University Jena, 07743 Jena, Germany
| | - Eva Perlt
- Otto Schott Institute of Materials Research, Friedrich Schiller University Jena, 07743 Jena, Germany
| | - Torsten Fritz
- Institute of Solid State Physics, Friedrich Schiller University Jena, 07743 Jena, Germany
| | - Marek Sierka
- Otto Schott Institute of Materials Research, Friedrich Schiller University Jena, 07743 Jena, Germany
| |
Collapse
|
6
|
Díaz Mirón JEZ, Stein M. A benchmark for non-covalent interactions in organometallic crystals. Phys Chem Chem Phys 2022; 24:29338-29349. [PMID: 36448535 DOI: 10.1039/d2cp04160j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
Organometallic complexes are the basis for homogeneous catalysis, have applications in materials science and are also active pharmaceutical ingredients. The interaction between transition metal complexes in the solid state is determining their thermodynamics and bio-availability. Non-covalent interactions such as hydrogen bonding and van der Waals are stabilizing crystals of transition metal complexes. The variation of ligand field, central metal atoms and their oxidation and spin states are determinants of the magnitude of their inter-molecular interactions. A comparison of a set of 43 manually curated experimental heats of sublimation (the new XTMC43 set) and results from periodic DFT calculations shows that an agreement to within 9% can be achieved using GGA or mGGA functionals with atom-centred Gaussian-type basis functions. The need for careful assessments of consistency, calibration and reproducibility of experimental and computational data is discussed. Results regarding the new XTMC43 benchmark set are suggested to serve as a starting point for further method development, systematic screening and crystal engineering.
Collapse
Affiliation(s)
- José Eduardo Zamudio Díaz Mirón
- Molecular Simulations and Design Group, Max Planck Institute for Dynamics of Complex Technical Systems, Sandtorstrasse 1, 39106 Magdeburg, Germany.
| | - Matthias Stein
- Molecular Simulations and Design Group, Max Planck Institute for Dynamics of Complex Technical Systems, Sandtorstrasse 1, 39106 Magdeburg, Germany.
| |
Collapse
|
7
|
Chen YF, Schroeder C, Lew CM, Zones SI, Koller H, Sierka M. Cooperativity of silanol defect chemistry in zeolites. Phys Chem Chem Phys 2022; 25:478-485. [PMID: 36477757 DOI: 10.1039/d2cp05218k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022]
Abstract
Deboronation treatment of zeolite B-SSZ-55 can generate vacancy defects consisting of four silanol groups (silanol nests). However, 1H solid-state NMR spectroscopy indicates the prevalence of two silanol groups (silanol dyads) instead of four silanol groups. Such silanol dyads must be formed by the silanol condensation of two silanol groups at the silanol nests. Yet, the exact mechanism of this condensation and detailed structure of the silanol defect are not known. Here, the structure and formation mechanism of silanol dyads in the SSZ-55 zeolite have been investigated by both cluster and periodic density functional theory calculations. The calculated 1H NMR chemical shifts agree with the experimental values, showing that the silanol dyads are indeed commonly present at the vacancies and the vacancy density plays a role in the relaxation of the zeolite framework. The nature (size) of the silanol clusters influences their acidity.
Collapse
Affiliation(s)
- Ya-Fan Chen
- Friedrich Schiller University Jena, Otto Schott Institute of Materials Research, Löbdergraben 32, 07743 Jena, Germany.
| | - Christian Schroeder
- Institute of Physical Chemistry, University of Münster, Corrensstr. 30, 48149 Münster, Germany
| | | | | | - Hubert Koller
- Institute of Physical Chemistry, University of Münster, Corrensstr. 30, 48149 Münster, Germany
| | - Marek Sierka
- Friedrich Schiller University Jena, Otto Schott Institute of Materials Research, Löbdergraben 32, 07743 Jena, Germany.
| |
Collapse
|
8
|
Sharma S, Beylkin G. Efficient Evaluation of Two-Center Gaussian Integrals in Periodic Systems. J Chem Theory Comput 2021; 17:3916-3922. [PMID: 34061523 DOI: 10.1021/acs.jctc.0c01195] [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/28/2022]
Abstract
By using Poisson's summation formula, we calculate periodic integrals over Gaussian basis functions by partitioning the lattice summations between the real and reciprocal space, where both sums converge exponentially fast with a large exponent. We demonstrate that the summation can be performed efficiently to calculate two-center Gaussian integrals over various kernels including overlap, kinetic, and Coulomb. The summation in real space is performed using an efficient flavor of the McMurchie-Davidson recurrence relation. The expressions for performing summation in the reciprocal space are also derived and implemented. The algorithm for reciprocal space summation allows us to reuse several terms and leads to a significant improvement in efficiency when highly contracted basis functions with large exponents are used. We find that the resulting algorithm is only between a factor of 5 and 15 slower than that for molecular integrals, indicating the very small number of terms needed in both the real and reciprocal space summations. An outline of the algorithm for calculating three-center Coulomb integrals is also provided.
Collapse
Affiliation(s)
- Sandeep Sharma
- Department of Chemistry, University of Colorado, Boulder, Boulder, Colorado 80302, United States
| | - Gregory Beylkin
- Department of Applied Mathematics, University of Colorado, Boulder, Boulder, Colorado 80309, United States
| |
Collapse
|
9
|
Mewes J, Hansen A, Grimme S. Comment on “The Nature of Chalcogen‐Bonding‐Type Tellurium–Nitrogen Interactions”: Fixing the Description of Finite‐Temperature Effects Restores the Agreement Between Experiment and Theory. Angew Chem Int Ed Engl 2021. [DOI: 10.1002/ange.202102679] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Jan‐Michael Mewes
- Mulliken Center for Theoretical Chemistry 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 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 Institut für Physikalische und Theoretische Chemie Rheinische Friedrich-Wilhelms Universität Bonn Beringstraße 4 53115 Bonn Germany
| |
Collapse
|
10
|
Mewes J, Hansen A, Grimme S. Comment on "The Nature of Chalcogen-Bonding-Type Tellurium-Nitrogen Interactions": Fixing the Description of Finite-Temperature Effects Restores the Agreement Between Experiment and Theory. Angew Chem Int Ed Engl 2021; 60:13144-13149. [PMID: 33960596 PMCID: PMC8252449 DOI: 10.1002/anie.202102679] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2021] [Indexed: 12/14/2022]
Abstract
Mitzel and co-workers recently presented an intriguing molecule displaying a tellurium-nitrogen interaction. Structural data obtained in the solid and in gas phase indicated a large increase of the Te-N equilibrium distance re from 2.64 to 2.92 Å, respectively. Although some DFT calculations appear to support the large re in gas phase, we argue that the lions share of the increase is due to an incomplete description of finite-temperature effects in the back-corrected experimental data. This hypothesis is based on high-level coupled-cluster (CC) and periodic DFT calculations, which consistently point towards a much smaller re in the isolated molecule. Further support comes through MD simulations with a tuned GFN2-xTB Hamiltonian: Calibrated against a CC reference, these show a six-times larger influence of temperature than with the originally used GFN1-xTB. Taking this into account, the back-corrected re in gas phase becomes 2.67±0.08 Å, in good agreement with high-level CC theory and most DFT methods.
Collapse
Affiliation(s)
- Jan‐Michael Mewes
- Mulliken Center for Theoretical ChemistryInstitut für Physikalische und Theoretische ChemieRheinische Friedrich-Wilhelms Universität BonnBeringstraße 453115BonnGermany
| | - Andreas Hansen
- Mulliken Center for Theoretical ChemistryInstitut für Physikalische und Theoretische ChemieRheinische Friedrich-Wilhelms Universität BonnBeringstraße 453115BonnGermany
| | - Stefan Grimme
- Mulliken Center for Theoretical ChemistryInstitut für Physikalische und Theoretische ChemieRheinische Friedrich-Wilhelms Universität BonnBeringstraße 453115BonnGermany
| |
Collapse
|
11
|
Ye HZ, Berkelbach TC. Fast periodic Gaussian density fitting by range separation. J Chem Phys 2021; 154:131104. [PMID: 33832275 DOI: 10.1063/5.0046617] [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/15/2022] Open
Abstract
We present an efficient implementation of periodic Gaussian density fitting (GDF) using the Coulomb metric. The three-center integrals are divided into two parts by range-separating the Coulomb kernel, with the short-range part evaluated in real space and the long-range part in reciprocal space. With a few algorithmic optimizations, we show that this new method-which we call range-separated GDF (RSGDF)-scales sublinearly to linearly with the number of k-points for small to medium-sized k-point meshes that are commonly used in periodic calculations with electron correlation. Numerical results on a few three-dimensional solids show about ten-fold speedups over the previously developed GDF with little precision loss. The error introduced by RSGDF is about 10-5Eh in the converged Hartree-Fock energy with default auxiliary basis sets and can be systematically reduced by increasing the size of the auxiliary basis with little extra work.
Collapse
Affiliation(s)
- Hong-Zhou Ye
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | | |
Collapse
|
12
|
Fink K, Höfener S. Combining wavefunction frozen-density embedding with one-dimensional periodicity. J Chem Phys 2021; 154:104114. [PMID: 33722017 DOI: 10.1063/5.0041501] [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/14/2022] Open
Abstract
We present the combination of wavefunction frozen-density embedding (FDE) with a periodic repetition in one dimension (1D) for molecular systems in the KOALA program. In this periodic orbital-uncoupled FDE ansatz, no wavefunction overlap is taken into account, and only the electron density of the active subsystem is computed explicitly. This density is relaxed in the presence of the environment potential, which is obtained by translating the updated active subsystem density, yielding a fully self-consistent solution at convergence. Treating only one subsystem explicitly, the method allows for the calculation of local properties in condensed molecular systems, while no orbital band structure is obtained preventing the application, e.g., to systems with metallic bonding. In order to illustrate possible applications of the new implementation, selected case studies are presented, ranging from ground-state dipole moments using configuration interaction methods via excitation energies using time-dependent density-functional theory to ionization potentials obtained from equation-of-motion correlation methods. Different levels of approximations are assessed, revealing that an active subsystem consisting of two or three molecules leads to results that are converged with respect to the environment contributions.
Collapse
Affiliation(s)
- Karin Fink
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), P.O. Box 3630, 76021 Karlsruhe, Germany
| | - Sebastian Höfener
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), P.O. Box 6980, 76049 Karlsruhe, Germany
| |
Collapse
|
13
|
Grimme S, Hansen A, Ehlert S, Mewes JM. r 2SCAN-3c: A "Swiss army knife" composite electronic-structure method. J Chem Phys 2021; 154:064103. [PMID: 33588555 DOI: 10.1063/5.0040021] [Citation(s) in RCA: 280] [Impact Index Per Article: 93.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
The recently proposed r2SCAN meta-generalized-gradient approximation (mGGA) of Furness and co-workers is used to construct an efficient composite electronic-structure method termed r2SCAN-3c. To this end, the unaltered r2SCAN functional is combined with a tailor-made triple-ζ Gaussian atomic orbital basis set as well as with refitted D4 and geometrical counter-poise corrections for London-dispersion and basis set superposition error. The performance of the new method is evaluated for the GMTKN55 database covering large parts of chemical space with about 1500 data points, as well as additional benchmarks for non-covalent interactions, organometallic reactions, and lattice energies of organic molecules and ices, as well as for the adsorption on polar salt and non-polar coinage-metal surfaces. These comprehensive tests reveal a spectacular performance and robustness of r2SCAN-3c: It by far surpasses its predecessor B97-3c at only twice the cost and provides one of the best results of all semi-local density-functional theory (DFT)/QZ methods ever tested for the GMTKN55 database at one-tenth of the cost. Specifically, for reaction and conformational energies as well as non-covalent interactions, it outperforms prominent hybrid-DFT/QZ approaches at two to three orders of magnitude lower cost. Perhaps, the most relevant remaining issue of r2SCAN-3c is self-interaction error (SIE), owing to its mGGA nature. However, SIE is slightly reduced compared to other (m)GGAs, as is demonstrated in two examples. After all, this remarkably efficient and robust method is chosen as our new group default, replacing previous composite DFT and partially even expensive high-level methods in most standard applications for systems with up to several hundreds of atoms.
Collapse
Affiliation(s)
- Stefan Grimme
- Mulliken Center for Theoretical Chemistry, 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, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Sebastian Ehlert
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Jan-Michael Mewes
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| |
Collapse
|
14
|
Maschita J, Banerjee T, Savasci G, Haase F, Ochsenfeld C, Lotsch BV. Ionothermal Synthesis of Imide-Linked Covalent Organic Frameworks. Angew Chem Int Ed Engl 2020; 59:15750-15758. [PMID: 32573890 PMCID: PMC7497034 DOI: 10.1002/anie.202007372] [Citation(s) in RCA: 93] [Impact Index Per Article: 23.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2020] [Indexed: 11/26/2022]
Abstract
Covalent organic frameworks (COFs) are an extensively studied class of porous materials, which distinguish themselves from other porous polymers in their crystallinity and high degree of modularity, enabling a wide range of applications. COFs are most commonly synthesized solvothermally, which is often a time-consuming process and restricted to well-soluble precursor molecules. Synthesis of polyimide-linked COFs (PI-COFs) is further complicated by the poor reversibility of the ring-closing reaction under solvothermal conditions. Herein, we report the ionothermal synthesis of crystalline and porous PI-COFs in zinc chloride and eutectic salt mixtures. This synthesis does not require soluble precursors and the reaction time is significantly reduced as compared to standard solvothermal synthesis methods. In addition to applying the synthesis to previously reported imide COFs, a new perylene-based COF was also synthesized, which could not be obtained by the classical solvothermal route. In situ high-temperature XRPD analysis hints to the formation of precursor-salt adducts as crystalline intermediates, which then react with each other to form the COF.
Collapse
Affiliation(s)
- Johannes Maschita
- Nanochemistry DepartmentMax Planck Institute for Solid State ResearchHeisenbergstraße 170569StuttgartGermany
- Department of ChemistryUniversity of Munich (LMU)Butenandtstraße 5–1381377MünchenGermany
| | - Tanmay Banerjee
- Nanochemistry DepartmentMax Planck Institute for Solid State ResearchHeisenbergstraße 170569StuttgartGermany
| | - Gökcen Savasci
- Nanochemistry DepartmentMax Planck Institute for Solid State ResearchHeisenbergstraße 170569StuttgartGermany
- Department of ChemistryUniversity of Munich (LMU)Butenandtstraße 5–1381377MünchenGermany
| | - Frederik Haase
- Nanochemistry DepartmentMax Planck Institute for Solid State ResearchHeisenbergstraße 170569StuttgartGermany
- Department of ChemistryUniversity of Munich (LMU)Butenandtstraße 5–1381377MünchenGermany
- Current address: Institute for Functional InterfacesKarlsruhe Institute of Technology (KIT)Hermann-von-Helmholtz-Platz 176344Eggenstein-LeopoldshafenGermany
| | - Christian Ochsenfeld
- Nanochemistry DepartmentMax Planck Institute for Solid State ResearchHeisenbergstraße 170569StuttgartGermany
- Department of ChemistryUniversity of Munich (LMU)Butenandtstraße 5–1381377MünchenGermany
- E-conversion and Center for NanoscienceLichtenbergstraße 4a85748Garching bei MünchenGermany
| | - Bettina V. Lotsch
- Nanochemistry DepartmentMax Planck Institute for Solid State ResearchHeisenbergstraße 170569StuttgartGermany
- Department of ChemistryUniversity of Munich (LMU)Butenandtstraße 5–1381377MünchenGermany
- E-conversion and Center for NanoscienceLichtenbergstraße 4a85748Garching bei MünchenGermany
| |
Collapse
|
15
|
Xing X, Huang H, Chow E. A linear scaling hierarchical block low-rank representation of the electron repulsion integral tensor. J Chem Phys 2020; 153:084119. [DOI: 10.1063/5.0010732] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Xin Xing
- School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-4017, USA
| | - Hua Huang
- School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-4017, USA
| | - Edmond Chow
- School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-4017, USA
| |
Collapse
|
16
|
|
17
|
Gottschling K, Savasci G, Vignolo-González H, Schmidt S, Mauker P, Banerjee T, Rovó P, Ochsenfeld C, Lotsch BV. Rational Design of Covalent Cobaloxime-Covalent Organic Framework Hybrids for Enhanced Photocatalytic Hydrogen Evolution. J Am Chem Soc 2020; 142:12146-12156. [PMID: 32564604 PMCID: PMC7366382 DOI: 10.1021/jacs.0c02155] [Citation(s) in RCA: 71] [Impact Index Per Article: 17.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
![]()
Covalent
organic frameworks (COFs) display a unique combination
of chemical tunability, structural diversity, high porosity, nanoscale
regularity, and thermal stability. Recent efforts are directed at
using such frameworks as tunable scaffolds for chemical reactions.
In particular, COFs have emerged as viable platforms for mimicking
natural photosynthesis. However, there is an indisputable need for
efficient, stable, and economical alternatives for the traditional
platinum-based cocatalysts for light-driven hydrogen evolution. Here,
we present azide-functionalized chloro(pyridine)cobaloxime hydrogen-evolution
cocatalysts immobilized on a hydrazone-based COF-42 backbone that
show improved and prolonged photocatalytic activity with respect to
equivalent physisorbed systems. Advanced solid-state NMR and quantum-chemical
methods allow us to elucidate details of the improved photoreactivity
and the structural composition of the involved active site. We found
that a genuine interaction between the COF backbone and the cobaloxime
facilitates recoordination of the cocatalyst during the photoreaction,
thereby improving the reactivity and hindering degradation of the
catalyst. The excellent stability and prolonged reactivity make the
herein reported cobaloxime-tethered COF materials promising hydrogen
evolution catalysts for future solar fuel technologies.
Collapse
Affiliation(s)
- Kerstin Gottschling
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany.,Department of Chemistry, University of Munich (LMU), Butenandtstrasse 5-13, 81377 Munich, Germany.,Cluster of Excellence e-conversion, Lichtenbergstrasse 4a, 85748 Garching, Germany.,Center for Nanoscience (CeNS), Schellingstrasse 4, 80799 Munich, Germany
| | - Gökcen Savasci
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany.,Department of Chemistry, University of Munich (LMU), Butenandtstrasse 5-13, 81377 Munich, Germany.,Cluster of Excellence e-conversion, Lichtenbergstrasse 4a, 85748 Garching, Germany.,Center for Nanoscience (CeNS), Schellingstrasse 4, 80799 Munich, Germany
| | - Hugo Vignolo-González
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany
| | - Sandra Schmidt
- Department of Chemistry, University of Munich (LMU), Butenandtstrasse 5-13, 81377 Munich, Germany
| | - Philipp Mauker
- Department of Chemistry, University of Munich (LMU), Butenandtstrasse 5-13, 81377 Munich, Germany.,Center for Nanoscience (CeNS), Schellingstrasse 4, 80799 Munich, Germany
| | - Tanmay Banerjee
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany
| | - Petra Rovó
- Department of Chemistry, University of Munich (LMU), Butenandtstrasse 5-13, 81377 Munich, Germany.,Center for Nanoscience (CeNS), Schellingstrasse 4, 80799 Munich, Germany
| | - Christian Ochsenfeld
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany.,Department of Chemistry, University of Munich (LMU), Butenandtstrasse 5-13, 81377 Munich, Germany.,Cluster of Excellence e-conversion, Lichtenbergstrasse 4a, 85748 Garching, Germany.,Center for Nanoscience (CeNS), Schellingstrasse 4, 80799 Munich, Germany
| | - Bettina V Lotsch
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany.,Department of Chemistry, University of Munich (LMU), Butenandtstrasse 5-13, 81377 Munich, Germany.,Cluster of Excellence e-conversion, Lichtenbergstrasse 4a, 85748 Garching, Germany.,Center for Nanoscience (CeNS), Schellingstrasse 4, 80799 Munich, Germany
| |
Collapse
|
18
|
Balasubramani SG, Chen GP, Coriani S, Diedenhofen M, Frank MS, Franzke YJ, Furche F, Grotjahn R, Harding ME, Hättig C, Hellweg A, Helmich-Paris B, Holzer C, Huniar U, Kaupp M, Marefat Khah A, Karbalaei Khani S, Müller T, Mack F, Nguyen BD, Parker SM, Perlt E, Rappoport D, Reiter K, Roy S, Rückert M, Schmitz G, Sierka M, Tapavicza E, Tew DP, van Wüllen C, Voora VK, Weigend F, Wodyński A, Yu JM. TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations. J Chem Phys 2020; 152:184107. [PMID: 32414256 PMCID: PMC7228783 DOI: 10.1063/5.0004635] [Citation(s) in RCA: 563] [Impact Index Per Article: 140.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2020] [Accepted: 04/07/2020] [Indexed: 01/30/2023] Open
Abstract
TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy-cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe-Salpeter methods, second-order Møller-Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE's functionality, including excited-state methods, RPA and Green's function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE's current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE's development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.
Collapse
Affiliation(s)
- Sree Ganesh Balasubramani
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Guo P Chen
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Sonia Coriani
- DTU Chemistry, Technical University of Denmark, Kemitorvet Build. 207, DK-2800 Kongens Lyngby, Denmark
| | - Michael Diedenhofen
- Dassault Systèmes Deutschland GmbH, Imbacher Weg 46, 51379 Leverkusen, Germany
| | - Marius S Frank
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - Yannick J Franzke
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), KIT Campus South, P.O. Box 6980, 76049 Karlsruhe, Germany
| | - Filipp Furche
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Robin Grotjahn
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, 10623 Berlin, Germany
| | | | - Christof Hättig
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - Arnim Hellweg
- Dassault Systèmes Deutschland GmbH, Imbacher Weg 46, 51379 Leverkusen, Germany
| | - Benjamin Helmich-Paris
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Christof Holzer
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), KIT Campus South, P.O. Box 6980, 76049 Karlsruhe, Germany
| | - Uwe Huniar
- Dassault Systèmes Deutschland GmbH, Imbacher Weg 46, 51379 Leverkusen, Germany
| | - Martin Kaupp
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Alireza Marefat Khah
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | | | - Thomas Müller
- Forschungszentrum Jülich, Jülich Supercomputer Centre, Wilhelm-Jonen Straße, 52425 Jülich, Germany
| | - Fabian Mack
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), KIT Campus South, P.O. Box 6980, 76049 Karlsruhe, Germany
| | - Brian D Nguyen
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Shane M Parker
- Department of Chemistry, Case Western Reserve University, 10900 Euclid Avenue, Cleveland, Ohio 44106, USA
| | - Eva Perlt
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Dmitrij Rappoport
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, USA
| | - Kevin Reiter
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), KIT Campus North, P.O. Box 3640, 76021 Karlsruhe, Germany
| | - Saswata Roy
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Matthias Rückert
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - Gunnar Schmitz
- Department of Chemistry, Aarhus Universitet, Langelandsgade 140, DK-8000 Aarhus, Denmark
| | - Marek Sierka
- TURBOMOLE GmbH, Litzenhardtstraße 19, 76135 Karlsruhe, Germany
| | - Enrico Tapavicza
- Department of Chemistry and Biochemistry, California State University, Long Beach, 1250 Bellflower Boulevard, Long Beach, California 90840, USA
| | - David P Tew
- Max Planck Institute for Solid State Research, Heisenbergstaße 1, 70569 Stuttgart, Germany
| | - Christoph van Wüllen
- Fachbereich Chemie and Forschungszentrum OPTIMAS, Technische Universität Kaiserslautern, Erwin-Schrödinger-Staße 52, 67663 Kaiserslautern, Germany
| | - Vamsee K Voora
- Department of Chemical Sciences, Tata Institute of Fundamental Research, Homi Bhabha Road, Colaba, Mumbai 400005, India
| | - Florian Weigend
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT), KIT Campus North, P.O. Box 3640, 76021 Karlsruhe, Germany
| | - Artur Wodyński
- Institut für Chemie, Theoretische Chemie/Quantenchemie, Technische Universität Berlin, Sekr. C7, Straße des 17. Juni 135, 10623 Berlin, Germany
| | - Jason M Yu
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| |
Collapse
|
19
|
Intermolecular Interactions in Molecular Organic Crystals upon Relaxation of Lattice Parameters. CRYSTALS 2019. [DOI: 10.3390/cryst9120665] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
Crystal structure prediction is based on the assumption that the most thermodynamically stable structure will crystallize first. The existence of other structures such as polymorphs or from counterenantiomers requires an accurate calculation of the electronic energy. Using atom-centered Gaussian basis functions in periodic Density Functional Theory (DFT) calculations in Turbomole, the performance of two dispersion-corrected functionals, PBE-D3 and B97-D, is assessed for molecular organic crystals of the X23 benchmark set. B97-D shows a MAE (mean absolute error) of 4 kJ/mol, compared to 9 kJ/mol for PBE-D3. A strategy for the convergence of lattice energies towards the basis set limit is outlined. A simultaneous minimization of molecular structures and lattice parameters shows that both methods are able to reproduce experimental unit cell parameters to within 4–5%. Calculated lattice energies, however, deviate slightly more from the experiment, i.e., by 0.4 kJ/mol after unit cell optimization for PBE-D3 and 0.5 kJ/mol for B97-D. The accuracy of the calculated lattice energies compared to the experimental values demonstrates the ability of current DFT methods to assist in the quest for possible polymorphs and enantioselective crystallization processes.
Collapse
|
20
|
Matczak P, Mlostoń G, Hamera-Fałdyga R, Görls H, Weigand W. Structure of Diferrocenyl Thioketone: From Molecule to Crystal. MOLECULES (BASEL, SWITZERLAND) 2019; 24:molecules24213950. [PMID: 31683693 PMCID: PMC6864675 DOI: 10.3390/molecules24213950] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/12/2019] [Revised: 10/28/2019] [Accepted: 10/29/2019] [Indexed: 01/07/2023]
Abstract
Ferrocenyl-functionalized thioketones have recently been recognized as useful building blocks for sulfur-containing compounds with potential applications in materials chemistry. This work is devoted to a single representative of such thioketones, namely diferrocenyl thioketone (Fc2CS), whose structure has been determined here for the first time. Both X-ray crystallography and a wide variety of quantum-chemical methods were used to explore the structure of Fc2CS. In addition to the X-ray structure determination, intermolecular interactions occurring in the crystal structure of Fc2CS were examined in detail by quantum-chemical methods. These methods were also an invaluable tool in studying the molecular structure of Fc2CS, from the gas phase to solutions and to its crystal. Intramolecular interactions governing the conformational behavior of an isolated Fc2CS molecule were deduced from quantum-chemical analyses carried out in orbital space and real space. Our experimental and theoretical results indicate that the main structural features of an isolated Fc2CS molecule in its lowest-energy geometry are retained both upon solvation and in the crystal. The tilt of ferrocenyl groups is only slightly affected by crystal packing forces that are dominated by dispersion. Nonetheless, a network of intermolecular interactions, such as H···H, C···H and S···H, was detected in the Fc2CS crystal but each of them is fairly weak.
Collapse
Affiliation(s)
- Piotr Matczak
- Department of Physical Chemistry, Faculty of Chemistry, University of Lodz, Pomorska 163/165, 90236 Lodz, Poland.
| | - Grzegorz Mlostoń
- Department of Organic and Applied Chemistry, Faculty of Chemistry, University of Lodz, Tamka 12, 91403 Lodz, Poland.
| | - Róża Hamera-Fałdyga
- Department of Organic and Applied Chemistry, Faculty of Chemistry, University of Lodz, Tamka 12, 91403 Lodz, Poland.
| | - Helmar Görls
- Institute of Inorganic and Analytical Chemistry, Friedrich-Schiller-University Jena, Humboldtstrasse 8, 07743 Jena, Germany.
| | - Wolfgang Weigand
- Institute of Inorganic and Analytical Chemistry, Friedrich-Schiller-University Jena, Humboldtstrasse 8, 07743 Jena, Germany.
| |
Collapse
|
21
|
Atomistic Simulations of Plasma-Enhanced Atomic Layer Deposition. MATERIALS 2019; 12:ma12162605. [PMID: 31443331 PMCID: PMC6719897 DOI: 10.3390/ma12162605] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/25/2019] [Revised: 08/07/2019] [Accepted: 08/12/2019] [Indexed: 11/17/2022]
Abstract
Plasma-enhanced atomic layer deposition (PEALD) is a widely used, powerful layer-by-layer coating technology. Here, we present an atomistic simulation scheme for PEALD processes, combining the Monte Carlo deposition algorithm and structure relaxation using molecular dynamics. In contrast to previous implementations, our approach employs a real, atomistic model of the precursor. This allows us to account for steric hindrance and overlap restrictions at the surface corresponding to the real precursor deposition step. In addition, our scheme takes various process parameters into account, employing predefined probabilities for precursor products at each Monte Carlo deposition step. The new simulation protocol was applied to investigate PEALD synthesis of SiO2 thin films using the bis-diethylaminosilane precursor. It revealed that increasing the probability for precursor binding to one surface oxygen atom favors amorphous layer growth, a large number of –OH impurities, and the formation of voids. In contrast, a higher probability for precursor binding to two surface oxygen atoms leads to dense SiO2 film growth and a reduction of –OH impurities. Increasing the probability for the formation of doubly bonded precursor sites is therefore the key factor for the formation of dense SiO2 PEALD thin films with reduced amounts of voids and –OH impurities.
Collapse
|
22
|
Goncalves TJ, Plessow PN, Studt F. On the Accuracy of Density Functional Theory in Zeolite Catalysis. ChemCatChem 2019. [DOI: 10.1002/cctc.201900791] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Affiliation(s)
- Tiago J. Goncalves
- Institute of Catalysis Research and TechnologyKarlsruhe Institute of Technology Hermann-von-Helmholtz Platz 1 76344 Eggenstein-Leopoldshafen Germany
| | - Philipp N. Plessow
- Institute of Catalysis Research and TechnologyKarlsruhe Institute of Technology Hermann-von-Helmholtz Platz 1 76344 Eggenstein-Leopoldshafen Germany
| | - Felix Studt
- Institute of Catalysis Research and TechnologyKarlsruhe Institute of Technology Hermann-von-Helmholtz Platz 1 76344 Eggenstein-Leopoldshafen Germany
- Institute for Chemical Technology and Polymer ChemistryKarlsruhe Institute of Technology Engesserstrasse 18 76131 Karlsruhe Germany
| |
Collapse
|
23
|
Becker M, Sierka M. Density functional theory for molecular and periodic systems using density fitting and continuous fast multipole method: Stress tensor. J Comput Chem 2019; 40:2563-2570. [PMID: 31322769 DOI: 10.1002/jcc.26033] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2019] [Revised: 06/24/2019] [Accepted: 06/27/2019] [Indexed: 11/06/2022]
Abstract
A full implementation of the analytical stress tensor for periodic systems is reported in the TURBOMOLE program package within the framework of Kohn-Sham density functional theory using Gaussian-type orbitals as basis functions. It is the extension of the implementation of analytical energy gradients (Lazarski et al., Journal of Computational Chemistry 2016, 37, 2518-2526) to the stress tensor for the purpose of optimization of lattice vectors. Its key component is the efficient calculation of the Coulomb contribution by combining density fitting approximation and continuous fast multipole method. For the exchange-correlation (XC) part the hierarchical numerical integration scheme (Burow and Sierka, Journal of Chemical Theory and Computation 2011, 7, 3097-3104) is extended to XC weight derivatives and stress tensor. The computational efficiency and favorable scaling behavior of the stress tensor implementation are demonstrated for various model systems. The overall computational effort for energy gradient and stress tensor for the largest systems investigated is shown to be at most two and a half times the computational effort for the Kohn-Sham matrix formation. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Martin Becker
- Otto-Schott-Institut für Materialforschung, Friedrich-Schiller-Universität Jena, Löbdergraben 32, Jena, D-07743, Germany
| | - Marek Sierka
- Otto-Schott-Institut für Materialforschung, Friedrich-Schiller-Universität Jena, Löbdergraben 32, Jena, D-07743, Germany
| |
Collapse
|
24
|
Topochemical conversion of an imine- into a thiazole-linked covalent organic framework enabling real structure analysis. Nat Commun 2018; 9:2600. [PMID: 29968723 PMCID: PMC6030076 DOI: 10.1038/s41467-018-04979-y] [Citation(s) in RCA: 158] [Impact Index Per Article: 26.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2017] [Accepted: 05/08/2018] [Indexed: 11/08/2022] Open
Abstract
Stabilization of covalent organic frameworks (COFs) by post-synthetic locking strategies is a powerful tool to push the limits of COF utilization, which are imposed by the reversible COF linkage. Here we introduce a sulfur-assisted chemical conversion of a two-dimensional imine-linked COF into a thiazole-linked COF, with full retention of crystallinity and porosity. This post-synthetic modification entails significantly enhanced chemical and electron beam stability, enabling investigation of the real framework structure at a high level of detail. An in-depth study by electron diffraction and transmission electron microscopy reveals a myriad of previously unknown or unverified structural features such as grain boundaries and edge dislocations, which are likely generic to the in-plane structure of 2D COFs. The visualization of such real structural features is key to understand, design and control structure-property relationships in COFs, which can have major implications for adsorption, catalytic, and transport properties of such crystalline porous polymers.
Collapse
|
25
|
Caldeweyher E, Brandenburg JG. Simplified DFT methods for consistent structures and energies of large systems. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:213001. [PMID: 29633964 DOI: 10.1088/1361-648x/aabcfb] [Citation(s) in RCA: 36] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
Kohn-Sham density functional theory (DFT) is routinely used for the fast electronic structure computation of large systems and will most likely continue to be the method of choice for the generation of reliable geometries in the foreseeable future. Here, we present a hierarchy of simplified DFT methods designed for consistent structures and non-covalent interactions of large systems with particular focus on molecular crystals. The covered methods are a minimal basis set Hartree-Fock (HF-3c), a small basis set screened exchange hybrid functional (HSE-3c), and a generalized gradient approximated functional evaluated in a medium-sized basis set (B97-3c), all augmented with semi-classical correction potentials. We give an overview on the methods design, a comprehensive evaluation on established benchmark sets for geometries and lattice energies of molecular crystals, and highlight some realistic applications on large organic crystals with several hundreds of atoms in the primitive unit cell.
Collapse
Affiliation(s)
- Eike Caldeweyher
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | | |
Collapse
|
26
|
Sachse T, Martínez TJ, Dietzek B, Presselt M. A program for automatically predicting supramolecular aggregates and its application to urea and porphin. J Comput Chem 2018; 39:763-772. [PMID: 29297589 DOI: 10.1002/jcc.25151] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2017] [Revised: 12/04/2017] [Accepted: 12/07/2017] [Indexed: 11/08/2022]
Abstract
Not only the molecular structure but also the presence or absence of aggregates determines many properties of organic materials. Theoretical investigation of such aggregates requires the prediction of a suitable set of diverse structures. Here, we present the open-source program EnergyScan for the unbiased prediction of geometrically diverse sets of small aggregates. Its bottom-up approach is complementary to existing ones by performing a detailed scan of an aggregate's potential energy surface, from which diverse local energy minima are selected. We crossvalidate this approach by predicting both literature-known and heretofore unreported geometries of the urea dimer. We also predict a diverse set of dimers of the less intensely studied case of porphin, which we investigate further using quantum chemistry. For several dimers, we find strong deviations from a reference absorption spectrum, which we explain using computed transition densities. This proof of principle clearly shows that EnergyScan successfully predicts aggregates exhibiting large structural and spectral diversity. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Torsten Sachse
- Friedrich Schiller University, Institute of Physical Chemistry, Helmholtzweg 4, 07743, Jena, Germany.,Leibniz Institute of Photonic Technology Jena (IPHT), Research Department Functional Interfaces, Albert-Einstein-Straße 9, Jena, 07745, Germany
| | - Todd J Martínez
- Stanford University, Department of Chemistry and the PULSE Institute, 333 Campus Drive, Stanford, California 94305.,SLAC National Accelerator Laboratory, 2575 Sand Hill Rd, Menlo Park, California, 94025
| | - Benjamin Dietzek
- Friedrich Schiller University, Institute of Physical Chemistry, Helmholtzweg 4, 07743, Jena, Germany.,Center for Energy and Environmental Chemistry Jena, Humboldtstraße 10, Jena, 07743, Germany
| | - Martin Presselt
- Leibniz Institute of Photonic Technology Jena (IPHT), Research Department Functional Interfaces, Albert-Einstein-Straße 9, Jena, 07745, Germany.,SciClus GmbH & Co. KG, Moritz-von-Rohr-Straße 1a, Jena, 07745, Germany
| |
Collapse
|
27
|
Buchholz HK, Stein M. Accurate lattice energies of organic molecular crystals from periodic turbomole calculations. J Comput Chem 2018; 39:1335-1343. [DOI: 10.1002/jcc.25205] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2017] [Revised: 02/14/2018] [Accepted: 02/15/2018] [Indexed: 01/16/2023]
Affiliation(s)
- Hannes Konrad Buchholz
- Physical and Chemical Foundations Group; Max Planck Institute for Dynamics of Complex Technical Systems, Sandtorstrasse 1; Magdeburg 39106 Germany
- Molecular Simulations and Design Group; Max Planck Institute for Dynamics of Complex Technical Systems, Sandtorstrasse 1; Magdeburg 39106 Germany
| | - Matthias Stein
- Molecular Simulations and Design Group; Max Planck Institute for Dynamics of Complex Technical Systems, Sandtorstrasse 1; Magdeburg 39106 Germany
| |
Collapse
|