Significant advances have been made in the study of quantum advantage both in theory and experiment, although these have mostly been limited to artificial setups. In this work, we extend the scope to address quantum advantage in tasks relevant to chemistry and physics. Specifically, we consider the unitary cluster Jastrow (UCJ) ansatz-a variant of the unitary coupled cluster ansatz, which is widely used to solve the electronic structure problem on quantum computers-to show that sampling from the output distributions of quantum circuits implementing the UCJ ansatz is likely to be classically hard. More specifically, we show that there exist UCJ circuits for which classical simulation of sampling cannot be performed in polynomial time, under a reasonable complexity-theoretical assumption that the polynomial hierarchy does not collapse. Our main contribution is to show that a class of UCJ circuits can be used to perform arbitrary instantaneous quantum polynomial-time (IQP) computations, which are already known to be classically hard to simulate under the same complexity assumption. As a side result, we also show that UCJ equipped with post-selection can generate the class post-BQP. Our demonstration, worst-case nonsimulatability of UCJ, would potentially imply quantum advantage in quantum algorithms for chemistry and physics using unitary coupled cluster type ansatzes, such as the variational quantum eigensolver and quantum-selected configuration interaction.

Quasiparticle band structures are fundamental for understanding strongly correlated electron systems. While solving these structures accurately on classical computers is challenging, quantum computing offers a promising alternative. Specifically, the quantum subspace expansion (QSE) method, combined with the variational quantum eigensolver (VQE), provides a quantum algorithm for calculating quasiparticle band structures. However, optimizing the variational parameters in VQE becomes increasingly difficult as the system size grows, due to device noise, statistical noise, and the barren plateau problem. To address these challenges, we propose a hybrid approach that combines QSE with the quantum-selected configuration interaction (QSCI) method for calculating quasiparticle band structures. QSCI may leverage the VQE ansatz as an input state but, unlike the standard VQE, it does not require full optimization of the variational parameters, making it more scalable for larger quantum systems. Based on this approach, we demonstrate the quantum computation of the quasiparticle band structure of a silicon using 16 qubits on an IBM quantum processor.

Quantum-selected configuration interaction (QSCI) is a novel quantum-classical hybrid algorithm for quantum chemistry calculations. This method identifies electron configurations having large weights for the target state using quantum devices and allows CI calculations to be performed with the selected configurations on classical computers. In principle, the QSCI algorithm can take advantage of the ability to handle large configuration spaces while reducing the negative effects of noise on the calculated values. At present, QSCI calculations are limited by qubit noise during the input state preparation and measurement process, restricting them to small active spaces. These limitations make it difficult to perform calculations with quantitative accuracy. The present study demonstrates a computational scheme based on multireference perturbation theory calculations on a classical computer, using the QSCI wavefunction as a reference. This method was applied to ground and excited state calculations for two typical aromatic molecules, naphthalene and tetracene. The incorporation of the perturbation treatment was found to provide improved accuracy. Extension of the reference space based on the QSCI-selected configurations as a means of further improvement was also investigated.

We study the Hamiltonian algebra of Jordan-Wigner-transformed interacting fermion models and its fast-forwardability. We prove that the dimension of the Hamiltonian algebra of the fermion model with single-site Coulomb interaction is bounded from below by the exponential function of the number of sites, and the circuit depth of the Cartan-based fast-forwarding method for such model also exhibits the same scaling. We apply this proposition to the Anderson impurity model and the Hubbard model and show that the dimension of the Hamiltonian algebra of these models scales exponentially with the number of sites. These behaviors of the Hamiltonian algebras imply that the qubit models obtained by the Jordan-Wigner transformation of these fermion models cannot be efficiently simulated using the Cartan-based fast-forwarding method.

We present a comprehensive study of the commute time kernel method via the effective resistance framework analyzing the quantum complexity of the originally classical approach. Our study reveals that while there is a trade-off between accuracy and computational complexity, significant improvements can be achieved in terms of runtime efficiency without substantially compromising on precision. Our investigation highlights a notable quantum speedup in calculating the kernel, which offers a quadratic improvement in time complexity over classical approaches in certain instances. In addition, we introduce methodical improvements over the original work on the commute time kernel and provide empirical evidence suggesting the potential reduction of kernel queries without significant impact on result accuracy. Benchmarking our method on several chemistry-based datasets: AIDS, NCL1, PTC−MR, MUTAG, PROTEINS - data points previously unexplored in existing literature, shows that while not always the most accurate, it excels in time efficiency. This makes it a compelling alternative for applications where computational speed is crucial. Our results highlight the balance between accuracy, computational complexity, and speedup offered by quantum computing, promoting further research into efficient algorithms for kernel methods and their applications in chemistry-based datasets.

