⚛️ quantum physics

Two-dimensional quantum lattice gas algorithm for anisotropic Burger-like equations

This paper presents a refined two-dimensional quantum lattice gas algorithm that corrects predicted viscosity and utilizes a minimal two-velocity model to simulate anisotropic Burger-like equations, offering a promising quantum-native pathway toward momentum-conserving Navier-Stokes dynamics.

Original authors: Niccoló Fonio, Pierre Sagaut, Giuseppe Di Molfetta

Published 2026-02-20

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Niccoló Fonio, Pierre Sagaut, Giuseppe Di Molfetta

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to simulate how water flows, smoke drifts, or traffic jams form. In the world of classical computers, we do this using a method called "Lattice Gas," which is like a giant digital chessboard. On this board, little particles (or "chess pieces") hop from square to square, bumping into each other, and following rules that mimic real-world physics like viscosity (thickness) and momentum.

Now, imagine trying to run this simulation on a quantum computer. Quantum computers are incredibly powerful but also very fragile; they operate on "qubits" (quantum bits) that can be in two states at once. The challenge is: How do we make these delicate quantum particles behave like a fluid without breaking the rules of quantum mechanics?

This paper, by Niccolò Fonio and colleagues, is like a blueprint for a new, more efficient way to do exactly that. Here is the breakdown in simple terms:

1. The Problem with the Old Blueprint

Previously, scientists (like Jeffrey Yepez) had a quantum model that worked well in one dimension (like a single line of traffic). It was "unitary," meaning it followed the strict, reversible laws of quantum physics perfectly. However, when they tried to expand this to two dimensions (a real 2D grid), things got messy.

The Viscosity Glitch: The old math predicted the fluid would be "thicker" (more viscous) than it actually was in the simulation. It was like a recipe that said "add 2 cups of flour," but the cake turned out dense because the math was slightly off.
The Complexity Trap: To simulate 2D fluids on a quantum computer, most methods require a huge number of qubits (quantum bits), which is like trying to build a skyscraper with a pile of bricks that is too heavy to lift.

2. The New Solution: A "Minimalist" Quantum Kitchen

The authors decided to revisit the old recipe and fix the math. They introduced a minimalist approach:

The "Two-Speed" Trick: Instead of using a complex grid with many different speeds for particles, they found a way to simulate 2D fluid dynamics using just two qubits per grid cell and only two possible directions for movement.
The Analogy: Imagine a busy kitchen. The old way required a chef to juggle 10 different pots, pans, and ingredients to make a soup. The new way shows you can make the exact same soup using just a pot and a spoon, provided you know the exact right timing and temperature.

3. Fixing the "Thickness" (Viscosity)

One of the paper's biggest achievements is fixing the viscosity calculation.

The Discovery: They derived a new mathematical correction. Think of it as realizing that the "friction" in their quantum fluid wasn't just a simple number; it depended on the angle of the "collision" between particles in a specific way.
The Result: With this new formula, the quantum simulation now matches the real physics perfectly. They proved that this model can simulate fluids with arbitrarily low viscosity (making the fluid almost frictionless) without the simulation crashing. This is a huge deal because other methods usually become unstable when you try to make the fluid too "slippery."

4. The 2D Breakthrough: Anisotropic Burgers

The authors took their 1D model and stretched it into 2D.

The Result: They created a model that simulates "Burger-like" equations. In plain English, this describes how waves (like shockwaves in a sonic boom or a traffic jam) form and move.
The Catch: Because they only used two directions (like a car that can only go North or East, but not diagonally), the fluid behaves slightly differently depending on the direction. The authors call this "anisotropic."
The Metaphor: Imagine walking through a forest. If you can only walk North or East, your path looks different than if you could walk in any direction. The fluid "flows" differently along the grid lines. However, the authors showed that by choosing the right grid layout (like a triangular lattice vs. a square one), you can control these quirks.

