Implementing time-evolution operators in shallow quantum circuits is important for quantum simulations. The standard method of Trotterization requires a large number of gates to achieve practical accuracy. Variational Quantum Simulation (VQS) is an algorithm that calculates the time evolution of a quantum state and can be executed with shallower circuits than Trotterization. However, the operator obtained by VQS evolves only a fixed initial state and is not the time evolution operator itself. In this paper, we propose Variational Quantum Operator Simulation (VQOS), a method to realize time evolution operators in shallow quantum circuits. This method is based on the variational principle for operators and does not require the implementation of the desired Trotter decomposition of the time evolution operator. We performed numerical simulations of the VQOS algorithm and successfully implemented the time evolution operator for closed systems in a quantum circuit that is up to 5 times shallower than the Trotterization. By providing a more practical way to implement time evolution operators, VQOS increases the applicability of near-term quantum computers.
Implementing time-evolution operators in shallow quantum circuits is important for quantum simulations. The standard method of Trotterization requires a large number of gates to achieve practical accuracy. Variational Quantum Simulation (VQS) is an algorithm that calculates the time evolution of a quantum state and can be executed with shallower circuits than Trotterization. However, the operator obtained by VQS evolves only a fixed initial state and is not the time evolution operator itself. In this paper, we propose Variational Quantum Operator Simulation (VQOS), a method to realize time evolution operators in shallow quantum circuits. This method is based on the variational principle for operators and does not require the implementation of the desired Trotter decomposition of the time evolution operator. We performed numerical simulations of the VQOS algorithm and successfully implemented the time evolution operator for closed systems in a quantum circuit that is up to 5 times shallower than the Trotterization. By providing a more practical way to implement time evolution operators, VQOS increases the applicability of near-term quantum computers.
Dynamical simulation of quantum many-body systems is a central task in quantum chemistry and requires efficient wavefunction representations. Although tensor-network and neural-network quantum states have achieved considerable success in ground-state calculations, entanglement growth hinders their application to quantum dynamics. Quantum computing may offer a route to quantum advantage, but algorithms such as Trotterization generally require fault-tolerant quantum computers, whereas near-term hardware supports only limited circuit sizes. Variational quantum simulation (VQS) instead represents the evolving state with a single shallow, fixed-size parameterized quantum circuit (PQC). Here, we introduce nonorthogonal variational quantum simulation (NOVQS), which applies linear combinations of parameterized quantum states to real- and imaginary-time evolutions. We design a shallow, hardware-friendly ansatz tailored to second-quantized electronic-structure Hamiltonians, together with resource-efficient protocols for measuring the matrices and vectors in the parameter equations of motion. Error analysis and resource estimation are also provided. Numerical simulations of hydrogen chains and the nitrogen molecule demonstrate that a collection of shallow, or even single-layer, PQCs can match or outperform a much deeper PQC in VQS. In particular, we identify a trade-off between circuit number and depth. NOVQS therefore provides a flexible route to enhancing wavefunction expressivity under circuit-depth constraints, making it a promising framework for quantum dynamics simulations on near-term quantum processors.
Abstract We present an empirical analysis of the scaling of the minimal quantum circuit depth required for a variational quantum simulation (VQS) method to obtain a solution to the time evolution of a quantum system within a predefined error tolerance. In a comparison against a non-variational method based on Trotterized time evolution, we observe similar scaling behaviours of the depth requirements of VQS and Trotter when the size of the system and the simulated time grow linearly. Results are also put into perspective by discussing the corresponding classical complexity required for VQS. Our results allow us to identify a possible advantage region for VQS over Trotterization.
Simulating time evolution of quantum systems is one of the most promising applications of quantum computing and also appears as a subroutine in many applications such as Green's function methods. In the current era of NISQ machines we assess the state of algorithms for simulating time dynamics with limited resources. We propose the Jaynes-Cummings model and extensions to it as useful toy models to investigate time evolution algorithms on near-term quantum computers. Using these simple models, direct Trotterisation of the time evolution operator produces deep circuits, requiring coherence times out of reach on current NISQ hardware. Therefore we test two alternative responses to this problem: variational compilation of the time evolution operator, and variational quantum simulation of the wavefunction ansatz. We demonstrate numerically to what extent these methods are successful in time evolving this system. The costs in terms of circuit depth and number of measurements are compared quantitatively, along with other drawbacks and advantages of each method. We find that computational requirements for both methods make them suitable for performing time evolution simulations of our models on NISQ hardware. Our results also indicate that variational quantum compilation produces more accurate results than variational quantum simulation, at the cost of a larger number of measurements.
Implementing time evolution operators on quantum circuits is important for quantum simulation. However, the standard way, Trotterization, requires a huge numbers of gates to achieve desirable accuracy. Here, we propose a local variational quantum compilation (LVQC) algorithm, which allows to accurately and efficiently compile a time evolution operators on a large-scale quantum system by the optimization with smaller-size quantum systems. LVQC utilizes a subsystem cost function, which approximates the fidelity of the whole circuit, defined for each subsystem as large as approximate causal cones brought by the Lieb-Robinson (LR) bound. We rigorously derive its scaling property with respect to the subsystem size, and show that the optimization conducted on the subsystem size leads to the compilation of whole-system time evolution operators. As a result, LVQC runs with limited-size quantum computers or classical simulators that can handle such smaller quantum systems. For instance, finite-ranged and short-ranged interacting L-size systems can be compiled with O(L0)- or O(logL)-size quantum systems depending on observables of interest. Furthermore, since this formalism relies only on the LR bound, it can efficiently construct time evolution operators of various systems in generic dimension involving finite-, short-, and long-ranged interactions. We also numerically demonstrate the LVQC algorithm for one-dimensional systems. Employing classical simulation by time-evolving block decimation, we succeed in compressing the depth of a time evolution operators up to 40 qubits by the compilation for 20 qubits. LVQC not only provides classical protocols for designing large-scale quantum circuits, but also will shed light on applications of intermediate-scale quantum devices in implementing algorithms in larger-scale quantum devices.
We propose a variational alternative to the Trotter-Suzuki decomposition that provides greater control over errors while preserving the unitary structure of time evolution. The variational parameters in our ansatz are derived from a global action principle, where Euler-Lagrange equations govern their optimal dynamics. Unlike conventional wavefunction-based variational methods, our approach specifically targets the time evolution operation and this allows a single set of optimized parameters to be applied to any initial state for a fixed Hamiltonian avoiding costly optimization procedures. Our method outperforms the standard Trotter-Suzuki formulas, typically achieving higher accuracy than higher-order Suzuki schemes. This translates directly to quantum computing applications, where it enables the design of quantum circuits with fewer gates which reduces noise and improves precision. Although we focus on quantum dynamics, the method is broadly applicable to problems involving general time-evolution operators. Applied to various model Hamiltonians, our approach reduces errors by factors of 2 to 5 compared to Trotter-Suzuki decompositions, demonstrating its promise for accurate quantum simulation with improved efficiency. In certain cases, the variational ansatz achieves higher accuracy than more complex higher-order Suzuki formulas while reducing the gate count by nearly half within a single circuit layer. Furthermore, we derive approximate analytical expressions for the variational parameters up to cubic order in time, valid for generic Hamiltonians. These approximations enable long-time quantum simulations with improved accuracy over equivalent Suzuki decompositions, providing ready-to-use evolution formulas that match Suzuki's gate complexity while delivering better performance.
Simulation of quantum many-body systems is one of the most promising applications of quantum computers. It is crucial to efficiently implement the time-evolution operator as a quantum circuit to execute such simulations on near-term quantum computing devices with limited computational resources. However, standard approaches such as Trotterization sometimes require a deep quantum circuit, which is hard to implement on near-term quantum computers. Here, we propose a hybrid quantum-classical algorithm, called local subspace variational quantum compilation (LSVQC), for compiling the time-evolution operator of quantum many-body systems. The LSVQC performs a variational optimization to reproduce the action of the target time-evolution operator within a physically reasonable subspace. The optimization is performed for small local subsystems based on the Lieb-Robinson bound, which allows us to execute the cost function evaluation using small-scale quantum devices and/or classical computers. We demonstrate the validity of the LSVQC algorithm through numerical simulations of a simple spin-lattice model and an effective model of a parent compound of cuprate superconductors, Sr 2 CuO 3 , constructed by the downfolding method. It is shown that the LSVQC achieves a 95% reduction of the circuit depth for simulating quantum many-body dynamics compared to the Trotterization at best while maintaining the same computational accuracy. We also demonstrate that the restriction to a subspace leads to a substantial reduction of required resources and improved accuracy compared to the case of considering the entire Hilbert space. 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 about one order of magnitude larger than that for the Trotterization.