1
|
Giovannini T, Marrazzini G, Scavino M, Koch H, Cappelli C. Integrated Multiscale Multilevel Approach to Open Shell Molecular Systems. J Chem Theory Comput 2023; 19:1446-1456. [PMID: 36780359 PMCID: PMC10018740 DOI: 10.1021/acs.jctc.2c00805] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/14/2023]
Abstract
We present a novel multiscale approach to study the electronic structure of open shell molecular systems embedded in an external environment. The method is based on the coupling of multilevel Hartree-Fock (MLHF) and Density Functional Theory (MLDFT), suitably extended to the unrestricted formalism, to Molecular Mechanics (MM) force fields (FF). Within the ML region, the system is divided into active and inactive parts, thus describing the most relevant interactions (electrostatic, polarization, and Pauli repulsion) at the quantum level. The surrounding MM part, which is formulated in terms of nonpolarizable or polarizable FFs, permits a physically consistent treatment of long-range electrostatics and polarization effects. The approach is extended to the calculation of hyperfine coupling constants and applied to selected nitroxyl radicals in an aqueous solution.
Collapse
Affiliation(s)
| | - Gioia Marrazzini
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Marco Scavino
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Henrik Koch
- Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy.,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
|
2
|
Sen S, Senjean B, Visscher L. Characterization of excited states in time-dependent density functional theory using localized molecular orbitals. J Chem Phys 2023; 158:054115. [PMID: 36754801 DOI: 10.1063/5.0137729] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
Abstract
Localized molecular orbitals are often used for the analysis of chemical bonds, but they can also serve to efficiently and comprehensibly compute linear response properties. While conventional canonical molecular orbitals provide an adequate basis for the treatment of excited states, a chemically meaningful identification of the different excited-state processes is difficult within such a delocalized orbital basis. In this work, starting from an initial set of supermolecular canonical molecular orbitals, we provide a simple one-step top-down embedding procedure for generating a set of orbitals, which are localized in terms of the supermolecule but delocalized over each subsystem composing the supermolecule. Using an orbital partitioning scheme based on such sets of localized orbitals, we further present a procedure for the construction of local excitations and charge-transfer states within the linear response framework of time-dependent density functional theory (TDDFT). This procedure provides direct access to approximate diabatic excitation energies and, under the Tamm-Dancoff approximation, also their corresponding electronic couplings-quantities that are of primary importance in modeling energy transfer processes in complex biological systems. Our approach is compared with a recently developed diabatization procedure based on subsystem TDDFT using projection operators, which leads to a similar set of working equations. Although both of these methods differ in the general localization strategies adopted and the type of basis functions (Slaters vs Gaussians) employed, an overall decent agreement is obtained.
Collapse
Affiliation(s)
- Souloke Sen
- Division of Theoretical Chemistry, Faculty of Sciences, Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - Bruno Senjean
- ICGM, Université de Montpellier, CNRS, ENSCM, Montpellier, France
| | - Lucas Visscher
- Division of Theoretical Chemistry, Faculty of Sciences, Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| |
Collapse
|
3
|
Chen M, Baer R, Rabani E. Structure optimization with stochastic density functional theory. J Chem Phys 2023; 158:024111. [PMID: 36641385 DOI: 10.1063/5.0126169] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
Linear-scaling techniques for Kohn-Sham density functional theory are essential to describe the ground state properties of extended systems. Still, these techniques often rely on the localization of the density matrix or accurate embedding approaches, limiting their applicability. In contrast, stochastic density functional theory (sDFT) achieves linear- and sub-linear scaling by statistically sampling the ground state density without relying on embedding or imposing localization. In return, ground state observables, such as the forces on the nuclei, fluctuate in sDFT, making optimizing the nuclear structure a highly non-trivial problem. In this work, we combine the most recent noise-reduction schemes for sDFT with stochastic optimization algorithms to perform structure optimization within sDFT. We compare the performance of the stochastic gradient descent approach and its variations (stochastic gradient descent with momentum) with stochastic optimization techniques that rely on the Hessian, such as the stochastic Broyden-Fletcher-Goldfarb-Shanno algorithm. We further provide a detailed assessment of the computational efficiency and its dependence on the optimization parameters of each method for determining the ground state structure of bulk silicon with varying supercell dimensions.
Collapse
Affiliation(s)
- Ming Chen
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | - Roi Baer
- Fritz Haber Center of Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Eran Rabani
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| |
Collapse
|
4
|
Initial Maximum Overlap Method Embedded with Extremely Localized Molecular Orbitals for Core-Ionized States of Large Systems. MOLECULES (BASEL, SWITZERLAND) 2022; 28:molecules28010136. [PMID: 36615331 PMCID: PMC9822432 DOI: 10.3390/molecules28010136] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/29/2022] [Revised: 12/21/2022] [Accepted: 12/22/2022] [Indexed: 12/28/2022]
Abstract
Despite great advances in X-ray absorption spectroscopy for the investigation of small molecule electronic structure, the application to biosystems of experimental techniques developed within this research field remains a challenge. To partially circumvent the problem, users resort to theoretical methods to interpret or predict the X-ray absorption spectra of large molecules. To accomplish this task, only low-cost computational strategies can be exploited. For this reason, some of them are single Slater determinant wavefunction approaches coupled with multiscale embedding techniques designed to treat large systems of biological interest. Therefore, in this work, we propose to apply the recently developed IMOM/ELMO embedding method to the determination of core-ionized states. The IMOM/ELMO technique resulted from the combination of the single Slater determinant Δself-consistent-field-initial maximum overlap approach (ΔSCF-IMOM) with the QM/ELMO (quantum mechanics/extremely localized molecular orbital) embedding strategy, a method where only the chemically relevant region of the examined system is treated at fully quantum chemical level, while the rest is described through transferred and frozen extremely localized molecular orbitals (ELMOs). The IMOM/ELMO technique was initially validated by computing core-ionization energies for small molecules, and it was afterwards exploited to study larger biosystems. The obtained results are in line with those reported in previous studies that applied alternative ΔSCF approaches. This makes us envisage a possible future application of the proposed method to the interpretation of X-ray absorption spectra of large molecules.
Collapse
|
5
|
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
|
6
|
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] [Key Words] [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@( H 2 O ) 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
|
7
|
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
|
8
|
Macetti G, Genoni A. Initial Maximum Overlap Method for Large Systems by the Quantum Mechanics/Extremely Localized Molecular Orbital Embedding Technique. J Chem Theory Comput 2021; 17:4169-4182. [PMID: 34196174 DOI: 10.1021/acs.jctc.1c00388] [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/30/2022]
Abstract
Quantum chemistry offers a large variety of methods to treat excited states. Many of them are based on a multireference wave function ansatz and are therefore characterized by an intrinsic complexity and high computational costs. To overcome these drawbacks and also some limitations of simpler single-reference approaches (e.g., configuration interaction singles and time-dependent density functional theory), the single-determinant Δself-consistent field-initial maximum overlap method (ΔSCF-IMOM) has been proposed. This strategy substitutes the aufbau principle with a criterion that occupies molecular orbitals at successive SCF iterations on the basis of their maximum overlap with a proper set of guess orbitals for the target excited state. In this way, it prevents the SCF process to collapse to the ground state wave function and provides excited state single Slater determinant solutions to the SCF equations. Here, we propose to extend the applicability of the IMOM to the treatment of localized excited states of large systems. To accomplish this task, we coupled it with the QM/ELMO (quantum mechanics/extremely localized molecular orbitals) strategy, a quantum mechanical embedding method in which the most chemically relevant part of the system is treated with traditional quantum chemical approaches, while the rest is described by extremely localized molecular orbitals transferred from recently constructed libraries or proper model molecules. After presenting the theoretical foundations of the new IMOM/ELMO technique, in this paper, we will show and discuss the results of preliminary test calculations carried out on both model systems (i.e., decanoic acid, decene, decapentaene, and solvated acrolein) and a system of biological interest (flavin mononucleotide in the flavodoxin protein). We observed that, for localized excited states, the new IMOM/ELMO method provides reliable results, and it reproduces the outcomes of fully IMOM calculations within the chemical accuracy threshold (i.e., 0.043 eV) by including only a limited number of atoms in the QM region. Furthermore, the first application of our embedding technique to a larger biological system gave completely plausible results in line with those obtained through more traditional quantum mechanical methods, thus opening the possibility of using the new approach in future investigations of photobiology problems.
Collapse
Affiliation(s)
- Giovanni Macetti
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| | - Alessandro Genoni
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| |
Collapse
|
9
|
Chen M, Baer R, Neuhauser D, Rabani E. Stochastic density functional theory: Real- and energy-space fragmentation for noise reduction. J Chem Phys 2021; 154:204108. [PMID: 34241170 DOI: 10.1063/5.0044163] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Stochastic density functional theory (sDFT) is becoming a valuable tool for studying ground-state properties of extended materials. The computational complexity of describing the Kohn-Sham orbitals is replaced by introducing a set of random (stochastic) orbitals leading to linear and often sub-linear scaling of certain ground-state observables at the account of introducing a statistical error. Schemes to reduce the noise are essential, for example, for determining the structure using the forces obtained from sDFT. Recently, we have introduced two embedding schemes to mitigate the statistical fluctuations in the electron density and resultant forces on the nuclei. Both techniques were based on fragmenting the system either in real space or slicing the occupied space into energy windows, allowing for a significant reduction in the statistical fluctuations. For chemical accuracy, further reduction of the noise is required, which could be achieved by increasing the number of stochastic orbitals. However, the convergence is relatively slow as the statistical error scales as 1/Nχ according to the central limit theorem, where Nχ is the number of random orbitals. In this paper, we combined the embedding schemes mentioned above and introduced a new approach that builds on overlapped fragments and energy windows. The new approach significantly lowers the noise for ground-state properties, such as the electron density, total energy, and forces on the nuclei, as demonstrated for a G-center in bulk silicon.
Collapse
Affiliation(s)
- Ming Chen
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Roi Baer
- Fritz Haber Center of Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Daniel Neuhauser
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095, USA
| | - Eran Rabani
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| |
Collapse
|
10
|
Macetti G, Wieduwilt EK, Genoni A. QM/ELMO: A Multi-Purpose Fully Quantum Mechanical Embedding Scheme Based on Extremely Localized Molecular Orbitals. J Phys Chem A 2021; 125:2709-2726. [DOI: 10.1021/acs.jpca.0c11450] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Affiliation(s)
- Giovanni Macetti
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| | - Erna K. Wieduwilt
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| | - Alessandro Genoni
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| |
Collapse
|
11
|
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
|
12
|
Zhang K, Ren S, Caricato M. Multistate QM/QM Extrapolation of UV/Vis Absorption Spectra with Point Charge Embedding. J Chem Theory Comput 2020; 16:4361-4372. [DOI: 10.1021/acs.jctc.0c00339] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Kaihua Zhang
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Sijin Ren
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Marco Caricato
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| |
Collapse
|
13
|
Macetti G, Wieduwilt EK, Assfeld X, Genoni A. Localized Molecular Orbital-Based Embedding Scheme for Correlated Methods. J Chem Theory Comput 2020; 16:3578-3596. [PMID: 32369363 DOI: 10.1021/acs.jctc.0c00084] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Embedding strategies currently provide the best compromise between accuracy and computational cost in modeling chemical properties and processes of large and complex systems. In this framework, different methods have been proposed all over the years, from the very popular QM/MM approaches to the more recent and very promising density matrix and density functional embedding techniques. Here, we present a further development of the quantum mechanics/extremely localized molecular orbital technique (QM/ELMO) method, a recently proposed multiscale embedding strategy in which the chemically active region of the investigated system is treated at a fully quantum mechanical level, while the rest is described by frozen extremely localized molecular orbitals previously transferred from proper libraries or tailor-made model molecules. In particular, in this work we discuss and assess in detail the extension of the QM/ELMO approach to density functional theory and post-Hartree-Fock techniques by evaluating its performances when it is used to describe chemical reactions, bond dissociations, and intermolecular interactions. The preliminary test calculations have shown that, in the investigated cases, the new embedding strategy enables the results of the corresponding fully quantum mechanical computations to be reproduced within chemical accuracy in almost all the cases but with a significantly reduced computational cost, especially when correlated post-Hartree-Fock strategies are used to describe the quantum mechanical subsystem. In light of the obtained results, we already envisage the future application of the new correlated QM/ELMO techniques to the investigation of more challenging problems, such as the modeling of enzyme catalysis, the study of excited states of biomolecules, and the refinement of macromolecular X-ray crystal structures.
Collapse
Affiliation(s)
- Giovanni Macetti
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| | - Erna K Wieduwilt
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| | - Xavier Assfeld
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, Boulevard des Aiguilletes, BP 70239, F-54506 Vandoeuvre-lès-Nancy, France
| | - Alessandro Genoni
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| |
Collapse
|
14
|
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
|
15
|
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
|
16
|
Chen M, Baer R, Neuhauser D, Rabani E. Energy window stochastic density functional theory. J Chem Phys 2019; 151:114116. [PMID: 31542024 DOI: 10.1063/1.5114984] [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
Linear scaling density functional theory is important for understanding electronic structure properties of nanometer scale systems. Recently developed stochastic density functional theory can achieve linear or even sublinear scaling for various electronic properties without relying on the sparsity of the density matrix. The basic idea relies on projecting stochastic orbitals onto the occupied space by expanding the Fermi-Dirac operator and repeating this for Nχ stochastic orbitals. Often, a large number of stochastic orbitals are required to reduce the statistical fluctuations (which scale as Nχ -1/2) below a tolerable threshold. In this work, we introduce a new stochastic density functional theory that can efficiently reduce the statistical fluctuations for certain observable and can also be integrated with an embedded fragmentation scheme. The approach is based on dividing the occupied space into energy windows and projecting the stochastic orbitals with a single expansion onto all windows simultaneously. This allows for a significant reduction of the noise as illustrated for bulk silicon with a large supercell. We also provide theoretical analysis to rationalize why the noise can be reduced only for a certain class of ground state properties, such as the forces and electron density.
Collapse
Affiliation(s)
- Ming Chen
- Department of Chemistry, University of California and Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Roi Baer
- Fritz Haber Center of Molecular Dynamics and 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
| | - Eran Rabani
- Department of Chemistry, University of California and Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
17
|
Lee SJR, Ding F, Manby FR, Miller TF. Analytical gradients for projection-based wavefunction-in-DFT embedding. J Chem Phys 2019. [DOI: 10.1063/1.5109882] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Affiliation(s)
- Sebastian J. R. Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Feizhi Ding
- 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
|
18
|
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
|
19
|
Vogiatzis KD, Polynski MV, Kirkland JK, Townsend J, Hashemi A, Liu C, Pidko EA. Computational Approach to Molecular Catalysis by 3d Transition Metals: Challenges and Opportunities. Chem Rev 2019; 119:2453-2523. [PMID: 30376310 PMCID: PMC6396130 DOI: 10.1021/acs.chemrev.8b00361] [Citation(s) in RCA: 214] [Impact Index Per Article: 42.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2018] [Indexed: 12/28/2022]
Abstract
Computational chemistry provides a versatile toolbox for studying mechanistic details of catalytic reactions and holds promise to deliver practical strategies to enable the rational in silico catalyst design. The versatile reactivity and nontrivial electronic structure effects, common for systems based on 3d transition metals, introduce additional complexity that may represent a particular challenge to the standard computational strategies. In this review, we discuss the challenges and capabilities of modern electronic structure methods for studying the reaction mechanisms promoted by 3d transition metal molecular catalysts. Particular focus will be placed on the ways of addressing the multiconfigurational problem in electronic structure calculations and the role of expert bias in the practical utilization of the available methods. The development of density functionals designed to address transition metals is also discussed. Special emphasis is placed on the methods that account for solvation effects and the multicomponent nature of practical catalytic systems. This is followed by an overview of recent computational studies addressing the mechanistic complexity of catalytic processes by molecular catalysts based on 3d metals. Cases that involve noninnocent ligands, multicomponent reaction systems, metal-ligand and metal-metal cooperativity, as well as modeling complex catalytic systems such as metal-organic frameworks are presented. Conventionally, computational studies on catalytic mechanisms are heavily dependent on the chemical intuition and expert input of the researcher. Recent developments in advanced automated methods for reaction path analysis hold promise for eliminating such human-bias from computational catalysis studies. A brief overview of these approaches is presented in the final section of the review. The paper is closed with general concluding remarks.
Collapse
Affiliation(s)
| | | | - Justin K. Kirkland
- Department
of Chemistry, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Jacob Townsend
- Department
of Chemistry, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Ali Hashemi
- Inorganic
Systems Engineering group, Department of Chemical Engineering, Delft University of Technology, Van der Maasweg 9, 2629 HZ Delft, The Netherlands
| | - Chong Liu
- Inorganic
Systems Engineering group, Department of Chemical Engineering, Delft University of Technology, Van der Maasweg 9, 2629 HZ Delft, The Netherlands
| | - Evgeny A. Pidko
- TheoMAT
group, ITMO University, Lomonosova 9, St. Petersburg 191002, Russia
- Inorganic
Systems Engineering group, Department of Chemical Engineering, Delft University of Technology, Van der Maasweg 9, 2629 HZ Delft, The Netherlands
| |
Collapse
|
20
|
Chen M, Baer R, Neuhauser D, Rabani E. Overlapped embedded fragment stochastic density functional theory for covalently-bonded materials. J Chem Phys 2019; 150:034106. [DOI: 10.1063/1.5064472] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ming Chen
- Department of Chemistry, University of California, and Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Roi Baer
- Fritz Haber Center for Molecular Dynamics and 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
| | - Eran Rabani
- Department of Chemistry, University of California, and Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
- The Raymond and Beverly Sackler Center of Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv 69978, Israel
| |
Collapse
|
21
|
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
|
22
|
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
|
23
|
Hégely B, Nagy PR, Kállay M. Dual Basis Set Approach for Density Functional and Wave Function Embedding Schemes. J Chem Theory Comput 2018; 14:4600-4615. [DOI: 10.1021/acs.jctc.8b00350] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Affiliation(s)
- Bence Hégely
- MTA-BME Lendület Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Péter R. Nagy
- MTA-BME Lendület Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Mihály Kállay
- MTA-BME Lendület Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| |
Collapse
|
24
|
Culpitt T, Brorsen KR, Hammes-Schiffer S. Communication: Density functional theory embedding with the orthogonality constrained basis set expansion procedure. J Chem Phys 2018; 146:211101. [PMID: 28576084 DOI: 10.1063/1.4984777] [Citation(s) in RCA: 33] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Density functional theory (DFT) embedding approaches have generated considerable interest in the field of computational chemistry because they enable calculations on larger systems by treating subsystems at different levels of theory. To circumvent the calculation of the non-additive kinetic potential, various projector methods have been developed to ensure the orthogonality of molecular orbitals between subsystems. Herein the orthogonality constrained basis set expansion (OCBSE) procedure is implemented to enforce this subsystem orbital orthogonality without requiring a level shifting parameter. This scheme is a simple alternative to existing parameter-free projector-based schemes, such as the Huzinaga equation. The main advantage of the OCBSE procedure is that excellent convergence behavior is attained for DFT-in-DFT embedding without freezing any of the subsystem densities. For the three chemical systems studied, the level of accuracy is comparable to or higher than that obtained with the Huzinaga scheme with frozen subsystem densities. Allowing both the high-level and low-level DFT densities to respond to each other during DFT-in-DFT embedding calculations provides more flexibility and renders this approach more generally applicable to chemical systems. It could also be useful for future extensions to embedding approaches combining wavefunction theories and DFT.
Collapse
Affiliation(s)
- Tanner Culpitt
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
| | - Kurt R Brorsen
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
| | - Sharon Hammes-Schiffer
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
| |
Collapse
|
25
|
Huang C, Chi YC. Directly patching high-level exchange-correlation potential based on fully determined optimized effective potentials. J Chem Phys 2017; 147:244111. [PMID: 29289130 DOI: 10.1063/1.5003663] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The key element in Kohn-Sham (KS) density functional theory is the exchange-correlation (XC) potential. We recently proposed the exchange-correlation potential patching (XCPP) method with the aim of directly constructing high-level XC potential in a large system by patching the locally computed, high-level XC potentials throughout the system. In this work, we investigate the patching of the exact exchange (EXX) and the random phase approximation (RPA) correlation potentials. A major challenge of XCPP is that a cluster's XC potential, obtained by solving the optimized effective potential equation, is only determined up to an unknown constant. Without fully determining the clusters' XC potentials, the patched system's XC potential is "uneven" in the real space and may cause non-physical results. Here, we developed a simple method to determine this unknown constant. The performance of XCPP-RPA is investigated on three one-dimensional systems: H20, H10Li8, and the stretching of the H19-H bond. We investigated two definitions of EXX: (i) the definition based on the adiabatic connection and fluctuation dissipation theorem (ACFDT) and (ii) the Hartree-Fock (HF) definition. With ACFDT-type EXX, effective error cancellations were observed between the patched EXX and the patched RPA correlation potentials. Such error cancellations were absent for the HF-type EXX, which was attributed to the fact that for systems with fractional occupation numbers, the integral of the HF-type EXX hole is not -1. The KS spectra and band gaps from XCPP agree reasonably well with the benchmarks as we make the clusters large.
Collapse
Affiliation(s)
- Chen Huang
- Department of Scientific Computing, Florida State University, Tallahassee, Florida 32306-4120, USA
| | - Yu-Chieh Chi
- Department of Scientific Computing, Florida State University, Tallahassee, Florida 32306-4120, USA
| |
Collapse
|
26
|
Lee SJ, Miyamoto K, Ding F, Manby FR, Miller TF. Density-based errors in mixed-basis mean-field electronic structure, with implications for embedding and QM/MM methods. Chem Phys Lett 2017. [DOI: 10.1016/j.cplett.2017.04.059] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
|
27
|
Ding F, Tsuchiya T, Manby FR, Miller TF. Linear-Response Time-Dependent Embedded Mean-Field Theory. J Chem Theory Comput 2017; 13:4216-4227. [DOI: 10.1021/acs.jctc.7b00666] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Feizhi Ding
- Division
of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Takashi Tsuchiya
- Centre
for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - 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
|
28
|
Koh KJ, Nguyen-Beck TS, Parkhill J. Accelerating Realtime TDDFT with Block-Orthogonalized Manby–Miller Embedding Theory. J Chem Theory Comput 2017; 13:4173-4178. [DOI: 10.1021/acs.jctc.7b00494] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Kevin J. Koh
- Department of Chemistry and
Biochemistry, The University of Notre Dame, 251 Nieuwland Science Hall, Notre Dame, Indiana 46556, United States
| | - Triet S. Nguyen-Beck
- Department of Chemistry and
Biochemistry, The University of Notre Dame, 251 Nieuwland Science Hall, Notre Dame, Indiana 46556, United States
| | - John Parkhill
- Department of Chemistry and
Biochemistry, The University of Notre Dame, 251 Nieuwland Science Hall, Notre Dame, Indiana 46556, United States
| |
Collapse
|
29
|
Libisch F, Marsman M, Burgdörfer J, Kresse G. Embedding for bulk systems using localized atomic orbitals. J Chem Phys 2017; 147:034110. [DOI: 10.1063/1.4993795] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- F. Libisch
- Institute for Theoretical Physics, Vienna University of Technology, Wiedner Hauptstraße 8-10/136, A-1040 Vienna, Austria, EU
| | - M. Marsman
- Faculty of Physics and Center for Computational Materials Science, University of Vienna, Sensengasse 8/12, A-1090 Vienna, Austria, EU
| | - J. Burgdörfer
- Institute for Theoretical Physics, Vienna University of Technology, Wiedner Hauptstraße 8-10/136, A-1040 Vienna, Austria, EU
| | - G. Kresse
- Faculty of Physics and Center for Computational Materials Science, University of Vienna, Sensengasse 8/12, A-1090 Vienna, Austria, EU
| |
Collapse
|
30
|
Hršak D, Olsen JMH, Kongsted J. Optimization and transferability of non-electrostatic repulsion in the polarizable density embedding model. J Comput Chem 2017. [DOI: 10.1002/jcc.24859] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Dalibor Hršak
- Department of Physics, Chemistry and Pharmacy; University of Southern Denmark; Campusvej 55 Odense M 5230 Denmark
| | - Jógvan Magnus Haugaard Olsen
- Department of Physics, Chemistry and Pharmacy; University of Southern Denmark; Campusvej 55 Odense M 5230 Denmark
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy; University of Southern Denmark; Campusvej 55 Odense M 5230 Denmark
| |
Collapse
|
31
|
Ding F, Manby FR, Miller TF. Embedded Mean-Field Theory with Block-Orthogonalized Partitioning. J Chem Theory Comput 2017; 13:1605-1615. [DOI: 10.1021/acs.jctc.6b01065] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Feizhi Ding
- 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, U.K
| | - Thomas F. Miller
- Division
of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
32
|
Pennifold RCR, Bennie SJ, Miller TF, Manby FR. Correcting density-driven errors in projection-based embedding. J Chem Phys 2017; 146:084113. [PMID: 28249446 DOI: 10.1063/1.4974929] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Projection-based embedding provides a simple and numerically robust framework for multiscale wavefunction-in-density-functional-theory (WF-in-DFT) calculations. The approach works well when the approximate DFT is sufficiently accurate to describe the energetics of the low-level subsystem and the coupling between subsystems. It is also necessary that the low-level DFT produces a qualitatively reasonable description of the total density, and in this work, we study model systems where delocalization error prevents this from being the case. We find substantial errors in embedding calculations on open-shell doublet systems in which self-interaction errors cause spurious delocalization of the singly occupied orbital. We propose a solution to this error by evaluating the DFT energy using a more accurate self-consistent density, such as that of Hartree-Fock (HF) theory. These so-called WF-in-(HF-DFT) calculations show excellent convergence towards full-system wavefunction calculations.
Collapse
Affiliation(s)
- Robert C R Pennifold
- Center for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Simon J Bennie
- Center 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
| | - Frederick R Manby
- Center for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
33
|
Gómez S, Nafziger J, Restrepo A, Wasserman A. Partition-DFT on the water dimer. J Chem Phys 2017; 146:074106. [DOI: 10.1063/1.4976306] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sara Gómez
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia
| | - Jonathan Nafziger
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, USA
| | - Albeiro Restrepo
- Instituto de Química, Universidad de Antioquia UdeA, Calle 70 No. 52-21, Medellín, Colombia
| | - Adam Wasserman
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, USA
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Avenue, West Lafayette, Indiana 47907, USA
| |
Collapse
|
34
|
Pijeau S, Hohenstein EG. Improved Complete Active Space Configuration Interaction Energies with a Simple Correction from Density Functional Theory. J Chem Theory Comput 2017; 13:1130-1146. [DOI: 10.1021/acs.jctc.6b00893] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Shiela Pijeau
- Department
of Chemistry and Biochemistry, The City College of New York, New York, New York 10031, United States
| | - Edward G. Hohenstein
- Department
of Chemistry and Biochemistry, The City College of New York, New York, New York 10031, United States
- Ph.D.
Program in Chemistry, The City University of New York, New York, New York 10016, United States
| |
Collapse
|
35
|
Biancardi A, Barnes J, Caricato M. Point charge embedding for ONIOM excited states calculations. J Chem Phys 2017; 145:224109. [PMID: 27984901 DOI: 10.1063/1.4972000] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Hybrid quantum mechanical methods can assist in the interpretation and prediction of the electronic spectra of large molecular structures. In this work, we study the performance of the ONIOM (Our own N-layered Integrated molecular Orbital molecular Mechanics) hybrid method for the calculation of transition energies and oscillator strengths by embedding the core region in a field of fixed point charges. These charges introduce polarization effects from the substituent groups to the core region. We test various charge definitions, with particular attention to the issue of overpolarization near the boundary between layers. To minimize this issue, we fit the charges on the electrostatic potential of the entire structure in the presence of the link atoms used to cap dangling bonds. We propose two constrained fitting strategies: one that produces an average set of charges common to both model system calculations, EE(L1), and one that produces two separate sets of embedding charges, EE(L2). The results from our tests show that indeed electronic embedding with constrained-fitted charges tends to improve the performance of ONIOM compared to non-embedded calculations. However, the EE(L2) charges work best for transition energies, and the EE(L1) charges work best for oscillator strengths. This may be an indication that fixed point charges do not have enough flexibility to adapt to each system, and other effects (e.g., polarization of the embedding field) may be necessary.
Collapse
Affiliation(s)
- Alessandro Biancardi
- Department of Chemistry, University of Kansas, 1251 Wescoe Hall Drive, Lawrence, Kansas 66045, USA
| | - Jeremy Barnes
- Department of Chemistry, University of Kansas, 1251 Wescoe Hall Drive, Lawrence, Kansas 66045, USA
| | - Marco Caricato
- Department of Chemistry, University of Kansas, 1251 Wescoe Hall Drive, Lawrence, Kansas 66045, USA
| |
Collapse
|
36
|
Tamukong PK, Khait YG, Hoffmann MR. Accurate Dissociation of Chemical Bonds Using DFT-in-DFT Embedding Theory with External Orbital Orthogonality. J Phys Chem A 2016; 121:256-264. [DOI: 10.1021/acs.jpca.6b09909] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Patrick K. Tamukong
- Chemistry Department, University of North Dakota, Grand
Forks, North Dakota 58202, United States
| | - Yuriy G. Khait
- Chemistry Department, University of North Dakota, Grand
Forks, North Dakota 58202, United States
| | - Mark R. Hoffmann
- Chemistry Department, University of North Dakota, Grand
Forks, North Dakota 58202, United States
| |
Collapse
|
37
|
Miyamoto K, Miller TF, Manby FR. Fock-Matrix Corrections in Density Functional Theory and Use in Embedded Mean-Field Theory. J Chem Theory Comput 2016; 12:5811-5822. [DOI: 10.1021/acs.jctc.6b00685] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Kaito Miyamoto
- 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
| | - Frederick R. Manby
- Centre
for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
38
|
Mills JD, Ben-Nun M, Rollin K, Bromley MWJ, Li J, Hinde RJ, Winstead CL, Sheehy JA, Boatz JA, Langhoff PW. Atomic Spectral Methods for Ab Initio Molecular Electronic Energy Surfaces: Transitioning From Small-Molecule to Biomolecular-Suitable Approaches. J Phys Chem B 2016; 120:8321-37. [PMID: 27232159 DOI: 10.1021/acs.jpcb.6b02021] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Continuing attention has addressed incorportation of the electronically dynamical attributes of biomolecules in the largely static first-generation molecular-mechanical force fields commonly employed in molecular-dynamics simulations. We describe here a universal quantum-mechanical approach to calculations of the electronic energy surfaces of both small molecules and large aggregates on a common basis which can include such electronic attributes, and which also seems well-suited to adaptation in ab initio molecular-dynamics applications. In contrast to the more familiar orbital-product-based methodologies employed in traditional small-molecule computational quantum chemistry, the present approach is based on an "ex-post-facto" method in which Hamiltonian matrices are evaluated prior to wave function antisymmetrization, implemented here in the support of a Hilbert space of orthonormal products of many-electron atomic spectral eigenstates familiar from the van der Waals theory of long-range interactions. The general theory in its various forms incorporates the early semiempirical atoms- and diatomics-in-molecules approaches of Moffitt, Ellison, Tully, Kuntz, and others in a comprehensive mathematical setting, and generalizes the developments of Eisenschitz, London, Claverie, and others addressing electron permutation symmetry adaptation issues, completing these early attempts to treat van der Waals and chemical forces on a common basis. Exact expressions are obtained for molecular Hamiltonian matrices and for associated energy eigenvalues as sums of separate atomic and interaction-energy terms, similar in this respect to the forms of classical force fields. The latter representation is seen to also provide a long-missing general definition of the energies of individual atoms and of their interactions within molecules and matter free from subjective additional constraints. A computer code suite is described for calculations of the many-electron atomic eigenspectra and the pairwise-atomic Hamiltonian matrices required for practical applications. These matrices can be retained as functions of scalar atomic-pair separations and employed in assembling aggregate Hamiltonian matrices, with Wigner rotation matrices providing analytical representations of their angular degrees of freedom. In this way, ab initio potential energy surfaces are obtained in the complete absence of repeated evaluations and transformations of the one- and two-electron integrals at different molecular geometries required in most ab inito molecular calculations, with large Hamiltonian matrix assembly simplified and explicit diagonalizations avoided employing partitioning and Brillouin-Wigner or Rayleigh-Schrödinger perturbation theory. Illustrative applications of the important components of the formalism, selected aspects of the scaling of the approach, and aspects of "on-the-fly" interfaces with Monte Carlo and molecular-dynamics methods are described in anticipation of subsequent applications to biomolecules and other large aggregates.
Collapse
Affiliation(s)
- Jeffrey D Mills
- Air Force Research Laboratory , 10 East Saturn Boulevard, Edwards AFB, California 93524-7680, United States
| | - Michal Ben-Nun
- Predictive Science, Inc. , 9990 Mesa Rim Road #170, San Diego, California 92121, United States
| | - Kyle Rollin
- Northrup Grumman Corporation , 1 Rancho Carmel Drive, San Diego, California 92128, United States
| | - Michael W J Bromley
- School of Mathematics and Physics, The University of Queensland , Brisbane, Queensland 4072, Australia
| | - Jiabo Li
- Accelrys Inc. , 10188 Telesis Court #100, San Diego, California 92121-4779, United States
| | - Robert J Hinde
- Department of Chemistry, University of Tennessee , Knoxville, Tennessee 37996-1600, United States
| | - Carl L Winstead
- A.A. Noyes Laboratory of Chemical Physics, California Institute of Technology , Pasadena, California 91125, United States
| | - Jeffrey A Sheehy
- NASA Headquarters , 300 E Street SW, Suite 5R30, Washington, DC 201546, United States
| | - Jerry A Boatz
- Air Force Research Laboratory , 10 East Saturn Boulevard, Edwards AFB, California 93524-7680, United States
| | - Peter W Langhoff
- Department of Chemistry and Biochemistry, University of California, San Diego , 9500 Gilman Drive, MS 0365, La Jolla, California 92093-0365, United States
| |
Collapse
|
39
|
Hégely B, Nagy PR, Ferenczy GG, Kállay M. Exact density functional and wave function embedding schemes based on orbital localization. J Chem Phys 2016. [DOI: 10.1063/1.4960177] [Citation(s) in RCA: 64] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
40
|
Culpitt T, Brorsen KR, Pak MV, Hammes-Schiffer S. Multicomponent density functional theory embedding formulation. J Chem Phys 2016; 145:044106. [DOI: 10.1063/1.4958952] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Affiliation(s)
- Tanner Culpitt
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Ave, Urbana, Illinois 61801, USA
| | - Kurt R. Brorsen
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Ave, Urbana, Illinois 61801, USA
| | - Michael V. Pak
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Ave, Urbana, Illinois 61801, USA
| | - Sharon Hammes-Schiffer
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Ave, Urbana, Illinois 61801, USA
| |
Collapse
|
41
|
Bennie SJ, van der Kamp MW, Pennifold RCR, Stella M, Manby FR, Mulholland AJ. A Projector-Embedding Approach for Multiscale Coupled-Cluster Calculations Applied to Citrate Synthase. J Chem Theory Comput 2016; 12:2689-97. [PMID: 27159381 DOI: 10.1021/acs.jctc.6b00285] [Citation(s) in RCA: 50] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Projector-based embedding has recently emerged as a robust multiscale method for the calculation of various electronic molecular properties. We present the coupling of projector embedding with quantum mechanics/molecular mechanics modeling and apply it for the first time to an enzyme-catalyzed reaction. Using projector-based embedding, we combine coupled-cluster theory, density-functional theory (DFT), and molecular mechanics to compute energies for the proton abstraction from acetyl-coenzyme A by citrate synthase. By embedding correlated ab initio methods in DFT we eliminate functional sensitivity and obtain high-accuracy profiles in a procedure that is straightforward to apply.
Collapse
Affiliation(s)
- Simon J Bennie
- Center for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, U.K
| | - Marc W van der Kamp
- Center for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, U.K
- School of Biochemistry, Biomedical Sciences Building, University Walk, University of Bristol , Bristol BS8 1TD, U.K
| | - Robert C R Pennifold
- Center for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, U.K
| | - Martina Stella
- Center for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, U.K
| | - Frederick R Manby
- Center for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, U.K
| | - Adrian J Mulholland
- Center for Computational Chemistry, School of Chemistry, University of Bristol , Bristol BS8 1TS, U.K
| |
Collapse
|
42
|
Ren S, Caricato M. Multi-state extrapolation of UV/Vis absorption spectra with QM/QM hybrid methods. J Chem Phys 2016; 144:184102. [PMID: 27179466 DOI: 10.1063/1.4948471] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
In this work, we present a simple approach to simulate absorption spectra from hybrid QM/QM calculations. The goal is to obtain reliable spectra for compounds that are too large to be treated efficiently at a high level of theory. The present approach is based on the extrapolation of the entire absorption spectrum obtained by individual subcalculations. Our program locates the main spectral features in each subcalculation, e.g., band peaks and shoulders, and fits them to Gaussian functions. Each Gaussian is then extrapolated with a formula similar to that of ONIOM (Our own N-layered Integrated molecular Orbital molecular Mechanics). However, information about individual excitations is not necessary so that difficult state-matching across subcalculations is avoided. This multi-state extrapolation thus requires relatively low implementation effort while affording maximum flexibility in the choice of methods to be combined in the hybrid approach. The test calculations show the efficacy and robustness of this methodology in reproducing the spectrum computed for the entire molecule at a high level of theory.
Collapse
Affiliation(s)
- Sijin Ren
- Department of Chemistry, University of Kansas, 1251 Wescoe Hall Dr., Lawrence, Kansas 66045, USA
| | - Marco Caricato
- Department of Chemistry, University of Kansas, 1251 Wescoe Hall Dr., Lawrence, Kansas 66045, USA
| |
Collapse
|
43
|
Bennie SJ, Stella M, Miller TF, Manby FR. Accelerating wavefunction in density-functional-theory embedding by truncating the active basis set. J Chem Phys 2015; 143:024105. [DOI: 10.1063/1.4923367] [Citation(s) in RCA: 55] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Simon J. Bennie
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Martina Stella
- 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
| | - Frederick R. Manby
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
44
|
Artiukhin DG, Jacob CR, Neugebauer J. Excitation energies from frozen-density embedding with accurate embedding potentials. J Chem Phys 2015; 142:234101. [DOI: 10.1063/1.4922429] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Denis G. Artiukhin
- 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
| | - Christoph R. Jacob
- Institute of Physical and Theoretical Chemistry, TU Braunschweig, Hans-Sommer-Straße 10, 38106 Braunschweig, 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
|
45
|
Śmiga S, Fabiano E, Laricchia S, Constantin LA, Della Sala F. Subsystem density functional theory with meta-generalized gradient approximation exchange-correlation functionals. J Chem Phys 2015; 142:154121. [DOI: 10.1063/1.4917257] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Szymon Śmiga
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziadzka 5, 87-100 Torun, Poland
| | - Eduardo Fabiano
- Istituto Nanoscienze-CNR, Euromediterranean Center for Nanomaterial Modelling and Technology (ECMT), Via per Arnesano 16, I-73100 Lecce, Italy
- Center for Biomolecular Nanotechnologies@UNILE, Istituto Italiano di Tecnologia (IIT), Via Barsanti, 73010 Arnesano, LE, Italy
| | - Savio Laricchia
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, USA
| | - Lucian A. Constantin
- Center for Biomolecular Nanotechnologies@UNILE, Istituto Italiano di Tecnologia (IIT), Via Barsanti, 73010 Arnesano, LE, Italy
| | - Fabio Della Sala
- Istituto Nanoscienze-CNR, Euromediterranean Center for Nanomaterial Modelling and Technology (ECMT), Via per Arnesano 16, I-73100 Lecce, Italy
- Center for Biomolecular Nanotechnologies@UNILE, Istituto Italiano di Tecnologia (IIT), Via Barsanti, 73010 Arnesano, LE, Italy
| |
Collapse
|
46
|
Odoh SO, Cramer CJ, Truhlar DG, Gagliardi L. Quantum-Chemical Characterization of the Properties and Reactivities of Metal–Organic Frameworks. Chem Rev 2015; 115:6051-111. [DOI: 10.1021/cr500551h] [Citation(s) in RCA: 206] [Impact Index Per Article: 22.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Samuel O. Odoh
- Department of Chemistry,
Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Christopher J. Cramer
- Department of Chemistry,
Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Donald G. Truhlar
- Department of Chemistry,
Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| | - Laura Gagliardi
- Department of Chemistry,
Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United States
| |
Collapse
|
47
|
Akimov AV, Prezhdo OV. Large-Scale Computations in Chemistry: A Bird’s Eye View of a Vibrant Field. Chem Rev 2015; 115:5797-890. [DOI: 10.1021/cr500524c] [Citation(s) in RCA: 159] [Impact Index Per Article: 17.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Alexey V. Akimov
- Department
of Chemistry, University of South California, Los Angeles, California 90089, United States
| | - Oleg V. Prezhdo
- Department
of Chemistry, University of South California, Los Angeles, California 90089, United States
| |
Collapse
|
48
|
Affiliation(s)
- Martina Stella
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, UK
| | - Simon J. Bennie
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, UK
| | - Frederick R. Manby
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, UK
| |
Collapse
|
49
|
Luber S. Local electric dipole moments for periodic systems via density functional theory embedding. J Chem Phys 2014; 141:234110. [DOI: 10.1063/1.4903828] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
|
50
|
Rusakov AA, Phillips JJ, Zgid D. Local Hamiltonians for quantitative Green's function embedding methods. J Chem Phys 2014; 141:194105. [DOI: 10.1063/1.4901432] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
| | - Jordan J. Phillips
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Dominika Zgid
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| |
Collapse
|