The Atomic Superfluid Quantum Interference Device with tunable Josephson Junctions

Boyang Liu; Jiatao Tan

REVIEW 3 major objections 4 minor 24 references

The Atomic Superfluid Quantum Interference Device with tunable Josephson Junctions

T0 review · 3 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash

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

Pith's one-line read The paper claims that the critical population bias and the critical time of a ring-shaped atomic superfluid interferometer both respond to rotation with a fixed period, providing practical rotation-sensing observables.

desk verdict The rotation-sensing observables are worth a thought, but Eq. (7) does not follow from the paper's own Lagrangian, so the numeric results are not established. read the letter →

arxiv 2411.17211 v1 pith:C67CSW27 submitted 2024-11-26 cond-mat.quant-gas

Jiatao Tan , Boyang Liu This is my paper · ORCID

classification cond-mat.quant-gas

keywords atomtronicsASQUIDJosephsonjunctionsrotationsensingcriticalpopulationbiasself-trappingtunnelingHamiltonianringBose-Einsteincondensate

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

This paper develops an analytical theory for an atomic superfluid quantum interference device (ASQUID), a ring-shaped Bose-Einstein condensate divided by two weak links, and uses it to find observables for rotation sensing. The central result is that the critical population bias $Z_c$, the initial imbalance at which the system switches between self-trapped and Josephson-oscillating behavior, oscillates periodically with the rotation angular velocity, with adjacent peaks separated by $\Omega_0/2$. The same period appears in a second observable, the critical time $t_c$, when the junctions are ramped up adiabatically. A reader should care because an ASQUID is an isolated neutral-atom system in which the usual DC-SQUID current readout is hard to observe; these two quantities give practical readouts tied directly to the fundamental rotation scale $\Omega_0$.

What carries the argument

The load-bearing object is the time-dependent plane-wave ansatz of Eq. (5), $\psi_i(\theta,t)=\sqrt{n_i(t)}e^{i(m_i(t)\theta+\alpha_i(t))}$, which assumes that turning on the weak links only lets the density, winding number, and phase of each half-ring condensate evolve without changing the functional form of the wave function. Substituted into the Lagrangian for a uniformly rotating ring and varied with respect to the six collective variables, it produces the coupled equations of motion (7). These equations generate both the critical population bias and the critical time; in the limit of very weak symmetric links and frozen imbalance they reduce to $I=I_c\sin\Phi$ with $I_c\propto\cos(2\pi\Omega/\Omega_0)$, the superfluid interference current relation that is the direct analogue of the conventional SQUID.

What would settle it

A direct check is to simulate the same ring and barriers with the full one-dimensional Gross-Pitaevskii equation and locate the self-trapping-to-Josephson-oscillation boundary as a function of angular velocity; if the peak spacing of that boundary is not $\Omega_0/2$, or if the boundary shifts with ramp history, the ansatz-based reduction is wrong. Experimentally, one could measure $Z_c$ for initial phase near $\pi/2$ over $\Omega$ from $0$ to $2\Omega_0$ in a ring Bose-Einstein condensate and look for the predicted periodic modulation.

Watch

Extended reading notes

Core claim

The paper claims that in an isolated ring-trap ASQUID with two Josephson junctions, rotation imprints itself in the threshold between dynamical regimes rather than in a steady current. Starting from the tunneling Hamiltonian and a plane-wave ansatz for each half-ring condensate, the authors derive a closed set of equations of motion for the population imbalance, relative phase, and winding numbers. Solving these, they find that the critical population bias $Z_c$ is periodically modulated by the angular velocity with period $\Omega_0/2$, and that the same period appears in the critical time $t_c$ defined by an adiabatic ramp of the junction strengths. They also show that symmetric junctions preserve the clean two-peak structure better than asymmetric ones, that initial relative phases near $\pi/2$ or $3\pi/2$ give the most sensitive response, and that the model reduces to the familiar sinusoidal superfluid interference current in the weak-link limit.

Load-bearing premise

