{"id":"db1ee297-6634-43ab-9206-5e99f0199b75","arxiv_id":"2412.16811","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A Magnus-expansion analysis of product formula errors for quantum phase estimation yields custom product formulas with up to quartic energy-error scaling and a low-energy bound with up to quadratic speedup in the target error.","lead":"This paper derives tighter error bounds for product formulas used in quantum phase estimation and designs new product formulas whose energy-measurement error scales as the fourth power of the step size. The result could make quantum energy calculations, such as in chemistry, require substantially fewer simulation steps.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Supplementary Eq. (22) does not enforce the claimed equal-coefficient condition for W2, so the printed construction fails to yield the asserted O(s^3) and O(s^4) energy errors.","rationale":"The reader's weakest assumption was the spectral-gap condition in the perturbation expansion. That is a standard nondegenerate-perturbation assumption and, while worth stating, is not the most load-bearing point for this paper. The more concrete and central issue is the derivation of the custom product formula coefficients: the printed Eq. (22) does not imply the claimed leading error form. I verified this by expanding the exponent operator from Eq. (14) for the five-term alternating formula, computing the coefficients of [A,[A,B]] and [B,[A,B]], and comparing with Eq. (22) after imposing the first-order and second-order conditions. The mismatch is numerical, not a stylistic concern: for a4=-0.3, the Eq. (22) solution leaves a nonzero non-H-outer component in the O(s^2) effective-Hamiltonian correction, which would make the QPE error O(s^2) for generic H=A+B. The numerical plots in Fig. 1 are encouraging and suggest the authors may have used a different, correct coefficient condition, but the manuscript provides no explicit coefficients or code, so the central existence claim is not reproducible from the printed equations. The recommended verdict remains CONDITIONAL, but the condition should explicitly require correcting Eq. (22) (or equivalently providing the correct coefficient equations and the values used in the numerics) and re-verifying the W2 and W2W2' scalings. I am not claiming the result is false; I am claiming the printed proof of the central construction is internally inconsistent and needs repair.","tokens_in":9597,"tokens_out":52658,"duration_ms":373222,"concrete_test":"Impose Eqs. (20)-(21) together with the actual equal-coefficient condition f_AA=f_BA computed from Eq. (14), and solve for a4=-0.3; the roots are a3 approximately 1.181 (or -0.181), not the Eq. (22) roots 0.021 and 6.908. Then rerun the Fig. 1 (right) experiment for W2(s/2)W2'(s/2) on a random 9-qubit Hamiltonian with both the printed Eq. (22) solution and the corrected f_AA=f_BA solution. If the printed solution gives deltaE scaling as s^2 while the corrected solution gives s^4, the error is confirmed. Alternatively, expand log W2(s) to O(s^3) for the printed solution and verify that the O(s^2) effective-Hamiltonian term contains [A,[A,B]]-[B,[A,B]], which is inconsistent with Eq. (16).","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central custom-formula result rests on the claim that the five-exponential W2 satisfies W2(s)-e^{-iHs}=alpha3 s^3[H,[H,B]]+O(s^4) (Eq. 16), equivalently F(sigma)=3alpha3 sigma^2[H,[H,B]]+O(sigma^3) in Eq. (19). Expanding F(sigma) from Eq. (14) in the basis [A,[A,B]] and [B,[A,B]] gives coefficients f_AA and f_BA; the stated form requires exactly f_AA=f_BA. The paper's Eq. (22) is not algebraically equivalent to f_AA=f_BA after imposing Eqs. (20)-(21). For the quoted example a4=-0.3, solving Eqs. (20)-(22) gives a3 approximately 0.021 (or 6.908), whereas f_AA-f_BA is approximately 0.137 (or -23.8), not zero. Consequently the effective-Hamiltonian correction contains a nonzero component [A,[A,B]]-[B,[A,B]] at O(s^2), whose expectation in a generic eigenstate of H=A+B need not vanish. Thus deltaE would be O(s^2), not the claimed O(s^3); the same O(s^2) contamination destroys the O(s^4) claim for W2(s/2)W2'(s/2). The numerics cannot be checked against the printed equations because explicit coefficient values are not given. This is an internal inconsistency in the derivation, not a disagreement with consensus.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a Magnus-based perturbative analysis of the error in energy estimation (e.g., QPE) when the Hamiltonian evolution is implemented with product formulas. The main theoretical result is Eq. (9), which expresses the leading energy error as an expectation value of the exponent operator of the product formula, and generalizes previous first-order results to arbitrary order. Based on this, the authors propose a family of five-exponential second-order product formulas W2 for H=A+B, and a nine-exponential symmetrized combination W2(s/2)W2'(s/2), claiming energy errors O(s^3) and O(s^4), respectively. They also connect their analysis to low-energy bounds from prior work [12] to argue for improved cost scaling for low-energy states of certain k-local Hamiltonians.","tokens_in":9911,"tokens_out":21398,"duration_ms":144287,"significance":"If the custom product formula result were correct, it would be a significant advance: a task-specific (rather than worst-case) error analysis for product formulas, and a concrete formula that outperforms second-order Trotter-Suzuki in energy estimation at equal exponentials count. The general framework in Eq. (9) and the low-energy bound are useful and likely correct. The paper states its spectral-gap assumption explicitly and provides a supplementary derivation of the perturbation bound. However, the central custom-formula construction contains an algebraic error (see major comments), and the numerics are not reproducible as reported. The general contributions do not compensate for the collapse of the paper's headline claim.","major_comments":[{"comment":"Supplementary Eq. (22) does not enforce the equal-coefficient condition needed for the central claim in Eq. (16). The condition that the sigma^2 term in F(sigma) (Eq. (19)) be proportional to [H,[H,B]] is f_AA,B = f_BA,B. Expanding Eq. (15) for the five exponentials in Eq. (17) gives f_AA,B = (1/2)a2 a3^2 + a2 a3 a5 + (1/2)a2 a5^2 + (1/2)a4 a5^2 - (a1 a2 a3 + a1 a2 a5 + a1 a4 a5 + a3 a4 a5) and f_BA,B = a2 a3 a4 - ((1/2)a1 a2^2 + a1 a2 a4 + (1/2)a1 a4^2 + (1/2)a3 a4^2). The printed Eq. (22) is not algebraically equivalent to f_AA,B = f_BA,B after imposing Eqs. (20)-(21). For the paper's stated example a4=-0.3, solving Eqs. (20)-(22) gives a3 approx 0.021 or a3 approx 6.908, and in both cases f_AA,B - f_BA,B is nonzero (approx 0.137 and approx -23.8, respectively). Consequently, the effective Hamiltonian contains a nonzero component [A,[A,B]]-[B,[A,B]] at O(s^2), whose expectation in a generic eigenstate of H=A+B need not vanish, so deltaE is O(s^2), not O(s^3). This invalidates Eq. (16) and, through Eq. (13), the O(s^4) claim for W2(s/2)W2'(s/2) and Table I.","section":"Supplementary material, Eq. (22)"},{"comment":"The numerical verification in Fig. 1 cannot be checked against the printed construction because the coefficient values a1, a2, a3, a5 used for a4=-0.3 are not reported. Given that the printed equations do not produce the claimed exponent operator, the observed scalings in Fig. 1 do not provide evidence for the validity of the construction as stated. The authors should provide the actual coefficient values and verify explicitly that they satisfy the correct equal-coefficient condition f_AA,B = f_BA,B, not Eq. (22).","section":"Section 'Customized product formulas'; Fig. 1"},{"comment":"The statement that 'one must choose the free variable outside of the range 0 < a4 < 1' is a property of the incorrect Eq. (22). The correct condition f_AA,B = f_BA,B leads to a different solution set; the existence of a real solution family and the admissible range of a4 must be re-established. Without a correct derivation, the claimed existence of the custom W2 family is unsupported.","section":"Main text after Eq. (10); Fig. 1 (left)"}],"minor_comments":[{"comment":"The figure legend labels the W2 curve as '~ s4', while the text says 'We expect an s3 behavior for W2, but the s4 behavior appearing here is likely due to the structure of the model.' This is confusing; please clarify whether the plotted scaling is the observed numerical slope or the leading-order theoretical prediction, and reconcile the two statements.","section":"Fig. 1 (middle) and text after Eq. (16)"},{"comment":"The phrase 'by inspection' is too terse for the claimed expansion W2(s/2)W2'(s/2) - e^{-iHs} = alpha3 s^3 [H,[H,B]] + alpha4 s^4 [H,[H,[H,B]]] + O(s^5); provide a short derivation or at least specify how the coefficients alpha3 and alpha4 are obtained from the W2 solution.","section":"Eq. (13)"},{"comment":"The symbol E(sigma) is used in Eq. (9) and in the text as the 'exponent operator', but it is not defined in the main text; please define it there or refer explicitly to the supplementary definition, and avoid any confusion with the eigenvalue E.","section":"Main text, Eq. (9) and surrounding text"},{"comment":"The nonvanishing spectral gap assumption is stated for all step sizes between 0 and s, but it is not checked for the numerical examples (XY model and random Hamiltonian). Please comment on whether this assumption is satisfied in the simulations, or how violation would affect the reported scalings.","section":"After Eq. (4)"}],"recommendation":"major_revision","confidential_remarks":"The algebraic error in Supplementary Eq. (22) is load-bearing and, as written, invalidates the paper's main advertised result. I have verified the counterexample numerically using the full expansion of Eq. (15): for a4=-0.3, the solutions of Eqs. (20)-(22) give f_AA,B - f_BA,B approx 0.137 (or approx -23.8), so the claimed O(s^3) energy error for W2 does not follow. This is not a circularity or consensus issue; it is an internal inconsistency. I recommend major revision rather than rejection because the general Magnus-based framework (Eq. (9)) and the low-energy analysis appear sound and are separable from the custom-formula construction. The authors should re-derive the coefficient conditions, re-check existence of solutions, and resubmit with reproducible coefficient values. If the corrected condition admits no real solution family, then the paper should be rejected."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe headline result—the 9-term formula with quartic energy error—does not survive contact with the supplementary algebra. The stress-test note is correct. Equation (22) is not algebraically equivalent to requiring f_AA = f_BA in Eq. (19), which is the necessary and sufficient condition for the leading error to be proportional to [H,[H,B]]. I checked with the a4 = -0.3 branch: solving (20)-(22) gives a3 ≈ 0.021, and then f_AA - f_BA ≈ 0.137, not zero. So the printed construction does not yield the claimed O(s^3) energy error, and the O(s^4) claim for W2(s/2)W2'(s/2) inherits the problem. The numerics cannot be verified because explicit coefficients are not given.\n\nThat said, the paper isn't empty. The Magnus-based analysis leading to Eq. (9) is a plausible generalization of Yi-Crosson to arbitrary order, and the argument that nested commutators with H on the outermost slot don't contribute to the QPE energy error is clean and useful. That part is worth taking seriously, and it may have independent value. The low-energy section is more speculative; it leans on the authors' unpublished [12] and contains a step-size expression that looks like a typo, but the idea of using low-energy bounds for energy estimation is a reasonable extension.\n\nThe problem is that the custom product formula is the paper's main new contribution, and the printed algebra doesn't support it. Either the coefficients were actually computed from a different (correct) condition, in which case the supplementary should be fixed and the coefficients provided, or the claim is simply wrong. As written, the paper is internally inconsistent.\n\nIf I were the editor, I would send it to a referee, but I'd flag the supplementary equation as the main thing to check. The Magnus-based QPE error analysis deserves a serious look, and if the custom formula can be repaired, the paper could be valuable. But in its current form, the central claim is not established.","headline":"The Magnus-based QPE error analysis is worth a look, but the custom product formula construction fails as printed: the supplementary condition does not enforce the claimed error scaling.","tokens_in":10452,"tokens_out":16681,"would_cite":false,"duration_ms":112973,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"For a gapped eigenstate, the energy-estimation error of an order-$p$ product formula equals the average of a nested-commutator expectation plus $O(s^{2p})$, and a custom 9-term formula attains fourth-order energy accuracy.","keywords":["quantum phase estimation","product formulas","Trotter-Suzuki formula","Magnus expansion","nested commutators","energy estimation error","low-energy bounds","Hamiltonian simulation"],"falsifier":"Compute the exact QPE energy error $\\delta E(s)$ for a two-level Hamiltonian $H=A+B$ whose effective Hamiltonian $\\tilde{H}(s)$ has a gap that closes at some step size inside $[0,s]$. If the measured error deviates from $\\frac{1}{s}\\int_0^s \\langle\\psi|E(\\sigma)|\\psi\\rangle\\,d\\sigma$ by a term that is not $O(s^{2p})$, the perturbation expansion fails. Alternatively, on a random multi-qubit Hamiltonian, if the 9-term formula $W_2(s/2)W_2'(s/2)$ shows energy-error scaling $s^2$ or $s^3$ rather than $s^4$ at small step sizes, the central design claim is wrong.","tokens_in":9414,"feed_emoji":"⚛️","tokens_out":11400,"duration_ms":80353,"temperature":0.7,"pith_summary":"The paper claims that when product formulas are used inside quantum phase estimation, the error in the estimated energy can be much smaller than the worst-case spectral error of the time-evolution operator. For an order-$p$ product formula on a gapped eigenstate, it derives a closed-form expansion: $\\delta E = \\frac{1}{s}\\int_0^s \\langle\\psi|E(\\sigma)|\\psi\\rangle\\,d\\sigma + O(s^{2p})$, where $E(\\sigma)$ is a sum of nested commutators of Hamiltonian terms. This formula turns the design of product formulas into a cancellation problem: any nested commutator with the full Hamiltonian on its outermost slot contributes nothing to the energy error. The paper exploits this to construct a 9-term second-order formula whose energy error is $O(s^4)$, matching a fourth-order Trotter-Suzuki formula that uses 11 terms.","feed_headline":"Nine exponentials match eleven for quantum energy estimation","feed_subtitle":"A tailored 9-term formula achieves 4th-order energy accuracy with fewer exponentials.","key_machinery":"The central object is the effective Hamiltonian $\\tilde{H}(s)$ defined by the single-step exponential identity $e^{-i\\tilde{H}(s)s}=W(s)$. The product formula is first written as a time-ordered exponential $W(s)=T\\exp(-i\\int_0^s F(\\sigma)\\,d\\sigma)$ with $F(\\sigma)=H+E(\\sigma)$, where $E(\\sigma)$ is an exponent operator built from nested commutators of the Hamiltonian terms, each weighted by product-formula coefficients. A Magnus expansion then converts the time-ordered exponential into the single exponential $e^{-i\\tilde{H}(s)s}$, giving $\\tilde{H}-H = \\frac{1}{s}\\int_0^s E(\\sigma)\\,d\\sigma$ plus higher nested-commutator terms. The identity $\\delta E=\\frac{1}{s}\\int\\langle\\psi|E(\\sigma)|\\psi\\rangle\\,d\\sigma+O(s^{2p})$ follows because terms with $H$ in the outermost commutator have zero expectation on an eigenstate. The custom product formula $W_2(s)=e^{-iAa_1s}e^{-iBa_2s}e^{-iAa_3s}e^{-iBa_4s}e^{-iAa_5s}$ is constructed so that the leading surviving term in $E(\\sigma)$ is proportional to $[H,[H,B]]$.","core_discovery":"The central discovery is a perturbation-theoretic identity for the error in energy estimates obtained from product formulas. For any product formula $W(s)$ of order $p$, define the effective Hamiltonian $\\tilde{H}(s)$ by $e^{-i\\tilde{H}(s)s}=W(s)$. Under a non-vanishing spectral gap, the shift in the target eigenvalue is, up to $O(s^{2p})$, the time-average of the expectation value of the exponent operator $E(\\sigma)$: $\\delta E = \\frac{1}{s}\\int_0^s \\langle\\psi|E(\\sigma)|\\psi\\rangle\\,d\\sigma + O(s^{2p})$ (Eq. (9)). Because $H$ appears as a factor in many of the nested commutators inside $E(\\sigma)$, and $\\langle\\psi|[H,\\cdot]|\\psi\\rangle=0$ for an eigenstate, many leading contributions vanish. For $H=A+B$, the paper exhibits a five-exponential second-order formula $W_2(s)$ with $E(\\sigma)=3\\alpha_3\\sigma^2[H,[H,B]]+O(\\sigma^3)$, making the energy error $O(s^3)$, and its symmetrized combination $W_2(s/2)W_2'(s/2)$ yields an error $O(s^4)$. For low-lying eigenstates of a class of $k$-local positive Hamiltonians, the coefficients in the expansion grow only polylogarithmically in system size up to order $2p-1$, which leads to a step-size scaling $s=O(\\epsilon^{1/p}/\\log(N)^{1+1/p}+\\epsilon^{1/(2p)}\\operatorname{poly}(N))$ and an up-to-quadratic asymptotic speedup in the target error.","pith_inferences":["The same cancellation criterion is generic: any order-$p$ formula whose nested commutators of orders $p$ through $2p-1$ all have $H$ on the outermost slot should inherit an $O(s^{2p})$ energy error; constructing such formulas at higher orders is left for future work.","Because the analysis depends only on the input/output structure of the energy-estimation routine, the step-size savings should carry over to other energy-estimation schemes that query simulated time evolution, such as amplitude-estimation-based approaches.","The observed $s^4$ scaling for $W_2$ on the XY model, instead of the generic $O(s^3)$, suggests that model-specific commutator structure can further improve the energy error beyond the universal family; this could be tested systematically on other lattice models.","If the paper's conjecture that the positivity assumption in the low-energy bound can be relaxed is correct, the polylogarithmic improvement would extend to fermionic Hamiltonians used in electronic structure simulation."],"forward_implications":["For a gapped eigenstate, an order-$p$ product formula with exponent operator $E(\\sigma)$ has energy-estimation error $\\delta E = \\frac{1}{s}\\int_0^s \\langle\\psi|E(\\sigma)|\\psi\\rangle\\,d\\sigma + O(s^{2p})$, so any nested-commutator term with $H$ on the outermost slot contributes nothing at leading order.","For $H=A+B$, the five-exponential formula $W_2(s)$ achieves $O(s^3)$ energy error, improving the required step size from $s=O(\\epsilon^{1/2})$ for second-order Trotter-Suzuki to $s=O(\\epsilon^{1/3})$.","The symmetrized 9-term formula $W_2(s/2)W_2'(s/2)$ achieves $O(s^4)$ energy error, matching a fourth-order Trotter-Suzuki formula while using 9 exponentials per step instead of 11.","For low-energy eigenstates of $k$-local positive-semidefinite fast-forwardable Hamiltonians, the step size can scale as $s=O(\\epsilon^{1/p}/\\log(N)^{1+1/p}+\\epsilon^{1/(2p)}\\operatorname{poly}(N))$, giving up to a quadratic asymptotic speedup in the target error $\\epsilon$.","The results apply to any energy-estimation routine that consumes simulated time evolution, not only to quantum phase estimation's internal protocol."],"supporting_citations":[{"why":"Supplies the nested-commutator, time-ordered-exponential expansion of product-formula error that defines the exponent operator E(σ) and F(σ)=H+E(σ).","marker":"[8]"},{"why":"Establishes the earlier first-order QPE error analysis and the s=O(ε^{1/2}) baseline that this paper generalizes to arbitrary order.","marker":"[10]"},{"why":"Provides the low-energy dynamics bounds for product formulas on k-local positive Hamiltonians that the paper adapts to energy-estimation error.","marker":"[12]"},{"why":"Gives the Magnus expansion used to convert the time-ordered exponential W(s) into the single effective exponential e^{-iH̃(s)s}.","marker":"[13]"},{"why":"Demonstrates the Magnus-based effective-Hamiltonian approach for low-order product formulas that is extended here to a perturbative expansion of arbitrary order.","marker":"[14]"},{"why":"Provides the explicit form of the Magnus expansion terms used to write each H̃_n as integrals over nested commutators.","marker":"[15]"},{"why":"Supplies the spectral-gap lemma bounding the projection difference, which controls the higher-order terms in the perturbation theory.","marker":"[17]"}],"fun_headline_variants":["Nine exponential terms yield quadratic energy error improvement","Quadratic speedup in phase estimation with 9-term formula","Better than Trotter-Suzuki: 9-term formula for quantum energies","Energy estimation errors drop quadratically with 9 exponentials"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analysis assumes the target energy level stays separated from all other levels of the simulated effective Hamiltonian at every step size up to the chosen one; if that gap closes, the predicted error formula no longer holds.","fun_headline_variants_meta":{"raw":{"variants":["Nine exponential terms yield quadratic energy error improvement","Quadratic speedup in phase estimation with 9-term formula","Better than Trotter-Suzuki: 9-term formula for quantum energies","Energy estimation errors drop quadratically with 9 exponentials"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000449,"raw_usage":{"total_tokens":2338,"prompt_tokens":1092,"completion_tokens":1246,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":708,"completion_tokens_details":{"reasoning_tokens":1176}},"tokens_in":708,"tokens_out":1246,"duration_ms":10977,"temperature":1.0,"reasoning_tokens":1176,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T10:16:10.863880+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the exact QPE energy error $\\delta E(s)$ for a two-level Hamiltonian $H=A+B$ whose effective Hamiltonian $\\tilde{H}(s)$ has a gap that closes at some step size inside $[0,s]$. If the measured error deviates from $\\frac{1}{s}\\int_0^s \\langle\\psi|E(\\sigma)|\\psi\\rangle\\,d\\sigma$ by a term that is not $O(s^{2p})$, the perturbation expansion fails. Alternatively, on a random multi-qubit Hamiltonian, if the 9-term formula $W_2(s/2)W_2'(s/2)$ shows energy-error scaling $s^2$ or $s^3$ rather than $s^4$ at small step sizes, the central design claim is wrong.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the nested-commutator, time-ordered-exponential expansion of product-formula error that defines the exponent operator E(σ) and F(σ)=H+E(σ)."},{"cited_title":"Yi and E","cited_arxiv_id":null,"evidence_quote":"Demonstrates the Magnus-based effective-Hamiltonian approach for low-order product formulas that is extended here to a perturbative expansion of arbitrary order."},{"cited_title":"S ¸ahino˘ glu and R","cited_arxiv_id":null,"evidence_quote":"Provides the explicit form of the Magnus expansion terms used to write each H̃_n as integrals over nested commutators."}],"review_version":1}