Optimal Linear Combination of Hamiltonian Simulation for Dissipative Dynamics
Linear combination of Hamiltonian simulation (LCHS) simulates dissipative dynamics on quantum computers. Its query cost depends on the choice of kernel. The open question was whether the Ni–Ying kernel family is optimal for the cost λR.
The minimax error over measures with support in [-R,R] and total variation at most λ is characterized exactly. For the lower bound, weak duality is combined with a realization lemma that approximates extremal dual functions by matrix elements of dissipative evolutions. The realization uses polynomial approximation, time-ordered evolutions, a clock register and clock discretization. The upper bound comes from the Ni–Ying local Wiener norm estimates. Theorem 1.1 was formalized in Lean 4 with Mathlib.
The minimal error is V_R(λ), given in closed form. The optimal cost is (2e/π)(ln(1/ε) − ln ln(1/ε)) + o(1), so the leading constant 2e/π is optimal. The formal Lean proof uses only Lean's standard axioms.
