Pith. sign in

REVIEW 3 major objections 6 minor 19 references

For amplitude-embedded quantum classifiers, input gradients reduce to overlaps that a Hadamard-test circuit can measure directly, enabling integrated-gradients attribution on quantum hardware.

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 →

HattriQ computes input-feature attributions for amplitude-encoded quantum classifiers by estimating amplitude gradients with Hadamard-test circuits and integrating them from a baseline image.

T0 review reviewed 2026-08-04 challenge →

load-bearing objection A useful Hadamard-test gradient circuit for amplitude-embedded QML, but the pixel-level attribution claims don't follow from the equations because the chain rule through normalization and overflow encoding is never derived. the 3 major comments →

arxiv 2510.02497 v2 pith:BQBIEPQQ submitted 2025-10-02 quant-ph

HattriQ: Designing Integrated Gradients for Feature Attribution in Quantum Machine Learning

classification quant-ph PACS 03.67.Lx
keywords integrated gradientsquantum machine learningamplitude embeddingHadamard testfeature attributionquantum gradientsparameterized quantum circuits
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper tries to establish that integrated-gradients feature attribution—a standard classical interpretability tool—can be computed directly on quantum hardware for quantum classifiers that encode inputs as amplitudes. It shows that the gradient of the model output with respect to an input amplitude equals twice the real part of a simple inner product, and that a Hadamard-test circuit can estimate that inner product by measuring a single ancilla qubit. If correct, HattriQ gives every amplitude-embedded QML model a hardware-native explanation method: one circuit per input feature, with no need to simulate the state vector or access hidden quantum states. The paper validates the method on several image-classification benchmarks, showing attributions that highlight class-relevant pixels and remain stable even at low measurement-shot counts.

Core claim

The central claim is that, for a quantum model F(x;θ)=⟨x|U†(θ)OU(θ)|x⟩ with amplitude-embedded input |x⟩=∑ x_k |b_k⟩, the input gradient takes the closed form ∂F/∂c_k = 2 Re[⟨b_k|U†OU|x⟩] (and similarly for imaginary parts). Because this is an overlap between a computational basis state and the evolved input state, it can be evaluated by a Hadamard-test circuit: prepare a controlled superposition of |x⟩ and |b_k⟩, apply U†OU under control, and measure the ancilla. The paper further shows how to parallelize across features using multiple ancilla qubits and applies the construction to compute integrated-gradients attributions on standard handwritten-digit and fashion benchmarks. The authors fr

What carries the argument

The load-bearing object is the Hadamard-test gradient circuit. It entangles an ancilla with the data register in a superposition of the encoded input V(x)|0⟩ and a basis state V(b_k)|0⟩, then applies U†OU conditioned on the ancilla. The probability of measuring the ancilla in |0⟩ is 1/2(1+Re[⟨b_k|U†OU|x⟩]), which, by Lemma 3.1, is an affine rescaling of the k-th feature gradient. A multi-ancilla generalization encodes several gradients into the joint ancilla measurement probabilities, which are then recovered by solving a linear system; this is the paper's parallelization mechanism.

Load-bearing premise

The pixel-level attributions shown assume that the gradients computed for normalized amplitudes are also the gradients with respect to the raw pixel values, even though the overflow-encoding map between them is nonlinear and its chain rule is never derived.

What would settle it

