REVIEW 3 major objections 6 minor 2 cited by
Repulsive barricade pushes tweezer loading past 80%
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · glm-5.2
2026-07-04 17:43 UTC pith:OOOYZ7MJ
load-bearing objection dark-state shelving approaches in Refs. 39 and 41. It's species-agnostic, which is the real selling point: molecules like CaF are stuck at ~35% fill fractions and don't have obvious shelving states, so this fills a genuine gap. the 3 major comments →
Near-deterministic loading of optical tweezer arrays via repulsive barricade potentials
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The paper's central result is that a repulsive optical barricade, formed by overlapping a blue-detuned beam with a standard red-detuned tweezer, can suppress the flux of particles entering an already-occupied trap by over three orders of magnitude. Using 3D Monte Carlo simulations calibrated to experimental parameters for 87Rb and CaF, the authors show that this suppression corresponds to collision-limited lifetimes of about 1.3 seconds in both cases. With cycle times of 75 ms (Rb) and 15 ms (CaF), four cumulative loading cycles would yield filling fractions of 87 percent and 81 percent respectively, approaching the theoretical maxima of 94 percent and 82 percent for four cycles in the limit
What carries the argument
The barricade tweezer potential is the sum of an attractive light shift (red-detuned beam, diffraction-limited waist) and a repulsive light shift (blue-detuned beam, larger waist). The repulsive beam creates a shell-like barrier around the trapping core. Flux suppression is modeled by 3D semi-classical Monte Carlo trajectory simulations of particles drawn from a Maxwell-Boltzmann distribution evolving in the combined potential. An analytical isotropic barrier model (complementary error function plus exponential term) is shown to fit the Monte Carlo results well once the barrier is sufficiently isotropic. The cumulative fill fraction follows a linear recurrence relation accounting for three f
Load-bearing premise
The simulation assumes the repulsive barricade can be ramped on adiabatically without heating or losing the already-trapped particle. If intensity noise, SLM switching artifacts, or non-adiabatic trap-frequency changes during the ramp cause significant heating, the protected particle would be lost before the next loading cycle, and the cumulative fill fraction model would overestimate performance.
What would settle it
Experimentally, if the measured collision-limited lifetime of a protected particle is substantially shorter than 1.3 s — for instance, because the ramp-on of the repulsive beam causes parametric heating or because off-resonant scattering from the blue-detuned light is higher than estimated — then the cumulative fill fractions would fall below the predicted 87 percent and 81 percent, and the advantage of the barricade scheme over existing dark-state protection methods would diminish.
If this is right
- If the barricade scheme works experimentally as simulated, molecule tweezer arrays — currently limited to 35 percent filling — could reach over 80 percent filling, removing a major bottleneck for molecular quantum simulation and computation.
- The scheme is species-independent and requires no internal-state shelving, so it could be applied to any atom or molecule for which both attractive and repulsive optical potentials can be engineered at compatible wavelengths.
- Combined with existing parallel rearrangement techniques, near-unity filling of large defect-free arrays becomes feasible, directly enabling scaling to larger quantum technology platforms.
- The collision-limited lifetime of 1.3 s sets a ceiling on the number of useful loading cycles; extending this lifetime (via higher barriers or lower-density reservoirs) would push filling fractions closer to unity.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript proposes a 'barricade tweezer' scheme for near-deterministic loading of optical tweezer arrays via multiple loading cycles. The core idea is to superimpose a blue-detuned repulsive beam with a red-detuned attractive tweezer to create a potential barrier around already-loaded sites, preventing new particles from entering and causing light-assisted collisions. The authors perform 3D semi-classical Monte Carlo trajectory simulations to compute the flux suppression through the trapping core for 87Rb and CaF, and use an analytical isotropic barrier model (Eq. 8) to fit the results. They report collision-limited lifetimes of ~1.3 s for both species, yielding cumulative fill fractions of 87% (Rb) and 81% (CaF) after four loading cycles. The fill fraction model (Eq. 1) is derived from a recursion relation (Eq. 4) in the Appendix. The scheme is general and does not require long-lived shelving states, distinguishing it from prior dark-state protection proposals.
Significance. The paper addresses a practical bottleneck in tweezer array scaling — the stochastic loading limit — with a conceptually simple and broadly applicable protocol. The Monte Carlo simulations use realistic experimental parameters and the code is publicly available (GitHub), which is a notable strength for reproducibility. The conservative choices for CaF (R_core = r_FWHM, worst-case m_F state) are appropriate and lend credibility to the predictions. The isotropic barrier model (Eq. 8) provides a useful analytical framework. The proposal is timely given the rapid scaling of tweezer arrays to thousands of sites and the difficulty of extending collisional shielding to molecules.
major comments (3)
- The central claim that cumulative fill fractions exceed 80% after four cycles depends on the barricade being ramped on without loss or heating of the already-trapped particle. The manuscript states this will be done 'adiabatically' (Fig. 1b, main text) but provides no analysis of the ramp dynamics: no estimate of the required ramp time, no simulation of the time-dependent potential during ramping, and no discussion of whether the trap frequencies change discontinuously as the repulsive beam is introduced. For typical tweezer parameters this is likely feasible, but the simultaneous power increase of the attractive beam to maintain trap depth creates a time-dependent potential landscape that could excite or eject the trapped particle. A brief analysis — even an order-of-magnitude estimate of the adiabaticity condition (ramp time vs. inverse trap frequency) and a discussion of potential SLM
- switching artifacts or intensity noise — would substantially strengthen the paper, as this is the single most critical unvalidated assumption underpinning the cumulative fill fraction model.
- Section 'CaF', paragraph 4: The manuscript notes that 'MHz-rate photon scattering during imaging and cooling causes a molecule to traverse the full m_F manifold many times while crossing the barrier,' justifying the use of the worst-case m_F state. However, the time scale for m_F mixing relative to the barrier transit time is not quantified. If the molecule scatters photons at MHz rates but the barrier transit time is sub-microsecond, the molecule may not fully explore the m_F manifold during a single crossing attempt. A brief estimate of the ratio of scattering rate to transit rate would make this conservative assumption more convincing.
minor comments (6)
- Eq. (1): The notation η(τ) is used both for the reload probability and appears in the same equation as η (the cumulative fill fraction). Consider renaming the reload probability (e.g., to p_reload) to avoid confusion.
- Appendix, Eq. (5): The statement 'which recovers the saturation limit of η_0 = 0.5 in the limit of t/τ ≫ 1' should be clarified — Eq. (5) gives η(τ) → 1/2 in that limit, but η_0 is defined separately as the single-cycle loading probability. The connection should be stated more precisely.
- Fig. 3: The axis labels in panels (c) and (d) are difficult to read in the rendered version. Consider increasing font sizes or simplifying the layout.
- Section 'CaF': The footnote defining the F = 1⁻ state (footnote 3) uses a superscript that could be confused with a reference. Consider using a different notation (e.g., F = 1_lower).
- The manuscript uses both 'filling fraction' and 'fill fraction' interchangeably. Standardizing on one term would improve consistency.
- Reference [51] (Humphreys et al., CPC 2025) appears to be by co-authors of this manuscript; this relationship should be disclosed if the journal requires it.
Simulated Author's Rebuttal
We thank the referee for their careful reading and positive assessment. The referee raises two substantive points: (1) the need for an analysis of the adiabatic ramp-on dynamics of the barricade potential, and (2) a quantitative justification for the assumption that molecules explore the full m_F manifold during barrier transit. Both points are well-taken, and we will address them in a revised manuscript by adding an order-of-magnitude adiabaticity estimate and a scattering-rate-to-transit-rate comparison for CaF.
read point-by-point responses
-
Referee: Major Comment 1: The central claim depends on the barricade being ramped on without loss or heating. No analysis of ramp dynamics, adiabaticity condition, trap frequency changes, or SLM switching artifacts is provided.
Authors: The referee is correct that the ramp-on dynamics of the barricade potential are not analyzed in the current manuscript, and that this is a critical assumption underpinning the cumulative fill fraction model. We agree that an explicit discussion is needed. In the revised manuscript, we will add an order-of-magnitude adiabaticity analysis. The key consideration is that the trapped particle sits in the ground state of the attractive tweezer (radial trap frequency ~100 kHz for both Rb and CaF with the parameters used), so the adiabaticity condition requires ramp times much longer than ~10 microseconds. SLM switching times are typically ~1 ms or faster with modern devices (as cited in Ref. 42), and intensity ramps over ~1-10 ms are standard experimental practice, comfortably satisfying the adiabaticity condition. We will also note that the trap depth is maintained constant during the ramp (as described in the manuscript), so the trap frequency changes are modest rather than discontinuous. We will add a brief discussion of potential parametric heating from intensity noise on the repulsive beam, noting that the heating rate from off-resonant scattering (already estimated at 3-4 microK/s) dominates over any parametric heating contribution at these ramp times. This will be added as a short subsection or paragraph in the Methods section. revision: partial
-
Referee: Major Comment 2: The m_F mixing argument for CaF is not quantified. If barrier transit time is sub-microsecond but scattering is MHz-rate, the molecule may not fully explore the m_F manifold during a single crossing attempt.
Authors: This is a fair point. We will add a quantitative estimate. For CaF, the photon scattering rate during imaging/cooling is on the order of 1 MHz (consistent with the manuscript's statement). The barrier transit time can be estimated as the barrier width divided by the thermal velocity: for a barrier width of ~1 micron and a thermal velocity of ~0.3 m/s (at 9 microK), the transit time is ~3 microseconds. This gives roughly 3 photon scattering events per transit, which is indeed marginal for full m_F manifold exploration during a single crossing. However, the relevant comparison is not for a single crossing attempt but over the full loading cycle (10 ms of cooling), during which a molecule near the barrier scatters ~10,000 photons and traverses the m_F manifold many times. The worst-case m_F state assumption is therefore conservative in the sense that it represents a lower bound on barrier effectiveness for any individual encounter, and over many encounters the molecule samples all m_F states. We will add this estimate to the CaF section to make the conservatism of the assumption explicit, and we will soften the claim from 'traverse the full m_F manifold many times while crossing the barrier' to a more precise statement about the number of scattering events per transit and the cumulative exploration over the loading cycle. revision: partial
Circularity Check
No significant circularity found; derivation chain is self-contained
full rationale
The paper's derivation chain proceeds in two independent stages: (1) a first-principles optical potential (Eqs. 2-3, vectorial non-paraxial model) feeds a 3D Monte Carlo trajectory simulation that extracts particle flux through the trapping core, from which the collision-limited lifetime τ = 1/flux is obtained; (2) τ is fed into the fill fraction model (Eq. 1), which is independently derived from a linear recurrence relation (Eq. 4) in the Appendix. No step defines its output in terms of itself. The isotropic barrier model (Eq. 8) is fitted to the Monte Carlo data with j₀ and T as free parameters, but it is used only as an analytical interpolation of the simulation results (shown as dashed lines in Fig. 3), not presented as an independent prediction. The fill fraction recursion (Eq. 4) is solved correctly as a linear recurrence with constant coefficients, and the η(τ) term (Eqs. 5-6) is negligible at the operating parameters (t/τ ≈ 0.06 for Rb, 0.01 for CaF), so even its heuristic treatment for molecules is not load-bearing. The only self-citation is Ref [51] (Humphreys, Matthies, Williams — two co-authors of this paper), cited for tensor polarizability computation methods for molecules; this is a methods/tooling citation, not a load-bearing theoretical premise. The central results (τ ≈ 1.3 s, η > 80%) follow from the Monte Carlo simulation and the independently derived fill fraction formula, with no reduction to inputs by construction.
Axiom & Free-Parameter Ledger
free parameters (3)
- j0 (unperturbed flux)
- T (effective temperature)
- R_core (trapping core radius) =
50 nm (Rb), 0.63 µm (CaF)
axioms (3)
- domain assumption Light-assisted collisions cause pairwise ejection, saturating loading at 50% for atoms.
- domain assumption The barricade potential can be ramped on adiabatically without parametric heating.
- ad hoc to paper The same pairwise ejection mechanism operates for molecules, with an additional unknown loss mechanism.
read the original abstract
Optical tweezers are a powerful tool for creating defect-free arrays of atoms and molecules, enabling advances in quantum simulation, computation, and precision metrology. However, the achievable array size is limited by the initial loading fraction, typically $50\,\%$ for atoms and $35\,\%$ for molecules. Here, we propose a general scheme for enabling multiple loading cycles by protecting trapped particles using a repulsive barrier. We show that collision-limited lifetimes of particles in protected tweezers can exceed one second, leading to filling fractions of over $80\%$ after four loading cycles. Combined with existing rearrangement techniques, this approach enables efficient unity filling of tweezer arrays and provides a scalable pathway towards larger quantum technology platforms.
Figures
Forward citations
Cited by 2 Pith papers
-
Deterministic loading of molecular arrays by microwave-assisted collisions
Microwave-assisted collisions with repulsive van der Waals shelving are predicted to load single molecules into tweezers at 87-96% filling fractions, up from the current 30-40%.
-
Deterministic loading of molecular arrays by microwave-assisted collisions
Shelving molecules in rotationally or vibrationally excited states and driving microwave-assisted collisions can suppress loss, eject one molecule deterministically, and yield predicted filling fractions up to 96%.
Reference graph
Works this paper leans on
- [1]
-
[2]
A. M. Kaufman and K.-K. Ni, Nature Physics17, 1324 (2021)
work page 2021
- [3]
-
[4]
D. Barredo, S. de L´ es´ eleuc, V. Lienhard, T. Lahaye, and A. Browaeys, Science354, 1021 (2016)
work page 2016
-
[5]
W. Lee, H. Kim, and J. Ahn, Opt. Express24, 9816 (2016)
work page 2016
-
[6]
M. A. Norcia, A. W. Young, and A. M. Kaufman, Phys. Rev. X8, 041054 (2018)
work page 2018
- [7]
-
[8]
N. C. Jackson, R. K. Hanley, M. Hill, F. Leroux, C. S. Adams, and M. P. A. Jones, SciPost Phys.8, 038 (2020)
work page 2020
- [9]
-
[10]
L. Anderegg, L. W. Cheuk, Y. Bao, S. Burchesky, W. Ketterle, K.-K. Ni, and J. M. Doyle, Science365, 1156 (2019)
work page 2019
-
[11]
Y. Bao, S. S. Yu, L. Anderegg, E. Chae, W. Ketterle, K.-K. Ni, and J. M. Doyle, Science382, 1138 (2023)
work page 2023
-
[12]
C. M. Holland, Y. Lu, and L. W. Cheuk, Phys. Rev. Lett. 131, 053202 (2023)
work page 2023
-
[13]
A. J. Park, L. R. B. Picard, G. E. Patenotte, J. T. Zhang, T. Rosenband, and K.-K. Ni, Phys. Rev. Lett. 131, 183401 (2023)
work page 2023
-
[14]
N. B. Vilas, P. Robichaud, C. Hallas, G. K. Li, L. An- deregg, and J. M. Doyle, Nature628, 282 (2024)
work page 2024
-
[15]
D. K. Ruttley, T. R. Hepworth, A. Guttridge, and S. L. Cornish, Nature637, 827 (2025)
work page 2025
-
[16]
F. Nogrette, H. Labuhn, S. Ravets, D. Barredo, L. B´ eguin, A. Vernier, T. Lahaye, and A. Browaeys, Phys. Rev. X4, 021034 (2014). 6
work page 2014
-
[17]
H. Kim, W. Lee, H. G. Lee, H. Jo, Y. Song, and J. Ahn, Nature Communications7, 13317 (2016)
work page 2016
-
[18]
K.-N. Schymik, V. Lienhard, D. Barredo, P. Scholl, H. Williams, A. Browaeys, and T. Lahaye, Phys. Rev. A102, 063107 (2020)
work page 2020
-
[19]
D. Bluvstein, H. Levine, G. Semeghini, T. T. Wang, S. Ebadi, M. Kalinowski, A. Keesling, N. Maskara, H. Pichler, M. Greiner, V. Vuleti´ c, and M. D. Lukin, Nature604, 451 (2022)
work page 2022
-
[20]
C. M. Holland, Y. Lu, and L. W. Cheuk, Science382, 1143 (2023)
work page 2023
-
[21]
H. J. Manetsch, G. Nomura, E. Bataille, X. Lv, K. H. Leung, and M. Endres, Nature647, 60 (2025)
work page 2025
-
[22]
G. Pichard, D. Lim, E. Bloch, J. Vaneecloo, L. Boura- chot, G.-J. Both, G. M´ eriaux, S. Dutartre, R. Hostein, J. Paris, B. Ximenez, A. Signoles, A. Browaeys, T. La- haye, and D. Dreon, Phys. Rev. Appl.22, 024073 (2024)
work page 2024
-
[23]
N.-C. Chiu, E. C. Trapp, J. Guo, M. H. Abobeih, L. M. Stewart, S. Hollerith, P. L. Stroganov, M. Kalinowski, A. A. Geim, S. J. Evered, S. H. Li, X. Lyu, L. M. Peters, D. Bluvstein, T. T. Wang, M. Greiner, V. Vuleti´ c, and M. D. Lukin, Nature646, 1075 (2025)
work page 2025
- [24]
- [25]
-
[26]
T. M. Graham, Y. Song, J. Scott, C. Poole, L. Phut- titarn, P. Jooya, K.and Eichler, X. Jiang, A. Marra, B. Grinkemeyer, M. Kwon, M. Ebert, J. Cherek, M. T. Lichtman, M. Gillette, J. Gilbert, D. Bowman, T. Bal- lance, C. Campbell, E. D. Dahl, O. Crawford, N. S. Blunt, B. Rogers, T. Noel, and M. Saffman, Nature604, 457 (2022)
work page 2022
-
[27]
S. J. Evered, D. Bluvstein, M. Kalinowski, S. Ebadi, T. Manovitz, H. Zhou, S. H. Li, A. A. Geim, T. T. Wang, N. Maskara, H. Levine, G. Semeghini, M. Greiner, V. Vuleti´ c, and M. D. Lukin, Nature622, 268 (2023)
work page 2023
-
[28]
D. Bluvstein, A. A. Geim, S. H. Li, S. J. Evered, J. P. B. Ataides, G. Baranes, A. Gu, T. Manovitz, M. Xu, M. Kalinowski, S. Majidy, C. Kokail, N. Maskara, E. C. Trapp, L. M. Stewart, S. Hollerith, H. Zhou, M. J. Gul- lans, S. F. Yelin, M. Greiner, V. Vuleti´ c, M. Cain, and M. D. Lukin, Nature649, 39 (2026)
work page 2026
-
[29]
I. S. Madjarov, A. Cooper, A. L. Shaw, J. P. Covey, V. Schkolnik, T. H. Yoon, J. R. Williams, and M. En- dres, Phys. Rev. X9, 041052 (2019)
work page 2019
-
[30]
A. W. Young, W. J. Eckner, W. R. Milner, D. Kedar, M. A. Norcia, E. Oelker, N. Schine, J. Ye, and A. M. Kaufman, Nature588, 408 (2020)
work page 2020
-
[31]
N. Schlosser, G. Reymond, and P. Grangier, Phys. Rev. Lett.89, 023005 (2002)
work page 2002
-
[32]
N. Schlosser, G. Reymond, I. Protsenko, and P. Grangier, Nature411, 1024 (2001)
work page 2001
-
[33]
T. Gr¨ unzweig, A. Hilliard, M. McGovern, and M. F. An- dersen, Nat. Phys.6, 951 (2010)
work page 2010
-
[34]
B. J. Lester, N. Luick, A. M. Kaufman, C. M. Reynolds, and C. A. Regal, Phys. Rev. Lett.115, 073003 (2015)
work page 2015
-
[35]
M. M. Aliyu, L. Zhao, X. Q. Quek, K. C. Yellapragada, and H. Loh, Phys. Rev. Res.3, 043059 (2021)
work page 2021
-
[36]
M. O. Brown, T. Thiele, C. Kiehl, T.-W. Hsu, and C. A. Regal, Phys. Rev. X9, 011057 (2019)
work page 2019
-
[37]
J. Ang’ong’a, C. Huang, J. P. Covey, and B. Gadway, Phys. Rev. Res.4, 013240 (2022)
work page 2022
-
[38]
A. Jenkins, J. W. Lis, A. Senoo, W. F. McGrew, and A. M. Kaufman, Phys. Rev. X12, 021027 (2022)
work page 2022
-
[39]
E. F. Walraven, M. R. Tarbutt, and T. Karman, Phys. Rev. Lett.132, 183401 (2024)
work page 2024
-
[40]
I. H. A. Knottnerus, Y. C. Tseng, A. Urech, R. J. C. Spreeuw, and F. Schreck, SciPost Phys.19, 118 (2025)
work page 2025
-
[41]
A. L. Shaw, P. Scholl, R. Finklestein, I. S. Madjarov, B. Grinkemeyer, and M. Endres, Phys. Rev. Lett.130, 193402 (2023)
work page 2023
-
[42]
T. Vanackere, A. Hermans, I. Christen, C. Panuski, M. Dong, M. Zimmermann, H. Raniwala, A. J. Leen- heer, M. Eichenfield, G. Gilbert, and D. Englund, Nature Communications16, 11637 (2025)
work page 2025
-
[43]
C. S. Adams and I. G. Hughes,Optics f2f: From Fourier to Fresnel(Oxford University Press, 2018)
work page 2018
-
[44]
Y. H. Fung, P. Sompet, and M. F. Andersen, Technologies 4, 4 (2016)
work page 2016
-
[45]
A. Fuhrmanek, R. Bourgain, Y. R. P. Sortais, and A. Browaeys, Phys. Rev. A85, 062708 (2012)
work page 2012
- [46]
-
[47]
E. B. Norrgard, D. J. McCarron, M. H. Steinecker, M. R. Tarbutt, and D. DeMille, Phys. Rev. Lett.116, 063004 (2016)
work page 2016
-
[48]
H. J. Williams, S. Truppe, M. Hambach, L. Caldwell, N. J. Fitch, E. A. Hinds, B. E. Sauer, and M. R. Tarbutt, New J. of Phys.19, 113035 (2017)
work page 2017
-
[49]
L. Anderegg, B. L. Augenbraun, E. Chae, B. Hemmer- ling, N. R. Hutzler, A. Ravi, A. Collopy, J. Ye, W. Ket- terle, and J. M. Doyle, Phys. Rev. Lett.119, 103201 (2017)
work page 2017
-
[50]
A. L. Collopy, S. Ding, Y. Wu, I. A. Finneran, L. An- deregg, B. L. Augenbraun, J. M. Doyle, and J. Ye, Phys. Rev. Lett.121, 213201 (2018)
work page 2018
-
[51]
B. Humphreys, A. J. Matthies, and H. J. Williams, Com- puter Physics Communications316, 109813 (2025). Appendix: Fill F raction Model We derive the cumulative fill fractionηafterNload- ing cycles by considering the probabilityP N that a given tweezer site is occupied at the end of theN-th loading cycle. In the first cycle, no barricade is present and thus...
work page 2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.