Pith. sign in

REVIEW 2 major objections 4 minor 36 references

A trapped-ion hybrid protocol simulates intense-field QED pair production with polynomial resources and noise mitigation that recovers the key signals.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · grok-4.5

2026-07-14 15:05 UTC pith:NKNA2WTF

load-bearing objection Solid hybrid mapping of Furry-picture IFQED onto ion spin-phonon gates, with a clean single-mode Breit-Wheeler benchmark; multimode noise remains untested. the 2 major comments →

arxiv 2607.09844 v1 pith:NKNA2WTF submitted 2026-07-10 quant-ph hep-lathep-ph

Hardware-efficient quantum simulation of intense-field QED

classification quant-ph hep-lathep-ph
keywords intense-field QEDFurry picturetrapped ionshybrid analog-digital simulationnonlinear Breit-WheelerVolkov stateszero-noise extrapolationspin-phonon gates
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

Strong laser fields turn quantum electrodynamics into a real-time nonperturbative problem: fermions are dressed by the background while photons remain dynamical. This paper proposes a concrete way to put that dynamics on trapped-ion hardware. Photons live in collective phonon modes and Volkov-dressed fermions live in ion spins; Clifford circuits compress the nonlocal strings that arise from the Jordan–Wigner encoding so that native spin–phonon gates can finish the interaction. For the nonlinear Breit–Wheeler process the gate count grows only polynomially with the number of momentum modes. A single-mode benchmark matches exact evolution once Trotter error is controlled, and experimentally realistic phonon heating and dephasing produce visible errors that zero-noise extrapolation largely removes from both photon survival and pair-production probabilities. The result is a hardware-efficient route to particle-production dynamics that sit beyond ordinary perturbation theory or static-field approximations.

Core claim

The authors show that intense-field QED in 3+1 dimensions, formulated in the Furry picture, can be mapped onto a trapped-ion hybrid analog–digital circuit whose resources scale polynomially with the number of retained momentum modes, and that the resulting circuit, when restricted to a single resonant mode of nonlinear Breit–Wheeler pair production, reproduces exact dynamics with controlled Trotter error and recovers the target observables under realistic noise via zero-noise extrapolation.

What carries the argument

Hybrid analog–digital compilation: Clifford circuits (nearest-neighbor CNOTs, Hadamard and phase gates) recursively compress nonlocal Jordan–Wigner strings into local operators so that native spin–phonon gates can implement the residual boson–fermion couplings of the Furry-picture interaction Hamiltonian.

Load-bearing premise

That the single-mode resonant truncation and the experimentally quoted noise rates continue to keep errors under control once the circuit depth reaches the full multimode polynomial scaling required for a genuine 3+1-dimensional lattice.

What would settle it

Run the hybrid circuit for the single-mode Hamiltonian of Eq. (14) on a trapped-ion device with the stated heating and dephasing rates; if zero-noise extrapolation fails to restore both the photon-survival probability and the pair-production signal to within the reported error bars of the exact evolution, the central feasibility claim is falsified.

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

If this is right

  • Nonlinear Breit–Wheeler pair production becomes accessible on near-term trapped-ion hardware without truncating photon occupation into qubits.
  • The same encoding immediately extends to nonlinear Compton scattering and to modest multimode truncations that include higher harmonics.
  • Resource counts remain polynomial (O(N_p^{3}) CNOTs and O(N_p^{2}) analog gates per Trotter step), so larger momentum lattices stay in principle reachable.
  • Zero-noise extrapolation of the dominant phonon-heating and dephasing channels recovers both survival and production observables under the noise model used.

Where Pith is reading between the lines

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

  • If the hybrid compilation generalizes cleanly, other strong-field processes that mix dressed fermions with dynamical bosons (e.g., multiphoton emission in magnetar magnetospheres) become natural targets for the same platform.
  • The separation of Clifford digital overhead from analog spin–phonon dynamics suggests that further gains may come from optimizing only the analog layer, a route not available to fully digital encodings.
  • Successful multimode runs would supply real-time correlation functions that classical lattice methods currently cannot reach in the nonperturbative regime.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 4 minor

