Curator's Take
AI Commentary
This article introduces an ancilla‑efficient quantum algorithm for simulating forced variable‑coefficient advection‑diffusion equations by exploiting a skew‑symmetric discretization and a Gray‑code Trotter sequence, keeping the logical qubit count to just n_q + 1. By avoiding register‑expanding dilations and heavy postselection, it tackles a key scalability bottleneck that has limited previous quantum PDE solvers, making near‑term hardware more plausible for fluid‑flow and transport applications. The authors demonstrate second‑order temporal convergence (up to fourth order with Richardson extrapolation) on modest lattice sizes, showing that realistic scientific problems can achieve useful accuracy without prohibitive gate depth. Nonetheless, the per‑step postselection failure probability grows with viscous dissipation, so error mitigation will be essential for noisy devices.
— Mark Eatherly
Summary
Quantum time-marching algorithms for transport PDEs often represent variable coefficients and forcing through register-expanding dilations, block-encoding oracles, or repeated postselection. We present an alternative algorithm for a forced variable-coefficient advection-diffusion equation in flow-inspired skew-symmetric form that incorporates spatially varying velocity, viscous dissipation, and persistent source injection with a peak logical requirement of $n_q+1$ qubits. A centered skew-symmetric discretization makes the advection operator strictly skew-Hermitian for arbitrary velocity profiles, enabling an ancilla-free unitary realization using a Gray-code Trotter sequence of controlled-$R_y$ rotations. Diffusion is applied in the Fourier basis through a uniformly controlled rotation on one postselected ancilla, which is measured, reset, and reused between the two diffusion half-steps, while the source is incorporated classically through second-order Strang splitting. Statevector simulations for $N=16$ and $32$ recover second-order temporal convergence against high-accuracy classical solutions, while Richardson extrapolation gives fourth-order accuracy and reduces kernel calls by factors of four to fourteen. Independent tests through $N=256$ confirm second-order spatial consistency. We further show that the per-step ancilla failure probability is proportional to the instantaneous viscous dissipation rate, making postselection cost self-regulating over a fifty-fold viscosity range. Stable evolution is demonstrated for $5\times10^4$ time steps without observable secular error growth, while Gray-code advection accounts for $71$--$95\%$ of transpiled controlled-NOT gates. The fixed-width kernel provides a qubit-efficient building block for near-term hardware studies, although classical readout and state re-preparation remain the main obstacles to coherent multistep evolution.