Generalize LCHS analysis to time-dependent or inhomogeneous differential equations

Extend the analysis of the linear combination of Hamiltonian simulation (LCHS) method from the homogeneous, time-independent setting du(t)/dt = −A u(t) with constant matrix A to time-dependent matrices A(t) and to inhomogeneous linear differential equations. Determine the appropriate LCU-based propagator representations, specify the truncation and discretization schemes that yield rigorous error bounds, and characterize the resulting sampling overhead via the reduction factor R within the proposed hybrid LCU framework.

Background

The paper develops a hybrid linear combination of unitaries (LCU) strategy that interpolates between coherent and randomized implementations, introducing a reduction factor R to quantify sampling overhead. In the application to differential equations, the authors analyze a specific LCHS construction for the homogeneous, time-independent case du/dt = −A u(t), deriving an integral representation of e{−AT} and practical truncation and discretization procedures, together with bounds on R.

While the presented analysis is for constant A and homogeneous dynamics, many practical problems involve time-dependent generators A(t) or inhomogeneous terms. Extending the LCHS framework to these settings would require deriving suitable LCU decompositions and adapting the error and resource analysis, thereby enabling the same tradeoff between circuit size and sampling overhead for a broader class of dynamical systems.

References

We here consider a specific example for solving differential equations ${\rm d}u(t)/{\rm d}t = -Au(t)$. Here, $A$ is a general time-independent matrix. A generalization to time-dependent cases or inhomogeneous cases of the following analysis would be possible, but we leave it for future work.

Tradeoffs between quantum and classical resources in linear combination of unitaries  (2512.06260 - Wada et al., 6 Dec 2025) in Section 4.1 (Application: Linear combination of Hamiltonian simulation)