Mean-field Pulse Adaptation for the Circularization of Interacting Rydberg Atoms

Eloisa Cuestas; Felix Motzoi; Matthias H\"uls; Tommaso Calarco

REVIEW 2 major objections 4 minor 16 references

Mean-field Pulse Adaptation for the Circularization of Interacting Rydberg Atoms

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

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

Pith's one-line read For two interacting 87Rb atoms at 7 µm, a mean-field adaptation of a single-atom pulse restores circularization fidelity from 66% to 98.3%.

desk verdict The mean-field pulse adaptation works for the two-atom case as claimed; just demand the missing derivation for the phase coefficient in Eq. (13) before publication. read the letter →

arxiv 2607.26978 v2 pith:Y5CB2ODA submitted 2026-07-29 quant-ph physics.atom-ph

Matthias Hüls , Felix Motzoi , Tommaso Calarco , Eloisa Cuestas This is my paper · ORCID

classification quant-phphysics.atom-ph MSC 81V8081P68

keywords mean-fieldapproximationRydbergcircularizationpulseadaptationdipole-dipoleinteractionquantumcontrol87RbarraysHartree

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper introduces a mean-field (semi-classical) model that simulates the circularization of interacting Rydberg atoms in a Hilbert space of dimension N*d instead of d^N, where d is the single-atom state space and N the atom count. The model treats the dipole-dipole interactions of other atoms as a classical time-dependent field and propagates each atom in parallel. From the same single evolution, the method produces an adapted pulse that adds amplitude and phase corrections canceling the interaction field. For pairs of 87Rb atoms in n=52 states separated by ≥7 µm, the effective model's final infidelity stays below 1% and the adapted pulse restores the mean circular-state probability from about 66% to 98.3% in an exact simulation.

What carries the argument

The load-bearing object is the single-atom effective Hamiltonian h_eff,i(t) = h_s - d·[F(t) + I_i(t)], where I_i(t) is a classical interaction field built from the time-dependent mean dipole moments of all other atoms (Eqs. 4–5). The pulse adaptation is carried by the complex-pulse equation f_ad(t_k) = exp(-i φ_k) f_sa(t_k) - I_x - i I_y with φ_k = (3n ea0/(2ℏ)) Δt Σ I_z, which subtracts the transverse field as an amplitude offset and integrates the longitudinal field into a chirp. This turns the exponential Hilbert-space problem into N*d parallel single-atom propagations and yields the adapted pulse as a by-product of one time evolution.

What would settle it

Perform an exact two-body simulation at R=5 µm for θ=0 and compare the effective model's predicted final overlap and adapted-pulse fidelity against the exact numbers; the paper's own scaling predicts the <1% error claim to break down, and a measured mismatch of the adapted pulse's exact fidelity from the predicted 95% would confirm the entanglement-driven limit. A more direct test is to compute the von Neumann entropy of the single-atom reduced state during the exact evolution — nonzero entropy quantifies the entanglement the mean-field model drops.

Watch

Extended reading notes

Core claim

The central discovery is that a Hartree product-state approximation — each atom evolving in the classical field generated by the mean dipole moments of the others — reproduces the exact two-body circularization dynamics to better than 99% overlap for interatomic distances R≥7 µm, and that the same calculation yields a pulse adapted to the interactions in one pass. The adaptation consists of subtracting the transverse interaction field components (amplitude compensation) and integrating the longitudinal component as a time-dependent phase chirp. With this adapted pulse, the mean circular-state probability recovers from 65.9% (unadapted) to 98.3% in an exact two-atom simulation at R=7 µm, θ=0.

Load-bearing premise

