1
|
Singh H, Majumder S, Mishra S. Hückel molecular orbital theory on a quantum computer: A scalable system-agnostic variational implementation with compact encoding. J Chem Phys 2024; 160:194106. [PMID: 38767256 DOI: 10.1063/5.0210597] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2024] [Accepted: 05/01/2024] [Indexed: 05/22/2024] Open
Abstract
Hückel molecular orbital (HMO) theory provides a semi-empirical treatment of the electronic structure in conjugated π-electronic systems. A scalable system-agnostic execution of HMO theory on a quantum computer is reported here based on a variational quantum deflation (VQD) algorithm for excited state quantum simulation. A compact encoding scheme is proposed here that provides an exponential advantage over the direct mapping and allows for quantum simulation of the HMO model for systems with up to 2n conjugated centers with n qubits. The transformation of the Hückel Hamiltonian to qubit space is achieved by two different strategies: an iterative refinement transformation and the Frobenius-inner-product-based transformation. These methods are tested on a series of linear, cyclic, and hetero-nuclear conjugated π-electronic systems. The molecular orbital energy levels and wavefunctions from the quantum simulation are in excellent agreement with the exact classical results. However, the higher excited states of large systems are found to suffer from error accumulation in the VQD simulation. This is mitigated by formulating a variant of VQD that exploits the symmetry of the Hamiltonian. This strategy has been successfully demonstrated for the quantum simulation of C60 fullerene containing 680 Pauli strings encoded on six qubits. The methods developed in this work are easily adaptable to similar problems of different complexity in other fields of research.
Collapse
Affiliation(s)
- Harshdeep Singh
- Center of Computational and Data Sciences, Indian Institute of Technology, Kharagpur, India
| | - Sonjoy Majumder
- Department of Physics, Indian Institute of Technology, Kharagpur, India
| | - Sabyashachi Mishra
- Department of Chemistry, Indian Institute of Technology, Kharagpur, India
| |
Collapse
|
2
|
Prasad VK, Cheng F, Fekl U, Jacobsen HA. Applications of noisy quantum computing and quantum error mitigation to "adamantaneland": a benchmarking study for quantum chemistry. Phys Chem Chem Phys 2024; 26:4071-4082. [PMID: 38225897 DOI: 10.1039/d3cp03523a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2024]
Abstract
The field of quantum computing has the potential to transform quantum chemistry. The variational quantum eigensolver (VQE) algorithm has allowed quantum computing to be applied to chemical problems in the noisy intermediate-scale quantum (NISQ) era. Applications of VQE have generally focused on predicting absolute energies instead of chemical properties that are relative energy differences and that are most interesting to chemists studying a chemical problem. We address this shortcoming by constructing a molecular benchmark data set in this work containing isomers of C10H16 and carbocationic rearrangements of C10H15+, calculated at a high-level of theory. Using the data set, we compared noiseless VQE simulations to conventionally performed density functional and wavefunction theory-based methods to understand the quality of results. We also investigated the effectiveness of a quantum state tomography-based error mitigation technique in applications of VQE under noise (simulated and real). Our findings reveal that the use of quantum error mitigation is crucial in the NISQ era and advantageous to yield almost noiseless quality results.
Collapse
Affiliation(s)
- Viki Kumar Prasad
- The Edward S. Rogers Sr. Department of Electrical and Computer Engineering, University of Toronto, 10 Kings College Road, Toronto, Ontario, Canada, M5S 3G4. arno,
- Department of Chemical and Physical Sciences, University of Toronto Mississauga, 3359 Mississauga Road, Mississauga, Ontario, Canada, L5L 1C6.
| | - Freeman Cheng
- Department of Computer Science, University of Toronto, 40 St. George Street, Toronto, Ontario, Canada, M5S 2E4
| | - Ulrich Fekl
- Department of Chemical and Physical Sciences, University of Toronto Mississauga, 3359 Mississauga Road, Mississauga, Ontario, Canada, L5L 1C6.
| | - Hans-Arno Jacobsen
- The Edward S. Rogers Sr. Department of Electrical and Computer Engineering, University of Toronto, 10 Kings College Road, Toronto, Ontario, Canada, M5S 3G4. arno,
- Department of Computer Science, University of Toronto, 40 St. George Street, Toronto, Ontario, Canada, M5S 2E4
| |
Collapse
|
3
|
Wang Y, Xue S, Wang Y, Ding J, Shi W, Wang D, Liu Y, Liu Y, Fu X, Huang G, Huang A, Deng M, Wu J. Experimental quantum natural gradient optimization in photonics. OPTICS LETTERS 2023; 48:3745-3748. [PMID: 37450740 DOI: 10.1364/ol.494560] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/18/2023] [Accepted: 06/14/2023] [Indexed: 07/18/2023]
Abstract
Variational quantum algorithms (VQAs) combining the advantages of parameterized quantum circuits and classical optimizers, promise practical quantum applications in the noisy intermediate-scale quantum era. The performance of VQAs heavily depends on the optimization method. Compared with gradient-free and ordinary gradient descent methods, the quantum natural gradient (QNG), which mirrors the geometric structure of the parameter space, can achieve faster convergence and avoid local minima more easily, thereby reducing the cost of circuit executions. We utilized a fully programmable photonic chip to experimentally estimate the QNG in photonics for the first time, to the best of our knowledge. We obtained the dissociation curve of the He-H+ cation and achieved chemical accuracy, verifying the outperformance of QNG optimization on a photonic device. Our work opens up a vista of utilizing QNG in photonics to implement practical near-term quantum applications.
Collapse
|
4
|
Xie QX, Song Y, Zhao Y. Power of the Sine Hamiltonian Operator for Estimating the Eigenstate Energies on Quantum Computers. J Chem Theory Comput 2022; 18:7586-7602. [PMID: 36449783 DOI: 10.1021/acs.jctc.2c00759] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/02/2022]
Abstract
Quantum computers have been shown to have tremendous potential in solving difficult problems in quantum chemistry. In this paper, we propose a new classical-quantum hybrid method, named as power of sine Hamiltonian operator (PSHO), to evaluate the eigenvalues of a given Hamiltonian (Ĥ). In PSHO, for any reference state |φ0⟩, the normalized energy of the sinn(H^τ)|φ0⟩ state can be determined. With the increase of the power, the initial reference state can converge to the eigenstate with the largest |sin(Eiτ)| value in the coefficients of the expansion of |φ0⟩, and the normalized energy of the sinn(H^τ)|φ0⟩ state converges to Ei. The ground- and excited-state energies of a Hamiltonian can be determined by taking different τ values. The performance of the PSHO method is demonstrated by numerical calculations of the H4 and LiH molecules. Compared with the current popular variational quantum eigensolver method, PSHO does not need to design the ansatz circuits and avoids the complex nonlinear optimization problems. PSHO has great application potential in near-term quantum devices.
Collapse
Affiliation(s)
- Qing-Xing Xie
- The Institute of Technological Sciences, Wuhan University, Wuhan430072, People's Republic of China
| | - Yi Song
- The Institute of Technological Sciences, Wuhan University, Wuhan430072, People's Republic of China
| | - Yan Zhao
- The Institute of Technological Sciences, Wuhan University, Wuhan430072, People's Republic of China
| |
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
|
Xie QX, Liu S, Zhao Y. Orthogonal State Reduction Variational Eigensolver for the Excited-State Calculations on Quantum Computers. J Chem Theory Comput 2022; 18:3737-3746. [PMID: 35621354 DOI: 10.1021/acs.jctc.2c00159] [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/30/2022]
Abstract
Variational quantum eigensolver (VQE) is a promising method for ground-state calculations on current noisy intermediate-scale quantum computers. However, the research progress of excited-state calculations on quantum computers is relatively slow. In order to extend the framework of VQE for excited-state calculations, we propose a new algorithm, orthogonal state reduction variational eigensolver (OSRVE), to determine the energies of excited states. Theoretical derivations prove that the optimized state in the OSRVE method can ensure the energy minimum and orthogonality constraint simultaneously, and OSRVE is also applicable for the degenerate state. The performance of OSRVE is demonstrated by numerical calculations of the H4 and H2O molecules. Compared with other excited-state calculation algorithms, OSRVE has obvious advantages in calculating lower-order excited states. This work extends the VQE algorithm to excited-state calculations, and OSRVE can be implemented on near-term quantum computers.
Collapse
Affiliation(s)
- Qing-Xing Xie
- The Institute of Technological Sciences, Wuhan University, Wuhan 430072, People's Republic of China
| | - Sheng Liu
- The Institute of Technological Sciences, Wuhan University, Wuhan 430072, People's Republic of China
| | - Yan Zhao
- The Institute of Technological Sciences, Wuhan University, Wuhan 430072, People's Republic of China
| |
Collapse
|
7
|
Yoshida R, Lötstedt E, Yamanouchi K. Quantum computing of Hückel molecular orbitals of π-electron systems. J Chem Phys 2022; 156:184117. [PMID: 35568559 DOI: 10.1063/5.0086489] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In order to demonstrate an applicability of quantum computing to fundamental electronic structure problems of molecules, we describe the Hückel Hamiltonian matrix in terms of quantum gates and obtain the orbital energies of fundamental π-electron molecules (C2H4, C3H4, C4H4, C4H6, and C6H6) using a superconducting-qubit-type quantum computer (ibm_kawasaki) with a post-selection error mitigation method. We show that the orbital energies are obtained with sufficiently high accuracy and small uncertainties and that characteristic features of the electronic structure of the π-electron molecules can be extracted by quantum computing in a straightforward manner.
Collapse
Affiliation(s)
- Ryuhei Yoshida
- Department of Chemistry, School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan
| | - Erik Lötstedt
- Department of Chemistry, School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan
| | - Kaoru Yamanouchi
- Department of Chemistry, School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan
| |
Collapse
|
8
|
Sugisaki K, Kato T, Minato Y, Okuwaki K, Mochizuki Y. Variational quantum eigensolver simulations with the multireference unitary coupled cluster ansatz: a case study of the C2v quasi-reaction pathway of beryllium insertion to H 2 molecule. Phys Chem Chem Phys 2022; 24:8439-8452. [DOI: 10.1039/d1cp04318h] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Variational quantum eigensolver (VQE)-based quantum chemical calculations have been extensively studied as a computational model using noisy intermediate-scale quantum devices. VQE uses a parametrized quantum circuit defined through an “ansatz”...
Collapse
|
9
|
Sugisaki K, Sakai C, Toyota K, Sato K, Shiomi D, Takui T. Bayesian phase difference estimation: a general quantum algorithm for the direct calculation of energy gaps. Phys Chem Chem Phys 2021; 23:20152-20162. [PMID: 34551045 DOI: 10.1039/d1cp03156b] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Abstract
Quantum computers can perform full configuration interaction (full-CI) calculations by utilising the quantum phase estimation (QPE) algorithms including Bayesian phase estimation (BPE) and iterative quantum phase estimation (IQPE). In these quantum algorithms, the time evolution of wave functions for atoms and molecules is simulated conditionally with an ancillary qubit as the control, which make implementation to real quantum devices difficult. Also, most of the problems in chemistry discuss energy differences between two electronic states rather than total energies themselves, and thus direct calculations of energy gaps are promising for future applications of quantum computers to real chemistry problems. In the race of finding efficient quantum algorithms to solve quantum chemistry problems, we test a Bayesian phase difference estimation (BPDE) algorithm, which is a general algorithm to calculate the difference of two eigenphases of unitary operators in the several cases of the direct calculations of energy gaps between two electronic states on quantum computers, including vertical ionisation energies, singlet-triplet energy gaps, and vertical excitation energies. In the BPDE algorithm, state preparation is carried out conditionally on the ancillary qubit, and the time evolution of the wave functions in superposition of two electronic states are executed unconditionally. Based on our test, we conclude that BPDE is capable of computing the energy gap with an accuracy similar to BPE without controlled-time evolution simulations and with the smaller number of iterations in Bayesian optimisations.
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 Centre, University Administration Division, Osaka City University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka 558-8585, Japan
| |
Collapse
|
10
|
Sugisaki K, Toyota K, Sato K, Shiomi D, Takui T. A quantum algorithm for spin chemistry: a Bayesian exchange coupling parameter calculator with broken-symmetry wave functions. Chem Sci 2020; 12:2121-2132. [PMID: 34163976 PMCID: PMC8179312 DOI: 10.1039/d0sc04847j] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2020] [Accepted: 12/13/2020] [Indexed: 01/03/2023] Open
Abstract
The Heisenberg exchange coupling parameter J (H = -2J S i · S j ) characterises the isotropic magnetic interaction between unpaired electrons, and it is one of the most important spin Hamiltonian parameters of multi-spin open shell systems. The J value is related to the energy difference between high-spin and low-spin states, and thus computing the energies of individual spin states are necessary to obtain the J values from quantum chemical calculations. Here, we propose a quantum algorithm, B̲ayesian ex̲change coupling parameter calculator with b̲roken-symmetry wave functions (BxB), which is capable of computing the J value directly, without calculating the energies of individual spin states. The BxB algorithm is composed of the quantum simulations of the time evolution of a broken-symmetry wave function under the Hamiltonian with an additional term j S 2, the wave function overlap estimation with the SWAP test, and Bayesian optimisation of the parameter j. Numerical quantum circuit simulations for H2 under a covalent bond dissociation, C, O, Si, NH, OH+, CH2, NF, O2, and triple bond dissociated N2 molecule revealed that the BxB can compute the J value within 1 kcal mol-1 of errors with less computational costs than conventional quantum phase estimation-based approaches.
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
| | - 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 Centre, University Administration Division, Osaka City University 3-3-138 Sugimoto, Sumiyoshi-ku Osaka 558-8585 Japan
| |
Collapse
|
11
|
Teplukhin A, Kendrick BK, Tretiak S, Dub PA. Electronic structure with direct diagonalization on a D-wave quantum annealer. Sci Rep 2020; 10:20753. [PMID: 33247201 PMCID: PMC7695747 DOI: 10.1038/s41598-020-77315-4] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2020] [Accepted: 11/06/2020] [Indexed: 11/22/2022] Open
Abstract
Quantum chemistry is regarded to be one of the first disciplines that will be revolutionized by quantum computing. Although universal quantum computers of practical scale may be years away, various approaches are currently being pursued to solve quantum chemistry problems on near-term gate-based quantum computers and quantum annealers by developing the appropriate algorithm and software base. This work implements the general Quantum Annealer Eigensolver (QAE) algorithm to solve the molecular electronic Hamiltonian eigenvalue-eigenvector problem on a D-Wave 2000Q quantum annealer. The approach is based on the matrix formulation, efficiently uses qubit resources based on a power-of-two encoding scheme and is hardware-dominant relying on only one classically optimized parameter. We demonstrate the use of D-Wave hardware for obtaining ground and excited electronic states across a variety of small molecular systems. The approach can be adapted for use by a vast majority of electronic structure methods currently implemented in conventional quantum-chemical packages. The results of this work will encourage further development of software such as qbsolv which has promising applications in emerging quantum information processing hardware and has expectation to address large and complex optimization problems intractable for classical computers.
Collapse
Affiliation(s)
- Alexander Teplukhin
- Theoretical Division (T-1, MS B221), Los Alamos National Laboratory, Los Alamos, NM, 87545, USA
| | - Brian K Kendrick
- Theoretical Division (T-1, MS B221), Los Alamos National Laboratory, Los Alamos, NM, 87545, USA.
| | - Sergei Tretiak
- Theoretical Division (T-1, MS B221), Los Alamos National Laboratory, Los Alamos, NM, 87545, USA
| | - Pavel A Dub
- Chemistry Division (C-IIAC, MS K558), Los Alamos National Laboratory, Los Alamos, NM, 87545, USA.
| |
Collapse
|
12
|
Wu B, Cheng B, Jia F, Zhang J, Yung MH, Sun X. Speedup in classical simulation of Gaussian boson sampling. Sci Bull (Beijing) 2020; 65:832-841. [PMID: 36659202 DOI: 10.1016/j.scib.2020.02.012] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2019] [Revised: 11/28/2019] [Accepted: 12/29/2019] [Indexed: 01/21/2023]
Abstract
Gaussian boson sampling is an alternative model for demonstrating quantum computational supremacy, where squeezed states are injected into every input mode, instead of applying single photons as in the case of standard boson sampling. Here by analyzing numerically the computational costs, we establish a lower bound for achieving quantum computational supremacy for a class of Gaussian boson-sampling problems. Specifically, we propose a more efficient method for calculating the transition probabilities, leading to a significant reduction of the simulation costs. Particularly, our numerical results indicate that one can simulate up to 18 photons for Gaussian boson sampling at the output subspace on a normal laptop, 20 photons on a commercial workstation with 256 cores, and about 30 photons for supercomputers. These numbers are significantly smaller than those in standard boson sampling, suggesting that Gaussian boson sampling could be experimentally-friendly for demonstrating quantum computational supremacy.
Collapse
Affiliation(s)
- Bujiao Wu
- Institute of Computing Technology, Chinese Academy of Sciences, Beijing 100190, China; University of Chinese Academy of Sciences, Beijing 100049, China
| | - Bin Cheng
- Institute for Quantum Science and Engineering, Department of Physics, Southern University of Science and Technology, Shenzhen 518055, China
| | - Fei Jia
- Central Research Institute, Huawei Technologies, Shenzhen 518129, China
| | - Jialin Zhang
- Institute of Computing Technology, Chinese Academy of Sciences, Beijing 100190, China; University of Chinese Academy of Sciences, Beijing 100049, China
| | - Man-Hong Yung
- Institute for Quantum Science and Engineering, Department of Physics, Southern University of Science and Technology, Shenzhen 518055, China; Shenzhen Key Laboratory of Quantum Science and Engineering, Southern University of Science and Technology, Shenzhen 518055, China; Central Research Institute, Huawei Technologies, Shenzhen 518129, China.
| | - Xiaoming Sun
- Institute of Computing Technology, Chinese Academy of Sciences, Beijing 100190, China; University of Chinese Academy of Sciences, Beijing 100049, China; CAS Center of Excellence in Topological Quantum Computation, Beijing 100190, China.
| |
Collapse
|
13
|
Sugisaki K, Toyota K, Sato K, Shiomi D, Takui T. A probabilistic spin annihilation method for quantum chemical calculations on quantum computers. Phys Chem Chem Phys 2020; 22:20990-20994. [PMID: 32940301 DOI: 10.1039/d0cp03745a] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
A probabilistic spin annihilation method based on the quantum phase estimation algorithm is presented for quantum chemical calculations on quantum computers. This approach can eliminate more than one spin component from the spin contaminated wave functions by single operation. Comparison with the spin annihilation operation on classical computers is given.
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.
| | | | | | | | | |
Collapse
|
14
|
Zhang T, Li C, Evangelista FA. Improving the Efficiency of the Multireference Driven Similarity Renormalization Group via Sequential Transformation, Density Fitting, and the Noninteracting Virtual Orbital Approximation. J Chem Theory Comput 2019; 15:4399-4414. [PMID: 31268704 DOI: 10.1021/acs.jctc.9b00353] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
This study examines several techniques to improve the efficiency of the linearized multireference driven similarity renormalization group truncated to one- and two-body operators [MR-LDSRG(2)]. We propose a sequential MR-LDSRG(2) [sq-MR-LDSRG(2)] scheme, in which one-body substitutions are folded exactly into the Hamiltonian. This new approach is combined with density fitting (DF) to reduce the storage cost of two-electron integrals. To further avoid storage of large four-index intermediates, we propose a noninteracting virtual orbital (NIVO) approximation of the Baker-Campbell-Hausdorff series that neglects commutators terms with three and four virtual indices. The NIVO approximation reduces the computational prefactor of the MR-LDSRG(2), bringing it closer to that of coupled cluster with singles and doubles (CCSD). We test the effect of the DF and NIVO approximations on the MR-LDSRG(2) and sq-MR-LDSRG(2) methods by computing properties of eight diatomic molecules. The diatomic constants obtained by DF-sq-MR-LDSRG(2)+NIVO are found to be as accurate as those from the original MR-LDSRG(2) and coupled cluster theory with singles, doubles, and perturbative triples. Finally, we demonstrate that the DF-sq-MR-LDSRG(2)+NIVO scheme can be applied to chemical systems with more than 550 basis functions by computing the automerization energy of cyclobutadiene with a quintuple-ζ basis set. The predicted automerization energy is found to be similar to the value computed with Mukherjee's state-specific multireference coupled cluster theory with singles and doubles.
Collapse
Affiliation(s)
- Tianyuan Zhang
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation , Emory University , Atlanta , Georgia 30322 , United States
| | - Chenyang Li
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation , Emory University , Atlanta , Georgia 30322 , United States
| | - Francesco A Evangelista
- Department of Chemistry and Cherry L. Emerson Center for Scientific Computation , Emory University , Atlanta , Georgia 30322 , United States
| |
Collapse
|
15
|
Yung MH, Gao X, Huh J. Universal bound on sampling bosons in linear optics and its computational implications. Natl Sci Rev 2019; 6:719-729. [PMID: 34691927 PMCID: PMC8291458 DOI: 10.1093/nsr/nwz048] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2018] [Revised: 03/24/2019] [Accepted: 04/06/2019] [Indexed: 12/12/2022] Open
Abstract
In linear optics, photons are scattered in a network through passive optical elements including beam splitters and phase shifters, leading to many intriguing applications in physics, such as Mach-Zehnder interferometry, the Hong-Ou-Mandel effect, and tests of fundamental quantum mechanics. Here we present the fundamental limit in the transition amplitudes of bosons, applicable to all physical linear optical networks. Apart from boson sampling, this transition bound results in many other interesting applications, including behaviors of Bose-Einstein condensates (BEC) in optical networks, counterparts of Hong-Ou-Mandel effects for multiple photons, and approximating permanents of matrices. In addition, this general bound implies the existence of a polynomial-time randomized algorithm for estimating the transition amplitudes of bosons, which represents a solution to an open problem raised by Aaronson and Hance (Quantum Inf Comput 2012; 14: 541-59). Consequently, this bound implies that computational decision problems encoded in linear optics, prepared and detected in the Fock basis, can be solved efficiently by classical computers within additive errors. Furthermore, our result also leads to a classical sampling algorithm that can be applied to calculate the many-body wave functions and the S-matrix of bosonic particles.
Collapse
Affiliation(s)
- Man-Hong Yung
- Shenzhen Institute for Quantum Science and Engineering and Department of Physics, Southern University of Science and Technology, Shenzhen 518055, China
- Shenzhen Key Laboratory of Quantum Science and Engineering, Southern University of Science and Technology, Shenzhen 518055, China
- Central Research Institute, Huawei Technologies, Shenzhen 518129, China
| | - Xun Gao
- Center for Quantum Information, Institute for Interdisciplinary Information Sciences, Tsinghua University, Beijing 100084, China
| | - Joonsuk Huh
- Department of Chemistry, Sungkyunkwan University, Suwon 440-746, Korea
| |
Collapse
|
16
|
Sugisaki K, Nakazawa S, Toyota K, Sato K, Shiomi D, Takui T. Quantum chemistry on quantum computers: quantum simulations of the time evolution of wave functions under the S 2 operator and determination of the spin quantum number S. Phys Chem Chem Phys 2019; 21:15356-15361. [PMID: 31270515 DOI: 10.1039/c9cp02546d] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Quantum computers have an enormous impact on quantum chemical calculations. Approaches to calculate the energies of atoms and molecules on quantum computers by utilizing quantum phase estimation (QPE) and the variational quantum eigensolver (VQE) have been well documented, and dozens of methodological improvements to decrease computational costs and to mitigate errors have been reported until recently. However, the possible methodological implementation of observables on quantum computers such as calculating the spin quantum numbers of arbitrary wave functions, which is a crucial issue in quantum chemistry, has been discussed less. Here, we propose a quantum circuit to simulate the time evolution of wave functions under an S2 operator, exp(-iS2t)|Ψ, and integrate it into the QPE circuit enabling us to determine the spin quantum number of the arbitrary wave functions. We demonstrate that the spin quantum numbers of up to three spins can be determined by only one qubit measurement in QPE.
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.
| | - Shigeaki Nakazawa
- 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. and Research Support Department/University Research Administrator Center, Unviersity Administration Division, Osaka City University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka 558-8585, Japan
| |
Collapse
|