The load-bearing premise is the ansatz of Eq. (5): when the weak links are switched on, the wave function in each half-ring keeps the plane-wave form $\sqrt{n_i(t)}e^{i(m_i(t)\theta+\alpha_i(t))}$, so six collective variables capture all the dynamics; if the condensate deforms, nucleates vortices, or excites higher modes, the equations of motion and the predicted $\Omega_0/2$ periodicity fail.

Editorial extensions

If this is right

Symmetric junctions preserve the two adjacent peaks separated by $\Omega_0/2$, while junction asymmetry blurs them and lowers the rotation-sensing sensitivity.
The sensitivity is set by the fundamental rotation rate $\Omega_0=h/(m\pi r_0^2)$, so increasing the ring radius $r_0$ improves the detectable angular-velocity scale.
Initial relative phases near $\pi/2$ or $3\pi/2$ give clear peak structure across the full rotation range, whereas phases near $0$ or $\pi$ produce flat, insensitive regions.
When the junction strengths are ramped adiabatically, the critical time $t_c$ shows the same $\Omega_0/2$ periodicity and can be found in a single experimental run, unlike $Z_c$ which requires many runs with different initial imbalances.
The model reduces to the standard superfluid interference current $I=I_c\sin\Phi$ in the weak symmetric limit, connecting the ASQUID rotation response to the established SQUID analogue.

Reading between the lines

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

Beyond the paper's symmetric-versus-asymmetric comparison, the gap that opens at $\Omega=(n+0.75)\Omega_0$ when one junction is weaker suggests that deliberately detuning the two barriers could isolate a single junction's phase response, a testable route to local junction diagnostics.
Because $\Omega_0$ scales inversely with trap area, the predicted periodicity could double as an in-situ calibration of the ring radius or the atomic mass, not only as a rotation readout.
The critical-time protocol is the more promising candidate for a practical atomtronic gyroscope, since a single adiabatic run returns a rotation value, whereas finding $Z_c$ requires repeated preparation of different initial conditions.

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

I read the Tan–Liu ASQUID paper. The new idea is to use the critical population bias Zc and the critical time tc as rotation readouts for an isolated two-junction ring BEC, and the claim that both show periodic modulation with period Ω0/2. That is a sensible extension of two-mode Josephson physics, and the paper does a nice job of distinguishing the isolated-superfluid case from a conventional DC SQUID. The comparison to the Ryu et al. experiment for Ω0 (15% agreement) is a useful sanity check.

The problem is Eq. (7). The paper says "straightforward calculation yields" those equations, but they are not the Euler–Lagrange equations of the ansatz in Eq. (5). Varying the Lagrangian with respect to m1 and m2 gives two algebraic constraints. In the paper's own variables, the correct constraints are

m+ + Z m− = 4Ω/Ω0 + (8√(1−Z²)K2/(ℏΩ0)) sin(Φ + π m+/2) + 2π Ż/Ω0, m− + Z m+ = 4ΩZ/Ω0.

Solving these for m± gives the K2 term divided by √(1−Z²); the printed Eq. (7) has it multiplied by √(1−Z²). Substituting the printed expressions back into the first constraint leaves a residual proportional to K2 Z² sin(...). So the dynamical system actually integrated in Figs. 2–5 is not the model's dynamics. Every quantitative result—the Ω0/2 period of Zc, the critical-time modulation, the asymmetric-junction gaps—is therefore unsupported.

I don't think this is a matter of interpretation. The paper claims a derivation and provides no intermediate steps, which is exactly why this slipped through. The operational definition of Zc is also fuzzy: it is picked from a threshold between Z(0)=0.14 and 0.139, with no tolerance specified. That is minor next to the equation bug, but it will need tightening.

Who should care: atomtronics experimentalists interested in rotation sensing will find the proposal intuitive, but they should not trust the numerics until the equations are corrected and the scans re-run. The paper is not a desk reject: the idea is new and the framework is close to something usable. Send it to referees, but the referee should be asked to verify Eq. (7) from first principles. I would expect the authors to fix the algebraic error and re-examine whether the periodic modulation survives.

Referee Report

3 major / 4 minor

