1
|
Majee K, Chakraborty S, Mukhopadhyay T, Nayak MK, Dutta AK. A reduced cost four-component relativistic unitary coupled cluster method for atoms and molecules. J Chem Phys 2024; 161:034101. [PMID: 39007370 DOI: 10.1063/5.0207091] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2024] [Accepted: 06/12/2024] [Indexed: 07/16/2024] Open
Abstract
We present a four-component relativistic unitary coupled cluster method for atoms and molecules. We have used commutator-based non-perturbative approximation using the "Bernoulli expansion" to derive an approximation to the relativistic unitary coupled cluster method. The performance of the full quadratic unitary coupled-cluster singles and doubles method (qUCCSD), as well as a perturbative approximation variant (UCC3), has been reported for both energies and properties. It can be seen that both methods give results comparable to those of the standard relativistic coupled cluster method. The qUCCSD method shows better agreement with experimental results due to the better inclusion of relaxation effects. The relativistic UCC3 and qUCCSD methods can simulate the spin-forbidden transition with easy access to transition properties. A natural spinor-based scheme to reduce the computational cost of relativistic UCC3 and qUCCSD methods has been discussed.
Collapse
Affiliation(s)
- Kamal Majee
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Sudipta Chakraborty
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Tamoghna Mukhopadhyay
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Malaya K Nayak
- Theoretical Chemistry Section, Bhabha Atomic Research Centre, Trombay, Mumbai 400085, India
- Homi Bhabha National Institute, BARC Training School Complex, Anushakti Nagar, Mumbai 400094, India
| | - Achintya Kumar Dutta
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| |
Collapse
|
2
|
Windom ZW, Claudino D, Bartlett RJ. A new "gold standard": Perturbative triples corrections in unitary coupled cluster theory and prospects for quantum computing. J Chem Phys 2024; 160:214113. [PMID: 38832905 DOI: 10.1063/5.0202567] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2024] [Accepted: 05/16/2024] [Indexed: 06/06/2024] Open
Abstract
A major difficulty in quantum simulation is the adequate treatment of a large collection of entangled particles, synonymous with electron correlation in electronic structure theory, with coupled cluster (CC) theory being the leading framework for dealing with this problem. Augmenting computationally affordable low-rank approximations in CC theory with a perturbative account of higher-rank excitations is a tractable and effective way of accounting for the missing electron correlation in those approximations. This is perhaps best exemplified by the "gold standard" CCSD(T) method, which bolsters the baseline CCSD with the effects of triple excitations using considerations from many-body perturbation theory (MBPT). Despite this established success, such a synergy between MBPT and the unitary analog of CC theory (UCC) has not been explored. In this work, we propose a similar approach wherein converged UCCSD amplitudes are leveraged to evaluate energy corrections associated with triple excitations, leading to the UCCSD[T] method. In terms of quantum computing, this correction represents an entirely classical post-processing step that improves the energy estimate by accounting for triple excitation effects without necessitating new quantum algorithm developments or increasing demand for quantum resources. The rationale behind this choice is shown to be rigorous by studying the properties of finite-order UCC energy functionals, and our efforts do not support the addition of the fifth-order contributions as in the (T) correction. We assess the performance of these approaches on a collection of small molecules and demonstrate the benefits of harnessing the inherent synergy between MBPT and UCC theories.
Collapse
Affiliation(s)
- Zachary W Windom
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
- Quantum Information Science Section, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Daniel Claudino
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Rodney J Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
3
|
Stein F, Hutter J. Massively parallel implementation of gradients within the random phase approximation: Application to the polymorphs of benzene. J Chem Phys 2024; 160:024120. [PMID: 38214385 DOI: 10.1063/5.0180704] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2023] [Accepted: 12/15/2023] [Indexed: 01/13/2024] Open
Abstract
The Random-Phase approximation (RPA) provides an appealing framework for semi-local density functional theory. In its Resolution-of-the-Identity (RI) approach, it is a very accurate and more cost-effective method than most other wavefunction-based correlation methods. For widespread applications, efficient implementations of nuclear gradients for structure optimizations and data sampling of machine learning approaches are required. We report a well scaling implementation of RI-RPA nuclear gradients on massively parallel computers. The approach is applied to two polymorphs of the benzene crystal obtaining very good cohesive and relative energies. Different correction and extrapolation schemes are investigated for further improvement of the results and estimations of error bars.
Collapse
Affiliation(s)
- Frederick Stein
- Center for Advanced Systems Understanding (CASUS), Helmholtz-Zentrum Dresden, Rossendorf (HZDR), Untermarkt 20, 02826 Görlitz, Germany
| | - Jürg Hutter
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, 8057 Zurich, Switzerland
| |
Collapse
|
4
|
Kim Y, Krylov AI. Two Algorithms for Excited-State Quantum Solvers: Theory and Application to EOM-UCCSD. J Phys Chem A 2023; 127:6552-6566. [PMID: 37505075 DOI: 10.1021/acs.jpca.3c02480] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/29/2023]
Abstract
Near-term quantum devices promise to revolutionize quantum chemistry, but simulations using the current noisy intermediate-scale quantum (NISQ) devices are not practical due to their high susceptibility to errors. This motivated the design of NISQ algorithms leveraging classical and quantum resources. While several developments have shown promising results for ground-state simulations, extending the algorithms to excited states remains challenging. This paper presents two cost-efficient excited-state algorithms inspired by the classical Davidson algorithm. We implemented the Davidson method into the quantum self-consistent equation-of-motion unitary coupled-cluster (q-sc-EOM-UCC) excited-state method adapted for quantum hardware. The circuit strategies for generating desired excited states are discussed, implemented, and tested. We demonstrate the performance and accuracy of the proposed algorithms (q-sc-EOM-UCC/Davidson and its variational variant) by simulations of H2, H4, LiH, and H2O molecules. Similar to the classical Davidson scheme, q-sc-EOM-UCC/Davidson algorithms are capable of targeting a small number of excited states of the desired character.
Collapse
Affiliation(s)
- Yongbin Kim
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, United States
| | - Anna I Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, United States
| |
Collapse
|
5
|
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
|
6
|
Wang M, Fang WH, Li C. Assessment of State-Averaged Driven Similarity Renormalization Group on Vertical Excitation Energies: Optimal Flow Parameters and Applications to Nucleobases. J Chem Theory Comput 2023; 19:122-136. [PMID: 36534617 DOI: 10.1021/acs.jctc.2c00966] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We present a comprehensive excited-state benchmark for the state-averaged (SA) driven similarity renormalization group (DSRG) [Li, C.; Evangelista, F. A. J. Chem. Phys. 2018, 148, 124106]. Following the QUEST database [Véril, M.; Scemama, A.; Caffarel, M.; Lipparini, F.; Boggio-Pasqua, M.; Jacquemin, D.; Loos, P.-F. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2021, 11, e1517], 280 vertical transition energies of 35 medium-sized molecules are computed using the SA-DSRG derived second- and third-order perturbation theories (PT2/PT3) along with a nonperturbative approach [sq-LDSRG(2)]. Comparing to the theoretical best estimates, the optimal flow parameter is found to be 0.35 and 2.0 Eh-2 for SA-DSRG-PT2 and SA-DSRG-PT3, respectively. For SA-sq-LDSRG(2), a flow parameter of 1.5 Eh-2 provides converged equations without compromising the accuracy. We then assess the accuracy of the SA-DSRG hierarchy using these parameters. The SA-DSRG-PT2 scheme outperforms the level-shifted CASPT2 by 0.10 eV in mean absolute error (MAE), yet this accuracy is slightly inferior than that of CASPT2 with the ionization-potential-electron-affinity shift. Both SA-DSRG-PT3 and SA-sq-LDSRG(2) yield a MAE of 0.10 eV, which is comparable to that of CASPT3 (0.09 eV). Finally, we compute vertical excitation energies of several low-lying singlet states of nucleobases. The SA-sq-LDSRG(2) approach provides highly accurate results for π → π* excitations, while n → π* transitions are better described by SA-DSRG-PT3.
Collapse
Affiliation(s)
- Meng Wang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Chenyang Li
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|