Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2408.00387 · Europhysics Letters (EPL)

Decomposition of nonlinear collision operator in quantum Lattice Boltzmann algorithm

E. Dinesh Kumar, Steven H. Frankel

Published Oct 28, 2024• 3 citations
arXivHTMLPDF

Abstract

Abstract We propose a quantum algorithm to tackle the quadratic nonlinearity in the Lattice Boltzmann (LB) collision operator. The key idea is to build the quantum gates based on the particle distribution functions (PDF) within the coherence time for qubits. Thus, both the operator and a state vector are linear functions of PDFs, and upon quantum state evolution, the resulting PDFs will have quadraticity. To this end, we decompose the collision operator for a DmQn lattice model into a product of operators, where n is the number of lattice velocity directions. After decomposition, the operators with constant entries remain unchanged throughout the simulation, whereas the remaining will be built based on the statevector of the previous time step. Also, we show that such a decomposition is not unique. Compared to the second-order Carleman-linearized LB, the present approach reduces the circuit width by half and circuit depth by exponential order. The proposed algorithm has been verified through the one-dimensional flow discontinuity and two-dimensional Kolmogrov-like flow test cases.

DOI
10.1209/0295-5075/ad8bef

Related papers

Quantum algorithm for lattice Boltzmann (QALB) simulation of incompressible fluids with a nonlinear collision term

Wael Itani, Katepalli R. Sreenivasan, Sauro Succi A

We present a full quantum algorithm for the lattice Boltzmann method for simulating fluid flows, the only such algorithm to implement both the streaming and collision substeps as unitary operators using an efficient number of qubits. We use Hamiltonian simulation as the main route underlying the algorithm, and show that it uses qubits that scale logarithmically in the flow Reynolds number but gates that scale only polynomially. We motivate the discussion by a brief overview of existing attempts at simulating classical fluids on quantum computers and present a pedagogical discussion on assigning quantum operators to classical variables in both the streaming and the collision substeps, after highlighting the incompatibility of the latter with the implementations of the streaming step that exists in the literature. We use the Bhatnagar–Gross–Krook ansatz for the collision term, representing the relaxation toward an equilibrium distribution. For nonlinear collisions, we use Kowalski's framework that links the nonlinear dynamics of a system to the evolution of bosonic modes, assigning a Carleman linearization order to the truncation in the Fock space of the bosons. We present the qubit and gate complexities, in terms of the chosen accuracy and the Reynolds number. In the Appendix, we work out the details of implementing the operators of the truncated bosonic Fock space in terms of single-qubit gates as well as the error scaling for a general polynomial driving function.

arXiv:2304.05915•Physics of Fluids•Jan 1, 2024•54 citations

Quantum lattice Boltzmann method for several time steps: A local Carleman linearization algorithm

Antonio David Bastida Zamora, Ljubomir Budinski, Valtteri Lahtinen +1 more

This article presents a novel encoding for quantum Lattice Boltzmann method algorithm using Carleman linearization. In contrast to previous articles \cite{Sanavio2024LatticeBC,sanavio2025carleman}, the encoding used allows for local collision rules while keeping a higher probability to obtain the right result, which is of the order of $10^{-2}$. The algorithm scales as $O(log_2^2(N)+Q^3)$ each time step with $N$ the number of lattice sites of the 2D lattice and $Q$ the number of channels with a constant number of qubits when using dynamical circuits.

quant-pharXiv:2511.13072•Physical review. E•Nov 18, 2025•2 citations

Lattice Boltzmann–Carleman quantum algorithm and circuit for fluid flows at moderate Reynolds number

Claudio Sanavio, Sauro Succi A

We present a quantum computing algorithm for fluid flows based on the Carleman-linearization of the Lattice Boltzmann (LB) method. First, we demonstrate the convergence of the classical Carleman procedure at moderate Reynolds numbers, namely, for Kolmogorov-like flows. Then we proceed to formulate the corresponding quantum algorithm, including the quantum circuit layout, and analyze its computational viability. We show that, at least for moderate Reynolds numbers between 10 and 100, the Carleman–LB procedure can be successfully truncated at second order, which is a very encouraging result. We also show that the quantum circuit implementing the single time-step collision operator has a fixed depth, regardless of the number of lattice sites. However, such depth is of the order of ten thousands quantum gates, meaning that quantum advantage over classical computing is not attainable today, but could be achieved in the near or mid-term future. The same goal for the multi-step version remains, however, an open topic for future research.