Summary. The manuscript proposes a hybrid analog–digital trapped-ion protocol for real-time simulation of intense-field QED in 3+1 dimensions in the Furry picture. Photon modes are encoded in collective phonons and Volkov-dressed fermion modes in ion spins via Jordan–Wigner; Clifford circuits compress the nonlocal strings so that native spin-phonon gates realize the local boson-fermion couplings. Resource counts are given as O(N_p^{3}) CNOTs and O(N_p^{2}) analog gates per Trotter step (with overall depth O(N_p^{5} t^{2} C/ε) for first-order Trotter). The construction is benchmarked on a single-mode (and l_max=2) resonant truncation of nonlinear Breit–Wheeler pair production: the interaction Hamiltonian is reduced to Eq. (14), exact diagonalization is compared with noiseless Trotter circuits, and a Lindblad model of phonon heating plus spin/phonon dephasing is mitigated by third-order polynomial zero-noise extrapolation, recovering both photon-survival and pair-production signals within fitting uncertainty.

Significance. If the protocol and its error-mitigation performance hold under the stated assumptions, the work supplies a concrete, hardware-native route from the Furry-picture IFQED Hamiltonian to near-term trapped-ion operations. Encoding dynamical photons directly in phonons avoids the usual bosonic truncation overhead of purely digital encodings, while the Clifford compression of Jordan–Wigner strings keeps the digital layer efficient. The single-mode benchmark is carefully cross-checked against exact dynamics and uses experimentally quoted noise rates rather than fitted parameters, giving a credible proof-of-principle that error-mitigated hybrid circuits can access nonperturbative pair-production signals. The result is therefore a useful bridge between high-intensity laser QED and quantum-simulation hardware, even though genuine multimode 3+1D scaling remains untested.

major comments (2)
  1. The central claim of a “hardware-efficient route to intense-field particle-production dynamics” rests on the asymptotic gate counts given after Eq. (9) (O(N_p^{3}) CNOTs and O(N_p^{2}) analog gates per Trotter step, overall depth O(N_p^{5} t^{2} C/ε)). All numerical evidence, however, is confined to the single-mode (or l_max=2) resonant truncation of Eq. (14) and the corresponding circuits of Fig. 1(d). No noisy multimode circuit is simulated, so it is not shown that the quoted phonon-heating and dephasing rates, together with third-order polynomial ZNE, continue to control errors once circuit depth reaches the full multimode scaling. The multimode extrapolation is therefore an unproven premise rather than a demonstrated result; either a modest multimode noisy benchmark or a clear statement that the claim is limited to the single-mode setting is needed.
  2. Energy-momentum conservation is enforced by retaining only resonant (Q_0=0) terms, so that the phase factor e^{i Q_0 t} in Eq. (3) becomes unity. Off-resonant contributions and the continuous time dependence they generate are discarded without a quantitative estimate of the truncation error for the chosen laser and photon parameters (ξ=1, ω=1.55 eV, ω′=1 TeV). Because the subsequent Trotter and noise analyses inherit this truncation, a bound or numerical check on the size of the neglected terms would strengthen the claim that the reduced Hamiltonian of Eq. (14) faithfully represents the target IFQED process.
minor comments (4)
  1. Fig. 1(e) caption states that error bars are enlarged by a factor of fifteen for visibility; the actual (unenlarged) fitting uncertainties should also be reported so that the quality of the ZNE recovery can be judged quantitatively.
  2. The coefficient values in Eq. (13) and the multimode extension (C9) are given to two decimal places without an indication of numerical precision or of the Bessel-function truncation used to obtain them; a brief statement of the computational procedure would aid reproducibility.
  3. Notation for the Volkov modes switches between ψ^((±)s)_p, U^s_p / V^s_p and the quadruplet B̂^m_p; a short glossary or consistent usage would improve readability of the Supplemental Material.
  4. The gate-duration scalings (100 µs CNOT, 10 µs single-qubit, τ_gate = 10×(t/10^7 eV^{-1}) µs) are stated without reference to a specific ion species or trap frequency beyond the later 40Ca+ example; clarifying the assumed platform would help experimental groups assess feasibility.

Circularity Check

0 steps flagged

No circularity: single-mode benchmarks compare the hybrid circuit to exact diagonalization of the same truncated Hamiltonian; noise rates and ZNE are external, not fitted to the target signals.

