1
|
Mukhopadhyay T, Jangid B, Dutta AK. State-specific frozen natural orbital for reduced-cost algebraic diagrammatic construction calculations: The application to ionization problem. J Chem Phys 2023; 159:084113. [PMID: 37638624 DOI: 10.1063/5.0160024] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2023] [Accepted: 08/03/2023] [Indexed: 08/29/2023] Open
Abstract
We have developed a reduced-cost algebraic diagrammatic construction (ADC) method based on state-specific frozen natural orbital and natural auxiliary functions. The newly developed method has been benchmarked on the GW100 test set for the ionization problem. The use of state-specific natural orbitals drastically reduces the size of the virtual space with a systematically controllable accuracy and offers a significant speedup over the standard ionization potential (IP)-ADC(3) method. The accuracy of the method can be controlled by two thresholds and nearly a black box to use. The inclusion of the perturbative correction significantly improves the accuracy of the calculated IP values, and the efficiency of the method has been demonstrated by calculating the IP of a molecule with 60 atoms and more than 2216 basis functions.
Collapse
Affiliation(s)
- Tamoghna Mukhopadhyay
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Bhavnesh Jangid
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Achintya Kumar Dutta
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| |
Collapse
|
2
|
Dreuw A, Papapostolou A, Dempwolff AL. Algebraic Diagrammatic Construction Schemes Employing the Intermediate State Formalism: Theory, Capabilities, and Interpretation. J Phys Chem A 2023; 127:6635-6646. [PMID: 37498297 DOI: 10.1021/acs.jpca.3c02761] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/28/2023]
Abstract
Algebraic diagrammatic construction (ADC) schemes represent a family of ab initio methods for the calculation of excited electronic states and electron-detached and -attached states. All ADC methods have been demonstrated to possess great potential for molecular applications, e.g., for the calculation of absorption or photoelectron spectra or electron attachment processes. ADC originates from Green's function or propagator theory; however, most recent ADC developments heavily rely on the intermediate state representation or effective Liouvillian formalisms, which comprise new ADC methods and computational schemes for high-order properties. The different approaches for the calculation of excitation energies, ionization potentials, and electron affinities are intimately related, and they provide a coherent description of these quantities at equivalent levels of theory and with comparable errors. Most quantum chemical program packages contain ADC methods; however, the most complete ADC suite of methods can be found in the recent release of Q-Chem.
Collapse
Affiliation(s)
- Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, 69120 Heidelberg, Germany
| | - Antonia Papapostolou
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, 69120 Heidelberg, Germany
| | - Adrian L Dempwolff
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, 69120 Heidelberg, Germany
| |
Collapse
|
3
|
Mester D, Kállay M. Vertical Ionization Potentials and Electron Affinities at the Double-Hybrid Density Functional Level. J Chem Theory Comput 2023. [PMID: 37326360 DOI: 10.1021/acs.jctc.3c00363] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
The double-hybrid (DH) time-dependent density functional theory is extended to vertical ionization potentials (VIPs) and electron affinities (VEAs). Utilizing the density fitting approximation, efficient implementations are presented for the genuine DH ansatz relying on the perturbative second-order correction, while an iterative analogue is also elaborated using our second-order algebraic-diagrammatic construction [ADC(2)]-based DH approach. The favorable computational requirements of the present schemes are discussed in detail. The performance of the recently proposed spin-component-scaled and spin-opposite-scaled (SOS) range-separated (RS) and long-range corrected (LC) DH functionals is comprehensively assessed, while popular hybrid and global DH approaches are also discussed. For the benchmark calculations, up-to-date test sets are selected with high-level coupled-cluster references. Our results show that the ADC(2)-based SOS-RS-PBE-P86 approach is the most accurate and robust functional. This method consistently outperforms the excellent SOS-ADC(2) approach for VIPs, although the results are somewhat less satisfactory for VEAs. Among the genuine DH functionals, the SOS-ωPBEPP86 approach is also recommended for describing ionization processes, but its performance is even less reliable for electron-attached states. In addition, surprisingly good results are attained by the LC hybrid ωB97X-D functional, where the corresponding occupied (unoccupied) orbital energies are retrieved as VIPs (VEAs) within the present formalism.
Collapse
Affiliation(s)
- Dávid Mester
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- MTA-BME Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Faculty of Chemical Technology and Biotechnology, Budapest University of Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- MTA-BME Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| |
Collapse
|
4
|
Asthana A, Kumar A, Abraham V, Grimsley H, Zhang Y, Cincio L, Tretiak S, Dub PA, Economou SE, Barnes E, Mayhall NJ. Quantum self-consistent equation-of-motion method for computing molecular excitation energies, ionization potentials, and electron affinities on a quantum computer. Chem Sci 2023; 14:2405-2418. [PMID: 36873839 PMCID: PMC9977410 DOI: 10.1039/d2sc05371c] [Citation(s) in RCA: 14] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2022] [Accepted: 01/26/2023] [Indexed: 01/30/2023] Open
Abstract
Near-term quantum computers are expected to facilitate material and chemical research through accurate molecular simulations. Several developments have already shown that accurate ground-state energies for small molecules can be evaluated on present-day quantum devices. Although electronically excited states play a vital role in chemical processes and applications, the search for a reliable and practical approach for routine excited-state calculations on near-term quantum devices is ongoing. Inspired by excited-state methods developed for the unitary coupled-cluster theory in quantum chemistry, we present an equation-of-motion-based method to compute excitation energies following the variational quantum eigensolver algorithm for ground-state calculations on a quantum computer. We perform numerical simulations on H2, H4, H2O, and LiH molecules to test our quantum self-consistent equation-of-motion (q-sc-EOM) method and compare it to other current state-of-the-art methods. q-sc-EOM makes use of self-consistent operators to satisfy the vacuum annihilation condition, a critical property for accurate calculations. It provides real and size-intensive energy differences corresponding to vertical excitation energies, ionization potentials and electron affinities. We also find that q-sc-EOM is more suitable for implementation on NISQ devices as it is expected to be more resilient to noise compared with the currently available methods.
Collapse
Affiliation(s)
- Ayush Asthana
- Department of Chemistry, Virginia Tech Blacksburg 24061 VA USA
- Virginia Tech Center for Quantum Information Science and Engineering Blacksburg 24061 VA USA
| | - Ashutosh Kumar
- Theoretical Division, Los Alamos National Laboratory Los Alamos 87545 NM USA
| | - Vibin Abraham
- Department of Chemistry, University of Michigan Ann Arbor 48109 MI USA
| | - Harper Grimsley
- Department of Chemistry, Virginia Tech Blacksburg 24061 VA USA
- Virginia Tech Center for Quantum Information Science and Engineering Blacksburg 24061 VA USA
| | - Yu Zhang
- Theoretical Division, Los Alamos National Laboratory Los Alamos 87545 NM USA
| | - Lukasz Cincio
- Theoretical Division, Los Alamos National Laboratory Los Alamos 87545 NM USA
| | - Sergei Tretiak
- Theoretical Division, Los Alamos National Laboratory Los Alamos 87545 NM USA
- Center for Integrated Nanotechnologies, Los Alamos National Laboratory Los Alamos 87545 NM USA
| | - Pavel A Dub
- Chemistry Division, Los Alamos National Laboratory Los Alamos 87545 NM USA
| | - Sophia E Economou
- Department of Physics, Virginia Tech Blacksburg 24061 VA USA
- Virginia Tech Center for Quantum Information Science and Engineering Blacksburg 24061 VA USA
| | - Edwin Barnes
- Department of Physics, Virginia Tech Blacksburg 24061 VA USA
- Virginia Tech Center for Quantum Information Science and Engineering Blacksburg 24061 VA USA
| | - Nicholas J Mayhall
- Department of Chemistry, Virginia Tech Blacksburg 24061 VA USA
- Virginia Tech Center for Quantum Information Science and Engineering Blacksburg 24061 VA USA
| |
Collapse
|
5
|
Banerjee S, Sokolov AY. Non-Dyson Algebraic Diagrammatic Construction Theory for Charged Excitations in Solids. J Chem Theory Comput 2022; 18:5337-5348. [PMID: 35976918 DOI: 10.1021/acs.jctc.2c00565] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present the first implementation and applications of non-Dyson algebraic diagrammatic construction theory for charged excitations in three-dimensional periodic solids (EA/IP-ADC). The EA/IP-ADC approach has a computational cost similar to the ground-state Møller-Plesset perturbation theory, enabling efficient calculations of a variety of crystalline excited-state properties (e.g., band structure, band gap, density of states) sampled in the Brillouin zone. We use EA/IP-ADC to compute the quasiparticle band structures and band gaps of several materials (from large-gap atomic and ionic solids to small-gap semiconductors) and analyze the errors of EA/IP-ADC approximations up to the third order in perturbation theory. Our work also reports the first-ever calculations of ground-state properties (equation-of-state and lattice constants) of three-dimensional crystalline systems using a periodic implementation of third-order Møller-Plesset perturbation theory (MP3).
Collapse
Affiliation(s)
- Samragni Banerjee
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Alexander Yu Sokolov
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
6
|
Dempwolff AL, Hodecker M, Dreuw A. Vertical ionization potential benchmark for unitary coupled-cluster and algebraic-diagrammatic construction methods. J Chem Phys 2022; 156:054114. [DOI: 10.1063/5.0079047] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Adrian L. Dempwolff
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Manuel Hodecker
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
- Department of Theoretical Chemistry and Biology, KTH Royal Institute of Technology, Malvinas Väg 10, 114 28 Stockholm, Sweden
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| |
Collapse
|