{"id":"a5bd0239-308b-4500-a0b7-ee10d60b3da5","arxiv_id":"1908.08222","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A GRAPE-optimized single gate on a simulated transmon propagates and spectroscopically resolves a frozen two-neutron spin Hamiltonian under realistic noise.","lead":"This paper proposes putting the spin part of the force between two neutrons into the four energy levels of a single superconducting chip, then using one specially designed microwave pulse to evolve it in time. The authors' computer simulation of the noisy chip shows the spin states oscillating for several cycles, and the oscillation frequencies reveal the energy levels of the encoded nuclear interaction.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The single fixed-r gate cannot implement the r-dependent spin propagator required in Eq. (4); full-dynamics and classical-quantum co-processing claims are unsupported, though the frozen-spin simulation itself is coherent.","rationale":"I agree with the reader that the decisive weakness is the frozen-separation scope. My pass adds a sharper technical reason: even the proposed classical-quantum co-processing extension is not immediate, because the spin propagator U_dt(r) is itself r-dependent. Applying one fixed gate U(r0) to a spin register does not reproduce the r-resolved evolution in Eq. (4) unless the spatial wavefunction is a delta function at r0. The frozen-spin part of the paper is, as far as I can see, internally consistent; the optimized-pulse fidelity and the V_SD^3 simulation are not fully documented, but those are data/reproducibility weaknesses rather than demonstrated mathematical errors. Since the reader already judged the paper CONDITIONAL on essentially this scope issue, my critique does not change the verdict; it strengthens the reason why the condition is needed. A concrete two-bin fidelity check would settle whether the single-gate approach can be extended to a full H_LO Trotter step, and I expect it would fail except at artificially close separations where V_SD is approximately constant.","tokens_in":14072,"tokens_out":10850,"duration_ms":116045,"concrete_test":"Use a two-bin position superposition, |Psi0> = (|r1> + |r2>)/sqrt(2) (x) |down up>, with r1 and r2 separated by roughly 1 fm so V_SD(r1) differs significantly from V_SD(r2), using the expressions in Appendix A. Compute the exact one-Trotter-step state exp[-i(H_SI + V_SD) dt]|Psi0> and compare it with the proposed scheme: apply the single fixed transmon gate U(r0 = 3.5 fm) to the spin state, then classically propagate positions according to T + V_SI. If the trace fidelity between the two states is not close to 1, the single-gate protocol cannot implement Eq. (4) and the full-dynamics / co-processing claims must be restricted to the frozen-separation problem.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The frozen-separation spin simulation is internally coherent: the eigenvalue algebra in Appendix A checks out, the V_SD^3 extraction in Eqs. (13)-(14) is algebraically consistent, and the Lindblad noise model is reasonable. The load-bearing gap is the bridge from that simulation to the paper's 'nuclear dynamics' claim. In the Trotterized full propagator, Eq. (4) requires applying U_dt(r) = exp[-i V_SD(r) dt / hbar] to the spin part for each spatial point r independently, followed by spatial propagation with T + V_SI. V_SD(r) is strongly r-dependent through the one-pion-exchange Yukawa and tensor factors in Appendix A. A single transmon pulse implementing one fixed 4x4 unitary U_dt(r0) cannot implement this r-dependent family of unitaries on a wavefunction that is a superposition over r. The classical-quantum co-processing idea in Sec. VI provides no mechanism, such as a position-controlled gate, for applying different U_dt(r) to different spatial components. Classically tabulating U(r) and running a separate quantum gate per r bin would discard spin-position entanglement and is not equivalent to Eq. (4). The paper's own limitation statement in Sec. VI confines the demonstration to fixed r, but the title and abstract and the proposed co-processing pathway go beyond that. Thus the central claim should be read as a proof-of-principle for spin dynamics at fixed separation, not as an implementation of general two-nucleon H_LO dynamics.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a single-gate, single-qudit implementation of the short-time propagator exp(-i V_SD(r0)Δt/ħ) for the spin-dependent part of the two-neutron interaction at leading order of chiral EFT. The authors fix the internuclear separation r=3.5 fm and encode the four neutron spin states in the lowest four transmon levels. A GRAPE-optimized pulse is computed in QuTIP, and repeated application of the pulse is modeled with a Lindblad master equation including T1=30 μs and Tφ=50 μs. From the time-dependent occupation probabilities they extract all pairwise eigenenergy differences of V_SD (Table I), and by a second propagation driven by V_SD^3 they infer absolute eigenvalues (Table II). The appendices provide the analytic spin-eigenvalue decomposition and the Fourier-peak fitting procedure.","tokens_in":14433,"tokens_out":10377,"duration_ms":103156,"significance":"If read as a proof-of-principle for encoding the fixed-r spin part of the chiral-EFT NN interaction in a multi-level transmon and for noise-robust spectroscopy of the resulting four-level system, the work is a solid and useful contribution. The numerical experiments are internally consistent: the GRAPE optimization is run to an infidelity threshold below 10^-4; the Lindblad simulations show several oscillation cycles; and the extracted peak locations agree with the analytic eigenvalues within the quoted uncertainties. The main weakness is that the title and abstract advertise 'nuclear dynamics' of two interacting neutrons, whereas the simulated object is only the 4x4 spin propagator at one fixed separation, with the spatial part of H_LO not implemented. With an appropriate reframing, the paper would be a meaningful step toward qudit-based simulation of spin-isospin nuclear interactions.","major_comments":[{"comment":"The central simulation is not the two-nucleon Hamiltonian of Eq. (1). Equation (4) shows that the full Trotterized propagator requires applying exp[-i V_SD(x)Δt/ħ] separately for each spatial point x while the spin state is entangled with the spatial wavefunction. The implemented pulse realizes one fixed unitary (Eq. (9)) for a single value r0=3.5 fm; a single fixed pulse cannot implement the x-dependent family of unitaries required by Eq. (4). The classical-quantum co-processing paragraph in Sec. VI is only a 'possibility' and does not provide a concrete mechanism (such as a position-controlled gate or a sum over r-bins that preserves spin-position entanglement) that would realize Eq. (4). The title, abstract, and several introductory statements therefore overstate the scope of the demonstration. Please either reframe the claims to the frozen-spin spin-propagator proof-of-principle, or provide and validate a concrete mechanism for the full H_LO dynamics.","section":"Sec. II, Eq. (4); Sec. IV, Eq. (9); Sec. VI"},{"comment":"The extraction of absolute eigenvalues uses the trace of V_SD, but the text says that because this matrix is traceless, a constant diagonal matrix is added to induce a nonzero trace. Adding c times the identity shifts every eigenvalue by c and changes no observable probabilities. The manuscript does not specify the value of c or how the reported eigenvalues in Table II are obtained from the shifted ones. Without this relation, the trace condition is underdetermined and the absolute calibration is not established. Please state the shift explicitly and show how the eigenvalues in Table II follow from it.","section":"Sec. V, Eqs. (13)-(14) and the trace paragraph"},{"comment":"The relationship between the plotted nuclear evolution time (in MeV^-1), the Trotter step Δt=0.30 MeV^-1, and the physical duration of each 100 ns GRAPE pulse is not given. This conversion determines how many pulses are applied and hence the total physical time over which T1 and Tφ act; without it the reader cannot check whether the observed attenuation is consistent with the stated coherence times. Please state the mapping between device-time and nuclear-time units and the total physical duration represented by the plots.","section":"Sec. V, Eq. (10) and Fig. 4"}],"minor_comments":[{"comment":"There are typos: 'seperation' in Sec. II and 'Linblad' in Sec. V should be 'separation' and 'Lindblad', respectively.","section":"Sec. II and Sec. V"},{"comment":"The phrase 'trace trace' is duplicated in the paragraph on absolute eigenvalue extraction; please correct.","section":"Sec. V"},{"comment":"The numeric labels '45 95 144 156 324' in the figure are unexplained; please clarify what these numbers denote.","section":"Fig. 4 caption"},{"comment":"References [30] and [49] are the same paper (Paik et al.), as are [43] and [50] (Rigetti et al.); please consolidate or cite once.","section":"References"},{"comment":"The axis label 'Detuned Frequency from QPU Ground State (GHz)' is unclear for a Fourier transform of the drive amplitudes; please clarify what is displayed.","section":"Fig. 3(b)"}],"recommendation":"major_revision","confidential_remarks":"To the editor: the underlying frozen-spin simulation appears sound and could be a legitimate proof-of-principle once the claims are rescaled. The main issue is the gap between the advertised 'nuclear dynamics' and the fixed-r spin propagator actually simulated; this is fixable by rewriting the title and abstract and adding the requested clarifications. The trace/calibration issue in the eigenvalue section should be resolved before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The useful core here is a single 4-level transmon gate, found by GRAPE, that implements the fixed-separation spin propagator of the LO chiral-EFT two-neutron interaction. The Lindblad simulation with realistic T1 and Tphi shows several oscillation cycles before decoherence, and the Fourier peaks reproduce the analytic eigenenergy differences within the quoted errors. The V^3 trick to get absolute eigenvalues without phase estimation is neat, and the algebra in Eqs. (13)-(14) is consistent. At the level of what is actually simulated, the work is sound: the spin Hamiltonian in Appendix A is correct, the pulse optimization reaches 1e-4 infidelity, and the noise model is reasonable. The paper deserves credit for a careful, internally coherent emulation.\n\nThe main weakness is a scope mismatch. Equation (4) requires applying exp[-i V_SD(x) dt] independently at each spatial point x on a wavefunction that is a superposition over r. A single fixed-r pulse cannot implement that family of unitaries. The classical-quantum co-processing idea in Sec. VI has no mechanism, such as a position-controlled gate, to apply different U(r) to different spatial components. Tabulating U(r) per bin would discard spin-position entanglement, so it is not equivalent to the Trotterized H_LO propagation. The paper's own limitation statement in Sec. VI does confine the demonstration to fixed r, which is honest, but the title, abstract, and conclusion go beyond that. That overreach should be fixed by reframing the claims as a proof-of-principle for frozen-spin dynamics.\n\nOther soft spots are minor by comparison: there is no hardware data, the optimized pulse data and numerical LEC values are not given, so the numerics are not independently reproducible, and the time-correlated Gaussian noise model in Appendix B is plausible but unvalidated. None of this undermines the emulated claim.\n\nThis paper is for people working on qudit encodings for NISQ simulation and nuclear theorists interested in hardware-efficient simulation of spin Hamiltonians. It deserves a serious referee. I would recommend conditional acceptance with major revision on scope, or at least a clear editorial direction to align the title and abstract with the actual frozen-r result.","headline":"A coherent fixed-separation spin-dynamics emulation with a real eigenvalue-extraction trick, wrapped in a title and abstract that overstate the scope to full nuclear dynamics.","tokens_in":614,"tokens_out":838,"would_cite":false,"duration_ms":25193,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A single optimized microwave gate on a 3D transmon can simulate the spin dynamics of two interacting neutrons, and with realistic noise the surviving oscillation signal yields the interaction's complete energy spectrum without quantum…","keywords":["quantum simulation","chiral effective field theory","two-neutron spin dynamics","3D transmon","gradient ascent pulse engineering","Fourier spectroscopy","Lindblad master equation","near-term quantum hardware"],"falsifier":"On real hardware, prepare a 3D transmon in the $|\\downarrow\\uparrow\\rangle$ state, apply the 100 ns optimized pulse repeatedly, and take the discrete Fourier transform of the four occupation probabilities: the central claim fails if the three dominant peaks do not appear at the time-converted values 2.5254, 3.3951, and 5.9205 MeV under the stated $T_1 = 30\\,\\mu\\mathrm{s}$, $T_\\phi = 50\\,\\mu\\mathrm{s}$ noise, or if the signal decays before completing one full oscillation cycle. A direct check of the encoding step is quantum process tomography of the implemented gate against $\\exp(-i\\hat{V}_{\\mathrm{SD}}\\Delta t/\\hbar)$, which should show infidelity below the claimed $10^{-4}$ threshold.","tokens_in":13901,"feed_emoji":"⚛️","tokens_out":16446,"duration_ms":134380,"temperature":0.7,"pith_summary":"This paper claims that the real-time evolution of two interacting neutrons can be enacted on a single multi-level superconducting device, a 3D transmon, using one numerically optimized microwave pulse instead of a long cascade of elementary gates. The simulated object is the spin-dependent part of the leading-order chiral effective field theory potential, $\\hat{V}_{\\mathrm{SD}}$, evaluated at a fixed internuclear separation. Simulating the driven device with realistic relaxation and dephasing, the authors show that the four spin-state occupation probabilities keep oscillating for multiple cycles, and the discrete Fourier transform of that signal exhibits peaks at every distinct pairwise eigenenergy difference of $\\hat{V}_{\\mathrm{SD}}$. A second simulation driven by the cube of the interaction then recovers the absolute eigenvalues without quantum phase estimation. If hardware reproduces these master-equation results, the approach offers a noise-tolerant route to real-time nuclear dynamics on near-term quantum processors.","feed_headline":"One optimized gate reads a neutron pair's energy spectrum","feed_subtitle":"Under realistic transmon noise the spin signal survives long enough to recover the nuclear eigenenergies.","key_machinery":"The load-bearing object is a single dense multi-level gate: a 100 ns microwave drive, optimized by gradient ascent pulse engineering (GRAPE), whose time-ordered evolution under the transmon-plus-drive Hamiltonian reproduces the target unitary $\\exp(-i\\hat{V}_{\\mathrm{SD}}\\Delta t/\\hbar)$ to an infidelity below $10^{-4}$. The gate works by mapping the four uncoupled two-neutron spin states onto the lowest four Fock levels of the transmon and by driving the device's own transition frequencies, so no population leaks out of the computational manifold. On the analysis side, the argument is carried by the identity $|\\langle \\xi_i | \\Psi(t) \\rangle|^2 = \\sum_{j,k} c_j b^i_j c_k^* b^{i*}_k e^{-i\\Delta\\lambda_{jk} t/\\hbar}$, which turns the discrete Fourier transform of the measured occupation probabilities into a direct readout of every pairwise eigenenergy difference $\\Delta\\lambda_{jk}$ of $\\hat{V}_{\\mathrm{SD}}$. A generalized least-squares fit with a time-correlated Gaussian covariance model locates the spectral peaks between Fourier grid points, and a second propagation with $\\hat{V}_{\\mathrm{SD}}^3$ supplies the nonlinear pair $\\lambda_3 - \\lambda_0 = \\alpha$, $\\lambda_3^3 - \\lambda_0^3 = \\beta$ whose solution, combined with the Hamiltonian trace, fixes the absolute eigenvalues without quantum phase estimation.","core_discovery":"The paper establishes that the short-time spin propagator $\\exp(-i\\hat{V}_{\\mathrm{SD}}\\Delta t/\\hbar)$ of the neutron–neutron interaction at leading order of chiral EFT can be embedded as a single dense gate on the lowest four levels of a 3D transmon, with the four uncoupled spin states $|\\downarrow\\downarrow\\rangle$, $|\\downarrow\\uparrow\\rangle$, $|\\uparrow\\downarrow\\rangle$, $|\\uparrow\\uparrow\\rangle$ mapped onto the Fock states $|0\\rangle$ through $|3\\rangle$. The gate is a 100 ns drive pulse found by gradient ascent pulse engineering, optimized to an infidelity below $10^{-4}$, and it acts by driving the device's own transition frequencies so that none of the population leaves the computational manifold. Lindblad master-equation simulations at $T_1 = 30\\,\\mu\\mathrm{s}$ and $T_\\phi = 50\\,\\mu\\mathrm{s}$ show the occupation probabilities undergoing multiple full oscillation cycles, and the discrete Fourier transform of the $|\\downarrow\\uparrow\\rangle$ probability resolves all three distinct pairwise eigenenergy differences of $\\hat{V}_{\\mathrm{SD}}$ — 2.5254, 3.3951, and 5.9205 MeV — to within roughly 0.01–0.03 MeV. Propagating $\\hat{V}_{\\mathrm{SD}}^3$ in a second simulation, and solving the two-equation pair relating the largest eigenvalue difference $\\alpha$ and the largest cubed difference $\\beta$, yields the absolute eigenvalues ($-2.3(2)$, $0.9(6)$, $3.6(2)$ MeV versus exact $-2.329$, $1.066$, $3.592$ MeV) without quantum phase estimation.","pith_inferences":["My inference: the cubed-interaction trick is a general spectroscopic lever — propagating a low-degree polynomial of the interaction and matching extremal eigenvalue differences between the polynomial and the linear operator would let other few-level simulations convert eigenvalue differences into absolute eigenvalues without quantum phase estimation.","My inference: scaling beyond two neutrons is the open question the paper leaves unanswered; the GRAPE search and the single-device encoding both grow with the number of spin configurations, and whether the method remains practical for three or more nucleons requires a dedicated study.","My inference: a direct hardware test on a real transmon at the stated T1 and Tφ values, comparing measured Fourier linewidths with the master-equation prediction, would show whether Markovian relaxation and dephasing capture the dominant noise in this regime.","My inference: the frozen-separation approximation could be lifted in a concrete hybrid algorithm by discretizing the internuclear separation, running one optimized gate per grid point, and interleaving classical updates of the spatial wavefunction between time steps."],"forward_implications":["Energy-difference spectroscopy becomes available on near-term quantum hardware: the eigenenergies of a simulated nuclear interaction can be read from the Fourier transform of a short time signal, as long as the propagator survives a few oscillation periods.","The qudit encoding sidesteps the compilation overhead of breaking a nuclear propagator into many elementary one- and two-qubit gates, because the multi-level structure of the device itself carries the computation.","The scheme enables a classical–quantum co-processing protocol in which the quantum processor advances the spin wavefunction while a classical computer handles the spatial propagation, potentially bypassing the exponential growth of spin configurations that limits quantum Monte Carlo methods.","Because the pulse is derived from the device's own Hamiltonian, the single-gate recipe transfers to any short-time propagator expressible in the lowest levels of a multilevel superconducting circuit, reaching beyond nuclear physics."],"supporting_citations":[{"why":"Defines the chiral EFT expansion whose leading-order spin-dependent potential is the Hamiltonian the scheme simulates.","marker":"[28, 29]"},{"why":"Supplies the explicit coordinate-space forms of the interaction coefficients used to build the simulated potential and its analytic eigenvalues.","marker":"[38]"},{"why":"Introduces the 3D transmon architecture whose long coherence times make the single-gate approach practical.","marker":"[30]"},{"why":"Gives the transmon-plus-readout Hamiltonian assumed by the pulse optimization.","marker":"[44]"},{"why":"Demonstrates numerically optimized multi-level control in superconducting circuits and provides the drive Hamiltonian form the authors adopt.","marker":"[26]"},{"why":"The gradient-ascent pulse-engineering algorithm used to find the control pulse realizing the nuclear propagator.","marker":"[31]"},{"why":"The open-source quantum optics toolbox used for pulse optimization and for the Lindblad master-equation simulations with realistic noise.","marker":"[32]"}],"fun_headline_variants":["One dense gate on a transmon reads out nuclear eigenenergies","Neutron pair dynamics and spectrum from one optimized pulse","A single gate replaces quantum phase estimation for nuclear energies","Optimal control extracts neutron pair spectrum without QPE","Four-level transmon gate simulates neutron interaction and yields its spectrum"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The scheme rests on treating the two neutrons as frozen at one fixed internuclear separation, so that the single gate simulates only the spin propagator of the interaction: kinetic energy and the spin-independent potential are set aside for a classical co-processor that the paper neither implements nor tests. If that separation of spin and spatial dynamics is not valid for the physics being simulated, the gate does not implement the full nuclear dynamics promised by the paper's title.","fun_headline_variants_meta":{"raw":{"variants":["One dense gate on a transmon reads out nuclear eigenenergies","Neutron pair dynamics and spectrum from one optimized pulse","A single gate replaces quantum phase estimation for nuclear energies","Optimal control extracts neutron pair spectrum without QPE","Four-level transmon gate simulates neutron interaction and yields its spectrum"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000259,"raw_usage":{"total_tokens":1626,"prompt_tokens":1029,"completion_tokens":597,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":645,"completion_tokens_details":{"reasoning_tokens":514}},"tokens_in":645,"tokens_out":597,"duration_ms":5867,"temperature":1.0,"reasoning_tokens":514,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:46:20.777411+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"On real hardware, prepare a 3D transmon in the $|\\downarrow\\uparrow\\rangle$ state, apply the 100 ns optimized pulse repeatedly, and take the discrete Fourier transform of the four occupation probabilities: the central claim fails if the three dominant peaks do not appear at the time-converted values 2.5254, 3.3951, and 5.9205 MeV under the stated $T_1 = 30\\,\\mu\\mathrm{s}$, $T_\\phi = 50\\,\\mu\\mathrm{s}$ noise, or if the signal decays before completing one full oscillation cycle. A direct check of the encoding step is quantum process tomography of the implemented gate against $\\exp(-i\\hat{V}_{\\mathrm{SD}}\\Delta t/\\hbar)$, which should show infidelity below the claimed $10^{-4}$ threshold.","supporting_citations":[{"cited_title":"Gezerlis, I","cited_arxiv_id":null,"evidence_quote":"Supplies the explicit coordinate-space forms of the interaction coefficients used to build the simulated potential and its analytic eigenvalues."},{"cited_title":"Rigetti, J","cited_arxiv_id":null,"evidence_quote":"Gives the transmon-plus-readout Hamiltonian assumed by the pulse optimization."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The gradient-ascent pulse-engineering algorithm used to find the control pulse realizing the nuclear propagator."},{"cited_title":"Khaneja, T","cited_arxiv_id":null,"evidence_quote":"The open-source quantum optics toolbox used for pulse optimization and for the Lindblad master-equation simulations with realistic noise."}],"review_version":1}