1
|
Demel O, Lecours MJ, Nooijen M. Further investigations into a Laplace MP2 method using range separated Coulomb potential and orbital selective virtuals: Multipole correction, OSV extrapolation, and critical assessment. J Chem Phys 2023; 158:114120. [PMID: 36948803 DOI: 10.1063/5.0135113] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/03/2023] Open
Abstract
We report further investigations to aid the development of a Laplace MP2 (second-order Møller Plesset) method with a range separated Coulomb potential partitioned into short- and long-range parts. The implementation of the method extensively uses sparse matrix algebra, density fitting techniques for the short-range part, and a Fourier transformation in spherical coordinates for the long-range part of the potential. Localized molecular orbitals are employed for the occupied space, whereas virtual space is described by orbital specific virtual orbitals (OSVs) associated with localized molecular orbitals. The Fourier transform is deficient for very large distances between localized occupied orbitals, and a multipole expansion for widely separated pairs is introduced for the direct MP2 contribution, which is applicable also to non-Coulombic potentials that do not satisfy the Laplace equation. For the exchange contribution, an efficient screening of contributing localized occupied pairs is employed, which is discussed more completely here. To mitigate errors due to the truncation of OSVs, a simple and efficient extrapolation procedure is used to obtain results close to MP2 for the full basis set of atomic orbitals Using a suitable set of default parameters, the accuracy of the approach is demonstrated. The current implementation of the approach is not very efficient, and the aim of this paper is to introduce and critically discuss ideas that can have more general applicability beyond MP2 calculations for large molecules.
Collapse
Affiliation(s)
- Ondřej Demel
- J. Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Michael J Lecours
- University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
| | - Marcel Nooijen
- University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
| |
Collapse
|
2
|
Demel O, Lecours MJ, Habrovský R, Nooijen M. Toward Laplace MP2 method using range separated Coulomb potential and orbital selective virtuals. J Chem Phys 2021; 155:154104. [PMID: 34686052 DOI: 10.1063/5.0060099] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We report the development of a new Laplace MP2 (second-order Møller-Plesset) implementation using a range separated Coulomb potential, partitioned into short- and long-range parts. The implementation heavily relies on the use of sparse matrix algebra, density fitting techniques for the short-range Coulomb interactions, while a Fourier transformation in spherical coordinates is used for the long-range part of the potential. Localized molecular orbitals are employed for the occupied space, whereas orbital specific virtual orbitals associated with localized molecular orbitals are obtained from the exchange matrix associated with specific localized occupied orbitals. The range separated potential is crucial to achieve efficient treatment of the direct term in the MP2, while extensive screening is employed to reduce the expense of the exchange contribution in MP2. The focus of this paper is on controllable accuracy and linear scaling of the data entering the algorithm.
Collapse
Affiliation(s)
- Ondřej Demel
- J. Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Michael J Lecours
- University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
| | - Richard Habrovský
- J. Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Marcel Nooijen
- University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
| |
Collapse
|
3
|
Antalík A, Veis L, Brabec J, Demel O, Legeza Ö, Pittner J. Toward the efficient local tailored coupled cluster approximation and the peculiar case of oxo-Mn(Salen). J Chem Phys 2019; 151:084112. [PMID: 31470730 DOI: 10.1063/1.5110477] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Abstract
We introduce a new implementation of the coupled cluster method with single and double excitations tailored by the matrix product state wave functions (DMRG-TCCSD), which employs the local pair natural orbital (LPNO) approach. By exploiting locality in the coupled cluster stage of the calculation, we were able to remove some of the limitations that hindered the application of the canonical version of the method to larger systems and/or with larger basis sets. We assessed the accuracy of the approximation using two systems: tetramethyleneethane (TME) and oxo-Mn(Salen). Using the default cut-off parameters, we were able to recover over 99.7% and 99.8% of the canonical correlation energy for the triplet and singlet state of TME, respectively. In the case of oxo-Mn(Salen), we found that the amount of retrieved canonical correlation energy depends on the size of the complete active space (CAS)-we retrieved over 99.6% for the larger 27 orbital CAS and over 99.8% for the smaller 22 orbital CAS. The use of LPNO-TCCSD allowed us to perform these calculations up to quadruple-ζ basis set, amounting to 1178 basis functions. Moreover, we examined dependence of the ground state of oxo-Mn(Salen) on the CAS composition. We found that the inclusion of 4dxy orbital plays an important role in stabilizing the singlet state at the DMRG-CASSCF level via double-shell effect. However, by including dynamic correlation, the ground state was found to be triplet regardless of the size of the basis set or the composition of CAS, which is in agreement with previous findings by canonical DMRG-TCCSD in smaller basis.
Collapse
Affiliation(s)
- Andrej Antalík
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - 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
| | - Jiří Brabec
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Ondřej Demel
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Örs Legeza
- Strongly Correlated Systems "Lendület" Research Group, Wigner Research Centre for Physics, H-1525 Budapest, Hungary
| | - Jiří Pittner
- 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
|
4
|
Lang J, Brabec J, Saitow M, Pittner J, Neese F, Demel O. Perturbative triples correction to domain-based local pair natural orbital variants of Mukherjee's state specific coupled cluster method. Phys Chem Chem Phys 2019; 21:5022-5038. [PMID: 30762044 DOI: 10.1039/c8cp03577f] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In this article we report an implementation of the perturbative triples correction to Mukherjee's state-specific multireference coupled cluster method based on the domain-based pair natural orbital approach (DLPNO-MkCC). We tested the performance of DLPNO-MkCCSD(T) in calculations involving tetramethyleneethane and isomers of naphthynes. These tests show that more than 97% of triples energy was recovered with respect to the canonical MkCCSD(T) method, which together with the DLPNO-MkCCSD part accounts for about 99.70-99.85% of the total correlation energy. The applicability of the method was demonstrated on calculations of singlet-triplet gaps for several large systems: triangulene, dynemicin A, and a beryllium complex.
Collapse
Affiliation(s)
- Jakub Lang
- J. Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic.
| | | | | | | | | | | |
Collapse
|
5
|
Brabec J, Lang J, Saitow M, Pittner J, Neese F, Demel O. Domain-Based Local Pair Natural Orbital Version of Mukherjee’s State-Specific Coupled Cluster Method. J Chem Theory Comput 2018; 14:1370-1382. [DOI: 10.1021/acs.jctc.7b01184] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jiri Brabec
- J. Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Jakub Lang
- J. Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
- Department of Physical and Macromolecular Chemistry, Faculty of Science, Charles University in Prague, Hlavova 2030, 12840 Prague 2, Czech Republic
| | - Masaaki Saitow
- Max Planck Institute for Chemical Energy Conversion, 45470 Mülheim an der Ruhr, Germany
| | - Jiří Pittner
- J. Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Frank Neese
- Max Planck Institute for Chemical Energy Conversion, 45470 Mülheim an der Ruhr, Germany
| | - Ondřej Demel
- J. Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| |
Collapse
|
6
|
Affiliation(s)
- Jakub Lang
- J. Heyrovsky Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Prague 8, Czech Republic
- Faculty of Science, Charles University in Prague, Prague 2, Czech Republic
| | - Matej Švaňa
- J. Heyrovsky Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Prague 8, Czech Republic
| | - Ondřej Demel
- J. Heyrovsky Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Prague 8, Czech Republic
| | - Jiri Brabec
- J. Heyrovsky Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Prague 8, Czech Republic
| | - Stanislav Kedžuch
- Institute of Inorganic Chemistry, Slovak Academy of Sciences, Bratislava, Slovakia
| | - Jozef Noga
- Institute of Inorganic Chemistry, Slovak Academy of Sciences, Bratislava, Slovakia
- Department of Inorganic Chemistry, Faculty of Natural Sciences, Comenius University, Bratislava, Slovakia
| | - Karol Kowalski
- William R. Wiley Environmental Molecular Sciences Laboratory, Battelle, Pacific Northwest National Laboratory, Richland, WA 99352, USA
| | - Jiří Pittner
- J. Heyrovsky Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Prague 8, Czech Republic
| |
Collapse
|
7
|
Huntington LMJ, Demel O, Nooijen M. Benchmark Applications of Variations of Multireference Equation of Motion Coupled-Cluster Theory. J Chem Theory Comput 2015; 12:114-32. [DOI: 10.1021/acs.jctc.5b00799] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Lee M. J. Huntington
- Department
of Chemistry, University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
| | - Ondřej Demel
- J.
Heyrovský Institute of Physical Chemistry, v. v. i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Marcel Nooijen
- Department
of Chemistry, University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
| |
Collapse
|
8
|
Affiliation(s)
- Ondřej Demel
- J.
Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Jiří Pittner
- J.
Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Frank Neese
- Max Planck Institute
of Chemical Energy Conversion, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
9
|
Nooijen M, Demel O, Datta D, Kong L, Shamasundar KR, Lotrich V, Huntington LM, Neese F. Communication: Multireference equation of motion coupled cluster: A transform and diagonalize approach to electronic structure. J Chem Phys 2014; 140:081102. [DOI: 10.1063/1.4866795] [Citation(s) in RCA: 41] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
|
10
|
Demel O, Kedžuch S, Noga J, Pittner J. Perturbative triples correction for explicitly correlated Mukherjee’s state-specific coupled cluster method. Mol Phys 2013. [DOI: 10.1080/00268976.2013.809488] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
11
|
Demel O, Datta D, Nooijen M. Additional global internal contraction in variations of multireference equation of motion coupled cluster theory. J Chem Phys 2013; 138:134108. [DOI: 10.1063/1.4796523] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
|
12
|
|
13
|
Demel O, Kedžuch S, Švaňa M, Ten-no S, Pittner J, Noga J. An explicitly correlated Mukherjee's state specific coupled cluster method: development and pilot applications. Phys Chem Chem Phys 2012; 14:4753-62. [DOI: 10.1039/c2cp23198k] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
14
|
Kedžuch S, Demel O, Pittner J, Ten-no S, Noga J. Multireference F12 coupled cluster theory: The Brillouin-Wigner approach with single and double excitations. Chem Phys Lett 2011. [DOI: 10.1016/j.cplett.2011.06.023] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
15
|
Bhaskaran-Nair K, Demel O, Šmydke J, Pittner J. Multireference state-specific Mukherjee's coupled cluster method with noniterative triexcitations using uncoupled approximation. J Chem Phys 2011; 134:154106. [DOI: 10.1063/1.3573373] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
16
|
Demel O, Bhaskaran-Nair K, Pittner J. Uncoupled multireference state-specific Mukherjee’s coupled cluster method with triexcitations. J Chem Phys 2010; 133:134106. [DOI: 10.1063/1.3495679] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
17
|
Bhaskaran-Nair K, Demel O, Pittner J. Multireference Mukherjee’s coupled cluster method with triexcitations in the linked formulation: Efficient implementation and applications. J Chem Phys 2010; 132:154105. [DOI: 10.1063/1.3376053] [Citation(s) in RCA: 47] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
|
18
|
Bhaskaran-Nair K, Demel O, Pittner J. Multireference state-specific Mukherjee’s coupled cluster method with noniterative triexcitations. J Chem Phys 2008; 129:184105. [DOI: 10.1063/1.3006401] [Citation(s) in RCA: 98] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
|
19
|
Demel O, Shamasundar KR, Kong L, Nooijen M. Application of Double Ionization State-Specific Equation of Motion Coupled Cluster Method to Organic Diradicals. J Phys Chem A 2008; 112:11895-902. [PMID: 18956856 DOI: 10.1021/jp800577q] [Citation(s) in RCA: 66] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ondřej Demel
- Department of Chemistry, University of Waterloo, 200 University Avenue West, Waterloo, Ontario, N2L 3G1 Canada
| | - K. R. Shamasundar
- Department of Chemistry, University of Waterloo, 200 University Avenue West, Waterloo, Ontario, N2L 3G1 Canada
| | - Liguo Kong
- Department of Chemistry, University of Waterloo, 200 University Avenue West, Waterloo, Ontario, N2L 3G1 Canada
| | - Marcel Nooijen
- Department of Chemistry, University of Waterloo, 200 University Avenue West, Waterloo, Ontario, N2L 3G1 Canada
| |
Collapse
|
20
|
Demel O, Pittner J. Multireference Brillouin–Wigner coupled cluster method with singles, doubles, and triples: Efficient implementation and comparison with approximate approaches. J Chem Phys 2008; 128:104108. [DOI: 10.1063/1.2832865] [Citation(s) in RCA: 57] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
21
|
Pittner J, Demel O. Publisher’s Note: “Towards the multireference Brillouin-Wigner coupled-clusters method with iterative connected triples: MR BWCCSDT-α approximation” [J. Chem. Phys. 122, 181101 (2005)]. J Chem Phys 2005. [DOI: 10.1063/1.2006670] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jiři Pittner
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Ondřej Demel
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| |
Collapse
|
22
|
Demel O, Pittner J, Čársky P, Hubač I. Multireference Brillouin−Wigner Coupled Cluster Singles and Doubles Study of the Singlet−Triplet Separation in Alkylcarbenes. J Phys Chem A 2004. [DOI: 10.1021/jp037135m] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Ondřej Demel
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic, and Department of Chemical Physics, Faculty of Mathematics and Physics, Comenius University, 84215 Bratislava, Slovakia
| | - Jiří Pittner
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic, and Department of Chemical Physics, Faculty of Mathematics and Physics, Comenius University, 84215 Bratislava, Slovakia
| | - Petr Čársky
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic, and Department of Chemical Physics, Faculty of Mathematics and Physics, Comenius University, 84215 Bratislava, Slovakia
| | - Ivan Hubač
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic, and Department of Chemical Physics, Faculty of Mathematics and Physics, Comenius University, 84215 Bratislava, Slovakia
| |
Collapse
|