{"id":"c66a59f5-fc61-4556-aeb7-cefd4979ae0e","arxiv_id":"1908.01353","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"Vortex states in two coupled discrete Bose-Hubbard rings undergo coherent Josephson-like oscillations between the rings at low interaction, and become self-trapped and fragmented at strong interaction.","lead":"Two tiny rings of ultracold atoms can swap a swirling vortex current back and forth between them, like a pendulum, when the interaction between atoms is weak. This coherent swapping is the same physics as a Josephson junction and could become a building block for atomtronic flux qubits.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Mean-field equations (11)–(12) have an interaction term with an extra factor of ring size M; the correct two-mode reduction of Hamiltonian (1) gives (Λ/2)(1+z²), not (MΛ/2)(1+z²).","rationale":"The paper's most valuable evidence is the exact many-body dynamics in Sec. III, which is a legitimate numerical experiment and would survive even if the mean-field model were flawed. However, the abstract and conclusions explicitly assert agreement with a mean-field approximation, and the paper uses Λ and Eqs. (12) to define the Josephson and self-trapping regimes. A direct calculation of the interaction energy for ansatz (10) with Hamiltonian (1) yields a factor of 1/M relative to Eq. (11): the vortex modes are delocalized over M sites, so the on-site interaction energy per particle is U(N−1)/(2M) for a fully imbalanced state, not U(N−1)/2. The manuscript's MΛ/2 term effectively treats each ring as a single site. This changes the predicted Rabi frequency by a factor √((1+MΛ)/(1+Λ)) and shifts the self-trapping boundary by M, so the displayed mean-field comparison cannot validate the equations as written. The concern is concrete and easily checked by recomputing the expectation value and rerunning the comparison with the corrected coefficient. I therefore keep the CONDITIONAL verdict: the qualitative picture may still hold, but the quantitative mean-field support and Λ-based regime boundaries require correction. I only partially agree with the reader's identified weakest assumption, because the short-time validity of the two-mode ansatz is a separate limitation; the more elementary problem is that the closed two-mode equations themselves are not the correct reduction of the stated Hamiltonian.","tokens_in":13188,"tokens_out":16959,"duration_ms":177643,"concrete_test":"Recompute the expectation value ⟨H⟩ by substituting ansatz (10) into Hamiltonian (1), using ⟨a†_{l,j}a_{l,j}⟩ = N|α_{l,j}|² and ⟨n_{l,j}(n_{l,j}−1)⟩ = N(N−1)|α_{l,j}|⁴ with |α_{l,↑}|² = (1+z)/(2M) and |α_{l,↓}|² = (1−z)/(2M). Verify whether the interaction coefficient is Λ/2 or MΛ/2. Then rerun the mean-field curves of Figs. 2 and 3 with dφ/dt̃ = Λ z + z/√(1−z²) cosφ, keeping Λ = (N−1)U/(2MJ⊥), and compare to the exact many-body data, especially the N=24, Λ=0.2 panels where the period difference between the two equations is about 13% (0.913 t_R vs 0.791 t_R). If the corrected curves no longer track the exact simulations, the paper's claimed quantitative agreement fails.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim's quantitative support rests on the Josephson-type mean-field equations (12), derived from ansatz (10) and Hamiltonian (1). Evaluating the interaction term explicitly for N bosons in the single-particle orbital w† = c↑Ψ†_{q,↑} + c↓Ψ†_{q,↓}, with c↑ = √((1+z)/2)e^{iφ/2} and c↓ = √((1−z)/2)e^{−iφ/2}, gives per-site amplitudes |α_{l,j}|² = (1±z)/(2M) and hence ⟨H⟩/(NJ⊥) = ϵ_q/J⊥ − √(1−z²) cosφ + [U(N−1)/(4MJ⊥)](1+z²). Using Λ = (N−1)U/(2MJ⊥), the interaction term is (Λ/2)(1+z²), not the (MΛ/2)(1+z²) appearing in Eq. (11). Correspondingly, the second of Eqs. (12) should contain Λ z, not MΛ z. This is not a cosmetic typo: it changes the linearized oscillation frequency from √(1+Λ) to √(1+MΛ) and shifts the predicted self-trapping threshold by a factor M. Because the abstract claims the vortex-flux transfer is 'in good agreement with a mean-field approximation', the quantitative agreement shown in Figs. 2 and 3 is not established until this coefficient is corrected and the comparison recomputed. The exact many-body simulations may still support the qualitative short-time coherent transfer, but the specific mean-field validation is internally inconsistent as written.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the quantum dynamics of vortices in two linearly coupled Bose-Hubbard rings. The authors prepare either a population-imbalanced vortex in one ring or a population-balanced fractional vortex state, then quench on the inter-ring tunneling and on-site interactions, and evolve the many-body state exactly using the SciPy implementation of the Al-Mohy-Higham exponential algorithm. For low interactions they observe coherent, Rabi-period oscillations of the population imbalance and of the chiral current, which they interpret as Josephson-like vortex-flux transfer. They introduce a two-mode mean-field ansatz leading to a short Josephson junction model and compare its predictions with the exact many-body dynamics across the Josephson and self-trapping regimes. The paper concludes that low interactions allow coherent vortex tunneling relevant for flux-qubit-inspired atomtronic devices, while strong interactions suppress chiral currents and lead to fragmentation.","tokens_in":13577,"tokens_out":10018,"duration_ms":96891,"significance":"If the central claim holds, the results are a useful step toward coherent vortex transfer in atomtronic circuits and connect few-site Bose-Hubbard dynamics with flux-qubit phenomenology. The exact many-body calculations are carefully specified and reproducible: the Hamiltonian, parameters, and numerical method are fully stated, and the figures show honest comparisons between exact and mean-field results. The qualitative short-time coherent oscillations seen in the exact dynamics are independent of the mean-field derivation and are supported by the data. However, the quantitative agreement with the mean-field approximation that the abstract claims is undermined by a factor-of-M error in the mean-field energy and equations of motion, so the main quantitative support requires correction and recomputation.","major_comments":[{"comment":"The interaction term in the mean-field energy has an extra factor of M. For the two-mode ansatz (10), the per-site amplitudes are |α_{l,j}|^2 = (1±z)/(2M), so the interaction energy per particle in units of NJ⊥ is U(N−1)/(4M J⊥)(1+z²) = (Λ/2)(1+z²) with Λ defined in Eq. (14). The manuscript's Eq. (11) instead contains (MΛ/2)(1+z²), which is independent of M and unphysical for fixed N and U. Consequently the second of Eqs. (12) should be dφ/d˜t = Λ z + (z/√(1−z²)) cosφ, not MΛ z + (z/√(1−z²)) cosφ. This error changes the linearized oscillation frequency from √(1+Λ) to √(1+MΛ) and shifts the self-trapping threshold by a factor M. Because the dashed mean-field curves in Figs. 2 and 3(a) are obtained from the incorrect Eqs. (12), the quantitative claim of 'good agreement with a mean-field approximation' in the abstract is not established as written. Please correct Eqs. (11) and (12), recompute the mean-field curves, and revisit the corresponding text; the exact many-body results may still support the qualitative short-time coherent transfer, but the current comparison is internally inconsistent.","section":"Section II.2, Eqs. (11)-(12)"}],"minor_comments":[{"comment":"The sentence 'Note that (12) are the typical equations of a single, short bosonic Josephson junction' is somewhat imprecise: the ansatz (10) is a single-orbital two-mode state, which is a stronger assumption than a purely two-site junction. Please rephrase to clarify the level of approximation.","section":"Section II.2"},{"comment":"The phrase 'dually follows' is unclear; consider 'correspondingly follows' or 'in a dual manner'.","section":"Abstract"},{"comment":"The citation for L. Amico, Scientific Reports 86, 153 (2014) appears to have an incorrect volume number; Scientific Reports volumes are numbered consecutively from 4 onward.","section":"Reference [18]"},{"comment":"The many-body transition probability P(t)=[sin(J⊥t/ℏ)]^{2N} is correct, but the difference from the single-particle result in Eq. (18) is a finite-N effect worth a brief parenthetical remark.","section":"Equation (24)"},{"comment":"The sentence 'the many-body dynamics approaches more to the mean field solution for increasing number of particles' is ungrammatical; suggest 'approaches the mean-field solution more closely'.","section":"Section III"},{"comment":"The legend uses (q,q′) notation, while the text refers to 'half-vortex states (1,0)' and 'vortex states (0,2)'; please define the notation in the caption.","section":"Figure 6 caption"}],"recommendation":"major_revision","confidential_remarks":"The factor-of-M error in Eqs. (11)–(12) is a genuine technical flaw in the mean-field comparison, but because the exact many-body results are independent and the error is straightforward to correct, I do not recommend rejection. The authors should recompute all mean-field curves and re-evaluate the abstract's phrasing about quantitative agreement. I would also suggest checking that the corrected mean-field curves still agree with the exact numerics before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The reader's take is close, but I'd put the stress-test note front and center: it's correct. Equation (11) has an interaction term MΛ/2 (1+z²), and the correct two-mode reduction of Hamiltonian (1) gives Λ/2 (1+z²). The second of Eqs. (12) should have Λ z, not MΛ z. This is not cosmetic: the linearized oscillation frequency changes from √(1+Λ) to √(1+MΛ), and the self-trapping threshold shifts by a factor M. Since Fig. 2 compares exact many-body curves against the erroneous mean-field equations, the abstract's \"good agreement\" is not established as written. The exact numerics may still match the corrected mean-field at low Λ, but the curves have to be recomputed before we know.\n\nWhat the paper does well: the exact diagonalization study of N bosons on two coupled M-site rings is genuine and useful. For M=3,4, N up to 24, the authors show that coherent vortex transfer survives at low interaction, that strong interaction causes fragmentation and self-trapping, and that revivals appear at longer times. The single-particle Rabi physics is textbook, but the many-body results—especially the one-body density matrix eigenvalues and the short coherence window—are the actual contribution. Parameters are fully specified, and the numerical method is standard.\n\nSofter spots beyond the mean-field error: the coherence window is short even in the friendly case U/J⊥=0.1, with the sinusoidal transfer degrading around t≈10 t_R; the abstract's \"feasible to experimental realization\" overstates that window. No code or data are provided, which is a real gap for a purely numerical paper. The two-mode ansatz's neglect of other quasimomenta is validated only at short times and low interaction, and the paper acknowledges as much.\n\nWho this is for: people working on atomtronic circuits, persistent currents, and Bose-Hubbard ladders. It deserves a serious referee, because the exact many-body data are a useful reference point even if the mean-field comparison needs correction. My recommendation: send it to peer review, but require the corrected Eq. (12), recomputed figures, and a more careful abstract.","headline":"The exact many-body vortex-transfer results are worth a look, but the mean-field comparison has a factor-of-M interaction error that the authors need to fix before the quantitative claims hold.","tokens_in":14054,"tokens_out":2512,"would_cite":false,"duration_ms":27097,"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":"Vortex tunneling between linearly coupled Bose-Hubbard rings follows a Josephson-like sinusoidal current for low interactions.","keywords":["Bose-Hubbard model","double ring","vortex tunneling","chiral current","Josephson effect","fragmentation","atomtronics","persistent currents"],"falsifier":"Prepare a double ring with $M=3$ sites per ring, $N=6$ atoms, $J/J_\\perp=1$ and $U/J_\\perp=0.1$, starting from a fully imbalanced $q=1$ vortex in the upper ring, and measure the chiral current and population imbalance as a function of hold time. The paper predicts $z(t)\\approx\\cos(2\\pi t/t_R)$ and a chiral current following the same cosine for several periods; if instead the imbalance decays within the first period or the chiral current does not oscillate at the expected frequency, the central claim is falsified.","tokens_in":13006,"feed_emoji":"🌀","tokens_out":9094,"duration_ms":81700,"temperature":0.7,"pith_summary":"This paper asks whether a quantized vortex can tunnel coherently between two weakly connected Bose-Hubbard rings, the atomic analogue of a persistent-current switch in a superconducting circuit. It shows that for low on-site interaction strength, a vortex initially localized in one ring oscillates back and forth into the other ring at the inter-ring oscillation period, and the associated chiral (angular-momentum) current follows the same sinusoidal law as the particle current of a Josephson junction. The same coherent transfer occurs for balanced pairs of fractional vortices, one in each ring. If correct, this gives a concrete few-site bosonic circuit whose persistent-current switching resembles a flux qubit and is within reach of current ring-lattice experiments. At high interaction, by contrast, the chiral current is suppressed and the system falls into fragmented condensation.","feed_headline":"Vortex tunneling between rings mirrors Josephson current","feed_subtitle":"Low-interaction double-ring Bose-Hubbard system switches persistent currents coherently, for atomtronic flux qubits.","key_machinery":"The load-bearing object is the two-mode mean-field ansatz of Eq. (10), which writes the $N$-particle state as a coherent superposition of the single-ring vortex modes $|\\Psi_{q,\\uparrow}\\rangle$ and $|\\Psi_{q,\\downarrow}\\rangle$ with imbalance $z$ and relative phase $\\varphi$. This ansatz reduces the Hamiltonian expectation value to that of a short bosonic Josephson junction and yields the equations of motion (12) for $(z,\\varphi)$, whose Josephson and self-trapping regimes organize the dynamics. The observable that carries the vorticity transfer is the chiral current operator $\\hat{L}_{\\mathrm{chi}}=\\hat{L}_\\uparrow-\\hat{L}_\\downarrow$, whose normalized expectation value is shown to oscillate as $\\cos(2J_\\perp t/\\hbar)$ in the low-interaction regime, dual to the population imbalance. The dimensionless interaction parameter $\\Lambda=(N-1)U/(2M J_\\perp)$ and the period $t_R=\\pi\\hbar/J_\\perp$ are the natural scales of the problem.","core_discovery":"At low on-site interaction, the exact many-body dynamics of a double-ring Bose-Hubbard system shows that an initially imbalanced vortex state $|\\Psi_{q,\\uparrow}\\rangle$ performs coherent oscillations into $|\\Psi_{q,\\downarrow}\\rangle$, with population imbalance $z(t)=\\cos(2J_\\perp t/\\hbar)$ and normalized chiral current following the same cosine, matching the two-mode mean-field equations for a bosonic Josephson junction. The same statement holds for population-balanced fractional-vortex states $|\\Psi_{q,\\uparrow}\\rangle\\otimes|\\Psi_{q',\\downarrow}\\rangle$, whose vortex charges are exchanged between rings while the population balance stays near zero or oscillates weakly depending on the energy difference of the charges. The paper claims that this duality between vortex-flux transfer and sinusoidal particle current establishes a coherent tunneling regime connecting current states with chiral symmetry, and that the loss of coherence around $t\\approx 10t_R$ at low interaction, followed by revivals, is the signature of the underlying many-body spectrum. Strong interactions, in contrast, drive the system into a self-trapped, fragmented state.","pith_inferences":["Beyond the paper, the vortex-charge oscillation suggests a mapping of the double ring to a two-level system in which the vortex charge plays the role of the flux degree of freedom; this could be tested by measuring the visibility of the oscillations as a function of $J_\\perp/J$.","The collapse and revival seen around $t\\approx 10t_R$ imply that a single-tone sinusoidal model is only an envelope; one could fit the revival time as a function of $U$ to extract the nonlinear energy spacing and verify the Bose-Hubbard description.","A natural next experiment is to start from the interacting stationary vortex (as in Fig. 3b) and vary the quench strength $J_\\perp$; the paper predicts almost no difference from the non-interacting preparation, which would confirm that the phase profile, not the initial interaction, controls the coherent transfer.","The result for $q=2$ (no azimuthal current) suggests that even current-free edge states can transfer population coherently, which implies the mechanism is tied to the phase winding rather than to the magnitude of the circulating current."],"forward_implications":["A two-ring Bose-Hubbard circuit can act as a coherent switch of persistent currents: the vortex charge moves back and forth between rings at the inter-ring period while population imbalance oscillates sinusoidally.","The chiral current provides a direct, observable readout of vortex transfer, since it is proportional to the imbalance in the low-interaction regime.","Increasing particle number lengthens the coherent Josephson regime, making the proposed dynamics easier to see in experiments with larger condensates.","The similarity between configurations A and B means the coherent transfer is insensitive to whether the initial vortex is fully imbalanced or a balanced pair of fractional vortices, as long as interactions stay low.","For strong interactions, self-trapping and fragmentation set in, so the same device would fail as a coherent current switch; this marks the boundary of the usable parameter regime."],"supporting_citations":[{"why":"Provides the static properties and vortex solutions of the double-ring Bose-Hubbard system that this work builds on for the dynamical quench.","marker":"[33]"},{"why":"Supplies the two-mode bosonic Josephson junction equations and the flux-qubit analogy used for the mean-field comparison.","marker":"[12]"},{"why":"Previous mean-field demonstration of coherent half-vortex oscillations in a two-component condensate, used as comparison for the Josephson-regime results.","marker":"[19]"},{"why":"Provides the algorithm used to compute the exact many-body time evolution of the initial vortex states.","marker":"[34]"},{"why":"Defines the one-body density matrix whose eigenvalues characterize fragmentation of the many-body state.","marker":"[35]"},{"why":"Identifies the Josephson and self-trapping dynamical regimes of the nonlinear two-mode model used to interpret the results.","marker":"[39]"},{"why":"Justifies the macroscopic superposition mean-field ansatz in the large-particle limit.","marker":"[38]"}],"fun_headline_variants":["Vortex flux tunneling mirrors Josephson particle current","Coherent vortex oscillations in coupled Bose-Hubbard rings","Fractional vortices switch persistent currents like flux qubits","Double-ring vortex dynamics follow Josephson sinusoid","Low-interaction vortex tunneling revives chiral currents"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central claim depends on the assumption that the coupled-ring dynamics remains confined to the two vortex modes used in the mean-field ansatz, with no significant scattering into other quasimomenta; the paper itself shows this coherence fades around ten inter-ring periods even at low interaction.","fun_headline_variants_meta":{"raw":{"variants":["Vortex flux tunneling mirrors Josephson particle current","Coherent vortex oscillations in coupled Bose-Hubbard rings","Fractional vortices switch persistent currents like flux qubits","Double-ring vortex dynamics follow Josephson sinusoid","Low-interaction vortex tunneling revives chiral currents"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000313,"raw_usage":{"total_tokens":1755,"prompt_tokens":895,"completion_tokens":860,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":511,"completion_tokens_details":{"reasoning_tokens":785}},"tokens_in":511,"tokens_out":860,"duration_ms":8783,"temperature":1.0,"reasoning_tokens":785,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:15:15.199801+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Prepare a double ring with $M=3$ sites per ring, $N=6$ atoms, $J/J_\\perp=1$ and $U/J_\\perp=0.1$, starting from a fully imbalanced $q=1$ vortex in the upper ring, and measure the chiral current and population imbalance as a function of hold time. The paper predicts $z(t)\\approx\\cos(2\\pi t/t_R)$ and a chiral current following the same cosine for several periods; if instead the imbalance decays within the first period or the chiral current does not oscillate at the expected frequency, the central claim is falsified.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the static properties and vortex solutions of the double-ring Bose-Hubbard system that this work builds on for the dynamical quench."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Previous mean-field demonstration of coherent half-vortex oscillations in a two-component condensate, used as comparison for the Josephson-regime results."},{"cited_title":"Atala, M","cited_arxiv_id":null,"evidence_quote":"Provides the algorithm used to compute the exact many-body time evolution of the initial vortex states."},{"cited_title":"Static properties of two linearly coupled discrete circuits","cited_arxiv_id":"1807.03838","evidence_quote":"Defines the one-body density matrix whose eigenvalues characterize fragmentation of the many-body state."},{"cited_title":"Son and M","cited_arxiv_id":null,"evidence_quote":"Identifies the Josephson and self-trapping dynamical regimes of the nonlinear two-mode model used to interpret the results."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Justifies the macroscopic superposition mean-field ansatz in the large-particle limit."}],"review_version":1}