{"id":"4d06416a-d2b5-4e99-ab8a-1a15c18dd950","arxiv_id":"2411.19070","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Near an engineered conical intersection in two trapped Rydberg ions, spontaneous decay damps the spin and phonon dynamics but not before several clear oscillations occur.","lead":"This paper simulates what happens to the motion and electronic state of two laser-excited Rydberg ions when the excited state decays, near a special crossing of energy levels called a conical intersection. It finds that the tell-tale oscillations survive long enough that the conical-intersection effect should be observable before the Rydberg state decays.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The model predicts that after Rydberg decay the y phonon is a free oscillator, so the reported ⟨Ny⟩→0 steady state cannot follow from the stated Lindblad equation; this points to an unmodeled phonon-loss channel or a numerical artifact.","rationale":"The reader's weakest assumption was that the single-channel Markovian decay γS=0.13 µs⁻¹ captures the relevant dissipation. That is a fair physical concern, but the most load-bearing problem I find is internal: the stated master equation cannot produce phonon relaxation of ⟨Ny⟩ to zero. Because the dissipator does not act on phonons and the post-decay spin sector is decoupled from the phonons, the y-phonon occupation is frozen once the Rydberg state has decayed. Since the coherent dynamics generate nonzero y-phonon occupation before decay, the long-time value must be positive. The paper's claim that ⟨Ny⟩=0 in the steady state therefore signals a numerical or implementation artifact, not a consequence of the model. This undermines the quantitative reliability of the dissipative phonon dynamics, though the central early-time visibility claim may still survive. The reader's conditional verdict remains appropriate, but the required revision is more specific than the reader's stated concern: the authors should either identify the extra phonon damping or correct the master-equation numerics. Because the final verdict is unchanged, I set verdict_should_be to UNCHANGED.","tokens_in":14771,"tokens_out":27816,"duration_ms":289216,"concrete_test":"Run the same master-equation simulation with Gx=0, Gy=2π×0.86 MHz, ωy=2π×1.6 MHz, γS=0.13 µs⁻¹, initial state |01⟩⊗|0_y⟩, and no phonon collapse operators. Track ⟨Ny⟩ after the Rydberg population has decayed. In the post-decay sector the Hamiltonian is a free oscillator, so ⟨Ny⟩ must become constant. If it instead decays to zero, there is an undocumented phonon-loss mechanism or a solver artifact. A positive, nonzero constant would confirm the reported ⟨Ny⟩→0 in Fig. 4b is inconsistent with the stated model.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Sec. 3.2 the paper states that although there is no phonon loss term, ⟨Ny⟩ decays and reaches ⟨Ny⟩=0 in the steady state, with the parity symmetry broken by Rydberg decay. This is internally inconsistent with the master equation (6)-(7). The dissipator (7) changes only the internal state |0⟩→|g⟩; it does not act on the phonon Fock state. The Hamiltonian (3) in the sector where one ion is in |g⟩ and the other in |1⟩ is just the free oscillator Hamiltonian, since Sx and Sz vanish on |g⟩. Therefore, once the Rydberg population has decayed, ⟨Ny⟩ is conserved. Moreover, the coherent dynamics generated by Gy(a†+a)Sx starting from |01,0⟩ (Eq. 4) populates y-phonon states with Ny>0 (see Figs. 2b and 3b); the no-jump y-phonon state associated with the |0⟩ component has conditional occupation |α|² tanh(|α|²)>0 for t>0. A quantum jump leaves that Fock state unchanged. Thus the steady-state value of ⟨Ny⟩ is a positive decay-time average, not zero. The observation ⟨Ny⟩=0 therefore indicates either an undocumented phonon damping in the simulation or a steady-state solver artifact. This does not by itself disprove the early-time oscillatory visibility claim, but it means the quantitative phonon relaxation results in Fig. 4b are not supported by the stated model.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies how the finite lifetime of a Rydberg state affects vibronic dynamics near a conical intersection in a system of two trapped Rydberg ions. The authors take the CI Hamiltonian from their earlier PRL, quantize the phonon modes, and add a Lindblad dissipator describing spontaneous decay from the Rydberg state |0> to a low-lying state |g>. Using mean-field equations and QuTiP master-equation simulations, they report damping of spin and phonon observables and claim that oscillatory spin-phonon dynamics remain visible within the Rydberg lifetime. The stated parameters are taken from prior work and atomic lifetime data without fitting to the target result.","tokens_in":15002,"tokens_out":8884,"duration_ms":82849,"significance":"The qualitative conclusion—that CI-induced spin-phonon oscillations survive for several microseconds before Rydberg decay destroys the collective spin sector—is interesting and experimentally relevant for the emerging trapped-Rydberg-ion platform. The paper has the merit of starting from a concrete Hamiltonian with parameters from previous experiments and independent atomic data, and the master-equation setup is clearly stated. If the numerical results are correct, the collective spin observables provide a plausible observable target. However, the reported phonon relaxation to zero in the absence of any phonon-loss term is internally inconsistent with the stated model, and the absence of convergence checks makes the quantitative claims unreliable. These issues directly affect the central 'relaxation dynamics' part of the paper.","major_comments":[{"comment":"The reported steady state ⟨Ny⟩=0 in Fig. 4b is inconsistent with the stated master equation. The dissipator in Eq. (7) acts only on the electronic states via σ^{g0} and leaves the phonon Fock state unchanged. Once the population has decayed to |g⟩, the Hamiltonian (3) in that sector reduces to the free oscillator ω_y a_y† a_y, since Sx=Sz=0 on |g⟩. Consequently the y-phonon number distribution in the |g⟩ sector is conserved, and the no-jump path starting from Eq. (4) populates y-phonon states with ⟨Ny⟩>0 (see Fig. 2b). The steady-state value of ⟨Ny⟩ is therefore a positive decay-time average, not zero. This means the phonon relaxation and steady-state statements in Sec. 3.2 and Appendix B (A=B=0) do not follow from Eqs. (6)-(7). The authors should either include an explicit phonon-loss mechanism or revise the phonon dynamics claims.","section":"Section 3.2, Eqs. (6)-(7)"},{"comment":"No Fock-basis truncation dimension or convergence test is reported for the QuTiP simulations. The stated phonon decay could be an artifact of an insufficient Fock cutoff, which removes population that would otherwise remain in high-|n⟩ states. Given the unexpected decay of ⟨Nx⟩ and ⟨Ny⟩ in Fig. 4b, the authors should specify the truncation used and demonstrate convergence with respect to the cutoff before the quantitative phonon relaxation results are accepted.","section":"Section 3.2, Fig. 4"},{"comment":"The statement 'Hence, (⟨Sz⟩+⟨Ny⟩) is a conserved quantity' does not follow from the parity symmetry P=Sz e^{iπNy}. A direct calculation gives [H, Sz+Ny] = -2i Gy(a_y†+a_y)Sy + Gy Sx(a_y†-a_y), which is generically non-zero. The parity P is indeed conserved, but the sum ⟨Sz⟩+⟨Ny⟩ is not. This invalidates the explanation of the small ⟨Ny⟩ in Fig. 2b and should be corrected.","section":"Section 2.2"},{"comment":"The mean-field fixed point A=B=0 is not the long-time limit of the mean-field equations. When the spin expectation values have decayed to zero, the equations reduce to ˙A=-iωx A and ˙B=-iωy B, whose solutions are undamped oscillations; A=B=0 is reached only if the initial phonon amplitudes vanish. Therefore the claimed steady state with vanishing expectation values is inconsistent with the initial coherent state used in the paper, and the agreement with the exact numerics is not established.","section":"Appendix B"}],"minor_comments":[{"comment":"There is a typo: 'Rybderg' should be 'Rydberg'.","section":"Section 2.1"},{"comment":"The decay rate γS = 0.13 μs⁻¹ corresponds to a lifetime of about 7.7 μs, while the text states 7.2 μs for the 50S state; please reconcile the numbers.","section":"Section 3.1"},{"comment":"The dotted black line representing e^{-γS t} is not identified in the legend; specify which curve it corresponds to.","section":"Figure 6"},{"comment":"The mean-field equations should explicitly state the factorization assumption ⟨AB⟩≈⟨A⟩⟨B⟩ used to close the equations, since this is essential for reproducing the derivation.","section":"Appendix B"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is closely based on the authors' earlier PRL model [30]; the incremental contribution here is the dissipative master-equation study. The central visibility claim is plausible and worth publishing after the phonon-relaxation inconsistency and convergence issues are resolved. The journal should verify that the numerical code or data are made available, since the current 'available from the authors upon reasonable request' statement is not independently checkable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this paper adds spontaneous decay to the trapped-Rydberg-ion conical intersection proposal and claims oscillatory vibronic dynamics survive a few microseconds, giving experimentalists a parameter window. That early-time visibility claim is plausible and useful. But the phonon relaxation results contain a real internal inconsistency: the stated master equation cannot produce ⟨Ny⟩→0 in the steady state. Without a phonon loss term, decay |0⟩→|g⟩ leaves the phonon Fock state untouched, and in the |g⟩ sector the y phonon is a free oscillator. The final ⟨Ny⟩ should be a positive decay-time average, not zero. This suggests either an unmodeled phonon damping channel in the simulation or a numerical artifact. The same logic applies to ⟨Nx⟩, though the paper is less explicit there. So the quantitative phonon relaxation claims in Fig. 4b and the surrounding text are not supported by the stated model.\n\nThe early-time oscillations themselves are probably fine—they occur within the first few microseconds, before the steady state matters—so the central visibility claim is credible. The setup is a direct extension of the authors’ PRL, adding quantized vibrations and a single-channel Lindblad decay. Parameters come from prior trap experiments and independent lifetime data, with no fitting to the target result. That makes it a legitimate feasibility case study, not a fundamental new effect.\n\nSoft spots beyond the inconsistency: Section 2.2 claims parity implies ⟨Sz⟩+⟨Ny⟩ is conserved, but the parity operator P = S_z e^{iπN_y} commutes with H; it does not make that sum conserved. The transition |↓,0⟩↔|↑,1⟩ changes the sum by 2, so that remark is simply wrong. Also, no convergence checks or code/data are provided; for a numerical case study that is a real issue. The simplifications—single-channel decay, and neglect of micromotion and motional heating—are acknowledged and bound the feasibility claim, though they are not fatal.\n\nThis paper is for experimentalists in trapped Rydberg ions and theorists estimating whether the CI signal is observable. It deserves a serious referee: the early-time visibility claim is worth checking, but the phonon relaxation analysis needs to be redone or heavily qualified. I would send it to review, expecting major revision.","headline":"Useful feasibility case study for observing conical-intersection dynamics in trapped Rydberg ions, but the phonon relaxation results are internally inconsistent with the stated master equation.","tokens_in":15629,"tokens_out":5026,"would_cite":false,"duration_ms":57169,"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":"This paper claims that spin–phonon oscillations encoding a conical intersection in two trapped Rydberg ions remain clearly visible within the microsecond Rydberg lifetime, so the conical effect can be observed despite dissipation.","keywords":["conical intersection","trapped Rydberg ions","spin-phonon dynamics","master equation","dissipative dynamics","Rydberg lifetime","vibronic dynamics","quantum simulation"],"falsifier":"A concrete falsifier would be a two-Sr+-ion experiment prepared in $|0,1\\rangle=|50S,50P\\rangle$ with the $x$-phonon mode in a coherent state: the model predicts that $\\langle xS_z\\rangle$ and $\\langle N_x\\rangle$ oscillate and damp on the ~7 $\\mu$s 50S lifetime. If the measured spin–phonon correlations disappear within about one microsecond, or if the $|0\\rangle$ population decays appreciably faster than $e^{-\\gamma_S t}$, then the single-channel decay model omits a dominant decoherence mechanism and the visibility claim fails.","tokens_in":14449,"feed_emoji":"⚛️","tokens_out":8018,"duration_ms":68301,"temperature":0.7,"pith_summary":"This paper asks whether the finite lifetime of Rydberg states prevents a trapped-ion simulator from showing the dynamics of a conical intersection, a point where two electronic energy surfaces meet and the usual separation between electronic and nuclear motion breaks down. Using a master equation for two trapped strontium ions, it finds that spontaneous decay of the nS Rydberg state damps both the collective spin and the phonon populations, yet the oscillatory spin–phonon response survives for several microseconds, which is within the Rydberg lifetime. The authors conclude that the conical-intersection dynamics can be observed before the Rydberg state decays, and that collective spin observables are the clearest place to look. This matters because trapped Rydberg ions have been proposed as a slow, controllable quantum simulator for conical-intersection physics, and the decay of the Rydberg states was an open threat to that programme.","feed_headline":"Conical intersection oscillations survive Rydberg decay","feed_subtitle":"A trapped-ion simulation shows the conical-intersection signal lasts several microseconds before the Rydberg state dies.","key_machinery":"The central object is the quantized vibronic Hamiltonian $H = \\omega_x(a_x^\\dagger a_x+\\tfrac12) + \\omega_y(a_y^\\dagger a_y+\\tfrac12) + G_x(a_x^\\dagger+a_x)S_z + G_y(a_y^\\dagger+a_y)S_x$, where $S_z$ and $S_x$ are collective spin operators of the two-ion Rydberg pair and $a_x$, $a_y$ are the two phonon modes; the perpendicular couplings $G_x$ and $G_y$ create the conical intersection in the collective spin sector. The dissipative part is a Lindblad master equation with a single decay channel $|0\\rangle\\to|g\\rangle$ at rate $\\gamma_S$. The coherent dynamics is constrained by the parity symmetry $P = S_z e^{i\\pi N_y}$, which leaves the Hamiltonian invariant and makes $\\langle S_z\\rangle+\\langle N_y\\rangle$ conserved; the decay breaks this symmetry, which is why the $y$-phonon number relaxes in the dissipative case. The argument proceeds by exact numerical integration of the master equation in a truncated Fock basis, with a mean-field steady-state analysis used to confirm that the long-time state has vanishing spin and phonon expectation values.","core_discovery":"The paper's central claim is that the conical-intersection signal in a trapped Rydberg ion pair is not erased by spontaneous decay of the Rydberg state forming it. Starting from the state $|0,1\\rangle$ with the $x$-phonon mode in a coherent state, the exact master-equation evolution shows $\\langle x\\rangle$ and the spin–phonon correlations $\\langle xS_z\\rangle$ and $\\langle yS_x\\rangle$ oscillating with microsecond periodicity while decaying on the ~7.2 $\\mu$s lifetime of the 50S state; the individual-ion populations oscillate only weakly. The decay $|0\\rangle\\to|g\\rangle$ at rate $\\gamma_S=0.13\\ \\mu\\text{s}^{-1}$ destroys the collective spin sector in which the conical intersection is defined, and because the spin and phonon modes are coupled, the phonon populations relax even though there is no direct phonon loss. The paper therefore frames the conical-intersection dynamics as a transient, observable phenomenon, with collective spin measurements as the recommended probe.","pith_inferences":["One could test how much of the geometric-phase suppression of tunneling survives the decay by extracting a tunneling probability from the dissipative wavepacket and comparing it with the coherent case; the paper reports the localization but does not quantify this visibility in the presence of decay.","Varying the principal quantum number $n$ would change both the spin–phonon couplings and the Rydberg lifetime together, so a natural extension is to search for an $n$ that maximizes the number of resolvable oscillation periods before the decay erases the signal.","The contrast with the superconducting-circuit experiment, where dephasing enhanced wavepacket branching at a conical intersection, suggests that different decoherence mechanisms act very differently on conical-intersection dynamics; a comparative study of amplitude damping versus dephasing in this trapped-ion setting could clarify when the conical signature survives.","Because the authors model only one decay channel, the inclusion of cascade decay through intermediate states could introduce additional dephasing of the spin–phonon correlations; a master equation with a multi-level decay ladder is a testable extension that would show whether the visibility window shrinks."],"forward_implications":["Within the Rydberg lifetime, the phonon position and the spin–phonon correlations still show clear oscillations, so an experimental run can resolve conical-intersection dynamics before decay erases the collective spin sector.","Collective spin observables such as $\\langle S_z\\rangle$ and $\\langle xS_z\\rangle$ are the better measurement channel: their oscillation amplitudes are much larger than those of single-ion Rydberg populations, which oscillate only weakly.","The phonon populations relax to their steady values even though the dissipator contains no phonon loss, because the spin–phonon coupling converts the electronic decay into damping of the vibrational modes.","The decay of the nP state $|1\\rangle$ can be neglected for the conical-intersection dynamics, since its lifetime is one to two orders of magnitude longer than that of the nS state $|0\\rangle$; including it would not change the relevant timescale.","The conical-intersection signal is transient by nature: once the $|0\\rangle$ population has decayed, the collective spin sector no longer exists, so experiments should target the first several microseconds of evolution."],"supporting_citations":[{"why":"Supplies the trapped-Rydberg-ion scheme and the Hamiltonian that creates the conical intersection, which this work quantizes by promoting the vibrations to phonon operators.","marker":"[30]"},{"why":"The superconducting-circuit conical-intersection experiment with dephasing that this paper contrasts with its amplitude-damping approach.","marker":"[26]"},{"why":"Model-potential calculation of Sr+ Rydberg properties used to obtain the lifetimes and level structure in the simulation.","marker":"[59]"},{"why":"Provides the coherent-state definition and the master-equation formalism used for the dissipative dynamics.","marker":"[61]"},{"why":"Standard open-systems master-equation theory underlying the Lindblad dissipator for spontaneous decay.","marker":"[65]"},{"why":"Numerical master-equation solver used to compute the exact dissipative time evolution.","marker":"[66]"},{"why":"Source of Rydberg-state lifetimes in group-II ions used to set the decay rate $\\gamma_S=0.13\\ \\mu\\text{s}^{-1}$.","marker":"[67]"}],"fun_headline_variants":["Conical-intersection dynamics survive Rydberg decay","Trapped Rydberg ions keep conical signal visible","Oscillatory conical dynamics within Rydberg lifetime","Rydberg decay fails to erase conical intersection trace","Conical intersection observed before Rydberg state dies"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the relevant dissipation on the microsecond timescale is a single Markovian decay channel from the nS Rydberg state $|0\\rangle$ to a passive ground state $|g\\rangle$ at rate $\\gamma_S=0.13\\ \\mu\\text{s}^{-1}$, with intermediate decay steps, micromotion, and motional heating neglected; if those produce comparable decoherence, the predicted visibility window for the conical-intersection oscillations could close.","fun_headline_variants_meta":{"raw":{"variants":["Conical-intersection dynamics survive Rydberg decay","Trapped Rydberg ions keep conical signal visible","Oscillatory conical dynamics within Rydberg lifetime","Rydberg decay fails to erase conical intersection trace","Conical intersection observed before Rydberg state dies"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000369,"raw_usage":{"total_tokens":1960,"prompt_tokens":910,"completion_tokens":1050,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":526,"completion_tokens_details":{"reasoning_tokens":972}},"tokens_in":526,"tokens_out":1050,"duration_ms":9263,"temperature":1.0,"reasoning_tokens":972,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T10:34:47.702244+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete falsifier would be a two-Sr+-ion experiment prepared in $|0,1\\rangle=|50S,50P\\rangle$ with the $x$-phonon mode in a coherent state: the model predicts that $\\langle xS_z\\rangle$ and $\\langle N_x\\rangle$ oscillate and damp on the ~7 $\\mu$s 50S lifetime. If the measured spin–phonon correlations disappear within about one microsecond, or if the $|0\\rangle$ population decays appreciably faster than $e^{-\\gamma_S t}$, then the single-channel decay model omits a dominant decoherence mechanism and the visibility claim fails.","supporting_citations":[{"cited_title":"Exploring the Many-Body Dynamics Near a Conical Intersection with Trapped Rydberg Ions","cited_arxiv_id":null,"evidence_quote":"Supplies the trapped-Rydberg-ion scheme and the Hamiltonian that creates the conical intersection, which this work quantizes by promoting the vibrations to phonon operators."},{"cited_title":"Observation of Wave-Packet Branching through an Engineered Conical Intersection","cited_arxiv_id":null,"evidence_quote":"The superconducting-circuit conical-intersection experiment with dephasing that this paper contrasts with its amplitude-damping approach."},{"cited_title":"Qutip Documentation and Coding","cited_arxiv_id":null,"evidence_quote":"Numerical master-equation solver used to compute the exact dissipative time evolution."},{"cited_title":"Lifetimes of Rydberg States in Ions of the Group II Elements.Opt","cited_arxiv_id":null,"evidence_quote":"Source of Rydberg-state lifetimes in group-II ions used to set the decay rate $\\gamma_S=0.13\\ \\mu\\text{s}^{-1}$."}],"review_version":1}