arXiv:2310.17973•AVS Quantum Science•Apr 22, 2024•40 citations

Quantum Lattice Boltzmann Method based on linear equilibrium distribution functions

Zhengliang Liu, Benzi John, David R. Emerson

In this paper, we propose a complete formulation of the Lattice Boltzmann Method adapted for quantum computing. The classical collision, based on linear equilibrium distribution functions and streaming steps, are reformulated as linear algebraic operations. The inherently non-unitary collision operator is decomposed using Singular Value Decomposition and the Linear Combination of Unitaries technique. Bounce-back boundary conditions are incorporated directly into the collision matrix, while the streaming step is realized through conditional unitary shift operations on spatial registers, controlled by lattice velocity indices encoded in the distribution function register. This formulation ensures that the streaming step remains purely unitary. The resulting quantum circuit is implemented using Qiskit and validated against Couette flow and Poiseuille flow benchmarks. The simulation accurately reproduces the expected velocity profile, with relative errors below 1 0 − 4 . This work establishes a foundational framework for quantum fluid solvers and provides a pathway toward quantum computational fluid dynamics.

Frontiers in Mechanical Engineering•Jan 16, 2026•0 citations

Quantum Lattice Boltzmann with Denoising Collision Operators

Trong Duong, Matthias Möller, Norbert Hosters

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.

quant-pharXiv:2604.09997•Apr 14, 2026•0 citations

Quantum unitary matrix representation of the lattice Boltzmann model for low Reynolds fluid flow simulation

E. Dinesh Kumar, Steven H. Frankel

We propose a quantum algorithm for the lattice Boltzmann (LB) method to simulate fluid flows in the low Reynolds number regime. First, we encode the particle distribution functions (PDFs) as probability amplitudes of the quantum state and demonstrate the need to control the state of the ancilla qubit during the initial state preparation. Second, we express the LB algorithm as a matrix-vector product by neglecting the quadratic non-linearity in the equilibrium distribution function, wherein the vector represents the PDFs, and the matrix represents the collision and streaming operators. Third, we employ classical singular value decomposition to decompose the non-unitary collision and streaming operators into a product of unitary matrices. Finally, we show the importance of having a Hadamard gate between the collision and the streaming operations. Our approach has been tested on linear/linearized flow problems such as the advection-diffusion of a Gaussian hill, Poiseuille flow, Couette flow, and lid-driven cavity problems. We provide counts for two-qubit controlled-NOT and single-qubit U gates for test cases involving 9–12 qubits with grid sizes ranging from 24 to 216 points. While the gate count aligns closely with theoretical limits, the high number of two-qubit gates on the order of 107 necessitates careful attention to circuit synthesis.

