Pith. sign in

REVIEW 4 major objections 4 minor 1 cited by

Fast simulations of X-ray absorption spectroscopy for battery materials on a quantum computer

T0 review · 4 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read An optimized time-domain algorithm simulates X-ray absorption spectra of a battery cathode cluster with 100 logical qubits and $3.11 \times 10^8$ Toffoli gates per circuit.

desk verdict The Trotter step-selection idea is a real advance, but the headline resource numbers are conditional on an unvalidated Y3 bound and the abstract overstates what was verified. read the letter →

arxiv 2506.15784 v2 pith:KOP4SX7R submitted 2025-06-18 quant-ph

Stepan Fomichev , Pablo A. M. Casares , Jay Soni , Utkarsh Azad , Alexander Kunitsa , Arne-Christian Voigt , Jonathan E. Mueller , Juan Miguel Arrazola This is my paper · ORCID
classification quant-ph MSC 81P6881V55
keywords quantumchemistryX-rayabsorptionspectroscopyTrotterproductformulascompresseddoublefactorizationfault-tolerantcomputingresourceestimationtime-domainsimulationbatterycathodematerials
verification ladder T0 review T1 audit T2 compute T3 formal

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 aims to show that X-ray absorption spectroscopy (XAS) of industrially relevant battery cathode materials can be simulated on early fault-tolerant quantum computers at a practical cost. It optimizes the time-domain XAS algorithm for a Li4Mn2O cluster in a CAS(22e,18o) active space, reporting resource estimates of 100 logical qubits and $3.11 \times 10^8$ Toffoli gates for the longest circuit, a reduction of about five orders of magnitude compared with the unoptimized implementation. The central reason the cost is low is that spectroscopy only requires the eigenvalues of the simulated Hamiltonian to be accurate, not the full time-evolution operator, which allows the Trotter step size to be fixed independently of total evolution time. The paper backs the resource estimates with simulator reconstructions of the absorption spectrum against classical references.

What carries the argument

The load-bearing objects are the leading Trotter-error operator $Y_3$, obtained from the Baker–Campbell–Hausdorff expansion of the second-order product formula, and the compressed double factorization (CDF) of the electronic Hamiltonian. $Y_3$ controls the eigenvalue error through its diagonal matrix elements: the paper fixes the Trotter step by $\Delta^2|\langle E_l|Y_3|E_l\rangle| \leq 0.05$ Ha, and the perturbative estimate of this quantity is what decouples the step size from the total simulation time. CDF writes the two-electron integrals as a sum of $L$ fragments, each diagonal in its own single-particle basis, reducing the number of fragments from $O(N^4)$ in a direct Jordan-Wigner picture to $O(N^3)$ with $L = N$; the basis rotations between fragments are implemented with Givens rotations. The third machinery element is the Lorentzian-kernel sampling distribution, which allocates shots to evolution times proportionally to $e^{-\eta\tau|j|}$, replacing uniform sampling and saving a factor of about 18 in total shots.

What would settle it

Compute the diagonal matrix elements $\langle E_l|Y_3|E_l\rangle$ for the Li4Mn2O CAS(22e,18o) Hamiltonian directly, for example with a matrix-product-state approximation at increasing bond dimension, and check whether the largest value stays at or below 1 Ha; if it exceeds 1 Ha, the quoted Trotter step $\Delta = \sqrt{\eta/1\ \mathrm{Ha}}$ is too large and the reported Toffoli counts are underestimates.

Watch

Extended reading notes

Core claim

On the paper's own terms, the discovery is that a time-domain XAS calculation can be reorganized so that its cost is governed by spectral accuracy rather than by generic Hamiltonian simulation error. Using perturbation theory on the effective Hamiltonian produced by a second-order Trotter product formula, the eigenvalue shift is $E'_l = E_l - \Delta^2 \langle E_l|Y_3|E_l\rangle$, so choosing the step $\Delta$ with $\Delta^2|\langle Y_3\rangle| \leq 0.05$ Ha keeps peak positions within the 1 eV target broadening. Because this choice depends only on the Hamiltonian and not on how long the evolution is run, each Hadamard-test circuit costs linearly in time instead of superlinearly. Combined with the compressed double factorization of the Hamiltonian, a Lorentzian-weighted sampling distribution, and circuit-level optimizations, this brings the Li4Mn2O CAS(22e,18o) system to 100 logical qubits and $3.11 \times 10^8$ Toffoli gates per circuit ($4.58 \times 10^{11}$ including sampling), which the authors estimate runs in under a day on a MHz-clock fault-tolerant machine.

Load-bearing premise

The resource counts rest on the extrapolated bound $|\langle E_l|Y_3|E_l\rangle| \leq 1$ Ha for the CAS(22e,18o) Li4Mn2O cluster, obtained from small-system data in Fig. 13; the authors explicitly call for future work to validate this bound more robustly, and if the true value is larger the Trotter step must shrink and the gate counts rise.

Editorial extensions

If this is right

  • For the Li4Mn2O CAS(22e,18o) target, the longest circuit needs only 100 logical qubits and $3.11 \times 10^8$ Toffoli gates, and the full spectrum with sampling costs $4.58 \times 10^{11}$ Toffoli gates.
  • Because the Trotter step no longer needs to shrink as evolution time grows, the cost of each XAS circuit scales linearly with time, which is what makes long-time, fine-resolution spectra affordable.
  • The simulator benchmarks reproduce the valence absorption spectrum of N2 against a classical full-configuration-interaction reference and show convergence of the core-excited LiMnO spectrum with maximal evolution time.
  • Combined with an active-volume compilation on 350 logical qubits, the algorithm runs in about $5.6 \times 10^7$ logical cycles for the longest circuit, under a day at a 1 MHz logical clock rate, versus an estimated month-scale runtime for a restricted-active-space classical approach on the same target.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • As an extension beyond the paper, the eigenvalue-error criterion for Trotter steps should apply to any spectroscopy computed by time evolution under a product formula, not just XAS, provided the observable of interest is dominated by eigenenergy differences.
  • The $|\langle Y_3\rangle| \leq 1$ Ha bound could be tested classically on intermediate-size clusters between the small systems in Fig. 13 and the full CAS(22e,18o) target; a confirmation would turn the resource estimate from an extrapolation into a computed input.
  • The kernel-weighted sampling idea is generic: any damped spectral kernel (Lorentzian, Gaussian, or lifetime-broadened line shapes with known decay) could be used to reallocate shots, potentially reducing costs in other time-domain spectroscopy algorithms.
  • If the method scales as $O(N^3)$ as argued, the crossover against classical state-by-state methods should only widen with active-space size, making K-edge XAS of larger transition-metal oxide clusters a natural next target.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 4 minor

