Quantum algorithms

Distributed Variational Quantum Linear Solver

A hybrid quantum-classical method that distributes a block-partitioned linear system across networked near-term quantum processors.

Year
2026
Status
Submitted
Topics
Quantum computing · Distributed optimization · Linear systems

The project in brief

Can several small near-term quantum computers cooperate to solve a linear system too large for any one of them?

Can several small quantum processors represent a solution to a linear system too large for one of them? Here the matrix is partitioned across agents, and only classical messages pass between processors.

The central connection is a least-squares identity: a graph Laplacian removes disagreement between local residuals while preserving their average. This gives a distributed objective; quantum circuits then parameterize its unknown vectors. The identity is exact, whereas training those circuits remains a nonconvex optimization problem.

A block matrix distributed across a grid of quantum computers

A conceptual illustration, not a plot of experimental results.

Objective

One global problem, smaller local state spaces

AA is a real 2n×2n2^n\times2^n matrix and b∈R2nb\in\mathbb R^{2^n}. Each agent holds one matrix block and uses qq system qubits to encode a local vector. Here nn counts global system qubits; ancillary qubits for measurement circuits are additional.

x⋆∈arg min⁡x∈R2n∥Ax−b∥22,A=[Aij]i,j=1m,Aij∈R2q×2q,m=2n−q,number of agents=m2.\begin{aligned} x^\star &\in \operatorname*{arg\,min}_{x\in\mathbb R^{2^n}}\|Ax-b\|_2^2,\\ A&=[A_{ij}]_{i,j=1}^{m},\qquad A_{ij}\in\mathbb R^{2^q\times2^q},\\ m&=2^{n-q},\qquad \text{number of agents}=m^2. \end{aligned}

Manuscript · Problem Formulation

Question
Recover a global least-squares solution from block-local data, normalized quantum states, and classical communication.
Method
Laplacian residual splitting, quantum direction-plus-scale representations, and column-wise gradient tracking with Adam.
Progress
An exact reformulation lemma and three numerical studies: representation size, partitioning, and optimizer ablation. No hardware or quantum-speedup claim.

Explore the project

01 / Problem

One linear system, many local views

The mathematical target is ordinary least squares. The constraints concern where the matrix lives, how vectors are represented, and which agents can communicate.

Read chapter
02 / Method

A graph Laplacian makes the objective local

The Laplacian of a connected undirected graph has every zero-sum vector in its image. Local auxiliary variables can remove residual disagreement, but cannot change the component representing the original equation.

Read chapter
03 / Progress & next steps

What the numerical studies establish

Three studies separate representation capacity, the cost of finer partitions, and optimizer design. All are classical simulations of quantum algorithms, averaged over ten random initializations—not hardware experiments.

Read chapter
↳

My contribution

  • Co-developed the distributed problem formulation and algorithm.
  • Connected local quantum cost functions with networked classical coordination.
  • Evaluated the construction in numerical simulations, including partition-size and optimizer comparisons.
P

Papers & resources

Distributed Variational Quantum Linear Solver

Tong Shen, Zeru Zhu, and Ji Liu

arXiv:2604.01426 · Submitted to IEEE CDC 2026 · Preprint / Submitted

Last updated