REVIEW 3 major objections 8 minor 33 references
Optimization by VarQITE on Adaptive Variational Quantum Kolmogorov-Arnold Network
T0 review · 3 major / 8 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read The paper proposes training an Adaptive Variational Quantum Kolmogorov-Arnold Network through Variational Quantum Imaginary Time Evolution (VarQITE), and reports that the pure damping of quantum states—no separate gradient-based…
desk verdict New combination of VarQITE and Adaptive VQKAN, but the target-independent evolution dynamics leave the central claim unsupported. 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 load-bearing object is the VarQITE loss (Eq. 5, summed in Eq. 8), which measures, for each sample, how faithfully the current variational state reproduces one step of imaginary-time damping: $$l_{\mathrm{VarQITE}}^m = \left|1 - \frac{\langle \Psi_\tau | $e^{{-\delta\tau H}}$ | \Psi_{\tau-1} \rangle}{\sqrt{\$\cosh$^2\delta\tau - \$\sinh$^2\delta\tau\,\langle\Psi_{\tau-1}|H|\Psi_{\tau-1}\rangle}}\right|,$$ with the propagator Taylor-expanded (Eq. 6) and reduced to the $\cosh\delta\tau - H\sinh\delta\tau$ identity (Eq. 7) when $H$ is a single Pauli operator. The ansatz being trained is the Adaptive VQKAN of Eqs. (1)–(3), whose per-layer operator strengths $\varphi^{np}_{jk}(n x_m)$ are built from B-spline functions and a Fermi–Dirac-like activation $E_f$; the supervised target $L$ (Eq. 4) is the absolute distance between the Hamiltonian expectation $\langle H\rangle$ and the aimed value. The purpose of the machinery is to replace an external optimizer with a physical contraction: imaginary-time damping drives the state toward lower-energy configurations, and the paper claims this contraction is what moves the readout toward the target function.
What would settle it
Run the same Adaptive VQKAN, ansatz-growth rule, initialization, and COBYLA budget on the five fitting targets, but minimize the supervised loss $L$ of Eq. (4) directly instead of the VarQITE loss of Eq. (8), and compare test-point sums of absolute distances. If the direct-$L$ version matches or beats VarQITE, the claim that imaginary-time damping is the operative training mechanism is falsified; alternatively, exhibiting a single sample where $L_{\mathrm{VarQITE}}$ decreases while $L$ increases would directly contradict the assumed link.
Extended reading notes
Core claim
The paper's claim, stated in its concluding section, is that VarQITE can optimize Adaptive VQKAN, and it presents this as the first practical application of QITE-family methods to quantum machine learning over multiple input samples. In the method, each training step propagates the variational state under the damping operator $e^{-\delta\tau H}$ with $\delta\tau = 0.1$, growing the adaptive ansatz every $\tau = 2.0$ and increasing the spline grid every $\tau = 0.1$, and the parameters are tuned so that the state best emulates one imaginary-time step; the readout Hamiltonian expectation is then compared with the target value. The reported results are that the supervised loss and the sum of absolute test distances decrease as $\tau$ grows, that the minimum distances on the target function of Eq. (9), the logarithmic function, and the fractional function are smaller than the average of a quantum neural network baseline, and that the fractional-function minimum is also smaller than the Enhanced VQKAN comparison quoted from the prior paper. The paper concedes that convergence is slow, that the distances rebound for $\tau > 6$, and that classification errors remain larger than those of ordinary Adaptive VQKAN, attributing the rebound to the difficulty of emulating imaginary-time propagators accurately.
Load-bearing premise
Everything rests on the assumption that pushing the VarQITE fidelity loss (Eq. 8) down also pushes the actual prediction loss $L$ (Eq. 4) down; the paper never proves or analyzes the relationship between these two different objectives, so the fitting claim is supported only by the empirical correlation seen in a few toy runs.
Editorial extensions
If this is right
- VarQITE becomes a viable training mechanism for the Adaptive VQKAN family, meaning quantum machine learning can proceed by imaginary-time damping without special training techniques.
- The same training rule transfers to other VQKAN variants and quantum neural networks, since the damping update does not depend on the specific ansatz form.
- On the target function of Eq. (9), the logarithm, and the fractional function, the minimum test distances are smaller than the average quantum neural network baseline; the fractional function also beats the Enhanced VQKAN comparison quoted from the prior paper.
- Slow convergence and the rebound of absolute distances for $\tau > 6$ are the identified open problems, pointing to the fidelity of emulating imaginary-time propagators as the limiting factor.
- In the 2-D classification setting the damping still reduces the losses, but the errors stay above ordinary Adaptive VQKAN, so more imaginary time is expected to be needed for practical accuracy.
Reading between the lines
- The paper leaves the link between $L_{\mathrm{VarQITE}}$ and the supervised loss $L$ unexamined; a testable extension is to compute both gradients on a single sample and check whether the VarQITE descent direction ever opposes the $L$ descent direction.
- If the claimed mechanism holds, the natural conjecture it implies is that any quantum machine learning ansatz whose loss is a Hamiltonian expectation could be trained by imaginary-time damping, and that conjecture is testable by swapping the Adaptive VQKAN ansatz for a plain hardware-efficient ansatz on the same five targets.
- The observed rebound at $\tau > 6$ is attributed to imperfect emulation of the propagator; flipping the sign of $H$ or replacing VarQITE with probabilistic imaginary time evolution on the same schedule would isolate whether the rebound is a property of the damping dynamics or of the variational emulation.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes training an Adaptive Variational Quantum Kolmogorov-Arnold Network (VQKAN) by Variational Quantum Imaginary Time Evolution (VarQITE). The method updates ansatz parameters by minimizing a VarQITE fidelity loss under a fixed Hamiltonian, with the readout given by the expectation value of the same Hamiltonian. Simulations are reported for regression of five elementary functions and for a two-dimensional classification task. The authors observe decreasing supervised losses and conclude that VarQITE can optimize Adaptive VQKAN, while acknowledging slow convergence, rebounds, poor classification results, and a possible sign problem. The claimed accuracy advantage over a quantum neural network is based on a previous paper rather than on comparisons in this manuscript.
Significance. The idea of using imaginary-time damping as a training signal for quantum machine learning is interesting and, if it worked as described, would offer a non-gradient optimization alternative for small-parameter ansatze. The paper is candid about its limitations and provides enough protocol detail to identify the key equations. However, the manuscript does not establish that minimizing the VarQITE loss reduces the supervised loss, and Eq. (5) is algebraically incorrect for the multi-term Hamiltonian actually used. The empirical results are suggestive but lack error bars for most curves, direct baseline comparisons, and code. As written, the central claim is therefore not yet supported.
major comments (3)
- [Section 2, Eqs. (4)-(8)] The central claim that VarQITE optimizes Adaptive VQKAN is not supported because no relationship is established between L_VarQITE and the supervised loss L. With the Hamiltonian H=Z0Z1+Z2Z3 used in Sec. 3.1, exact imaginary-time evolution obeys d<H>/dτ = -2 Var_Ψ(H) ≤ 0, so every input state's readout expectation moves monotonically toward the ground-state expectation of the same fixed H. Since the targets in Eq. (9) range up to e^2≈7.39 while <H>∈[-2,2], the supervised loss |<H>-f_aim| must increase for target values above the current expectation. The sample weights a_m depend on f_aim, but they only rescale the per-sample VarQITE losses; they do not make the evolution target-directed. The observed decrease of L in Sec. 3 therefore requires an additional mechanism (e.g., trainable θ_j in H, a target-dependent H, or optimizing L directly instead of L_VarQITE) that is not described. Please provide an analysis of when minimizing L_VarQITE reduces L, or modify the algorithm so that the evolution is guided by the targets.
- [Eq. (5) and (7)] The normalization in Eq. (5) is incorrect for the Hamiltonian used in this paper. The identity e^{-δτH}=cosh δτ - H sinh δτ and the denominator sqrt(cosh^2 δτ - sinh^2 δτ <H>) apply only when H^2=I. For H=Z0Z1+Z2Z3, H^2=2I+2Z0Z1Z2Z3, so Eq. (5) does not compute the fidelity between Ψ_τ and the imaginary-time-propagated state. Consequently the quantity actually minimized in the simulations is not the stated VarQITE loss. The authors should either restrict H to a single Pauli operator or derive the correct normalization for the multi-term Hamiltonian.
- [Section 3.1 and Abstract] The abstract's claim that 'our method was more accurate than a quantum neural network in solving some problems' is not supported by any comparison in this manuscript. The only evidence is the statement that, according to the previous paper [28], the minimum test losses for Eq. (9), the logarithmic function, and the fractional function are smaller than the average of a Quantum Neural Network. No baseline results, experimental settings, or statistical tests are shown here, and the classification experiments are explicitly reported to be worse than ordinary Adaptive VQKAN. Please include direct comparisons in the present paper or qualify the accuracy claim to what is actually demonstrated.
minor comments (8)
- [Section 1] The sentence 'VQKAN is the quantum version of VQKAN' should read 'quantum version of KAN'.
- [Section 2] The sentence 'We use blueqat SDK30 for numerical simulation of quantum calculations and COBYLA of scipy to optimize parameters but to declare the use of others' is incomplete and should be rephrased.
- [Eq. (8)] The phrase 'the same as that of eq. 5' is incorrect because Eq. (5) does not define a_m; the cross-reference should be to Eq. (4) or to the surrounding text.
- [Section 2] The notation 'nxm' and '1xm' is confusing; the dimension of the input vector and the meaning of the superscript m should be defined once in the text.
- [Figures 1-3] Most curves have no error bars and it is not stated whether they are single runs or averages; please specify this and label the y-axes clearly.
- [Section 3] The 'sum of absolute distances' on test points is used as a metric but never formally defined; please define this test loss in Section 2 or 3.
- [Eq. (2)] The argument of 'acos' may leave the interval [-1,1] for arbitrary parameter values; please explain how this is handled in the implementation.
- [Section 1] The paper cites [25] as the Adaptive VQKAN paper but does not clearly state what is new in the present work beyond replacing the optimizer; please add an explicit statement of the incremental contribution.
Circularity Check
No circularity in the core VarQITE-VQKAN derivation; the weak links are missing proof and technical normalization errors, not by-construction equivalence.
full rationale
The derivation chain is not circular. In Section 2, Eq. (4) defines the supervised loss L as the absolute distance between the readout expectation <H> and the target f_aim, while Eqs. (5)-(8) define a separate VarQITE objective that tracks imaginary-time propagation under H. Parameters are optimized with COBYLA against L_VarQITE; the predicted values are not obtained by fitting or renaming the training targets. The target enters only through the sample weights a_m in Eq. (8), and L_VarQITE is not algebraically identical to L, so the reported test loss is not forced by construction. The fixed-H concern is a real missing-mechanism/correctness gap rather than circularity: the paper does not define H in terms of f_aim or relabel a fitted parameter as a prediction. Eq. (5)'s denominator assumes H^2=I and is invalid for H=Z0Z1+Z2Z3; this is a technical flaw, not a circular reduction. Self-citations [25] and [28] supply the adaptive ansatz and baseline numbers, but they are not used as an unverified uniqueness theorem, and the baseline comparisons are externally checkable. The manuscript itself concedes limitations (weak classification, possible sign problem, rebound for tau>6), which weaken the empirical claim but do not show that the prediction is equivalent to the input by construction.
Assumptions & free parameters
free parameters (7)
- imaginary time step delta_tau =
0.1
- initial ansatz hyperparameters =
N=10, Nl=1, Nq=4, Ng=8, Ns=8
- ansatz growth and grid schedule =
grow ansatz per 2.0 tau, increase Ns by 4 per 0.1 tau
- initial random perturbation =
uniform 0 to 0.1
- Hamiltonian for fitting =
Z0Z1 + Z2Z3
- convergence threshold =
1e-16
- sample weighting coefficient a_m =
(1 - f_aim + 1/N)/(2 + 1/N)
assumptions (5)
- standard math The Taylor expansion of the imaginary-time propagator converges, and for a single Pauli H, e^{-delta_tau H} = cosh(delta_tau) - H sinh(delta_tau).
- domain assumption The parameterized quantum state can represent the target functions with sufficient accuracy.
- ad hoc to paper Minimizing the VarQITE loss L_VarQITE reduces the supervised loss L.
- domain assumption Swap test and Hadamard test can estimate the overlap in Eq. (5) without error, and COBYLA converges to the optimum.
- ad hoc to paper The adaptive ansatz growth rule preserves the validity of the VarQITE update.
Cite this review
Pith. "Pith review of Optimization by VarQITE on Adaptive Variational Quantum Kolmogorov-Arnold Network." pith.science (2026). https://pith.science/paper/BSSWKLWZ
@misc{pith2026250622801,
author = {Pith},
title = {Pith review of: Optimization by VarQITE on Adaptive Variational Quantum Kolmogorov-Arnold Network},
year = {2026},
howpublished = {\url{https://pith.science/paper/BSSWKLWZ}},
note = {Machine review of arXiv:2506.22801}
}
read the original abstract
Quantum imaginary time evolution (QITE) is a powerful method to derive the ground states of the systems. Only the damping of quantum states leads it; hence, reaching the ground state is guaranteed by nature without any external manipulation. Numerous QITE methods by many groups are used to improve speed and accuracy, derive excited states, and solve combined optimization problems. However, the QITE methods have not been used for quantum machine learning to predict the ideal values for multiple input values. Therefore, we propose a method for applying QITE methods for quantum machine learning and demonstrate fitting problems of elementary functions and classification problems on a 2-D plane. As a result, we confirmed that our method was more accurate than a quantum neural network in solving some problems. Our method can be used for other quantum machine learning algorithms; hence, it may be the milestone for applying QITE to quantum machine learning.
Figures
Figures from the paper (5 more)
Reference graph
Works this paper leans on
-
[28]
Wakaura, H., Mulyawan, R. & Suksmono, A. B. Enhanced Variational Quantum Kolmogorov-Arnold Network. arXiv e-prints arXiv:2503.22604, DOI: 10.48550/arXiv.2503.22604 (2025). 2503.22604
-
[1]
Feynman, R. P. Simulating physics with computers. Int. J. Theor. Phys. 21, 467–488, DOI: 10.1007/BF02650179 (1982)
- [2]
-
[3]
Algorithms for quantum computation: discrete logarithms and factoring
Shor, P. Algorithms for quantum computation: discrete logarithms and factoring. In Proceedings 35th Annual Symposium on F oundations of Computer Science, 124–134, DOI: 10.1109/SFCS.1994.365700 (1994)
arXiv 1994
-
[4]
Mezzacapo, A., Casanova, J., Lamata, L. & Solano, E. Topological qubits with majorana fermions in trapped ions. New J. Phys. 15, 033005, DOI: 10.1088/1367-2630/15/3/033005 (2013)
-
[5]
Dutt, M. V . G.et al. Quantum register based on individual electronic and nuclear spin qubits in diamond. Science 316, 1312–1316, DOI: 10.1126/science.1139831 (2007). https://www.science.org/doi/pdf/10.1126/science.1139831
-
[6]
Wu, X. et al. A concise review of Rydberg atom based quantum computation and quantum simulation. arXiv e-prints arXiv:2012.10614 (2020). 2012.10614. 7/9
arXiv 2020
-
[7]
D., Perdomo-Ortiz, A., Yung, M.-H
Kassal, I., Whitfield, J. D., Perdomo-Ortiz, A., Yung, M.-H. & Aspuru-Guzik, A. Simulating chemistry using quantum computers. Annu. Rev. Phys. Chem. 62, 185–207, DOI: 10.1146/annurev-physchem-032210-103512 (2011). https: //doi.org/10.1146/annurev-physchem-032210-103512
Show all 33 references
-
[8]
R., Romero, J., Babbush, R
McClean, J. R., Romero, J., Babbush, R. & Aspuru-Guzik, A. The theory of variational hybrid quantum-classical algorithms. New J. Phys. 18, 023023, DOI: 10.1088/1367-2630/18/2/023023 (2016)
2016 doi
-
[9]
R., Economou, S
Grimsley, H. R., Economou, S. E., Barnes, E. & Mayhall, N. J. An adaptive variational algorithm for exact molecular simulations on a quantum computer. Nat. Commun. 10, 3007, DOI: 10.1038/s41467-019-10988-2 (2019). 1812.11173
2019 arXiv
-
[10]
M., Hohenstein, E
Parrish, R. M., Hohenstein, E. G., McMahon, P. L. & Martinez, T. J. Hybrid Quantum/Classical Derivative Theory: Analytical Gradients and Excited-State Dynamics for the Multistate Contracted Variational Quantum Eigensolver. arXiv e-prints arXiv:1906.08728 (2019). 1906.08728
2019 arXiv
-
[11]
Khoshaman, A. et al. Quantum variational autoencoder. Quantum Sci. Technol. 4, 014001, DOI: 10.1088/2058-9565/aada1f (2019). 1802.05779
2019 arXiv
-
[12]
Supervised learning with quantum-enhanced feature spaces
Havlíˇcek, V .et al. Supervised learning with quantum-enhanced feature spaces. Nature 567, 209–212, DOI: 10.1038/ s41586-019-0980-2 (2019). 1804.11326
2019 arXiv
-
[13]
& Rosenkranz, M
Benedetti, M., Coyle, B., Fiorentini, M., Lubasch, M. & Rosenkranz, M. Variational Inference with a Quantum Computer. Phys. Rev. Appl. 16, 044057, DOI: 10.1103/PhysRevApplied.16.044057 (2021). 2103.06720
2021 arXiv
-
[14]
Abel, S., Criado, J. C. & Spannowsky, M. Completely quantum neural networks. Phys. Rev. A 106, 022601, DOI: 10.1103/PhysRevA.106.022601 (2022). 2202.11727
2022 arXiv
-
[15]
T., Ashida, Y
Wang, Z. T., Ashida, Y . & Ueda, M. Deep Reinforcement Learning Control of Quantum Cartpoles.Phys. Rev. Lett. 125, 100401, DOI: 10.1103/PhysRevLett.125.100401 (2020). 1910.09200
2020 arXiv
-
[16]
Quantum Distributed Deep Learning Architectures: Models, Discussions, and Applications
Kwak, Y .et al. Quantum Distributed Deep Learning Architectures: Models, Discussions, and Applications. arXiv e-prints arXiv:2202.11200 (2022). 2202.11200
2022 arXiv
-
[17]
Yang, D., Xiao, Z. & Yu, W. Boosting the Adversarial Transferability of Surrogate Model with Dark Knowledge.arXiv e-prints arXiv:2206.08316 (2022). 2206.08316
2022 arXiv
-
[18]
& Fujii, K
Mitarai, K., Negoro, M., Kitagawa, M. & Fujii, K. Quantum circuit learning. Phys. Rev. A 98, 032309, DOI: 10.1103/ PhysRevA.98.032309 (2018)
2018
-
[19]
Motta, M. et al. Determining eigenstates and thermal states on a quantum computer using quantum imaginary time evolution. Nat. Phys. 16, 205–210, DOI: 10.1038/s41567-019-0704-4 (2019)
2019 doi
-
[20]
McArdle, S. et al. Variational ansatz-based quantum simulation of imaginary time evolution. npj Quantum Inf. 5, 75, DOI: 10.1038/s41534-019-0187-2 (2019)
2019 doi
-
[21]
& Anisimovas, E
Ledinauskas, E. & Anisimovas, E. Scalable imaginary time evolution with neural network quantum states. SciPost Phys. 15, DOI: 10.21468/scipostphys.15.6.229 (2023)
2023 doi
-
[22]
& Fiorentini, M
Amaro, D., Rosenkranz, M., Fitzpatrick, N., Hirano, K. & Fiorentini, M. A case study of variational quantum algorithms for a job shop scheduling problem. EPJ Quantum Technol. 9, DOI: 10.1140/epjqt/s40507-022-00123-4 (2022)
2022 doi
-
[23]
& Matsushita, Y .-i
Kosugi, T., Nishiya, Y ., Nishi, H. & Matsushita, Y .-i. Imaginary-time evolution using forward and backward real-time evolution with a single ancilla: First-quantized eigensolver algorithm for quantum chemistry. Phys. Rev. Res. 4, 033121, DOI: 10.1103/PhysRevResearch.4.033121...
2022 arXiv
-
[24]
& Shingu, Y
Kuji, H., Nikuni, T. & Shingu, Y . Variational quantum-neural hybrid imaginary time evolution. arXiv e-prints arXiv:2503.22570, DOI: 10.48550/arXiv.2503.22570 (2025). 2503.22570
2025 doi
-
[25]
& Suksmono, A
Wakaura, H., Mulyawan, R. & Suksmono, A. B. Adaptive Variational Quantum Kolmogorov-Arnold Network. arXiv e-prints arXiv:2503.21336 (2025). 2503.21336
2025 arXiv
-
[26]
& Mulyawan, R
Wakaura, H., Bayu Suksmono, A. & Mulyawan, R. Variational quantum kolmogorov-arnold network. Res. Sq. DOI: 10.21203/rs.3.rs-4504342/v3 (2024). PREPRINT (Version 3)
2024 doi
- [27]
-
[29]
Dalton, K. et al. Quantifying the effect of gate errors on variational quantum eigensolvers for quantum chemistry. npj Quantum Inf. 10, 18, DOI: 10.1038/s41534-024-00808-x (2024). 2211.04505. 8/9
2024 arXiv
-
[30]
https://github.com/Qaqarot
Kato, T. https://github.com/Qaqarot. Opensource software development kit (2018)
2018
-
[31]
Huggins, W. J. et al. Unbiasing fermionic quantum Monte Carlo with a quantum computer. Nature 603, 416–420, DOI: 10.1038/s41586-021-04351-z (2022). 2106.16235
2022 arXiv
- [32]
-
[33]
& Matsushita, Y .-i
Huang, X., Nishi, H., Kosugi, T., Kawada, Y . & Matsushita, Y .-i. A probabilistic imaginary-time evolution quantum algorithm for advection-diffusion equation: Explicit gate-level implementation and comparisons to quantum linear system algorithms. arXiv e-prints arXiv:2409.185...
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.