Summary. The manuscript presents a heavily optimized quantum algorithm for time-domain X-ray absorption spectroscopy (XAS), building on earlier work by Fomichev et al. The main algorithmic contributions are: (i) use of compressed double factorization (CDF) for product-formula time evolution; (ii) a perturbative analysis showing that, for spectroscopy, it is sufficient to control the eigenvalue error of the approximate Hamiltonian, allowing the Trotter step size to be fixed independently of total simulation time; and (iii) a Lorentzian-kernel-aware sampling distribution that drastically reduces the shot budget. Additional circuit-level optimizations include the double-phase trick, combining consecutive basis rotations, reduced rotation precision, BLISS symmetry shifts, and a double-measurement scheme. The paper reports constant-factor resource estimates for the Li4Mn2O cluster with CAS(22e,18o): 100 logical qubits and 3.11e8 Toffoli gates for the longest circuit, with 4.58e11 Toffoli gates for the full spectrum, a claimed five-orders-of-magnitude reduction compared to the unoptimized algorithm. The algorithm is benchmarked on a simulator: it reproduces the N2 valence spectrum against full configuration interaction, and shows convergence for a smaller LiMnO CAS(14e,11o) core-excited spectrum under the core-valence separation approximation.

Significance. If the resource estimates are taken at face value, this is a significant step toward making XAS simulation of industrially relevant battery cathode models practical on early fault-tolerant quantum computers. The paper is commendable for its transparency: the cost formulas in Sec. IV are explicit, all free parameters are enumerated, the error sources are discussed one by one, and the simulator benchmarks do validate the core time-domain algorithm against an independent classical reference for N2. The CDF-based Trotter implementation and the sampling-distribution analysis are careful and appear internally consistent. The main weakness is that the headline resource numbers for CAS(22e,18o) rest on the assumed bound |<Y3>| <= 1 Ha, which is extrapolated from small systems and is explicitly flagged by the authors as needing future validation. A second concern is that the 'eigenvalue-only' error criterion neglects eigenvector corrections that can distort peak intensities. Both issues are load-bearing for the central claims, so the paper needs revision before the resource estimates can be considered reliable.

major comments (4)
  1. [Sec. IV.B.2, Eq. (42)] The displayed inequality Δ ≤ (|⟨Y_{2k+1}⟩| / 1 eV)^{1/(2k)} is inverted with respect to Eq. (41). Given ϵ_trot = Δ^{2k} ⟨Y_{2k+1}⟩ and the requirement ϵ_trot ≤ 1 eV, the correct bound is Δ ≤ (1 eV / |⟨Y_{2k+1}⟩|)^{1/(2k)}. As written, Eq. (42) would force Δ to shrink as the error matrix element decreases and is dimensionally inconsistent. The actual choice in Sec. V.C.6, Δ^2 |⟨Y_3⟩| ≤ 0.05 Ha, follows the corrected form, so this is likely a typographical error; nevertheless Eq. (42) is the formal basis for the Trotter step and must be corrected.
  2. [Sec. V.C.6, Fig. 13, Table I] The resource estimates in Table I rely on the input |⟨E_l|Y_3|E_l⟩| ≤ 1 Ha for the CAS(22e,18o) Li4Mn2O cluster. This bound is not evaluated on the target system: Fig. 13 only covers active spaces up to N = 14, and the data are for the first five eigenvalues of small systems, not specifically for the core-excited states that dominate the XAS response. Since the Trotter step is set by Δ^2 |⟨Y_3⟩| ≤ 0.05 Ha and the number of Trotter steps scales as τ/Δ, a true bound B would multiply the longest-circuit Toffoli count (and the total sampling cost) by roughly sqrt(B / 1 Ha). The manuscript itself says 'Future work should validate this thesis more robustly.' I request either a more direct estimate of the bound—for example, an MPO/tensor-network evaluation of Y_3 on approximate core-excited eigenstates of the actual CAS(22e,18o) Hamiltonian—or a sensitivity table for B = 1, 10, 100 Ha, so the reader can assess how the headline 3.11e8 and the five-orders-of-magnitude claim degrade with the true value of the bound.
  3. [Sec. III.B and Appendix B] The argument that controlling the eigenvalue error is sufficient for spectroscopy ignores the first-order eigenvector corrections in Eq. (26) and Eq. (B6). The spectral intensity depends on |⟨mρ|E_l⟩|^2, and the eigenvector correction is proportional to Δ^2 ⟨E_k|Y_3|E_l⟩ / (E_k - E_l). In the core-excited manifold of a transition-metal cluster, energy denominators can be much smaller than the 1 eV eigenvalue tolerance; near-degeneracies would amplify the eigenvector error and distort peak intensities even when the eigenvalue shifts are within tolerance. The manuscript should either bound the eigenvector correction and its effect on the spectral function, or provide numerical evidence (on the smaller N2 and LiMnO benchmarks where exact or converged spectra are available) that peak-height errors remain below the target under the same Δ selection.
  4. [Sec. III.B vs. Sec. V.B and Table I] The perturbative analysis in Sec. III.B, Eqs. (22)-(27), is developed for the deterministic second-order Trotter formula, but the benchmarks and resource estimates use a 'randomized second-order product formula' (Sec. V.B and Table I caption). The manuscript does not explain how randomization is incorporated into the effective-Hamiltonian analysis, nor whether Eq. (41) remains a valid spectral-error bound for the randomized circuit. Please clarify the randomized construction and either extend the analysis or provide numerical verification that the randomized formula obeys the same (or better) eigenvalue-error scaling.
minor comments (4)
  1. [Abstract and Table I] The abstract states 'less than 4×10^8 T gates per circuit', but Table I and the body report Toffoli gates (3.11e8). Since a Toffoli gate is not the same as a T gate, please harmonize the terminology throughout.
  2. [Sec. IV.B and Sec. V.C] The target error is given as '1 eV' in Sec. IV.B.2, but Sec. V.C.6 sets η = 0.05 Ha (approximately 1.36 eV) and the text elsewhere quotes 1 eV ≈ 0.039 Ha. Please clarify whether the target eigenvalue tolerance and the line-broadening parameter are meant to be identical, and use a consistent conversion between atomic units and eV.
  3. [Sec. IV.B.3, Eq. (45)] The expression jmax = π/(2ητ) log(1/ϵ_trunc) appears to have an extra factor of π and a missing factor of π inside the logarithm relative to the bound derived in Appendix B.3, ϵ_trunc ≤ e^{-2jmax τη}/π. Since the discrepancy is conservative, it does not invalidate the resource numbers, but the equations should be reconciled.
  4. [Sec. V.C.6] The sentence 'For τ = π/4 Ha, this is satisfied for τ = 4Δ' is confusing because τ and Δ both have units of inverse energy; it would be clearer to write Δ = τ/4 and to give the value of Δ explicitly (Δ ≈ 0.196 Ha^{-1}).

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the resource estimate is conditional on an explicitly flagged extrapolated Trotter-error bound, which is an input, not a re-statement of the claimed outputs.

full rationale