5. Why Does This Matter?

This paper is a bridge between two worlds:

Quantum Computing: It shows a way to simulate complex physics (fluids) using very few quantum resources (qubits). This is crucial because current quantum computers are small and noisy; we need algorithms that don't need thousands of qubits to work.
Fluid Dynamics: It offers a new way to simulate turbulence and shockwaves that is more stable than current classical methods when dealing with very low viscosity.

The Bottom Line

The authors have built a lean, mean, quantum fluid simulator. They fixed the math to get the "thickness" of the fluid right, proved it works in 2D, and showed that you can do it with very few quantum bits.

It's like taking a heavy, clunky steam engine (old methods) and replacing it with a sleek, high-efficiency electric motor (this new model). It runs smoother, uses less energy (fewer qubits), and can handle the tricky parts of the road (low viscosity) without stalling. This paves the way for future quantum computers to simulate weather patterns, aerodynamics, and ocean currents in ways we couldn't do before.

1. Problem Statement

The application of Quantum Computing (QC) to Computational Fluid Dynamics (CFD) faces challenges in balancing quantum advantage with physical accuracy and numerical stability. Existing approaches often rely on encoding classical lattice gas models (like Lattice Boltzmann Method - LBM) into quantum states, which can suffer from non-unitary collision operators, requiring complex multi-time-step implementations or linearization techniques that struggle with turbulence.

Specifically, the authors revisit the Yepez Quantum Lattice Gas (QLG) model (Phys. Rev. A, 2006). While this model is unitary and conceptually complete, previous formulations contained an approximation in the derivation of the macroscopic viscosity. Furthermore, extending this model to two dimensions (2D) while maintaining analytical tractability and minimizing resource usage (qubits) remains an open challenge. The goal is to derive a corrected viscosity, generalize the model to 2D with minimal resources (2 qubits per site), and determine the resulting macroscopic partial differential equations (PDEs).

2. Methodology

A. Theoretical Framework: Q-D1Q2 Model

The authors analyze a 1D Quantum Lattice Gas with 2 qubits per site (Q-D1Q2).

State Representation: Each lattice site $x$ holds a state $|\psi(x,t)\rangle = |q_0\rangle \otimes |q_1\rangle$ , where qubits represent particle populations moving with discrete velocities.
Algorithm: The evolution consists of three steps:
1. Initialization: Qubits are initialized based on classical distribution functions $f_i$ .
2. Collision: A unitary operator $\hat{C}$ (parametrized by Euler angles $\theta, \zeta, \xi$ ) is applied to conserve mass.
3. Streaming & Measurement: Populations are measured (collapsing the state) to update distributions, followed by classical streaming.
Derivation: The authors perform a Chapman-Enskog expansion on the lattice Boltzmann equation derived from the quantum evolution. They expand the finite difference equation up to second order in space ( $\delta x$ ) and first order in time ( $\delta t$ ), assuming the diffusive scaling regime ( $\delta x \sim \epsilon, \delta t \sim \epsilon^2$ ).

B. Viscosity Correction

A key methodological step is the rigorous calculation of the term $J_1 - J_0$ (the derivative of the collision term with respect to populations at equilibrium).

Previous work (Yepez) approximated this term, leading to a viscosity expression $\nu \propto \cot^2(\theta)$ .
The authors derive an exact closed-form solution for the equilibrium populations and the collision term, resulting in a corrected viscosity formula that includes higher-order terms in the low-Mach limit.

C. 2D Generalization (Q-D2Q2)

The model is extended to 2D using a square lattice with only 2 qubits per site (minimal resource usage).

Streaming: Populations stream in both $x$ and $y$ directions.
Collision: The same 1D collision operator is applied locally.
Analysis: The authors derive the macroscopic PDE for various discrete velocity sets (e.g., square lattice, triangular lattice) to observe how anisotropy and advection terms emerge.

3. Key Contributions

