1
|
Drwal D, Pernal K, Pastorczak E. Multireference Correlated Oscillator Strengths from Adiabatic Connection Approaches Based on Extended Random Phase Approximation. J Chem Theory Comput 2024; 20:3659-3668. [PMID: 38669448 PMCID: PMC11099974 DOI: 10.1021/acs.jctc.4c00103] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2024] [Revised: 03/23/2024] [Accepted: 04/05/2024] [Indexed: 04/28/2024]
Abstract
We show that accurate oscillator strengths can be obtained from adiabatic connection (AC) approaches based on the extended random phase approximation (ERPA) combined with multireference (complete active space, CAS) wave functions. The oscillator strengths calculated using the perturbation-corrected ERPA transition density matrices, proposed in this work, and the excitation energies calculated with recently introduced AC correlation energy methods, AC0 and AC0D, compete with accuracy in the perturbational CASPT2 approach and require less computational effort. AC0 and AC0D methods scale more favorably with the number of active orbitals than multiconfigurational perturbation approaches like CASPT2 and NEVPT2 thanks to their dependence on reduced density matrices up to the order of 2. Importantly, the newly developed approach for computing correlated transition dipole moments does not entail any additional costs, as all intermediate quantities become available when AC0 energies are being computed. We also test the performance of the recently proposed AC method corrected for the negative-transition contributions to the correlation energy, AC0D, for triplet excitation energies. Similarly, as for the singlet excitations, the correction improves the performance of the AC0 method, particularly for the low-lying excited states.
Collapse
Affiliation(s)
- Daria Drwal
- 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
| | - Ewa Pastorczak
- Institute of Physics, Lodz
University of Technology, ul. Wolczanska 217/221, 93-005 Lodz, Poland
| |
Collapse
|
2
|
Tecmer P, Gałyńska M, Szczuczko L, Boguslawski K. Geminal-Based Strategies for Modeling Large Building Blocks of Organic Electronic Materials. J Phys Chem Lett 2023; 14:9909-9917. [PMID: 37903084 PMCID: PMC10641881 DOI: 10.1021/acs.jpclett.3c02434] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2023] [Revised: 10/09/2023] [Accepted: 10/24/2023] [Indexed: 11/01/2023]
Abstract
We elaborate on unconventional electronic structure methods based on geminals and their potential to advance the rapidly developing field of organic photovoltaics (OPVs). Specifically, we focus on the computational advantages of geminal-based methods over standard approaches and identify the critical aspects of OPV development. Examples are reliable and efficient computations of orbital energies, electronic spectra, and van der Waals interactions. Geminal-based models can also be combined with quantum embedding techniques and a quantum information analysis of orbital interactions to gain a fundamental understanding of the electronic structures and properties of realistic OPV building blocks. Furthermore, other organic components present in, for instance, dye-sensitized solar cells (DSSCs) represent another promising scope of application. Finally, we provide numerical examples predicting the properties of a small building block of OPV components and two carbazole-based dyes proposed as possible DSSC sensitizers.
Collapse
Affiliation(s)
- Paweł Tecmer
- Institute of Physics, Faculty of Physics,
Astronomy, and Informatics, Nicolaus Copernicus
University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland
| | - Marta Gałyńska
- Institute of Physics, Faculty of Physics,
Astronomy, and Informatics, Nicolaus Copernicus
University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland
| | - Lena Szczuczko
- Institute of Physics, Faculty of Physics,
Astronomy, and Informatics, Nicolaus Copernicus
University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland
| | - Katharina Boguslawski
- Institute of Physics, Faculty of Physics,
Astronomy, and Informatics, Nicolaus Copernicus
University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland
| |
Collapse
|
3
|
Hapka M, Krzemińska A, Modrzejewski M, Przybytek M, Pernal K. Efficient Calculation of the Dispersion Energy for Multireference Systems with Cholesky Decomposition: Application to Excited-State Interactions. J Phys Chem Lett 2023; 14:6895-6903. [PMID: 37494637 PMCID: PMC10405273 DOI: 10.1021/acs.jpclett.3c01568] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2023] [Accepted: 07/17/2023] [Indexed: 07/28/2023]
Abstract
Accurate and efficient prediction of dispersion interactions in excited-state complexes poses a challenge due to the complex nature of electron correlation effects that need to be simultaneously considered. We propose an algorithm for computing the dispersion energy in nondegenerate ground- or excited-state complexes with arbitrary spin. The algorithm scales with the fifth power of the system size due to employing Cholesky decomposition of Coulomb integrals and a recently developed recursive formula for density response functions of the monomers. As a numerical illustration, we apply the new algorithm in the framework of multiconfigurational symmetry adapted perturbation theory, SAPT(MC), to study interactions in dimers with localized excitons. The SAPT(MC) analysis reveals that the dispersion energy may be the main force stabilizing excited-state dimers.
Collapse
Affiliation(s)
- Michał Hapka
- Faculty
of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Agnieszka Krzemińska
- Institute
of Physics, Lodz University of Technology, ul. Wolczanska 217/221, 93-005 Lodz, Poland
| | - Marcin Modrzejewski
- Faculty
of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Michał Przybytek
- 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 217/221, 93-005 Lodz, Poland
| |
Collapse
|
4
|
Żuchowski PS, Moszynski R. Dispersion Energy from the Time-Independent Coupled-Cluster Polarization Propagator. J Chem Theory Comput 2023; 19:1177-1185. [PMID: 36734970 PMCID: PMC9979596 DOI: 10.1021/acs.jctc.2c00902] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
We present a new method of calculation of the dispersion energy in the second-order symmetry-adapted perturbation theory. Using the Longuet-Higgins integral and time-independent coupled-cluster response theory, one shows that the general expression for the dispersion energy can be written in terms of cluster amplitudes and the excitation operators σ, which can be obtained by solving a linear equation. We introduced an approximate scheme dubbed CCPP2(T) for the dispersion energy accurate to the second order of intramonomer correlation, which includes certain classes to be summed to infinity. Assessment of the accuracy of the CCPP2(T) dispersion energy against the FCI dispersion for He2 demonstrates its high accuracy. For more complex systems, CCPP2(T) matches the accuracy of the best methods introduced for calculations of dispersion so far. The method can be extended to higher-order levels of excitations, providing a systematically improvable theory of dispersion interaction.
Collapse
Affiliation(s)
- Piotr S. Żuchowski
- Faculty
of Physics, Astronomy and Applied Informatics, Institute of Physics, Nicolas Copernicus University in Torun, Grudziadzka 5/7, Torun87-100, Poland,E-mail:
| | - Robert Moszynski
- Quantum
Chemistry Laboratory, Faculty of Chemistry, University of Warsaw, Pasteura 1, Warsaw02-093, Poland,E-mail:
| |
Collapse
|
5
|
Tecmer P, Boguslawski K. Geminal-based electronic structure methods in quantum chemistry. Toward a geminal model chemistry. Phys Chem Chem Phys 2022; 24:23026-23048. [PMID: 36149376 DOI: 10.1039/d2cp02528k] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this review, we discuss the recent progress in developing geminal-based theories for challenging problems in quantum chemistry. Specifically, we focus on the antisymmetrized geminal power, generalized valence bond, antisymmetrized product of strongly orthogonal geminals, singlet-type orthogonal geminals, the antisymmetric product of 1-reference orbital geminal, also known as the pair coupled cluster doubles ansatz, and geminals constructed from Richardson-Gaudin states. Furthermore, we review various corrections to account for the missing dynamical correlation effects in geminal models and possible extensions to target electronically excited states and open-shell species. Finally, we discuss some numerical examples and present-day challenges for geminal-based models, including a quantitative and qualitative analysis of wave functions, and software availability.
Collapse
Affiliation(s)
- Paweł Tecmer
- Institute of Physics, Faculty of Physics, Astronomy, and Informatics, Nicolaus Copernicus University in Toruń, Grudziądzka 5, 87-100 Toruń, Poland.
| | - Katharina Boguslawski
- Institute of Physics, Faculty of Physics, Astronomy, and Informatics, Nicolaus Copernicus University in Toruń, Grudziądzka 5, 87-100 Toruń, Poland.
| |
Collapse
|
6
|
Szalewicz K, Jeziorski B. Physical mechanisms of intermolecular interactions from symmetry-adapted perturbation theory. J Mol Model 2022; 28:273. [PMID: 36006512 DOI: 10.1007/s00894-022-05190-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2022] [Accepted: 05/12/2022] [Indexed: 10/15/2022]
Abstract
Symmetry-adapted perturbation theory (SAPT) is a method for computational studies of noncovalent interactions between molecules. This method will be discussed here from the perspective of establishing the paradigm for understanding mechanisms of intermolecular interactions. SAPT interaction energies are obtained as sums of several contributions. Each contribution possesses a clear physical interpretation as it results from some specific physical process. It also exhibits a specific dependence on the intermolecular separation R. The four major contributions are the electrostatic, induction, dispersion, and exchange energies, each due to a different mechanism, valid at any R. In addition, at large R, SAPT interaction energies are seamlessly connected with the corresponding terms in the asymptotic multipole expansion of interaction energy in inverse powers of R. Since such expansion explicitly depends on monomers' multipole moments and polarizabilities, this connection provides additional insights by rigorously relating interaction energies to monomers' properties.
Collapse
Affiliation(s)
- Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, DE, 19716, USA.
| | - Bogumił Jeziorski
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02093, Warsaw, Poland
| |
Collapse
|
7
|
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
|
8
|
Hapka M, Pernal K, Jensen HJA. An efficient implementation of time-dependent linear-response theory for strongly orthogonal geminal wave function models. J Chem Phys 2022; 156:174102. [DOI: 10.1063/5.0082155] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present an implementation of time-dependent linear-response equations for strongly orthogonal geminal wave function models: the generalized valence bond perfect-pairing theory (TD-GVB) and the antisymmetrized product of strongly orthogonal geminals (TD-APSG). The geminal wave functions are optimized using a restricted-step second-order algorithm suitable for handling many geminals, and the linear-response equations are solved in an efficient way using a direct iterative approach. The wave function optimization algorithm features an original scheme to create initial orbitals for the geminal functions in a black-box fashion. The implementation is employed to examine the accuracy of the geminal linear response for singlet excitation energies of small and medium-sized molecules. In systems dominated by dynamic correlation, geminal models constitute only a minor improvement with respect to time-dependent Hartree-Fock. Compared to linear-response complete active space self-consistent field (LR-CASSCF), TD-GVB either misses or gives large errors for states dominated by double excitations.
Collapse
Affiliation(s)
- Michal Hapka
- Faculty of Chemistry, Warsaw University Faculty of Chemistry, Poland
| | - Katarzyna Pernal
- Institute of Physics, Lodz University of Technology Institute of Physics, Poland
| | | |
Collapse
|
9
|
Kowalski PH, Krzemińska A, Pernal K, Pastorczak E. Dispersion Interactions between Molecules in and out of Equilibrium Geometry: Visualization and Analysis. J Phys Chem A 2022; 126:1312-1319. [PMID: 35166552 PMCID: PMC8883464 DOI: 10.1021/acs.jpca.2c00004] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The London dispersion interactions between systems undergoing bond breaking, twisting, or compression are not well studied due to the scarcity and the high computational cost of methods being able to describe both the dynamic correlation and the multireference character of the system. Recently developed methods based on the Generalized Valence Bond wave function, such as EERPA-GVB and SAPT(GVB) (SAPT = symmetry-adapted perturbation theory), allow one to accurately compute and analyze noncovalent interactions between multireference systems. Here, we augment this analysis by introducing a local indicator for dispersion interactions inspired by Mata and Wuttke's Dispersion Interaction Density [ J. Comput. Chem. 2017, 38, 15-23] applied on top of an EERPA-GVB computation. Using a few model systems, we show what insights into the nature and evolution of the dispersion interaction during bond breaking and twisting such an approach is able to offer. The new indicator can be used at a minimal cost additional to an EERPA-GVB computation and can be complemented by an energy decomposition employing the SAPT(GVB) method. We explain the physics behind the initial increase, followed by a decrease in the interaction of linear molecules upon bond stretching. Namely, the elongation of covalent bonds leads to the enhancement of attractive dispersion interactions. For even larger bond lengths, this effect is canceled by the increase of the repulsive exchange forces resulting in a suppression of the interaction and finally leading to repulsion between monomers.
Collapse
Affiliation(s)
- Piotr H Kowalski
- 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
| | - Katarzyna Pernal
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 217/221, 93-005 Lodz, Poland
| | - Ewa Pastorczak
- Institute of Physics, Lodz University of Technology, ul. Wolczanska 217/221, 93-005 Lodz, Poland
| |
Collapse
|
10
|
Hapka M, Przybytek M, Pernal K. Symmetry-Adapted Perturbation Theory Based on Multiconfigurational Wave Function Description of Monomers. J Chem Theory Comput 2021; 17:5538-5555. [PMID: 34517707 PMCID: PMC8444344 DOI: 10.1021/acs.jctc.1c00344] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
We present a formulation
of the multiconfigurational (MC) wave
function symmetry-adapted perturbation theory (SAPT). The method is
applicable to noncovalent interactions between monomers which require
a multiconfigurational description, in particular when the interacting
system is strongly correlated or in an electronically excited state.
SAPT(MC) is based on one- and two-particle reduced density matrices
of the monomers and assumes the single-exchange approximation for
the exchange energy contributions. Second-order terms are expressed
through response properties from extended random phase approximation
(ERPA). The dispersion components of SAPT(MC) have been introduced
in our previous works [HapkaM.J. Chem. Theory Comput.2019, 15, 1016−102730525591; HapkaM.J. Chem. Theory Comput.2019, 15, 6712–672331670950]. SAPT(MC) is applied either with generalized valence
bond perfect pairing (GVB) or with complete active space self-consistent
field (CASSCF) treatment of the monomers. We discuss two model multireference
systems: the H2 ··· H2 dimer
in out-of-equilibrium geometries and interaction between the argon
atom and excited state of ethylene. Using the C2H4* ··· Ar complex as an example, we examine second-order
terms arising from negative transitions in the linear response function
of an excited monomer. We demonstrate that the negative-transition
terms must be accounted for to ensure qualitative prediction of induction
and dispersion energies and develop a procedure allowing for their
computation. Factors limiting the accuracy of SAPT(MC) are discussed
in comparison with other second-order SAPT schemes on a data set of
small single-reference dimers.
Collapse
Affiliation(s)
- 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
| | - Michał Przybytek
- 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
| |
Collapse
|
11
|
Kodrycka M, Patkowski K. Efficient Density-Fitted Explicitly Correlated Dispersion and Exchange Dispersion Energies. J Chem Theory Comput 2021; 17:1435-1456. [PMID: 33606539 DOI: 10.1021/acs.jctc.0c01158] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The leading-order dispersion and exchange-dispersion terms in symmetry-adapted perturbation theory (SAPT), Edisp(20) and Eexch-disp(20), suffer from slow convergence to the complete basis set limit. To alleviate this problem, explicitly correlated variants of these corrections, Edisp(20)-F12 and Eexch-disp(20)-F12, have been proposed recently. However, the original formalism (M., Kodrycka , J. Chem. Theory Comput. 2019, 15, 5965-5986), while highly successful in terms of improving convergence, was not competitive to conventional orbital-based SAPT in terms of computational efficiency due to the need to manipulate several kinds of two-electron integrals. In this work, we eliminate this need by decomposing all types of two-electron integrals using robust density fitting. We demonstrate that the error of the density fitting approximation is negligible when standard auxiliary bases such as aug-cc-pVXZ/MP2FIT are employed. The new implementation allowed us to study all complexes in the A24 database in basis sets up to aug-cc-pV5Z, and the Edisp(20)-F12 and Eexch-disp(20)-F12 values exhibit vastly improved basis set convergence over their conventional counterparts. The well-converged Edisp(20)-F12 and Eexch-disp(20)-F12 numbers can be substituted for conventional Edisp(20) and Eexch-disp(20) ones in a calculation of the total SAPT interaction energy at any level (SAPT0, SAPT2+3, ...). We show that the addition of F12 terms does not improve the accuracy of low-level SAPT treatments. However, when the theory errors are minimized in high-level SAPT approaches such as SAPT2+3(CCD)δMP2, the reduction of basis set incompleteness errors thanks to the F12 treatment substantially improves the accuracy of small-basis calculations.
Collapse
Affiliation(s)
- Monika Kodrycka
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, United States
| | - Konrad Patkowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, Alabama 36849, United States
| |
Collapse
|
12
|
Hapka M, Krzemińska A, Pernal K. How Much Dispersion Energy Is Included in the Multiconfigurational Interaction Energy? J Chem Theory Comput 2020; 16:6280-6293. [PMID: 32877179 PMCID: PMC7586340 DOI: 10.1021/acs.jctc.0c00681] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2020] [Indexed: 11/30/2022]
Abstract
We demonstrate how to quantify the amount of dispersion interaction recovered by supermolecular calculations with the multiconfigurational self-consistent field (MCSCF) wave functions. For this purpose, we present a rigorous derivation which connects the portion of dispersion interaction captured by the assumed wave function model-the residual dispersion interaction-with the size of the active space. Based on the obtained expression for the residual dispersion contribution, we propose a dispersion correction for the MCSCF that avoids correlation double counting. Numerical demonstration for model four-electron dimers in both ground and excited states described with the complete active space self-consistent field (CASSCF) reference serves as a proof-of-concept for the method. Accurate results, largely independent of the size of the active space, are obtained. For many-electron systems, routine CASSCF interaction energy calculations recover a tiny fraction of the full second-order dispersion energy. We found that the residual dispersion is non-negligible only for purely dispersion-bound complexes.
Collapse
Affiliation(s)
- 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
| | - Agnieszka Krzemińska
- 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
|
13
|
Garcia J, Podeszwa R, Szalewicz K. SAPT codes for calculations of intermolecular interaction energies. J Chem Phys 2020; 152:184109. [PMID: 32414261 DOI: 10.1063/5.0005093] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Symmetry-adapted perturbation theory (SAPT) is a method for calculations of intermolecular (noncovalent) interaction energies. The set of SAPT codes that is described here, the current version named SAPT2020, includes virtually all variants of SAPT developed so far, among them two-body SAPT based on perturbative, coupled cluster, and density functional theory descriptions of monomers, three-body SAPT, and two-body SAPT for some classes of open-shell monomers. The properties of systems governed by noncovalent interactions can be predicted only if potential energy surfaces (force fields) are available. SAPT is the preferred approach for generating such surfaces since it is seamlessly connected to the asymptotic expansion of interaction energy. SAPT2020 includes codes for automatic development of such surfaces, enabling generation of complete dimer surfaces with a rigid monomer approximation for dimers containing about one hundred atoms. These codes can also be used to obtain surfaces including internal degrees of freedom of monomers.
Collapse
Affiliation(s)
- Javier Garcia
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| | - Rafał Podeszwa
- Institute of Chemistry, University of Silesia at Katowice, Szkolna 9, Katowice, Poland
| | - Krzysztof Szalewicz
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, USA
| |
Collapse
|
14
|
Abstract
This work examines the suitability of meta-GGA functionals for symmetry-adapted perturbation theory (SAPT) calculations. The assessment is based on the term-by-term comparison with the benchmark SAPT variant based on coupled-cluster singles and doubles description of monomers, SAPT(CCSD). Testing systems include molecular complexes ranging from strong to weak and the He dimer. The following nonempirical meta-GGAs are examined: TPSS, revTPSS, MVS, SCAN, and SCAN0 with and without the asymptotic correction (AC) of the exchange-correlation potential. One range-separated meta-GGA functional, LC-PBETPSS, is also included. The AC-corrected pure meta-GGAs (with the exception of MVS) represent a definite progress in SAPT(DFT) compared to pure GGA, such as PBEAC, with their more consistent predictions of energy components. However, none of the meta-GGAs is better than the hybrid GGA approach SAPT(PBE0AC). The SAPT(DFT) electrostatic energy offers the most sensitive probe of the quality of the underlying DFT density. Both SCAN- and TPSS-based electrostatic energies agree with reference to within 5% or better which is an excellent result. We find that SCAN0 can be used in SAPT without the AC correction. The long-range corrected LC-PBETPSS is a reliable performer both for the components and total interaction energies.
Collapse
|
15
|
Hapka M, Przybytek M, Pernal K. Second-Order Exchange-Dispersion Energy Based on a Multireference Description of Monomers. J Chem Theory Comput 2019; 15:6712-6723. [PMID: 31670950 DOI: 10.1021/acs.jctc.9b00925] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a method for calculation of the second-order exchange-dispersion energy in the framework of the symmetry-adapted perturbation theory (SAPT) for weakly interacting monomers described with multiconfigurational wave functions. The proposed formalism is based on response properties obtained from extended random phase approximation (ERPA) equations and assumes the single-exchange (S2) approximation. The approach is applicable to closed shell systems where static correlation cannot be neglected or to systems in nondegenerate excited states. We examine the new method in combination with either generalized valence bond perfect pairing (GVB) or complete active space self-consistent field (CASSCF) description of the interacting monomers. For model multireference dimers in ground states (H2···H2, Be···Be, He···H2), exchange-dispersion energies are reproduced accurately. For the interaction between the excited hydrogen molecule and the helium atom we found unacceptably large errors which is attributed to the neglect of diagonal double excitations in the employed approximation to the linear response function.
Collapse
Affiliation(s)
- 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
| | - Michał Przybytek
- 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
| |
Collapse
|
16
|
Patkowski K. Recent developments in symmetry‐adapted perturbation theory. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1452] [Citation(s) in RCA: 53] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Konrad Patkowski
- Department of Chemistry and Biochemistry Auburn University Auburn Alabama
| |
Collapse
|
17
|
Brzęk F, Boguslawski K, Tecmer P, Żuchowski PS. Benchmarking the Accuracy of Seniority-Zero Wave Function Methods for Noncovalent Interactions. J Chem Theory Comput 2019; 15:4021-4035. [PMID: 31136703 DOI: 10.1021/acs.jctc.9b00189] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Filip Brzęk
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Torun, Grudziadzka 5, 87-100 Torun, Poland
| | - Katharina Boguslawski
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Torun, Grudziadzka 5, 87-100 Torun, Poland
- Faculty of Chemistry, Nicolaus Copernicus University in Torun, Gagarina 7, 87-100 Torun, Poland
| | - Paweł Tecmer
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Torun, Grudziadzka 5, 87-100 Torun, Poland
| | - Piotr Szymon Żuchowski
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Torun, Grudziadzka 5, 87-100 Torun, Poland
| |
Collapse
|