While the main thrust of quantum computing research in materials science is to accurately measure the classically intractable electron correlation effects due to Coulomb repulsion, designing optimal quantum algorithms for simpler problems with well-understood solutions is a useful tactic to advance our quantum "toolbox". With this in mind, we consider the quantum calculation of a periodic system's single-electron band structure over a path through reciprocal space. Previous efforts have used the Variational Quantum Eigensolver algorithm to solve the energy of each band, which involves numerically optimizing the parameters of a variational quantum circuit to minimize a cost function, constructed as the expectation value of a Hamiltonian operator. Traditionally, a unique Hamiltonian operator is constructed for each k-point, so that many cost functions, each with their own parameter space, must be optimized to generate a single band. Similarly, calculating higher bands than the first has traditionally involved modifying the cost function with additional overlap terms to ensure higher-energy eigenstates are orthogonal to those of lower bands. In this paper, we adopt a direct space approach, using a novel hybrid first/second-quantized qubit mapping which allows us to construct a single Hamiltonian, and a single cost-function, suitable for solving the entire electronic band structure. In contrast to previous approaches, the k-point and the band index are selected by additional parameters in our quantum circuit, rather than through modifications to the cost function. The result is a technically and conceptually simpler approach to band structure calculations on a quantum computer. Moreover, we expect that the tools developed herein will motivate new strategies for tackling highly-correlated materials beyond the grasp of classical computing.
While the main thrust of quantum computing research in materials science is to accurately measure the classically intractable electron correlation effects due to Coulomb repulsion, designing optimal quantum algorithms for simpler problems with well-understood solutions is a useful tactic to advance our quantum "toolbox". With this in mind, we consider the quantum calculation of a periodic system's single-electron band structure over a path through reciprocal space. Previous efforts have used the Variational Quantum Eigensolver algorithm to solve the energy of each band, which involves numerically optimizing the parameters of a variational quantum circuit to minimize a cost function, constructed as the expectation value of a Hamiltonian operator. Traditionally, a unique Hamiltonian operator is constructed for each k-point, so that many cost functions, each with their own parameter space, must be optimized to generate a single band. Similarly, calculating higher bands than the first has traditionally involved modifying the cost function with additional overlap terms to ensure higher-energy eigenstates are orthogonal to those of lower bands. In this paper, we adopt a direct space approach, using a novel hybrid first/second-quantized qubit mapping which allows us to construct a single Hamiltonian, and a single cost-function, suitable for solving the entire electronic band structure. In contrast to previous approaches, the k-point and the band index are selected by additional parameters in our quantum circuit, rather than through modifications to the cost function. The result is a technically and conceptually simpler approach to band structure calculations on a quantum computer. Moreover, we expect that the tools developed herein will motivate new strategies for tackling highly-correlated materials beyond the grasp of classical computing.
Abstract While the main thrust of quantum computing research in materials science is to accurately measure the classically intractable electron correlation effects due to Coulomb repulsion, designing optimal quantum algorithms for simpler problems with well-understood solutions is a useful tactic to advance our quantum ``toolbox''. With this in mind, we consider the quantum calculation of a periodic system's single-electron band structure over a path through reciprocal space. Previous efforts have used the Variational Quantum Eigensolver algorithm to solve the energy of each band, which involves numerically optimizing the parameters of a variational quantum circuit to minimize a cost function, constructed as the expectation value of a Hamiltonian operator. Traditionally, a unique Hamiltonian operator is constructed for each k-point, so that many cost functions, each with their own parameter space, must be optimized to generate a single band. Similarly, calculating higher bands than the first has traditionally involved modifying the cost function with additional overlap terms to ensure higher-energy eigenstates are orthogonal to those of lower bands. In this paper, we adopt a direct space approach, using a novel hybrid first/second-quantized qubit mapping which allows us to construct a single Hamiltonian, and a single cost-function, suitable for solving the entire band-structure. In contrast to previous approaches, the k-point and the band index are selected by additional parameters in our quantum circuit, rather than through modifications to the cost function. The result is a technically and conceptually simpler approach to band structure calculations on a quantum computer. Moreover, we expect that the tools developed herein will motivate new strategies for tackling highly-correlated materials beyond the grasp of classical computing.
Development of quantum architectures during the last decade has inspired hybrid classical-quantum algorithms in physics and quantum chemistry that promise simulations of fermionic systems beyond the capability of modern classical computers, even before the era of quantum computing fully arrives. Strong research efforts have been recently made to obtain minimal depth quantum circuits which could accurately represent chemical systems. Here, we show that unprecedented methods used in quantum chemistry, designed to simulate molecules on quantum processors, can be extended to calculate properties of periodic solids. In particular, we present minimal depth circuits implementing the variational quantum eigensolver algorithm and successfully use it to compute the band structure of silicon on a quantum machine for the first time. We are convinced that the presented quantum experiments performed on cloud-based platforms will stimulate more intense studies towards scalable electronic structure computation of advanced quantum materials.
The variational quantum eigensolver (VQE) typically optimizes variational parameters in a quantum circuit to prepare eigenstates for a quantum system. Its application to many problems may involve a group of Hamiltonians, e.g., a Hamiltonian of a molecule is a function of the nuclear configurations. In this paper we incorporate derivatives of the Hamiltonian into the VQE and develop some hybrid quantum-classical algorithms, which explore both Hamiltonian and wave-function spaces for solving quantum chemistry problems more efficiently. For geometry optimization, we propose a mutual gradient-descent algorithm that updates the parameters of the Hamiltonian and wave function alternately, which can give a rapid convergence towards equilibrium structures of molecules as demonstrated in numerical simulations. Moreover, to speed up the calculation of the energy potential surface, we establish differential equations that govern how the optimized variational parameters of the wave function change with the intrinsic parameters of the Hamiltonian. Our study suggests a direction for a hybrid quantum-classical algorithm for solving quantum systems more efficiently by considering spaces of both the Hamiltonian and the wave function.
Accurately predicting the properties of molecules and materials from first principles requires solving the electronic structure problem. However, the computational cost of solving the electronic structure problem scales prohibitively in the number of atoms, and the usual approximate methods designed to reduce the overhead are often inaccurate for strongly correlated systems. For this reason, quantum algorithms have been developed as a promising alternative. In this work, we develop algebraic tools to improve the efficiency of three quantum algorithms for the electronic structure problem, addressing key bottlenecks for both near-term and fault-tolerant quantum computers. The first is the variational quantum eigensolver (VQE), which suffers from a large measurement overhead needed to estimate the Hamiltonian expectation value. To this end, we introduce a general class of measurable fragments of the electronic Hamiltonian, which enable partitioning schemes that significantly reduce the sampling cost for energy estimation. The second are quantum subspace methods (QSMs). We introduce a new QSM, the quantum seniority-based subspace expansion algorithm, which is a scalable alternative to VQE that alleviates VQE's state preparation, measurement, and optimization challenges while reducing qubit requirements. The third is quantum phase estimation (QPE), which requires an encoding of the electronic Hamiltonian directly on the quantum computer, the cost of which scales indirectly with the Hamiltonian spectral range. To address this, we develop linear programming approaches to implement the Block Invariant Symmetry Shift protocol, a method which reduces the spectral range of the Hamiltonian, improving its block-encoding cost for QPE. Collectively, these methods demonstrate that one can exploit the Hamiltonian's algebraic structure to mitigate computational costs, accelerating the path toward practical quantum advantage for electronic structure.
Abstract Variational quantum algorithms such as the variational quantum eigensolver (VQE) and variational quantum deflation (VQD) have emerged as leading candidates for computing electronic structure on near-term quantum hardware. For periodic solids, however, their reach has been limited not by the quantum algorithm itself but by the cost of constructing a faithful second-quantized Hamiltonian for each material, typically via density-functional theory (DFT) followed by Wannierization or by hand-fitted empirical tight-binding parameters. SlaKoNet, a recently developed method, significantly reduces the computational cost of electronic Hamiltonian construction with a small accuracy tradeoff. It integrates modern deep-learning methods and the Slater-Koster tight-binding formalism to fit thousands of onsite and distance-dependent hopping and overlap tight-binding parameters across 65 elements of the periodic table. After fitting, these parameters can be used to deterministically construct an electronic Hamiltonian for any crystal made up of any of the 65 elements in the training set, and the crystal's band structure can subsequently be computed. In this work, we combine a SlaKoNet model trained on the JARVIS-TBmBJ dataset with a Qiskit-based variational quantum deflation algorithm, replacing the traditional and costly Hamiltonian construction step with a universal neural Hamiltonian generator. Within the 65-element set, the resulting SlaKoNet-VQD workflow is structure-agnostic, fully differentiable, fast, and suited to high-throughput screening of band structures. We benchmark the workflow on silicon, recovering the full eight-band band structure along the standard high-symmetry k-path with a mean absolute deviation of 1.78 meV from exact diagonalization on a 3-qubit statevector simulator. We further compute the full band structure for five additional conventional superconductors, Al, Ta, Nb, V, and ZrN, achieving a mean absolute deviation within approximately 1.8 meV, and demonstrate full execution of SlaKoNet-VQD on IBM Quantum hardware for a representative k-point ground-state calculation on aluminum with a mean absolute error of roughly 0.37 eV. Beyond the single-particle level, we promote the same SlaKoNet Hamiltonian to a correlated Hubbard model and solve it with dynamical mean-field theory, recovering the weakening of correlations across the group-5 metals and strong quasiparticle renormalization in the cuprate La₂CuO₄, and we identify the resulting impurity problem as the natural target for a quantum solver. The pipeline opens a path to high-throughput benchmarking of variational quantum algorithms across the periodic table and to gradient-based co-optimization of ansatz and Hamiltonian for materials discovery. An interactive web app is available at https://atomgpt.org/quantum.
High Resolution Image Download MS PowerPoint Slide We propose a cost-efficient measurement scheme of the variational quantum eigensolver (VQE) for atomistic simulations of electronic structures based on a tight-binding (TB) theory. Leveraging the lattice geometry of a material domain, the sparse TB Hamiltonian is constructed in a bottom-up manner and is represented as a linear combination of the standard-basis (SB) operators. The cost function is evaluated with an extended version of the Bell measurement circuit that can simultaneously measure multiple SB operators and therefore reduces the number of circuits required by the evaluation process. The proposed VQE scheme is applied to find band gap energies of metal-halide-perovskite supercells that have finite dimensions with closed boundaries and are described with an sp 3 TB model. Experimental results confirm that the proposed scheme gives solutions that follow well the accurate ones but, more importantly, has computing efficiency that is obviously superior to the commutativity-based Pauli grouping methods. Extending the application scope of VQE to three-dimensional confined atomic structures, this work can serve as a practical guideline for handling TB simulations and, more generally, for calculations of sparse Hermitian matrices in the noise-intermediate-scale quantum devices.