Run the method on a trained amplitude-embedded model with overflow encoding, then independently perturb a single raw pixel by a small ε and measure the change in the model output; if (F(x+εe_i)−F(x))/ε does not match the attribution gradient the method assigns to that pixel (up to shot noise), the pixel-level attribution claim fails. Equivalently, check whether the attributions satisfy the integrated-gradients completeness identity, sum_i IG_i(x)=F(x)−F(x′), on the actual encoding used; a violation would falsify the 'provably faithful' claim.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • Amplitude-embedded QML classifiers can be explained with per-feature attribution scores using only circuit evaluations, not classical state-vector simulation.
  • The same Hadamard-test construction works for any state-preparation circuit V(x), since the gradient identity is independent of how |x⟩ is prepared.
  • Attribution cost scales linearly with the number of features and can be reduced by a factor of 2^m−1 using m ancilla qubits.
  • Because the gradient is exact up to shot noise, attributions remain stable at very low shot counts (the paper reports 10 shots), making the method usable on near-term hardware.
  • The parameter-shift rule covers angle-encoded models, so HattriQ provides a unified attribution scheme for both common encoding strategies.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the missing chain-rule link between normalized amplitudes and raw pixel values is supplied, the same inner-product identity would extend to any differentiable normalization, making HattriQ a general quantum attribution primitive rather than a dataset-specific recipe.
  • The integrated-gradients completeness property could be turned into an on-device sanity check: sum the attributions and compare to F(x)−F(x′); a mismatch would flag shot-noise bias or implementation errors without any classical simulation. The paper does not propose this test.
  • The multi-ancilla parallelization trades qubit count for circuit depth; on devices with limited connectivity, the controlled preparations may dominate the overhead, so a depth–qubit tradeoff analysis would be a natural next step.
  • The same Hadamard-test formalism might apply to parameter gradients, which would connect attribution to training dynamics and could yield a hardware-native version of parameter-shift rules for non-rotation gates.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. This paper proposes HattriQ, a method for computing integrated-gradients feature attributions for quantum machine learning models with amplitude embedding. The authors derive a closed-form expression for the gradient of the expectation value F(x)=⟨x|U†(θ)OU(θ)|x⟩ with respect to the real and imaginary parts of the amplitude-encoding coefficients (Lemma 3.1), and show that each gradient component can be estimated by a Hadamard-test circuit whose ancilla measurement probability encodes Re[⟨b_k|U†OU|x⟩] (Theorem 4.2). A multi-ancilla extension is proposed to compute several gradient components in parallel. The method is evaluated on binary classification tasks on Bars and Stripes, NIST, MNIST, and FashionMNIST, with attributions visualized for amplitude- and angle-encoded models and a shot-noise study.

Significance. The core idea—computing amplitude-coordinate gradients via a Hadamard test rather than via parameter-shift or state tomography—is appealing, and the algebra in Lemma 3.1 and Theorem 4.2 is correct for the unnormalized amplitude coordinates. If the gap described below is closed, the resulting tool would be a useful interpretability primitive for QML: it is hardware-compatible, avoids accessing internal quantum states, and the parallelization proposal is a reasonable resource-optimization idea. The manuscript also ships open-source code and gives detailed circuit constructions. However, the paper's advertised 'provably faithful' pixel-level attributions are not supported by the equations as written.

major comments (3)
  1. [§5.2, Eq. (2), Lemma 3.1] The paper's headline claim is 'provably faithful' attributions for input features (pixels). Lemma 3.1 gives the gradient with respect to amplitude coordinates c_k of a unit-norm state, not with respect to raw pixel values p_i. In §5.2 the pixels are mapped to amplitudes by the overflow encoding a_i = p_i/s and a_overflow = sqrt(1−Σ a_i^2); this map is nonlinear, and every pixel affects the overflow amplitude. The Jacobian ∂a/∂p is never written down; the sentence 'All of our discussion from before still applies upon simple modification using the chain rule' is not a derivation. Hence the quantities in Figs. 1–3 are not shown to be the integrated gradients of Eq. (2), and the completeness axiom need not hold. The same applies to the tanh output composition in §5.2.
  2. [Def. 2.2, §3.2] The IG path x′+α(x−x′) lies in raw feature space. For an amplitude-encoded model, the intermediate states are generally not normalized; the circuit model F in Eq. (1) is only defined for unit-norm states. The paper does not specify how to evaluate the model at non-normalized points (renormalization? overflow encoding applied at each α?), and Lemma 3.1's derivative is for the homogeneous extension F(c)=⟨c|O|c⟩, not the physical normalized expectation. The integral in Eq. (2) is therefore not formally defined for the amplitude-embedded model as described.
  3. [§5.4] The shot-noise study compares the sampled Hadamard-test estimates to exact evaluation of the same inner-product formula; this demonstrates sampling convergence, not attribution correctness. Because of the missing chain rule, the pixel-level attributions are never checked against the completeness axiom Σ_i IG_i = F(x)−F(x′) or any other ground truth. A direct completeness check on the actual encoded model would be a meaningful validation.
