1
|
Cioslowski J, Strasburger K. Constraints upon Functionals of the 1-Matrix, Universal Properties of Natural Orbitals, and the Fallacy of the Collins "Conjecture". J Phys Chem Lett 2024; 15:1328-1337. [PMID: 38285733 PMCID: PMC10860149 DOI: 10.1021/acs.jpclett.3c03118] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2023] [Revised: 01/02/2024] [Accepted: 01/04/2024] [Indexed: 01/31/2024]
Abstract
Reliability of quantum-chemical calculations based upon the density functional theory and its 1-matrix counterpart hinges upon minimizing the extent of empirical parameterization in the approximate energy expressions of these formalisms while imposing as many rigorous constraints upon them as possible. The recently uncovered universal properties of the natural orbitals facilitate the construction of such constraints for the 1-matrix functionals. The benefits of their employment in the validation of these functionals are vividly demonstrated by a critical review of the three incarnations of the so-called Collins conjecture. Although the incorporation of rigorous definitions of the correlation energy and entropy, and the identification of individual potential energy hypersurfaces as probable domains of its applicability turn the originally published unsubstantiated claim into a proper conjecture, the resulting formalism is found to be merely a conduit for incorporation of static correlation effects in electronic structure calculations that is unlikely to allow attaining chemical accuracy.
Collapse
Affiliation(s)
- Jerzy Cioslowski
- Institute
of Physics, University of Szczecin, Wielkopolska 15, 70-451 Szczecin, Poland
- Max-Planck-Institut
für Physik komplexer Systeme, Nöthnitzer Straße 38, 01187 Dresden, Germany
| | - Krzysztof Strasburger
- Department
of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology, Wybrzeże Wyspiańskiego
27, 50-370 Wrocław, Poland
| |
Collapse
|
2
|
Cioslowski J, Strasburger K. Symmetry Equiincidence of Natural Orbitals. J Phys Chem Lett 2023; 14:9296-9303. [PMID: 37815811 DOI: 10.1021/acs.jpclett.3c01738] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/11/2023]
Abstract
The symmetry equiincidence principle quantifies the apportionment of the natural orbitals (NOs), ordered according to their nonascending occupation numbers, among the irreducible representations (irreps) of the point group pertaining to the underlying on-top two-electron density. This principle, which is rigorously proven for the resolvable Cs, C2v, C3v, C4v, C6v, D2h, D3h, D4h, D6h, and Oh point groups, states that the symmetry incidences, i.e., the asymptotic probabilities with which the NOs belonging to different irreps occur, are proportional to the squares of irreps' dimensions. Since its proof hinges upon a sufficient number of planes of symmetry among the elements of a given point group, it yields only linear combinations of the symmetry incidences for the quasiresolvable groups with too few such planes and fails for the unresolvable C1, Ci, Cn, Dn, S2n, T, O, and I groups whose nontrivial elements comprise only symmetry axes and/or the center of inversion.
Collapse
Affiliation(s)
- Jerzy Cioslowski
- Institute of Physics, University of Szczecin, Wielkopolska 15, 70-451 Szczecin, Poland
| | - Krzysztof Strasburger
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology, Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland
| |
Collapse
|
3
|
Cioslowski J, Englert BG, Trappe MI, Hue JH. Contactium: A strongly correlated model system. J Chem Phys 2023; 158:2890212. [PMID: 37171198 DOI: 10.1063/5.0150523] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2023] [Accepted: 04/21/2023] [Indexed: 05/13/2023] Open
Abstract
At the limit of an infinite confinement strength ω, the ground state of a system that comprises two fermions or bosons in harmonic confinement interacting through the Fermi-Huang pseudopotential remains strongly correlated. A detailed analysis of the one-particle description of this "contactium" reveals several peculiarities that are not encountered in conventional model systems (such as the two-electron harmonium atom, ballium, and spherium) involving Coulombic interparticle interactions. First of all, none of the natural orbitals (NOs) {ψn(ω;r)} of the contactium is unoccupied, which implies nonzero collective occupancies for all the angular momenta. Second, the NOs and their non-ascendingly ordered occupation numbers {νn} turn out to be related to the eigenfunctions and eigenvalues of a zero-energy Schrödinger equation with an attractive Gaussian potential. This observation enables the derivation of their properties, such as the n-4/3 asymptotic decay of νn at the n→∞ limit (which differs from that of n-8/3 in the Coulombic systems), the independence of the confinement energy vn=⟨ψn(ω;r)|12ω2r2|ψn(ω;r)⟩ of n, and the n-2/3 asymptotic decay of the respective contribution νntn to the kinetic energy. Upon suitable scaling, the weakly occupied NOs of the contactium turn out to be virtually identical to those of the two-electron harmonium atom at the ω → ∞ limit, despite the entirely different interparticle interactions in these systems.
Collapse
Affiliation(s)
- Jerzy Cioslowski
- Institute of Physics, University of Szczecin, Wielkopolska 15, 70-451 Szczecin, Poland
| | - Berthold-Georg Englert
- Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, Singapore 117543, Singapore
- Department of Physics, National University of Singapore, 2 Science Drive 3, Singapore 117542, Singapore
- MajuLab, CNRS-UCA-SU-NUS-NTU International, Joint Research Unit, Singapore
| | - Martin-Isbjörn Trappe
- Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, Singapore 117543, Singapore
| | - Jun Hao Hue
- Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, Singapore 117543, Singapore
| |
Collapse
|
4
|
Cioslowski J, Schilling C, Schilling R. 1-Matrix functional for long-range interaction energy of two hydrogen atoms. J Chem Phys 2023; 158:084106. [PMID: 36859076 DOI: 10.1063/5.0139897] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
Abstract
The leading terms in the large-R asymptotics of the functional of the one-electron reduced density matrix for the ground-state energy of the H2 molecule with the internuclear separation R are derived thanks to the solution of the phase dilemma at the R → ∞ limit. At this limit, the respective natural orbitals (NOs) are given by symmetric and antisymmetric combinations of "half-space" orbitals with the corresponding natural amplitudes having the same amplitudes but opposite signs. Minimization of the resulting explicit functional yields the large-R asymptotics for the occupation numbers of the weakly occupied NOs and the C6 dispersion coefficient. The highly accurate approximates for the radial components of the p-type "half-space" orbitals and the corresponding occupation numbers (that decay like R-6), which are available for the first time thanks to the development of the present formalism, have some unexpected properties.
Collapse
Affiliation(s)
- Jerzy Cioslowski
- Institute of Physics, University of Szczecin, Wielkopolska 15, 70-451 Szczecin, Poland
| | - Christian Schilling
- Department of Physics, Arnold Sommerfeld Center for Theoretical Physics, Ludwig-Maximilians-Universität München, Theresienstrasse 37, 80333 München, Germany
| | - Rolf Schilling
- Institute of Physics, Johannes Gutenberg University, Staudinger Weg 9, 55099 Mainz, Germany
| |
Collapse
|
5
|
Savin A, Karwowski J. Correcting Models with Long-Range Electron Interaction Using Generalized Cusp Conditions. J Phys Chem A 2023; 127:1377-1385. [PMID: 36720050 DOI: 10.1021/acs.jpca.2c08426] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Sources of energy errors resulting from the replacement of the physical Coulomb interaction by its long-range erfc(μr)/r approximation are explored. It is demonstrated that the results can be dramatically improved and the range of μ giving energies within chemical accuracy limits significantly extended if the generalized cusp conditions are used to represent the wave function at small r. The numerical results for two-electron harmonium are presented and discussed.
Collapse
Affiliation(s)
- Andreas Savin
- Laboratoire de Chimie Théorique, CNRS and Sorbonne University4 Place Jussieu, 75252Paris cedex 05, France
| | - Jacek Karwowski
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Grudziadzka 5, 87-100Toruń, Poland
| |
Collapse
|
6
|
Cioslowski J, Strasburger K. A Universal Power Law Governing the Accuracy of Wave Function-Based Electronic Structure Calculations. J Phys Chem Lett 2022; 13:8055-8061. [PMID: 35994623 DOI: 10.1021/acs.jpclett.2c01987] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
A universal power law governing the accuracy of wave function-based electronic structure calculations is derived from first principles. The resulting expression ΔE(N,N)N≳19π2gNN, where g is a system-specific factor assuming values between zero and one and ≳ stands for asymptotic inequality at the limit of N→∞, allows facile estimation of the error ΔE(N,N) in the electronic energy of a singlet state of an N-electron system computed with a basis set of N one-electron functions. Several approaches to the estimation of the factor g, which depends on the on-top two-electron density, are presented.
Collapse
Affiliation(s)
- Jerzy Cioslowski
- Institute of Physics, University of Szczecin, Wielkopolska 15, 70-451 Szczecin, Poland
| | - Krzysztof Strasburger
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology, Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland
| |
Collapse
|
7
|
Cioslowski J. ERRORS IN APPROXIMATE IONIZATION ENERGIES DUE TO THE ONE-ELECTRON SPACE TRUNCATION OF THE EKT EIGENPROBLEM. J Chem Phys 2022; 157:064109. [DOI: 10.1063/5.0096145] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Unless the approximate wavefunction of the parent system is expressed in terms of explicitly correlated basis functions, the finite size of the generalized Fock matrix is unlikely to be the leading source of the truncation error in the ionization energy E produced by the EKT (extended Koopmans' theorem) formalism. This conclusion is drawn from a rigorous analysis that involves error partitioning into the parent- and ionized-system contributions, the former being governed by asymptotic power laws when the underlying wavefunction is assembled from a large number of spinorbitals and the latter arising from the truncation of the infinite-dimensional matrix V whose elements involve the 1-, 2- and 3-matrices of the parent system. Quite surprisingly, the decay of the second contribution with the number $n$ of the natural spinorbitals (NOs) employed in the construction of the truncated V turns out to be strongly system-dependent even in the simplest case of the 1S states of two-electron systems, following the n-5 power law for the helium atom while exhibiting an erratic behavior for the H- anion. This phenomenon, which stems from the presence of the so-called solitonic natural spinorbitals among the NOs, renders the extrapolation of the EKT approximates of E to the complete-basis-set limit generally unfeasible. However, attaining that limit is not contingent upon attempted reproduction of the ill-defined one-electron function known as 'the removal orbital', which does not have to be invoked in the derivation of EKT and whose expansion in terms of the NOs diverges.
Collapse
|
8
|
Karwowski J, Savin A. Two-particle coalescence conditions revisited. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2069055] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Jacek Karwowski
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University, Torun, Poland
| | - Andreas Savin
- Laboratoire de Chimie Théorique, CNRS and Sorbonne University, Paris cedex 05, France
| |
Collapse
|
9
|
Cioslowski J, Prątnicki F, Strasburger K. Solitonic natural orbitals in Coulombic systems. J Chem Phys 2022; 156:034108. [DOI: 10.1063/5.0075986] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Jerzy Cioslowski
- Institute of Physics, University of Szczecin, Wielkopolska 15, 70-451 Szczecin, Poland
| | - Filip Prątnicki
- Institute of Physics, University of Szczecin, Wielkopolska 15, 70-451 Szczecin, Poland
| | - Krzysztof Strasburger
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology, Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland
| |
Collapse
|