Advances in quantum simulator technology are increasingly required because research on quantum algorithms is becoming more sophisticated and complex. State vector simulation utilizes CPU and memory resources in computing nodes exponentially with respect to the number of qubits; furthermore, in a variational quantum algorithm, the large number of repeated runs by classical optimization is also a heavy load. This problem has been addressed by preparing numerous computing nodes or simulation frameworks that operate effectively. This study aimed to accelerate quantum simulation using two newly proposed methods: to efficiently utilize limited computational resources by adjusting the ratio of the MPI and distributed processing parallelism corresponding to the target problem settings and to slim down the Hamiltonian by considering the effect of accuracy on the calculation result. Ground-state energy calculations of a fermionic model were performed using a variational quantum eigensolver (VQE) on an HPC cluster with up to 1024 FUJITSU Processor A64FX connected to each other by InfiniBand; the processor is also used on supercomputer Fugaku. We achieved 200 times higher speed over VQE simulations and demonstrated 32 qubits ground-state energy calculations in acceptable time. This result indicates that > 30 qubit state vector simulations can be realistically utilized for further research on variational quantum algorithms.
Advances in quantum simulator technology are increasingly required because research on quantum algorithms is becoming more sophisticated and complex. State vector simulation utilizes CPU and memory resources in computing nodes exponentially with respect to the number of qubits; furthermore, in a variational quantum algorithm, the large number of repeated runs by classical optimization is also a heavy load. This problem has been addressed by preparing numerous computing nodes or simulation frameworks that operate effectively. This study aimed to accelerate quantum simulation using two newly proposed methods: to efficiently utilize limited computational resources by adjusting the ratio of the MPI and distributed processing parallelism corresponding to the target problem settings and to slim down the Hamiltonian by considering the effect of accuracy on the calculation result. Ground-state energy calculations of a fermionic model were performed using a variational quantum eigensolver (VQE) on an HPC cluster with up to 1024 FUJITSU Processor A64FX connected to each other by InfiniBand; the processor is also used on supercomputer Fugaku. We achieved 200 times higher speed over VQE simulations and demonstrated 32 qubits ground-state energy calculations in acceptable time. This result indicates that > 30 qubit state vector simulations can be realistically utilized for further research on variational quantum algorithms.
Quantum computer simulators running on classical computers are essential for developing real quantum computers and emerging quantum applications. In particular, state vector simulators, which store a full state vector in memory and update it in every quantum operation, can be used to simulate an arbitrary form of quantum circuits, debug quantum applications, and validate future quantum computers. Since the time and space complexity grows exponentially with the number of qubits and easily exceeds the capability of a single machine, various distributed state vector simulators are developed. However, these are tuned for common CPUs and GPUs, not for ARM CPUs. Therefore, we develop a distributed state vector simulator, mpiQulacs, that is optimized for large-scale simulation on A64FX-based cluster systems. A64FX is a 48-core ARM-based CPU having high-bandwidth memory (HBM2) and also used in one of the world's top supercomputers, Fugaku. As computational optimization, we vectorize qubit operations with Scalable Vector Extension (SVE) to fully utilize its HBM2. To reduce the amount of inter-node communication for scalable simulation, mpiQulacs provides a fused-swap gate that collectively performs multiple swap operations as one operation and a quantum circuit optimizer that automatically adds fused-swap gates to a circuit. We evaluate weak and strong scaling of mpiQulacs with up to 39 qubits on our 512-node A64FX-based cluster system. The result shows that mpiQulacs achieves high-scalability on hundreds of nodes in benchmark and practical quantum circuits, and its performance is comparable with the existing well-known distributed state vector simulators.
Quantum computer simulators running on classical computers are essential for developing real quantum computers and emerging quantum applications. In particular, state vector simulators, which store a full state vector in memory and update it in every quantum operation, are available to simulate an arbitrary form of quantum circuits, debug quantum applications, and validate future quantum computers. However, the time and space complexity grows exponentially with the number of qubits and easily exceeds the capability of a single machine. Therefore, we develop a distributed state vector simulator, mpiQulacs, that is optimized for large-scale simulation on A64FX-based cluster systems. A64FX is an ARM-based CPU that is also equipped in the world's top Fugaku supercomputer. We evaluate weak and strong scaling of mpiQulacs with up to 36 qubits on a new 64-node A64FX-based cluster system named Todoroki. By comparing mpiQulacs with existing distributed state vector simulators, we show that mpiQulacs achieves the highest performance for large-scale simulation on tens of nodes while sustaining a nearly ideal scalability. Besides, we define a new metric, quantumB/Fratio, and use it to demonstrate that mpiQulacs running on Todoroki fits the requirements of distributed state vector simulation rather than the existing simulators running on general purpose CPU-based or GPU-based cluster systems.
Quantum simulators have been actively developed that faithfully simulate the quantum computer operations on a classical computer. Since the amount of memory and operations required by quantum simulators increases exponentially with the number of qubits, MPI (Message Passing Interface) and distributed processing, which are methods to effectively use a large number of computing resources, can be used to achieve higher computation speed. In this paper, we profile the performance of the variational quantum eigensolver (VQE) algorithm by applying it to a multi-node quantum simulator to calculate the energy values of carbon dioxide (CO2) molecule. We also discuss how the size of CI (configuration interaction) ansatz and the convergence threshold of the classical optimization change the energy values. Experimental results demonstrate that the computational time is reduced as the number of computation nodes increases. The speed-up ratio is increased as the number of qubits increases. When the qubit size is 30, a maximum of 12.2 times speed-up is achieved if the computing node size is increased from 1 to 128.
The Variational Quantum Eigensolver (VQE) is a hybrid quantum-classical algorithm for ground-state energies of molecular systems. We compute the potential energy surface of the hydrogen molecule (H2) at 100 bond lengths with PennyLane on the ERAU Vega cluster (AMD EPYC 9654 CPUs, 4 NVIDIA H100 GPUs) and compare a serial baseline, an Optax and Catalyst JIT version, a lightning.gpu version, and an MPI version. Early stopping, a larger learning rate, and warm starts reduce the iteration count from 30,000 to 7,033 and the runtime from 593.95 s to 143.80 s, while the time per iteration stays between 0.0198 s and 0.0204 s. At 4 qubits, lightning.gpu lowers the time per iteration to 0.0145 s. MPI runs with 2 to 32 ranks finish in 8.45 s to 5.04 s; the serial and MPI jobs ran under different core allocations, so we do not report their ratio as a parallel speedup. In a 4 to 26 qubit study with different CPU and GPU software stacks, the CPU/GPU time ratio ranges from 3.5 to 81.5. One H100 runs up to 29 qubits.
Exponential scaling presents a significant challenge in electronic structure calculations performed on classical computers. This paper explores how quantum computer algorithms can accurately represent quantum systems. Variational Quantum Eigensolver (VQE) algorithm is used to compute the ground state energy of hydrogen and helium sequences by implementing variational principle and quantum gates as trial wavefunction. This technique combines classical optimization with quantum computing calculations to simulate quantum systems on noisy and resource-limited computers. The resulting calculated energy is highly consistent to the corresponding exact values and Hartree-Fock calculations with a trend of when the number of atoms increases the calculated energy becomes more negative, leading to a decrease in the percentage error. Moreover, the convergence of the ground state energy of hydrogen and helium atoms was effectively optimized. The desired energy was reached, proven by adjusting the expectation value, and gradually achieving unity in state overlap. These findings demonstrate the VQE method's accuracy in calculating simple quantum systems and its scalability for larger atomic and molecular system, such as those in quantum chemistry and material science. However, challenges in quantum computer simulations, such as limited in qubit numbers and the presence of noise, require further advancements. Therefore, implementing a larger basis sets, advanced qubit mapping, specific chemistry ansatz, and flexible optimization techniques is one way to improve overall calculation.
The variational quantum eigensolver (or VQE), first developed by Peruzzo et al. (2014), has received significant attention from the research community in recent years. It uses the variational principle to compute the ground state energy of a Hamiltonian, a problem that is central to quantum chemistry and condensed matter physics. Conventional computing methods are constrained in their accuracy due to the computational limits facing exact modeling of the exponentially growing electronic wavefunction for these many-electron systems. The VQE may be used to model these complex wavefunctions in polynomial time, making it one of the most promising near-term applications for quantum computing. One important advantage is that variational algorithms have been shown to present some degree of resilience to the noise in the quantum hardware. Finding a path to navigate the relevant literature has rapidly become an overwhelming task, with many methods promising to improve different parts of the algorithm, but without clear descriptions of how the diverse parts fit together. The potential practical advantages of the algorithm are also widely discussed in the literature, but with varying conclusions. Despite strong theoretical underpinnings suggesting excellent scaling of individual VQE components, studies have pointed out that their various pre-factors could be too large to reach a quantum computing advantage over conventional methods. This review aims at disentangling the relevant literature to provide a comprehensive overview of the progress that has been made on the different parts of the algorithm, and to discuss future areas of research that are fundamental for the VQE to deliver on its promises. All the different components of the algorithm are reviewed in detail. These include the representation of Hamiltonians and wavefunctions on a quantum computer, the optimization process to find ground state energies, the post processing mitigation of quantum errors, and suggested best practices. We identify four main areas of future research: (1) optimal measurement schemes for reduction of circuit repetitions required; (2) large scale parallelization across many quantum computers; (3) ways to overcome the potential appearance of vanishing gradients in the optimization process for large systems, and how the number of iterations required for the optimization scales with system size; (4) the extent to which VQE suffers for quantum noise, and whether this noise can be mitigated in a tractable manner. The answers to these open research questions will determine the routes for the VQE to achieve quantum advantage as the quantum computing hardware scales up and as the noise levels are reduced.