minor comments (6)
  1. [Abstract] The abstract mentions validation on 'TFIM quantum data', but no TFIM experiment appears anywhere in the paper.
  2. [Table 1] The class-pair lists include entries such as '(0,0)' that are not binary classification tasks, and the number of accuracy entries does not match the number of listed pairs for MNIST and NIST. Please correct.
  3. [Appendix A] The notation in the rewritten F(x) uses '⟨bi|x∗i⟩' in a confusing way; it should be x_i^* ⟨b_i|.
  4. [Appendix B] The expression for P(A=0) mixes vector norm and absolute value ('| 1/2 (|bk⟩+U†OU|x⟩)|^2'); clarify that this is the squared norm of the data-register component.
  5. [§4.3 / Appendix C] The control conditions for the parallel multi-ancilla circuit are not shown in the circuit diagram, and the general claim of 2^m−1 gradient components is only illustrated for m=2 without a proof for general m.
  6. [References] The reference for the COBYLA optimizer (page 10) is incomplete; author names are missing.

Circularity Check

0 steps flagged

No significant circularity: HattriQ's gradient and Hadamard-test derivations are self-contained; the overflow-encoding chain-rule gap is an omitted derivation, not a circular reduction.

full rationale

The central derivation is self-contained. Lemma 3.1 (Sec. 3.2, App. A) differentiates F(x)=⟨x|U†OU|x⟩ with the product rule and obtains ∂F/∂c_k = 2Re[⟨b_k|U†OU|x⟩]; this is a direct calculation from the definition of F, not an assumption of the result. Theorem 4.2 (Sec. 4.2, App. B) derives the ancilla probability P(A=0)=1/2(1+Re[⟨b_k|U†OU|x⟩]) from the Hadamard-test circuit; the circuit construction and proof are standard and do not invoke fitting or self-citation. The multi-ancilla extension is likewise derived by explicit calculation (App. C). Integrated gradients are then the integral of these gradients along the path of Eq. 2, following Sundararajan et al.; no parameter is fitted and then renamed as a prediction. Self-citations (DiBrita et al. 2024; Han et al. 2025; Luo et al. 2024; Cho et al. 2025) appear in the introduction/related work and are not load-bearing for Lemma 3.1 or Theorem 4.2; there is no imported uniqueness theorem and no ansatz smuggled through a self-citation. The real weakness is not circularity: Sec. 5.2 introduces the overflow amplitude encoding (and a tanh output) and asserts "All of our discussion from before still applies upon simple modification using the chain rule," but the Jacobian connecting raw pixel features to the overflow-normalized amplitudes is never written down. This is an omitted derivation in the path from amplitude gradients to the pixel attributions of Figs. 1–3; it undermines the empirical claim but is not an equivalence-by-construction of the kind counted as circular. For the mathematical object actually treated in the paper's lemmas (the homogeneous quadratic function of the amplitude vector), the derivation does not reduce to its inputs.

Axiom & Free-Parameter Ledger

4 free parameters · 4 axioms · 0 invented entities

The central equations are self-contained, but the mapping from raw pixel features to unit-norm amplitude states introduces a nonlinear normalization that is not accounted for in the attribution derivation. No quantities are fitted to make the formulas work; the listed free parameters are modeling choices that set the scale and meaning of the attributions.

