Surrogate Quantum Circuit Design for the Lattice Boltzmann Collision Operator

TL;DR

This paper develops a low-depth surrogate quantum circuit to efficiently approximate the non-unitary collision operator in the lattice Boltzmann method, enabling quantum simulations of fluid dynamics with preserved physical properties.

Contribution

It introduces a novel quantum circuit design that models dissipative nonlinear dynamics without ancilla qubits or post-selection, tailored for quantum fluid simulations.

Findings

Accurately reproduces vortex decay and flow recirculation in benchmark tests.

Requires only 724 native gates on IBM quantum hardware, independent of lattice size.

Preserves key physical properties like mass and scale conservation in the quantum circuit.

Abstract

This study introduces a framework for learning a low-depth surrogate quantum circuit (SQC) that approximates the nonlinear, dissipative, and hence non-unitary Bhatnagar-Gross-Krook (BGK) collision operator in the lattice Boltzmann method (LBM) for the D2Q9 lattice. By appropriately selecting the quantum state encoding, circuit architecture, and measurement protocol, non-unitary dynamics emerge naturally within the physical population space. This approach removes the need for probabilistic algorithms relying on any ancilla qubits and post-selection to reproduce dissipation, or for multiple state copies to capture nonlinearity. The SQC is designed to preserve key physical properties of the BGK operator, including mass conservation, scale equivariance, and D8 equivariance, while momentum conservation is encouraged through penalization in the training loss. When compiled to the IBM Heron quantum processor's native gate set, assuming all-to-all qubit connectivity, the circuit requires only 724 native gates and operates locally on the velocity register, making it independent of the lattice size. The learned SQC is validated on two benchmark cases, the Taylor-Green vortex decay and the lid-driven cavity, showing accurate reproduction of vortex decay and flow recirculation. While integration of the SQC into a quantum LBM framework presently requires measurement and re-initialization at each timestep, the necessary steps towards a measurement-free formulation are outlined.