full rationale

The paper's central claims are a hardware mapping (phonons for photons, Jordan–Wigner spins for Volkov modes, Clifford compression of strings, native spin-phonon gates) and a resource-scaling count, followed by a controlled numerical benchmark of the single-mode nonlinear Breit–Wheeler Hamiltonian (Eq. 14) against its own exact time evolution. The noiseless Trotterized circuit is shown to track the exact unitary of that same Hamiltonian; the noisy circuit uses literature values for heating and dephasing (not fitted to the survival or pair-production curves); and ZNE is a standard unitary-folding extrapolation whose zero-noise limit is again compared to the same exact dynamics. Nothing is predicted from a parameter that was itself fitted to the quantity being predicted, no uniqueness theorem is imported from overlapping authors to force the construction, and the Volkov/Furry starting point is standard textbook material rather than a self-citation that closes a loop. The multimode scaling argument is asymptotic and untested under noise, but that is an unproven extrapolation, not a circular reduction of a claimed result to its inputs. Score 0 is therefore the correct assessment.

Axiom & Free-Parameter Ledger

5 free parameters · 5 axioms · 0 invented entities

The central claim rests on standard QED and ion-trap machinery plus a handful of numerical choices that fix the single-mode benchmark and the noise model. No new physical entities are postulated; the free parameters are experimental or lattice cut-offs that control the size of the demonstrated circuit rather than the form of the Hamiltonian.

free parameters (5)
  • laser intensity parameter ξ and photon energies = ξ=1, ω=1.55 eV, ω′=1 TeV
    ξ=1, ω=1.55 eV, ω′=1 TeV (hence χ≈11.87) are chosen by hand to place the process in a nonperturbative regime; they fix the numerical coefficients in Eq. (13).
  • momentum lattice spacing Δp = 0.03 MeV
    Δp=0.03 MeV (L=4.12×10^{-11} m) is selected so that MeV-scale transverse momenta remain distinguishable while the TeV energy scale is retained; it sets the overall prefactor of Hint.
  • phonon heating rate and coherence times = 3 s^{-1}, 3 ms, 10 ms
    Γ_heat=3 s^{-1}, phonon T2=3 ms, spin T2=10 ms are taken from conservative experimental reports and directly determine the size of the noisy deviations that ZNE must correct.
  • gate durations and Trotter step schedule = 100 µs / 10 µs / 1-or-2 steps
    CNOT 100 µs, single-qubit 10 µs, analog gate duration scaled as 10×(t/10^7 eV^{-1}) µs, and the switch from 1 to 2 Trotter steps at t=4×10^7 eV^{-1} are free numerical choices that control both circuit depth and Trotter error.
  • ZNE folding factors and polynomial order = Ns=1..9, order 3
    Ns∈{1,3,5,7,9} and a third-order polynomial fit are chosen to extrapolate to zero noise; the order is not derived from a noise model but selected for empirical recovery.
axioms (5)
  • domain assumption Furry-picture Volkov states resum the classical background field to all orders, leaving only the quantized photon field dynamical.
    Invoked throughout Appendix A and the main-text IFQED framework section; standard in intense-field QED literature.
  • standard math Jordan-Wigner transformation correctly encodes the anticommutation relations of the Volkov fermion modes on a 1-D spin chain.
    Eq. (4) and the subsequent string compression; textbook encoding.
  • domain assumption The native ion-laser interaction in the Lamb-Dicke and rotating-wave approximations yields the spin-phonon Hamiltonian of Eq. (10).
    Appendix B; taken from the trapped-ion literature (Davoudi et al., Monroe et al.).
  • domain assumption Phonon heating and dephasing are adequately described by the Lindblad operators L[a], L[a†] and L[a†a] (and spin dephasing by L[σz]).
    Supplemental Appendix D; standard open-system model for ion traps.
  • standard math First-order Trotter-Suzuki decomposition converges with the stated O(t² N_p² C/ε) step count.
    Used for the resource scaling after Eq. (9).

pith-pipeline@v1.1.0-grok45 · 27180 in / 3579 out tokens · 40205 ms · 2026-07-14T15:05:32.457968+00:00 · methodology