The paper does not derive its central claims from those claims themselves. The resource estimate in Table I is an input-conditional calculation: given the empirically extrapolated bound |<Y3>| <= 1 Ha (Sec. V.C.6 and Fig. 13), Eq. (33) fixes the Trotter step count, and the cost formulas in Sec. IV.A translate that into Toffoli counts. That bound is an input to the estimate, not a restatement of the target spectrum or the gate count. The perturbation-theory Trotter error analysis in Sec. III.B (Eqs. 26-27) follows from the Baker-Campbell-Hausdorff expansion and is not defined in terms of the spectrum it later predicts. The algorithm is benchmarked against an independent classical full-configuration-interaction reference for N2 (Sec. V.B, Fig. 8), providing an external check of the time-domain method. The LiMnO simulation shows convergence with jmax but is explicitly presented without a classical reference. The paper itself flags the weakest link: 'Extrapolating from the growth of this error for small systems, we conclude that it is reasonable to expect |<El|Y3|El>| <= 1 Ha for the CAS(22e,18o) Li4Mn2O cluster. Future work should validate this thesis more robustly.' This is a correctness risk and an unvalidated assumption, not circular reasoning: a larger true bound would change the input Delta and grow the Toffoli count, but it would not make the derivation equivalent to its inputs. Self-citations (Ref. [11] for the base time-domain algorithm and Ref. [46] for initial-state preparation) supply components of the construction but are not invoked as uniqueness theorems or as the sole justification for any forbidden alternative; the central optimization claims (CDF, eigenvalue-error-based Trotter step selection, Lorentzian sampling, double measurement) are derived and analyzed in this paper. A likely typo in Eq. (42) (inverted ratio in the displayed inequality) does not affect the resource estimate because Table I and Sec. V.C.6 use the correct Delta = sqrt(eta/Y3). Overall, no step reduces by construction to its own inputs.

Assumptions & free parameters 9 free parameters · 7 assumptions · 0 invented entities

The central resource estimate depends on several assumed constants: the Trotter error matrix element bound (1 Ha), the CDF rank N, the shot budget, and the Slater count. None of these are fitted to the target XAS spectrum itself, but they are extrapolated from smaller systems or chosen heuristically. The paper contains no new physical entities.

free parameters (9)
  • Assumed Trotter error matrix element |<Y3>| = 1 Ha (assumed upper bound)
    Sec. V.C.6: extrapolated from small-system data in Fig. 13; sets the Trotter step Delta via Delta^2 * |<Y3>| <= 0.05 Ha. Not verified for the target CAS(22e,18o) system.
  • Number of CDF fragments L = 18 (= N)
    Sec. V.B and Fig. 11: L=N chosen to meet ~1 eV target error, extrapolated from orbital counts N=6 to 14.
  • Shot budget S = 2500
    Sec. V.C: set by measurement-noise analysis in Appendix B.2.
  • Maximal evolution time j_max = 100
    Sec. V.C: from truncation error bound Eq. (45) for eta=0.05 Ha, tau=pi/4 Ha^-1.
  • Time step tau = pi/4 Ha^-1
    Sec. V.C: from discretization error bound and ||H||_omega = 2 Ha.
  • Sampling exponent alpha = 1.3384
    Sec. V.C and Appendix B.2.a: optimized to minimize total shot cost; gives 1.178x saving over alpha=1.
  • Rotation precision epsilon_rot = 1e-3 Ha
    Sec. V.C.5 and Fig. 12: empirical trade-off; rotation synthesis via Gidney adder with ~10 Toffoli per rotation.
  • State preparation Slater count D = 10^4
    Sec. V.C: assumed sufficient to capture the state m|I>; no convergence evidence shown.
  • Double measurement savings factor = 1.49
    Appendix C.2: average savings over the unit circle; actual savings depend on the expectation value <rho|U|rho>.
assumptions (7)
  • domain assumption Fermi's golden rule and dipole approximation for XAS cross-section
    Eq. (1), Sec. II.A: physics model for absorption cross-section.
  • domain assumption Core-valence separation (CVS) approximation
    Sec. V.B: zeroing small two-electron integrals to decouple core-excited and valence-excited manifolds; error quantified in cited literature [45] but not for this system.
  • domain assumption Compressed double factorization with L ~ O(N) gives sufficient accuracy
    Sec. III.A and Appendix A; validated on small systems, extrapolated to N=18.
  • standard math First-order perturbation theory for eigenvalue shift is sufficient
    Sec. III.B and Appendix B.1; assumes the leading-order Trotter error dominates.
  • domain assumption Initial state |rho> = m|I>/||...|| can be prepared by sum-of-Slaters with D=10^4
    Sec. V.C and Ref [46]; no fidelity check for Li4Mn2O.
  • domain assumption Active volume architecture model from Ref [40] applies with quoted cycle costs
    Sec. IV.C, used to convert gate counts to runtime.
  • domain assumption Randomized second-order Trotter formula is optimal for the precision regime
    Sec. III.A.2: empirical comparison, not proven optimal.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Fast simulations of X-ray absorption spectroscopy for battery materials on a quantum computer." pith.science (2026). https://pith.science/paper/KOP4SX7R

@misc{pith2026250615784,
  author       = {Pith},
  title        = {Pith review of: Fast simulations of X-ray absorption spectroscopy for battery materials on a quantum computer},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/KOP4SX7R}},
  note         = {Machine review of arXiv:2506.15784}
}
abstract

X-ray absorption spectroscopy (XAS) is a leading technique for understanding structural changes in advanced battery materials such as lithium-excess cathodes. However, extracting critical information like oxidation states from the experimental spectra requires expensive and time-consuming simulations. Building upon a recent proposal to simulate XAS using quantum computers, this work proposes a highly-optimized implementation of the time-domain algorithm for X-ray absorption. Among a host of improvements to Hamiltonian representation, circuit implementation, and measurement strategies, three optimizations are key to the efficiency of the algorithm. The first is the use of product formulas with the compressed double factorized form of the Hamiltonian. The second is recognizing that for spectroscopy applications, it is sufficient to control the error in the eigenvalues of the (approximate) Hamiltonian being implemented by the product formula, rather than the generic error on the full time evolution operator. Using perturbation theory to estimate this eigenvalue error, we find that significantly fewer Trotter steps are needed than expected from the time evolution error bound. The third is the choice of an optimized distribution of samples that takes advantage of the exponentially decaying Lorentzian kernel. Through constant factor resource estimates, we show that a challenging model Li$_4$Mn$_2$O cluster system with 18 spatial orbitals and 22 electrons in the active space can be simulated with 100 logical qubits and less than $4 \times 10^8$ T gates per circuit. Finally, the algorithm is implemented on a simulator, and the reconstructed spectrum is verified against a classical computational reference. The low cost of our algorithm makes it attractive to use on fault-tolerant quantum devices to accelerate the development and commercialization of high-capacity battery cathodes.

Figures

Figures reproduced from arXiv: 2506.15784 by the authors.

