1
|
Sugisaki K. Projective Measurement-Based Quantum Phase Difference Estimation Algorithm for the Direct Computation of Eigenenergy Differences on a Quantum Computer. J Chem Theory Comput 2023; 19:7617-7625. [PMID: 37874368 PMCID: PMC10653105 DOI: 10.1021/acs.jctc.3c00784] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2023] [Indexed: 10/25/2023]
Abstract
Quantum computers are capable of calculating the energy difference of two electronic states using the quantum phase difference estimation (QPDE) algorithm. The Bayesian inference-based implementations for the QPDE have been reported so far, but in this approach, the quality of the calculated energy difference depends on the input wave functions being used. Here, we report the inverse quantum Fourier transformation-based QPDE with Na of ancillary qubits, which allows us to compute the difference of eigenenergies based on the single-shot projective measurement. As proof-of-concept demonstrations, we report numerical experiments for the singlet-triplet energy difference of the hydrogen molecule and the vertical excitation energies of halogen-substituted methylenes (CHF, CHCl, CF2, CFCl, and CCl2) and formaldehyde (HCHO).
Collapse
Affiliation(s)
- Kenji Sugisaki
- Graduate
School of Science and Technology, Keio University, 7-1 Shinkawasaki, Saiwai-ku, Kawasaki, Kanagawa 212-0032, Japan
- Quantum
Computing Center, Keio University, 3-14-1 Hiyoshi,
Kohoku-ku Yokohama, Kanagawa 223-8522, Japan
- Centre
for Quantum Engineering, Research and Education
TCG Centres for Research and Education in Science and Technology, Sector V, Salt Lake, Kolkata 700091, India
| |
Collapse
|
2
|
Shirai S, Iwakiri H, Kanno K, Horiba T, Omiya K, Hirai H, Koh S. Computational Analysis of Chemical Reactions Using a Variational Quantum Eigensolver Algorithm without Specifying Spin Multiplicity. ACS OMEGA 2023; 8:19917-19925. [PMID: 37305284 PMCID: PMC10249088 DOI: 10.1021/acsomega.3c01875] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 03/20/2023] [Accepted: 05/10/2023] [Indexed: 06/13/2023]
Abstract
The analysis of a chemical reaction along the ground-state potential energy surface in conjunction with an unknown spin state is challenging because electronic states must be separately computed several times using different spin multiplicities to find the lowest energy state. However, in principle, the ground state could be obtained with just a single calculation using a quantum computer without specifying the spin multiplicity in advance. In the present work, ground-state potential energy curves for PtCO were calculated as a proof-of-concept using a variational quantum eigensolver (VQE) algorithm. This system exhibits a singlet-triplet crossover as a consequence of the interaction between Pt and CO. VQE calculations using a statevector simulator were found to converge to a singlet state in the bonding region, while a triplet state was obtained at the dissociation limit. Calculations performed using an actual quantum device provided potential energies within ±2 kcal/mol of the simulated energies after error mitigation techniques were adopted. The spin multiplicities in the bonding and dissociation regions could be clearly distinguished even in the case of a small number of shots. The results of this study suggest that quantum computing can be a powerful tool for the analysis of the chemical reactions of systems for which the spin multiplicity of the ground state and variations in this parameter are not known in advance.
Collapse
Affiliation(s)
- Soichi Shirai
- Toyota
Central Research and Development Laboratories, Inc., 41-1 Yokomichi, Nagakute, Aichi 480-1192, Japan
| | - Hokuto Iwakiri
- QunaSys
Inc., Aqua Hakusan Building
9F, 1-13-7 Hakusan, Bunkyo, Tokyo 113-0001, Japan
| | - Keita Kanno
- QunaSys
Inc., Aqua Hakusan Building
9F, 1-13-7 Hakusan, Bunkyo, Tokyo 113-0001, Japan
| | - Takahiro Horiba
- Toyota
Central Research and Development Laboratories, Inc., 41-1 Yokomichi, Nagakute, Aichi 480-1192, Japan
| | - Keita Omiya
- QunaSys
Inc., Aqua Hakusan Building
9F, 1-13-7 Hakusan, Bunkyo, Tokyo 113-0001, Japan
| | - Hirotoshi Hirai
- Toyota
Central Research and Development Laboratories, Inc., 41-1 Yokomichi, Nagakute, Aichi 480-1192, Japan
| | - Sho Koh
- QunaSys
Inc., Aqua Hakusan Building
9F, 1-13-7 Hakusan, Bunkyo, Tokyo 113-0001, Japan
| |
Collapse
|
3
|
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
|
4
|
Sugisaki K, Wakimoto H, Toyota K, Sato K, Shiomi D, Takui T. Quantum Algorithm for Numerical Energy Gradient Calculations at the Full Configuration Interaction Level of Theory. J Phys Chem Lett 2022; 13:11105-11111. [PMID: 36444985 PMCID: PMC9743205 DOI: 10.1021/acs.jpclett.2c02737] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2022] [Accepted: 11/07/2022] [Indexed: 06/16/2023]
Abstract
A Bayesian phase difference estimation (BPDE) algorithm allows us to compute the energy gap of two electronic states of a given Hamiltonian directly by utilizing the quantum superposition of their wave functions. Here we report an extension of the BPDE algorithm to the direct calculation of the energy difference of two molecular geometries. We apply the BPDE algorithm for the calculation of numerical energy gradients based on the two-point finite-difference method, enabling us to execute geometry optimization of one-dimensional molecules at the full-CI level on a quantum computer. Results of numerical quantum circuit simulations of the geometry optimization of the H2 molecule with the STO-3G and 6-31G basis sets, the LiH and BeH2 molecules at the full-CI/STO-3G level, and the N2 molecule at the CASCI(6e,6o)/6-311G* level are given.
Collapse
Affiliation(s)
- Kenji Sugisaki
- Department
of Chemistry, Graduate School of Science, Osaka Metropolitan University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka558-8585, Japan
- JSTPRESTO, 4-1-8 Honcho, Kawaguchi, Saitama, 332-0012, Japan
- Centre
for Quantum Engineering, Research and Education (CQuERE), TCG Centres for Research and Education in Science
and Technology (TCG CREST), Sector V,
Salt Lake, Kolkata700091, India
| | - Hiroyuki Wakimoto
- Department
of Chemistry, Graduate School of Science, Osaka Metropolitan University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka558-8585, Japan
| | - Kazuo Toyota
- Department
of Chemistry, Graduate School of Science, Osaka Metropolitan University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka558-8585, Japan
| | - Kazunobu Sato
- Department
of Chemistry, Graduate School of Science, Osaka Metropolitan University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka558-8585, Japan
| | - Daisuke Shiomi
- Department
of Chemistry, Graduate School of Science, Osaka Metropolitan University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka558-8585, Japan
| | - Takeji Takui
- Department
of Chemistry, Graduate School of Science, Osaka Metropolitan University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka558-8585, Japan
- Research
Support Department/University Research Administrator Center, University
Administration Division, Osaka Metropolitan
University, 3-3-138 Sugimoto,
Sumiyoshi-ku, Osaka558-8585, Japan
| |
Collapse
|
5
|
Sugisaki K, Toyota K, Sato K, Shiomi D, Takui T. Adiabatic state preparation of correlated wave functions with nonlinear scheduling functions and broken-symmetry wave functions. Commun Chem 2022; 5:84. [PMID: 36698020 PMCID: PMC9814591 DOI: 10.1038/s42004-022-00701-8] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2022] [Accepted: 06/30/2022] [Indexed: 01/28/2023] Open
Abstract
Adiabatic state preparation (ASP) can generate the correlated wave function by simulating the time evolution of wave function under the time-dependent Hamiltonian that interpolates the Fock operator and the full electronic Hamiltonian. However, ASP is inherently unsuitable for studying strongly correlated systems, and furthermore practical computational conditions for ASP are unknown. In quest for the suitable computational conditions for practical applications of ASP, we performed numerical simulations of ASP in the potential energy curves of N2, BeH2, and in the C2v quasi-reaction pathway of the Be atom insertion to the H2 molecule, examining the effect of nonlinear scheduling functions and the ASP with broken-symmetry wave functions with the S2 operator as the penalty term, contributing to practical applications of quantum computing to quantum chemistry. Eventually, computational guidelines to generate the correlated wave functions having the square overlap with the complete-active space self-consistent field wave function close to unity are discussed.
Collapse
Affiliation(s)
- Kenji Sugisaki
- grid.261445.00000 0001 1009 6411Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka 558-8585 Japan ,grid.419082.60000 0004 1754 9200JST PRESTO, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012 Japan ,grid.510650.7Centre for Quantum Engineering, Research and Education (CQuERE), TCG Centres for Research and Education in Science and Technology (TCG CREST), Sector V, Salt Lake, Kolkata, 700091 India
| | - Kazuo Toyota
- grid.261445.00000 0001 1009 6411Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka 558-8585 Japan
| | - Kazunobu Sato
- grid.261445.00000 0001 1009 6411Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka 558-8585 Japan
| | - Daisuke Shiomi
- grid.261445.00000 0001 1009 6411Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka 558-8585 Japan
| | - Takeji Takui
- grid.261445.00000 0001 1009 6411Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka 558-8585 Japan ,grid.261445.00000 0001 1009 6411Research Support Department/University Research Administrator Center, University Administration Division, Osaka City University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka, 558-8585 Japan
| |
Collapse
|
6
|
Shirai S, Horiba T, Hirai H. Calculation of Core-Excited and Core-Ionized States Using Variational Quantum Deflation Method and Applications to Photocatalyst Modeling. ACS OMEGA 2022; 7:10840-10853. [PMID: 35382310 PMCID: PMC8973155 DOI: 10.1021/acsomega.2c01053] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/21/2022] [Accepted: 03/03/2022] [Indexed: 06/14/2023]
Abstract
The possibility of performing quantum-chemical calculations using quantum computers has attracted much interest. Variational quantum deflation (VQD) is a quantum-classical hybrid algorithm for the calculation of excited states with noisy intermediate-scale quantum devices. Although the validity of this method has been demonstrated, there have been few practical applications, primarily because of the uncertain effect of calculation conditions on the results. In the present study, calculations of the core-excited and core-ionized states for common molecules based on the VQD method were examined using a classical computer, focusing on the effects of the weighting coefficients applied in the penalty terms of the cost function. Adopting a simplified procedure for estimating the weighting coefficients based on molecular orbital levels allowed these core-level states to be successfully calculated. The O 1s core-ionized state for a water molecule was calculated with various weighting coefficients, and the resulting ansatz states were systematically examined. The application of this technique to functional materials was demonstrated by calculating the core-level states for titanium dioxide (TiO2) and nitrogen-doped TiO2 models. The results demonstrate that VQD calculations employing an appropriate cost function can be applied to the analysis of functional materials in conjunction with an experimental approach.
Collapse
|
7
|
Sugisaki K, Sakai C, Toyota K, Sato K, Shiomi D, Takui T. Quantum Algorithm for Full Configuration Interaction Calculations without Controlled Time Evolutions. J Phys Chem Lett 2021; 12:11085-11089. [PMID: 34749498 DOI: 10.1021/acs.jpclett.1c03214] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
A quantum phase estimation algorithm allows us to perform full configuration interaction (full-CI) calculations on quantum computers with polynomial costs against the system size under study, but it requires quantum simulation of the time evolution of the wave function conditional on an ancillary qubit, which makes the algorithm implementation on real quantum devices difficult. Here, we discuss an application of the Bayesian phase difference estimation algorithm that is free from controlled time evolution operations to the full-CI calculations.
Collapse
Affiliation(s)
- Kenji Sugisaki
- Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka 558-8585, Japan
- JST PRESTO, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
- Centre for Quantum Engineering, Research and Education (CQuERE), TCG Centres for Research and Education in Science and Technology (TCG CREST), 16th Floor, Omega, BIPL Building, Blocks EP & GP, Sector V, Salt Lake, Kolkata 700091, India
| | - Chikako Sakai
- Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka 558-8585, Japan
| | - Kazuo Toyota
- Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka 558-8585, Japan
| | - Kazunobu Sato
- Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka 558-8585, Japan
| | - Daisuke Shiomi
- Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka 558-8585, Japan
| | - Takeji Takui
- Department of Chemistry and Molecular Materials Science, Graduate School of Science, Osaka City University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka 558-8585, Japan
- Research Support Department, University Research Administrator Center, University Administration Division, Osaka City University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka 558-8585, Japan
| |
Collapse
|