Expectation value estimation is ubiquitous in quantum algorithms. The expectation value of a Hamiltonian, which is essential in various practical applications, is often estimated by measuring a large number of Pauli strings on quantum computers and performing classical post-processing. In the case of n-qubit molecular Hamiltonians in quantum chemistry calculations, it is necessary to evaluate O(n⁴) Pauli strings, requiring a large number of measurements for accurate estimation. To reduce the measurement cost, we assess an existing idea that uses two copies of an n-qubit quantum state of interest and coherently measures them in the Bell basis, which enables the simultaneous estimation of the absolute values of expectation values of all the n-qubit Pauli strings. We numerically investigate the efficiency of energy estimation for molecular Hamiltonians of up to 12 qubits. The results show that, when the target precision is no smaller than tens of milli-Hartree, this method requires fewer measurements than conventional sampling methods. This suggests that the method may be useful for many applications that rely on expectation value estimation of Hamiltonians and other observables as well when moderate precision is sufficient.

Quantum-selected configuration interaction (QSCI) utilizes an input quantum state on a quantum device to select important bases (electron configurations in quantum chemistry) which define a subspace where we diagonalize a target Hamiltonian, i.e., perform selected configuration interaction, on classical computers. Previous proposals for preparing a good input state, which is crucial for the quality of QSCI, based on optimization of quantum circuits may suffer from optimization difficulty and require many runs of the quantum device. Here we propose using a time-evolved state by the target Hamiltonian (for some initial state) as an input of QSCI. Our proposal is based on the intuition that the time evolution by the Hamiltonian creates electron excitations of various orders when applied to the initial state. We numerically investigate the accuracy of the energy obtained by the proposed method for quantum chemistry Hamiltonians describing electronic states of small molecules. Numerical results reveal that our method can yield sufficiently accurate ground-state energies for the investigated molecules. Our proposal provides a systematic and optimization-free method to prepare the input state of QSCI and could contribute to practical applications of quantum computers to quantum chemistry calculations.

In this review, we give an overview of the proposed applications in the early-FTQC (EFTQC) era. Starting from the error correction architecture for EFTQC device, we first review the recently developed space-time efficient analogue rotation (STAR) architecture, which is a partially fault-tolerant error correction architecture. Then, we review the requirements of an EFTQC algorithm. In particular, the class of ground state energy estimation (GSEE) algorithm known as the statistical phase estimation algorithm (SPE) is studied. We especially cast our attention on two SPE-type algorithms, the step-function filter-based variant by Lin and Tong (LT22) and Gaussian Filter. Based on the latter, we introduce the Gaussian Fitting algorithm, which uses an alternative post-processing procedure compared to the previous study. Finally, we systematically simulate the aforementioned algorithms and Variational Quantum Eigensolver (VQE) using the 1-uCJ ansatz with different shot counts. Most importantly, we perform noisy simulations based on the STAR architecture. We find that for estimating the ground state energy of the 4-qubit H2 Hamiltonian in the STO-3G basis, SPE becomes more advantageous over VQE when the physical error rate is sufficiently low.

Simulating thermal-equilibrium properties at finite temperature plays a crucial role in studying the nature of quantum many-body systems. In particular, implementing a finite-temperature simulation on a quantum computer is expected to overcome the difficulty in simulating large-sized systems, for which the quantum Monte-Carlo technique on a classical computer suffers from the sign problem in general. While several methods suitable for fault-tolerant quantum computing (FTQC) devices are expected to be useful in studying large-scale quantum many-body systems, those proposed so far involve a large number of ancilla qubits and a deep quantum circuit with many basic gates, making them unsuitable for the early-FTQC era, i.e., the early stage of the FTQC era, at which only a limited number of qubits and quantum gates are available. In this paper, we propose a method suitable for quantum devices in this early stage to calculate the thermal-equilibrium expectation value of an observable at finite temperature. Our proposal, named the Markov-chain Monte Carlo with sampled-pairs of unitaries (MCMC-SPU) algorithm, is based on sampling simple quantum circuits and generating the corresponding statistical ensembles, and overcomes the difficulties in the resource requirements and the decay in probability associated with postselection of measurement outcomes on ancilla qubits. We demonstrate the validity of our proposal by performing a numerical simulation of the MCMC-SPU algorithm on the one-dimensional transverse-field Ising model as an illustrative example.

Simulation of quantum many-body systems is a promising application of quantum computers. However, implementing the time-evolution operator as a quantum circuit efficiently on near-term devices with limited resources is challenging. Standard approaches like Trotterization often require deep circuits, making them impractical. This paper proposes a hybrid quantum-classical algorithm called Local Subspace Variational Quantum Compilation (LSVQC) for compiling the time-evolution operator. The LSVQC uses variational optimization to reproduce the action of the target time-evolution operator within a physically reasonable subspace. Optimization is performed on small local subsystems based on the Lieb-Robinson bound, allowing for cost function evaluation using small-scale quantum devices or classical computers. Numerical simulations on a spin-lattice model and an ab initio effective model of strongly correlated material Sr2CuO3 demonstrate the algorithm's effectiveness. It is shown that the LSVQC achieves a 95% reduction in circuit depth compared to Trotterization while maintaining accuracy. The subspace restriction also reduces resource requirements and improves accuracy. Furthermore, we estimate the gate count needed to execute the quantum simulations using the LSVQC on near-term quantum computing architectures in the noisy intermediate-scale or early fault-tolerant quantum computing era. Our estimation suggests that the acceptable physical gate error rate for the LSVQC can be significantly larger than for Trotterization.