Figure 1
Figure 1. FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Schematic of electronic excitations relevant in XAS [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Circuit evaluating the time evolution expectation [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figures from the paper (9 more)
Figure 6
Figure 6. Figure 6: FIG. 6. Double phase trick: we can use cheap CNOT gates [PITH_FULL_IMAGE:figures/full_fig_p006_6.png]
Figure 5
Figure 5. Figure 5: FIG. 5. Mapping a two-qubit [PITH_FULL_IMAGE:figures/full_fig_p006_5.png]
Figure 7
Figure 7. Figure 7: FIG. 7. Double measurement circuit, replacing Fig. 3, where [PITH_FULL_IMAGE:figures/full_fig_p008_7.png]
Figure 8
Figure 8. Figure 8: FIG. 8. Benchmarking the XAS algorithm: valence-excited [PITH_FULL_IMAGE:figures/full_fig_p012_8.png]
Figure 9
Figure 9. Figure 9: FIG. 9. Benchmarking the XAS algorithm: core-excited [PITH_FULL_IMAGE:figures/full_fig_p013_9.png]
Figure 10
Figure 10. Figure 10: FIG. 10. Comparing expected runtimes of our proposed quan [PITH_FULL_IMAGE:figures/full_fig_p015_10.png]
Figure 11
Figure 11. Figure 11: FIG. 11. Mean and standard deviation of the error of the [PITH_FULL_IMAGE:figures/full_fig_p022_11.png]
Figure 13
Figure 13. Figure 13: FIG. 13. Error of the second and fourth order formulas in the first five eigenvalues of CDF Hamiltonians for active spaces [PITH_FULL_IMAGE:figures/full_fig_p023_13.png]
Figure 14
Figure 14. Figure 14: FIG. 14. Double measurement Hadamard circuit for non-unitary state preparation. It performs a [PITH_FULL_IMAGE:figures/full_fig_p027_14.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Unified Uncertainty Quantification Framework Bridging Noisy Quantum Backends Across Variational Quantum Algorithms and Quantum Signal Processing

    cs.ET 2026-07 conditional novelty 6.0 of 10

    A unified Bayesian-optimization and uncertainty-quantification workflow benchmarks variational and QSVT workloads on four noisy backends, showing backend quality is workload-dependent.

Reference graph

Works this paper leans on

114 extracted references · 61 canonical work pages · cited by 1 Pith paper

  1. [1]

    Our empirical results, shown in Fig

    Number of fragments in the CDF:The number of fragments L needed to achieve a fixed accuracy scales approximately linearly with the number of orbitals L = ˜O(N ) [52–54]. Our empirical results, shown in Fig. 11 in Appendix A, also validate this scaling. Based on these results, we estimate that we need to choose rank L = N in Eq. (18) to achieve the target ...

  2. [2]

    As derived in Eqs

    T rotter time step:The error source determining the choice of the Trotter time step ∆ is the Trotter error. As derived in Eqs. (26) and (27) the error in the eigenvalues is ϵtrot = ∆2k ⟨El|Y2k+1|El⟩ , (41) to first order in perturbation theory, for a 2 k or- der product formula whose leading order term in the effective Hamiltonian is Y2k+1, see Eqs. (24) ...

  3. [3]

    They are associated to variables jmax and τ respectively

    Evolution time and time step:Two more er- ror sources are related to the truncation and dis- cretization of the Fourier transform integral. They are associated to variables jmax and τ respectively. The error ϵtrunc arises from truncating the formally infinite boundaries of the Fourier transform inte- gral to a maximal evolution time jmax: it can be bounde...

  4. [4]

    Let us distribute the shot budget according to Eq

    Shots budget:The last source of error is the mea- surement noise in each Hadamard test, determining the required shot budget S. Let us distribute the shot budget according to Eq. (31). Then, the num- ber of shots needed to achieve a measurement error ϵmeas is S = ητ 2πϵmeas 2 X j e(α−2)|j|τ ηX k e−α|k|τ η, (47) see Appendix B 2. Meanwhile the uniform dist...

  5. [5]

    CDF: Using compressed double factorization with L = N unlocks significant savings. A sparsi- fied Jordan-Wigner mapping of the CAS(22e, 18o) Hamiltonian results in approximately 5 .7 × 104 Pauli string rotations with angles larger than 10 −2 for a single first-order Trotter step with error close to 1 eV. By contrast, the CDF Hamiltonian with rank L = N on...

  6. [6]

    Symmetries (BLISS):For the active space with N = 18, the benefit in cost to shaving off one frag- ment from the CDF is ( L + 1)/L = 19/18 ≈ 1.1

  7. [7]

    Double phase: As argued in Sec. III, replac- ing controlled rotations by un-controlled rotations saves a factor of 2 at the cost of two Clifford gates (CNOTs), since a controlled rotation can be imple- mented via 2 uncontrolled rotations; moreover, the double phase approach doubles the effective num- ber of steps being simulated at no extra cost. Thus the...

  8. [8]

    This allows to estimate the cost savings from having one less separate unitary to apply per Trot- ter step, resulting in a multiplicative saving factor of ×1.46

    Combining consecutive rotations: Using the CDF with L = N for the CAS(22e, 18o) Li 4Mn2O cluster, we evaluate the relative costs of implement- ing the unitaries U (ℓ) versus the Z (ℓ) kl Pauli Z rota- tions. This allows to estimate the cost savings from having one less separate unitary to apply per Trot- ter step, resulting in a multiplicative saving fact...

Show all 114 references
  1. [9]

    12, by studying a range of active spaces for the Li4Mn2O cluster, we empirically find that ϵrot = 10−3 results in an eigenvalue error on the order of 10 −2 Ha

    Rotation precision: As demonstrated in Fig. 12, by studying a range of active spaces for the Li4Mn2O cluster, we empirically find that ϵrot = 10−3 results in an eigenvalue error on the order of 10 −2 Ha. Considering our target accuracy and other error sources, we thus choose ϵ...

  2. [10]

    13 shows the second and fourth order Trotter error for several chemical systems

    Perturbative analysis of the spectrum:Fig. 13 shows the second and fourth order Trotter error for several chemical systems. We estimated this error by evaluating the difference between time evolu- tions for the same t = 1, but different number of Trotter steps. Using the secon...

  3. [11]

    IV B we find a multiplicative savings factor of ×17.9 relative to uniform sampling

    Sampling distribution: For the choice of param- eters jmax = 100 of the maximal evolution time, τ = π/4 Ha , η = 0.05 Ha as above, and N = 18, using the expressions derived in Sec. IV B we find a multiplicative savings factor of ×17.9 relative to uniform sampling

  4. [12]

    III D, we can extract further informa- tion from the residual state of a Hadamard test using this double measurement approach

    Double measurement As we explained in Sec. III D, we can extract further informa- tion from the residual state of a Hadamard test using this double measurement approach. On average, this provides a ×1.49 savings factor. D. Classical vs. quantum In this concluding section, we u...

  5. [13]

    K. Luo, M. R. Roberts, R. Hao, N. Guerrini, D. M. Pickup, Y.-S. Liu, K. Edstr¨ om, J. Guo, A. V. Chad- wick, L. C. Duda, et al., Charge-compensation in 3 d-transition-metal-oxide intercalation cathodes through the generation of localized electron holes on oxygen, Na- ture chem...

    2016

  6. [14]

    Lu and J

    Z. Lu and J. R. Dahn, Understanding the anomalous ca- pacity of li/li [ni x li (1/3- 2x/3) mn (2/3- x/3)] o 2 cells using in situ x-ray diffraction and electrochemical stud- ies, Journal of the Electrochemical Society 149, A815 (2002)

    2002

  7. [15]

    Z. Lu, L. Beaulieu, R. Donaberger, C. Thomas, and J. Dahn, Synthesis, structure, and electrochemical be- havior of li [ni x li1/3- 2x/3mn2/3- x/3] o 2, Journal of the Electrochemical Society 149, A778 (2002)

    2002

  8. [16]

    Sauban` ere, E

    M. Sauban` ere, E. McCalla, J.-M. Tarascon, and M.-L. Doublet, The intriguing question of anionic redox in high- energy density cathodes for li-ion batteries, Energy & Environmental Science 9, 984 (2016)

    2016

  9. [17]

    N. Tran, L. Croguennec, M. M´ en´ etrier, F. Weill, P. Bi- ensan, C. Jordy, and C. Delmas, Mechanisms associated with the “plateau” observed at high voltage for the over- lithiated li1. 12 (ni0. 425mn0. 425co0. 15) 0.88 o2 system, Chemistry of Materials 20, 4815 (2008)

    2008

  10. [18]

    Zhang, D

    M. Zhang, D. A. Kitchaev, Z. Lebens-Higgins, J. Vinck- eviciute, M. Zuba, P. J. Reeves, C. P. Grey, M. S. Whit- tingham, L. F. Piper, A. Van der Ven, et al., Pushing the limit of 3 d transition metal-based layered oxides that use both cation and anion redox for energy storage,...

    2022

  11. [19]

    J. Hong, H. Gwon, S.-K. Jung, K. Ku, and K. Kang, Lithium-excess layered cathodes for lithium rechargeable batteries, Journal of The Electrochemical Society 162, A2447 (2015)

    2015

  12. [20]

    S. Hy, H. Liu, M. Zhang, D. Qian, B.-J. Hwang, and Y. S. Meng, Performance and design considerations for lithium excess layered oxide positive electrode materials for lithium ion batteries, Energy & Environmental Sci- ence 9, 1931 (2016)

    2016

  13. [21]

    M. D. Radin, S. Hy, M. Sina, C. Fang, H. Liu, J. Vincke- viciute, M. Zhang, M. S. Whittingham, Y. S. Meng, and A. Van der Ven, Narrowing the gap between theoretical and practical capacities in li-ion layered oxide cathode materials, Advanced Energy Materials7, 1602888 (2017)

    2017

  14. [22]

    C. S. Johnson, N. Li, C. Lefief, J. T. Vaughey, and M. M. Thackeray, Synthesis, characterization and elec- trochemistry of lithium battery electrodes: x li2mno3 ·(1- x) limn0. 333ni0. 333co0. 333o2 (0 ≤ x≤ 0.7), Chemistry of Materials 20, 6095 (2008)

    2008

  15. [23]

    Bettge, Y

    M. Bettge, Y. Li, K. Gallagher, Y. Zhu, Q. Wu, W. Lu, I. Bloom, and D. P. Abraham, Voltage fade of layered oxides: its measurement and impact on energy den- sity, Journal of The Electrochemical Society 160, A2046 (2013)

    2013

  16. [24]

    Fomichev, K

    S. Fomichev, K. Hejazi, I. Loaiza, M. S. Zini, A. Delgado, A.-C. Voigt, J. E. Mueller, and J. M. Arrazola, Simulat- ing x-ray absorption spectroscopy of battery materials on a quantum computer, arXiv preprint arXiv:2405.11015 (2024). 17

    2024 arXiv

  17. [25]

    D.-H. Seo, J. Lee, A. Urban, R. Malik, S. Kang, and G. Ceder, The structural and chemical origin of the oxy- gen redox activity in layered and cation-disordered li- excess cathode materials, Nature chemistry8, 692 (2016)

    2016

  18. [26]

    P. J. Lestrange, P. D. Nguyen, and X. Li, Calibration of energy-specific tddft for modeling k-edge xas spectra of light elements, Journal of Chemical Theory and Compu- tation 11, 2994 (2015)

    2015

  19. [27]

    H. Koga, L. Croguennec, M. M´ en´ etrier, P. Mannessiez, F. Weill, and C. Delmas, Different oxygen redox partic- ipation for bulk and surface: A possible global expla- nation for the cycling mechanism of li1. 20mn0. 54co0. 13ni0. 13o2, Journal of Power Sources 236, 250 (2013)

    2013

  20. [28]

    H. Koga, L. Croguennec, M. M´ en´ etrier, P. Mannessiez, F. Weill, C. Delmas, and S. Belin, Operando x-ray ab- sorption study of the redox processes involved upon cy- cling of the li-rich layered oxide li1. 20mn0. 54co0. 13ni0. 13o2 in li ion batteries, The Journal of Physical...

    2014

  21. [29]

    R. A. House, G. J. Rees, K. McColl, J.-J. Marie, M. Garcia-Fernandez, A. Nag, K.-J. Zhou, S. Cassidy, B. J. Morgan, M. Saiful Islam, et al., Delocalized electron holes on oxygen in a battery cathode, Nature Energy 8, 351 (2023)

    2023

  22. [30]

    Kalyani, S

    P. Kalyani, S. Chitra, T. Mohan, and S. Gopukumar, Lithium metal rechargeable cells using li2mno3 as the positive electrode, Journal of power sources 80, 103 (1999)

    1999

  23. [31]

    Ohzuku, M

    T. Ohzuku, M. Nagayama, K. Tsuji, and K. Ariyoshi, High-capacity lithium insertion materials of lithium nickel manganese oxides for advanced lithium-ion bat- teries: toward rechargeable capacity more than 300 ma hg- 1, Journal of Materials chemistry 21, 10179 (2011)

    2011

  24. [32]

    M. D. Radin, J. Vinckeviciute, R. Seshadri, and A. Van der Ven, Manganese oxidation as the origin of the anomalous capacity of mn-containing li-excess cath- ode materials, Nature Energy 4, 639 (2019)

    2019

  25. [33]

    De Groot, High-resolution x-ray emission and x-ray absorption spectroscopy, Chemical Reviews 101, 1779 (2001)

    F. De Groot, High-resolution x-ray emission and x-ray absorption spectroscopy, Chemical Reviews 101, 1779 (2001)

    2001

  26. [34]

    Stener, G

    M. Stener, G. Fronzoni, and M. de Simone, Time depen- dent density functional theory of core electrons excita- tions, Chemical physics letters 373, 115 (2003)

    2003

  27. [35]

    K. Ray, S. DeBeer George, E. I. Solomon, K. Wieghardt, and F. Neese, Description of the ground-state covalencies of the bis (dithiolato) transition-metal complexes from x- ray absorption spectroscopy and time-dependent density- functional calculations, Chemistry–A European Jou...

    2007

  28. [36]

    N. A. Besley and F. A. Asmuruf, Time-dependent den- sity functional theory calculations of the spectroscopy of core electrons, Physical Chemistry Chemical Physics 12, 12024 (2010)

    2010

  29. [37]

    Liang, S

    W. Liang, S. A. Fischer, M. J. Frisch, and X. Li, Energy- specific linear response tdhf/tddft for calculating high- energy excited states, Journal of chemical theory and computation 7, 3540 (2011)

    2011

  30. [38]

    Lopata, B

    K. Lopata, B. E. Van Kuiken, M. Khalil, and N. Govind, Linear-response and real-time time-dependent density functional theory studies of core-level near-edge x-ray ab- sorption, Journal of chemical theory and computation 8, 3284 (2012)

    2012

  31. [39]

    D. Casanova, Restricted active space configuration inter- action methods for strong correlation: Recent develop- ments, Wiley Interdisciplinary Reviews: Computational Molecular Science 12, e1561 (2022). 18

    2022

  32. [40]

    T. F. Stetina, J. M. Kasper, and X. Li, Modeling l2, 3- edge x-ray absorption spectroscopy with linear response exact two-component relativistic time-dependent density functional theory, The Journal of Chemical Physics 150 (2019)

    2019

  33. [41]

    B. O. Roos, P. R. Taylor, and P. E. Sigbahn, A complete active space scf method (casscf) using a density matrix formulated super-ci approach, Chemical Physics 48, 157 (1980)

    1980

  34. [42]

    Casanova and M

    D. Casanova and M. Head-Gordon, Restricted active space spin-flip configuration interaction approach: the- ory, implementation and examples, Physical Chemistry Chemical Physics 11, 9779 (2009)

    2009

  35. [43]

    Maganas, P

    D. Maganas, P. Kristiansen, L.-C. Duda, A. Knop- Gericke, S. DeBeer, R. Schlo¨ ogl, and F. Neese, Combined experimental and ab initio multireference configuration interaction study of the resonant inelastic x-ray scatter- ing spectrum of co2, The Journal of Physical Chemistry ...

    2014

  36. [44]

    Y. Yang, M. A. Ratner, and G. C. Schatz, Multirefer- ence ab initio study of ligand field d–d transitions in oc- tahedral transition-metal oxide clusters, The Journal of Physical Chemistry C 118, 29196 (2014)

    2014

  37. [45]

    R. V. Pinjari, M. G. Delcey, M. Guo, M. Odelius, and M. Lundberg, Restricted active space calculations of l- edge x-ray absorption spectra: From molecular orbitals to multiplet states, The Journal of chemical physics 141 (2014)

    2014

  38. [46]

    M. Guo, L. K. Sørensen, M. G. Delcey, R. V. Pinjari, and M. Lundberg, Simulations of iron k pre-edge x-ray ab- sorption spectra using the restricted active space method, Physical Chemistry Chemical Physics 18, 3250 (2016)

    2016

  39. [47]

    Sassi, C

    M. Sassi, C. I. Pearce, P. S. Bagus, E. Arenholz, and K. M. Rosso, First-principles fe l2, 3-edge and o k-edge xanes and xmcd spectra for iron oxides, The Journal of Physical Chemistry A 121, 7613 (2017)

    2017

  40. [48]

    M. G. Delcey, L. K. Sørensen, M. Vacher, R. C. Couto, and M. Lundberg, Efficient calculations of a large num- ber of highly excited states for multiconfigurational wave- functions, Journal of computational chemistry 40, 1789 (2019)

    2019

  41. [49]

    Maganas, J

    D. Maganas, J. K. Kowalska, M. Nooijen, S. DeBeer, and F. Neese, Comparison of multireference ab initio wave- function methodologies for x-ray absorption edges: A case study on [fe (ii/iii) cl4] 2–/1–molecules, The Journal of Chemical Physics 150 (2019)

    2019

  42. [50]

    M. Guo, X. Liu, and R. He, Restricted active space sim- ulations of the metal l-edge x-ray absorption spectra and resonant inelastic x-ray scattering: revisiting [co ii/iii (bpy) 3] 2+/3+ complexes, Inorganic Chemistry Fron- tiers 7, 1927 (2020)

    2020

  43. [51]

    S. Lee, H. Zhai, and G. K.-L. Chan, An ab initio correc- tion vector restricted active space approach to the l-edge xas and 2p3d rixs spectra of transition metal complexes, Journal of Chemical Theory and Computation 19, 7753 (2023)

    2023

  44. [52]

    Motta, E

    M. Motta, E. Ye, J. R. McClean, Z. Li, A. J. Minnich, R. Babbush, and G. K.-L. Chan, Low rank representa- tions for quantum simulation of electronic structure, npj Quantum Information 7, 83 (2021)

    2021

  45. [53]

    Litinski and N

    D. Litinski and N. Nickerson, Active volume: An ar- chitecture for efficient fault-tolerant quantum comput- ers with limited non-local connections, arXiv preprint arXiv:2211.15465 (2022)

    2022 arXiv

  46. [54]

    J. Cohn, M. Motta, and R. M. Parrish, Quantum filter di- agonalization with compressed double-factorized hamil- tonians, PRX Quantum 2, 040352 (2021)

    2021

  47. [55]

    L. S. Cederbaum, W. Domcke, and J. Schirmer, Many- body theory of core holes, Physical Review A 22, 206 (1980)

    1980

  48. [56]

    Barth and L

    A. Barth and L. Cederbaum, Many-body theory of core- valence excitations, Physical Review A 23, 1038 (1981)

    1981

  49. [57]

    Norman and A

    P. Norman and A. Dreuw, Simulating x-ray spectro- scopies and calculating core-excited states of molecules, Chemical reviews 118, 7208 (2018)

    2018

  50. [58]

    M. F. Herbst and T. Fransson, Quantifying the error of the core–valence separation approximation, The Journal of Chemical Physics 153 (2020)

    2020

  51. [59]

    Fomichev, K

    S. Fomichev, K. Hejazi, M. S. Zini, M. Kiser, J. F. Morales, P. A. M. Casares, A. Delgado, J. Huh, A.- C. Voigt, J. E. Mueller, et al., Initial state preparation for quantum chemistry on quantum computers, arXiv preprint arXiv:2310.18410 (2023)

    2023

  52. [60]

    Yen and A

    T.-C. Yen and A. F. Izmaylov, Cartan subalgebra ap- proach to efficient measurements of quantum observables, PRX quantum 2, 040320 (2021)

    2021

  53. [61]

    Oumarou, M

    O. Oumarou, M. Scheurer, R. M. Parrish, E. G. Hohen- stein, and C. Gogolin, Accelerating quantum computa- tions of chemistry through regularized compressed dou- ble factorization, Quantum 8, 1371 (2024)

    2024

  54. [62]

    I. D. Kivlichan, J. McClean, N. Wiebe, C. Gidney, A. Aspuru-Guzik, G. K.-L. Chan, and R. Babbush, Quantum simulation of electronic structure with linear depth and connectivity, Physical Review Letters 120, 110501 (2018)

    2018

  55. [63]

    D. J. Thouless, Stability conditions and nuclear rota- tions in the hartree-fock theory, Nuclear Physics 21, 225 (1960)

    1960

  56. [64]

    J. M. Arrazola, O. Di Matteo, N. Quesada, S. Jahangiri, A. Delgado, and N. Killoran, Universal quantum circuits for quantum chemistry, Quantum 6, 742 (2022)

    2022

  57. [65]

    A. C. Vazquez, D. J. Egger, D. Ochsner, and S. Woerner, Well-conditioned multi-product formulas for hardware- friendly hamiltonian simulation, Quantum 7, 1067 (2023)

    2023

  58. [66]

    von Burg, G

    V. von Burg, G. H. Low, T. H¨ aner, D. S. Steiger, M. Rei- her, M. Roetteler, and M. Troyer, Quantum computing enhanced computational catalysis, Physical Review Re- search 3, 033055 (2021)

    2021

  59. [67]

    J. Lee, D. W. Berry, C. Gidney, W. J. Huggins, J. R. Mc- Clean, N. Wiebe, and R. Babbush, Even more efficient quantum computations of chemistry through tensor hy- percontraction, PRX Quantum 2, 030305 (2021)

    2021

  60. [68]

    P. G. Anastasiou, Y. Chen, N. J. Mayhall, E. Barnes, and S. E. Economou, Tetris-adapt-vqe: an adaptive algo- rithm that yields shallower, denser circuit ans¨ atze, Phys- ical Review Research 6, 013254 (2024)

    2024

  61. [69]

    Loaiza and A

    I. Loaiza and A. F. Izmaylov, Block-invariant symme- try shift: Preprocessing technique for second-quantized hamiltonians to improve their decompositions to linear combination of unitaries, Journal of Chemical Theory and Computation 19, 8201 (2023)

    2023

  62. [70]

    Rocca, C

    D. Rocca, C. L. Cortes, J. F. Gonthier, P. J. Ollitrault, R. M. Parrish, G.-L. Anselmetti, M. Degroote, N. Moll, R. Santagati, and M. Streif, Reducing the runtime of fault-tolerant quantum simulations in chemistry through symmetry-compressed double factorization, Journal of Ch...

    2024

  63. [71]

    Patel, A

    S. Patel, A. S. Brahmachari, J. T. Cantin, L. Wang, and A. F. Izmaylov, Global minimization of electronic hamil- tonian 1-norm via linear programming in the block invari- ant symmetry shift (bliss) method, Journal of Chemical Theory and Computation (2025)

    2025

  64. [72]

    Caesura, C

    A. Caesura, C. L. Cortes, W. Poll, S. Sim, M. Steudtner, G.-L. R. Anselmetti, M. Degroote, N. Moll, R. Santa- gati, M. Streif, et al., Faster quantum chemistry simula- tions on a quantum computer with improved tensor fac- torization and active volume compilation, arXiv preprin...

    2025 arXiv

  65. [73]

    Suzuki, General theory of higher-order decomposi- tion of exponential operators and symplectic integrators, Physics Letters A 165, 387 (1992)

    M. Suzuki, General theory of higher-order decomposi- tion of exponential operators and symplectic integrators, Physics Letters A 165, 387 (1992)

    1992

  66. [74]

    G. H. Low, V. Kliuchnikov, and N. Wiebe, Well- conditioned multiproduct hamiltonian simulation, arXiv preprint arXiv:1907.11679 (2019)

    2019 arXiv

  67. [75]

    Blanes, F

    S. Blanes, F. Casas, and L. Shaw, Generalized extrapo- lation methods based on compositions of a basic 2nd- order scheme, Applied Mathematics and Computation 473, 128663 (2024)

    2024

  68. [76]

    S. Endo, Q. Zhao, Y. Li, S. Benjamin, and X. Yuan, Mit- igating algorithmic errors in a hamiltonian simulation, Physical Review A 99, 012334 (2019)

    2019

  69. [77]

    A. C. Vazquez, R. Hiptmair, and S. Woerner, Enhancing the quantum linear systems algorithm using richardson extrapolation, ACM Transactions on Quantum Comput- ing 3, 1 (2022)

    2022

  70. [78]

    H¨ aner, D

    T. H¨ aner, D. S. Steiger, K. Svore, and M. Troyer, A software methodology for compiling quantum programs, Quantum Science and Technology 3, 020501 (2018)

    2018

  71. [79]

    S. Zhuk, N. Robertson, and S. Bravyi, Trotter error bounds and dynamic multi-product formulas for hamilto- nian simulation, arXiv preprint arXiv:2306.12569 (2023)

    2023 arXiv

  72. [80]

    N. F. Robertson, B. Pokharel, B. Fuller, E. Switzer, O. Shtanko, M. Amico, A. Byrne, A. D’Urbano, S. Hayes- Shuptar, A. Akhriev, N. Keenan, S. Bravyi, and S. Zhuk, Tensor network enhanced dynamic multiproduct formu- las (2024), arXiv:2407.17405 [quant-ph]

    2024 arXiv

  73. [81]

    Rendon, J

    G. Rendon, J. Watkins, and N. Wiebe, Improved accu- racy for trotter simulations using chebyshev interpola- tion, Quantum 8, 1266 (2024)

    2024

  74. [82]

    Campbell, Random compiler for fast hamiltonian sim- ulation, Physical review letters 123, 070503 (2019)

    E. Campbell, Random compiler for fast hamiltonian sim- ulation, Physical review letters 123, 070503 (2019)

    2019

  75. [83]

    Hagan and N

    M. Hagan and N. Wiebe, Composite quantum simula- tions, Quantum 7, 1181 (2023)

    2023

  76. [84]

    Pocrnic, M

    M. Pocrnic, M. Hagan, J. Carrasquilla, D. Segal, and N. Wiebe, Composite qdrift-product formulas for quan- tum and classical simulations in real and imaginary time, Physical Review Research 6, 013224 (2024)

    2024

  77. [85]

    Blanes, F

    S. Blanes, F. Casas, and J. Ros, Processing symplectic methods for near-integrable hamiltonian systems, Celes- tial Mechanics and Dynamical Astronomy 77, 17 (2000)

    2000

  78. [86]

    Blanes, F

    S. Blanes, F. Casas, and A. Murua, On the numerical in- tegration of ordinary differential equations by processed methods, SIAM Journal on Numerical Analysis 42, 531 (2004)

    2004

  79. [87]

    Sharma and M

    K. Sharma and M. C. Tran, Hamiltonian simulation in the interaction picture using the magnus expansion, 19 arXiv preprint arXiv:2404.02966 (2024)

    2024 arXiv

  80. [88]

    J. L. Bosse, A. M. Childs, C. Derby, F. M. Gambetta, A. Montanaro, and R. A. Santos, Efficient and practi- cal hamiltonian simulation from time-dependent product formulas, arXiv preprint arXiv:2403.08729 (2024)

    2024 arXiv

  81. [89]

    Strang, On the construction and comparison of dif- ference schemes, SIAM journal on numerical analysis 5, 506 (1968)

    G. Strang, On the construction and comparison of dif- ference schemes, SIAM journal on numerical analysis 5, 506 (1968)

    1968

  82. [90]

    A. M. Childs, A. Ostrander, and Y. Su, Faster quantum simulation by randomization, Quantum 3, 182 (2019)

    2019

  83. [91]

    Asadi, A

    A. Asadi, A. Dusko, C.-Y. Park, V. Michaud-Rioux, I. Schoch, S. Shu, T. Vincent, and L. J. O’Riordan, Hy- brid quantum programming with PennyLane Lightning on HPC platforms (2024), arXiv:2403.02512 [quant-ph]

    2024 arXiv

  84. [92]

    Babbush, C

    R. Babbush, C. Gidney, D. W. Berry, N. Wiebe, J. Mc- Clean, A. Paler, A. Fowler, and H. Neven, Encoding elec- tronic spectra in quantum circuits with linear t complex- ity, Physical Review X 8, 041015 (2018)

    2018

  85. [93]

    Lin and Y

    L. Lin and Y. Tong, Heisenberg-limited ground-state en- ergy estimation for early fault-tolerant quantum comput- ers, PRX Quantum 3, 010318 (2022)

    2022

  86. [94]

    Y. Dong, L. Lin, and Y. Tong, Ground-state prepara- tion and energy estimation on early fault-tolerant quan- tum computers via quantum eigenvalue transformation of unitary matrices, PRX Quantum 3, 040305 (2022)

    2022

  87. [95]

    Trenev, P

    D. Trenev, P. J. Ollitrault, S. M. Harwood, T. P. Gu- jarati, S. Raman, A. Mezzacapo, and S. Mostame, Refin- ing resource estimation for the quantum computation of vibrational molecular spectra through trotter error anal- ysis, Quantum 9, 1630 (2025)

    2025

  88. [96]

    Blanes and F

    S. Blanes and F. Casas, On the convergence and opti- mization of the baker–campbell–hausdorff formula, Lin- ear algebra and its applications 378, 135 (2004)

    2004

  89. [97]

    Gidney, Halving the cost of quantum addition, Quan- tum 2, 74 (2018)

    C. Gidney, Halving the cost of quantum addition, Quan- tum 2, 74 (2018)

    2018

  90. [98]

    Litinski, A game of surface codes: Large-scale quan- tum computing with lattice surgery, Quantum 3, 128 (2019)

    D. Litinski, A game of surface codes: Large-scale quan- tum computing with lattice surgery, Quantum 3, 128 (2019)

    2019

  91. [99]

    E. R. Sayfutyarova, Q. Sun, G. K.-L. Chan, and G. Knizia, Automated construction of molecular active spaces from atomic valence orbitals, Journal of chemical theory and computation 13, 4063 (2017)

    2017

  92. [100]

    Sun, Libcint: An efficient general integral library for g aussian basis functions, Journal of computational chem- istry 36, 1664 (2015)

    Q. Sun, Libcint: An efficient general integral library for g aussian basis functions, Journal of computational chem- istry 36, 1664 (2015)

    2015

  93. [101]

    Q. Sun, T. C. Berkelbach, N. S. Blunt, G. H. Booth, S. Guo, Z. Li, J. Liu, J. D. McClain, E. R. Sayfutyarova, S. Sharma, et al., Pyscf: the python-based simulations of chemistry framework, Wiley Interdisciplinary Reviews: Computational Molecular Science 8, e1340 (2018)

    2018

  94. [102]

    Q. Sun, X. Zhang, S. Banerjee, P. Bao, M. Barbry, N. S. Blunt, N. A. Bogdanov, G. H. Booth, J. Chen, Z.-H. Cui, et al., Recent developments in the pyscf program package, The Journal of chemical physics 153 (2020)

    2020

  95. [103]

    Bergholm, J

    V. Bergholm, J. Izaac, M. Schuld, C. Gogolin, S. Ahmed, V. Ajith, M. S. Alam, G. Alonso-Linaje, B. Akash- Narayanan, A. Asadi, et al., Pennylane: Automatic dif- ferentiation of hybrid quantum-classical computations, arXiv preprint arXiv:1811.04968 (2018)

    2018 arXiv

  96. [105]

    H. Gao, S. Imamura, A. Kasagi, and E. Yoshida, Dis- tributed implementation of full configuration interaction for one trillion determinants, Journal of Chemical Theory and Computation 20, 1185 (2024)

    2024

  97. [106]

    Zwiebach, Perturbation theory (2018)

    B. Zwiebach, Perturbation theory (2018)

    2018

  98. [107]

    A. M. Childs, Y. Su, M. C. Tran, N. Wiebe, and S. Zhu, Theory of trotter error with commutator scaling, Physi- cal Review X 11, 011020 (2021)

    2021

  99. [108]

    N. J. Ross and P. Selinger, Optimal ancilla-free clifford+t approximation of z-rotations, Quantum Info. Comput. 16, 901–953 (2016)

    2016

  100. [109]

    Kliuchnikov, K

    V. Kliuchnikov, K. Lauter, R. Minko, A. Paetznick, and C. Petit, Shorter quantum circuits via single-qubit gate approximation, Quantum 7, 1208 (2023). 20 Appendix A: Compressed Double F actorization Here we present the mathematical description of the Compressed Double Factoriz...

    2023

  101. [110]

    Here we make the case for taking a constant number of Trotter steps with each j

    T rotter error In the sections above we have considered that we want to simulate each j samples with a constant amount of error (or increasing with the sampling error). Here we make the case for taking a constant number of Trotter steps with each j. In other words, we are not ...

  102. [111]

    Measurement noise Now, we want to understand what is the effect of the sampling error in the result. Recall that the variance of a binomial distribution is given by σ2 = p(1 − p) N ≤ 1 4N (B15) The noise in the signal can be characterized by measuring the real and imaginary co...

  103. [112]

    The first one is related to truncating the Fourier transform integral to the interval τ j∈ (−τ jmax, +τ jmax) instead of τ j∈ (−∞, +∞)

    Modelling discrete-time F ourier transform error There are two more sources of error. The first one is related to truncating the Fourier transform integral to the interval τ j∈ (−τ jmax, +τ jmax) instead of τ j∈ (−∞, +∞). The error can be bounded by considering the norm of all...

  104. [113]

    Non-unitary state preparation Instead of separately measuring the real and imaginary components in a Hadamard test, we can measure the cost of one almost directly from the other, as shown in Fig. 14. The trick is to recycle the output of the Hadamard test, either |ψℜ⟩ = (1 ± U...

  105. [114]

    7, the probability of measuring a 0 in the second register is P±(0) = | ⟨ψℜ|ρ⟩ |2 = |1 ± ⟨ρ|U |ρ⟩ |2 2(1 ± ℜ ⟨ρ|U |ρ⟩)

    Unitary state preparation In the circuit depicted in the main text, Fig. 7, the probability of measuring a 0 in the second register is P±(0) = | ⟨ψℜ|ρ⟩ |2 = |1 ± ⟨ρ|U |ρ⟩ |2 2(1 ± ℜ ⟨ρ|U |ρ⟩) . (C22) As a consequence, |1 ± ⟨ρ|U |ρ⟩ |2 = 2P±(0)(1 ± ℜ(⟨ρ|U |ρ⟩)) (C23) ⇒ |1 ± (ℜ ...

  106. [1000]

    Combining all of this, we get a classical runtime = 0 .76 · 1000 hours ∼ 32 days

    Finally, we make the assumption that CVS, RAS and a projection-like approach is employed to eliminate intermediate states, so that the solver is directly com- puting the (core-excited) states of interest. Combining all of this, we get a classical runtime = 0 .76 · 1000 hours ∼...

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.