REVIEW 4 major objections 4 minor 19 references
The paper claims that a single-layer complex-valued neural network trained on the Stiefel manifold can learn the unitary of a quantum circuit and reach fidelity 1 on the tested 2-qubit and 5-qubit examples.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
A single-layer complex-valued neural network constrained to the Stiefel manifold via Cayley updates is used to learn and transpile quantum circuit unitaries, with fidelity reaching 1 on reported toy examples.
T0 review reviewed 2026-08-05 challenge →
load-bearing objection Standard optimizer applied to a new domain, with a local-decrease proof overclaimed as global and an experiment too underspecified to support the central claim. the 4 major comments →
Quantum Circuit Design using Complex valued Neural Network in Stiefel Manifold
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
Core claim
The paper's central claim is that solving the constrained problem min f(W) subject to W†W=I by gradient descent on the Stiefel manifold preserves unitarity throughout training and produces a learned unitary that matches the target circuit. With W as the network weight matrix and G the Euclidean gradient, the paper constructs A=GW†−WG†, which is skew-Hermitian, and updates W(k+1)=(I+λA/2)^{-1}(I−λA/2)W(k). The appendix proves that each update stays on the Stiefel manifold and that, for small enough λ, the loss strictly decreases unless the Riemannian gradient vanishes. On a 5-qubit sample circuit the unitary error ‖WW†−I‖² reaches about 10⁻¹¹, and the fidelity F=(1/d)|Tr(U_true† U_learned)| r
What carries the argument
The central object is the Stiefel manifold V_N(C^N)={W∈C^{N×N}: W†W=I}, the set of unitary matrices. The carrying mechanism is the Cayley update W(k+1)=(I+λA/2)^{-1}(I−λA/2)W(k) with A=GW†−WG†, where G is the Euclidean gradient of the loss. Because A is skew-Hermitian, the Cayley factor is unitary and maps one unitary matrix to another; the proof shows the curve's tangent at λ=0 is −AW and the directional derivative of f is −‖AW‖²_c, so the loss descends along the manifold. This is what guarantees unitarity at every step rather than only after intermittent re-orthogonalization.
Load-bearing premise
The paper assumes the input-output state pairs used for training uniquely determine the target unitary, so that reaching fidelity 1 on those pairs means the learned circuit is the target circuit.
What would settle it
Take a fixed 2-qubit circuit U_true and provide only one state pair, say |00⟩→U_true|00⟩, to the network. Train until fidelity F=1 on that pair, then compute F on the held-out basis states |01⟩, |10⟩, and |11⟩. If the held-out fidelity is measurably below 1 because many unitaries agree on |00⟩, the paper's claim that F=1 implies U_learned=U_true fails in that setting.
If this is right
- Because every updated weight matrix is unitary, the learned matrix is always physically valid as a quantum gate and can be handed to a transpiler to produce an executable circuit.
- The loss decreases monotonically for sufficiently small learning rates, so training avoids the loss spikes that intermittent Gram-Schmidt orthogonalization introduces.
- If the fidelity result generalizes, a user would only need known input-output quantum state pairs to obtain a circuit, bypassing manual gate decomposition.
- The same update rule could be applied to any complex-valued neural network that needs a unitary weight matrix, not just quantum circuit synthesis.
Where Pith is reading between the lines
- The paper leaves open how many input-output state pairs are used and does not test on held-out states; the claim that reaching fidelity 1 means reproducing the target circuit depends on the training set uniquely determining the unitary, which is not demonstrated.
- A natural extension is to train on a single basis-state pair for a known random unitary and then evaluate fidelity on the remaining basis states; if many unitaries fit the one pair, the learned matrix will not match the target on the held-out states.
- The method could be applied to learning unknown processes: given experimentally measured input-output state pairs, the same Stiefel-constrained optimization would return the unitary closest to the data, with unitarity acting as a built-in regularizer.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes learning a quantum circuit's unitary transformation by training a single-layer complex-valued neural network whose weight matrix is constrained to the Stiefel manifold (i.e., unitary) via the Cayley update of Eq. (4). The authors claim this update preserves unitarity at every step and guarantees a globally decreasing loss, and they report numerical experiments on a 5-qubit circuit and a 2-qubit entangling operation, with training loss decreasing and fidelity (Eq. 5) reaching 1, after which the learned unitary is transpiled into a quantum circuit. The appendix contains proofs of unitarity preservation and of a local decrease property for sufficiently small learning rates.
Significance. If fully validated, the proposal would be a simple baseline for learning unitary maps from input--output quantum state pairs, with the useful property that the learned weight matrix is exactly unitary throughout training. The mathematical appendix correctly proves the key feasibility of the Cayley/Stiefel update: the update (A7) stays on the Stiefel manifold, and the directional derivative at λ=0 is non-positive. The method also has the practical advantage of being directly expressible in existing transpilation toolchains. However, the central empirical claim—that the learned Wfinal implements the target quantum circuit—is not established by the reported experiments, and the claimed global monotone decrease is stronger than what the proof supports. These issues are load-bearing for the paper's main contribution.
major comments (4)
- [§4 and Appendix, Claim 2] The paper states as a guaranteed property that "the loss function decreases globally, i.e., f(W(k+1)) ≤ f(W(k))" for the update (4). The appendix proof, Eqs. (A12)–(A15), only shows that if the Riemannian gradient AW is nonzero, then for sufficiently small λ one has f(Y(λ)) < f(W). At a stationary point (AW=0) the derivative is zero and no decrease is proved. Moreover, no rule is given for choosing λ that satisfies the bound (A15). The monotone-decrease claim used to motivate the method is therefore not established; please replace it with the local statement and provide a concrete step-size procedure or line search.
- [§2, §5, Fig. 2] The empirical claim that the trained Wfinal reproduces the target unitary Utrue is not supported by the experiments as reported. A training set of k pairs (|ψ_i>, Utrue|ψ_i>) determines a unitary only on the span of the input states; for k < 2^n or for linearly dependent inputs, infinitely many unitaries fit the training data. Section 2 only says "many input and output quantum states are known," and Section 5 reports a 5-qubit experiment without stating the number of training pairs, their distribution, or their linear independence. Figure 2b shows fidelity reaching 1, but no held-out fidelity on states not used in training is reported. Please specify the training data and add held-out evaluation before the fidelity claim can be assessed.
- [§5, Eqs. (1)–(4)] The loss f(W) and its Euclidean gradient G are never concretely specified for the reported experiments. Equation (1) defines the constrained optimization abstractly, but the experiments require a concrete loss—mean squared error over output states, trace fidelity, or another—to compute G in Eq. (3) and the update (4). The paper also omits the learning-rate schedule, initialization, number of epochs, and the exact 5-qubit target circuit. Without these details, the numerical results are not reproducible and the fidelity curves cannot be interpreted.
- [Appendix, 'Additional Result'] The 2-qubit example in the appendix is also under-specified: the target entangling operation is not named, the training data are not described, and Figure 4 shows a transpiled circuit without any comparison to the target or a correctness metric beyond the stated unitary error. This example should either be integrated into a fully specified experimental section or removed.
minor comments (4)
- [References] Reference [1] is listed as "Anonymous" and reference [6] contains "and 1 others"; these are incomplete citations and should be corrected.
- [Throughout] There are numerous typos: "varius," "ouput," "Caley" (should be Cayley), "uniatry," and "Steifel" in figure captions. Please proofread carefully.
- [§5, Fig. 2] The axes in Figures 2a and 2b are not fully described. In particular, the fidelity plot should state whether the fidelity is evaluated on the training data or on a separate test set, and the loss plot should state which loss function is displayed.
- [§1, Literature Study] The claim that this is "one of the first applications of Cayley-transform-based manifold optimization in quantum computing" is vague, especially given reference [17] on unitary neural networks for quantum circuit design. Please clarify the precise novelty with a more systematic comparison to prior work.
Circularity Check
No construction-level circularity; the Cayley/Stiefel derivation is self-contained and the reported F=1 is an in-sample training metric, not a circular reduction.
full rationale
The paper's central derivation is the Cayley-transform update on the Stiefel manifold (Eq. 4) with the proof in the Appendix that W remains unitary (Claim 1) and that the loss decreases for small learning rates (Claim 2). These proofs are self-contained and rely only on standard Riemannian optimization facts; they do not depend on the paper's conclusions and there are no self-citations used as load-bearing premises. The fidelity F=1 reported in Fig. 2b is computed between the learned unitary U_learned (the optimized weight matrix W) and the target unitary U_true (Eq. 5). This is an in-sample training evaluation, not an out-of-sample prediction: the paper does not claim to predict held-out states, and it does not present F as an independent confirmation. The main weakness is that Sections 2 and 5 do not specify the number of training input-output pairs, the loss function, or a train/test split, so the identifiability of the target unitary from the training data is not demonstrated. This is a correctness/validation gap, not a circular reduction: no equation equates F to the training objective by construction, and no load-bearing step reduces to its own inputs. Therefore no significant circularity is identified.
Axiom & Free-Parameter Ledger
free parameters (2)
- Learning rate lambda =
unspecified (positive, sufficiently small)
- Initial unitary W(0) =
random, seed unspecified
axioms (4)
- standard math The loss f(W) is a smooth real-valued function of the complex weight matrix with a well-defined Euclidean gradient.
- domain assumption A = G W dagger - W G dagger gives the Riemannian gradient direction under the canonical inner product (A3).
- ad hoc to paper The input-output state pairs used in training are sufficient to identify the target unitary U.
- domain assumption Qiskit's transpile function can convert the learned unitary into an executable quantum circuit.
Cite this review
Pith. "Pith review of Quantum Circuit Design using Complex valued Neural Network in Stiefel Manifold." pith.science (2026). https://pith.science/paper/4TB3FQFO
@misc{pith2026250902374,
author = {Pith},
title = {Pith review of: Quantum Circuit Design using Complex valued Neural Network in Stiefel Manifold},
year = {2026},
howpublished = {\url{https://pith.science/paper/4TB3FQFO}},
note = {Machine review of arXiv:2509.02374}
}
read the original abstract
Quantum algorithms operate on quantum states through unitary transformations in high dimensional complex Hilbert space. In this work, we propose a machine learning approach to create the quantum circuit using a single-layer complex-valued neural network. The input and ouput quantum states are provided to the network, which is trained to approximate the output state of a given quantum algorithm. To ensure that the fundamental property of unitarity is preserved throughout the training process, we employ optimization in Stiefel Manifold.
Figures
Reference graph
Works this paper leans on
-
[1]
Anonymous. 2021. The bitter truth about gate-based quantum algorithms in the nisq era. Quantum In- formation Processing (IOP) . Discusses how lim- ited qubit connectivity forces additional gates and increases circuit depth
work page 2021
-
[2]
Martin Arjovsky, Amar Shah, and Yoshua Bengio
-
[3]
Bennett, Richard Cleve, David P
Adriano Barenco, Charles H. Bennett, Richard Cleve, David P. DiVincenzo, Norman Margolus, Peter Shor, Tycho Sleator, John A. Smolin, and Harald Wein- furter. 1995. Elementary gates for quantum computa- tion. Physical Review A, 52(5):3457
work page 1995
-
[4]
Olivier Coulaud, Luc Giraud, and Martina Iannacito. 2022. On some orthogonalization schemes in tensor train format. arXiv preprint arXiv:2211.08770. Available at https://arxiv. org/abs/2211.08770
work page internal anchor Pith review Pith/arXiv arXiv 2022
-
[5]
David Deutsch and Richard Jozsa. 1992. Rapid so- lution of problems by quantum computation. Pro- ceedings of the Royal Society of London. Series A: Mathematical and Physical Sciences, 439(1907):553– 558
work page 1992
-
[6]
Yan Ge, Wu Wenjie, Chen Yuheng, and 1 others
-
[7]
Houssein, Zainab Abohashima, Mohamed Elhoseny, and Waleed M
Essam H. Houssein, Zainab Abohashima, Mohamed Elhoseny, and Waleed M. Mohamed. 2022. Machine learning in the quantum realm: The state-of-the-art, challenges, and future vision. Expert Systems with Applications, 194:116512
work page 2022
-
[8]
Zakaria Mhammedi, Andrew Hellicar, Ashfaqur Rahman, and James Bailey. 2016. Efficient or- thogonal parametrisation of recurrent neural net- works using householder reflections. arXiv preprint arXiv:1612.00188
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[9]
Michael A. Nielsen and Isaac L. Chuang. 2002. Quantum Computation and Quantum Information . Cambridge University Press
work page 2002
-
[10]
Edward Parker. 2025. Unitary time evolution in quantum mechanics is a stronger physical postulate than linear time evolution. F oundations of Physics, 55:1–13
work page 2025
-
[11]
Rohit Sarma Sarkar and Bibhas Adhikari. 2024. A quantum neural network framework for scalable quantum circuit approximation of unitary matrices. Accessible at arXiv:2405.00012
work page internal anchor Pith review Pith/arXiv arXiv 2024
-
[12]
Vivek V . Shende, Stephen S. Bullock, and Igor L. Markov. 2006. Synthesis of quantum-logic circuits. IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems, 25(6):1000–1010
work page 2006
-
[13]
Peter W. Shor. 1994. Polynomial-time algorithms for prime factorization and discrete logarithms on a quantum computer. In Proceedings of the 35th Annual Symposium on F oundations of Computer Sci- ence, pages 124–134. IEEE
work page 1994
-
[14]
Hemant D. Tagare. 2011. Notes on optimization on stiefel manifolds. Technical report, Yale University
work page 2011
-
[15]
Kwok Ho Wan, Oscar Dahlsten, Hlér Kristjánsson, Robert Gardner, and M. S. Kim. 2016. Quantum generalisation of feedforward neural networks. arXiv preprint arXiv:1612.01045
work page internal anchor Pith review Pith/arXiv arXiv 2016
-
[16]
Hershey, Jonathan Le Roux, and Les Atlas
Scott Wisdom, Thomas Powers, John R. Hershey, Jonathan Le Roux, and Les Atlas. 2016. Full- capacity unitary recurrent neural networks. In Ad- vances in Neural Information Processing Systems , pages 4880–4888
work page 2016
-
[17]
M. Zomorodi, H. Amini, M. Abbaszadeh, J. Sohrabi, V . Salari, and P. Plawiak. 2024. Optimal quantum circuit design via unitary neural networks. arXiv preprint arXiv:2408.13211. Appendix Optimization in Stiefel Manifold We consider the Stiefel manifold of all N × N complex-valued matrices whose columns are N orthonormal vectors in CN. It is defined as VN (...
Pith/arXiv arXiv 2024
-
[2016]
In Proceedings of the 33rd International Conference on Machine Learning (ICML) , pages 1120–1128
Unitary evolution recurrent neural networks. In Proceedings of the 33rd International Conference on Machine Learning (ICML) , pages 1120–1128
-
[2024]
Quantum circuit synthesis and compilation optimization: Overview and prospects. Preprint (arXiv). Provides a survey on translating algorithms into optimized gate sequences under NISQ con- straints
This paper was first reviewed by deepseek-v4-flash on August 5, 2026.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.