The most widely used approach for simulating the dynamics of time-dependent Hamiltonians via quantum computation depends on the quantum-classical hybrid variational quantum time evolution algorithm, in which ordinary differential equations of the variational coefficients for determining time evolution are solved via classical simulations with a time discretization method. We here present a full-quantum approach, in which ordinary differential equations of the variational coefficients are transformed into static linear equations via the Chebyshev spectral discretization method and then solved via the quantum singular value transformation algorithm. Our full quantum algorithm avoids classical feedback, achieves exponential convergence for smooth Hamiltonians, and yields a quantum circuit depth that is independent of the number of time steps. We demonstrate two implementation strategies, with a global formulation designed for fault-tolerant architectures and a sequential formulation tailored to near-term devices, and validate the approach through numerical simulations of proton-hydrogen charge-transfer dynamics, a prototypical time-dependent quantum chemistry problem. This work establishes a systematic pathway from quantum-classical hybrid variational quantum algorithms to full-quantum solvers for general time-dependent Hamiltonians, particularly those whose dynamics admit compact variational descriptions, opening a route toward full quantum computational advantages in time-dependent simulations.
The most widely used approach for simulating the dynamics of time-dependent Hamiltonians via quantum computation depends on the quantum-classical hybrid variational quantum time evolution algorithm, in which ordinary differential equations of the variational coefficients for determining time evolution are solved via classical simulations with a time discretization method. We here present a full-quantum approach, in which ordinary differential equations of the variational coefficients are transformed into static linear equations via the Chebyshev spectral discretization method and then solved via the quantum singular value transformation algorithm. Our full quantum algorithm avoids classical feedback, achieves exponential convergence for smooth Hamiltonians, and yields a quantum circuit depth that is independent of the number of time steps. We demonstrate two implementation strategies, with a global formulation designed for fault-tolerant architectures and a sequential formulation tailored to near-term devices, and validate the approach through numerical simulations of proton-hydrogen charge-transfer dynamics, a prototypical time-dependent quantum chemistry problem. This work establishes a systematic pathway from quantum-classical hybrid variational quantum algorithms to full-quantum solvers for general time-dependent Hamiltonians, particularly those whose dynamics admit compact variational descriptions, opening a route toward full quantum computational advantages in time-dependent simulations.
The simulation of ion-atom collisions remains a formidable challenge due to the complex interplay between electronic and nuclear degrees of freedom. We present a hybrid quantum-classical computing framework for simulating time-dependent ion-atom collision dynamics, within which two variational quantum time evolution algorithms are implemented. To validate our framework, we simulate the charge-transfer dynamics and compute the corresponding cross sections for the proton-hydrogen collision system across an energy range 1--25 keV. Our results accurately reproduce the charge-transfer dynamics with high fidelity and exhibit very good agreement with available experimental and theoretical cross-section data across the entire energy range. These results highlight the accuracy and applicability of our hybrid quantum-classical framework for scattering cross-section calculations. Our work demonstrates an effective approach for mapping time-dependent many-body collision problems onto near-term quantum computing devices, and also provides promising directions for practical applications of universal quantum computing in the noisy intermediate-scale quantum era.
The simulation of many-body open quantum systems is key to solving numerous outstanding problems in physics, chemistry, material science, and in the development of quantum technologies. Near-term quantum computers may bring considerable advantage for the efficient simulation of their static and dynamical properties, thanks to hybrid quantum-classical variational algorithms to approximate the dynamics of the density matrix describing the quantum state in terms of an ensemble average. Here, a variational quantum algorithm is developed to simulate the real-time evolution of the density matrix governed by the Lindblad master equation, under the assumption that the quantum state has a bounded entropy along the dynamics, entailing a low-rank representation of its density matrix. The algorithm encodes each pure state of the statistical mixture as a parametrized quantum circuit, and the associated probabilities as additional variational parameters stored classically, thereby requiring a significantly lower number of qubits than algorithms where the full density matrix is encoded in the quantum memory. Two variational ansatze are proposed, and their effectiveness is assessed in the simulation of the dynamics of a 2D dissipative transverse field Ising model. The results underscore the algorithm's efficiency in simulating the dynamics of open quantum systems in the low-rank regime with limited quantum resources on a near-term quantum device.
Classical simulations of real-space quantum dynamics are challenging due to the exponential scaling of computational cost with system dimensions. Quantum computers offer the potential to simulate quantum dynamics with polynomial complexity; however, existing quantum algorithms based on the split-operator techniques require large-scale fault-tolerant quantum computers that remain elusive in the near future. Here, we present variational simulations of real-space quantum dynamics suitable for implementation in noisy intermediate-scale quantum (NISQ) devices. The Hamiltonian is first encoded onto qubits using a discrete variable representation and binary encoding scheme. We show that direct application of a real-time variational quantum algorithm based on the McLachlan's principle is inefficient as the measurement cost grows exponentially with the qubit number for general potential energy, and an extremely small time-step size is required to achieve accurate results. Motivated by the insights that many chemical dynamics occur in the low-energy subspace, we propose a subspace expansion method by projecting the total Hamiltonian, including the time-dependent driving field, onto the system low-energy eigenstate subspace using quantum computers, and the exact quantum dynamics within the subspace can then be solved classically. We show that the measurement cost of the subspace approach grows polynomially with dimensionality for general potential energy. Our numerical examples demonstrate the capability of our approach, even under intense laser fields. Our work opens the possibility of simulating chemical dynamics with NISQ hardware.
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.
Differential equations underpin a vast range of scientific and engineering domains, including structural mechanics, fluid dynamics, and financial modelling. Their intricate nature often renders classical numerical methods computationally demanding, requiring fine discretization or large basis expansions, which can lead to prohibitive costs. Quantum computing offers a fundamentally different paradigm, leveraging superposition and entanglement with the promise of computational speedup. However, in the current Noisy Intermediate-Scale Quantum (NISQ) era, fully fault-tolerant quantum computers remain decades away, motivating the development of hybrid strategies such as Variational Quantum Algorithms (VQAs), which train parametrised quantum circuits using classical optimizers. This work investigates spectral VQAs for solving differential equations, extending the Chebyshev-encoded VQA of Kyriienko et al. and introducing a novel quantum spectral approach: the Lagrange-encoded VQA. This approach shifts the paradigm by incorporating entanglement directly within the encoding block while maintaining a shallow variational ansatz. Inspired by the Hadamard test structure, the design enables an efficient representation of the complete set of Lagrange polynomial basis functions using a limited number of gates and facilitates the use of the Hadamard test differentiation method, thereby reducing the computational load. VQA performance is analysed through the accuracy of the result in combination with the gate count and circuit count metrics, which serve as proxies for runtime and noise accumulation on NISQ devices. The proposed VQA is benchmarked on multiple problems: (i) a first-order linear ODEs and (ii) the second-order ODEs damped mass-spring system, both comparing Lagrange- and Chebyshev-encoded VQAs, (iii) the Poisson equation, comparing the Lagrange-encoded VQA with the discretized approach of Sato et al., and (iv) a first-order ODE with nonlinear coefficients addressed with the Lagrange-encoded VQA. Across these cases, the Lagrange-encoded VQA consistently achieves comparable or improved accuracy while significantly reducing circuit depth and computational overload, making it a promising candidate for NISQ devices.
Abstract Non-autonomous dynamical systems appear in a very wide range of interesting applications, both in classical and quantum dynamics, where in the latter case, it corresponds to having a time-dependent Hamiltonian. Designing quantum algorithms for time-dependent Hamiltonian is generally much more complicated than time-independent Hamiltonian problems, due to the challenge in handling the time-ordering; moreover, almost all existing such algorithms are developed for qubit-based quantum devices. Here we propose an alternative formalism based on Sambe-Howland’s clock that turns any non-autonomous unitary dynamical system into an autonomous unitary system, i.e. quantum system with a time-independent Hamiltonian, in one higher dimension, while keeping time continuous. This makes the simulation with time-dependent Hamiltonians not much more difficult than that of time-independent Hamiltonians, and can also be framed in terms of an analogue quantum system evolving continuously in time. We show how our new quantum protocol for time-dependent Hamiltonians can be performed in a resource-efficient way and without measurements, and can be made possible on either continuous-variable, qubit or hybrid systems. Combined with a technique called Schrödingerisation, this dilation technique can be applied to the quantum simulation of any linear ODEs and PDEs, and nonlinear ODEs and certain nonlinear PDEs, with time-dependent coefficients.