Saved in:
Bibliographic Details
Main Authors: Das, Elin Ranjan, Zheng, Muqing, Dutta, Rishab, Li, Ang, Stavenger, Timothy, Liu, Yuan
Format: Preprint
Published: 2026
Subjects:
Online Access:https://arxiv.org/abs/2605.10708
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866918494658887680
author Das, Elin Ranjan
Zheng, Muqing
Dutta, Rishab
Li, Ang
Stavenger, Timothy
Liu, Yuan
author_facet Das, Elin Ranjan
Zheng, Muqing
Dutta, Rishab
Li, Ang
Stavenger, Timothy
Liu, Yuan
contents We introduce a hybrid oscillator-qubit formulation of linear combination of Hamiltonian simulation (LCHS) for solving linear ordinary differential equations. Instead of representing the quadrature rule with a discrete-variable (DV) ancilla register in qubit-only LCHS, the method encodes the LCHS kernel in a continuous-variable (CV) ancillary mode, thereby eliminating the explicit $O(\log M_a)$ ancilla-qubit overhead, where $M_a$ is the number of discretized integral terms in the DV quadrature rule. We derive analytical error bounds for two main approximation mechanisms for the ideal kernel state preparation, showing superalgebraic convergence for Schwartz-class kernels in the truncation cutoff $N$. The required CV non-Gaussianity is captured by the finite squeezed-Fock kernel state, which generically has stellar rank $N-1$, identifying the truncation cutoff as a discrete measure of the oracle's non-Gaussian resource. For the hybrid oscillator-qubit evolution, we also obtain a product-formula bound showing that a $p$th-order formula requires $O(t^{1+1/p}(Γ_{p,N}/ε_t)^{1/p})$ Trotter steps to reach error $ε_t$, where $Γ_{p,N}$ collects Pauli commutator terms weighted by powers of the truncated position-operator norm $\|\hat{x}\|_N$. We further derive a perturbation bound for the probability of obtaining the required oscillator measurement outcome, showing that an $ε$-close implementation of the ideal LCHS oracle in operator norm induces only an $O(ε)$ perturbation in the postselection probability. In the heat-equation benchmarks, the Law--Eberly protocol achieves end-to-end solution fidelity at least 99.90%. A comparison with a matrix-product-state-based DV LCHS implementation further shows that, the hybrid construction uses a substantially more compact oracle description with reduce circuit cost.
format Preprint
id arxiv_https___arxiv_org_abs_2605_10708
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle Quantum Differential Equation Solver via Hybrid Oscillator-Qubit Linear Combination of Hamiltonian Simulations
Das, Elin Ranjan
Zheng, Muqing
Dutta, Rishab
Li, Ang
Stavenger, Timothy
Liu, Yuan
Quantum Physics
We introduce a hybrid oscillator-qubit formulation of linear combination of Hamiltonian simulation (LCHS) for solving linear ordinary differential equations. Instead of representing the quadrature rule with a discrete-variable (DV) ancilla register in qubit-only LCHS, the method encodes the LCHS kernel in a continuous-variable (CV) ancillary mode, thereby eliminating the explicit $O(\log M_a)$ ancilla-qubit overhead, where $M_a$ is the number of discretized integral terms in the DV quadrature rule. We derive analytical error bounds for two main approximation mechanisms for the ideal kernel state preparation, showing superalgebraic convergence for Schwartz-class kernels in the truncation cutoff $N$. The required CV non-Gaussianity is captured by the finite squeezed-Fock kernel state, which generically has stellar rank $N-1$, identifying the truncation cutoff as a discrete measure of the oracle's non-Gaussian resource. For the hybrid oscillator-qubit evolution, we also obtain a product-formula bound showing that a $p$th-order formula requires $O(t^{1+1/p}(Γ_{p,N}/ε_t)^{1/p})$ Trotter steps to reach error $ε_t$, where $Γ_{p,N}$ collects Pauli commutator terms weighted by powers of the truncated position-operator norm $\|\hat{x}\|_N$. We further derive a perturbation bound for the probability of obtaining the required oscillator measurement outcome, showing that an $ε$-close implementation of the ideal LCHS oracle in operator norm induces only an $O(ε)$ perturbation in the postselection probability. In the heat-equation benchmarks, the Law--Eberly protocol achieves end-to-end solution fidelity at least 99.90%. A comparison with a matrix-product-state-based DV LCHS implementation further shows that, the hybrid construction uses a substantially more compact oracle description with reduce circuit cost.
title Quantum Differential Equation Solver via Hybrid Oscillator-Qubit Linear Combination of Hamiltonian Simulations
topic Quantum Physics
url https://arxiv.org/abs/2605.10708