free parameters (4)
  • IG baseline = all-zero (blank) image
    Integrated gradients requires a reference input; all attributions are relative to this baseline and no sensitivity study is reported (§5.3).
  • Pixel scaling factor for amplitude encoding = features scaled into [0, (2^n−1)^(-1/2)]
    The overflow encoding is chosen to make amplitude-embedded states unit-norm; attribution units depend on this arbitrary scaling (§5.2).
  • Number of measurement shots = 10, 100, 500, and exact simulation
    Gradient components are estimated by sampling the Hadamard-test ancilla; shot count sets the attribution noise (§5.4).
  • Ansatz depth = 1–2× qubit count per dataset
    Hardware-efficient ansatz layer counts were chosen empirically; Table 1 reports chosen depths without a sensitivity analysis (§5.2, Appendix D).
axioms (4)
  • domain assumption Model output has the form F(x;θ) = ⟨x|U†(θ)OU(θ)|x⟩, with U unitary and O Hermitian.
    Eq. 1. If a nonlinear tanh output is used, the paper says the chain rule applies but does not derive the modified gradient (§5.2).
  • domain assumption The amplitude-encoded input state must have unit norm; the overflow encoding x → (x_i, sqrt(1−Σ|x_i|²)) maps pixel features to normalized states.
    §5.2: 'we cannot directly encode them as quantum states... must first apply some transformation... maintain the normalization condition of quantum states.'
  • ad hoc to paper The integrated-gradients path x′ + α(x−x′) is a valid domain for the model, and the amplitude-gradient formula ∂F/∂x_k = 2Re⟨b_k|U†OU|x⟩ applies along it.
    Eq. 2 and Lemma 3.1. The path between normalized amplitude-encoded states is generally not normalized, and the overflow component is a nonlinear function of the pixel features; no chain rule through normalization is given.
  • domain assumption The observable O must be both Hermitian and unitary so that U†OU can be implemented as a controlled gate.
    Remark 3.2 and Theorem 4.2 require this. The paper uses the single-qubit Z operator, but the 'general-purpose' framing depends on this restriction.

reviewed 2026-08-04 · how reviews work

0 comments
Cite this review

Pith. "Pith review of HattriQ: Designing Integrated Gradients for Feature Attribution in Quantum Machine Learning." pith.science (2026). https://pith.science/paper/BQBIEPQQ

@misc{pith2026251002497,
  author       = {Pith},
  title        = {Pith review of: HattriQ: Designing Integrated Gradients for Feature Attribution in Quantum Machine Learning},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/BQBIEPQQ}},
  note         = {Machine review of arXiv:2510.02497}
}
Share X Bluesky LinkedIn Reddit HN
read the original abstract

Quantum machine learning (QML) algorithms have demonstrated early promise across hardware platforms, but remain difficult to interpret due to the inherent opacity of quantum state evolution. Widely used classical interpretability methods, such as integrated gradients and surrogate-based sensitivity analysis, are not directly compatible with quantum circuits due to measurement collapse and the exponential complexity of simulating state evolution. In this work, we introduce HattriQ, a general-purpose framework for computing amplitude-based input-attribution scores in circuit-based QML models. HattriQ supports the widely-used input amplitude embedding feature encoding scheme and uses a Hadamard test-based construction to compute input gradients directly on quantum hardware to compute integrated gradient attributions. We validate HattriQ on classification tasks across several datasets (Bars and Stripes, MNIST, FashionMNIST, and TFIM quantum data).

Figures

Figures reproduced from arXiv: 2510.02497 by Hengrui Luo, Jason Han, Nicholas S. DiBrita, Tirthak Patel, Younghyun Cho.

