Abstract Classical computation of electronic properties in large-scale materials remains challenging. Quantum computation has the potential to offer advantages in memory footprint and computational scaling. However, general and viable quantum algorithms for simulating large-scale materials are still limited. We propose and implement random-state quantum algorithms to calculate electronic-structure properties of real materials. Using a random state circuit on a small number of qubits, we employ real-time evolution with first-order Trotter decomposition and Hadamard test to obtain electronic density of states, and we develop a modified quantum phase estimation algorithm to calculate real-space local density of states via direct quantum measurements. Furthermore, we validate these algorithms by numerically computing the density of states and spatial distributions of electronic states in graphene, twisted bilayer graphene quasicrystals, and fractal lattices, covering system sizes from hundreds to thousands of atoms. Our results manifest that the random-state quantum algorithms provide a general and qubit-efficient route to scalable simulations of electronic properties in large-scale periodic and aperiodic materials.
Abstract Classical computation of electronic properties in large-scale materials remains challenging. Quantum computation has the potential to offer advantages in memory footprint and computational scaling. However, general and viable quantum algorithms for simulating large-scale materials are still limited. We propose and implement random-state quantum algorithms to calculate electronic-structure properties of real materials. Using a random state circuit on a small number of qubits, we employ real-time evolution with first-order Trotter decomposition and Hadamard test to obtain electronic density of states, and we develop a modified quantum phase estimation algorithm to calculate real-space local density of states via direct quantum measurements. Furthermore, we validate these algorithms by numerically computing the density of states and spatial distributions of electronic states in graphene, twisted bilayer graphene quasicrystals, and fractal lattices, covering system sizes from hundreds to thousands of atoms. Our results manifest that the random-state quantum algorithms provide a general and qubit-efficient route to scalable simulations of electronic properties in large-scale periodic and aperiodic materials.
With quantum computers on the horizon, it is essential to understand what problems they can be used to solve. One potential application is simulating quantum systems, in particular, systems with many interacting electrons. However, there is a disparity between the computational resources required to study these problems and the capabilities of modern experimental quantum devices. In this thesis, we make progress toward closing this gap by reducing the cost of quantum simulation of quantum chemistry on a number of fronts. In the first half, we develop three improved quantum simulation algorithms by approximating time evolution under the Hamiltonian using a truncated Taylor series rather than the Trotter-Suzuki decomposition. We do this first for the molecular electronic structure Hamiltonian in two representations, obtaining asymptotic improvements in terms of both the number of orbitals and the required precision. We also apply the truncated Taylor series paradigm to quantum simulation in real space, where we consider simulation of general wave functions discretized on a uniform grid. In the second half, by working with Hamiltonians tailored specifically to quantum algorithms, we demonstrate that Trotter steps of the electronic structure Hamiltonian can be simulated in exactly N depth and with only linear connectivity, where N is the number of spin-orbitals. This resolves the problem of non-locality in fermion-to-qubit mappings, and has wide-ranging applications for variational and phase estimation-based quantum simulation of quantum chemistry. We further study this algorithm on an error-corrected quantum computer using the surface code. Benchmarking it alongside other product formula approaches, we show that quantum simulation of scientifically important instances of the electronic structure problem that are classically intractable can be performed with only a few hundred thousand physical qubits and fewer than one billion costly fault-tolerant gates. Finally, we study the effect of using randomized Hamiltonians in iterative phase estimation, and how this allows us to reduce the costs of simulation as well as to exploit knowledge of the ground state in how we design our simulation algorithm.
Abstract Determining the ground and excited state properties of materials is considered one of the most promising applications of quantum computers. On near-term hardware, the limiting constraint on such simulations is the requisite circuit depths and qubit numbers, which currently lie well beyond near-term capabilities. Here we develop a quantum algorithm which reduces the estimated cost of material simulations. For example, we obtain a circuit depth improvement by up to 6 orders of magnitude for a Trotter layer of time-dynamics simulation in the transition-metal oxide SrVO3compared with the best previous quantum algorithms. We achieve this by introducing a collection of connected techniques, including highly localised and physically compact representations of materials Hamiltonians in the Wannier basis, a hybrid fermion-to-qubit mapping, and an efficient circuit compiler. Combined together, these methods leverage locality of materials Hamiltonians and result in a design that generates quantum circuits with depth independent of the system’s size. Although the requisite resources for the quantum simulation of materials are still beyond current hardware, our results show that realistic simulation of specific properties may be feasible without necessarily requiring fully scalable, fault-tolerant quantum computers, providing quantum algorithm design incorporates deeper understanding of the target materials and applications.
At the heart of quantum chemistry and materials science lies the critical task of estimating ground-state properties. We present a quantum algorithm for this task by quantizing the density-functional theory (DFT). A key aspect of imple-menting DFT faithfully is the requirement for self-consistent calculations, which involve repeated di-agonalizations of the Hamiltonian. This procedure, however, creates a significant bottleneck, as a classical algorithm generally demands a computational complexity that grows cubically with the number of electrons, restricting the scalability of DFT for tackling large-scale problems that involve complex chemical environments and microstructures. This article presents the first quantum algorithm that has provided substantial speedup for the ground state computation, by improving the complexity to one with a linear scaling with the number of atoms. The algorithm leverages the exponential speedup by the quantum singular value transformation to generate a quantum circuit to encode the density-matrix, followed by an efficient estimation method for the output electron density, which constitutes a simple hybrid approach for achieving self-consistency. Moreover, the algorithm produces the ground state Hamiltonian, from which the ground state energy and band structures can be efficiently computed. The proposed framework is accompanied by a rigorous error analysis that establishes the convergence and quantifies various sources of error and the overall computational complexity. The combination of effi-ciency and precision opens new avenues for exploring large-scale physical systems.
We devise a quantum-circuit algorithm to solve the ground state and ground energy of artificial graphene. The algorithm implements a Trotterized adiabatic evolution from a purely tight-binding Hamiltonian to one including kinetic, spin-orbit, and Coulomb terms. The initial state is obtained efficiently using Gaussian-state preparation, while the readout of the ground energy is organized into seventeen sets of measurements, irrespective of the size of the problem. The total depth of the corresponding quantum circuit scales polynomially with the size of the system. A full simulation of the algorithm is performed and ground energies are obtained for lattices with up to four hexagons. Our results are benchmarked with exact diagonalization for systems with one and two hexagons. For larger systems we use the exact state vector and approximate matrix product state simulation techniques. The latter allows us to systematically trade off precision with memory and therefore to tackle larger systems. We analyze adiabatic and Trotterization errors, providing estimates for optimal periods and time discretizations given a finite accuracy. In the case of large systems we also study approximation errors.
Density-functional theory (DFT) has revolutionized computer simulations in chemistry and material science. A faithful implementation of the theory requires self-consistent calculations. However, this effort involves repeatedly diagonalizing the Hamiltonian, for which a classical algorithm typically requires a computational complexity that scales cubically with respect to the number of electrons. This limits DFT's applicability to large-scale problems with complex chemical environments and microstructures. This article presents a quantum algorithm that has a linear scaling with respect to the number of atoms, which is much smaller than the number of electrons. Our algorithm leverages the quantum singular value transformation (QSVT) to generate a quantum circuit to encode the density-matrix, and an estimation method for computing the output electron density. In addition, we present a randomized block coordinate fixed-point method to accelerate the self-consistent field calculations by reducing the number of components of the electron density that needs to be estimated. The proposed framework is accompanied by a rigorous error analysis that quantifies the function approximation error, the statistical fluctuation, and the iteration complexity. In particular, the analysis of our self-consistent iterations takes into account the measurement noise from the quantum circuit. These advancements offer a promising avenue for tackling large-scale DFT problems, enabling simulations of complex systems that were previously computationally infeasible.
Determining properties of molecules and materials is one of the premier applications of quantum computing. A major question in the thesis is: how might we use long- and near-term quantum computers to solve problems of chemical interest? First, we give an overview of the standard model in electronic structure theory. Electronic structure methods have been improved significantly in the last decades; however, it is still a non-trivial task to fit experimental results for complex molecular systems. Second, we introduce quantum computation and quantum algorithms as an alternative device to solve the electronic structure problem. In the remaining part of the thesis, we describe the research projects. The thesis covers two (rather different) research directions: software and hardware projects. For the software part, we contribute with two new algorithms: Phase-Estimation Interval Target Energy Readout (PHILTER) and Recursive Variational Series Estimation (RVSE), where both are designed to be run on a quantum digital computer and extract materials and molecular properties. The PHILTER algorithm is a modification of the Quantum Phase Estimation and is designed to compute excited state energies within a target energy-interval for cases where good approximations for the states in the target energy-interval are either unknown or hard to prepare. For the RVSE algorithm, we pose the question: how might we use imperfect near-term quantum computers to solve problems of practical value? In this work, we introduce a method to compute expectation values of general functions of a Hamiltonian with respect to an input state using low depth quantum circuits suitable for noisy intermediate-scale quantum computers. For the hardware part, we propose a new chemical quantum computer by using molecular electronics as a basis for quantum computing. In this work, we pose the question: can scattering of electrons in molecules produce quantum gates? We construct one-qubit gates using one-electron scattering in molecules, and two-qubit controlled-phase gates using electron-electron scattering along metallic leads. We show initial applications of the framework by illustrating one-qubit gates using the molecular electronic structure of molecular hydrogen as a baseline model.