1
|
Alcoba DR, Oña OB, Torre A, Lain L, Sierra G, Massaccesi GE. A variance-based optimization for determining ground and excited N-electron wave functions within the doubly occupied configuration interaction scheme. J Chem Phys 2024; 160:164107. [PMID: 38656446 DOI: 10.1063/5.0191857] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2023] [Accepted: 04/07/2024] [Indexed: 04/26/2024] Open
Abstract
This work describes optimizations of N-electron system wave functions by means of the simulated annealing technique within the doubly occupied configuration interaction framework. Using that technique, we minimize the energy variance of a Hamiltonian, providing determinations of wave functions corresponding to ground or excited states in an identical manner. The procedure that allows us to determine electronic spectra can be performed using treatments of restricted or unrestricted types. The results found in selected systems, described in terms of energy, spin, and wave function, are analyzed, showing the performance of each method. We also compare these results with those arising from more traditional approaches that minimize the energy, in both restricted and unrestricted versions, and with those obtained from the full configuration interaction treatment.
Collapse
Affiliation(s)
- Diego R Alcoba
- Universidad de Buenos Aires, Facultad de Ciencias Exactas y Naturales, Departamento de Física, Ciudad Universitaria, 1428 Buenos Aires, Argentina
- CONICET - Universidad de Buenos Aires, Instituto de Física de Buenos Aires (IFIBA), Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Ofelia B Oña
- Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas, Universidad Nacional de La Plata, Consejo Nacional de Investigaciones Cientıficas y Tócnicas, Diag. 113 y 64 (S/N), Sucursal 4, CC 16, 1900 La Plata, Argentina
| | - Alicia Torre
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco, Apdo. 644, E-48080 Bilbao, Spain
| | - Luis Lain
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco, Apdo. 644, E-48080 Bilbao, Spain
| | - Guadalupe Sierra
- Universidad de Buenos Aires, Facultad de Ciencias Exactas y Naturales, Departamento de Física, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Gustavo E Massaccesi
- Departamento de Ciencias Exactas, Ciclo Básico Común, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
- Instituto de Investigaciones Matemáticas Luis Santaló, Consejo Nacional de Investigaciones Científicas y Técnicas, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| |
Collapse
|
2
|
Garros A, Alcoba DR, Capuzzi P, Lain L, Torre A, Oña OB, Dukelsky J. Determination of electronic excitation energies within the doubly occupied configuration interaction space by means of the Hermitian operator method. J Chem Phys 2023; 159:124107. [PMID: 38127377 DOI: 10.1063/5.0168585] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2023] [Accepted: 09/06/2023] [Indexed: 12/23/2023] Open
Abstract
In this work, we formulate the equations of motion corresponding to the Hermitian operator method in the framework of the doubly occupied configuration interaction space. The resulting algorithms turn out to be considerably simpler than the equations provided by that method in more conventional spaces, enabling the determination of excitation energies in N-electron systems under an affordable polynomial computational cost. The implementation of this technique only requires to know the elements of low-order reduced density matrices of an N-electron reference state, which can be obtained from any approximate method. We contrast our procedure against the reduced Bardeen-Cooper-Schrieffer and Richardson-Gaudin-Kitaev integrable models, pointing out the reliability of our proposal.
Collapse
Affiliation(s)
- Adán Garros
- Universidad de Buenos Aires, Facultad de Ciencias Exactas y Naturales, Departamento de Física, Ciudad Universitaria, 1428 Buenos Aires, Argentina
- CONICET-Universidad de Buenos Aires, Instituto de Física de Buenos Aires (IFIBA), Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Diego R Alcoba
- Universidad de Buenos Aires, Facultad de Ciencias Exactas y Naturales, Departamento de Física, Ciudad Universitaria, 1428 Buenos Aires, Argentina
- CONICET-Universidad de Buenos Aires, Instituto de Física de Buenos Aires (IFIBA), Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Pablo Capuzzi
- Universidad de Buenos Aires, Facultad de Ciencias Exactas y Naturales, Departamento de Física, Ciudad Universitaria, 1428 Buenos Aires, Argentina
- CONICET-Universidad de Buenos Aires, Instituto de Física de Buenos Aires (IFIBA), Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Luis Lain
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco, Apdo. 644, E-48080 Bilbao, Spain
| | - Alicia Torre
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco, Apdo. 644, E-48080 Bilbao, Spain
| | - Ofelia B Oña
- Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas, Universidad Nacional de La Plata, Consejo Nacional de Investigaciones Científicas y Técnicas, Diag. 113 y 64 (S/N), Sucursal 4, CC 16, 1900 La Plata, Argentina
| | - Jorge Dukelsky
- Instituto de Estructura de la Materia, CSIC, Serrano 123, E-28006 Madrid, Spain
| |
Collapse
|
3
|
Jones GM, Li RR, DePrince AE, Vogiatzis KD. Data-Driven Refinement of Electronic Energies from Two-Electron Reduced-Density-Matrix Theory. J Phys Chem Lett 2023:6377-6385. [PMID: 37418691 DOI: 10.1021/acs.jpclett.3c01382] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/09/2023]
Abstract
The exponential computational cost of describing strongly correlated electrons can be mitigated by adopting a reduced-density matrix (RDM)-based description of the electronic structure. While variational two-electron RDM (v2RDM) methods can enable large-scale calculations on such systems, the quality of the solution is limited by the fact that only a subset of known necessary N-representability constraints can be applied to the 2RDM in practical calculations. Here, we demonstrate that violations of partial three-particle (T1 and T2) N-representability conditions, which can be evaluated with knowledge of only the 2RDM, can serve as physics-based features in a machine-learning (ML) protocol for improving energies from v2RDM calculations that consider only two-particle (PQG) conditions. Proof-of-principle calculations demonstrate that the model yields substantially improved energies relative to reference values from configuration-interaction-based calculations.
Collapse
Affiliation(s)
- Grier M Jones
- Department of Chemistry, University of Tennessee, Knoxville, Tennessee 37996, United States
| | - Run R Li
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | - A Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | | |
Collapse
|
4
|
Mazziotti DA. Quantum Many-Body Theory from a Solution of the N-Representability Problem. PHYSICAL REVIEW LETTERS 2023; 130:153001. [PMID: 37115895 DOI: 10.1103/physrevlett.130.153001] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/15/2022] [Revised: 02/13/2023] [Accepted: 03/16/2023] [Indexed: 06/19/2023]
Abstract
Here we present a many-body theory based on a solution of the N-representability problem in which the ground-state two-particle reduced density matrix (2-RDM) is determined directly without the many-particle wave function. We derive an equation that re-expresses physical constraints on higher-order RDMs to generate direct constraints on the 2-RDM, which are required for its derivation from an N-particle density matrix, known as N-representability conditions. The approach produces a complete hierarchy of 2-RDM constraints that do not depend explicitly upon the higher RDMs or the wave function. By using the two-particle part of a unitary decomposition of higher-order constraint matrices, we can solve the energy minimization by semidefinite programming in a form where the low-rank structure of these matrices can be potentially exploited. We illustrate by computing the ground-state electronic energy and properties of the H_{8} ring.
Collapse
Affiliation(s)
- David A Mazziotti
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
5
|
Sager-Smith LM, Mazziotti DA. Reducing the Quantum Many-Electron Problem to Two Electrons with Machine Learning. J Am Chem Soc 2022; 144:18959-18966. [PMID: 36194786 DOI: 10.1021/jacs.2c07112] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
An outstanding challenge in chemical computation is the many-electron problem where computational methodologies scale prohibitively with system size. The energy of any molecule can be expressed as a weighted sum of the energies of two-electron wave functions that are computable from only a two-electron calculation. Despite the physical elegance of this extended "aufbau" principle, the determination of the distribution of weights─geminal occupations─for general molecular systems has remained elusive. Here we introduce a new paradigm for electronic structure where approximate geminal-occupation distributions are "learned" via a convolutional neural network. We show that the neural network learns the N-representability conditions, constraints on the distribution for it to represent an N-electron system. By training on hydrocarbon isomers with only 2-7 carbon atoms, we are able to predict the energies for isomers of octane as well as hydrocarbons with 8-15 carbons. The present work demonstrates that machine learning can be used to reduce the many-electron problem to an effective two-electron problem, opening new opportunities for accurately predicting electronic structure.
Collapse
Affiliation(s)
- LeeAnn M Sager-Smith
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois60637, United States
| | - David A Mazziotti
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois60637, United States
| |
Collapse
|
6
|
Li RR, Rubin NC, DePrince AE. Challenges for Variational Reduced-Density-Matrix Theory: Total Angular Momentum Constraints. J Chem Theory Comput 2022; 18:5966-5977. [PMID: 36150206 DOI: 10.1021/acs.jctc.2c00402] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The variational two-electron reduced density matrix (v2RDM) method is generalized for the description of total angular momentum (J) and projection of total angular momentum (MJ) states in atomic systems described by nonrelativistic Hamiltonians, and it is shown that the approach exhibits serious deficiencies. Under ensemble N-representability constraints, v2RDM theory fails to retain the appropriate degeneracies among various J states for fixed spin (S) and orbital angular momentum (L), and for fixed L, S, and J, the manifold of MJ states is not necessarily degenerate. Moreover, a substantial energy error is observed for a system for which the two-electron reduced density matrix is exactly ensemble N-representable; in this case, the error stems from violations in pure-state N-representability conditions. Unfortunately, such violations do not appear to be good indicators of the reliability of energies from v2RDM theory in general. Several states are identified for which energy errors are near zero and yet pure-state conditions are clearly violated.
Collapse
Affiliation(s)
- Run R Li
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | | | - A Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| |
Collapse
|
7
|
Ríos E, Alcoba DR, Lain L, Torre A, Oña OB, Massaccesi GE. Variational determination of the two-electron reduced density matrix within the doubly occupied configuration interaction framework: Treatments of triplet N-electron systems. J Chem Phys 2022; 157:014102. [DOI: 10.1063/5.0088125] [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
This work performs variational determinations of two-electron reduced density matrices corresponding to open-shell N-electron systems within the framework of the doubly occupied configuration interaction treatment, traditionally limited to studies of closed-shell systems. The procedure has allowed us to describe satisfactorily molecular systems in triplet states following two methods. One of them adds hydrogen atoms at an infinite distance of the triplet system studied, constituting a singlet supersystem. Energies and reduced density matrices of the triplet system are obtained by removing the contributions of the added atoms from the singlet supersystem results. The second procedure determines variationally the two-electron reduced density matrices corresponding to the triplet systems by means of adequate couplings of basis-set functions. Both models have been managed by imposing N-representability conditions on the reduced density matrix calculations. Results obtained from these methods for molecular systems in triplet ground states are reported and compared with those provided by benchmark methods.
Collapse
Affiliation(s)
| | - Diego Ricardo Alcoba
- Departmento de Fisica, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Argentina
| | - Luis Lain
- Quimica Fisica, Universidad del Pais Vasco Facultad de Ciencia y Tecnologia, Spain
| | - Alicia Torre
- Quimica Fisica, Universidad del Pais Vasco, Spain
| | | | | |
Collapse
|
8
|
Elayan IA, Gupta R, Hollett JW. ΔNO and the complexities of electron correlation in simple hydrogen clusters. J Chem Phys 2022; 156:094102. [DOI: 10.1063/5.0073227] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
9
|
Li RR, Liebenthal MD, DePrince AE. Challenges for variational reduced-density-matrix theory with three-particle N-representability conditions. J Chem Phys 2021; 155:174110. [PMID: 34742213 DOI: 10.1063/5.0066404] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The direct variational optimization of the two-electron reduced density matrix (2RDM) can provide a reference-independent description of the electronic structure of many-electron systems that naturally capture strong or nondynamic correlation effects. Such variational 2RDM approaches can often provide a highly accurate description of strong electron correlation, provided that the 2RDMs satisfy at least partial three-particle N-representability conditions (e.g., the T2 condition). However, recent benchmark calculations on hydrogen clusters [N. H. Stair and F. A. Evangelista, J. Chem. Phys. 153, 104108 (2020)] suggest that even the T2 condition leads to unacceptably inaccurate results in the case of two- and three-dimensional clusters. We demonstrate that these failures persist under the application of full three-particle N-representability conditions (3POS). A variety of correlation metrics are explored in order to identify regimes under which 3POS calculations become unreliable, and we find that the relative squared magnitudes of the cumulant three- and two-particle reduced density matrices correlate reasonably well with the energy error in these systems. However, calculations on other molecular systems reveal that this metric is not a universal indicator for the reliability of the reduced-density-matrix theory with 3POS conditions.
Collapse
Affiliation(s)
- Run R Li
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, USA
| | - Marcus D Liebenthal
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, USA
| | - A Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, USA
| |
Collapse
|
10
|
Affiliation(s)
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, FL, USA
| |
Collapse
|
11
|
Alcoba DR, Oña OB, Lain L, Torre A, Capuzzi P, Massaccesi GE, Ríos E, Rubio-García A, Dukelsky J. Variational determination of ground and excited-state two-electron reduced density matrices in the doubly occupied configuration space: A dispersion operator approach. J Chem Phys 2021; 154:224104. [PMID: 34241224 DOI: 10.1063/5.0051793] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This work implements a variational determination of the elements of two-electron reduced density matrices corresponding to the ground and excited states of N-electron interacting systems based on the dispersion operator technique. The procedure extends the previously reported proposal [Nakata et al., J. Chem. Phys. 125, 244109 (2006)] to two-particle interaction Hamiltonians and N-representability conditions for the two-, three-, and four-particle reduced density matrices in the doubly occupied configuration interaction space. The treatment has been applied to describe electronic spectra using two benchmark exactly solvable pairing models: reduced Bardeen-Cooper-Schrieffer and Richardson-Gaudin-Kitaev Hamiltonians. The dispersion operator combined with N-representability conditions up to the four-particle reduced density matrices provides excellent results.
Collapse
Affiliation(s)
- Diego R Alcoba
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Ofelia B Oña
- Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas, Universidad de la Plata, Consejo Nacional de Investigaciones Científicas y Técnicas, Diag. 113 y 64 (S/N), Sucursal 4, CC 16, 1900 La Plata, Argentina
| | - Luis Lain
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco, Apdo. 644, E-48080 Bilbao, Spain
| | - Alicia Torre
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco, Apdo. 644, E-48080 Bilbao, Spain
| | - Pablo Capuzzi
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Gustavo E Massaccesi
- Departamento de Ciencias Exactas, Ciclo Básico Común, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Elías Ríos
- Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas, Universidad de la Plata, Consejo Nacional de Investigaciones Científicas y Técnicas, Diag. 113 y 64 (S/N), Sucursal 4, CC 16, 1900 La Plata, Argentina
| | - Alvaro Rubio-García
- Instituto de Estructura de la Materia, CSIC, Serrano 123, 28006 Madrid, Spain
| | - Jorge Dukelsky
- Instituto de Estructura de la Materia, CSIC, Serrano 123, 28006 Madrid, Spain
| |
Collapse
|
12
|
Head-Marsden K, Flick J, Ciccarino CJ, Narang P. Quantum Information and Algorithms for Correlated Quantum Matter. Chem Rev 2020; 121:3061-3120. [PMID: 33326218 DOI: 10.1021/acs.chemrev.0c00620] [Citation(s) in RCA: 37] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Discoveries in quantum materials, which are characterized by the strongly quantum-mechanical nature of electrons and atoms, have revealed exotic properties that arise from correlations. It is the promise of quantum materials for quantum information science superimposed with the potential of new computational quantum algorithms to discover new quantum materials that inspires this Review. We anticipate that quantum materials to be discovered and developed in the next years will transform the areas of quantum information processing including communication, storage, and computing. Simultaneously, efforts toward developing new quantum algorithmic approaches for quantum simulation and advanced calculation methods for many-body quantum systems enable major advances toward functional quantum materials and their deployment. The advent of quantum computing brings new possibilities for eliminating the exponential complexity that has stymied simulation of correlated quantum systems on high-performance classical computers. Here, we review new algorithms and computational approaches to predict and understand the behavior of correlated quantum matter. The strongly interdisciplinary nature of the topics covered necessitates a common language to integrate ideas from these fields. We aim to provide this common language while weaving together fields across electronic structure theory, quantum electrodynamics, algorithm design, and open quantum systems. Our Review is timely in presenting the state-of-the-art in the field toward algorithms with nonexponential complexity for correlated quantum matter with applications in grand-challenge problems. Looking to the future, at the intersection of quantum information science and algorithms for correlated quantum matter, we envision seminal advances in predicting many-body quantum states and describing excitonic quantum matter and large-scale entangled states, a better understanding of high-temperature superconductivity, and quantifying open quantum system dynamics.
Collapse
Affiliation(s)
- Kade Head-Marsden
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Johannes Flick
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
| | - Christopher J Ciccarino
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States.,Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Prineha Narang
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| |
Collapse
|
13
|
Oña OB, Torre A, Lain L, Alcoba DR, Ríos E, Massaccesi GE. Variational determination of the two-electron reduced density matrix within the doubly occupied configuration interaction scheme: An extension to the study of open-shell systems. J Chem Phys 2020; 153:084101. [DOI: 10.1063/5.0020581] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ofelia B. Oña
- Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas, Universidad Nacional de La Plata, CCT La Plata, Consejo Nacional de Investigaciones Científicas y Técnicas, Diag. 113 y 64 (S/N), Sucursal 4, CC 16, 1900 La Plata, Argentina
| | - Alicia Torre
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco, Apdo. 644, E-48080 Bilbao, Spain
| | - Luis Lain
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco, Apdo. 644, E-48080 Bilbao, Spain
| | - Diego R. Alcoba
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
- Instituto de Física de Buenos Aires, Consejo Nacional de Investigaciones Científicas y Técnicas, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Elías Ríos
- Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas, Universidad Nacional de La Plata, CCT La Plata, Consejo Nacional de Investigaciones Científicas y Técnicas, Diag. 113 y 64 (S/N), Sucursal 4, CC 16, 1900 La Plata, Argentina
| | - Gustavo E. Massaccesi
- Departamento de Ciencias Exactas, Ciclo Básico Común, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| |
Collapse
|
14
|
Maradzike E, Hapka M, Pernal K, DePrince AE. Reduced Density Matrix-Driven Complete Active Apace Self-Consistent Field Corrected for Dynamic Correlation from the Adiabatic Connection. J Chem Theory Comput 2020; 16:4351-4360. [PMID: 32538086 DOI: 10.1021/acs.jctc.0c00324] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
The recently proposed multireference adiabatic connection (AC) formalism [Pernal, Phys. Rev. Lett. 120, 013001 (2018)] is applied to recover dynamic electron correlation effects lacking in variational two-electron reduced density matrix (v2RDM)-driven complete active space self-consistent field theory (CASSCF). The AC approach is validated by computing potential energy curves for the dissociation of molecular nitrogen and the symmetric double dissociation of H2O while enforcing two sets of approximate N-representability conditions in the underlying v2RDM-driven CASSCF calculations (either two-particle or two-particle plus partial three-particle conditions). The AC yields smaller absolute errors than second-order N-electron perturbation theory (NEVPT2) at all molecular geometries for both sets of the N-representability conditions considered. The efficacy of the approach for thermochemistry is also assessed for a set of 31 small-molecule reactions. When imposing partial three-particle N-representability conditions, mean and maximum unsigned errors in reaction energies from the AC are superior to those from NEVPT2.
Collapse
Affiliation(s)
- Elvis Maradzike
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | - 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
| | - A Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| |
Collapse
|
15
|
Vu N, DePrince AE. Size-extensive seniority-zero energy functionals derived from configuration interaction with double excitations. J Chem Phys 2020; 152:244103. [DOI: 10.1063/5.0011637] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023] Open
Affiliation(s)
- Nam Vu
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, USA
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, USA
| |
Collapse
|
16
|
Head-Marsden K, Mazziotti DA. Active-Space Pair Two-Electron Reduced Density Matrix Theory for Strong Correlation. J Phys Chem A 2020; 124:4848-4854. [DOI: 10.1021/acs.jpca.0c01937] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Kade Head-Marsden
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637, United States
| | - David A. Mazziotti
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
17
|
Hollett JW, Loos PF. Capturing static and dynamic correlation with ΔNO-MP2 and ΔNO-CCSD. J Chem Phys 2020; 152:014101. [PMID: 31914756 DOI: 10.1063/1.5140669] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The ΔNO method for static correlation is combined with second-order Møller-Plesset perturbation theory (MP2) and coupled-cluster singles and doubles (CCSD) to account for dynamic correlation. The MP2 and CCSD expressions are adapted from finite-temperature CCSD, which includes orbital occupancies and vacancies, and expanded orbital summations. Correlation is partitioned with the aid of damping factors incorporated into the MP2 and CCSD residual equations. Potential energy curves for a selection of diatomics are in good agreement with extrapolated full configuration interaction results and on par with conventional multireference approaches.
Collapse
Affiliation(s)
- Joshua W Hollett
- Department of Chemistry, University of Winnipeg, Winnipeg, Manitoba R3B 2G3, Canada
| | - Pierre-François Loos
- Laboratoire de Chimie et Physique Quantiques (UMR 5626), Université de Toulouse, CNRS, UPS, Toulouse, France
| |
Collapse
|
18
|
Vu N, Mitxelena I, DePrince AE. An adiabatic connection for doubly-occupied configuration interaction wave functions. J Chem Phys 2019; 151:244121. [PMID: 31893894 DOI: 10.1063/1.5130660] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Abstract
An adiabatic connection (AC) is developed as an electron correlation correction for doubly occupied configuration interaction (DOCI) wave functions. Following the work of Pernal [Phys. Rev. Lett. 120, 013001 (2018)], the working equations of the approach, termed AC-DOCI, are rooted in the extended random phase approximation (ERPA) and require knowledge of only the ground-state two-electron reduced density matrix (2RDM) from the DOCI. As such, the AC is naturally suited to modeling electron correlation in variational 2RDM (v2RDM)-based approximations to the DOCI. The v2RDM-driven AC-DOCI is applied to the dissociation of molecular nitrogen and the double dissociation of water; the approach yields energies that are similar in quality to those from second-order multireference perturbation theory near equilibrium, but the quality of the AC-DOCI energy degrades at stretched geometries. The exact adiabatic connection path suggests the assumption that the one-electron reduced-density matrix is constant along the AC path is invalid at stretched geometries, but this deficiency alone cannot explain the observed behavior. Rather, it appears that the ERPA's single-particle-transition ansatz cannot, in general, provide good approximations to the 2RDM along the AC path. The AC-DOCI is also applied to a set of 45 reaction energies; for these systems, the approach has an average accuracy that is comparable to that of single-reference second-order many-body perturbation theory.
Collapse
Affiliation(s)
- Nam Vu
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, USA
| | - Ion Mitxelena
- Kimika Fakultatea, Euskal Herriko Unibertsitatea (UPV/EHU), P.K. 1072, 20080 Donostia, Euskadi, Spain
| | - A Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, USA
| |
Collapse
|
19
|
Rubio-García A, Dukelsky J, Alcoba DR, Capuzzi P, Oña OB, Ríos E, Torre A, Lain L. Variational reduced density matrix method in the doubly-occupied configuration interaction space using four-particle N-representability conditions: Application to the XXZ model of quantum magnetism. J Chem Phys 2019; 151:154104. [PMID: 31640387 DOI: 10.1063/1.5118899] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This work deals with the variational determination of the two-particle reduced density matrix (2-RDM) and the energy corresponding to the ground state of N-particle systems within the doubly occupied configuration interaction (DOCI) space. Here, we impose for the first time up to four-particle N-representability constraint conditions in the variational determination of the 2-RDM matrix elements using the standard semidefinite programming algorithms. The energies and 2-RDMs obtained from this treatment and the corresponding computational costs are compared with those arisen from previously reported less restrictive variational methods [D. R. Alcoba et al., J. Chem. Phys. 149, 194105 (2018)] as well as with the exact DOCI values. We apply the different approximations to the one-dimensional XXZ model of quantum magnetism, which has a rich phase diagram with one critical phase and constitutes a stringent test for the method. The numerical results show the usefulness of our treatment to achieve a high degree of accuracy.
Collapse
Affiliation(s)
- A Rubio-García
- Instituto de Estructura de la Materia, CSIC, Serrano 123, 28006 Madrid, Spain
| | - J Dukelsky
- Instituto de Estructura de la Materia, CSIC, Serrano 123, 28006 Madrid, Spain
| | - D R Alcoba
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - P Capuzzi
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - O B Oña
- Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas, Universidad Nacional de La Plata, CCT La Plata, Consejo Nacional de Investigaciones Científicas y Técnicas, Diag. 113 y 64 (S/N), Sucursal 4, CC 16, 1900 La Plata, Argentina
| | - E Ríos
- Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas, Universidad Nacional de La Plata, CCT La Plata, Consejo Nacional de Investigaciones Científicas y Técnicas, Diag. 113 y 64 (S/N), Sucursal 4, CC 16, 1900 La Plata, Argentina
| | - A Torre
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco, Apdo. 644, E-48080 Bilbao, Spain
| | - L Lain
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco, Apdo. 644, E-48080 Bilbao, Spain
| |
Collapse
|
20
|
Mullinax JW, Maradzike E, Koulias LN, Mostafanejad M, Epifanovsky E, Gidofalvi G, DePrince AE. Heterogeneous CPU + GPU Algorithm for Variational Two-Electron Reduced-Density Matrix-Driven Complete Active-Space Self-Consistent Field Theory. J Chem Theory Comput 2019; 15:6164-6178. [DOI: 10.1021/acs.jctc.9b00768] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- J. Wayne Mullinax
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, United States
| | - Elvis Maradzike
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, United States
| | - Lauren N. Koulias
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, United States
| | - Mohammad Mostafanejad
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, United States
| | - Evgeny Epifanovsky
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, United States
| | - Gergely Gidofalvi
- Department of Chemistry and Biochemistry, Gonzaga University, Spokane, Washington 99258, United States
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, United States
| |
Collapse
|
21
|
Alcoba DR, Torre A, Lain L, Massaccesi GE, Oña OB, Ríos E. Unrestricted treatment for the direct variational determination of the two-electron reduced density matrix for doubly occupied-configuration-interaction wave functions. J Chem Phys 2019; 150:164106. [PMID: 31042927 DOI: 10.1063/1.5092182] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This work extends to the unrestricted orbital approach the procedure described in our previous report [Alcoba et al., J. Chem. Phys. 148, 024105 (2018)] for determining variationally the two-electron reduced density matrix arising from doubly occupied-configuration-interaction wave functions by imposing two- and three-index N-representability conditions. An analysis of the numerical results obtained in selected systems, from both restricted and unrestricted treatments, allows one to assess the performance of these methodologies as well as to show the influence of the P, Q, G, T1, and T2 positivity conditions. We highlight the satisfactory results obtained within the unrestricted scheme.
Collapse
Affiliation(s)
- Diego R Alcoba
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires and Instituto de Física de Buenos Aires, Consejo Nacional de Investigaciones Científicas y Técnicas. Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Alicia Torre
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco. Apdo. 644, E-48080 Bilbao, Spain
| | - Luis Lain
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco. Apdo. 644, E-48080 Bilbao, Spain
| | - Gustavo E Massaccesi
- Departamento de Ciencias Exactas, Ciclo Básico Común, Universidad de Buenos Aires. Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Ofelia B Oña
- Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas, Universidad Nacional de La Plata, CCT La Plata, Consejo Nacional de Investigaciones Científicas y Técnicas. Diag. 113 y 64 (S/N), Sucursal 4, CC 16, 1900 La Plata, Argentina
| | - Elías Ríos
- Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas, Universidad Nacional de La Plata, CCT La Plata, Consejo Nacional de Investigaciones Científicas y Técnicas. Diag. 113 y 64 (S/N), Sucursal 4, CC 16, 1900 La Plata, Argentina
| |
Collapse
|
22
|
Alcoba DR, Capuzzi P, Rubio-García A, Dukelsky J, Massaccesi GE, Oña OB, Torre A, Lain L. Variational reduced density matrix method in the doubly occupied configuration interaction space using three-particleN-representability conditions. J Chem Phys 2018; 149:194105. [DOI: 10.1063/1.5056247] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Diego R. Alcoba
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
- Instituto de Física de Buenos Aires, Consejo Nacional de Investigaciones Científicas y Técnicas, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Pablo Capuzzi
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
- Instituto de Física de Buenos Aires, Consejo Nacional de Investigaciones Científicas y Técnicas, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Alvaro Rubio-García
- Instituto de Estructura de la Materia, CSIC, Serrano 123, 28006 Madrid, Spain
| | - Jorge Dukelsky
- Instituto de Estructura de la Materia, CSIC, Serrano 123, 28006 Madrid, Spain
| | - Gustavo E. Massaccesi
- Departamento de Ciencias Exactas, Ciclo Básico Común, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Ofelia B. Oña
- Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas, Universidad Nacional de La Plata, CCT La Plata, Consejo Nacional de Investigaciones Científicas y Técnicas, Diag. 113 y 64 (S/N), Sucursal 4, CC 16, 1900 La Plata, Argentina
| | - Alicia Torre
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco, Apdo. 644, E-48080 Bilbao, Spain
| | - Luis Lain
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco, Apdo. 644, E-48080 Bilbao, Spain
| |
Collapse
|
23
|
Bytautas L, Dukelsky J. Seniority based energy renormalization group (Ω-ERG) approach in quantum chemistry: Initial formulation and application to potential energy surfaces. COMPUT THEOR CHEM 2018. [DOI: 10.1016/j.comptc.2018.08.011] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
24
|
Rubio-García A, Alcoba DR, Capuzzi P, Dukelsky J. Benchmarking the Variational Reduced Density Matrix Theory in the Doubly Occupied Configuration Interaction Space with Integrable Pairing Models. J Chem Theory Comput 2018; 14:4183-4192. [PMID: 29906104 DOI: 10.1021/acs.jctc.8b00387] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The variational reduced density matrix theory has been recently applied with great success to models within the truncated doubly occupied configuration interaction space, which corresponds to the seniority zero subspace. Conservation of the seniority quantum number restricts the Hamiltonians to be based on the SU(2) algebra. Among them there is a whole family of exactly solvable Richardson-Gaudin pairing Hamiltonians. We benchmark the variational theory against two different exactly solvable models, the Richardson-Gaudin-Kitaev and the reduced BCS Hamiltonians. We obtain exact numerical results for the so-called [Formula: see text] N-representability conditions in both cases for systems that go from 10 to 100 particles. However, when random single-particle energies as appropriate for small superconducting grains are considered, the exactness is lost but still a high accuracy is obtained.
Collapse
Affiliation(s)
- A Rubio-García
- Instituto de Estructura de la Materia, CSIC, Serrano 123 , 28006 Madrid , Spain
| | - D R Alcoba
- Departamento de Física, Facultad de Ciencias Exactas y Naturales , Universidad de Buenos Aires, Ciudad Universitaria , 1428 Buenos Aires , Argentina.,Instituto de Física de Buenos Aires, Consejo Nacional de Investigaciones Científicas y Técnicas , Ciudad Universitaria , 1428 Buenos Aires , Argentina
| | - P Capuzzi
- Departamento de Física, Facultad de Ciencias Exactas y Naturales , Universidad de Buenos Aires, Ciudad Universitaria , 1428 Buenos Aires , Argentina.,Instituto de Física de Buenos Aires, Consejo Nacional de Investigaciones Científicas y Técnicas , Ciudad Universitaria , 1428 Buenos Aires , Argentina
| | - J Dukelsky
- Instituto de Estructura de la Materia, CSIC, Serrano 123 , 28006 Madrid , Spain
| |
Collapse
|
25
|
Lanssens C, Ayers PW, Van Neck D, De Baerdemacker S, Gunst K, Bultinck P. Method for making 2-electron response reduced density matrices approximately N-representable. J Chem Phys 2018; 148:084104. [PMID: 29495792 DOI: 10.1063/1.4994618] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In methods like geminal-based approaches or coupled cluster that are solved using the projected Schrödinger equation, direct computation of the 2-electron reduced density matrix (2-RDM) is impractical and one falls back to a 2-RDM based on response theory. However, the 2-RDMs from response theory are not N-representable. That is, the response 2-RDM does not correspond to an actual physical N-electron wave function. We present a new algorithm for making these non-N-representable 2-RDMs approximately N-representable, i.e., it has the right symmetry and normalization and it fulfills the P-, Q-, and G-conditions. Next to an algorithm which can be applied to any 2-RDM, we have also developed a 2-RDM optimization procedure specifically for seniority-zero 2-RDMs. We aim to find the 2-RDM with the right properties which is the closest (in the sense of the Frobenius norm) to the non-N-representable 2-RDM by minimizing the square norm of the difference between this initial response 2-RDM and the targeted 2-RDM under the constraint that the trace is normalized and the 2-RDM, Q-matrix, and G-matrix are positive semidefinite, i.e., their eigenvalues are non-negative. Our method is suitable for fixing non-N-representable 2-RDMs which are close to being N-representable. Through the N-representability optimization algorithm we add a small correction to the initial 2-RDM such that it fulfills the most important N-representability conditions.
Collapse
Affiliation(s)
- Caitlin Lanssens
- Center for Molecular Modeling, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| | - Paul W Ayers
- Department of Chemistry and Chemical Biology, McMaster University, 1280 Main Street West, Hamilton, Ontario L8S 4L8, Canada
| | - Dimitri Van Neck
- Center for Molecular Modeling, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| | - Stijn De Baerdemacker
- Center for Molecular Modeling, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| | - Klaas Gunst
- Center for Molecular Modeling, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| | - Patrick Bultinck
- Department of Inorganic and Physical Chemistry, Ghent University, Krijgslaan 281 (S3), 9000 Ghent, Belgium
| |
Collapse
|
26
|
Alcoba DR, Torre A, Lain L, Massaccesi GE, Oña OB, Honoré EM, Poelmans W, Van Neck D, Bultinck P, De Baerdemacker S. Direct variational determination of the two-electron reduced density matrix for doubly occupied-configuration-interaction wave functions: The influence of three-index N-representability conditions. J Chem Phys 2018; 148:024105. [PMID: 29331147 DOI: 10.1063/1.5008811] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This work proposes the variational determination of two-electron reduced density matrices corresponding to the ground state of N-electron systems within the doubly occupied-configuration-interaction methodology. The P, Q, and G two-index N-representability conditions have been extended to the T1 and T2 (T2') three-index ones and the resulting optimization problem has been addressed using a standard semidefinite program. We report results obtained from the doubly occupied-configuration-interaction method, from the two-index constraint variational procedure and from the two- and three-index constraint variational treatment. The discussion of these results along with a study of the computational cost demanded shows the usefulness of our proposal.
Collapse
Affiliation(s)
- Diego R Alcoba
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Alicia Torre
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco. Apdo. 644, E-48080 Bilbao, Spain
| | - Luis Lain
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco. Apdo. 644, E-48080 Bilbao, Spain
| | - Gustavo E Massaccesi
- Departamento de Ciencias Exactas, Ciclo Básico Común, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Ofelia B Oña
- Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas, Universidad Nacional de La Plata, CCT La Plata, Consejo Nacional de Investigaciones Científicas y Técnicas. Diag. 113 y 64 (S/N), Sucursal 4, CC 16, 1900 La Plata, Argentina
| | - Eduardo M Honoré
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Ward Poelmans
- Center for Molecular Modeling, Ghent University, Technologiepark 903, B-9052 Zwijnaarde, Belgium
| | - Dimitri Van Neck
- Center for Molecular Modeling, Ghent University, Technologiepark 903, B-9052 Zwijnaarde, Belgium
| | - Patrick Bultinck
- Department of Inorganic and Physical Chemistry, Ghent University, Krijgslaan 281 (S3), 9000 Gent, Belgium
| | - Stijn De Baerdemacker
- Center for Molecular Modeling, Ghent University, Technologiepark 903, B-9052 Zwijnaarde, Belgium
| |
Collapse
|
27
|
Hybrid Treatments Based on Determinant Seniority Numbers and Spatial Excitation Levels in the Configuration Interaction Framework. ADVANCES IN QUANTUM CHEMISTRY 2018. [DOI: 10.1016/bs.aiq.2017.05.003] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
28
|
Vogiatzis KD, Ma D, Olsen J, Gagliardi L, de Jong WA. Pushing configuration-interaction to the limit: Towards massively parallel MCSCF calculations. J Chem Phys 2017; 147:184111. [DOI: 10.1063/1.4989858] [Citation(s) in RCA: 88] [Impact Index Per Article: 12.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
Affiliation(s)
- Konstantinos D. Vogiatzis
- Department of Chemistry, Minnesota Supercomputing Institute, and Chemical Theory Center, University of Minnesota, 207 Pleasant Street Southeast, Minneapolis, Minnesota 55455-0431, USA
| | - Dongxia Ma
- Department of Chemistry, Minnesota Supercomputing Institute, and Chemical Theory Center, University of Minnesota, 207 Pleasant Street Southeast, Minneapolis, Minnesota 55455-0431, USA
| | - Jeppe Olsen
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Laura Gagliardi
- Department of Chemistry, Minnesota Supercomputing Institute, and Chemical Theory Center, University of Minnesota, 207 Pleasant Street Southeast, Minneapolis, Minnesota 55455-0431, USA
| | - Wibe A. de Jong
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
29
|
Johnson PA, Limacher PA, Kim TD, Richer M, Miranda-Quintana RA, Heidar-Zadeh F, Ayers PW, Bultinck P, De Baerdemacker S, Van Neck D. Strategies for extending geminal-based wavefunctions: Open shells and beyond. COMPUT THEOR CHEM 2017. [DOI: 10.1016/j.comptc.2017.05.010] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
30
|
Head-Marsden K, Mazziotti DA. Pair 2-electron reduced density matrix theory using localized orbitals. J Chem Phys 2017; 147:084101. [PMID: 28863513 DOI: 10.1063/1.4999423] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Full configuration interaction (FCI) restricted to a pairing space yields size-extensive correlation energies but its cost scales exponentially with molecular size. Restricting the variational two-electron reduced-density-matrix (2-RDM) method to represent the same pairing space yields an accurate lower bound to the pair FCI energy at a mean-field-like computational scaling of O(r3) where r is the number of orbitals. In this paper, we show that localized molecular orbitals can be employed to generate an efficient, approximately size-extensive pair 2-RDM method. The use of localized orbitals eliminates the substantial cost of optimizing iteratively the orbitals defining the pairing space without compromising accuracy. In contrast to the localized orbitals, the use of canonical Hartree-Fock molecular orbitals is shown to be both inaccurate and non-size-extensive. The pair 2-RDM has the flexibility to describe the spectra of one-electron RDM occupation numbers from all quantum states that are invariant to time-reversal symmetry. Applications are made to hydrogen chains and their dissociation, n-acene from naphthalene through octacene, and cadmium telluride 2-, 3-, and 4-unit polymers. For the hydrogen chains, the pair 2-RDM method recovers the majority of the energy obtained from similar calculations that iteratively optimize the orbitals. The localized-orbital pair 2-RDM method with its mean-field-like computational scaling and its ability to describe multi-reference correlation has important applications to a range of strongly correlated phenomena in chemistry and physics.
Collapse
Affiliation(s)
- Kade Head-Marsden
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| | - David A Mazziotti
- Department of Chemistry and The James Franck Institute, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
31
|
Affiliation(s)
- Peter A. Limacher
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany
| |
Collapse
|
32
|
Alcoba DR, Torre A, Lain L, Massaccesi GE, Oña OB, Capuzzi P. Spin contamination-free N-electron wave functions in the excitation-based configuration interaction treatment. J Chem Phys 2016; 145:014109. [PMID: 27394101 DOI: 10.1063/1.4955159] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
This work deals with the spin contamination in N-electron wave functions provided by the excitation-based configuration interaction methods. We propose a procedure to ensure a suitable selection of excited N-electron Slater determinants with respect to a given reference determinant, required in these schemes. The procedure guarantees the construction of N-electron wave functions which are eigenfunctions of the spin-squared operator Sˆ(2), avoiding any spin contamination. Our treatment is based on the evaluation of the excitation level of the determinants by means of the expectation value of an excitation operator formulated in terms of spin-free replacement operators. We report numerical determinations of energies and 〈Sˆ(2)〉 expectation values, arising from our proposal as well as from traditional configuration interaction methods, in selected open-shell systems, in order to compare the behavior of these procedures and their computational costs.
Collapse
Affiliation(s)
- Diego R Alcoba
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires and Instituto de Física de Buenos Aires, Consejo Nacional de Investigaciones Científicas y Técnicas, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Alicia Torre
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco, Apdo. 644 E-48080 Bilbao, Spain
| | - Luis Lain
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco, Apdo. 644 E-48080 Bilbao, Spain
| | - Gustavo E Massaccesi
- Departamento de Ciencias Exactas, Ciclo Básico Común, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Ofelia B Oña
- Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicadas, Universidad Nacional de La Plata, CCT La Plata, Consejo Nacional de Investigaciones Científicas y Técnicas, Diag. 113 y 64 (S/N), Sucursal 4, CC 16, 1900 La Plata, Argentina
| | - Pablo Capuzzi
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires and Instituto de Física de Buenos Aires, Consejo Nacional de Investigaciones Científicas y Técnicas, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| |
Collapse
|
33
|
Alcoba DR, Torre A, Lain L, Massaccesi GE, Oña OB, Ayers PW, Van Raemdonck M, Bultinck P, Van Neck D. Performance of Shannon-entropy compacted N-electron wave functions for configuration interaction methods. Theor Chem Acc 2016. [DOI: 10.1007/s00214-016-1905-x] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
34
|
Fosso-Tande J, Nguyen TS, Gidofalvi G, DePrince AE. Large-Scale Variational Two-Electron Reduced-Density-Matrix-Driven Complete Active Space Self-Consistent Field Methods. J Chem Theory Comput 2016; 12:2260-71. [DOI: 10.1021/acs.jctc.6b00190] [Citation(s) in RCA: 95] [Impact Index Per Article: 11.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jacob Fosso-Tande
- Department
of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | - Truong-Son Nguyen
- Department of Chemistry and Biochemistry, Gonzaga University, Spokane, Washington 99258-0089, United States
| | - Gergely Gidofalvi
- Department of Chemistry and Biochemistry, Gonzaga University, Spokane, Washington 99258-0089, United States
| | - A. Eugene DePrince
- Department
of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| |
Collapse
|
35
|
Van Raemdonck M, Alcoba DR, Poelmans W, De Baerdemacker S, Torre A, Lain L, Massaccesi GE, Van Neck D, Bultinck P. Polynomial scaling approximations and dynamic correlation corrections to doubly occupied configuration interaction wave functions. J Chem Phys 2015; 143:104106. [PMID: 26374017 DOI: 10.1063/1.4930260] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
A class of polynomial scaling methods that approximate Doubly Occupied Configuration Interaction (DOCI) wave functions and improve the description of dynamic correlation is introduced. The accuracy of the resulting wave functions is analysed by comparing energies and studying the overlap between the newly developed methods and full configuration interaction wave functions, showing that a low energy does not necessarily entail a good approximation of the exact wave function. Due to the dependence of DOCI wave functions on the single-particle basis chosen, several orbital optimisation algorithms are introduced. An energy-based algorithm using the simulated annealing method is used as a benchmark. As a computationally more affordable alternative, a seniority number minimising algorithm is developed and compared to the energy based one revealing that the seniority minimising orbital set performs well. Given a well-chosen orbital basis, it is shown that the newly developed DOCI based wave functions are especially suitable for the computationally efficient description of static correlation and to lesser extent dynamic correlation.
Collapse
Affiliation(s)
- Mario Van Raemdonck
- Department of Inorganic and Physical Chemistry, Ghent University, Krijgslaan 281 (S3), 9000 Gent, Belgium
| | - Diego R Alcoba
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires and Instituto de Física de Buenos Aires, Consejo Nacional de Investigaciones Científicas y Técnicas, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Ward Poelmans
- Center for Molecular Modeling, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| | - Stijn De Baerdemacker
- Department of Inorganic and Physical Chemistry, Ghent University, Krijgslaan 281 (S3), 9000 Gent, Belgium
| | - Alicia Torre
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco, Apdo. 644, E-48080 Bilbao, Spain
| | - Luis Lain
- Departamento de Química Física, Facultad de Ciencia y Tecnología, Universidad del País Vasco, Apdo. 644, E-48080 Bilbao, Spain
| | - Gustavo E Massaccesi
- Departamento de Ciencias Exactas, Ciclo Básico Común, Universidad de Buenos Aires, Ciudad Universitaria, 1428 Buenos Aires, Argentina
| | - Dimitri Van Neck
- Center for Molecular Modeling, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| | - Patrick Bultinck
- Department of Inorganic and Physical Chemistry, Ghent University, Krijgslaan 281 (S3), 9000 Gent, Belgium
| |
Collapse
|