{"id":"ac7d5159-97a1-4c78-8c04-24cf7ce487e4","arxiv_id":"2607.10772","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":3,"one_line_summary":"QPE-based quantum amplitude estimation recovers fission-neutron-yield and U-238 resonance-integral expectations with O(1/T²) squared-error scaling, confirming quadratic Monte Carlo speedup on nuclear-engineering examples.","lead":"Quantum amplitude estimation is shown to estimate fission-neutron yield and a U-238 resonance integral with squared error falling as 1/T², a quadratic improvement over classical Monte Carlo. The result matters because reactor-physics Monte Carlo is limited by slow statistical convergence and quantum methods could cut the cost of high-precision calculations.","discovery_kind":"new_application","skeptic_critique":{"model":"grok-4.5","headline":"No significant objection identified beyond the already-flagged surrogate limitation.","rationale":"The paper’s strongest claim is a transparent numerical confirmation of the known QAE scaling on two nuclear-engineering integrals. Both the gate-level toy circuit and the eigendecomposition surrogate produce squared-error slopes of −2 that match the analytical mode-error formula, so the ideal-case quadratic advantage is solidly demonstrated. The sole material caveat—the surrogate does not yield a hardware-ready circuit—is already stated by the authors and correctly elevated by the reader as the conditioning factor. No additional load-bearing flaw (incorrect encoding of the Breit-Wigner amplitudes, mis-decoding of the QPE register, or mismatch between claimed and observed error) appears in the text. Consequently the reader’s CONDITIONAL verdict, low correctness risk, and identification of the weakest assumption remain appropriate; no adjustment is warranted.","tokens_in":7156,"tokens_out":455,"duration_ms":6627,"concrete_test":"Independently recompute the m=14 mode of the QPE register from the analytical eigenvalues of Q (or from a small-n gate-level circuit for the toy problem) and verify that |μ̂−μ|^{2} still tracks the O(1/T^{2}) bound of Eq. (3) to within shot noise; any systematic deviation would indicate an encoding or decoding error.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is the ideal-case O(1/T^{2}) squared-error scaling of QPE-based QAE on two nuclear expectation values, with the U-238 resonance integral recovered to ~0.03% relative error at m=14. The toy problem is fully gate-level and matches the analytical mode-error formula; the resonance-integral case uses an exact eigendecomposition of the 512\times512 Grover operator that is mathematically identical to the circuit under ideal simulation (explicitly stated in §II and Conclusions). The reader already correctly identifies the hardware-readiness gap (state-preparation cost / noise). No further internal inconsistency, hidden assumption, or numerical discrepancy undermines the ideal-case scaling that is actually asserted.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The manuscript applies standard QPE-based quantum amplitude estimation (QAE) to two nuclear-engineering expectation values: a four-outcome mean fission-neutron yield and the U-238 6.674 eV resonance integral under a 1/E spectrum. The toy problem is realized as a gate-level Qiskit circuit; the resonance-integral case is obtained by exact eigendecomposition of the 512×512 Grover operator to avoid state-preparation decomposition cost. In both cases the squared error is shown to follow the theoretical O(1/T²) scaling with oracle calls T (Eq. 3), in contrast to classical Monte Carlo O(1/N). At m=14 the resonance integral is recovered to ~0.03% relative error, corresponding to a query-complexity speedup ratio N/T of order 10³ relative to classical sampling at the same MSE.","tokens_in":7307,"tokens_out":1307,"duration_ms":31997,"significance":"The work is a clean numerical confirmation of the known Brassard–Montanaro QAE convergence theorem on two physically motivated nuclear integrands. Strengths include (i) an explicit gate-level Qiskit construction for the discrete yield problem that matches the analytical mode-error formula, (ii) side-by-side classical MC baselines on the same log-log plots, and (iii) transparent disclosure that the resonance-integral demonstration uses an eigendecomposition surrogate. The paper does not claim a new algorithm or a hardware-ready circuit for the larger problem; its value is as a domain-specific demonstration that the ideal-case quadratic query advantage is attainable for reactor-physics expectation values. That is a useful, if incremental, contribution for the nuclear-engineering / quantum-computing interface.","major_comments":[{"comment":"§II (U-238) and Conclusions: the resonance-integral result that underpins the abstract’s 0.03% / N/T≈3016 claim is obtained from an exact eigendecomposition of Q, not from a compiled circuit. While the paper states that the surrogate is “mathematically identical” under ideal simulation, the reported speedup ratio compares classical sample count N only to the number of oracle applications T=2^m−1. It does not account for the O(2^{n_x}) multi-controlled Ry cost of the state-preparation operator A that is bypassed by the surrogate. Because the central claim is a “demonstration” of quadratic Monte Carlo speedup for a nuclear quantity, the manuscript should either (a) supply a gate-complexity accounting that shows the query advantage survives realistic A cost, or (b) rephrase the abstract, results paragraph, and speedup ratio so that they are explicitly limited to the ideal query model and do","section":"§II U-238 Resonance Integral / Conclusions"},{"comment":"Eq. (3) and Figs. 2–3: the plotted “QAE analytical mode error” and the circuit/eigendecomposition points coincide almost perfectly, which is expected only when |ψ⟩=A|0⟩ lies exactly in the two-dimensional eigenspace of Q and the mode of ŷ is used. The paper should state explicitly that the reported MSE is the squared mode-estimation error (infinite-shot limit of the most-probable register outcome), not the mean-squared error over the full QPE distribution. Without that clarification a reader may over-interpret the empirical points as shot-noise-limited performance rather than the ideal mode error already guaranteed by the theory.","section":"Eq. (3), Figs. 2–3"}],"minor_comments":[{"comment":"Abstract and Introduction: “quadratic speedup for Monte Carlo-type expectation values” is standard terminology, but a one-sentence reminder that the speedup is in the number of oracle calls (query complexity), not wall-clock time or total gate count, would prevent over-reading by a nuclear-engineering audience.","section":"Abstract / Introduction"},{"comment":"Fig. 1 caption: the circuit is shown for m=4; it would help to note that the controlled-Q^{2^k} blocks for m>4 are obtained by matrix powering rather than circuit repetition, as stated in the text.","section":"Fig. 1"},{"comment":"Eq. (2): the absolute-value form ˆθ=π|ŷ/2^m−1/2| is correct for the −e^{±2iθ} eigenvalues, but a brief remark that both peaks φ± map to the same µ would aid reproducibility.","section":"Eq. (2)"},{"comment":"Reference [6] is cited as “arXiv preprint (2024)” without an identifier; please supply the arXiv number or DOI.","section":"References"},{"comment":"Notation: the same symbol T is used for oracle calls and (implicitly) for temperature nowhere, but the classical sample count is N; a short glossary or consistent “T_QAE vs N_MC” labeling on the dual axes of Figs. 2–3 would improve readability.","section":"Figs. 2–3"},{"comment":"C_phys≈1578 b is given as a classically precomputed constant; stating the precise numerical value used (or the quadrature rule) would make the IR conversion fully reproducible.","section":"§II"}],"recommendation":"minor_revision","confidential_remarks":"The manuscript is a solid, honest demonstration paper rather than an algorithmic advance. The eigendecomposition surrogate is the only material caveat and is already flagged by the authors; once the abstract and speedup-ratio language are tightened to the query model, the paper is suitable for a methods / applications venue. Fit for a pure quant-ph theory journal is weaker; fit for a nuclear-engineering or quantum-for-science journal is good."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"This is a clean methods note that takes the known QAE quadratic speedup and shows it works, in ideal simulation, on two nuclear-engineering expectations: a four-outcome fission-yield mean and a U-238 resonance integral under a 1/E spectrum. That domain application is the actual novelty; the algorithm and the error bound are textbook Brassard/Montanaro.\n\nWhat they do well is straightforward. The toy problem is a full gate-level Qiskit circuit with an explicit Ry-tree state prep; the squared-error plots sit on the analytical mode-error curve and give the clean slope -2 versus classical -1. For the resonance integral they discretize into 256 bins, build the 512-dimensional Grover operator analytically, and run QPE via eigendecomposition + FFT. At m=14 they recover IR to ~0.03 % relative error and show that matching that MSE classically would take ~3 000\times more samples. The math is standard, the classical baselines are side-by-side, and they never hide the surrogate.\n\nThe soft spot is exactly the one they name: multi-controlled state preparation for n≥9 is already prohibitive on the simulator, so the more interesting example never produces a hardware-ready circuit. Noise and realistic prep cost could erase the practical edge at this scale; they point to future arithmetic circuits and iterative QAE. No code is shipped. Those are real limitations on impact, not on the ideal-case claim that is actually asserted.\n\nWho it is for: people already working quantum-for-nuclear or looking for concrete QAE encodings outside finance. It is not a theoretical advance and does not claim to be. The citation pattern is appropriate (Brassard, Montanaro, finance demos, the one prior nuclear preprint). I would send it to referees; it is honest enough and useful enough for that niche that a serious editor should not desk-reject it. Engage if the nuclear-quantum intersection is on your radar; otherwise file it as a clean confirmation note.","headline":"Solid, transparent numerical demo that standard QPE-QAE recovers the expected O(1/T^{2}) scaling on two nuclear expectation values; the resonance case uses an ideal eigendecomposition surrogate that the authors flag clearly.","tokens_in":7938,"tokens_out":517,"would_cite":false,"duration_ms":6907,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Quantum amplitude estimation delivers quadratic Monte Carlo speedup on nuclear engineering expectations, recovering a U-238 resonance integral to 0.03% error with far fewer oracle calls than classical sampling.","keywords":["quantum amplitude estimation","Monte Carlo speedup","quantum phase estimation","resonance integral","fission neutron yield","nuclear engineering","U-238","Grover iterate"],"falsifier":"Compile an explicit gate-level circuit for the same 256-bin Breit-Wigner resonance integral (or a smaller but still multi-controlled version) and measure whether the observed squared-error scaling remains O(1/T^{2}) once state-preparation cost and realistic noise are included.","tokens_in":8014,"feed_emoji":"⚛","tokens_out":972,"duration_ms":12066,"temperature":0.7,"pith_summary":"Classical Monte Carlo is the gold standard for high-fidelity reactor calculations, but its error falls only as one over the square root of the number of particle histories. This paper shows that quantum amplitude estimation can improve that scaling to one over the number of oracle calls, a quadratic advantage. Using quantum phase estimation, the authors encode two nuclear quantities as amplitudes: a simple mean fission-neutron yield and the U-238 resonance integral under a 1/E slowing-down spectrum. In both cases the squared error falls as 1/T squared, matching theory. For the resonance integral they recover 120.65 barns versus the true 120.69 barns (0.03 percent relative error) with 14 phase-estimation qubits, a precision that would classically require tens of millions of samples. The work therefore supplies a concrete demonstration that quantum speedup for Monte Carlo-type nuclear integrals is already visible on present-day simulators and points toward hardware-ready circuits for larger resonance problems.","feed_headline":"Quantum amplitude estimation cuts Monte Carlo cost for nuclear integrals","feed_subtitle":"U-238 resonance integral recovered to 0.03% error with three orders fewer samples than classical MC","key_machinery":"Quantum amplitude estimation via quantum phase estimation of the Grover iterate Q: the target expectation is encoded as the success probability μ = sin^{2}\theta of an ancilla, Q rotates the good/bad subspace by 2\theta, and QPE extracts \theta (hence μ) from the eigenphase after T = 2^m - 1 applications of Q.","core_discovery":"QPE-based quantum amplitude estimation applied to two nuclear-engineering expectation values produces squared error that scales as O(1/T^{2}) with oracle calls T, confirming the theoretical quadratic improvement over classical Monte Carlo's O(1/N). For the U-238 resonance integral the method reaches approximately 0.03 percent relative error at m=14 phase-estimation qubits.","pith_inferences":["If the state-preparation bottleneck can be removed by quantum arithmetic rather than multi-controlled rotations, the demonstrated three-order-of-magnitude sample-count advantage becomes a realistic target for early fault-tolerant machines.","The same QAE template can be reused for other continuum integrals in nuclear data evaluation (Doppler-broadened resonances, unresolved-resonance averages) without redesigning the amplification step.","Noise-aware iterative QAE would let current intermediate-scale devices test the quadratic scaling on the simpler fission-yield toy problem before full resonance circuits become available."],"forward_implications":["Nuclear Monte Carlo integrals that today require tens of millions of particle histories can in principle be estimated to the same precision with a few tens of thousands of coherent oracle calls.","Once fault-tolerant arithmetic circuits for Breit-Wigner amplitudes exist, the same QAE pipeline can be extended to n ≥ 18 energy bins and multiple resonances.","Iterative amplitude-estimation variants can be substituted for full QPE to reduce circuit depth on near-term noisy hardware while retaining the quadratic scaling.","The same encoding of a 1/E spectrum and energy-dependent cross section applies immediately to other resonance integrals and reaction-rate tallies used in reactor physics."],"fun_headline_variants":["QAE yields O(1/T²) Monte Carlo error for nuclear integrals","Quantum amplitude estimation hits 0.03% U-238 integral error","Quadratic speedup via QAE on fission yield and U-238 resonance","QPE-based QAE confirms quadratic Monte Carlo gain in nuclear cases","Nuclear expectation values scale as 1/T² with quantum amplitude estimation"],"cache_read_input_tokens":128,"weakest_assumption_plain":"The resonance-integral result rests on an exact matrix eigendecomposition of the Grover operator rather than a fully compiled gate-level circuit, because state-preparation decomposition already becomes prohibitive at nine system qubits.","fun_headline_variants_meta":{"raw":{"variants":["QAE yields O(1/T²) Monte Carlo error for nuclear integrals","Quantum amplitude estimation hits 0.03% U-238 integral error","Quadratic speedup via QAE on fission yield and U-238 resonance","QPE-based QAE confirms quadratic Monte Carlo gain in nuclear cases","Nuclear expectation values scale as 1/T² with quantum amplitude estimation"]},"model":"grok-4.5","effort":"low","cost_usd":0.002676,"raw_usage":{"total_tokens":984,"prompt_tokens":711,"num_sources_used":0,"completion_tokens":100,"cost_in_usd_ticks":26760000,"prompt_tokens_details":{"text_tokens":711,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":173,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":711,"tokens_out":100,"duration_ms":3503,"temperature":1.0,"reasoning_tokens":173,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-14T09:21:38.184806+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Compile an explicit gate-level circuit for the same 256-bin Breit-Wigner resonance integral (or a smaller but still multi-controlled version) and measure whether the observed squared-error scaling remains O(1/T^{2}) once state-preparation cost and realistic noise are included.","supporting_citations":[],"review_version":1}