Summary. The paper develops a mean-field model for an atomic superfluid quantum interference device (ASQUID) with two tunable Josephson junctions, using a six-parameter ansatz for the condensate wave functions in the two half-rings. From the resulting equations of motion, the authors numerically compute the critical population bias Zc that separates self-trapping from Josephson oscillations, and find that Zc is periodically modulated by the rotation rate Ω with period Ω0/2. They also study asymmetric junctions and time-dependent junctions, introducing a second observable, critical time tc, with the same period. The paper proposes Zc and tc as practical rotation sensors and argues that symmetric junctions are preferable to asymmetric ones.

Significance. If the central results were correct, the paper would provide a concrete and potentially useful scheme for rotation sensing in ring-geometry Bose–Einstein condensates, with two experimentally accessible observables (Zc and tc) and a simple analytical framework. The model is transparent and the numerics are straightforward to reproduce. However, the significance is strongly undercut by three issues: the equations of motion are not derived and appear inconsistent with the stated action principle; the central observable Zc is defined through a numerical threshold with no specified tolerance; and the only experimental comparison is a single value of Ω0 with about 15% agreement. These issues prevent the current manuscript from establishing the claimed periodic-modulation effect.

major comments (3)

[Sec. II, Eq. (7)] Equation (7) is not the Euler–Lagrange system derived from the ansatz (5) and the Lagrangian (6). Varying the Lagrangian with respect to m1 and m2 gives two algebraic constraints for m±. In the notation of Eq. (7), these constraints are m+ + Z m− = 4Ω/Ω0 + [8K2/(ℏΩ0)]√(1−Z²) sin(Φ + πm+/2) + 2πŻ/Ω0 and m− + Z m+ = 4ΩZ/Ω0. Substituting the printed expressions for m+ and m− from Eq. (7) into the first constraint leaves a residual 8K2 Z²/(ℏΩ0)√(1−Z²) sin(Φ + πm+/2). Thus the dynamical system actually integrated in Figs. 2–5 is not the model described in Sec. II, and all subsequent predictions for Zc and tc are not established. A corrected derivation of the equations of motion is required before the numerical results can be trusted.
[Sec. III, Fig. 2] The critical population bias Zc is defined operationally by a numerical threshold: Z(0)=0.14 gives self-trapping while Z(0)=0.139 gives Josephson oscillations, and Zc is then quoted as 0.139. The manuscript does not specify the tolerance or the precise criterion used to distinguish the two regimes. Because the periodic modulation in Fig. 3 is the central claim, it must be demonstrated that the period and peak positions are independent of the chosen tolerance rather than artifacts of the numerical definition.
[Sec. III, after Eq. (8)] The only quantitative comparison to experiment is a single measurement of the fundamental rotation rate Ω0, which agrees with Eq. (8) at the 15% level. No experimental or independent numerical comparison is provided for the predicted periodic modulation of Zc or tc. The claim that these quantities can serve as practical rotation sensors therefore remains unvalidated; the modulation is predicted entirely within the same model used to define the observables.

minor comments (4)

[Sec. IV, Fig. 4] Both panels in Fig. 4 are labeled (a); the second panel should be labeled (b). In addition, the text says 'In Fig.4 (a) we set K2 = 0.01E0' but the caption indicates that panel (a) fixes K1, so the text likely should refer to Fig. 4(b).
[References] References [14] and [20] are identical (both are G. D. Pace et al., Phys. Rev. X 12, 041037 (2022)); one of them should be removed or replaced.
[Title] The title contains a typo: 'tuna ble' should be 'tunable'.
[Sec. II, Eq. (7)] The text states that the equations of motion are obtained by 'Straight forward calculation' but provides no intermediate steps. Given the inconsistency noted above, a full derivation (or a supplementary file with the derivation) is essential for verification.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the derivation chain is self-contained, with no fitted parameters renamed as predictions and no load-bearing self-citations.

full rationale

