REVIEW 3 major objections 5 minor 1 cited by
Solving 1D Poisson problem with a Variational Quantum Linear Solver
T0 review · 3 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read This paper shows that tridiagonal matrices from 1D Poisson problems can be decomposed into $2^{n-1}+n$ unitary terms using SWAP and center-switch gates, enabling the first VQLS runs on such systems.
desk verdict A concrete small-case decomposition for tridiagonal VQLS, but the general-n claim is asserted, not proved. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the center-switch gate $\mathrm{CS}(n-2)$, a permutation matrix on $n$ qubits that swaps the computational basis states $|011\ldots1\rangle$ and $|100\ldots0\rangle$; for $n=2$ it reduces to the standard SWAP gate. The paper decomposes it as a sequence of $2n-1$ multi-controlled Toffoli gates and uses it, together with Pauli strings and $Z$-type tensor products, to reproduce the tridiagonal structure of the matrix $A$ in the VQLS cost function, reducing the number of terms in the expansion from $2^n$ to $2^{n-1}+n$.
What would settle it
For a $16\times 16$ tridiagonal matrix ($n=4$), explicitly construct the paper's claimed 12 unitary terms, sum them, and compare against the original matrix; any missing or extra entry would refute the counting claim, and the same check can be repeated for $n=5$ with the claimed 21 terms.
Extended reading notes
Core claim
The central claim is that the class of tridiagonal Toeplitz matrices used in 1D Poisson problems admits a unitary decomposition with roughly half as many terms as the standard Pauli-string expansion, by adding multi-qubit permutation gates. The key observation is that the center off-diagonal entries of a tridiagonal matrix are exactly the off-diagonal entries of a SWAP gate for the $4\times4$ case, and of a generalized center-switch gate $\mathrm{CS}(n-2)$ that swaps the bitstrings $011\ldots1$ and $100\ldots0$ for larger $n$. The paper constructs these gates explicitly in terms of Toffoli gates and compensates their diagonal entries with $n$ tensor products of $Z$ operators, leading to the count $2^{n-1}+n$. It then uses this decomposition inside VQLS to solve $2\times2$ and $4\times4$ linear systems, obtaining fidelities close to 1 in simulation and on hardware, and documents the tradeoff between fewer terms and deeper quantum circuits.
Load-bearing premise
The central efficiency claim depends on the unproved assumption that the Pauli strings, SWAP gate, and center-switch gates $\mathrm{CS}(n-2)$ span the space of tridiagonal Toeplitz matrices for every $n$ and that the term count $2^{n-1}+n$ is exact; if the count is off for larger $n$, the advantage over the Pauli decomposition disappears.
Editorial extensions
If this is right
- For a $2^n\times 2^n$ tridiagonal system, the number of unitary terms in the matrix expansion drops from $2^n$ to $2^{n-1}+n$, reducing the number of quantum circuits needed to estimate the VQLS cost function.
- The multi-qubit-gate decomposition extends VQLS beyond few-qubit Hamiltonians to a class of engineering matrices derived from finite-element discretizations.
- The depth penalty is real but not prohibitive at small scale: the $4\times4$ case runs in circuits of depth 99 versus 51, and still reaches high fidelity on hardware.
- The small hardware runs show that even when the cost estimate is biased by sampling noise and decoherence, the optimizer can land at parameters close to the analytical solution, a resilience the authors attribute to noise tolerance of variational algorithms.
Reading between the lines
- The $O(L^2)$ circuit count means the term reduction from $2^n$ to about $2^{n-1}$ should cut the number of cost-circuit evaluations by roughly a factor of four, a consequence the paper does not state explicitly.
- The same permutation-gate construction could be adapted to other banded or structured matrices that have off-diagonal blocks resembling SWAP-like permutations, though the diagonal compensation would need to be re-derived.
- For larger $n$, the center-switch gate requires $2n-1$ multi-controlled Toffoli gates, so one could test whether alternative decompositions of the permutation yield shallower circuits and whether the depth penalty stays acceptable on current hardware.
- If the term count holds at scale, the main obstacle to using this decomposition becomes circuit depth rather than circuit count, shifting the focus to depth-reduction techniques such as the Hadamard-Overlap Test.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper applies the variational quantum linear solver (VQLS) to the one-dimensional Poisson problem, whose discrete matrix is tridiagonal and Toeplitz. It proposes a decomposition of such a matrix into Pauli strings together with SWAP and newly introduced 'center-switch' multi-qubit gates, claiming a reduction in the number of unitary terms from 2^n to 2^{n-1}+n for a 2^n x 2^n matrix. The 4x4 and 8x8 cases are worked out explicitly, and simulated plus real-hardware VQLS results are reported for 2x2 and 4x4 systems.
Significance. If the general-n decomposition is correct, the paper offers a concrete, non-trivial instance of a VQLS decomposition for an engineering-relevant matrix class, with a term-count saving that is exponential in n (from 2^n to 2^{n-1}+n). The trade-off against circuit depth is honestly acknowledged. The 4x4 explicit decomposition is correct and the hardware demonstration, though small, addresses a gap in the VQLS literature, which has mostly focused on physics Hamiltonians. However, the central claim for arbitrary n is asserted rather than proved, and the text contains internal inconsistencies about the gate family and the number of diagonal correction terms. The significance is therefore conditional on a rigorous general construction, which the current manuscript does not supply.
major comments (3)
- [Section III.B.2, Eq. (20)] The text states that 'The set of Pauli gates and their tensor products, together with the SWAP and CS(n−2) gates are all the unitary terms required to decompose a tridiagonal 2n × 2n matrix' (Section III.B.2, final paragraph). This is not correct for n>3: the single gate CS(n−2) defined in Eq. (20) only swaps the central pair (2^{n-1}-1, 2^{n-1}). For n=4, Table II also needs I3CS(2−0), i.e., a center-switch gate on fewer qubits, covering the pairs (3,4) and (11,12). The general construction requires a hierarchy of center-switch gates on k+1 qubits for k=1,...,n−1 (with k=1 being SWAP and k=0 being X on the least significant qubit). This hierarchy is visible in Table II but never stated in the text. Please define the full gate family explicitly and prove that the n gates cover each of the 2^n−1 tridiagonal off-diagonal edges exactly once.
- [Section III.C] The claim that the diagonal correction 'only requires n additional unitaries' is contradicted by Table II: for n=4 the decomposition contains seven non-identity Z strings (I3I2Z1Z0, I3Z2I1Z0, I3Z2Z1I0, Z3I2I1Z0, Z3I2Z1I0, Z3Z2I1I0, Z3Z2Z1Z0). The correct count of non-identity Z strings is 2^{n-1}-1, so the total term count is n + 1 + (2^{n-1}-1) = 2^{n-1}+n, which matches the stated formula. The explanation in the text must be corrected so that it is consistent with Table II and with the final count.
- [Section III.B.2 and III.C] No proof is given that the proposed set of unitaries spans the space of tridiagonal matrices of the form (10) for general n, nor are the expansion coefficients specified for general n. The explicit 4x4 and 8x8 examples are not a substitute, because it is not demonstrated that the same pattern of coefficients and Z corrections works for every n. Please provide a constructive proof (for example, by induction, or by explicitly verifying that the n off-diagonal gates cover all edges of the tridiagonal graph and that the even-parity Z strings fix the diagonal) and state the coefficients in Eq. (3) for arbitrary n.
minor comments (5)
- [Abstract and Section I] Minor language issues: 'less terms' should be 'fewer terms', and 'Proabably' is a typo for 'Probably'.
- [Section III.C] The formula in the text '2(n−1) + n' should read '2^{n-1} + n'; the superscript formatting appears to have been lost in compilation.
- [Section IV] The 51/99 circuit-depth comparison in Section III.C is asserted without details of the compiled circuits; for reproducibility, please include the depth calculation method or a code repository.
- [Section IV] The hardware results for the 4x4 case consist of only two runs. Please state the number of shots and any error bars, and consider adding more repetitions or at least an explicit caveat about the limited statistics.
- [References] Reference [30] is a master's thesis cited as 'to be published'; please provide a stable reference or archive link.
Circularity Check
No significant circularity: the new decomposition is a direct mathematical construction, and the VQLS runs are demonstrations rather than fitted predictions.
full rationale
The paper's derivation chain is non-circular. The central contribution in Section III.B is a constructive linear-algebra statement: a tridiagonal 2^n by 2^n matrix is represented as a linear combination of Pauli strings plus SWAP and center-switch unitaries, with Table II giving explicit term lists for n = 2, 3, 4 and a stated general count of 2^(n-1) + n terms. This count is not obtained by fitting parameters to data and then predicting the same data, nor does any result reduce to its input by construction. The VQLS experiments in Section IV are demonstrations of an algorithm; the reported fidelities are measured outputs of a variational optimization, not theoretical predictions derived from fitted quantities. The only remotely self-referential citation is [30], used for the standard fact that a tridiagonal matrix has only N nonzero Pauli coefficients; that fact is elementary and independently checkable, so the citation is not load-bearing. There is a genuine rigor gap: Section III.B.2 asserts without proof that the set of Pauli strings, SWAP gates, and CS(n-2) gates spans the space of tridiagonal Toeplitz matrices for general n, and the text's description of the gate set is under-specified for n > 3, where Table II appears to require a hierarchy of center-switch gates rather than a single CS(n-2) gate. These are unproved mathematical claims and presentation issues, not circular reasoning. Therefore no circular step is present and the score is 0.
Assumptions & free parameters
assumptions (4)
- standard math Every N by N Hermitian matrix can be expressed as a linear combination of N^2 Pauli strings, and for tridiagonal Toeplitz matrices only N terms are non-zero.
- standard math The center-switch gate CS(n-2) is a unitary permutation matrix and can be decomposed into 2n-1 multi-controlled gates.
- domain assumption The QLSP and VQLS framework, including the Hadamard test and cost function, works as described in Bravo-Prieto et al. [14].
- domain assumption The RY single-qubit ansatz V(theta) can represent the solution state well enough for the solver to converge.
invented entities (1)
-
center-switch gate CS(n-2)
independent evidence
Cite this review
Pith. "Pith review of Solving 1D Poisson problem with a Variational Quantum Linear Solver." pith.science (2026). https://pith.science/paper/YKKN2SIP
@misc{pith2026241204938,
author = {Pith},
title = {Pith review of: Solving 1D Poisson problem with a Variational Quantum Linear Solver},
year = {2026},
howpublished = {\url{https://pith.science/paper/YKKN2SIP}},
note = {Machine review of arXiv:2412.04938}
}
read the original abstract
Different hybrid quantum-classical algorithms have recently been developed as a near-term way to solve linear systems of equations on quantum devices. However, the focus has so far been mostly on the methods, rather than the problems that they need to tackle. In fact, these algorithms have been run on real hardware only for problems in quantum physics, such as Hamiltonians of a few qubits systems. These problems are particularly favorable for quantum hardware, since their matrices are the sum of just a few unitary terms and since only shallow quantum circuits are required to estimate the cost function. However, for many interesting problems in linear algebra, it appears far less trivial to find an efficient decomposition and to trade it off with the depth of the cost quantum circuits. A first simple yet interesting instance to consider are tridiagonal systems of equations. These arise, for instance, in the discretization of one-dimensional finite element analyses. This work presents a method to solve a class of tridiagonal systems of equations with the variational quantum linear solver (VQLS), a recently proposed variational hybrid algorithm for solving linear systems. In particular, we present a new decomposition for this class of matrices based on both Pauli strings and multi--qubit gates, resulting in less terms than those obtained by just using Pauli gates. Based on this decomposition, we discuss the tradeoff between the number of terms and the near-term implementability of the quantum circuits. Furthermore, we present the first simulated and real-hardware results obtained by solving tridiagonal linear systems with VQLS, using the decomposition proposed.
Figures
Forward citations
Cited by 1 Pith paper
-
A generalized variational quantum linear solver on photonic platform
A photonic two-qubit experiment runs VQLS on 4D linear systems, adding Tikhonov regularization for singular cases and a modulo-2 cost function for finite-field equations.
Reference graph
Works this paper leans on
-
[1]
Polynomial-time algorithms for prime factorization and discrete logarithms on a quantum computer,
P. W. Shor, “Polynomial-time algorithms for prime factorization and discrete logarithms on a quantum computer,” SIAM Review , vol. 41, no. 2, pp. 303–332, 1999. [Online]. Available: http: //www.jstor.org/stable/2653075
-
[2]
A fast quantum mechanical algorithm for database search,
L. K. Grover, “A fast quantum mechanical algorithm for database search,” in Proceedings of the Twenty-Eighth Annual ACM Symposium on Theory of Computing , ser. STOC ’96. New York, NY , USA: Association for Computing Machinery, 1996, p. 212–219. [Online]. Available: https://doi.org/10.1145/237814.237866
arXiv 1996
-
[3]
Quantum Computing in the NISQ era and beyond,
J. Preskill, “Quantum Computing in the NISQ era and beyond,” Quan- tum, vol. 2, p. 79, Aug. 2018
2018
-
[4]
Quantum supremacy using a programmable super- conducting processor,
F. Arute et al. , “Quantum supremacy using a programmable super- conducting processor,” Nature, vol. 574, no. 7779, pp. 505–510, Oct
-
[5]
The theory of variational hybrid quantum-classical algorithms,
J. R. McClean, J. Romero, R. Babbush, and A. Aspuru-Guzik, “The theory of variational hybrid quantum-classical algorithms,” New Journal of Physics , vol. 18, no. 2, p. 023023, Feb. 2016. [Online]. Available: https://doi.org/10.1088/1367-2630/18/2/023023
-
[6]
A variational eigenvalue solver on a photonic quantum processor,
A. Peruzzo, J. McClean, P. Shadbolt, M.-H. Yung, X.-Q. Zhou, P. J. Love, A. Aspuru-Guzik, and J. L. O’Brien, “A variational eigenvalue solver on a photonic quantum processor,” Nature Communications , vol. 5, no. 1, Jul. 2014. [Online]. Available: https://doi.org/10.1038/ ncomms5213
work page 2014
-
[7]
Scalable quantum simulation of molecular energies,
P. O’Malley et al. , “Scalable quantum simulation of molecular energies,” Physical Review X , vol. 6, no. 3, Jul. 2016. [Online]. Available: https://doi.org/10.1103/physrevx.6.031007
-
[8]
Variational quantum algorithms for discovering hamiltonian spectra,
T. Jones, S. Endo, S. McArdle, X. Yuan, and S. C. Benjamin, “Variational quantum algorithms for discovering hamiltonian spectra,” Phys. Rev. A , vol. 99, p. 062304, Jun 2019. [Online]. Available: https://link.aps.org/doi/10.1103/PhysRevA.99.062304
Show all 37 references
-
[9]
Quantum machine learning in feature hilbert spaces,
M. Schuld and N. Killoran, “Quantum machine learning in feature hilbert spaces,” Phys. Rev. Lett., vol. 122, p. 040504, Feb 2019. [Online]. Available: https://link.aps.org/doi/10.1103/PhysRevLett.122.040504
2019 doi
-
[10]
Variational quantum circuits for deep reinforcement learning,
S. Y .-C. Chen, C.-H. H. Yang, J. Qi, P.-Y . Chen, X. Ma, and H.-S. Goan, “Variational quantum circuits for deep reinforcement learning,” IEEE Access , vol. 8, pp. 141 007–141 024, 2020
2020
-
[11]
Variational quantum generators: Generative adversarial quantum machine learning for continuous distributions,
J. Romero and A. Aspuru-Guzik, “Variational quantum generators: Generative adversarial quantum machine learning for continuous distributions,” Advanced Quantum Technologies , vol. 4, no. 1, p. 2000003, Dec. 2020. [Online]. Available: https://doi.org/10.1002/qute. 202000003
2020 doi
-
[12]
Near-term quantum algorithms for linear systems of equations,
H.-Y . Huang, K. Bharti, and P. Rebentrost, “Near-term quantum algorithms for linear systems of equations,” 2019. [Online]. Available: arXiv:1909.07344
2019 arXiv
-
[13]
Variational algorithms for linear algebra,
X. Xu, J. Sun, S. Endo, Y . Li, S. C. Benjamin, and X. Yuan, “Variational algorithms for linear algebra,” 2019. [Online]. Available: arXiv:1909.03898
2019 arXiv
-
[14]
Variational quantum linear solver,
C. Bravo-Prieto, R. LaRose, M. Cerezo, Y . Subasi, L. Cincio, and P. J. Coles, “Variational quantum linear solver,” 2020. [Online]. Available: arXiv:1909.05820v2
2020 arXiv
-
[15]
Quantum linear system solver based on time-optimal adiabatic quantum computing and quantum approximate optimization algorithm,
D. An and L. Lin, “Quantum linear system solver based on time-optimal adiabatic quantum computing and quantum approximate optimization algorithm,” 2020. [Online]. Available: arXiv:1909.05500
2020 arXiv
-
[16]
Quantum-classical algorithms for skewed linear systems with an optimized hadamard test,
B. Wu, M. Ray, L. Zhao, X. Sun, and P. Rebentrost, “Quantum-classical algorithms for skewed linear systems with an optimized hadamard test,” Physical Review A , vol. 103, no. 4, Apr. 2021. [Online]. Available: https://doi.org/10.1103/physreva.103.042422
2021 doi
-
[17]
Quantum algorithm for linear systems of equations,
A. W. Harrow, A. Hassidim, and S. Lloyd, “Quantum algorithm for linear systems of equations,” Phys. Rev. Lett. , vol. 103, p. 150502, Oct
-
[18]
Quarteroni, Numerical Models for Differential Problems
A. Quarteroni, Numerical Models for Differential Problems . Springer Milan, 2014. [Online]. Available: https://doi.org/10.1007/ 978-88-470-5522-3
2014
-
[19]
Power of one bit of quantum information,
E. Knill and R. Laflamme, “Power of one bit of quantum information,” Phys. Rev. Lett., vol. 81, pp. 5672–5675, Dec 1998. [Online]. Available: https://link.aps.org/doi/10.1103/PhysRevLett.81.5672
1998 doi
-
[20]
Hardness of classically sampling the one-clean- qubit model with constant total variation distance error,
T. Morimae, “Hardness of classically sampling the one-clean- qubit model with constant total variation distance error,” Phys. Rev. A , vol. 96, p. 040302, Oct 2017. [Online]. Available: https://link.aps.org/doi/10.1103/PhysRevA.96.040302
2017 doi
-
[21]
Impossibility of classically simulating one-clean-qubit model with multiplicative error,
K. Fujii, H. Kobayashi, T. Morimae, H. Nishimura, S. Tamate, and S. Tani, “Impossibility of classically simulating one-clean-qubit model with multiplicative error,” Phys. Rev. Lett. , vol. 120, p. 200502, May
-
[22]
Hardware-efficient variational quantum eigensolver for small molecules and quantum magnets,
A. Kandala, A. Mezzacapo, K. Temme, M. Takita, M. Brink, J. M. Chow, and J. M. Gambetta, “Hardware-efficient variational quantum eigensolver for small molecules and quantum magnets,” Nature, vol. 549, no. 7671, pp. 242–246, Sep. 2017. [Online]. Available: https://doi.org/10.10...
2017 doi
-
[23]
From the quantum approximate optimization algorithm to a quantum alternating operator ansatz,
S. Hadfield, Z. Wang, B. O'Gorman, E. Rieffel, D. Venturelli, and R. Biswas, “From the quantum approximate optimization algorithm to a quantum alternating operator ansatz,” Algorithms, vol. 12, no. 2, p. 34, Feb. 2019. [Online]. Available: https://doi.org/10.3390/a12020034
2019 doi
-
[24]
Quantum approximate optimization is computationally universal,
S. Lloyd, “Quantum approximate optimization is computationally universal,” 2018. [Online]. Available: arXiv:1812.11075
2018 arXiv
-
[25]
Stochastic gradient descent for hybrid quantum-classical optimization,
R. Sweke, F. Wilde, J. Meyer, M. Schuld, P. K. Faehrmann, B. Meynard-Piganeau, and J. Eisert, “Stochastic gradient descent for hybrid quantum-classical optimization,” Quantum, vol. 4, p. 314, Aug
-
[26]
Fletcher, Practical Methods of Optimization; (2nd Ed.)
R. Fletcher, Practical Methods of Optimization; (2nd Ed.) . USA: Wiley- Interscience, 1987
1987
-
[27]
A simplex method for function minimization,
J. A. Nelder and R. Mead, “A simplex method for function minimization,” The Computer Journal , vol. 7, no. 4, pp. 308–313, Jan
-
[28]
M. J. D. Powell, A Direct Search Optimization Method That Models the Objective and Constraint Functions by Linear Interpolation . Dordrecht: Springer Netherlands, 1994, pp. 51–67. [Online]. Available: https://doi.org/10.1007/978-94-015-8330-5 4
1994 doi
-
[29]
M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information: 10th Anniversary Edition . Cambridge University Press, 2010
2010
-
[30]
A variational linear solver for the poisson 1d matrix,
E. Cappanera, “A variational linear solver for the poisson 1d matrix,” Master’s thesis, 2021, to be published
2021
-
[31]
Qiskit: An open-source framework for quantum computing,
H. Abraham et al. , “Qiskit: An open-source framework for quantum computing,” 2019
2019
-
[32]
Ibm quantum,
“Ibm quantum,” 2021. [Online]. Available: https://quantum-computing. ibm.com/
2021
-
[33]
Noise resilience of variational quantum compiling,
K. Sharma, S. Khatri, M. Cerezo, and P. J. Coles, “Noise resilience of variational quantum compiling,” New Journal of Physics, vol. 22, no. 4, p. 043006, Apr. 2020. [Online]. Available: https://doi.org/10.1088/1367-2630/ab784c
2020 doi
-
[1965]
Available: https://doi.org/10.1093/comjnl/7.4.308
[Online]. Available: https://doi.org/10.1093/comjnl/7.4.308
-
[2018]
Available: https://link.aps.org/doi/10.1103/PhysRevLett
[Online]. Available: https://link.aps.org/doi/10.1103/PhysRevLett. 120.200502
-
[2019]
Available: https://doi.org/10.1038/s41586-019-1666-5
[Online]. Available: https://doi.org/10.1038/s41586-019-1666-5
-
[2020]
Available: https://doi.org/10.22331/q-2020-08-31-314
[Online]. Available: https://doi.org/10.22331/q-2020-08-31-314
2020 doi
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.