Reliable multistep quantum lattice-Boltzmann evolution requires controlling computational faults. When ideal collision preserves encoded mass and momentum exactly, changes in those charges can provide syndromes for selected faults at the post-collision checkpoint. A fault that shifts one population then leaves a unique mass-momentum residual labeled by the lattice velocity. A coherent reference register stores the expected charges, so the same test remains valid while the physical fluid charges vary across the lattice. Single bit flips on population number registers occupy disjoint syndrome sectors and satisfy the Knill-Laflamme condition. A state-vector demonstration on eighteen data qubits of the two-dimensional nine-velocity lattice (D2Q9) recovers that single-bit family to double-precision roundoff. Conservation still leaves a fifteen-dimensional kinetic nullspace on the three-dimensional nineteen-velocity lattice (D3Q19). Exact multiple-relaxation-time (MRT) streaming analysis ranks those charge-preserving modes by the order at which they return to density or momentum. A classical D3Q19 decaying-flow simulation at moderate Reynolds number then injects about $1.75\times10^4$ selected single-population shifts per realization. Recovery of that alphabet returns the trajectory to floating-point roundoff. The quantum statements assume an ideal charge-preserving collision and a reliable reference. The result is an inner recovery map for a stated charge-changing alphabet together with a kinetic classification of the unresolved sector.
Reliable multistep quantum lattice-Boltzmann evolution requires controlling computational faults. When ideal collision preserves encoded mass and momentum exactly, changes in those charges can provide syndromes for selected faults at the post-collision checkpoint. A fault that shifts one population then leaves a unique mass-momentum residual labeled by the lattice velocity. A coherent reference register stores the expected charges, so the same test remains valid while the physical fluid charges vary across the lattice. Single bit flips on population number registers occupy disjoint syndrome sectors and satisfy the Knill-Laflamme condition. A state-vector demonstration on eighteen data qubits of the two-dimensional nine-velocity lattice (D2Q9) recovers that single-bit family to double-precision roundoff. Conservation still leaves a fifteen-dimensional kinetic nullspace on the three-dimensional nineteen-velocity lattice (D3Q19). Exact multiple-relaxation-time (MRT) streaming analysis ranks those charge-preserving modes by the order at which they return to density or momentum. A classical D3Q19 decaying-flow simulation at moderate Reynolds number then injects about 1.75×104 selected single-population shifts per realization. Recovery of that alphabet returns the trajectory to floating-point roundoff. The quantum statements assume an ideal charge-preserving collision and a reliable reference. The result is an inner recovery map for a stated charge-changing alphabet together with a kinetic classification of the unresolved sector.
Turbulent fluid simulation is computationally demanding because nonlinear interactions couple a wide range of scales. Quantum computing offers another computational model, but fluid evolution is nonlinear and dissipative whereas closed gate-based dynamics is linear and reversible. The lattice Boltzmann method (LBM) is an attractive discrete target because streaming is a local permutation, yet conventional collision reconstructs a nonlinear equilibrium and applies dissipative relaxation. We present a digital quantum lattice Boltzmann (QLBM) algorithm that keeps that nonlinear timestep in computational-basis registers. Each step writes the updated populations into a clean destination bank, uncomputes the workspace, and resets the obsolete source. The resulting channel is completely positive and trace preserving (CPTP), composes without intermediate measurement, and applies to any lattice Boltzmann stencil that can be evaluated reversibly into a clean bank. Forced D3Q19 homogeneous isotropic turbulence, in direct numerical simulation (DNS) and Smagorinsky large-eddy simulation (LES) at Reynolds number 15000 over one million lattice steps, shows statistical agreement in mass, energy budgets, spectra, and intermittency. The three-dimensional calculations use a digital-register emulator of the circuit. A separate D2Q9 implementation compiles the same contract to reversible gates and agrees bit-exactly with its integer twin. Two-bank storage and serial lattice depth limit near-term feasibility.
Lattice Boltzmann (LB) on quantum devices must reconcile unitary gate evolution with the dissipative \emph{collision} step. In the multiple-relaxation-time (MRT) class, we work in the common setting of \emph{modewise diagonal} moment relaxation, δmr′=λrδmr with λr∈[−1,1] (overrelaxation if λr<0). Embedding that contraction in a unitary by block encoding or a linear combination of unitaries (LCU) typically yields subunitary success probability that decays multiplicatively across modes, sites, and time, a key bottleneck for quantum LB. \emph{For the dissipative MRT block alone} we give a \emph{block-encoding-free} construction: a signed \emph{two-rail} population encoding, then a completely positive trace-preserving (CPTP) map (per-rail amplitude damping with survival ∣λr∣ and, if λr<0, a rail SWAP) so that, after the decode, the map agrees with classical MRT relaxation exactly (expectations of the rail number operators, common encoding--decode scale). Trace preservation gives success probability 1 for that substage. The main result is the dissipative MRT block; construction of the equilibrium moment vector~meq=Mfeq (prescribed~feq, host moment matrix~M; notation as in Section~\ref{subsec:generic-mrt}), moment transforms, streaming, and boundaries are composed with it as in a standard host pipeline and lie outside the scope of the formal theorem. Hybrid and fully coherent encodings, adaptive scales, Carleman-based context, and a one-rail no-go in the same nonnegative population framework are in the main text. Audits of the open-channel map on a long LBM collide-stream simulation and on stencil-free inputs both match the target to machine precision.
The quantum lattice Boltzmann scheme of Succi and Dellar advances a four-component Dirac spinor on a lattice by a fixed sequence of local, exactly norm-preserving operations: a basis rotation, a collision, a streaming shift, and the inverse rotation. While this correspondence has been noted qualitatively, it has remained unclear whether it can be realized exactly and unitarily as a gate-model quantum circuit. This unitarity is a structural property of the scheme, which should, in principle, enable an operator-by-operator port onto quantum hardware. We give a gate-level construction of every operation of the three-dimensional Dirac QLB scheme: the fixed rotation gates, the collision gate, the streaming shift as a controlled increment on a position register, the position-dependent potential as a phase oracle, and periodic and reflecting boundary conditions as unitary circuits. Together, these constructions provide the first complete demonstration that the scheme is exactly representable in the gate model without any approximation beyond the classical method itself. We then compose them into single-axis, two- and three-dimensional time steps. On a state-vector emulator the resulting circuits reproduce the classical QLB solver to machine precision, thus verifying our quantum implementation. Our aim is representability and correctness rather than computational speedup; questions of state preparation, asymptotic cost, and hardware noise are identified as open avenues for future work. Appendices restore physical units, derive the scheme's Zitterbewegung in closed form and compare it with a trapped-ion simulation, and execute a reduced two-site instance on a Rigetti superconducting processor, where choosing the qubits from the live calibration lowers the state-preparation error. All operators, circuits, and figures are reproducible from the open-source quantumKineticMethods library.
The Carleman lattice Boltzmann method replaces the quadratic collision of the lattice Boltzmann equation by a linear map on a lifted state. Existing block encodings of this map succeed with probability 10−5 to 10−2 per step, so that the probability of a multi-step evolution decays exponentially with the number of steps. This study develops multi-relaxation-time Carleman lattice Boltzmann, a representation split in which the Carleman-lifted collision is performed in a moment basis on the amplitude-encoded state, while streaming remains a permutation in population space. The level-2 hierarchy closes under collision but not under streaming, which generates the work of pressure, advection and stress on the momentum flux. The state therefore remains the full pair array, and the moment representation is an orthogonal change of basis of the population one, with identical dynamics and truncation error. The operator to be block encoded changes with the basis, and with it the success probability of a quantum step. In the weighted Hermite basis with the amplitude encoding fi/wi, the relaxation is a diagonal operator of norm one, and the level-2 collision stage becomes this diagonal plus a sparse coupling from products of conserved moments. A structured block encoding of this form raises the single-step success probability on the D2Q9 lattice from 10−2 for the best population encoding to 0.98. With the pair sector in relative coordinates and scaled, one step requires O(log2N) Toffoli gates and is verified by statevector simulation. The encoding and readout overhead of a T-step run is then polynomial in T and N rather than accumulating a factor pT. For resolved flows, the remaining exponential loss in coherent multi-step Carleman lattice Boltzmann evolution is set by the non-equilibrium norm relaxed by dissipation rather than by the block encoding.
Quantum computing shows substantial potential in accelerating simulations and alleviating memory bottlenecks in computational fluid dynamics (CFD), owing to its inherent properties of superposition and entanglement. The lattice Boltzmann method (LBM), being largely algebraic in nature, has inspired the development of various quantum LBMs. However, most existing approaches fix the relaxation time at τ = 1, thereby confining a given mesh resolution to simulations at a single Reynolds number. Although our earlier quantum lattice kinetic scheme (LKS) lifted this restriction, it suffers from instability at high Reynolds numbers. To address this challenge, we propose a quantum fractional-step LBM (FS-LBM). In this framework, the predictor step is implemented on a quantum circuit using the standard LBM formulation, while the corrector step is performed classically. The relaxation time is retained at τ = 1 to ensure seamless compatibility with existing quantum LBMs. Benchmark simulations of representative two- and three-dimensional incompressible isothermal and thermal flows demonstrate that the quantum FS-LBM achieves accuracy and convergence orders consistent with its classical counterpart, while significantly outperforming the quantum LKS in both precision and stability. Notably, this work presents the first quantum LBM simulation of three-dimensional incompressible thermal flows.
The Lattice Boltzmann method (LBM) is a well-established mesoscopic approach for simulating fluid dynamics by evolving particle distribution functions on discrete lattices. While the LBM is highly parallelizable on classical hardware, its translation to quantum algorithms is impeded by the collision process, which is intrinsically nonlinear and irreversible. Several existing quantum formulations implement this process through repeated quantum tomography and state preparation at every timestep, leading to significant overheads. We introduce a quantum LBM based on a denoising-type collision operator that avoids tomography-based updates. The collision dynamics are reformulated as an orthogonal projection onto the linearized manifold of equilibrium distributions around a reference state. This geometric approach filters non-equilibrium components while preserving lattice symmetries and approximating nonlinear terms needed to recover hydrodynamic behavior. A complete pipeline is presented with efficient gate-level realizations, incorporating encoding of distributions, collision, streaming, boundary conditions, and measurement of physical quantities such as hydrodynamic forces. In addition, we outline an approach for implementing projector-based operators deterministically without postselection, paving the way to fully coherent multi-timestep LBM simulations. Numerical experiments for advection-diffusion and flow problems demonstrate that the method reproduces macroscopic behaviors with high accuracy, with performance depending on the choice of reference state.