Analytical Correction to Viscosity:
The authors derive a corrected expression for the kinematic viscosity $\nu$ in the 1D model. The new formula is:
$\nu = -\frac{\delta x^2}{\delta t} \frac{1}{2} \left( 1 - \frac{\sqrt{\alpha^2 + 1}}{\sin^2(\theta)} \right)$
where $\alpha = \cot(\theta)\cos(\zeta - \xi)$ . This corrects the previous $\cot^2(\theta)$ approximation, extending the validity of the Burgers-like equation to order $O(\epsilon^4)$ .
Minimal 2D Generalization:
They propose a Q-D2Q2 model (2 dimensions, 2 qubits per site). Unlike traditional 2D lattice gases (like FHP) which require 6 or 7 velocities, this model uses only 2.
- Result: The model successfully reproduces anisotropic Burgers-like dynamics.
- Equation: The derived macroscopic equation is:
  $\partial_t \rho + \mathbf{d} \cdot \nabla \rho - \nabla \cdot (\mathbf{D} \nabla \rho) = 0$
  where $\mathbf{d}$ is an advective vector and $\mathbf{D}$ is an anisotropic diffusion tensor dependent on the chosen velocity set.
Arbitrarily Small Viscosity with Stability:
The paper demonstrates that the QLG model can achieve arbitrarily small viscosities without the numerical instability typically found in single-relaxation-time LBM models (which become unstable as the relaxation time $\tau \to 0.5$ ). This is due to the model's inherent unitary structure and satisfaction of an H-theorem.

4. Results

Numerical Validation of Viscosity:
Simulations in 1D confirmed the corrected viscosity formula. The experimental viscosity matched the theoretical prediction (Eq. 22) significantly better than the previous approximation (Eq. 19), especially in the low-viscosity regime.
Shock Formation and Stability:
- The model successfully simulated shock formation in Burgers' equation.
- Trade-off: Achieving very low viscosity (sharp shocks) requires longer simulation times for thermalization (relaxation) and/or higher grid resolution, but the scheme remains unconditionally stable.
2D Anisotropic Behavior:
- Different velocity sets produced different macroscopic behaviors.
- Symmetric sets (e.g., $c_0 = -c_1$ ) resulted in 1D-like dynamics projected onto 2D.
- Asymmetric sets (e.g., $c_0$ horizontal, $c_1$ vertical) introduced cross-derivative terms ( $\partial_{xy}\rho$ ) and additional advective terms, confirming the anisotropic nature of the minimal 2D model.
Comparison with Classical Methods:
Numerical comparisons between the QLG simulation and a Finite Difference Method (FDM) solver for the derived 2D PDE showed excellent qualitative and quantitative agreement ( $L_2$ norm error analysis), validating the derived PDE.

5. Significance and Future Perspectives

Quantum-Native Simulation: The work establishes a path for "quantum-native" fluid simulations that do not merely encode classical algorithms but utilize unitary quantum dynamics to naturally embed dissipative phenomena (viscosity) within reversible microscopic rules.
Resource Efficiency: By achieving 2D anisotropic dynamics with only 2 qubits per site, the model offers a resource-efficient alternative to classical lattice gas automata (like FHP) and standard LBM, which require more states to capture similar physics.
Bridge to Navier-Stokes: While the current 2D model captures Burgers-like dynamics (scalar conservation), the authors identify that adding a third velocity (3 qubits) could enable momentum conservation, potentially leading to a full Navier-Stokes solver on a quantum lattice.
Hardware Relevance: The minimal qubit requirement and unconditional stability make this approach particularly promising for implementation on near-term quantum hardware (NISQ devices).

In conclusion, this paper refines the theoretical foundation of quantum lattice gas algorithms, corrects a long-standing approximation in viscosity derivation, and successfully demonstrates a minimal, resource-efficient 2D extension capable of simulating complex nonlinear fluid dynamics.