Figure 1
Figure 1. Figure 1: Sample images and the accompanying integrated gradients attribution for various samples [PITH_FULL_IMAGE:figures/full_fig_p007_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Sample images and gradient attribution for the amplitude-embedded MNIST dataset. We [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Sample images and attributions for the FashionMNIST dataset using amplitude encoding. [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Sample images and the accompanying integrated gradients attribution for the Bars and [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Integrated gradients computed using various amounts of measurement shots (samples). [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Sample images and the accompanying integrated gradients attribution for various samples [PITH_FULL_IMAGE:figures/full_fig_p017_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Sample images and the accompanying integrated gradients attribution for various samples [PITH_FULL_IMAGE:figures/full_fig_p018_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Sample images and the accompanying integrated gradients attribution for various samples [PITH_FULL_IMAGE:figures/full_fig_p018_8.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

19 extracted references · 1 canonical work pages

  1. [2]

    URL https://link.aps.org/doi/10.1103/RevModPhys.94.015004

    doi: 10.1103/RevModPhys.94.015004. URL https://link.aps.org/doi/10.1103/RevModPhys.94.015004. Joseph Bowles, Shahnawaz Ahmed, and Maria Schuld. Better than classical? the subtle art of benchmarking quantum machine learning models.arXiv preprint arXiv:2403.07059,

  2. [5]

    Scalable and interpretable quantum natural language processing: an implementation on trapped ions.arXiv preprint arXiv:2409.08777,

    Tiffany Duneau, Saskia Bruhn, Gabriel Matos, Tuomas Laakkonen, Katerina Saiti, Anna Pearson, Konstantinos Meichanetzidis, and Bob Coecke. Scalable and interpretable quantum natural language processing: an implementation on trapped ions.arXiv preprint arXiv:2409.08777,

  3. [6]

    Enqode: Fast amplitude embedding for quantum machine learning using classical data.arXiv preprint arXiv:2503.14473,

    11 Jason Han, Nicholas S DiBrita, Younghyun Cho, Hengrui Luo, and Tirthak Patel. Enqode: Fast amplitude embedding for quantum machine learning using classical data.arXiv preprint arXiv:2503.14473,

  4. [8]

    The disagreement problem in explainable machine learning: A practitioner’s perspective.arXiv preprint arXiv:2202.01602,

    Satyapriya Krishna, Tessa Han, Alex Gu, Steven Wu, Shahin Jabbari, and Himabindu Lakkaraju. The disagreement problem in explainable machine learning: A practitioner’s perspective.arXiv preprint arXiv:2202.01602,

  5. [9]

    Quantum embeddings for machine learning.arXiv preprint arXiv:2001.03622,

    Seth Lloyd, Maria Schuld, Aroosa Ijaz, Josh Izaac, and Nathan Killoran. Quantum embeddings for machine learning.arXiv preprint arXiv:2001.03622,

  6. [12]

    doi: 10.18653/v1/P18-1176

    Association for Computational Linguistics. doi: 10.18653/v1/P18-1176. URL https://aclanthology.org/P18-1176/. Kim Nicoli, Christopher J Anders, Lena Funcke, Tobias Hartung, Karl Jansen, Stefan Kühn, Klaus- Robert Müller, Paolo Stornati, Pan Kessel, and Shinichi Nakajima. Physics-informed bayesian optimization of variational quantum circuits.Advances in Ne...

  7. [13]

    Maria Schuld and Nathan Killoran

    URL https://proceedings.neurips.cc/ paper_files/paper/2024/file/386432c7534eec9a1cd7cbeea90d7e9f-Paper-Conference.pdf. Maria Schuld and Nathan Killoran. Quantum machine learning in feature hilbert spaces.Physical review letters, 122(4):040504,

  8. [14]

    Supervised learning with quantum computers.Quantum science and technology (Springer, 2018),

    Maria Schuld and Francesco Petruccione. Supervised learning with quantum computers.Quantum science and technology (Springer, 2018),

  9. [15]

    URL https://link.aps.org/doi/10.1103/PhysRevA.99.032331

    doi: 10.1103/ PhysRevA.99.032331. URL https://link.aps.org/doi/10.1103/PhysRevA.99.032331. 13 Maria Schuld, Alex Bocharov, Krysta M. Svore, and Nathan Wiebe. Circuit-centric quantum classifiers.Phys. Rev. A, 101:032308, Mar

  10. [17]

    Andrew G White, DFV James, William J Munro, and PG Kwiat

    doi: 10.1038/s41592-019-0686-2. Andrew G White, DFV James, William J Munro, and PG Kwiat. Exploring hilbert space: Accurate characterization of quantum information.Physical Review A, 65(1):012301,

  11. [18]

    Yusen Wu, Bujiao Wu, Yanqi Song, Xiao Yuan, and Jingbo Wang

    doi: 10.22331/q-2022-03-30-677. Yusen Wu, Bujiao Wu, Yanqi Song, Xiao Yuan, and Jingbo Wang. Learning the complexity of weakly noisy quantum states. InThe Thirteenth International Conference on Learning Representations,

  12. [2013]

    Pramod Kaushik Mudrakarta, Ankur Taly, Mukund Sundararajan, and Kedar Dhamdhere. Did the model understand the question? In Iryna Gurevych and Yusuke Miyao (eds.),Proceedings of the 56th Annual Meeting of the Association for Computational Linguistics (Volume 1: Long Papers), pp. 1896–1906, Melbourne, Australia, July

  13. [2018]

    URL https://link.aps.org/doi/10.1103/ PhysRevA.98.032309

    doi: 10.1103/PhysRevA.98.032309. URL https://link.aps.org/doi/10.1103/ PhysRevA.98.032309. 12 Ashley Montanaro and Ronald de Wolf. A survey of quantum property testing.arXiv preprint arXiv:1310.2035,

  14. [2019]

    What is my quantum computer good for? quantum capability learning with physics-aware neural networks.arXiv preprint arXiv:2406.05636,

    Daniel Hothem, Ashe Miller, and Timothy Proctor. What is my quantum computer good for? quantum capability learning with physics-aware neural networks.arXiv preprint arXiv:2406.05636,

  15. [2020]

    URL https://link.aps.org/doi/10.1103/PhysRevA.101.032308

    doi: 10.1103/PhysRevA.101.032308. URL https://link.aps.org/doi/10.1103/PhysRevA.101.032308. Daniel Silver, Tirthak Patel, William Cutler, Aditya Ranjan, Harshitta Gandhi, and Devesh Tiwari. Mosaiq: Quantum generative adversarial networks for image generation on nisq computers. In Proceedings of the IEEE/CVF International Conference on Computer Vision, pp....

  16. [2022]

    Kishor Bharti, Alba Cervera-Lierta, Thi Ha Kyaw, Tobias Haug, Sumner Alperin-Lea, Abhinav Anand, Matthias Degroote, Hermanni Heimonen, Jakob S

    URL https://arxiv.org/abs/1811.04968. Kishor Bharti, Alba Cervera-Lierta, Thi Ha Kyaw, Tobias Haug, Sumner Alperin-Lea, Abhinav Anand, Matthias Degroote, Hermanni Heimonen, Jakob S. Kottmann, Tim Menke, Wai-Keong Mok, Sukin Sim, Leong-Chuan Kwek, and Alán Aspuru-Guzik. Noisy intermediate-scale quantum algorithms.Rev. Mod. Phys., 94:015004, Feb

  17. [2023]

    Multi-qubit rydberg gates between distant atoms.arXiv preprint arXiv:2507.16602,

    Antonis Delakouras, Georgios Doultsinos, and David Petrosyan. Multi-qubit rydberg gates between distant atoms.arXiv preprint arXiv:2507.16602,

  18. [2024]

    Efficiently manipulating pauli strings with pauliarray.arXiv preprint arXiv:2405.19287,

    Maxime Dion, Tania Belabbas, and Nolan Bastien. Efficiently manipulating pauli strings with pauliarray.arXiv preprint arXiv:2405.19287,

  19. [2025]

    Fashion-mnist: a novel image dataset for benchmarking machine learning algorithms.arXiv preprint arXiv:1708.07747,

    Han Xiao, Kashif Rasul, and Roland V ollgraf. Fashion-mnist: a novel image dataset for benchmarking machine learning algorithms.arXiv preprint arXiv:1708.07747,

This paper was first reviewed by deepseek-v4-flash on August 4, 2026.