0 comments
read the original abstract

Strong electromagnetic backgrounds make quantum electrodynamics a real-time nonperturbative problem involving dressed fermions and dynamical photons. We propose a trapped-ion protocol for simulating intense-field QED in $3+1$ dimensions in the Furry picture. The construction encodes photon modes in collective phonons and Volkov-dressed fermion modes in ion spins, combining native spin-phonon couplings with Clifford circuits that compress nonlocal Jordan--Wigner strings. For nonlinear Breit--Wheeler pair production, the protocol has polynomial resource scaling and is benchmarked against exact single-mode dynamics with controlled Trotter errors. With experimentally motivated phonon heating and dephasing, zero-noise extrapolation substantially reduces deviations in photon-survival and pair-production signals. These results provide a hardware-efficient route to intense-field particle-production dynamics beyond perturbative or static-field descriptions.

Figures

Figures reproduced from arXiv: 2607.09844 by Bin Xu, Jing Shu, Yiheng Lin, Ying-Ying Li, Yuxiang Huang, Zhongtian Dong, Zhuoyi Li.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Time evolution of channel probabilities for an incident left-polarized photon. The probabilities for [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: The survival channel [PITH_FULL_IMAGE:figures/full_fig_p017_3.png] view at source ↗

discussion (0)

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

Reference graph

Works this paper leans on

36 extracted references · 11 linked inside Pith

  1. [1]

    Fedotov, A

    A. Fedotov, A. Ilderton, F. Karbstein, B. King, D. Seipt, H. Taya, and G. Torgrimsson, Advances in QED with intense background fields, Phys. Rept.1010, 1 (2023), arXiv:2203.00019 [hep-ph]

  2. [2]

    Di Piazza, C

    A. Di Piazza, C. Muller, K. Z. Hatsagortsyan, and C. H. Keitel, Extremely high-intensity laser interactions with fundamental quantum systems, Rev. Mod. Phys.84, 1177 (2012), arXiv:1111.3886 [hep-ph]

  3. [3]

    V. M. Kaspi and A. Beloborodov, Magnetars, Ann. Rev. Astron. Astrophys.55, 261 (2017), arXiv:1703.00068 [astro-ph.HE]

  4. [4]

    M. C. Ba˜ nuls, R. Blatt, J. Catani, A. Celi, J. I. Cirac, M. Dalmonte, L. Fallani, K. Jansen, M. Lewenstein, S. Montangero, C. A. Muschik, B. Reznik, E. Rico, L. Tagliacozzo, K. Van Acoleyen, F. Verstraete, U.-J. Wiese, M. Wingate, J. Zakrzewski, and P. Zoller, Simu- lating lattice gauge theories within quantum technologies, Eur. Phys. J. D74, 165 (2020)...

  5. [5]

    C. W. Baueret al., Quantum Simulation for High- Energy Physics, PRX Quantum4, 027001 (2023), arXiv:2204.03381 [quant-ph]

  6. [6]

    Funcke, T

    L. Funcke, T. Hartung, K. Jansen, and S. K¨ uhn, Review on quantum computing for lattice field theory (2023), arXiv:2302.00467 [hep-lat]

  7. [7]

    Di Meglioet al., Quantum Computing for High-Energy Physics: State of the Art and Challenges

    A. Di Meglioet al., Quantum Computing for High-Energy Physics: State of the Art and Challenges. Summary of the QC4HEP Working Group (2023), arXiv:2307.03236 [quant-ph]

  8. [8]

    Y. Fang, C. Gao, Y.-Y. Li, J. Shu, Y. Wu, H. Xing, B. Xu, L. Xu, and C. Zhou, Quantum frontiers in high energy physics, Sci. China Phys. Mech. Astron.68, 260301 (2025), arXiv:2411.11294 [hep-ph]

  9. [9]

    Hidalgo and P

    L. Hidalgo and P. Draper, Quantum simulations for strong-field QED, Phys. Rev. D109, 076004 (2024), arXiv:2311.18209 [hep-ph]

  10. [10]

    Draper, L

    P. Draper, L. Hidalgo, and A. Ilderton, Hamiltonian truncation and quantum simulation of strong-field QED beyond tree level, Phys. Rev. D113, 056010 (2026), arXiv:2509.15495 [hep-ph]

  11. [11]

    Davoudi, N

    Z. Davoudi, N. M. Linke, and G. Pagano, Toward simu- lating quantum field theories with controlled phonon-ion dynamics: A hybrid analog-digital approach, Phys. Rev. Res.3, 043072 (2021), arXiv:2104.09346 [quant-ph]

  12. [12]

    W. H. Furry, On Bound States and Scattering in Positron Theory, Phys. Rev.81, 115 (1951)

  13. [13]

    Jordan and E

    P. Jordan and E. P. Wigner, About the Pauli exclusion principle, Z. Phys.47, 631 (1928)

  14. [14]

    Davoudi, N

    Z. Davoudi, N. M. Linke, and G. Pagano, Toward simu- lating quantum field theories with controlled phonon-ion dynamics: A hybrid analog-digital approach, Phys. Rev. Res.3, 043072 (2021)

  15. [15]

    Monroe, W

    C. Monroe, W. C. Campbell, L.-M. Duan, Z.-X. Gong, A. V. Gorshkov, P. W. Hess, R. Islam, K. Kim, N. M. Linke, G. Pagano,et al., Programmable quantum sim- ulations of spin systems with trapped ions, Reviews of Modern Physics93, 025001 (2021)

  16. [16]

    Suzuki, Generalized Trotter’s Formula and System- atic Approximants of Exponential Operators and Inner Derivations with Applications to Many Body Problems, Commun

    M. Suzuki, Generalized Trotter’s Formula and System- atic Approximants of Exponential Operators and Inner Derivations with Applications to Many Body Problems, Commun. Math. Phys.51, 183 (1976)

  17. [17]

    C. R. Clark, H. N. Tinkey, B. C. Sawyer, A. M. Meier, K. A. Burkhardt, C. M. Seck, C. M. Shappert, N. D. Guise, C. E. Volin, S. D. Fallek, H. T. Hayden, W. G. Rellergert, and K. R. Brown, High-fidelity bell-state preparation with 40ca+ optical qubits, Phys. Rev. Lett.127, 130505 (2021)

  18. [18]

    X. Zhao, J. Bian, Y. Li, Y. Li, M. Zhang, and Y. Lin, High-fidelity two-qubit quantum logic gates in a trapped- ion chain using axial motional modes, Chinese Physics Letters42, 110601 (2025)

  19. [19]

    C. J. Ballance, T. P. Harty, N. M. Linke, M. A. Sepiol, and D. M. Lucas, High-fidelity quantum logic gates us- ing trapped-ion hyperfine qubits, Phys. Rev. Lett.117, 060504 (2016)

  20. [20]

    J. P. Gaebler, T. R. Tan, Y. Lin, Y. Wan, R. Bowler, A. C. Keith, S. Glancy, K. Coakley, E. Knill, D. Leibfried, and D. J. Wineland, High-fidelity universal gate set for 9Be + ion qubits, Phys. Rev. Lett.117, 060505 (2016)

  21. [21]

    C. D. Bruzewicz, J. Chiaverini, R. McConnell, and J. M. Sage, Trapped-ion quantum computing: Progress and challenges, Applied physics reviews6(2019)

  22. [22]

    S. Endo, S. C. Benjamin, and Y. Li, Practical quantum error mitigation for near-future applications, Phys. Rev. X8, 031027 (2018)

  23. [23]

    Temme, S

    K. Temme, S. Bravyi, and J. M. Gambetta, Error mitiga- tion for short-depth quantum circuits, Phys. Rev. Lett. 119, 180509 (2017)

  24. [24]

    Seipt, Volkov States and Non-linear Compton Scatter- ing in Short and Intense Laser Pulses (Verlag Deutsches Elektronen-Synchrotron, Hamburg, 2017) pp

    D. Seipt, Volkov States and Non-linear Compton Scatter- ing in Short and Intense Laser Pulses (Verlag Deutsches Elektronen-Synchrotron, Hamburg, 2017) pp. 24–43

  25. [25]

    Yakaboylu, Volkov wave function: its orthonormality and completeness (2015), arXiv:1505.02801 [quant-ph]

    E. Yakaboylu, Volkov wave function: its orthonormality and completeness (2015), arXiv:1505.02801 [quant-ph]

  26. [26]

    Q. A. Turchette, C. J. Myatt, B. E. King, C. A. Sack- ett, D. Kielpinski, W. M. Itano, C. Monroe, and D. J. Wineland, Decoherence and decay of motional quantum states of a trapped atom coupled to engineered reservoirs, Phys. Rev. A62, 053807 (2000)

  27. [27]

    Breuer and F

    H.-P. Breuer and F. Petruccione,The theory of open quantum systems(OUP Oxford, 2002)

  28. [28]

    M. A. Nielsen and I. L. Chuang,Quantum computation and quantum information(Cambridge university press, 2010). 7 Supplemental Material Appendix A: Intense-field QED in the Furry Picture

  29. [29]

    Furry picture

    Derivation of Volkov States In intense field QED (IFQED), the laser field is usually described as a background field. This back-ground field has to be treated to all orders, which can be achieved by going to the Furry interaction picture in which the background field is treated as part of the unperturbed system. When an electron interacts with a high-inte...

  30. [30]

    Quantization of Volkov state In IFQED, the free fermion field operators are linear combinations of the Volkov wave functions with the coefficient operators denoting electron annihilationb s p and positron creationd s† p : ψ(x) = ∫ d3p (2π)3 1√ 2Ep 2∑ s=1 ( bs pψ(+)s p (x) +ds† pψ(−)s p (x) ) (A13) ¯ψ(x) = ∫ d3p′ (2π)3 1√ 2Ep′ 2∑ s′=1 ( bs′† p′¯ψ(+)s′ p′ (...

  31. [31]

    The interaction Hamiltonian Unlike the electron case, photons free streaming in the strong laser due to its long mean free path. Thus the field operator for free photon is given by A(x) = ∫ d3k′ (2π)3 1√2Ek′ 2∑ r=1 ( ar k′ϵr k′e−ik′·x+ar† k′ϵr∗ k′eik′·x ) (A17) wherea r k′are the bosonic annihilation operators that follow the commutation rules [ar k′,ar′†...

  32. [32]

    and making the rotating- wave approximation (RWA) to neglect high frequency oscillating terms, we get the interaction Hamiltonian H′ ion-laser = 1 2 ℏΩjσ+ j eiηk,j(ar ke−iωkt+ar† k eiωkt)e−i(ωL j−ω0)t+iϕj + h.c.,(B3) where ωk and ω0, correspond respectively to the vibrational mode frequency, and the energy difference between the two internal levels of the...

  33. [33]

    We have simulated the situation of etaining only the l = 1 harmonic contribution

    truncatelto 2 As mentioned, the effect of the background field enters through the exchange of an effective harmonic number l. We have simulated the situation of etaining only the l = 1 harmonic contribution. To show the scalability of our framework, we should consider more background virtual photons participating in the interaction which will produce more...

  34. [34]

    For the case where l is truncated to 2, there are at least two set momentum modes, that is, we need to extend the spin chain to 8 ions

    Quantum circuit construction We order the fermionic modes sequentially along a one-dimensional spin chain, with the four internal states associated with each momentum mode occupying four consecutive qubits. For the case where l is truncated to 2, there are at least two set momentum modes, that is, we need to extend the spin chain to 8 ions. And the fermio...

  35. [35]

    In the present work, we focus on two main classes of errors

    Noise Model In trapped-ion quantum platforms, quantum coherence is affected by several distinct noise channels arising from environmental coupling and technical imperfections. In the present work, we focus on two main classes of errors. The first involves energy-exchange processes, including the anomalous heating of phonon modes caused by fluctuating ambi...

  36. [36]

    The Fock space of the phonons is truncated to five, which is sufficient to reduce truncation errors in the subspace encoding the photon occupation states under consideration

    Extrapolation to the zero noise limit To simulate the noisy dynamics of the quantum circuit shown in Fig.1(d), we represent the global state by a density matrix ρ∈HS1⊗HS2⊗HPs⊗HPa, whereHSi denotes the i-th qubit andHPs,a represent the storage and auxiliary phonon modes, respectively. The Fock space of the phonons is truncated to five, which is sufficient ...