1
|
Vu N, Mejia-Rodriguez D, Bauman NP, Panyala A, Mutlu E, Govind N, Foley JJ. Cavity Quantum Electrodynamics Complete Active Space Configuration Interaction Theory. J Chem Theory Comput 2024; 20:1214-1227. [PMID: 38291561 PMCID: PMC10876286 DOI: 10.1021/acs.jctc.3c01207] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2023] [Revised: 12/21/2023] [Accepted: 12/26/2023] [Indexed: 02/01/2024]
Abstract
Polariton chemistry has attracted great attention as a potential route to modify chemical structure, properties, and reactivity through strong interactions among molecular electronic, vibrational, or rovibrational degrees of freedom. A rigorous theoretical treatment of molecular polaritons requires the treatment of matter and photon degrees of freedom on equal quantum mechanical footing. In the limit of molecular electronic strong or ultrastrong coupling to one or a few molecules, it is desirable to treat the molecular electronic degrees of freedom using the tools of ab initio quantum chemistry, yielding an approach we refer to as ab initio cavity quantum electrodynamics, where the photon degrees of freedom are treated at the level of cavity quantum electrodynamics. Here, we present an approach called Cavity Quantum Electrodynamics Complete Active Space Configuration Interaction theory to provide ground- and excited-state polaritonic surfaces with a balanced description of strong correlation effects among electronic and photonic degrees of freedom. This method provides a platform for ab initio cavity quantum electrodynamics when both strong electron correlation and strong light-matter coupling are important and is an important step toward computational approaches that yield multiple polaritonic potential energy surfaces and couplings that can be leveraged for ab initio molecular dynamics simulations of polariton chemistry.
Collapse
Affiliation(s)
- Nam Vu
- Department
of Chemistry, University of North Carolina
Charlotte, 9201 University City Blvd., Charlotte, North Carolina 28223, United States
| | - Daniel Mejia-Rodriguez
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Nicholas P. Bauman
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Ajay Panyala
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Erdal Mutlu
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Niranjan Govind
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Jonathan J. Foley
- Department
of Chemistry, University of North Carolina
Charlotte, 9201 University City Blvd., Charlotte, North Carolina 28223, United States
| |
Collapse
|
2
|
Li S, Misiewicz JP, Evangelista FA. Intruder-free cumulant-truncated driven similarity renormalization group second-order multireference perturbation theory. J Chem Phys 2023; 159:114106. [PMID: 37712785 DOI: 10.1063/5.0159403] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2023] [Accepted: 08/14/2023] [Indexed: 09/16/2023] Open
Abstract
Accurate multireference electronic structure calculations are important for constructing potential energy surfaces. Still, even in the case of low-scaling methods, their routine use is limited by the steep growth of the computational and storage costs as the active space grows. This is primarily due to the occurrence of three- and higher-body density matrices or, equivalently, their cumulants. This work examines the effect of various cumulant truncation schemes on the accuracy of the driven similarity renormalization group second-order multireference perturbation theory. We test four different levels of three-body reduced density cumulant truncations that set different classes of cumulant elements to zero. Our test cases include the singlet-triplet gap of CH2, the potential energy curves of the XΣg+1 and AΣu+3 states of N2, and the singlet-triplet splittings of oligoacenes. Our results show that both relative and absolute errors introduced by these cumulant truncations can be as small as 0.5 kcal mol-1 or less. At the same time, the amount of memory required is reduced from O(NA6) to O(NA5), where NA is the number of active orbitals. No additional regularization is needed to prevent the intruder state problem in the cumulant-truncated second-order driven similarity renormalization group multireference perturbation theory methods.
Collapse
Affiliation(s)
- Shuhang Li
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Jonathon P Misiewicz
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
3
|
Jangrouei MR, Krzemińska A, Hapka M, Pastorczak E, Pernal K. Dispersion Interactions in Exciton-Localized States. Theory and Applications to π-π* and n-π* Excited States. J Chem Theory Comput 2022; 18:3497-3511. [PMID: 35587598 PMCID: PMC9202351 DOI: 10.1021/acs.jctc.2c00221] [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/29/2022]
Abstract
![]()
We address the problem
of intermolecular interaction energy calculations
in molecular complexes with localized excitons. Our focus is on the
correct representation of the dispersion energy. We derive an extended
Casimir-Polder formula for direct computation of this contribution
through second order in the intermolecular interaction operator V̂. An alternative formula, accurate to infinite order
in V̂, is derived within the framework of the
adiabatic connection (AC) theory. We also propose a new parametrization
of the VV10 nonlocal correlation density functional, so that it corrects
the CASSCF energy for the dispersion contribution and can be applied
to excited-state complexes. A numerical investigation is carried out
for benzene, pyridine, and peptide complexes with the local exciton
corresponding to the lowest π–π* or n– π*
states. The extended Casimir-Polder formula is implemented in the
framework of multiconfigurational symmetry-adapted perturbation theory,
SAPT(MC). A SAPT(MC) analysis shows that the creation of a localized
exciton affects mostly the electrostatic component of the interaction
energy of investigated complexes. Nevertheless, the changes in Pauli
repulsion and dispersion energies cannot be neglected. We verify the
performance of several perturbation- and AC-based methods. Best results
are obtained with a range-separated variant of an approximate AC approach
employing extended random phase approximation and CASSCF wave functions.
Collapse
Affiliation(s)
- Mohammad Reza Jangrouei
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 217/221, 93-005, Lodz, Poland
| | - Agnieszka Krzemińska
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 217/221, 93-005, Lodz, Poland
| | - Michał Hapka
- Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093, Warsaw, Poland
| | - Ewa Pastorczak
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 217/221, 93-005, Lodz, Poland
| | - Katarzyna Pernal
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 217/221, 93-005, Lodz, Poland
| |
Collapse
|
4
|
Kollmar C, Sivalingam K, Guo Y, Neese F. An efficient implementation of the NEVPT2 and CASPT2 methods avoiding higher-order density matrices. J Chem Phys 2021; 155:234104. [PMID: 34937355 DOI: 10.1063/5.0072129] [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
A factorization of the matrix elements of the Dyall Hamiltonian in N-electron valence state perturbation theory allowing their evaluation with a computational effort comparable to the one needed for the construction of the third-order reduced density matrix at the most is presented. Thus, the computational bottleneck arising from explicit evaluation of the fourth-order density matrix is avoided. It is also shown that the residual terms arising in the case of an approximate complete active space configuration interaction solution and containing even the fifth-order density matrix for two excitation classes can be evaluated with little additional effort by choosing again a favorable factorization of the corresponding matrix elements. An analogous argument is also provided for avoiding the fourth-order density matrix in complete active space second-order perturbation theory. Practical calculations indicate that such an approach leads to a considerable gain in computational efficiency without any compromise in numerical accuracy or stability.
Collapse
Affiliation(s)
- Christian Kollmar
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Kantharuban Sivalingam
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Yang Guo
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, China
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
5
|
Beran P, Matoušek M, Hapka M, Pernal K, Veis L. Density Matrix Renormalization Group with Dynamical Correlation via Adiabatic Connection. J Chem Theory Comput 2021; 17:7575-7585. [PMID: 34762423 DOI: 10.1021/acs.jctc.1c00896] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The quantum chemical version of the density matrix renormalization group (DMRG) method has established itself as one of the methods of choice for calculations of strongly correlated molecular systems. Despite its great ability to capture strong electronic correlation in large active spaces, it is not suitable for computations of dynamical electron correlation. In this work, we present a new approach to the electronic structure problem of strongly correlated molecules, in which DMRG is responsible for a proper description of the strong correlation, whereas dynamical correlation is computed via the recently developed adiabatic connection (AC) technique which requires only up to two-body active space reduced density matrices. We report the encouraging results of this approach on typical candidates for DMRG computations, namely, n-acenes (n = 2 → 7), Fe(II)-porphyrin, and the Fe3S4 cluster.
Collapse
Affiliation(s)
- Pavel Beran
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic.,Faculty of Mathematics and Physics, Charles University, 121 16 Prague, Czech Republic
| | - Mikuláš Matoušek
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic.,Faculty of Mathematics and Physics, Charles University, 121 16 Prague, Czech Republic
| | - Michał Hapka
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland.,Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Katarzyna Pernal
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| | - Libor Veis
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| |
Collapse
|
6
|
Golub P, Antalik A, Veis L, Brabec J. Machine Learning-Assisted Selection of Active Spaces for Strongly Correlated Transition Metal Systems. J Chem Theory Comput 2021; 17:6053-6072. [PMID: 34570505 DOI: 10.1021/acs.jctc.1c00235] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Active space quantum chemical methods could provide very accurate description of strongly correlated electronic systems, which is of tremendous value for natural sciences. The proper choice of the active space is crucial but a nontrivial task. In this article, we present a neural network-based approach for automatic selection of active spaces, focused on transition metal systems. The training set has been formed from artificial systems composed of one transition metal and various ligands, on which we have performed the density matrix renormalization group and calculated the single-site entropy. On the selected set of systems, ranging from small benchmark molecules up to larger challenging systems involving two metallic centers, we demonstrate that our machine learning models could predict the active space orbitals with reasonable accuracy. We also tested the transferability on out-of-the-model systems, including bimetallic complexes and complexes with ligands, which were not involved in the training set. Also, we tested the correctness of the automatically selected active spaces on a Fe(II)-porphyrin model, where we studied the lowest states at the DMRG level and compared the energy difference between spin states or the energy difference between conformations of ferrocene with recent studies.
Collapse
Affiliation(s)
- Pavlo Golub
- J. Heyrovsky Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejskova 3, 18223 Prague 8, Czech Republic
| | - Andrej Antalik
- J. Heyrovsky Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejskova 3, 18223 Prague 8, Czech Republic
| | - Libor Veis
- J. Heyrovsky Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejskova 3, 18223 Prague 8, Czech Republic
| | - Jiri Brabec
- J. Heyrovsky Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejskova 3, 18223 Prague 8, Czech Republic
| |
Collapse
|
7
|
Daas T, Fabiano E, Della Sala F, Gori-Giorgi P, Vuckovic S. Noncovalent Interactions from Models for the Møller-Plesset Adiabatic Connection. J Phys Chem Lett 2021; 12:4867-4875. [PMID: 34003655 PMCID: PMC8280728 DOI: 10.1021/acs.jpclett.1c01157] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2021] [Accepted: 05/13/2021] [Indexed: 05/08/2023]
Abstract
Given the omnipresence of noncovalent interactions (NCIs), their accurate simulations are of crucial importance across various scientific disciplines. Here we construct accurate models for the description of NCIs by an interpolation along the Møller-Plesset adiabatic connection (MP AC). Our interpolation approximates the correlation energy, by recovering MP2 at small coupling strengths and the correct large-coupling strength expansion of the MP AC, recently shown to be a functional of the Hartree-Fock density. Our models are size consistent for fragments with nondegenerate ground states, have the same cost as double hybrids, and require no dispersion corrections to capture NCIs accurately. These interpolations greatly reduce large MP2 errors for typical π-stacking complexes (e.g., benzene-pyridine dimers) and for the L7 data set. They are also competitive with state-of-the-art dispersion enhanced functionals and can even significantly outperform them for a variety of data sets, such as CT7 and L7.
Collapse
Affiliation(s)
- Timothy
J. Daas
- Department
of Chemistry & Pharmaceutical Sciences and Amsterdam Institute
of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Eduardo Fabiano
- Institute
for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, 73100 Lecce, Italy
- Center
for Biomolecular Nanotechnologies, Istituto
Italiano di Tecnologia, Via Barsanti 14, 73010 Arnesano (LE), Italy
| | - Fabio Della Sala
- Institute
for Microelectronics and Microsystems (CNR-IMM), Via Monteroni, Campus Unisalento, 73100 Lecce, Italy
- Center
for Biomolecular Nanotechnologies, Istituto
Italiano di Tecnologia, Via Barsanti 14, 73010 Arnesano (LE), Italy
| | - Paola Gori-Giorgi
- Department
of Chemistry & Pharmaceutical Sciences and Amsterdam Institute
of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
| | - Stefan Vuckovic
- Department
of Chemistry & Pharmaceutical Sciences and Amsterdam Institute
of Molecular and Life Sciences (AIMMS), Faculty of Science, Vrije Universiteit, De Boelelaan 1083, 1081HV Amsterdam, The Netherlands
- Physical
and Theoretical Chemistry, University of
Saarland, 66123 Saarbrücken, Germany
- Department
of Chemistry, University of California, Irvine, California 92697, United States
| |
Collapse
|
8
|
Drwal D, Pastorczak E, Pernal K. Excited states in the adiabatic connection fluctuation-dissipation theory: Recovering missing correlation energy from the negative part of the density response spectrum. J Chem Phys 2021; 154:164102. [PMID: 33940850 DOI: 10.1063/5.0046852] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023] Open
Abstract
The adiabatic connection (AC) theory offers an alternative to the perturbation theory methods for computing correlation energy in the multireference wavefunction framework. We show that the AC correlation energy formula can be expressed in terms of the density linear response function as a sum of components related to positive and negative parts of the transition energy spectrum. Consequently, generalization of the adiabatic connection fluctuation-dissipation theory to electronically excited states is obtained. The component of the linear response function related to the negative-transition energy enters the correlation energy expression with an opposite sign to that of the positive-transition part and is non-negligible in the description of excited states. To illustrate this, we analyze the approximate AC model in which the linear response function is obtained in the extended random phase approximation (ERPA). We demonstrate that AC can be successfully combined with the ERPA for excited states, provided that the negative-excitation component of the response function is rigorously accounted for. The resulting AC0D model, an extension of the AC0 scheme introduced in our earlier works, is applied to a benchmark set of singlet excitation energies of organic molecules. AC0D constitutes a significant improvement over AC0 by bringing the excitation energies of the lowest excited states to a satisfactory agreement with theoretical best estimates, which parallels or even exceeds the accuracy of the n-electron valence state perturbation theory method. For higher excitations, AC0D is less reliable due to the gradual deterioration of the underlying ERPA linear response.
Collapse
Affiliation(s)
- Daria Drwal
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| | - Ewa Pastorczak
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| | - Katarzyna Pernal
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 219, 90-924 Lodz, Poland
| |
Collapse
|