The paper proceeds from an explicit tunneling Hamiltonian and the stated plane-wave ansatz of Eq. (5) to the equations of motion (7) via the Euler–Lagrange equations. The subsequent quantities—critical population bias Zc and critical time tc—are obtained by numerical integration of this dynamical system and by operational threshold definitions (trajectory classification for Zc, first zero-crossing for tc). They are not fitted to data, and no parameter is adjusted to produce the Ω0/2 periodicity. The only external comparison is the calculation of Ω0 = 62.6 Hz versus the experimental value ≈72 Hz from Ref. [8], which is an independent benchmark rather than an input. There are no self-citations; the references are to experiments and standard results. The ansatz is explicitly presented as an assumption for weak links, not smuggled in through a citation. The conventional-SQUID current I = Ic sin Φ with Ic ∝ cos(2πΩ/Ω0) is derived from the model and checked against a known result, which is a consistency test rather than circular reasoning. The skeptic's claim that Eq. (7) violates the Euler–Lagrange equations would be a correctness defect if true, but it is not a circularity: the predictions would be invalid, not equivalent to the inputs by construction. Therefore, under the specified criteria, there is no identifiable circular step.

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

The central claim rests on a collective-variable ansatz and the weak-link approximation. These are standard modeling choices in two-mode BEC physics, but the ansatz is not derived from a more microscopic picture. No new entities, forces, or fitted parameters are introduced.

free parameters (1)

Zc detection tolerance = not specified explicitly (between 0.001 and 0.002 in Fig. 2)
The critical population bias is identified as the turning point between self-trapping and Josephson oscillation, but the precision threshold for this turning point is arbitrary. Changing the tolerance shifts the reported Zc value, though it does not affect the claimed period Ω0/2.

assumptions (5)

domain assumption Collective-coordinate plane-wave ansatz
Eq. (5) assumes the superfluid in each region stays in a plane-wave state with time-dependent n_i, m_i, α_i, reducing the field theory to six collective variables.
standard math Quantization of phase winding
The order parameter must be single-valued around the ring, giving integer winding numbers and the quantum of circulation used to define Ω0.
domain assumption Weak-link regime
The ansatz and the perturbative treatment of tunneling require K_i small compared to the intra-region energy; the paper states 'we only consider the case of weak link between the two condensates'.
domain assumption Uniform rotation and rotating frame
The rotating-frame Hamiltonian H' = H - ΩJz describes a uniformly rotating system with constant Ω, simulated experimentally by moving barriers.
domain assumption Adiabatic ramping
Sec. V assumes K1=K2=αt with small α so the evolution is adiabatic and the system tracks the instantaneous self-trapping boundary.

how reviews work

0 comments

Cite this review

Pith. "Pith review of The Atomic Superfluid Quantum Interference Device with tunable Josephson Junctions." pith.science (2026). https://pith.science/paper/C67CSW27

@misc{pith2026241117211,
  author       = {Pith},
  title        = {Pith review of: The Atomic Superfluid Quantum Interference Device with tunable Josephson Junctions},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/C67CSW27}},
  note         = {Machine review of arXiv:2411.17211}
}

read the original abstract

The atomic superfluid quantum interference device (ASQUID) with tunable Josephson junctions is theoretically investigated. ASQUID is a device that can be used for the detection of rotation. In this work we establish an analytical theory for the ASQUID using the tunneling Hamiltonian method and find two physical quantities that can be used for the rotation sensing. The first one is the critical population bias, which characterizes the transition between the self-trapping and the Josephson oscillation regimes and demonstrates a periodic modulation behavior due to the rotation of the system. We discuss the variation of critical population bias when the tunneling strengths of the junctions are tuned in different values, and find that the symmetric junctions are better choice than the asymmetric ones in terms of rotation sensing. Furthermore, how the initial phase difference between the two condensates affects the measurement of rotation is also discussed. Finally, we investigate the case of time-dependent junctions and find there is another physical quantity, named critical time, that can be used to detect the rotation.

Figures

Figures reproduced from arXiv: 2411.17211 by the authors.

