51
|
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
|
52
|
Zhou C, Chen Z, Wu W. Reciprocal transformation of seniority number restricted wave function. J Chem Phys 2018; 149:044111. [DOI: 10.1063/1.5032129] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
Affiliation(s)
- Chen Zhou
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Zhenhua Chen
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Wei Wu
- The State Key Laboratory of Physical Chemistry of Solid Surfaces, iChem and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| |
Collapse
|
53
|
van Meer R, Gritsenko OV, Baerends EJ. A non-JKL density matrix functional for intergeminal correlation between closed-shell geminals from analysis of natural orbital configuration interaction expansions. J Chem Phys 2018; 148:104102. [DOI: 10.1063/1.5018671] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- R. van Meer
- Computational Molecular Science Research Team, RIKEN Advanced Institute for Computational Science, Kobe, Japan
- Section Theoretical Chemistry, VU University, Amsterdam, The Netherlands
| | - O. V. Gritsenko
- Section Theoretical Chemistry, VU University, Amsterdam, The Netherlands
| | - E. J. Baerends
- Section Theoretical Chemistry, VU University, Amsterdam, The Netherlands
| |
Collapse
|
54
|
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
|
55
|
Benavides-Riveros CL, Lathiotakis NN, Marques MAL. Towards a formal definition of static and dynamic electronic correlations. Phys Chem Chem Phys 2018; 19:12655-12664. [PMID: 28474027 DOI: 10.1039/c7cp01137g] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Some of the most spectacular failures of density-functional and Hartree-Fock theories are related to an incorrect description of the so-called static electron correlation. Motivated by recent progress in the N-representability problem of the one-body density matrix for pure states, we propose a method to quantify the static contribution to the electronic correlation. By studying several molecular systems we show that our proposal correlates well with our intuition of static and dynamic electron correlation. Our results bring out the paramount importance of the occupancy of the highest occupied natural spin-orbital in such quantification.
Collapse
|
56
|
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
|
57
|
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]
|
58
|
Wahlen-Strothman JM, Henderson TM, Scuseria GE. Influence of broken-pair excitations on the exact pair wavefunction. Mol Phys 2017. [DOI: 10.1080/00268976.2017.1372646] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
| | - Thomas M. Henderson
- Department of Physics Astronomy, Rice University, Houston, TX, USA
- Department of Chemistry, Rice University, Houston, TX, USA
| | - Gustavo E. Scuseria
- Department of Physics Astronomy, Rice University, Houston, TX, USA
- Department of Chemistry, Rice University, Houston, TX, USA
| |
Collapse
|
59
|
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]
|
60
|
A comparative study on seniority-based MO and VB calculations of the singlet and triplet energy gaps of open-shell molecules. COMPUT THEOR CHEM 2017. [DOI: 10.1016/j.comptc.2017.03.031] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
61
|
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
|
62
|
Affiliation(s)
- Ágnes Szabados
- Laboratory of Theoretical Chemistry, Loránd Eötvös University, Budapest, Hungary
| | - Ádám Margócsy
- Laboratory of Theoretical Chemistry, Loránd Eötvös University, Budapest, Hungary
| |
Collapse
|
63
|
Affiliation(s)
| | - Markus Reiher
- Laboratorium für Physikalische Chemie, ETH Zürich, Zürich, Switzerland
| |
Collapse
|
64
|
Varsano D, Caprasecca S, Coccia E. Theoretical description of protein field effects on electronic excitations of biological chromophores. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:013002. [PMID: 27830666 DOI: 10.1088/0953-8984/29/1/013002] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Photoinitiated phenomena play a crucial role in many living organisms. Plants, algae, and bacteria absorb sunlight to perform photosynthesis, and convert water and carbon dioxide into molecular oxygen and carbohydrates, thus forming the basis for life on Earth. The vision of vertebrates is accomplished in the eye by a protein called rhodopsin, which upon photon absorption performs an ultrafast isomerisation of the retinal chromophore, triggering the signal cascade. Many other biological functions start with the photoexcitation of a protein-embedded pigment, followed by complex processes comprising, for example, electron or excitation energy transfer in photosynthetic complexes. The optical properties of chromophores in living systems are strongly dependent on the interaction with the surrounding environment (nearby protein residues, membrane, water), and the complexity of such interplay is, in most cases, at the origin of the functional diversity of the photoactive proteins. The specific interactions with the environment often lead to a significant shift of the chromophore excitation energies, compared with their absorption in solution or gas phase. The investigation of the optical response of chromophores is generally not straightforward, from both experimental and theoretical standpoints; this is due to the difficulty in understanding diverse behaviours and effects, occurring at different scales, with a single technique. In particular, the role played by ab initio calculations in assisting and guiding experiments, as well as in understanding the physics of photoactive proteins, is fundamental. At the same time, owing to the large size of the systems, more approximate strategies which take into account the environmental effects on the absorption spectra are also of paramount importance. Here we review the recent advances in the first-principle description of electronic and optical properties of biological chromophores embedded in a protein environment. We show their applications on paradigmatic systems, such as the light-harvesting complexes, rhodopsin and green fluorescent protein, emphasising the theoretical frameworks which are of common use in solid state physics, and emerging as promising tools for biomolecular systems.
Collapse
Affiliation(s)
- Daniele Varsano
- S3 Center, CNR Institute of Nanoscience, Via Campi 213/A, 41125 Modena, Italy
| | | | | |
Collapse
|
65
|
Gomez JA, Degroote M, Zhao J, Qiu Y, Scuseria GE. Spin polynomial similarity transformation for repulsive Hamiltonians: interpolating between coupled cluster and spin-projected unrestricted Hartree–Fock. Phys Chem Chem Phys 2017; 19:22385-22394. [DOI: 10.1039/c7cp04075j] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Our overarching goal is to be able to describe both weak and strong correlation with a single, computationally affordable method without sacrificing important qualities of the wavefunction, e.g. symmetries of the Hamiltonian.
Collapse
Affiliation(s)
- John A. Gomez
- Department of Chemistry
- Rice University
- Houston
- USA
- Applied Physics Program
| | | | - Jinmo Zhao
- Department of Chemistry
- Rice University
- Houston
- USA
| | - Yiheng Qiu
- Department of Chemistry
- Rice University
- Houston
- USA
| | | |
Collapse
|
66
|
Zhao L, Neuscamman E. Amplitude Determinant Coupled Cluster with Pairwise Doubles. J Chem Theory Comput 2016; 12:5841-5850. [DOI: 10.1021/acs.jctc.6b00812] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Luning Zhao
- Department
of Chemistry, University of California, Berkeley, California 94720, United States
| | - Eric Neuscamman
- Department
of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical
Science Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
67
|
Affiliation(s)
- Peter A. Limacher
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany
| |
Collapse
|
68
|
Abstract
We introduce an energy functional for ground-state electronic structure calculations. Its variables are the natural spin-orbitals of singlet many-body wave functions and their joint occupation probabilities deriving from controlled approximations to the two-particle density matrix that yield algebraic scaling in general, and Hartree-Fock scaling in its seniority-zero version. Results from the latter version for small molecular systems are compared with those of highly accurate quantum-chemical computations. The energies lie above full configuration interaction calculations, close to doubly occupied configuration interaction calculations. Their accuracy is considerably greater than that obtained from current density-functional theory approximations and from current functionals of the one-particle density matrix.
Collapse
Affiliation(s)
- Ralph Gebauer
- The Abdus Salam International Centre for Theoretical Physics, 34151 Trieste, Italy
| | - Morrel H Cohen
- Department of Physics and Astronomy, Rutgers University, Piscataway, NJ 08854;
- Department of Chemistry, Princeton University, Princeton, NJ 08544
| | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, NJ 08544;
- Department of Physics, Princeton University, Princeton, NJ 08544
| |
Collapse
|
69
|
Bhattacharya D, Shil S, Misra A, Bytautas L, Klein DJ. Toward Molecular Magnets of Organic Origin via Anion−π Interaction Involving m-Aminyl Diradical: A Theoretical Study. J Phys Chem A 2016; 120:9117-9130. [DOI: 10.1021/acs.jpca.6b09666] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Debojit Bhattacharya
- Department of Marine Sciences, Texas A&M University at Galveston, Galveston, Texas 77553, United States
| | - Suranjan Shil
- Center
for Atomic Scale Materials Design, Department of Physics, Technical University of Denmark, Kgs. Lyngby 2800, Denmark
| | - Anirban Misra
- Department
of Chemistry, University of North Bengal, Darjeeling, PIN. 734013, West Bengal, India
| | - Laimutis Bytautas
- Department
of Chemistry, Galveston College, 4015 Avenue Q, Galveston, Texas 77550, United States
| | - Douglas J. Klein
- Department of Marine Sciences, Texas A&M University at Galveston, Galveston, Texas 77553, United States
| |
Collapse
|
70
|
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
|
71
|
Zen A, Coccia E, Gozem S, Olivucci M, Guidoni L. Quantum Monte Carlo Treatment of the Charge Transfer and Diradical Electronic Character in a Retinal Chromophore Minimal Model. J Chem Theory Comput 2016; 11:992-1005. [PMID: 25821414 PMCID: PMC4357234 DOI: 10.1021/ct501122z] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2014] [Indexed: 01/22/2023]
Abstract
![]()
The
penta-2,4-dieniminium cation (PSB3) displays similar ground
state and first excited state potential energy features as those of
the retinal protonated Schiff base (RPSB) chromophore in rhodopsin.
Recently, PSB3 has been used to benchmark several electronic structure
methods, including highly correlated multireference wave function
approaches, highlighting the necessity to accurately describe the
electronic correlation in order to obtain reliable properties even
along the ground state (thermal) isomerization paths. In this work,
we apply two quantum Monte Carlo approaches, the variational Monte
Carlo and the lattice regularized diffusion Monte Carlo, to study
the energetics and electronic properties of PSB3 along representative
minimum energy paths and scans related to its thermal cis–trans isomerization. Quantum Monte Carlo
is used in combination with the Jastrow antisymmetrized geminal power
ansatz, which guarantees an accurate and balanced description of the
static electronic correlation thanks to the multiconfigurational nature
of the antisymmetrized geminal power term, and of the dynamical correlation,
due to the presence of the Jastrow factor explicitly depending on
electron–electron distances. Along the two ground state isomerization
minimum energy paths of PSB3, CASSCF calculations yield wave functions
having either charge transfer or diradical character in proximity
of the two transition state configurations. Here, we observe that
at the quantum Monte Carlo level of theory, only the transition state
with charge transfer character can be located. The conical intersection,
which becomes highly sloped, is observed only if the path connecting
the two original CASSCF transition states is extended beyond the diradical
one, namely by increasing the bond-length-alternation (BLA). These
findings are in good agreement with the results obtained by MRCISD+Q
calculations, and they demonstrate the importance of having an accurate
description of the static and dynamical correlation when studying
isomerization and transition states of conjugated systems.
Collapse
|
72
|
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]
|
73
|
Shepherd JJ, Henderson TM, Scuseria GE. Using full configuration interaction quantum Monte Carlo in a seniority zero space to investigate the correlation energy equivalence of pair coupled cluster doubles and doubly occupied configuration interaction. J Chem Phys 2016; 144:094112. [DOI: 10.1063/1.4942770] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Affiliation(s)
- James J. Shepherd
- Department of Chemistry and Department of Physics and Astronomy, Rice University, Houston, Texas 77005-1892, USA
| | - Thomas M. Henderson
- Department of Chemistry and Department of Physics and Astronomy, Rice University, Houston, Texas 77005-1892, USA
| | - Gustavo E. Scuseria
- Department of Chemistry and Department of Physics and Astronomy, Rice University, Houston, Texas 77005-1892, USA
| |
Collapse
|
74
|
Garza AJ, Bulik IW, Alencar AGS, Sun J, Perdew JP, Scuseria GE. Combinations of coupled cluster, density functionals, and the random phase approximation for describing static and dynamic correlation, and van der Waals interactions. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1123315] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Affiliation(s)
| | | | | | - Jianwei Sun
- Department of Physics, Temple University, Philadelphia, PA, USA
| | - John P. Perdew
- Department of Physics, Temple University, Philadelphia, PA, USA
- Department of Chemistry, Temple University, Philadelphia, PA, USA
| | - Gustavo E. Scuseria
- Department of Chemistry, Rice University, Houston, TX, USA
- Department of Physics and Astronomy, Rice University, Houston, TX, USA
| |
Collapse
|
75
|
Barborini M, Coccia E. Investigating Disjoint Non-Kekulé Diradicals with Quantum Monte Carlo: The Tetramethyleneethane Molecule through the Jastrow Antisymmetrized Geminal Power Wave Function. J Chem Theory Comput 2015; 11:5696-704. [DOI: 10.1021/acs.jctc.5b00819] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Matteo Barborini
- S3
Research Center, CNR-NANO, Via Campi 213/a, 41125 Modena, Modena, Italy
| | - Emanuele Coccia
- Laboratoire de Chimie Théorique, Sorbonne Universités, UPMC Univ Paris 06 & CNRS, UMR 7616, F-75005 Paris, France
| |
Collapse
|
76
|
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
|
77
|
Barborini M, Guidoni L. Ground State Geometries of Polyacetylene Chains from Many-Particle Quantum Mechanics. J Chem Theory Comput 2015; 11:4109-18. [PMID: 26405437 PMCID: PMC4570661 DOI: 10.1021/acs.jctc.5b00427] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2015] [Indexed: 11/28/2022]
Abstract
Due to the crucial role played by electron correlation, the accurate determination of ground state geometries of π-conjugated molecules is still a challenge for many quantum chemistry methods. Because of the high parallelism of the algorithms and their explicit treatment of electron correlation effects, Quantum Monte Carlo calculations can offer an accurate and reliable description of the electronic states and of the geometries of such systems, competing with traditional quantum chemistry approaches. Here, we report the structural properties of polyacetylene chains H-(C₂H₂)(N)-H up to N = 12 acetylene units, by means of Variational Monte Carlo (VMC) calculations based on the multi-determinant Jastrow Antisymmetrized Geminal Power (JAGP) wave function. This compact ansatz can provide for such systems an accurate description of the dynamical electronic correlation as recently detailed for the 1,3-butadiene molecule [J. Chem. Theory Comput. 2015 11 (2), 508-517]. The calculated Bond Length Alternation (BLA), namely the difference between the single and double carbon bonds, extrapolates, for N → ∞, to a value of 0.0910(7) Å, compatible with the experimental data. An accurate analysis was able to distinguish between the influence of the multi-determinantal AGP expansion and of the Jastrow factor on the geometrical properties of the fragments. Our size-extensive and self-interaction-free results provide new and accurate ab initio references for the structures of the ground state of polyenes.
Collapse
Affiliation(s)
- Matteo Barborini
- Dipartimento
di Ingegneria e Scienze dell’Informazione e Matematica, Università degli studi dell’Aquila, Via Vetoio (Coppito), 67100 L’Aquila, Italy
- Dipartimento
di Scienze Fisiche e Chimiche, Università
degli studi dell’Aquila, Via Vetoio (Coppito), 67100 L’Aquila, Italy
| | - Leonardo Guidoni
- Dipartimento
di Scienze Fisiche e Chimiche, Università
degli studi dell’Aquila, Via Vetoio (Coppito), 67100 L’Aquila, Italy
| |
Collapse
|
78
|
Bytautas L, Scuseria GE, Ruedenberg K. Seniority number description of potential energy surfaces: Symmetric dissociation of water, N2, C2, and Be2. J Chem Phys 2015; 143:094105. [DOI: 10.1063/1.4929904] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Laimutis Bytautas
- Department of Chemistry, Galveston College, 4015 Ave. Q, Galveston, Texas 77550, USA
| | - Gustavo E. Scuseria
- Department of Chemistry, Rice University, Houston, Texas 77005, USA
- Chemistry Department, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia
| | - Klaus Ruedenberg
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, USA
| |
Collapse
|
79
|
Poelmans W, Van Raemdonck M, Verstichel B, De Baerdemacker S, Torre A, Lain L, Massaccesi GE, Alcoba DR, Bultinck P, Van Neck D. Variational Optimization of the Second-Order Density Matrix Corresponding to a Seniority-Zero Configuration Interaction Wave Function. J Chem Theory Comput 2015; 11:4064-76. [DOI: 10.1021/acs.jctc.5b00378] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Ward Poelmans
- Center
for Molecular Modeling, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| | - Mario Van Raemdonck
- Department
of Inorganic and Physical Chemistry, Ghent University, Krijgslaan
281 (S3), 9000 Gent, Belgium
| | - Brecht Verstichel
- 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
- 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
| | - 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
| | - Patrick Bultinck
- Department
of Inorganic and Physical Chemistry, Ghent University, Krijgslaan
281 (S3), 9000 Gent, Belgium
| | - Dimitri Van Neck
- Center
for Molecular Modeling, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| |
Collapse
|
80
|
Abstract
A new single-pair operator for seniority-zero wave functions is introduced closely related to the single-particle Fock operator from conventional Hartree-Fock theory. This allows one to ascribe orbital energies in the context of model Hamiltonians, for which no single-particle operators exist. Several applications demonstrate the usefulness of these orbital energies. In analogy to Koopmans' theorem, atomic double ionization potentials are successfully predicted. A computationally efficient second-order perturbation scheme for seniority-zero wave functions scaling quadratically with system size is defined and applied to the dissociation of nitrogen and to strongly correlated two-dimensional Heisenberg lattices. An extension of the method for full seniority is presented leading to an intruder-free single reference perturbation theory with improved asymptotic convergence.
Collapse
Affiliation(s)
- Peter A Limacher
- Institute of Physical Chemistry, Karlsruhe Institute of Technology , D-76131 Karlsruhe, Germany
| |
Collapse
|
81
|
Garza AJ, Bulik IW, Henderson TM, Scuseria GE. Range separated hybrids of pair coupled cluster doubles and density functionals. Phys Chem Chem Phys 2015; 17:22412-22. [PMID: 26249820 DOI: 10.1039/c5cp02773j] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Pair coupled cluster doubles (pCCD) is a size-consistent, size-extensive, low-cost simplification of CCD that has been shown to be able to describe static correlation without breaking symmetry. We combine pCCD with Kohn-Sham functionals of the density and the local pair density in order to incorporate dynamic correlation in pCCD while maintaining its low cost. Double counting is eliminated by splitting the (interelectron) Coulomb operator into complementary short- and long-range parts, and evaluating the two-body energy with pCCD in the long-range and with density functionals in the short-range. This simultaneously suppresses self-interaction in the Hartree-exchange term of the functionals. Generalizations including a fraction of wavefunction two-body energy in the short-range are also derived and studied. The improvement of our pCCD+DFT hybrids over pCCD is demonstrated in calculations on benchmarks where both types of correlation are important.
Collapse
Affiliation(s)
- Alejandro J Garza
- Department of Chemistry, Rice University, Houston, Texas 77251-1892, USA.
| | | | | | | |
Collapse
|
82
|
Affiliation(s)
- Zhenhua Chen
- The State Key Laboratory
of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial
Key Laboratory of Theoretical and Computational Chemistry and College
of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Chen Zhou
- The State Key Laboratory
of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial
Key Laboratory of Theoretical and Computational Chemistry and College
of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| | - Wei Wu
- The State Key Laboratory
of Physical Chemistry of Solid Surfaces, iChEM, Fujian Provincial
Key Laboratory of Theoretical and Computational Chemistry and College
of Chemistry and Chemical Engineering, Xiamen University, Xiamen, Fujian 361005, China
| |
Collapse
|
83
|
A study of the compactness of wave functions based on Shannon entropy indices: a seniority number approach. Theor Chem Acc 2015. [DOI: 10.1007/s00214-015-1688-5] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
84
|
Henderson TM, Bulik IW, Scuseria GE. Pair extended coupled cluster doubles. J Chem Phys 2015; 142:214116. [DOI: 10.1063/1.4921986] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- Thomas M. Henderson
- Department of Chemistry, Rice University, Houston, Texas 77005-1892, USA
- Department of Physics and Astronomy, Rice University, Houston, Texas 77005-1892, USA
| | - Ireneusz W. Bulik
- Department of Chemistry, Rice University, Houston, Texas 77005-1892, USA
| | - Gustavo E. Scuseria
- Department of Chemistry, Rice University, Houston, Texas 77005-1892, USA
- Department of Physics and Astronomy, Rice University, Houston, Texas 77005-1892, USA
| |
Collapse
|
85
|
Small DW, Sundstrom EJ, Head-Gordon M. A simple way to test for collinearity in spin symmetry broken wave functions: General theory and application to generalized Hartree Fock. J Chem Phys 2015; 142:094112. [DOI: 10.1063/1.4913740] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- David W. Small
- Department of Chemistry, University of California, Berkeley, California 94720, USA and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Eric J. Sundstrom
- Department of Chemistry, University of California, Berkeley, California 94720, USA and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
86
|
Henderson TM, Bulik IW, Stein T, Scuseria GE. Seniority-based coupled cluster theory. J Chem Phys 2014; 141:244104. [DOI: 10.1063/1.4904384] [Citation(s) in RCA: 90] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Affiliation(s)
- Thomas M. Henderson
- Department of Chemistry, Rice University, Houston, Texas 77005-1892, USA
- Department of Physics and Astronomy, Rice University, Houston, Texas 77005-1892, USA
| | - Ireneusz W. Bulik
- Department of Chemistry, Rice University, Houston, Texas 77005-1892, USA
| | - Tamar Stein
- Department of Chemistry, Rice University, Houston, Texas 77005-1892, USA
| | - Gustavo E. Scuseria
- Department of Chemistry, Rice University, Houston, Texas 77005-1892, USA
- Department of Physics and Astronomy, Rice University, Houston, Texas 77005-1892, USA
| |
Collapse
|
87
|
Alcoba DR, Torre A, Lain L, Oña OB, Capuzzi P, Van Raemdonck M, Bultinck P, Van Neck D. A hybrid configuration interaction treatment based on seniority number and excitation schemes. J Chem Phys 2014; 141:244118. [DOI: 10.1063/1.4904755] [Citation(s) in RCA: 32] [Impact Index Per Article: 3.2] [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 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
| | - 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
| | - Mario Van Raemdonck
- Department of Inorganic and Physical Chemistry, Ghent University, Krijgslaan 281 (S3), 9000 Gent, Belgium
| | - Patrick Bultinck
- Department of Inorganic and Physical Chemistry, Ghent University, Krijgslaan 281 (S3), 9000 Gent, Belgium
| | - Dimitri Van Neck
- Center for Molecular Modelling, Ghent University, Technologiepark 903, 9052 Zwijnaarde, Belgium
| |
Collapse
|
88
|
Abstract
A method is suggested to build simple multiconfigurational wave functions specified uniquely by an energy cutoff Λ. These are constructed from a model space containing determinants with energy relative to that of the most stable determinant no greater than Λ. The resulting Λ-CI wave function is adaptive, being able to represent both single-reference and multireference electronic states. We also consider a more compact wave function parameterization (Λ+SD-CI), which is based on a small Λ-CI reference and adds a selection of all the singly and doubly excited determinants generated from it. We report two heuristic algorithms to build Λ-CI wave functions. The first is based on an approximate prescreening of the full configuration interaction space, while the second performs a breadth-first search coupled with pruning. The Λ-CI and Λ+SD-CI approaches are used to compute the dissociation curve of N2 and the potential energy curves for the first three singlet states of C2. Special attention is paid to the issue of energy discontinuities caused by changes in the size of the Λ-CI wave function along the potential energy curve. This problem is shown to be solvable by smoothing the matrix elements of the Hamiltonian. Our last example, involving the Cu2O2(2+) core, illustrates an alternative use of the Λ-CI method: as a tool to both estimate the multireference character of a wave function and to create a compact model space to be used in subsequent high-level multireference coupled cluster computations.
Collapse
Affiliation(s)
- Francesco A Evangelista
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
89
|
Boguslawski K, Tecmer P, Bultinck P, De Baerdemacker S, Van Neck D, Ayers PW. Nonvariational Orbital Optimization Techniques for the AP1roG Wave Function. J Chem Theory Comput 2014; 10:4873-82. [DOI: 10.1021/ct500759q] [Citation(s) in RCA: 45] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Katharina Boguslawski
- Department
of Chemistry and Chemical Biology, McMaster University, 1280 Main
Street West, Hamilton, L8S 4M1 Ontario, Canada
| | - Paweł Tecmer
- Department
of Chemistry and Chemical Biology, McMaster University, 1280 Main
Street West, Hamilton, L8S 4M1 Ontario, Canada
| | - Patrick Bultinck
- Department
of Inorganic and Physical Chemistry, Ghent University, Krijgslaan
281 (S3), 9000 Gent, East Flanders, Belgium
| | - Stijn De Baerdemacker
- Center
for Molecular Modelling, Ghent University, Technologiepark 903, 9052 Gent, East
Flanders, Belgium
| | - Dimitri Van Neck
- Center
for Molecular Modelling, Ghent University, Technologiepark 903, 9052 Gent, East
Flanders, Belgium
| | - Paul W. Ayers
- Department
of Chemistry and Chemical Biology, McMaster University, 1280 Main
Street West, Hamilton, L8S 4M1 Ontario, Canada
| |
Collapse
|
90
|
Alcoba DR, Torre A, Lain L, Massaccesi GE, Oña OB. Configuration interaction wave functions: A seniority number approach. J Chem Phys 2014; 140:234103. [DOI: 10.1063/1.4882881] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
|
91
|
Boguslawski K, Tecmer P, Limacher PA, Johnson PA, Ayers PW, Bultinck P, De Baerdemacker S, Van Neck D. Projected seniority-two orbital optimization of the antisymmetric product of one-reference orbital geminal. J Chem Phys 2014; 140:214114. [PMID: 24907997 DOI: 10.1063/1.4880820] [Citation(s) in RCA: 58] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
Affiliation(s)
- Katharina Boguslawski
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, 1280 Main Street West, Ontario L8S 4M1, Canada
| | - Paweł Tecmer
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, 1280 Main Street West, Ontario L8S 4M1, Canada
| | - Peter A. Limacher
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, 1280 Main Street West, Ontario L8S 4M1, Canada
| | - Paul A. Johnson
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, 1280 Main Street West, Ontario L8S 4M1, Canada
| | - Paul W. Ayers
- Department of Chemistry and Chemical Biology, McMaster University, Hamilton, 1280 Main Street West, Ontario L8S 4M1, Canada
| | - Patrick Bultinck
- Department of Inorganic and Physical Chemistry, Ghent University, Krijgslaan 281 (S3), 9000 Gent, Belgium
| | - Stijn De Baerdemacker
- Center for Molecular Modelling, Ghent University, Technologiepark 903, 9052 Gent, Belgium
| | - Dimitri Van Neck
- Center for Molecular Modelling, Ghent University, Technologiepark 903, 9052 Gent, Belgium
| |
Collapse
|
92
|
Stein T, Henderson TM, Scuseria GE. Seniority zero pair coupled cluster doubles theory. J Chem Phys 2014; 140:214113. [DOI: 10.1063/1.4880819] [Citation(s) in RCA: 125] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Tamar Stein
- Department of Chemistry, Rice University, Houston, Texas 77251, USA
| | - Thomas M. Henderson
- Department of Chemistry, Rice University, Houston, Texas 77251, USA
- Department of Physics and Astronomy, Rice University, Houston, Texas 77251, USA
| | - Gustavo E. Scuseria
- Department of Chemistry, Rice University, Houston, Texas 77251, USA
- Department of Physics and Astronomy, Rice University, Houston, Texas 77251, USA
| |
Collapse
|
93
|
Mentel ŁM, van Meer R, Gritsenko OV, Baerends EJ. The density matrix functional approach to electron correlation: Dynamic and nondynamic correlation along the full dissociation coordinate. J Chem Phys 2014; 140:214105. [DOI: 10.1063/1.4879776] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
94
|
|
95
|
Zen A, Coccia E, Luo Y, Sorella S, Guidoni L. Static and Dynamical Correlation in Diradical Molecules by Quantum Monte Carlo Using the Jastrow Antisymmetrized Geminal Power Ansatz. J Chem Theory Comput 2014; 10:1048-61. [DOI: 10.1021/ct401008s] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Affiliation(s)
- Andrea Zen
- Dipartimento
di Fisica, Sapienza-Università di Roma, Piazzale Aldo Moro
5, 00185 Rome, Italy
| | - Emanuele Coccia
- Dipartimento
di Scienze Fisiche e Chimiche, Università degli Studi dell’Aquila, Via Vetoio 2, 67010 L’Aquila, Italy
| | - Ye Luo
- Scuola
Internazionale Superiore di Studi Avanzati (SISSA) and Democritos
National Simulation Center, Istituto Officina dei Materiali del CNR, via Bonomea 265, 34136 Trieste, Italy
| | - Sandro Sorella
- Scuola
Internazionale Superiore di Studi Avanzati (SISSA) and Democritos
National Simulation Center, Istituto Officina dei Materiali del CNR, via Bonomea 265, 34136 Trieste, Italy
| | - Leonardo Guidoni
- Dipartimento
di Scienze Fisiche e Chimiche, Università degli Studi dell’Aquila, Via Vetoio 2, 67010 L’Aquila, Italy
| |
Collapse
|
96
|
Booth GH, Smart SD, Alavi A. Linear-scaling and parallelisable algorithms for stochastic quantum chemistry. Mol Phys 2014. [DOI: 10.1080/00268976.2013.877165] [Citation(s) in RCA: 78] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
97
|
Limacher PA, Kim TD, Ayers PW, Johnson PA, De Baerdemacker S, Van Neck D, Bultinck P. The influence of orbital rotation on the energy of closed-shell wavefunctions. Mol Phys 2014. [DOI: 10.1080/00268976.2013.874600] [Citation(s) in RCA: 61] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
98
|
Limacher PA, Ayers PW, Johnson PA, De Baerdemacker S, Neck DV, Bultinck P. Simple and inexpensive perturbative correction schemes for antisymmetric products of nonorthogonal geminals. Phys Chem Chem Phys 2014; 16:5061-5. [DOI: 10.1039/c3cp53301h] [Citation(s) in RCA: 37] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
|
99
|
Giesbertz K. Are natural orbitals useful for generating an efficient expansion of the wave function? Chem Phys Lett 2014. [DOI: 10.1016/j.cplett.2013.11.038] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
100
|
Alcoba DR, Torre A, Lain L, Massaccesi GE, Oña OB. Seniority number in spin-adapted spaces and compactness of configuration interaction wave functions. J Chem Phys 2013; 139:084103. [DOI: 10.1063/1.4818755] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|