1
|
Hu G, Liu P, Jensen L. Calculating Molecular Polarizabilities Using Exact Frozen Density Embedding with External Orthogonality. J Chem Theory Comput 2024. [PMID: 39105755 DOI: 10.1021/acs.jctc.4c00692] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/07/2024]
Abstract
Frozen density embedding (FDE) with freeze-thaw cycles is a formally exact embedding scheme. In practice, this method is limited to systems with small density overlaps when approximate nonadditive kinetic energy functionals are used. It has been shown that the use of approximate nonadditive kinetic energy functionals can be avoided when external orthogonality (EO) is enforced, and FDE can then generate exact results even for strongly overlapping subsystems. In this work, we present an implementation of exact FDEc-EO (coupled FDE TDDFT with EO) for the calculation of polarizabilities in the Amsterdam density functional program package. EO is enforced using the level-shift projection operator method, which ensures that orbitals between fragments are orthogonal. For pure functionals, we show that only the symmetric EO contributions to the induced density matrix are needed. This leads to a simplified implementation for the calculation of polarizability that can exactly reproduce the supermolecular TDDFT results. We further discuss the limitation of exact FDEc-EO in interpreting subsystem polarizabilities due to the nonunique partitioning of the total density. We show that this limitation is due to the fact that subsystem polarizability partitioning is dependent on how the subsystems are initially polarized. As supermolecular virtual orbitals are exactly reproduced, this dependence is attributed to the description of the occupied orbitals. In contrast, for excitations of subsystems that are localized within one subsystem, we show that the excitation energies are stable with respect to the order of polarization. This observation shows that impacts from the nonunique nature of exact FDE on subsystem properties can be minimized by better fragmentation of the supermolecular systems if the property is localized. For global properties like polarizability, this is not the case, and nonuniqueness remains independent of the fragmentation used.
Collapse
Affiliation(s)
- Gaohe Hu
- Department of Chemistry, The Pennsylvania State University, 104 Benkovic Building, University Park, Pennsylvania 16802, United States
| | - Pengchong Liu
- Department of Chemistry, The Pennsylvania State University, 104 Benkovic Building, University Park, Pennsylvania 16802, United States
| | - Lasse Jensen
- Department of Chemistry, The Pennsylvania State University, 104 Benkovic Building, University Park, Pennsylvania 16802, United States
| |
Collapse
|
2
|
Sundaram V, Baumeier B. Quantum-Quantum and Quantum-Quantum-Classical Schemes for Near-Gap Excitations with Projection-Based-Embedded GW-Bethe-Salpeter Equation. J Chem Theory Comput 2024; 20:5451-5465. [PMID: 38916411 PMCID: PMC11238541 DOI: 10.1021/acs.jctc.4c00163] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/26/2024]
Abstract
We present quantum-quantum and quantum-quantum-classical schemes based on many-body Green's functions theory in the GW approximation with the Bethe-Salpeter equation (GW-BSE) employing projection-based-embedding (PbE). Such approaches allow defining active and inactive subsystems of larger, complex molecular systems, with only the smaller active subsystem being explicitly treated by GW-BSE offering significant computational advantages. However, as PbE can modify the single-particle states in the Kohn-Sham (KS) ground state calculation and screening effects from the inactive region are not automatically included in GW-BSE, results from such PbE-GW-BSE calculations can deviate from a full-system reference. Here, we scrutinize in detail, e.g., the individual and combined effects of different choices of active regions, the influence of omitting the screening from the inactive region, and strategies for basis set truncation on frontier orbital and near-gap electron-hole excitation energies. As prototypical systems, we consider a diketopyrrolopyrrole bicyclic ring including side-chains, a polarity-sensitive dye (prodan) in aqueous environment, and a π-stacked dimer of benzene and tetracyanoethylene in water, respectively, covering a variety of excitation characters in molecular systems with complex chemical environments and photoinduced processes. Our results suggest that to obtain agreement of approximately 0.1 eV between near-gap excitation energies from embedded and full calculations, the active region should be chosen based on the Mulliken population of the full highest-occupied molecular orbital and that careful benchmarking should be done on the KS level before the actual GW-BSE steps when basis set truncation is used. We find that PbE-GW-BSE offers significant reductions in computation times and, more importantly, memory requirements, making calculations for considerably larger systems tractable.
Collapse
Affiliation(s)
- Vivek Sundaram
- Department of Mathematics and Computer Science, Eindhoven University of Technology, P.O. Box 513, 5600MB Eindhoven, The Netherlands
- Institute for Complex Molecular Systems, Eindhoven University of Technology, P.O. Box 513, 5600MB Eindhoven, The Netherlands
- Department of Applied Physics and Science Education, Eindhoven University of Technology, P.O. Box 513, 5600MB Eindhoven, The Netherlands
| | - Björn Baumeier
- Department of Mathematics and Computer Science, Eindhoven University of Technology, P.O. Box 513, 5600MB Eindhoven, The Netherlands
- Institute for Complex Molecular Systems, Eindhoven University of Technology, P.O. Box 513, 5600MB Eindhoven, The Netherlands
| |
Collapse
|
3
|
Giovannini T, Scavino M, Koch H. Time-Dependent Multilevel Density Functional Theory. J Chem Theory Comput 2024; 20:3601-3612. [PMID: 38648031 DOI: 10.1021/acs.jctc.4c00041] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2024]
Abstract
We present a novel three-layer approach based on multilevel density functional theory (MLDFT) and polarizable molecular mechanics to simulate the electronic excitations of chemical systems embedded in an external environment within the time-dependent DFT formalism. In our method, the electronic structure of a target system, the chromophore, is determined in the field of an embedded inactive layer, which is treated as frozen. Long-range interactions are described by employing the polarizable fluctuating charge (FQ) force field. The resulting MLDFT/FQ thus accurately describes both electrostatics (and polarization) and non-electrostatic target-environment interactions. The robustness and reliability of the approach are demonstrated by comparing our results with experimental data reported for various organic molecules in solution.
Collapse
Affiliation(s)
| | - Marco Scavino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Henrik Koch
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| |
Collapse
|
4
|
Lafiosca P, Rossi F, Egidi F, Giovannini T, Cappelli C. Multiscale Frozen Density Embedding/Molecular Mechanics Approach for Simulating Magnetic Response Properties of Solvated Systems. J Chem Theory Comput 2024; 20:266-279. [PMID: 38109486 PMCID: PMC10782454 DOI: 10.1021/acs.jctc.3c00850] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2023] [Revised: 11/06/2023] [Accepted: 11/10/2023] [Indexed: 12/20/2023]
Abstract
We present a three-layer hybrid quantum mechanical/quantum embedding/molecular mechanics approach for calculating nuclear magnetic resonance (NMR) shieldings and J-couplings of molecular systems in solution. The model is based on the frozen density embedding (FDE) and polarizable fluctuating charges (FQ) and fluctuating dipoles (FQFμ) force fields and permits the accurate ab initio description of short-range nonelectrostatic interactions by means of the FDE shell and cost-effective treatment of long-range electrostatic interactions through the polarizable force field FQ(Fμ). Our approach's accuracy and potential are demonstrated by studying NMR spectra of Brooker's merocyanine in aqueous and nonaqueous solutions.
Collapse
Affiliation(s)
- Piero Lafiosca
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Federico Rossi
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Franco Egidi
- Software
for Chemistry and Materials BV, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | | | - Chiara Cappelli
- Scuola
Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
5
|
Mi W, Luo K, Trickey SB, Pavanello M. Orbital-Free Density Functional Theory: An Attractive Electronic Structure Method for Large-Scale First-Principles Simulations. Chem Rev 2023; 123:12039-12104. [PMID: 37870767 DOI: 10.1021/acs.chemrev.2c00758] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2023]
Abstract
Kohn-Sham Density Functional Theory (KSDFT) is the most widely used electronic structure method in chemistry, physics, and materials science, with thousands of calculations cited annually. This ubiquity is rooted in the favorable accuracy vs cost balance of KSDFT. Nonetheless, the ambitions and expectations of researchers for use of KSDFT in predictive simulations of large, complicated molecular systems are confronted with an intrinsic computational cost-scaling challenge. Particularly evident in the context of first-principles molecular dynamics, the challenge is the high cost-scaling associated with the computation of the Kohn-Sham orbitals. Orbital-free DFT (OFDFT), as the name suggests, circumvents entirely the explicit use of those orbitals. Without them, the structural and algorithmic complexity of KSDFT simplifies dramatically and near-linear scaling with system size irrespective of system state is achievable. Thus, much larger system sizes and longer simulation time scales (compared to conventional KSDFT) become accessible; hence, new chemical phenomena and new materials can be explored. In this review, we introduce the historical contexts of OFDFT, its theoretical basis, and the challenge of realizing its promise via approximate kinetic energy density functionals (KEDFs). We review recent progress on that challenge for an array of KEDFs, such as one-point, two-point, and machine-learnt, as well as some less explored forms. We emphasize use of exact constraints and the inevitability of design choices. Then, we survey the associated numerical techniques and implemented algorithms specific to OFDFT. We conclude with an illustrative sample of applications to showcase the power of OFDFT in materials science, chemistry, and physics.
Collapse
Affiliation(s)
- Wenhui Mi
- Key Laboratory of Material Simulation Methods & Software of Ministry of Education, College of Physics, Jilin University, Changchun 130012, PR China
- State Key Laboratory of Superhard Materials, Jilin University, Changchun 130012, PR China
- International Center of Future Science, Jilin University, Changchun 130012, PR China
| | - Kai Luo
- Department of Applied Physics, Nanjing University of Science and Technology, Nanjing 210094, PR China
| | - S B Trickey
- Quantum Theory Project, Department of Physics and Department of Chemistry, University of Florida, Gainesville, Florida 32611, United States
| | - Michele Pavanello
- Department of Physics and Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| |
Collapse
|
6
|
Kolodzeiski E, Stein CJ. Automated, Consistent, and Even-Handed Selection of Active Orbital Spaces for Quantum Embedding. J Chem Theory Comput 2023; 19:6643-6655. [PMID: 37775093 PMCID: PMC10569175 DOI: 10.1021/acs.jctc.3c00653] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2023] [Indexed: 10/01/2023]
Abstract
A widely used strategy to reduce the computational cost of quantum-chemical calculations is to partition the system into an active subsystem, which is the focus of the computational efforts, and an environment that is treated at a lower computational level. The system partitioning is mostly based on localized molecular orbitals. When reaction paths or energy differences are to be calculated, it is crucial to keep the orbital space consistent for all structures. Inconsistencies in orbital space can lead to unpredictable errors on the potential energy surface. While successful strategies to ensure this consistency have been established for organic and even metal-organic systems, these methods often fail for metal clusters or nanoparticles with a high density of near-degenerate and delocalized molecular orbitals. However, such systems are highly relevant for catalysis. Accurate yet feasible quantum-mechanical ab initio calculations are therefore highly desired. In this work, we present an approach based on the subsystem projected atomic orbital decomposition algorithm that allows us to ensure automated and consistent partitioning even for systems with delocalized and near-degenerate molecular orbitals and demonstrate the validity of this method for the binding energies of small molecules on transition-metal clusters.
Collapse
Affiliation(s)
- Elena Kolodzeiski
- Technical University of Munich, TUM
School of Natural Sciences, Department of Chemistry, Lichtenbergstr. 4, Garching D-85748, Germany
| | - Christopher J. Stein
- Technical University of Munich, TUM
School of Natural Sciences, Department of Chemistry, Lichtenbergstr. 4, Garching D-85748, Germany
| |
Collapse
|
7
|
Pauls M, Schnieders D, Dronskowski R. Embedded Localized Molecular-Orbital Representations for Periodic Wave Functions. J Phys Chem A 2023; 127:6541-6551. [PMID: 37523466 DOI: 10.1021/acs.jpca.3c02283] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/02/2023]
Abstract
To more straightforwardly provide local chemical-bonding reasoning in crystalline matter, we introduce a new approach to generate a real-space analogue of periodic electronic structures using "exact" top-down frozen-density embedding calculations. Based on the obtained real-space electronic structure, we then construct localized molecular orbitals and evidence that our technique compares favorably against the commonly used Wannier method, both in terms of numerical efficiency and details of chemical bonding. The new method has been implemented into the LOBSTER software package and designed as a black-box approach, digesting any periodic electronic structure from the currently supported codes, i.e., VASP, Quantum ESPRESSO, and ABINIT.
Collapse
Affiliation(s)
- Mike Pauls
- Institute of Inorganic Chemistry, RWTH Aachen University, Landoltweg 1, D-52056 Aachen, Germany
| | - David Schnieders
- Institute of Inorganic Chemistry, RWTH Aachen University, Landoltweg 1, D-52056 Aachen, Germany
| | - Richard Dronskowski
- Institute of Inorganic Chemistry, RWTH Aachen University, Landoltweg 1, D-52056 Aachen, Germany
- Hoffmann Institute of Advanced Materials, Shenzhen Polytechnic, 7098 Liuxian Blvd, Nanshan District, 518055 Shenzhen, China
| |
Collapse
|
8
|
Polak E, Englert T, Gander MJ, Wesolowski TA. Symmetrized non-decomposable approximations of the non-additive kinetic energy functional. J Chem Phys 2023; 158:2887765. [PMID: 37129139 DOI: 10.1063/5.0143602] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2023] [Accepted: 04/13/2023] [Indexed: 05/03/2023] Open
Abstract
In subsystem density functional theory (DFT), the bottom-up strategy to approximate the multivariable functional of the non-additive kinetic energy (NAKE) makes it possible to impose exact properties on the corresponding NAKE potential (NAKEP). Such a construction might lead to a non-symmetric and non-homogeneous functional, which excludes the use of such approximations for the evaluation of the total energy. We propose a general formalism to construct a symmetric version based on a perturbation theory approach of the energy expression for the asymmetric part. This strategy is then applied to construct a symmetrized NAKE corresponding to the NAKEP developed recently [Polak et al., J. Chem. Phys. 156, 044103 (2022)], making it possible to evaluate consistently the energy. These functionals were used to evaluate the interaction energy in several model intermolecular complexes using the formal framework of subsystem DFT. The new symmetrized energy expression shows a superior qualitative performance over common decomposable models.
Collapse
Affiliation(s)
- Elias Polak
- Department of Physical Chemistry, University of Geneva, Quai Ernest-Ansermet 30, CH-1211 Genève 4, Switzerland
- Section of Mathematics, University of Geneva, Rue du Conseil-Général 7-9, CP 64, CH-1211 Genève 4, Switzerland
| | - Tanguy Englert
- Department of Physical Chemistry, University of Geneva, Quai Ernest-Ansermet 30, CH-1211 Genève 4, Switzerland
| | - Martin J Gander
- Section of Mathematics, University of Geneva, Rue du Conseil-Général 7-9, CP 64, CH-1211 Genève 4, Switzerland
| | - Tomasz A Wesolowski
- Department of Physical Chemistry, University of Geneva, Quai Ernest-Ansermet 30, CH-1211 Genève 4, Switzerland
| |
Collapse
|
9
|
Niemeyer N, Eschenbach P, Bensberg M, Tölle J, Hellmann L, Lampe L, Massolle A, Rikus A, Schnieders D, Unsleber JP, Neugebauer J. The subsystem quantum chemistry program
Serenity. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1647] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Niklas Niemeyer
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Patrick Eschenbach
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Moritz Bensberg
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Johannes Tölle
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Lars Hellmann
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Lukas Lampe
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Anja Massolle
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Anton Rikus
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - David Schnieders
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| | - Jan P. Unsleber
- Laboratorium für Physikalische Chemie ETH Zürich Zürich Switzerland
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch‐Chemisches Institut and Center for Multiscale Theory and Computation Westfälische Wilhelms‐Universität Münster Münster Germany
| |
Collapse
|
10
|
Sharma M, Sierka M. Efficient Implementation of Density Functional Theory Based Embedding for Molecular and Periodic Systems Using Gaussian Basis Functions. J Chem Theory Comput 2022; 18:6892-6904. [PMID: 36223886 DOI: 10.1021/acs.jctc.2c00380] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A practical and effective implementation of density functional theory based embedding is reported, which allows us to treat both periodic and aperiodic systems on an equal footing. Its essence is the expansion of orbitals and electron density of the periodic system using Gaussian basis functions, rather than plane-waves, which provides a unique all-electron direct-space representation, thus avoiding the need for pseudopotentials. This makes the construction of embedding potential for a molecular active subsystem due to a periodic environment quite convenient, as transformation between representations is far from trivial. The three flavors of embedding, molecule-in-molecule, molecule-in-periodic, and periodic-in-periodic embedding, are implemented using embedding potentials based on non-additive kinetic energy density functionals (approximate) and level-shift projection operator (exact). The embedding scheme is coupled with a variety of correlated wave function theory (WFT) methods, thereby providing an efficient way to study the ground and excited state properties of low-dimensional systems using high-level methods for the region of interest. Finally, an implementation of real time-time-dependent density functional embedding theory (RT-TDDFET) is presented that uses a projection operator-based embedding potential and provides accurate results compared to full RT-TDDFT for systems with uncoupled excitations. The embedding potential is calculated efficiently using a combination of density fitting and continuous fast multipole method for the Coulomb term. The applicability of (i) WFT-in-DFT embedding, in predicting the adsorption and excitation energies, and (ii) RT-TDDFET, in predicting the absorption spectra, is explored for various test systems.
Collapse
Affiliation(s)
- Manas Sharma
- Otto Schott Institute of Materials Research, Friedrich Schiller University Jena, Löbdergraben 32, 07743Jena, Germany
| | - Marek Sierka
- Otto Schott Institute of Materials Research, Friedrich Schiller University Jena, Löbdergraben 32, 07743Jena, Germany
| |
Collapse
|
11
|
Eschenbach P, Neugebauer J. Subsystem density-functional theory: A reliable tool for spin-density based properties. J Chem Phys 2022; 157:130902. [PMID: 36209003 DOI: 10.1063/5.0103091] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Subsystem density-functional theory compiles a set of features that allow for efficiently calculating properties of very large open-shell radical systems such as organic radical crystals, proteins, or deoxyribonucleic acid stacks. It is computationally less costly than correlated ab initio wave function approaches and can pragmatically avoid the overdelocalization problem of Kohn-Sham density-functional theory without employing hard constraints on the electron-density. Additionally, subsystem density-functional theory calculations commonly start from isolated fragment electron densities, pragmatically preserving a priori specified subsystem spin-patterns throughout the calculation. Methods based on subsystem density-functional theory have seen a rapid development over the past years and have become important tools for describing open-shell properties. In this Perspective, we address open questions and possible developments toward challenging future applications in connection with subsystem density-functional theory for spin-dependent properties.
Collapse
Affiliation(s)
- Patrick Eschenbach
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| |
Collapse
|
12
|
De Santis M, Sorbelli D, Vallet V, Gomes ASP, Storchi L, Belpassi L. Frozen-Density Embedding for Including Environmental Effects in the Dirac-Kohn-Sham Theory: An Implementation Based on Density Fitting and Prototyping Techniques. J Chem Theory Comput 2022; 18:5992-6009. [PMID: 36172757 PMCID: PMC9558305 DOI: 10.1021/acs.jctc.2c00499] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Frozen density embedding (FDE) represents an embedding scheme in which environmental effects are included from first-principles calculations by considering the surrounding system explicitly by means of its electron density. In the present paper, we extend the full four-component relativistic Dirac-Kohn-Sham (DKS) method, as implemented in the BERTHA code, to include environmental and confinement effects with the FDE scheme (DKS-in-DFT FDE). The implementation, based on the auxiliary density fitting techniques, has been enormously facilitated by BERTHA's python API (PyBERTHA), which facilitates the interoperability with other FDE implementations available through the PyADF framework. The accuracy and numerical stability of this new implementation, also using different auxiliary fitting basis sets, has been demonstrated on the simple NH3-H2O system, in comparison with a reference nonrelativistic implementation. The computational performance has been evaluated on a series of gold clusters (Aun, with n = 2, 4, 8) embedded into an increasing number of water molecules (5, 10, 20, 40, and 80 water molecules). We found that the procedure scales approximately linearly both with the size of the frozen surrounding environment (consistent with the underpinnings of the FDE approach) and with the size of the active system (in line with the use of density fitting). Finally, we applied the code to a series of heavy (Rn) and super-heavy elements (Cn, Fl, Og) embedded in a C60 cage to explore the confinement effect induced by C60 on their electronic structure. We compare the results from our simulations, with respect to more-approximate models employed in the atomic physics literature. Our results indicate that the specific interactions described by FDE are able to improve upon the cruder approximations currently employed, and, thus, they provide a basis from which to generate more-realistic radial potentials for confined atoms.
Collapse
Affiliation(s)
- Matteo De Santis
- Univ. Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers Atomes et Molécules, F-59000 Lille, France
| | - Diego Sorbelli
- Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy.,Istituto di Scienze e Tecnologie Chimiche (SCITEC), Consiglio Nazionale delle Ricerche c/o Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
| | - Valérie Vallet
- Univ. Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers Atomes et Molécules, F-59000 Lille, France
| | | | - Loriano Storchi
- Istituto di Scienze e Tecnologie Chimiche (SCITEC), Consiglio Nazionale delle Ricerche c/o Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy.,Dipartimento di Farmacia, Università degli Studi 'G. D'Annunzio', Via dei Vestini 31, 66100 Chieti, Italy
| | - Leonardo Belpassi
- Istituto di Scienze e Tecnologie Chimiche (SCITEC), Consiglio Nazionale delle Ricerche c/o Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
| |
Collapse
|
13
|
Wesolowski TA. Is the non-additive kinetic potential always equal to the difference of effective potentials from inverting the Kohn-Sham equation? J Chem Phys 2022; 157:081102. [DOI: 10.1063/5.0101791] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The relation used frequently in the literature {according to which the non-additive kinetic potential which is a functional depending on a pair of electron densities is equal (up to a constant) to} the difference of two potentials obtained from inverting two Kohn-Sham equations, is examined. <p>The relation is based on a silent assumption that the two densities can be obtained from two independent Kohn-Sham equations - i.e. are $v_s$-representable.</p> <p>It is shown that this assumption does not hold for pairs of densities: $\rho_{tot}$ - being the Kohn-Sham density in some system and $\rho_B$ obtained from such partitioning of $\rho_{tot}$ that the difference $\rho_{tot}-\rho_B$ vanish on a Lebesgue measurable volume element. The inversion procedure is still applicable for $\rho_{tot}-\rho_B$ but cannot be interpreted as the inversion of Kohn-Sham equation but rather inversion of Kohn-Sham equation with an additional constraint. The obtained effective potential comprises a "contaminant" that might even not be unique. It is shown, that the construction of the non-additive kinetic potential based on the examined relation is not applicable for such pairs.
Collapse
Affiliation(s)
- Tomasz A. Wesolowski
- Department of Physical Chemistry, University of Geneva Faculty of Science, Switzerland
| |
Collapse
|
14
|
Miranda-Quintana RA, Heidar-Zadeh F, Fias S, Chapman AEA, Liu S, Morell C, Gómez T, Cárdenas C, Ayers PW. Molecular Interactions From the Density Functional Theory for Chemical Reactivity: The Interaction Energy Between Two-Reagents. Front Chem 2022; 10:906674. [PMID: 35769444 PMCID: PMC9234655 DOI: 10.3389/fchem.2022.906674] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2022] [Accepted: 04/19/2022] [Indexed: 12/13/2022] Open
Abstract
Reactivity descriptors indicate where a reagent is most reactive and how it is most likely to react. However, a reaction will only occur when the reagent encounters a suitable reaction partner. Determining whether a pair of reagents is well-matched requires developing reactivity rules that depend on both reagents. This can be achieved using the expression for the minimum-interaction-energy obtained from the density functional reactivity theory. Different terms in this expression will be dominant in different circumstances; depending on which terms control the reactivity, different reactivity indicators will be preferred.
Collapse
Affiliation(s)
- Ramón Alain Miranda-Quintana
- Department of Chemistry and Quantum Theory Project, University of Florida, Gainesville, FL, United States
- *Correspondence: Ramón Alain Miranda-Quintana, ; Carlos Cárdenas, ; Paul W. Ayers, ; Tatiana Gómez,
| | | | - Stijn Fias
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, ON, Canada
| | - Allison E. A. Chapman
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, ON, Canada
| | - Shubin Liu
- Research Computing Center, University of North Carolina, Chapel Hill, NC, United States
| | - Christophe Morell
- Université de Lyon, Université Claude Bernard Lyon 1, Institut des Sciences Analytiques-UMR CNRS 5280, Villeurbanne, France
| | - Tatiana Gómez
- Theoretical and Computational Chemistry Center, Institute of Applied Chemical Sciences, Faculty of Engineering, Universidad Autonoma de Chile, Santiago, Chile
- *Correspondence: Ramón Alain Miranda-Quintana, ; Carlos Cárdenas, ; Paul W. Ayers, ; Tatiana Gómez,
| | - Carlos Cárdenas
- Departamento de Fisica, Facultad de Ciencias, Universidad de Chile, Santiago, Chile
- Centro para el desarrollo de la Nanociencias y Nanotecnologia, CEDENNA, Santiago, Chile
- *Correspondence: Ramón Alain Miranda-Quintana, ; Carlos Cárdenas, ; Paul W. Ayers, ; Tatiana Gómez,
| | - Paul W. Ayers
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, ON, Canada
- *Correspondence: Ramón Alain Miranda-Quintana, ; Carlos Cárdenas, ; Paul W. Ayers, ; Tatiana Gómez,
| |
Collapse
|
15
|
Erhard J, Trushin E, Görling A. Numerically stable inversion approach to construct Kohn-Sham potentials for given electron densities within a Gaussian basis set framework. J Chem Phys 2022; 156:204124. [PMID: 35649824 DOI: 10.1063/5.0087356] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a Kohn-Sham (KS) inversion approach to construct KS exchange-correlation potentials corresponding to given electron densities. This method is based on an iterative procedure using linear response to update potentials. All involved quantities, i.e., orbitals, potentials, and response functions, are represented by Gaussian basis functions. In contrast to previous KS inversion methods relying on Gaussian basis sets, the method presented here is numerically stable even for standard basis sets from basis set libraries due to a preprocessing of the auxiliary basis used to represent an exchange-correlation charge density that generates the exchange-correlation potential. The new KS inversion method is applied to reference densities of various atoms and molecules obtained by full configuration interaction or CCSD(T) (coupled cluster singles doubles perturbative triples). The considered examples encompass cases known to be difficult, such as stretched hydrogen or lithium hydride molecules or the beryllium isoelectronic series. For the stretched hydrogen molecule, potentials of benchmark quality are obtained by employing large basis sets. For the carbon monoxide molecule, we show that the correlation potential from the random phase approximation (RPA) is in excellent qualitative and quantitative agreement with the correlation potential from the KS inversion of a CCSD(T) reference density. This indicates that RPA correlation potentials, in contrast to those from semi-local density-functionals, resemble the exact correlation potential. Besides providing exchange-correlation potentials for benchmark purposes, the proposed KS inversion method may be used in density-partition-based quantum embedding and in subsystem density-functional methods because it combines numerical stability with computational efficiency.
Collapse
Affiliation(s)
- Jannis Erhard
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Egor Trushin
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| | - Andreas Görling
- Lehrstuhl für Theoretische Chemie, Universität Erlangen-Nürnberg, Egerlandstr. 3, D-91058 Erlangen, Germany
| |
Collapse
|
16
|
Shi Y, Chávez VH, Wasserman A. n2v
: A density‐to‐potential inversion suite. A sandbox for creating, testing, and benchmarking density functional theory inversion methods. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1617] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Yuming Shi
- Department of Physics and Astronomy Purdue University West Lafayette Indiana USA
| | - Victor H. Chávez
- Department of Chemistry Purdue University West Lafayette Indiana USA
| | - Adam Wasserman
- Department of Physics and Astronomy Purdue University West Lafayette Indiana USA
- Department of Chemistry Purdue University West Lafayette Indiana USA
| |
Collapse
|
17
|
De Santis M, Vallet V, Gomes ASP. Environment Effects on X-Ray Absorption Spectra With Quantum Embedded Real-Time Time-Dependent Density Functional Theory Approaches. Front Chem 2022; 10:823246. [PMID: 35295974 PMCID: PMC8919347 DOI: 10.3389/fchem.2022.823246] [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: 11/27/2021] [Accepted: 02/04/2022] [Indexed: 11/13/2022] Open
Abstract
In this work we implement the real-time time-dependent block-orthogonalized Manby-Miller embedding (rt-BOMME) approach alongside our previously developed real-time frozen density embedding time-dependent density functional theory (rt-TDDFT-in-DFT FDE) code, and investigate these methods’ performance in reproducing X-ray absorption spectra (XAS) obtained with standard rt-TDDFT simulations, for model systems comprised of solvated fluoride and chloride ions ([X@(H2O)8−, X = F, Cl). We observe that for ground-state quantities such as core orbital energies, the BOMME approach shows significantly better agreement with supermolecular results than FDE for the strongly interacting fluoride system, while for chloride the two embedding approaches show more similar results. For the excited states, we see that while FDE (constrained not to have the environment densities relaxed in the ground state) is in good agreement with the reference calculations for the region around the K and L1 edges, and is capable of reproducing the splitting of the 1s1 (n + 1)p1 final states (n + 1 being the lowest virtual p orbital of the halides), it by and large fails to properly reproduce the 1s1 (n + 2)p1 states and misses the electronic states arising from excitation to orbitals with important contributions from the solvent. The BOMME results, on the other hand, provide a faithful qualitative representation of the spectra in all energy regions considered, though its intrinsic approximation of employing a lower-accuracy exchange-correlation functional for the environment induces non-negligible shifts in peak positions for the excitations from the halide to the environment. Our results thus confirm that QM/QM embedding approaches are viable alternatives to standard real-time simulations of X-ray absorption spectra of species in complex or confined environments.
Collapse
|
18
|
Abstract
Quantum embedding schemes are a promising way to extend multireference computations to large molecules with strong correlation effects localized on a small number of atoms. This work introduces a second-order active-space embedding theory [ASET(2)] which improves upon mean-field frozen embedding by treating fragment-environment interactions via an approximate canonical transformation. The canonical transformation employed in ASET(2) is formulated using the driven similarity renormalization group. The ASET(2) scheme is benchmarked on the N═N bond dissociation in pentyldiazene, the S0 to S1 excitation in 1-octene, and the interaction energy of the O2-benzene complex. The ASET(2) explicit treatment of fragment-environment interactions beyond the mean-field level generally improves the accuracy of embedded computations, and it becomes necessary to achieve an accurate description of excitation energies of 1-octene and the singlet-triplet gap of the O2-benzene complex.
Collapse
Affiliation(s)
- Nan He
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Chenyang Li
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
19
|
Tölle J, Neugebauer J. The Seamless Connection of Local and Collective Excited States in Subsystem Time-Dependent Density Functional Theory. J Phys Chem Lett 2022; 13:1003-1018. [PMID: 35061387 DOI: 10.1021/acs.jpclett.1c04023] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/21/2023]
Abstract
The theoretical understanding of photoinduced processes in multichromophoric systems requires, as an essential ingredient, the possibility of accurately describing their electronically excited states. However, the size of these systems often prohibits the usage of conventional electronic-structure methods, so that often multiscale approaches based on phenomenologically motivated models are employed. In contrast, subsystem time-dependent density functional theory (sTDDFT) allows for a subsystem-based ab initio description of multichromophoric systems and therefore allows for, in principle, an exact description of photoinduced processes. This Perspective aims to outline the theoretical foundations and commonly used practical realizations as well as to illustrate benefits of recent developments and open issues in the field of sTDDFT. Prospective, potential future applications and possible methodological developments are discussed.
Collapse
Affiliation(s)
- Johannes Tölle
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| |
Collapse
|
20
|
Treß RS, Hättig C, Höfener S. Employing Pseudopotentials to Tackle Excited-State Electron Spill-Out in Frozen Density Embedding Calculations. J Chem Theory Comput 2022; 18:1737-1747. [PMID: 35107998 DOI: 10.1021/acs.jctc.1c00732] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
In frozen density embedding (FDE), the properties of a target molecule are computed in the presence of an effective embedding potential, which accounts for the attractive and repulsive contributions of the environment. The formally exact embedding potential, however, is in practice calculated using explicit kinetic-energy functionals for which the resulting potentials are in many cases not repulsive enough to account fully for Pauli repulsion by the electrons of the environment and to compensate thereby the strong electron-nuclear attraction. For the excited states on the target molecule, this leads to charge spill-out when diffuse basis functions are included, which allow that valence electrons are excited to those regions of the environment where the strong nuclear attraction is not sufficiently compensated by repulsive contributions. To reduce this insufficiency, we propose in the present work the inclusion of atomic all-electron pseudopotentials for all environment atoms on top of the conventional embedding potential. In the current work, the pseudopotentials are applied for computing vertical excitation energies of local excited states in complex systems employing the second-order algebraic diagrammatic construction (ADC(2)) scheme. The proposed approach leads to significantly reduced charge spill-out and an improved agreement of FDE and supermolecular calculations in the frozen solvent approximation. In particular, when diffuse functions are employed, the mean absolute deviation (MAD) is reduced from 0.27 to 0.05 eV for the investigated cases.
Collapse
Affiliation(s)
- Robert S Treß
- Department of Theoretical Chemistry, Ruhr University Bochum, Bochum 44801, Germany
| | - Christof Hättig
- Department of Theoretical Chemistry, Ruhr University Bochum, Bochum 44801, Germany
| | - Sebastian Höfener
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Karlsruhe 76131, Germany
| |
Collapse
|
21
|
Eschenbach P, Artiukhin DG, Neugebauer J. Multi-state formulation of the frozen-density embedding quasi-diabatization approach. J Chem Phys 2021; 155:174104. [PMID: 34749503 DOI: 10.1063/5.0070486] [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 a multi-state implementation of the recently developed frozen-density embedding diabatization (FDE-diab) methodology [D. G. Artiukhin and J. Neugebauer, J. Chem. Phys. 148, 214104 (2018)] in the Serenity program. The new framework extends the original approach such that any number of charge-localized quasi-diabatic states can be coupled, giving an access to calculations of ground and excited state spin-density distributions as well as to excitation energies. We show that it is possible to obtain results similar to those from correlated wave function approaches such as the complete active space self-consistent field method at much lower computational effort. Additionally, we present a series of approximate computational schemes, which further decrease the overall computational cost and systematically converge to the full FDE-diab solution. The proposed methodology enables computational studies on spin-density distributions and related properties for large molecular systems of biochemical interest.
Collapse
Affiliation(s)
- Patrick Eschenbach
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Denis G Artiukhin
- Department of Chemistry, Aarhus Universitet, DK-8000 Aarhus, Denmark
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| |
Collapse
|
22
|
Giannone G, Śmiga S, D'Agostino S, Fabiano E, Della Sala F. Plasmon Couplings from Subsystem Time-Dependent Density Functional Theory. J Phys Chem A 2021; 125:7246-7259. [PMID: 34403247 DOI: 10.1021/acs.jpca.1c05384] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Many applications in plasmonics are related to the coupling between metallic nanoparticles (MNPs) or between an emitter and a MNP. The theoretical analysis of such a coupling is thus of fundamental importance to analyze the plasmonic behavior and to design new systems. While classical methods neglect quantum and spill-out effects, time-dependent density functional theory (TD-DFT) considers all of them and with Kohn-Sham orbitals delocalized over the whole system. Thus, within TD-DFT, no definite separation of the subsystems (the single MNP or the emitter) and their couplings is directly available. This important feature is obtained here using the subsystem formulation of TD-DFT, which has been originally developed in the context of weakly interacting organic molecules. In subsystem TD-DFT, interacting MNPs are treated independently, thus allowing us to compute the plasmon couplings directly from the subsystem TD-DFT transition densities. We show that subsystem TD-DFT, as well as a simplified version of it in which kinetic contributions are neglected, can reproduce the reference TD-DFT calculations for gap distances greater than about 6 Å or even smaller in the case of hybrid plasmonic systems (i.e., molecules interacting with MNPs). We also show that the subsystem TD-DFT can be also used as a tool to analyze the impact of charge-transfer effects.
Collapse
Affiliation(s)
- Giulia Giannone
- Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, Arnesano (LE) 73010, Italy.,Department of Mathematics and Physics "E. De Giorgi", University of Salento, Via Arnesano, Lecce 73100, Italy
| | - Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudzia̧dzka 5, Toruń 87-100, Poland
| | - Stefania D'Agostino
- Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, Arnesano (LE) 73010, Italy.,Department of Mathematics and Physics "E. De Giorgi", University of Salento, Via Arnesano, Lecce 73100, Italy.,Institute of Nanotechnology, National Research Council (CNR-NANOTEC), c/o Campus Ecotekne, via Monteroni, Lecce 73100, Italy
| | - Eduardo Fabiano
- Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, Arnesano (LE) 73010, Italy.,Institute for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, Lecce 73100, Italy
| | - Fabio Della Sala
- Center for Biomolecular Nanotechnologies, Istituto Italiano di Tecnologia, Via Barsanti 14, Arnesano (LE) 73010, Italy.,Institute for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, Lecce 73100, Italy
| |
Collapse
|
23
|
Wang Z, Liu W. iOI: An Iterative Orbital Interaction Approach for Solving the Self-Consistent Field Problem. J Chem Theory Comput 2021; 17:4831-4845. [PMID: 34240856 DOI: 10.1021/acs.jctc.1c00445] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
An iterative orbital interaction (iOI) approach is proposed to solve, in a bottom-up fashion, the self-consistent field problem in quantum chemistry. While it belongs grossly to the family of fragment-based quantum chemical methods, iOI is distinctive in that (1) it divides and conquers not only the energy but also the wave function and that (2) the subsystem sizes are automatically determined by successively merging neighboring small subsystems until they are just enough for converging the wave function to a given accuracy. Orthonormal occupied and virtual localized molecular orbitals are obtained in a natural manner, which can be used for all post-SCF purposes.
Collapse
Affiliation(s)
- Zikuan Wang
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, P. R. China
| | - Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, Shandong 266237, P. R. China
| |
Collapse
|
24
|
Schmitt-Monreal D, Jacob CR. Density-Based Many-Body Expansion as an Efficient and Accurate Quantum-Chemical Fragmentation Method: Application to Water Clusters. J Chem Theory Comput 2021; 17:4144-4156. [PMID: 34196558 DOI: 10.1021/acs.jctc.1c00340] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Fragmentation methods based on the many-body expansion offer an attractive approach for the quantum-chemical treatment of large molecular systems, such as molecular clusters and crystals. Conventionally, the many-body expansion is performed for the total energy, but such an energy-based many-body expansion often suffers from a slow convergence with respect to the expansion order. For systems that show strong polarization effects such as water clusters, this can render the energy-based many-body expansion infeasible. Here, we establish a density-based many-body expansion as a promising alternative approach. By performing the many-body expansion for the electron density instead of the total energy and inserting the resulting total electron density into the total energy functional of density functional theory, one can derive a density-based energy correction, which in principle accounts for all higher-order polarization effects. Here, we systematically assess the accuracy of such a density-based many-body expansion for test sets of water clusters. We show that already a density-based two-body expansion is able to reproduce interaction energies per fragment within chemical accuracy and is able to accurately predict the energetic ordering as well as the relative interaction energies of different isomers of water clusters.
Collapse
Affiliation(s)
- Daniel Schmitt-Monreal
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstr. 17, 38106 Braunschweig, Germany
| | - Christoph R Jacob
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstr. 17, 38106 Braunschweig, Germany
| |
Collapse
|
25
|
Sen R, González-Espinoza CE, Zech A, Dreuw A, Wesolowski TA. Benchmark of the Extension of Frozen-Density Embedding Theory to Nonvariational Correlated Methods: The Embedded-MP2 Case. J Chem Theory Comput 2021; 17:4049-4062. [PMID: 34137597 DOI: 10.1021/acs.jctc.1c00228] [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/30/2022]
Abstract
The extension of the frozen-density embedding theory for nonvariational methods [J. Chem. Theory Comput. 2020, 16, 6880] was utilized to evaluate intermolecular interaction energies for complexes in the Zhao-Truhlar basis set. In the applied method (FDET-MP2-FAT-LDA), the same auxiliary system is used to evaluate the correlation energy by means of the second-order Møller-Plesset perturbation theory (MP2), as in our previous work [J. Chem. Phys. 2019, 150, 121101]. Local density approximation is used for ExcTnad[ρA,ρB] in both cases. Additionally, the contribution to the energy due to the neglected correlation potential was evaluated and analyzed. The domain of applicability of the local density approximation for ExcTnad[ρA,ρB] was determined based on deviations from the interaction energies from the conventional MP2 calculations. The local density approximation for ExcTnad[ρA,ρB] performs well for hydrogen- or dipole-bound complexes. The relative errors in the interaction energy lie within 3-30%. While for charge-transfer complexes, this approximation fails consistently, and for other types of complexes, the performance of this approximation is not systematic. The sources of error are discussed in detail.
Collapse
Affiliation(s)
- Reena Sen
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University Heidelberg, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | | | - Alexander Zech
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University Heidelberg, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Tomasz A Wesolowski
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| |
Collapse
|
26
|
Shao X, Mi W, Pavanello M. GGA-Level Subsystem DFT Achieves Sub-kcal/mol Accuracy Intermolecular Interactions by Mimicking Nonlocal Functionals. J Chem Theory Comput 2021; 17:3455-3461. [PMID: 33983729 DOI: 10.1021/acs.jctc.1c00283] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
The key feature of nonlocal kinetic energy functionals is their ability to reduce to the Thomas-Fermi functional in the regions of high density and to the von Weizsäcker functional in the region of low-density/high reduced density gradient. This behavior is crucial when these functionals are employed in subsystem DFT simulations to approximate the nonadditive kinetic energy. We propose a GGA nonadditive kinetic energy functional which mimics the good behavior of nonlocal functionals, retaining the computational complexity of typical semilocal functionals. Crucially, this functional depends on the inter-subsystem density overlap. The new functional reproduces Kohn-Sham DFT and benchmark CCSD(T) interaction energies of weakly interacting dimers in the S22-5 and S66 test sets with a mean absolute deviation well below 1 kcal/mol.
Collapse
Affiliation(s)
- Xuecheng Shao
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| | - Wenhui Mi
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States
| | - Michele Pavanello
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, United States.,Department of Physics, Rutgers University, Newark, New Jersey 07102, United States
| |
Collapse
|
27
|
Brunken C, Reiher M. Automated Construction of Quantum–Classical Hybrid Models. J Chem Theory Comput 2021; 17:3797-3813. [DOI: 10.1021/acs.jctc.1c00178] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Affiliation(s)
- Christoph Brunken
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Markus Reiher
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
28
|
Egidi F, Angelico S, Lafiosca P, Giovannini T, Cappelli C. A polarizable three-layer frozen density embedding/molecular mechanics approach. J Chem Phys 2021; 154:164107. [PMID: 33940798 DOI: 10.1063/5.0045574] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
We present a novel multilayer polarizable embedding approach in which the system is divided into three portions, two of which are treated using density functional theory and their interaction is based on frozen density embedding (FDE) theory, and both also mutually interact with a polarizable classical layer described using an atomistic model based on fluctuating charges (FQ). The efficacy of the model is demonstrated by extending the formalism to linear response properties and applying it to the simulation of the excitation energies of organic molecules in aqueous solution, where the solute and the first solvation shell are treated using FDE, while the rest of the solvent is modeled using FQ charges.
Collapse
Affiliation(s)
- Franco Egidi
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Sara Angelico
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Piero Lafiosca
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Tommaso Giovannini
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
29
|
Chi YC, Shaban Tameh M, Huang C. Efficient Embedded Cluster Density Approximation Calculations with an Orbital-Free Treatment of Environments. J Chem Theory Comput 2021; 17:2737-2751. [PMID: 33856795 DOI: 10.1021/acs.jctc.0c01133] [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/28/2022]
Abstract
The computational cost of the Kohn-Sham density functional theory (KS-DFT), employing advanced orbital-based exchange-correlation (XC) functionals, increases quickly for large systems. To tackle this problem, we recently developed a local correlation method in the framework of KS-DFT: the embedded cluster density approximation (ECDA). The aim of ECDA is to obtain accurate electronic structures in an entire system. With ECDA, for each atom in a system, we define a cluster to enclose that atom, with the rest atoms treated as the environment. The system's electron density is then partitioned among the cluster and the environment. The cluster's XC energy density is then calculated based on its electron density using an advanced orbital-based XC functional. The system's XC energy is obtained by patching all clusters' XC energy densities in an atom-by-atom manner. In our previous formulation of ECDA, environments were treated by KS-DFT, which makes the following two tasks computationally expensive for large systems. The first task is to partition the system's electron density among a cluster and its environment. The second task is to solve the environments' Sternheimer equations for calculating the system's XC potential. In this work, we remove these two computational bottlenecks by treating the environments with the orbital-free (OF) DFT. The new method is called ECDA-envOF. The performance of ECDA-envOF is examined in two systems: ester and Cl-tetracene, for which the exact exchange (EXX) is used as the advanced XC functional. We show that ECDA-envOF gives results that are very close to the previous formulation in which the environments were treated by KS-DFT. Therefore, ECDA-envOF can be used for future large-scale simulations. Another focus of this work is to examine ECDA-envOF's performance on systems having different bond types. With ECDA-envOF, we calculate the energy curves for stretching/compressing some covalent, metallic, and ionic systems. ECDA-envOF's predictions agree well with the benchmarks by using reasonably large clusters. These examples demonstrate that ECDA-envOF is nearly a black-box local correlation method for investigating heterogeneous materials in which different bond types exist.
Collapse
Affiliation(s)
- Yu-Chieh Chi
- Department of Scientific Computing, Florida State University, Tallahassee, Florida 32306, United States
| | - Maliheh Shaban Tameh
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Chen Huang
- Department of Scientific Computing, Materials Science and Engineering Program, National High Magnetic Field Laboratory, Florida State University, Tallahassee, Florida 32306, United States
| |
Collapse
|
30
|
Marrazzini G, Giovannini T, Scavino M, Egidi F, Cappelli C, Koch H. Multilevel Density Functional Theory. J Chem Theory Comput 2021; 17:791-803. [PMID: 33449681 PMCID: PMC7880574 DOI: 10.1021/acs.jctc.0c00940] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
![]()
Following recent
developments in multilevel embedding methods,
we introduce a novel density matrix-based multilevel approach within
the framework of density functional theory (DFT). In this multilevel
DFT, the system is partitioned in an active and an inactive fragment,
and all interactions are retained between the two parts. The decomposition
of the total system is performed upon the density matrix. The orthogonality
between the two parts is maintained by solving the Kohn–Sham
equations in the MO basis for the active part only, while keeping
the inactive density matrix frozen. This results in the reduction
of computational cost. We outline the theory and implementation and
discuss the differences and similarities with state-of-the-art DFT
embedding methods. We present applications to aqueous solutions of
methyloxirane and glycidol.
Collapse
Affiliation(s)
- Gioia Marrazzini
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Tommaso Giovannini
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | - Marco Scavino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Franco Egidi
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Chiara Cappelli
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Henrik Koch
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
31
|
Powell AD, Kroes GJ, Doblhoff-Dier K. Quantum Monte Carlo calculations on dissociative chemisorption of H2 + Al(110): Minimum barrier heights and their comparison to DFT values. J Chem Phys 2020; 153:224701. [DOI: 10.1063/5.0022919] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Andrew D. Powell
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, Netherlands
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, Netherlands
| | - Katharina Doblhoff-Dier
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, Netherlands
| |
Collapse
|
32
|
Massolle A, Neugebauer J. Subsystem density-functional theory for interacting open-shell systems: spin densities and magnetic exchange couplings. Faraday Discuss 2020; 224:201-226. [PMID: 33000819 DOI: 10.1039/d0fd00063a] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
We investigate the possibility of describing interacting open-shell systems in high-spin and broken-symmetry (BS) states with subsystem density-functional theory (sDFT). This subsystem method typically starts from the electronic-structure results obtained for individual systems, for which the spin states can be individually defined. Through the confining effect of the embedding potential and/or the use of monomer basis sets, these individual spin states can be preserved in sDFT calculations. This offers the possibility of easy convergence to broken-symmetry states with arbitrary local spin patterns. We show that the resulting spin densities are in very good agreement with successfully converged broken-symmetry Kohn-Sham density-functional theory (KS-DFT) calculations. Yet sDFT can even cure those BS cases where KS-DFT suffers from convergence problems or convergence to undesired spin states. In contrast to KS-DFT, the sDFT-results only show a mild exchange-correlation functional dependence. We also show that magnetic coupling constants from sDFT are not satisfactory with standard approximations for the non-additive kinetic energy. When this component is evaluated "exactly", i.e. based on potential reconstruction, however, the magnetic coupling constants derived from spin-state energy differences are greatly improved. Hence, the interacting radicals studied here represent cases where even (semi-)local approximations for the non-additive kinetic-energy potential work well, while the parent energy functionals do not yield satisfactory results for spin-state energy differences.
Collapse
Affiliation(s)
- Anja Massolle
- Theoretische Organische Chemie, Organisch-Chemisches Institut, Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany.
| | | |
Collapse
|
33
|
De Santis M, Belpassi L, Jacob CR, Severo Pereira Gomes A, Tarantelli F, Visscher L, Storchi L. Environmental Effects with Frozen-Density Embedding in Real-Time Time-Dependent Density Functional Theory Using Localized Basis Functions. J Chem Theory Comput 2020; 16:5695-5711. [PMID: 32786918 PMCID: PMC8009524 DOI: 10.1021/acs.jctc.0c00603] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2020] [Indexed: 12/14/2022]
Abstract
Frozen-density embedding (FDE) represents a versatile embedding scheme to describe the environmental effect on electron dynamics in molecular systems. The extension of the general theory of FDE to the real-time time-dependent Kohn-Sham method has previously been presented and implemented in plane waves and periodic boundary conditions [Pavanello, M.; J. Chem. Phys. 2015, 142, 154116]. In the current paper, we extend our recent formulation of the real-time time-dependent Kohn-Sham method based on localized basis set functions and developed within the Psi4NumPy framework to the FDE scheme. The latter has been implemented in its "uncoupled" flavor (in which the time evolution is only carried out for the active subsystem, while the environment subsystems remain at their ground state), using and adapting the FDE implementation already available in the PyEmbed module of the scripting framework PyADF. The implementation was facilitated by the fact that both Psi4NumPy and PyADF, being native Python API, provided an ideal framework of development using the Python advantages in terms of code readability and reusability. We employed this new implementation to investigate the stability of the time-propagation procedure, which is based on an efficient predictor/corrector second-order midpoint Magnus propagator employing an exact diagonalization, in combination with the FDE scheme. We demonstrate that the inclusion of the FDE potential does not introduce any numerical instability in time propagation of the density matrix of the active subsystem, and in the limit of the weak external field, the numerical results for low-lying transition energies are consistent with those obtained using the reference FDE calculations based on the linear-response TDDFT. The method is found to give stable numerical results also in the presence of a strong external field inducing nonlinear effects. Preliminary results are reported for high harmonic generation (HHG) of a water molecule embedded in a small water cluster. The effect of the embedding potential is evident in the HHG spectrum reducing the number of the well-resolved high harmonics at high energy with respect to the free water. This is consistent with a shift toward lower ionization energy passing from an isolated water molecule to a small water cluster. The computational burden for the propagation step increases approximately linearly with the size of the surrounding frozen environment. Furthermore, we have also shown that the updating frequency of the embedding potential may be significantly reduced, much less than one per time step, without jeopardizing the accuracy of the transition energies.
Collapse
Affiliation(s)
- Matteo De Santis
- Dipartimento di
Chimica, Biologia e Biotecnologie, Università
degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
- Istituto di Scienze
e Tecnologie Chimiche (SCITEC), Consiglio Nazionale delle Ricerche
c/o Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
| | - Leonardo Belpassi
- Istituto di Scienze
e Tecnologie Chimiche (SCITEC), Consiglio Nazionale delle Ricerche
c/o Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
| | - Christoph R. Jacob
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstr. 17, 38106 Braunschweig, Germany
| | | | - Francesco Tarantelli
- Dipartimento di
Chimica, Biologia e Biotecnologie, Università
degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
| | - Lucas Visscher
- Theoretical Chemistry, Faculty of Science, Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - Loriano Storchi
- Istituto di Scienze
e Tecnologie Chimiche (SCITEC), Consiglio Nazionale delle Ricerche
c/o Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
- Dipartimento di Farmacia, Università
degli Studi ‘G. D’Annunzio’, Via dei Vestini 31, 66100 Chieti, Italy
| |
Collapse
|
34
|
Hoyer CE, Li X. Relativistic two-component projection-based quantum embedding for open-shell systems. J Chem Phys 2020; 153:094113. [DOI: 10.1063/5.0012433] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Chad E. Hoyer
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| |
Collapse
|
35
|
Ricardi N, Ernst M, Macchi P, Wesolowski TA. Embedding-theory-based simulations using experimental electron densities for the environment. Acta Crystallogr A Found Adv 2020; 76:571-579. [PMID: 32869754 PMCID: PMC7459768 DOI: 10.1107/s2053273320008062] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2020] [Accepted: 06/16/2020] [Indexed: 11/26/2022] Open
Abstract
The basic idea of frozen-density embedding theory (FDET) is the constrained minimization of the Hohenberg-Kohn density functional EHK[ρ] performed using the auxiliary functional E_{v_{AB}}^{\rm FDET}[\Psi _A, \rho _B], where ΨA is the embedded NA-electron wavefunction and ρB(r) is a non-negative function in real space integrating to a given number of electrons NB. This choice of independent variables in the total energy functional E_{v_{AB}}^{\rm FDET}[\Psi _A, \rho _B] makes it possible to treat the corresponding two components of the total density using different methods in multi-level simulations. The application of FDET using ρB(r) reconstructed from X-ray diffraction data for a molecular crystal is demonstrated for the first time. For eight hydrogen-bonded clusters involving a chromophore (represented as ΨA) and the glycylglycine molecule [represented as ρB(r)], FDET is used to derive excitation energies. It is shown that experimental densities are suitable for use as ρB(r) in FDET-based simulations.
Collapse
Affiliation(s)
- Niccolò Ricardi
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| | - Michelle Ernst
- University of Bern, Freiestraße 3, 3012 Bern, Switzerland
| | - Piero Macchi
- Department of Chemistry, Materials and Chemical Engineering, Polytechnic of Milan, via Mancinelli 7, Milano 20131, Italy
| | - Tomasz Adam Wesolowski
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| |
Collapse
|
36
|
Bensberg M, Neugebauer J. Orbital Alignment for Accurate Projection-Based Embedding Calculations along Reaction Paths. J Chem Theory Comput 2020; 16:3607-3619. [DOI: 10.1021/acs.jctc.0c00104] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Moritz Bensberg
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| |
Collapse
|
37
|
Bensberg M, Neugebauer J. Density functional theory based embedding approaches for transition-metal complexes. Phys Chem Chem Phys 2020; 22:26093-26103. [PMID: 33201953 DOI: 10.1039/d0cp05188h] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
Transition metal species are commonly discussed by considering the metal atom embedded in a ligand environment. This apparently makes them interesting targets for modern embedding strategies based on Kohn-Sham density functional theory (DFT), which aim at modelling accurate predictions for large systems by combining different quantum chemical methods. In this perspective, we will focus on subsystem density functional theory and projection-based embedding. We review the developments in the field for transition metal species, demonstrate benefits, drawbacks and analyse error sources of the different strategies using the example of chromium hexacarbonyle, before giving a perspective where the field is currently heading.
Collapse
Affiliation(s)
- Moritz Bensberg
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster Corrensstraße 36, 48149 Münster, Germany.
| | | |
Collapse
|
38
|
Wen X, Graham DS, Chulhai DV, Goodpaster JD. Absolutely Localized Projection-Based Embedding for Excited States. J Chem Theory Comput 2019; 16:385-398. [DOI: 10.1021/acs.jctc.9b00959] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Xuelan Wen
- Department of Chemistry, University of Minnesota, 207 Pleasant Street Southeast, Minneapolis, Minnesota 55455, United States
| | - Daniel S. Graham
- Department of Chemistry, University of Minnesota, 207 Pleasant Street Southeast, Minneapolis, Minnesota 55455, United States
| | - Dhabih V. Chulhai
- Department of Chemistry, University of Minnesota, 207 Pleasant Street Southeast, Minneapolis, Minnesota 55455, United States
| | - Jason D. Goodpaster
- Department of Chemistry, University of Minnesota, 207 Pleasant Street Southeast, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
39
|
Li W, Chen M, Rabani E, Baer R, Neuhauser D. Stochastic embedding DFT: Theory and application to p-nitroaniline in water. J Chem Phys 2019; 151:174115. [PMID: 31703523 DOI: 10.1063/1.5110226] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Over this past decade, we combined the idea of stochastic resolution of identity with a variety of electronic structure methods. In our stochastic Kohn-Sham density functional theory (DFT) method, the density is an average over multiple stochastic samples, with stochastic errors that decrease as the inverse square root of the number of sampling orbitals. Here, we develop a stochastic embedding density functional theory method (se-DFT) that selectively reduces the stochastic error (specifically on the forces) for a selected subsystem(s). The motivation, similar to that of other quantum embedding methods, is that for many systems of practical interest, the properties are often determined by only a small subsystem. In stochastic embedding DFT, two sets of orbitals are used: a deterministic one associated with the embedded subspace and the rest, which is described by a stochastic set. The method agrees exactly with deterministic calculations in the limit of a large number of stochastic samples. We apply se-DFT to study a p-nitroaniline molecule in water, where the statistical errors in the forces on the system (the p-nitroaniline molecule) are reduced by an order of magnitude compared with nonembedding stochastic DFT.
Collapse
Affiliation(s)
- Wenfei Li
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, USA
| | - Ming Chen
- Department of Chemistry, University of California and Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Eran Rabani
- Department of Chemistry, University of California and Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Roi Baer
- Fritz Haber Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Daniel Neuhauser
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, USA
| |
Collapse
|
40
|
Schulz A, Jacob CR. Description of intermolecular charge transfer with subsystem density-functional theory. J Chem Phys 2019; 151:131103. [PMID: 31594348 DOI: 10.1063/1.5125218] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Efficient quantum-chemical methods that are able to describe intermolecular charge transfer are crucial for modeling organic semiconductors. However, the correct description of intermolecular charge transfer with density-functional theory (DFT) is hampered by the fractional charge error of approximate exchange-correlation (xc) functionals. Here, we investigate the charge transfer induced by an external electric field in a tetrathiafulvalene-tetracyanoquinodimethane (TTF-TCNQ) complex as a test case. For this seemingly simple model system, a supermolecular DFT treatment fails with most conventional xc functionals. Here, we present an extension of subsystem DFT to subsystems with a fractional number of electrons. We show that within such a framework, it becomes possible to overcome the fractional charge error by enforcing the correct dependence of each subsystem's total energy on the subsystem's fractional charge. Such a subsystem DFT approach allows for a correct description of the intermolecular charge transfer in the TTF-TCNQ model complex. The approach presented here can be generalized to larger molecular aggregates and will thus allow for modeling organic semiconductor materials accurately and efficiently.
Collapse
Affiliation(s)
- Anika Schulz
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry, Gaußstr. 17, 38106 Braunschweig, Germany
| | - Christoph R Jacob
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry, Gaußstr. 17, 38106 Braunschweig, Germany
| |
Collapse
|
41
|
Huang C. Analytical energy gradient for the embedded cluster density approximation. J Chem Phys 2019; 151:134101. [DOI: 10.1063/1.5112789] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Chen Huang
- Department of Scientific Computing, Materials Science and Engineering Program, and National High Magnetic Field Laboratory, Florida State University, Tallahassee, Florida 32306, USA
| |
Collapse
|
42
|
Bensberg M, Neugebauer J. Direct orbital selection for projection-based embedding. J Chem Phys 2019; 150:214106. [DOI: 10.1063/1.5099007] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Affiliation(s)
- Moritz Bensberg
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| |
Collapse
|
43
|
Abstract
Complex chemical systems present challenges to electronic structure theory stemming from large system sizes, subtle interactions, coupled dynamical time scales, and electronically nonadiabatic effects. New methods are needed to perform reliable, rigorous, and affordable electronic structure calculations for simulating the properties and dynamics of such systems. This Account reviews projection-based quantum embedding for electronic structure, which provides a formally exact method for density functional theory (DFT) embedding. The method also provides a rigorous and accurate approach for describing a small part of a chemical system at the level of a correlated wavefunction (WF) method while the remainder of the system is described at the level of DFT. A key advantage of projection-based embedding is that it can be formulated in terms of an extremely simple level-shift projection operator, which eliminates the need for any optimized effective potential calculation or kinetic energy functional approximation while simultaneously ensuring that no extra programming is needed to perform WF-in-DFT embedding with an arbitrary WF method. The current work presents the theoretical underpinnings of projection-based embedding, describes use of the method for combining wavefunction and density functional theories, and discusses technical refinements that have improved the applicability and robustness of the method. Applications of projection-based WF-in-DFT embedding are also reviewed, with particular focus on recent work on transition-metal catalysis, enzyme reactivity, and battery electrolyte decomposition. In particular, we review the application of projection-based embedding for the prediction of electrochemical potentials and reaction pathways in a Co-centered hydrogen evolution catalyst. Projection-based WF-in-DFT calculations are shown to provide quantitative accuracy while greatly reducing the computational cost compared with a reference coupled cluster calculation on the full system. Additionally, projection-based WF-in-DFT embedding is used to study the mechanism of citrate synthase; it is shown that projection-based WF-in-DFT largely eliminates the sensitivity of the potential energy landscape to the employed DFT exchange-correlation functional. Finally, we demonstrate the use of projection-based WF-in-DFT to study electron transfer reactions associated with battery electrolyte decomposition. Projection-based WF-in-DFT embedding is used to calculate the oxidation potentials of neat ethylene carbonate (EC), neat dimethyl carbonate (DMC), and 1:1 mixtures of EC and DMC in order to overcome qualitative inaccuracies in the electron densities and ionization energies obtained from conventional DFT methods. By further embedding the WF-in-DFT description in a molecular mechanics point-charge environment, this work enables an explicit description of the solvent and ensemble averaging of the solvent configurations. Looking forward, we anticipate continued refinement of the projection-based embedding methodology as well as its increasingly widespread application in diverse areas of chemistry, biology, and materials science.
Collapse
Affiliation(s)
- Sebastian J. R. Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Matthew Welborn
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Frederick R. Manby
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Thomas F. Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
44
|
Hoyer CE, Williams-Young DB, Huang C, Li X. Embedding non-collinear two-component electronic structure in a collinear quantum environment. J Chem Phys 2019; 150:174114. [DOI: 10.1063/1.5092628] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Affiliation(s)
- Chad E. Hoyer
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | | | - Chen Huang
- Department of Scientific Computing, Materials Science and Engineering Program, and National High Magnetic Field Laboratory, Florida State University, Tallahassee, Florida 32306, USA
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| |
Collapse
|
45
|
Haldar S, Riplinger C, Demoulin B, Neese F, Izsak R, Dutta AK. Multilayer Approach to the IP-EOM-DLPNO-CCSD Method: Theory, Implementation, and Application. J Chem Theory Comput 2019; 15:2265-2277. [DOI: 10.1021/acs.jctc.8b01263] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Soumi Haldar
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400 076, India
| | | | - Baptiste Demoulin
- Max-Planck-Institut
für Kohlenforschung, 1 Kaiser-Wilhelm Platz, Mülheim an der Ruhr 45470, Germany
| | - Frank Neese
- Max-Planck-Institut
für Kohlenforschung, 1 Kaiser-Wilhelm Platz, Mülheim an der Ruhr 45470, Germany
| | - Robert Izsak
- Max-Planck-Institut
für Kohlenforschung, 1 Kaiser-Wilhelm Platz, Mülheim an der Ruhr 45470, Germany
| | - Achintya Kumar Dutta
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400 076, India
| |
Collapse
|
46
|
Zhang X, Carter EA. Subspace Density Matrix Functional Embedding Theory: Theory, Implementation, and Applications to Molecular Systems. J Chem Theory Comput 2018; 15:949-960. [DOI: 10.1021/acs.jctc.8b00990] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Xing Zhang
- Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, New Jersey 08544-5263, United States
| | - Emily A. Carter
- School of Engineering and Applied Science, Princeton University, Princeton, New Jersey 08544-5263, United States
| |
Collapse
|
47
|
Mühlbach AH, Reiher M. Quantum system partitioning at the single-particle level. J Chem Phys 2018; 149:184104. [DOI: 10.1063/1.5055942] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Adrian H. Mühlbach
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, CH-8093 Zürich, Switzerland
| | - Markus Reiher
- Laboratorium für Physikalische Chemie, ETH Zürich, Vladimir-Prelog-Weg 2, CH-8093 Zürich, Switzerland
| |
Collapse
|
48
|
Huang C. Embedded Cluster Density Approximation for Exchange–Correlation Energy: A Natural Extension of the Local Density Approximation. J Chem Theory Comput 2018; 14:6211-6225. [DOI: 10.1021/acs.jctc.8b00471] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Chen Huang
- Department of Scientific Computing and National High Magnetic Field Laboratory, Florida State University, Tallahassee, Florida 32306, United States
| |
Collapse
|
49
|
Jiang K, Nafziger J, Wasserman A. Constructing a non-additive non-interacting kinetic energy functional approximation for covalent bonds from exact conditions. J Chem Phys 2018; 149:164112. [PMID: 30384751 DOI: 10.1063/1.5051455] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Kaili Jiang
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Ave., West Lafayette, Indiana 47907, USA
| | - Jonathan Nafziger
- Department of Chemistry, Purdue University, 560 Oval Dr., West Lafayette Indiana 47907, USA
| | - Adam Wasserman
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Ave., West Lafayette, Indiana 47907, USA
- Department of Chemistry, Purdue University, 560 Oval Dr., West Lafayette Indiana 47907, USA
| |
Collapse
|
50
|
Welborn M, Manby FR, Miller TF. Even-handed subsystem selection in projection-based embedding. J Chem Phys 2018; 149:144101. [DOI: 10.1063/1.5050533] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Affiliation(s)
- Matthew Welborn
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Frederick R. Manby
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Thomas F. Miller
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|