**Figure 2.** FIG. 2: (Color online) The temporal variation of the popu [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗

**Figure 3.** FIG. 3: (Color online) The critical population bias [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗

Figures from the paper (2 more)

**Figure 4.** Figure 4: FIG. 4: (Color online) The critical population bias [PITH_FULL_IMAGE:figures/full_fig_p004_4.png]

**Figure 5.** Figure 5: FIG. 5: (Color online) (a) The temporal variation of the pop [PITH_FULL_IMAGE:figures/full_fig_p005_5.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

24 extracted references · 19 canonical work pages

[1]

Amico, M

L. Amico, M. Boshier, G. Birkl, et al, A VS Quantum Sci. 3, 039201 (2021)

work page 2021
[2]

To compare their measurement with our calculations the measured value must be multiplied by 2. After do- 4 0.0 0.5 1.0 1.5 2.0 -0.1 0.0 0.1 0.2 0.3 Ω/Ω0 Zc Φ(0) = π Φ(0) = 3π/4 Φ(0) = π/2 Φ(0) = π/4 Φ(0) = 0 0.0 0.5 1.0 1.5 2.0 -0.1 0.0 0.1 0.2 0.3 Ω/Ω0 Zc Φ(0) = π Φ(0) = 3π/4,5π/4 Φ(0) = π/2, 3π/2 Φ(0) = π/4,7π/4 Φ(0) = 0, 2π FIG. 3: (Color online) The c...

work page
[3]

Amico, D

L. Amico, D. Anderson, M. Boshier, J.-P. Brantut, L.- C. Kwek, A. Minguzzi, and W. von Klitzing, Rev. Mod. Phys. 94, 041001 (2022)

work page 2022
[4]

Clarke and A

J. Clarke and A. I. Braginski, The SQUID Handbook (Wiley-VCH, Weinheim, 2004), Vols. 1–2

work page 2004
[5]

K. C. Wright, R. B. Blakestad, C. J. Lobb, W. D. Phillips, and G. K. Campbell, Phys. Rev. Lett. 110, 025302 (2013)

2013
[6]

Ryu, P.W

C. Ryu, P.W. Blackburn, A. A. Blinova, and M. G. Boshier, Phys. Rev. Lett. 111, 205301 (2013)

work page 2013
[7]

Jendrzejewski, S

F. Jendrzejewski, S. Eckel, N. Murray, C. Lanier, M. Ed- wards, C. J. Lobb, and G. K. Campbell, Phys. Rev. Lett. 113, 045305 (2014)

2014
[8]

Eckel, J

S. Eckel, J. G. Lee, F. Jendrzejewski, N. Murray, C. W. Clark, C. J. Lobb, W. D. Phillips, M. Edwards, and G. K. Campbell, Nature 506, 200 (2014)

2014

Show all 24 references

[9]

C. Ryu, E. C. Samson, and M. G. Boshier, Nature Com- munications 11, 3338 (2020)

2020
[10]

Gupta, K

S. Gupta, K. W. Murch, K. L. Moore, T. P. Purdy, and D. M. Stamper-Kurn, Phys. Rev. Lett. 95, 143201 (2005)

2005
[11]

A. S. Arnold, C. S. Garvie, and E. Riis, Phys. Rev. A 73, 041606 (2006)

2006
[12]

C. Ryu, M. F. Andersen, P. Clad´ e, V. Natarajan, K. Helmerson, and W. D. Phillips, Phys. Rev. Lett. 99, 6 260401 (2007)

2007
[13]

Ramanathan, K

A. Ramanathan, K. Wright, S. R. Muniz, M. Zelan, W. Hill III, C. Lobb, K. Helmerson, W. Phillips, and G. Campbell, Phys. Rev. Lett. 106, 130401 (2011)

2011
[14]

Y. Cai, D. G. Allman, P. Sabharwal, and K. C. Wright, Phys. Rev. Lett. 128, 150401 (2022)

2022
[15]

G. D. Pace, K. Xhani, A. Muzi Falconi, M. Fedrizzi, N. Grani, D. H. Rajkov, M. Inguscio, F. Scazza, W. J. Kwon, and G. Roati,Phys. Rev. X 12, 041037 (2022)

2022
[16]

Smerzi, S

A. Smerzi, S. Fantoni, S. Giovanazzi, and S. R. Shenoy Phys. Rev. Lett. 79, 4950 (1997)

1997
[17]

Zheng and J

Y. Zheng and J. Javanainen, Classical and Quantum Models for Phase Imprinting, Phys. Rev. A 67, 035602 (2003)

2003
[18]

Yefsah, A

T. Yefsah, A. T. Sommer, M. J. Ku, L. W. Cheuk, W. Ji, W. S. Bakr, and M. W. Zwierlein, Nature (London) 499, 426 (2013)

2013
[19]

Kumar, R

A. Kumar, R. Dubessy, T. Badr, C. De Rossi, M. de G¨ oer de Herve, L. Longchambon, and H. Perrin, Phys. Rev. A 97, 043615 (2018)

2018
[20]

Luick, L

N. Luick, L. Sobirey, M. Bohlen, V. P. Singh, L. Mathey, T. Lompe, and H. Moritz, Science 369, 89 (2020)

2020
[21]

Del Pace, K

G. Del Pace, K. Xhani, A. Muzi Falconi, M. Fedrizzi, N. Grani, D. Hernandez Rajkov, M. Inguscio, F. Scazza, W. J. Kwon, and G. Roati,Phys. Rev. X 12, 041037 (2022)

2022
[22]

J. F. Annett, Superconductivity, Superﬂuids andConden- sates, Oxford Master Series in Condensed MatterPhysics, Oxford, (2004)

2004
[23]

R. W. Simmonds, A. Marchenkov, E. Hoskinson, J. C. Davis, and R. E. Packard, Nature 412, 55 (2001)

2001
[24]

R. P. Feynman, R. B. Leighton, and M. Sands, The Feynman lectures on physics. 3: Quantum mechanics , Addison-Wesley Publishing Co., Inc., Reading, Mass, London (1965)

1965

Pith tools

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

[1] [1]

Amico, M

L. Amico, M. Boshier, G. Birkl, et al, A VS Quantum Sci. 3, 039201 (2021)

work page 2021

[2] [2]

To compare their measurement with our calculations the measured value must be multiplied by 2. After do- 4 0.0 0.5 1.0 1.5 2.0 -0.1 0.0 0.1 0.2 0.3 Ω/Ω0 Zc Φ(0) = π Φ(0) = 3π/4 Φ(0) = π/2 Φ(0) = π/4 Φ(0) = 0 0.0 0.5 1.0 1.5 2.0 -0.1 0.0 0.1 0.2 0.3 Ω/Ω0 Zc Φ(0) = π Φ(0) = 3π/4,5π/4 Φ(0) = π/2, 3π/2 Φ(0) = π/4,7π/4 Φ(0) = 0, 2π FIG. 3: (Color online) The c...

work page

[3] [3]

Amico, D

L. Amico, D. Anderson, M. Boshier, J.-P. Brantut, L.- C. Kwek, A. Minguzzi, and W. von Klitzing, Rev. Mod. Phys. 94, 041001 (2022)

work page 2022

[4] [4]

Clarke and A

J. Clarke and A. I. Braginski, The SQUID Handbook (Wiley-VCH, Weinheim, 2004), Vols. 1–2

work page 2004

[5] [5]

K. C. Wright, R. B. Blakestad, C. J. Lobb, W. D. Phillips, and G. K. Campbell, Phys. Rev. Lett. 110, 025302 (2013)

2013

[6] [6]

Ryu, P.W

C. Ryu, P.W. Blackburn, A. A. Blinova, and M. G. Boshier, Phys. Rev. Lett. 111, 205301 (2013)

work page 2013

[7] [7]

Jendrzejewski, S

F. Jendrzejewski, S. Eckel, N. Murray, C. Lanier, M. Ed- wards, C. J. Lobb, and G. K. Campbell, Phys. Rev. Lett. 113, 045305 (2014)

2014

[8] [8]

Eckel, J

S. Eckel, J. G. Lee, F. Jendrzejewski, N. Murray, C. W. Clark, C. J. Lobb, W. D. Phillips, M. Edwards, and G. K. Campbell, Nature 506, 200 (2014)

2014