The load-bearing premise is the Hartree product-state ansatz — the assertion that the multi-atom wavefunction stays factorizable at all times; when interatomic entanglement becomes non-negligible (as at R=4 µm, where the model's error reaches 46.7%), both the simulated evolution and the adapted pulse lose validity.

Editorial extensions

If this is right

Preparing circular Rydberg arrays beyond two atoms becomes numerically tractable: the simulation cost grows as N*d rather than d^N, enabling pulse design for large arrays.
Adapted pulses can be generated on the fly during a single mean-field evolution, avoiding iterative optimal-control loops for weakly to moderately interacting systems.
For atom pairs separated by ≥7 µm, the adapted pulses restore the single-atom performance in exact simulation to ≥98% mean circular probability.
At the magic angle θ=acos(1/√3), where the z-z interaction cancels, amplitude-only adaptation recovers 99% even at R=4 µm, showing that the detuning chirp is the main limiting mechanism.
In the strong-interaction regime, the adapted pulse remains useful as a warm start for a subsequent pulse-shaping optimization.

Reading between the lines

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

Extension: the same one-pass adaptation should work for moving atoms (flying qubits), since the mean-field computation does not assume fixed geometry; a testable prediction is that a time-dependent R(t) produces a smoothly chirped pulse that preserves about 99% fidelity for trajectories that keep R≥7 µm.
Testable extension: for geometrically asymmetric N>2 arrays, the spread of the local fields I_i(t) could serve as a quantitative figure of merit for how much the global-pulse adaptation will degrade; arrays arranged at the magic angle should show the smallest spread.
Neighbouring problem: the phase chirp breakdown at strong interactions resembles a breakdown of phase-locking in driven anharmonic systems; adapting the chirp rate dynamically from the instantaneous I_z rather than a fixed integral might extend the method to smaller R.
If the Hartree ansatz is replaced by a one-axis-twisting or Gaussian-state correction, the same I_i(t) machinery could incorporate the leading entanglement corrections, extending the <1% regime to shorter distances.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

read the letter

Core result holds up. The paper's central demonstration—effective model matches exact two-body evolution below 1% infidelity for R≥7 µm, and the adapted pulse recovers 98.3% circular probability in exact simulation—is credible and well presented. The Hartree approximation is the right tool for this weak-to-moderate interaction regime, and Fig. 1(d) honestly shows where it breaks down. This is genuinely new relative to Ref. [15]: instead of optimizing a couple of frequency parameters, they compute amplitude and phase corrections on the fly from a single mean-field evolution, and the O(Nd) scaling is the useful selling point.

The weak spot that will bite in review is Eq. (13). The coefficient 3n/2 is stated without derivation or citation, and the compensation is the limiting factor at strong interactions. I don't think it's wrong—the exact simulation at R=7 µm, θ=0.1π would catch a gross error—but it's an asserted effective average over many straddled Stark sublevels, and there's no sensitivity check. A referee should ask for one or two lines of derivation, or at minimum a scan over the coefficient showing the 98.3% is not a knife-edge. That's an omitted derivation, not a demonstrated flaw.

The other limitations are honestly flagged in the text: the model fails for strong interactions (46.7% infidelity at R=4 µm), and the N>2 extension rests on the heuristic average field in Eq. (14). The authors admit asymmetric arrays remain unsolved. That's fine for a first paper if the claims are scoped.

No code or data is included, which hurts reproducibility for a numerical methods paper. I'd like the authors to release the propagation and the pulse adaptation code, especially given the terse description of simulation details ("see [15]"). Not fatal, but it makes the sensitivity check harder for the reader.

Who gets value: groups working on circular Rydberg arrays and quantum control. It deserves a serious referee, not a desk reject. My recommendation: send to peer review with an explicit request to justify Eq. (13) and provide a sensitivity analysis, plus archive code. If the coefficient derivation is provided, the paper is in good shape.

Referee Report

2 major / 4 minor

Summary. The paper introduces a mean-field effective model for simulating the circularization of interacting Rydberg atoms. The model replaces the dipole-dipole interaction with classical fields generated by the mean dipole moments of the other atoms, reducing the state-space dimension from d^N to N d. Using this model, the authors adapt an existing single-atom optimal pulse by adding amplitude and phase compensations that cancel the interaction fields at each time step. They validate the effective model against exact two-body simulations for two 87Rb atoms in n=52 states, reporting that the final state overlap infidelity remains below 1% for interatomic distances R ≥ 7 µm (Fig. 1d). For a pair at R = 7 µm, θ = 0.1π, the adapted pulse F_ad(t) recovers a mean circular-state probability of 98.3% in the exact simulation, close to the 99% target of the single-atom pulse. The paper also discusses the extension to N > 2 atoms via an average-field adaptation, acknowledging that asymmetric arrays remain an open problem.

Significance. If the result holds, the method offers a scalable route to simulating and controlling interacting Rydberg arrays without exponential Hilbert-space growth, which is directly relevant to quantum simulation and computing platforms. A notable strength is that the adapted pulse is computed from the mean-field model and then tested against an independent exact two-body simulation, with no parameters fitted to the target results. The authors are also transparent about the Hartree product-state limitation and the heuristic nature of the N > 2 extension. However, the central numerical demonstration is only for N = 2, and the most critical element of the phase compensation—the coefficient in Eq. (13)—is asserted without derivation or a sensitivity analysis. This currently limits the confidence in the generality of the method.

major comments (2)

[Eq. (13)] The phase-modulation coefficient C = (3n/2)(e a0/ℏ) in Eq. (13) is stated without derivation or citation. This coefficient determines the integrated phase compensation that is crucial for the adapted-pulse performance; the exact-simulation value p̄C = 98.3% at R = 7 µm depends directly on its correctness. During the chirp, the population occupies many Stark sublevels, so the differential frequency shift per unit I_z is generally state-dependent. Please provide a derivation of C (e.g., from the linear Stark shifts of the relevant |n,m⟩ states) or demonstrate via a sensitivity scan (varying C by ±10% and recomputing p̄C in the exact simulation) that the result is robust to its precise value. Without this, the adaptation scheme is not fully specified.
[Abstract / Conclusion] The abstract claims the method 'enables the simulation of large atomic systems,' and the conclusion proposes simulating 'arbitrary arrays of N interacting atoms.' However, the only numerical demonstration is for N = 2, and it exploits permutation symmetry to reduce the simulation to a single atom. For N > 2, the average-field adaptation in Eq. (14) is explicitly heuristic, and the authors state that asymmetric arrays remain unsolved. To support the generality claimed in the title and abstract, please add at least one simulation for a small N > 2 system (e.g., a symmetric square array of four atoms), or alternatively restrict the stated performance claims to the two-atom case and clearly label the N > 2 portion as a proposal.

minor comments (4)

[Abstract / Sec. 'Mean-field pulse adaptation'] The claim that the error of the method remains below 1% refers to the state infidelity 1 − O(t_f) defined in Eq. (10), not directly to the error in the circular-state probability p̄C. The paper should clarify this in the abstract and results to avoid misinterpretation, especially since the p̄C differences in Fig. 2 can be larger than 1% in some regimes.
[Eq. (6)] The numerical parameters of the simulation—particularly the time step Δt and the size of the truncated Hilbert space—are not reported in the main text. Since the effective model propagates with a piecewise-constant Hamiltonian, the accuracy of the results depends on these parameters. Please state them or point to the specific section of Ref. [15] where they are given.
[Fig. 1(d) caption] The caption refers to 'the term ∝ ⟨ẑ_j⟩ of Eq. (9)', but the operator ⟨ẑ_j⟩ is not defined. This presumably refers to the ⟨d̂_zj⟩ term in Eq. (9). Please use a consistent notation.
[Section 'Following a semi-classical approach'] The statement that spin-orbit coupling is neglected is made without justification. For n = 52 Rb, this is likely a reasonable approximation, but a one-sentence justification or a reference would be helpful for the reader.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the adapted pulse is a feedforward compensation validated against an independent exact two-body simulation.

full rationale

The paper's central derivation is a mean-field (Hartree) effective model, Eqs. (4)-(6), which is benchmarked against an exact two-body simulation via the overlap O defined in Eq. (10). This comparison is not circular: the effective model is not fitted to the exact result, and the 0.8% infidelity is an independent numerical test. The pulse adaptation, Eqs. (11)-(13), is a feedforward scheme that subtracts the mean-field interaction field and adds a phase chirp; the effective-model recovery of 99.0% is indeed largely by construction, since F_ad is designed to cancel I(t) within the mean-field approximation. However, the paper does not stop there: it validates F_ad in an exact two-body simulation, finding p̄_C=98.3% versus 65.9% for the unadapted pulse. That exact-simulation test is the load-bearing evidence and is not self-referential. The coefficient (3n/2)(e a0/ℏ) in Eq. (13) is asserted without derivation, but it is a standard differential Stark-shift coefficient, not a parameter fitted to the target circular-state probability; at most this is a missing derivation/correctness risk, not circularity. The only self-citation is Ref. [15] (same group), used for the input pulse F_sa and simulation details; F_sa is an input to the method, not a predicted conclusion, so this self-citation is not load-bearing. Overall, the derivation chain is self-contained against an external benchmark and no circular step is exhibited.

Assumptions & free parameters 0 free parameters · 5 assumptions · 1 invented entities

No fitted parameters were introduced; static fields, distances, and angles are test scenarios inherited from Ref. [15], and F_sa(t) is a fixed input pulse. The main axioms are the Hartree approximation, the dipole-dipole Hamiltonian, neglection of spin-orbit coupling, permutation symmetry for the pair demonstration, and the unstated Stark-shift coefficient in Eq. (13).

assumptions (5)

domain assumption Hartree product-state approximation: |Ψ(t)> ≈ ⊗_i |ψ_i(t)> throughout the evolution
Core of the effective model, stated after Eq. (3). The paper tests its validity for N=2 but assumes it; Fig. 1d shows it breaks down for strong interactions.
domain assumption Classical dipole-dipole interaction Hamiltonian (Eq. 3) with uniform static fields and no propagation delays
Standard dipole coupling used to define mean fields in Eqs. (4)-(9); uniform-field/no-delay assumption is stated in the text after Eq. (2).
domain assumption Neglect of spin-orbit coupling
Stated after Eq. (1); simplifies the Hamiltonian to single-electron orbital dynamics.
domain assumption Permutation symmetry for two atoms: |ψ_1(t)> = |ψ_2(t)>
Used to reduce the two-atom validation to a single-atom simulation; exact symmetric initial state and identical atoms preserve exchange symmetry, but this is specific to the demonstrated configuration.
domain assumption Phase-modulation coefficient in Eq. (13) equals (3n/2)(e a0/ħ) I_z, converting interaction-induced detuning into a chirp
Asserted without derivation; relies on hydrogenic circular-state Stark shift. If this linear-response coefficient is inaccurate or amplitude-dependent, the adapted-pulse results change.

invented entities (1)

Mean-field interaction field I_i(t)
purpose: Effective classical field representing dipole-dipole interactions in the single-atom evolution Eq. (5)
A computational construct defined by expectation values of physical dipole operators (Eq. 4), not a new physical entity with independent falsifiable predictions.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Mean-field Pulse Adaptation for the Circularization of Interacting Rydberg Atoms." pith.science (2026). https://pith.science/paper/Y5CB2ODA

@misc{pith2026260726978,
  author       = {Pith},
  title        = {Pith review of: Mean-field Pulse Adaptation for the Circularization of Interacting Rydberg Atoms},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/Y5CB2ODA}},
  note         = {Machine review of arXiv:2607.26978}
}