AVS Quantum Science•Mar 1, 2025•15 citations
Abstract We propose a quantum algorithm to tackle the quadratic nonlinearity in the Lattice Boltzmann (LB) collision operator. The key idea is to build the quantum gates based on the particle distribution functions (PDF) within the coherence time for qubits. Thus, both the operator and a state vector are linear functions of PDFs, and upon quantum state evolution, the resulting PDFs will have quadraticity. To this end, we decompose the collision operator for a DmQn lattice model into a product of operators, where n is the number of lattice velocity directions. After decomposition, the operators with constant entries remain unchanged throughout the simulation, whereas the remaining will be built based on the statevector of the previous time step. Also, we show that such a decomposition is not unique. Compared to the second-order Carleman-linearized LB, the present approach reduces the circuit width by half and circuit depth by exponential order. The proposed algorithm has been verified through the one-dimensional flow discontinuity and two-dimensional Kolmogrov-like flow test cases.
We present a full quantum algorithm for the lattice Boltzmann method for simulating fluid flows, the only such algorithm to implement both the streaming and collision substeps as unitary operators using an efficient number of qubits. We use Hamiltonian simulation as the main route underlying the algorithm, and show that it uses qubits that scale logarithmically in the flow Reynolds number but gates that scale only polynomially. We motivate the discussion by a brief overview of existing attempts at simulating classical fluids on quantum computers and present a pedagogical discussion on assigning quantum operators to classical variables in both the streaming and the collision substeps, after highlighting the incompatibility of the latter with the implementations of the streaming step that exists in the literature. We use the Bhatnagar–Gross–Krook ansatz for the collision term, representing the relaxation toward an equilibrium distribution. For nonlinear collisions, we use Kowalski's framework that links the nonlinear dynamics of a system to the evolution of bosonic modes, assigning a Carleman linearization order to the truncation in the Fock space of the bosons. We present the qubit and gate complexities, in terms of the chosen accuracy and the Reynolds number. In the Appendix, we work out the details of implementing the operators of the truncated bosonic Fock space in terms of single-qubit gates as well as the error scaling for a general polynomial driving function.
This article presents a novel encoding for quantum Lattice Boltzmann method algorithm using Carleman linearization. In contrast to previous articles \cite{Sanavio2024LatticeBC,sanavio2025carleman}, the encoding used allows for local collision rules while keeping a higher probability to obtain the right result, which is of the order of 10−210^{-2}10−2. The algorithm scales as O(log22(N)+Q3)O(log_2^2(N)+Q^3)O(log22​(N)+Q3) each time step with NNN the number of lattice sites of the 2D lattice and QQQ the number of channels with a constant number of qubits when using dynamical circuits.
We present a quantum computing algorithm for fluid flows based on the Carleman-linearization of the Lattice Boltzmann (LB) method. First, we demonstrate the convergence of the classical Carleman procedure at moderate Reynolds numbers, namely, for Kolmogorov-like flows. Then we proceed to formulate the corresponding quantum algorithm, including the quantum circuit layout, and analyze its computational viability. We show that, at least for moderate Reynolds numbers between 10 and 100, the Carleman–LB procedure can be successfully truncated at second order, which is a very encouraging result. We also show that the quantum circuit implementing the single time-step collision operator has a fixed depth, regardless of the number of lattice sites. However, such depth is of the order of ten thousands quantum gates, meaning that quantum advantage over classical computing is not attainable today, but could be achieved in the near or mid-term future. The same goal for the multi-step version remains, however, an open topic for future research.
In this paper, we propose a complete formulation of the Lattice Boltzmann Method adapted for quantum computing. The classical collision, based on linear equilibrium distribution functions and streaming steps, are reformulated as linear algebraic operations. The inherently non-unitary collision operator is decomposed using Singular Value Decomposition and the Linear Combination of Unitaries technique. Bounce-back boundary conditions are incorporated directly into the collision matrix, while the streaming step is realized through conditional unitary shift operations on spatial registers, controlled by lattice velocity indices encoded in the distribution function register. This formulation ensures that the streaming step remains purely unitary. The resulting quantum circuit is implemented using Qiskit and validated against Couette flow and Poiseuille flow benchmarks. The simulation accurately reproduces the expected velocity profile, with relative errors below 1 0 − 4 . This work establishes a foundational framework for quantum fluid solvers and provides a pathway toward quantum computational fluid dynamics.
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.
We propose a quantum algorithm for the lattice Boltzmann (LB) method to simulate fluid flows in the low Reynolds number regime. First, we encode the particle distribution functions (PDFs) as probability amplitudes of the quantum state and demonstrate the need to control the state of the ancilla qubit during the initial state preparation. Second, we express the LB algorithm as a matrix-vector product by neglecting the quadratic non-linearity in the equilibrium distribution function, wherein the vector represents the PDFs, and the matrix represents the collision and streaming operators. Third, we employ classical singular value decomposition to decompose the non-unitary collision and streaming operators into a product of unitary matrices. Finally, we show the importance of having a Hadamard gate between the collision and the streaming operations. Our approach has been tested on linear/linearized flow problems such as the advection-diffusion of a Gaussian hill, Poiseuille flow, Couette flow, and lid-driven cavity problems. We provide counts for two-qubit controlled-NOT and single-qubit U gates for test cases involving 9–12 qubits with grid sizes ranging from 24 to 216 points. While the gate count aligns closely with theoretical limits, the high number of two-qubit gates on the order of 107 necessitates careful attention to circuit synthesis.