{"id":"c7753e42-4392-499a-a780-7b521309bd45","arxiv_id":"2607.05216","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":3,"one_line_summary":"Interaction-induced frequency shifts on two-atom Rydberg transitions are analytically derived and corrected with a two-parameter phase modulation, restoring circularization fidelity above 95% at 65 ns for atom pairs down to 6.5 μm separation.","lead":"This paper shows that dipole-dipole interactions between two Rydberg atoms disrupt fast circularization pulses by shifting transition frequencies, and provides an analytical two-parameter phase correction that restores fidelity above 95% for interatomic distances down to 6.5 μm. A smart generalist might read it because it removes a key bottleneck for using interacting circular Rydberg atoms in quantum computing and simulation.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"No significant objection identified. The M⁻=0 approximation is acknowledged, perturbatively bounded (App. E), and validated by full numerical simulation including M⁻≠0 states, with <0.1% gain from parameter optimization confirming near-optimality.","rationale":"The reader correctly identifies the M⁻=0 approximation as the weakest link in the argument. However, this concern does not undermine the central claim because: (1) the full numerical simulations include M⁻≠0 states, so the analytical parameters are validated in a setting that does not assume M⁻=0; (2) the <0.1% improvement from parameter optimization confirms near-optimality; (3) the perturbative bounds in App. E place the breakdown well below the claimed threshold. The paper is transparent about limitations, the derivation is parameter-free (not fitted), and the numerical verification is thorough across R, θ, and t_f. The minor distance threshold inconsistency (6.3 vs 6.5 μm) is cosmetic. The ACCEPT verdict with HIGH confidence is appropriate. The correctness risk remains 'unknown' only in the sense that experimental validation is still pending, but the theoretical framework is internally consistent and well-supported.","tokens_in":27616,"tokens_out":2949,"duration_ms":36446,"concrete_test":"To further validate the M⁻=0 approximation specifically for intra-ladder dynamics (as distinct from the inter-ladder check in App. E), compute the interaction-induced frequency shifts for pair states with M⁻≠0 (e.g., M⁻=2,4) using the same hydrogenic framework as Eq. (4.1), and compare them to the M⁻=0 shifts. If the relative difference exceeds ~10% for M⁻ values within the populated range (|M⁻| ≤ 28) at R = 6.5 μm, the analytical correction would be missing non-negligible physics even in the claimed regime. This is an analytical calculation that could be done without new simulations.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The reader identifies the synchronous-evolution / M⁻=0 approximation as the weakest assumption. This is a reasonable identification, but the concern does not land as load-bearing for the central claim (R ≥ 6.5 μm) for three reasons. First, the full numerical simulations used to verify p_CC include pair states with M⁻ up to 28 (Appendix F, Table F.I), so the analytical parameters α_int and β_int—derived assuming M⁻=0—are being tested in a simulation that does not assume M⁻=0. The fact that p_CC ≥ 95% is achieved in these simulations directly validates the approximation in the claimed regime. Second, the paper shows that numerically optimizing α_int and β_int yields <0.1% improvement over the analytical values for all R (Sec. V, Fig. 5), confirming that the two-parameter linear correction captures the dominant interaction effect. If the M⁻=0 approximation were causing significant missing physics, we would expect the optimized parameters to deviate substantially from the analytical ones. Third, Appendix E provides a perturbative bound on inter-ladder couplings, showing R_min ≤ 3.1 μm for all relevant terms—well below the 6.5 μm threshold. The paper is transparent about the breakdown regime (R < 6.5 μm, where p_CC drops to 62.7% at R = 4 μm) and does not claim validity there. One could note that the perturbative analysis in App. E addresses inter-ladder couplings rather than intra-ladder ΔM⁻≠0 transitions directly, but the numerical evidence (sparse M⁻ population in Fig. F.1, near-optimality of analytical parameters) compensates for this gap. No internally inconsistent or unsupported step was found in the derivation chain from Eq. (4.1) through Eq. (5.1) to the numerical verification.","agreement_with_reader":"agree"},"referee_report":{"model":"glm-5.2","summary":"This manuscript addresses the problem of fast circularization of two interacting Rydberg atoms, a key bottleneck for quantum simulation and computation with circular Rydberg states. The authors identify interaction-induced frequency shifts on the lowest-diagonal-ladder transitions as the dominant mechanism by which dipole-dipole interactions degrade single-atom circularization pulses. They derive analytically (in the hydrogenic approximation) a two-parameter linear correction to the transition frequencies (Eqs. 4.1–4.2), and use it to construct an adapted pulse (Eq. 5.1) by phase-modulating a pre-optimized single-atom pulse. The adapted pulse is shown numerically to restore circular pair-state probabilities ≥95% for two 87Rb atoms at n=52, for interatomic distances down to 6.5 μm and all angular configurations, within 65 ns and under experimental amplitude/frequency constraints. A multi-layer optimization scheme combining the analytical adaptation with Krotov's method is also proposed, extending high-fidelity operation to R > 5.9 μm.","tokens_in":27791,"tokens_out":2250,"duration_ms":55753,"significance":"The paper makes a valuable contribution by extending fast circularization protocols—previously developed only for single atoms—to the interacting two-atom regime, which is essential for quantum simulation and computation with circular Rydberg atoms. The analytical derivation of the correction parameters α_int and β_int (Eq. 4.2) is parameter-free, obtained from the dipole-dipole Hamiltonian and hydrogenic energy expressions without fitting to the target result. The key approximation (synchronous evolution on the M⁻=0 subspace) is validated by full numerical simulations that include pair states with M⁻ up to 28 (Appendix F, Table F.I), and the near-optimality of the analytical parameters is confirmed by the <0.1% improvement obtained from numerical optimization of α_int and β_int (Sec. V). The perturbative bounds in Appendix E provide additional support by showing that inter-ladder couplings are negligible for R well below the claimed threshold. The proposed multi-layer optimization scheme is practical and experimentally motivated, with explicit enforcement of amplitude and bandwidth constraints.","major_comments":[],"minor_comments":[{"comment":"There is an inconsistency in the stated distance threshold for the 95% fidelity claim: the abstract and Sec. V state 'down to 6.5 μm' / 'larger than 6.5 μm', while the conclusion (Sec. VII) states 'interatomic distances above 6.3 μm'. These should be reconciled.","section":null},{"comment":"Sec. V, Fig. 4(b) caption: the pulse is labeled 'F_aa^{ad.freq.}' in the text but the caption refers to 'F_a^{Krotov}' in green. The caption should clarify that the green curve is the single-atom optimized pulse F_a^{opt}, not a Krotov-optimized pair pulse.","section":null},{"comment":"Eq. (4.1): the expression for ΔE_{M+→M++2} contains terms from both the second-order Stark shift and the dipole-dipole interaction, but the separation into α_0, β_0 (non-interacting) and α_int, β_int (interaction) in Eq. (4.2) is only made explicit after the equation. A brief inline annotation of which terms in Eq. (4.1) correspond to which parameters would improve readability.","section":null},{"comment":"Sec. III, Fig. 2(b): the inset is mentioned but the scale of the interaction-induced population trapping is not easily discernible. Consider enlarging or adding a quantitative label for the trapped population fraction.","section":null},{"comment":"Appendix E, Table E.I: the term ∝ Ĵ^{m,1}_- Ĵ^{p(m),2}_z + Ĵ^{p(m),1}_z Ĵ^{m,2}_- lists R_min = 1.8 μm, but the text states 'Most of the terms lead to negligible couplings with the exception of the term ∝ Ĵ^{m,1}_- Ĵ^{p,2}_- + Ĵ^{p,1}_- Ĵ^{m,2}_-' with R_min = 3.1 μm. It would help to explicitly state in the table which term is the most restrictive.","section":null},{"comment":"Sec. VI, Fig. 5(a): the y-axis label 'Circular pair state prob. p_CC' could include the target line at p_CC = 0.98 (the non-interacting maximum) for reference, to make the gap to the theoretical maximum visible.","section":null},{"comment":"The conclusion mentions 'interatomic distances above 6.3 μm' for the adapted pulse, but Sec. VI states that Krotov's method with constraints yields differences below 1% for R > 5.7 μm. The distinction between the adapted-pulse-only threshold and the Krotov-refined threshold could be stated more sharply in the conclusion.","section":null},{"comment":"A few typographical issues: 'circulrization' in Fig. 2(b) caption; 'longe-range' in the Introduction; 'constrains' used in place of 'constraints' in several places (e.g., Sec. VI, Appendix D).","section":null}],"recommendation":"minor_revision","confidential_remarks":"The paper is a solid theoretical contribution with clean analytical work and thorough numerical validation. The central claim is well-supported and the approximation regime is transparently stated and bounded. I see no load-bearing issues. The minor inconsistencies in the distance threshold (6.3 vs 6.5 μm) and presentation issues should be addressed but do not affect the scientific conclusions. The fit with the journal's scope is good given the combination of analytical derivation, numerical simulation, and experimental constraint awareness."},"author_rebuttal":{"model":"glm-5.2","summary":"We thank the referee for a careful and constructive reading of the manuscript and for the positive assessment of our work. The referee's summary accurately captures the main contributions: the analytical derivation of the interaction-induced frequency shifts, the two-parameter phase-modulation correction, and the multi-layer optimization scheme. We note that the referee report contains no major comments, specific revision requests, or criticisms requiring point-by-point response. The recommendation is minor revision, but no specific minor points have been enumerated. We have carefully re-read the report and the manuscript in light of the referee's assessment and confirm that the manuscript already addresses the points the referee highlights as strengths: the perturbative bounds in Appendix E, the state-space convergence analysis in Appendix F (Table F.I), and the near-optimality check of the analytical parameters (Section V). We are not aware of any substantive changes the referee is requesting beyond what is already present. If the editor or referee has specific minor corrections in mind that were not included in the report as transmitted to us, we would be grateful to receive them and will incorporate them promptly.","responses":[],"tokens_in":26950,"tokens_out":232,"duration_ms":18447,"standing_objections":[]},"desk_editor":{"model":"glm-5.2","letter":"Bottom line: this paper takes the single-atom fast circularization pulses from Patsch et al. (2018) and extends them to the interacting two-atom regime. The key new result is an analytical, parameter-free derivation of interaction-induced frequency shifts on the relevant two-body transitions (Eq. 4.1–4.2), and a simple phase-modulation correction (Eq. 5.1) that restores the single-atom pulse fidelity to ≥95% for pairs at distances down to 6.5 μm, all angles, 65 ns duration. The 60× speedup over adiabatic passage is real and matters for the circular Rydberg experimental program—especially because the method works at arbitrary angles, not just the magic angle, which is a genuine limitation of current experimental approaches for 2D geometries. The analytical parameters α_int and β_int are derived from the dipole-dipole Hamiltonian and hydrogenic energy expressions without fitting to the target result. The paper then confirms near-optimality: numerical optimization of those two parameters yields <0.1% improvement. That is a strong internal consistency check. The full numerical simulations include pair states with M⁻ up to 28 (Appendix F), so the analytical correction—derived assuming M⁻ = 0—is being tested in a simulation that does not assume M⁻ = 0. The fact that p_CC ≥ 95% holds in these simulations validates the approximation in the claimed regime. The perturbative bound in Appendix E (R_min ≤ 3.1 μm for inter-ladder couplings) is consistent with the observed breakdown below 6.5 μm. The paper is transparent about where the method fails: p_CC drops to 62.7% at R = 4 μm, and the synchronous-evolution picture breaks down there. One minor inconsistency: the abstract says 6.5 μm, the conclusion says 6.3 μm. This does not affect the central claim. The work is purely theoretical/simulation-based. Experimental validation is stated as next step. The theoretical framework is sound and complete as presented. The multi-layer optimization scheme (Sec. VI) is a practical contribution: the adapted pulse serves as a good initial guess for Krotov optimization, reducing computational cost. This is useful for experimental groups even if they ultimately need the full optimization pipeline. The paper is for researchers working on circular Rydberg atom experiments and quantum control—theoretical groups designing pulses and experimental groups planning 2D arrays. It deserves a serious referee. The derivation chain from Eq. 4.1 through Eq. 5.1 to the numerical verification is clean, the limitations are honestly stated, and the result removes a real bottleneck for interacting circular Rydberg atoms.","headline":"Solid theoretical contribution: analytically derived two-parameter phase correction restores ≥95% circularization fidelity for interacting Rydberg pairs at 65 ns, with clean numerics and honest limitations.","tokens_in":28554,"tokens_out":650,"would_cite":true,"duration_ms":27190,"reading_group":"no","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"glm-5.2","headline":"Phase trick restores 95% circularization of interacting Rydberg atoms","keywords":[],"falsifier":"Experimental measurement of circular pair-state probability for two interacting 87Rb atoms at n=52, R ~ 7 micrometers, theta = 0, with the adapted pulse applied. If the measured p_CC falls substantially below the predicted ~97%, the analytical frequency-shift model or the synchronous-climbing approximation does not capture the dominant physics. A second falsifier: if the adapted pulse's performance does not follow the (1 - 3 cos^2 theta)/R^3 scaling across angles and distances, the identified mechanism is not the primary disturbance.","tokens_in":27789,"feed_emoji":"🌀","tokens_out":836,"duration_ms":39223,"temperature":0.7,"pith_summary":"Circular Rydberg atoms are prized for quantum simulation, computation, and sensing, but their preparation — circularization, the transfer of ~50 angular momentum quanta to the valence electron — has been a bottleneck. Prior fast-pulse methods worked for isolated atoms but failed when interatomic dipole-dipole interactions were present, exactly the regime quantum technologies require. This paper identifies the dominant disturbance as interaction-induced shifts to transition frequencies along the so-called lowest diagonal ladder of Rydberg states, derives those shifts analytically in a hydrogenic picture as two linear parameters (alpha_int, beta_int), and corrects for them by multiplying a pre-optimized single-atom RF pulse by a time-dependent phase factor that tracks the instantaneous frequency shift. The resulting adapted pulse, requiring no numerical optimization beyond the single-atom pulse already in hand, restores circular pair-state fidelity to at least 95% for two rubidium-87 atoms at n=52 separated by as little as 6.5 micrometers, in any angular orientation, at 65 ns duration — a 60-fold speedup over the standard adiabatic passage. Combined with Krotov's gradient-based pulse shaping as a refinement layer, the method reaches 99% fidelity for distances above 5.9 micrometers.","feed_headline":"Phase trick restores 95% circularization of interacting Rydberg atoms","feed_subtitle":"Analytical frequency-shift correction lets fast RF pulses prepare two interacting atoms in 65 ns, a 60x speedup over standard methods.","key_machinery":"The adapted pulse F_aa^{ad.freq.}(t) = F_a^{opt}(t) * exp(-i * integral of Delta_omega(tau) d tau), with Delta_omega(M+) = (alpha_int + beta_int * M+) / hbar, where alpha_int and beta_int are given in closed form by Eq. (4.2) and depend on (1 - 3 cos^2 theta)/R^3. The phase integral uses the single-atom trajectory m_bar_opt(t) as proxy for M+/2.","core_discovery":"The central mechanism is that dipole-dipole interactions between two Rydberg atoms shift the transition frequencies along the lowest diagonal ladder by an amount that is linear in the total magnetic quantum number M+ and scales as (1 - 3 cos^2 theta) / R^3. Because the single-atom pulse was designed for unshifted frequencies, these shifts drive it off resonance and destroy fidelity. The correction is a phase modulation F_aa(t) = F_a_opt(t) * exp(-i * integral of Delta_omega(t) dt), where Delta_omega depends on the analytically derived parameters alpha_int and beta_int and on the single-atom trajectory m_bar(t). This correction is constructed entirely from the single-atom pulse and the known氢","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["Interaction shifts break single-atom circularization; phase fix restores 95% fidelity","Fast RF pulses circularize interacting Rydberg atoms in 65 ns with 95% fidelity","Two-parameter phase fix cancels dipole shifts in interacting Rydberg circularization","Analytical phase modulation recovers 95% fidelity for two-atom Rydberg circularization","Correcting dipole-induced frequency shifts enables fast circularization of atom pairs"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The analytical correction assumes both atoms climb their respective ladders synchronously, populating only pair states with m1 = m2, which reduces the two-body dynamics to two simultaneous single-atom transitions. This breaks down at short distances where interaction terms coupling to states with m1 != m2 become significant, and indeed the fidelity drops to 62.7% at R = 4 micrometers.","fun_headline_variants_meta":{"raw":{"variants":["Interaction shifts break single-atom circularization; phase fix restores 95% fidelity","Fast RF pulses circularize interacting Rydberg atoms in 65 ns with 95% fidelity","Two-parameter phase fix cancels dipole shifts in interacting Rydberg circularization","Analytical phase modulation recovers 95% fidelity for two-atom Rydberg circularization","Correcting dipole-induced frequency shifts enables fast circularization of atom pairs"]},"model":"glm-5.2","effort":"high","cost_usd":0.0,"raw_usage":{"total_tokens":819,"prompt_tokens":710,"completion_tokens":109,"prompt_tokens_details":null},"tokens_in":710,"tokens_out":109,"duration_ms":8097,"temperature":1.0,"reasoning_tokens":null,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-07T23:01:51.252896+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"Experimental measurement of circular pair-state probability for two interacting 87Rb atoms at n=52, R ~ 7 micrometers, theta = 0, with the adapted pulse applied. If the measured p_CC falls substantially below the predicted ~97%, the analytical frequency-shift model or the synchronous-climbing approximation does not capture the dominant physics. A second falsifier: if the adapted pulse's performance does not follow the (1 - 3 cos^2 theta)/R^3 scaling across angles and distances, the identified mechanism is not the primary disturbance.","supporting_citations":[],"review_version":1}