abstract

Arrays of circular Rydberg atoms provide a promising platform for quantum simulation and computation; however, their preparation in the presence of interatomic interactions remains a major challenge. While optimal control methods have enabled the design of fast and accurate radio-frequency pulses for the circularization of a single atom and of an atom pair, the extension to more atoms is fundamentally limited by the exponential growth of the Hilbert space, which renders numerical simulations computationally infeasible. Here, we introduce an effective model that treats interactions within a mean-field approximation, thereby enabling the simulation of large atomic systems. Our model further enables the adaptation of pulses optimized for non-interacting atoms to interacting systems, based on the computation of a single time evolution. For two interacting $^{87}\mathrm{Rb}$ atoms, we demonstrate that the error of our method remains below $1 \, \%$ and that our adapted pulses recover the initial performance of optimal pulses in the regime of weak to moderate interaction strengths.

Figures

Figures reproduced from arXiv: 2607.26978 by the authors.

**Figure 2.** FIG. 2. Mean circular state probability ¯p [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

Review history (2 revisions) →

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

16 extracted references

[15]

H¨ uls, A

M. H¨ uls, A. A. Young, C. Sayrin, M. Brune, J.-M. Rai- mond, T. Calarco, F. Motzoi, R. Zeier, and E. Cuestas, Fast pulses for high-fidelity circularization of interacting Rydberg atoms (2026)

2026
[4]

S. R. Cohen and J. D. Thompson, PRX Quantum2, 10.1103/prxquantum.2.030322 (2021)

doi:10.1103/prxquantum.2.030322 2021
[6]

Meinert, C

F. Meinert, C. H¨ olzl, M. A. Nebioglu, A. D’Arnese, P. Karl, M. Dressel, and M. Scheffler, Physical Review Research2, 10.1103/physrevresearch.2.023192 (2020)

doi:10.1103/physrevresearch.2.023192 2020
[7]

Facon, E.-K

A. Facon, E.-K. Dietsche, D. Grosso, S. Haroche, J.- M. Raimond, M. Brune, and S. Gleyzes, Nature535, 262–265 (2016)

2016
[8]

E. K. Dietsche, A. Larrouy, S. Haroche, J. M. Raimond, M. Brune, and S. Gleyzes, Nature Physics15, 326–329 (2019)

2019
[9]

Nussenzveig, F

P. Nussenzveig, F. Bernardot, M. Brune, J. Hare, J. M. Raimond, S. Haroche, and W. Gawlik, Physical Review A48, 3991–3994 (1993)

1993
[10]

Signoles, E

A. Signoles, E. Dietsche, A. Facon, D. Grosso, S. Haroche, J. Raimond, M. Brune, and S. Gleyzes, Phys- ical Review Letters118, 10.1103/physrevlett.118.253603 (2017)

doi:10.1103/physrevlett.118.253603 2017
[11]

Patsch, D

S. Patsch, D. M. Reich, J.-M. Raimond, M. Brune, S. Gleyzes, and C. P. Koch, Physical Review A97, 10.1103/physreva.97.053418 (2018)

doi:10.1103/physreva.97.053418 2018

Show all 16 references

[12]

Larrouy, S

A. Larrouy, S. Patsch, R. Richaud, J.-M. Raimond, M. Brune, C. P. Koch, and S. Gleyzes, Physical Review X10, 10.1103/physrevx.10.021058 (2020)

2020 doi
[13]

M´ ehaignerie,Interactions entre atomes de Rydberg cir- culaires pi´ eg´ es pour la simulation quantique, Ph.D

P. M´ ehaignerie,Interactions entre atomes de Rydberg cir- culaires pi´ eg´ es pour la simulation quantique, Ph.D. thesis, Sorbonne universit´ e (2023)

2023
[14]

M´ ehaignerie, Y

P. M´ ehaignerie, Y. Machu, A. Dur´ an Hern´ andez, G. Creutzer, D. Papoular, J. Raimond, C. Sayrin, and M. Brune, PRX Quantum6, 10.1103/prxquan- tum.6.010353 (2025)

2025 doi
[16]

T. F. Gallagher,Rydberg Atoms, Cambridge Monographs on Atomic, Molecular and Chemical Physics (Cambridge University Press, 1994)

1994
[17]

Weber, C

S. Weber, C. Tresp, H. Menke, A. Urvoy, O. Firstenberg, H. P. B¨ uchler, and S. Hofferberth, Journal of Physics B: Atomic, Molecular and Optical Physics50, 133001 (2017)

2017
[18]

Ravets, H

S. Ravets, H. Labuhn, D. Barredo, T. Lahaye, and A. Browaeys, Physical Review A92, 10.1103/phys- reva.92.020701 (2015)

2015 doi
[19]

Hwang, A

H. Hwang, A. Byun, J. Park, S. de L´ es´ eleuc, and J. Ahn, Optica10, 401 (2023)

2023
[20]

Barth, L

I. Barth, L. Friedland, O. Gat, and A. G. Shagalov, Phys- ical Review A84, 10.1103/physreva.84.013837 (2011)

2011 doi

Pith tools

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

[1] [15]

H¨ uls, A

M. H¨ uls, A. A. Young, C. Sayrin, M. Brune, J.-M. Rai- mond, T. Calarco, F. Motzoi, R. Zeier, and E. Cuestas, Fast pulses for high-fidelity circularization of interacting Rydberg atoms (2026)

2026

[2] [4]

S. R. Cohen and J. D. Thompson, PRX Quantum2, 10.1103/prxquantum.2.030322 (2021)

doi:10.1103/prxquantum.2.030322 2021

[3] [6]

Meinert, C

F. Meinert, C. H¨ olzl, M. A. Nebioglu, A. D’Arnese, P. Karl, M. Dressel, and M. Scheffler, Physical Review Research2, 10.1103/physrevresearch.2.023192 (2020)

doi:10.1103/physrevresearch.2.023192 2020

[4] [7]

Facon, E.-K

A. Facon, E.-K. Dietsche, D. Grosso, S. Haroche, J.- M. Raimond, M. Brune, and S. Gleyzes, Nature535, 262–265 (2016)

2016

[5] [8]

E. K. Dietsche, A. Larrouy, S. Haroche, J. M. Raimond, M. Brune, and S. Gleyzes, Nature Physics15, 326–329 (2019)

2019

[6] [9]

Nussenzveig, F

P. Nussenzveig, F. Bernardot, M. Brune, J. Hare, J. M. Raimond, S. Haroche, and W. Gawlik, Physical Review A48, 3991–3994 (1993)

1993

[7] [10]

Signoles, E

A. Signoles, E. Dietsche, A. Facon, D. Grosso, S. Haroche, J. Raimond, M. Brune, and S. Gleyzes, Phys- ical Review Letters118, 10.1103/physrevlett.118.253603 (2017)

doi:10.1103/physrevlett.118.253603 2017

[8] [11]

Patsch, D

S. Patsch, D. M. Reich, J.-M. Raimond, M. Brune, S. Gleyzes, and C. P. Koch, Physical Review A97, 10.1103/physreva.97.053418 (2018)

doi:10.1103/physreva.97.053418 2018