{"id":"40281616-3e29-4825-be4d-7cadcc20965e","arxiv_id":"1908.02468","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In simulated mergers of two toroidal BECs, the final vorticity is set by the initially more populated ring above an imbalance threshold, while elongated traps support long-lived hybrid vortex complexes.","lead":"The authors simulate the merger of two stacked toroidal Bose-Einstein condensates with different circulation. They find that the more populated ring can drag the other into rotation, and that elongated traps support long-lived hybrid vortex structures.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The paper asserts γ-insensitivity without showing evidence; the P_cr threshold and hybrid stability are dissipation-driven, so a different or spatially non-uniform γ could shift the central predictions.","rationale":"I read the paper in good faith: it presents a plausible 3D numerical study of merging toroidal condensates, with a clear physical mechanism (dissipation-driven vortex drift and annihilation) and novel predictions about final vorticity selection. The central argument is not internally inconsistent, and the 3D treatment is a genuine step beyond prior 2D work. However, the most load-bearing condition for the central claims is that the weakly dissipative GPE with γ=0.03 (Eq. (5)) faithfully models the relaxation. The paper explicitly asserts insensitivity to γ but supplies no evidence; this is a missing-support passage (Section II). Since the threshold P_cr and hybrid lifetimes are determined by the dissipation rate and spatial distribution, an unverified γ assumption directly threatens the quantitative predictions. The reader's weakest_assumption identified exactly this issue, so I agree. The appropriate verdict remains CONDITIONAL: the paper's conclusions are plausible but should not be fully accepted until the γ-sensitivity of P_cr and of hybrid stability is demonstrated. A concrete numerical test—varying γ and introducing spatial dependence at fixed P values around the threshold—would settle whether the concern actually lands.","tokens_in":12628,"tokens_out":3714,"duration_ms":44627,"concrete_test":"Recompute the pancake merger for (m1,m2)=(1,0) at P=0.18, 0.20, 0.22, and 0.24 using γ=0.01, 0.03, and 0.1, and also with a spatially dependent dissipation γ(z)=0.03[1+0.5(n(z)/n_max−1)] to approximate thermal-cloud inhomogeneity. Record the final Lp after t=0.1 s for each case; if P_cr moves outside [0.20,0.22] or any final-state assignment flips relative to the γ=0.03 run, then the threshold claim requires qualification. As a secondary check, repeat the elongated A=0.16 hybrid run with γ=0.01 and γ=0.1 to see whether the hybrid remains long-lived and whether the residual-barrier state remains stationary over 2 s.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claims—the threshold initial imbalance P_cr≈0.21 for (m1,m2)=(1,0) and P_cr≈0.29 for (2,0), and the long-lived/stable 3D hybrids in elongated traps—are obtained from the weakly dissipative Gross-Pitaevskii equation (Eq. (1) / Eq. (5)) with a single spatially uniform dissipation parameter γ=0.03. The paper states in Section II: “we have verified that results reported below do not essentially depend on a specific value of γ≪1,” but no supporting data, figure, or scan over γ is presented. This is a load-bearing assumption because the relaxation mechanism itself—vortex drift to the periphery or central hole, fluxon bending and splitting, and antivortex annihilation—is driven by dissipation. Section III.B explicitly says that the vertical drift velocity of fluxons is determined by γ and that hybrid lifetimes depend on dissipation. If the true dissipation is spatially non-uniform (as expected for a thermal cloud in a trapped condensate) or if the effective γ differs from 0.03, the balance between vortex escape and annihilation that sets P_cr could shift outside the quoted 0.20–0.22 interval. Similarly, the claimed “completely stable” hybrid with a residual barrier (Fig. 8) is only demonstrated for one γ value; without quantifying how pinning competes with dissipation-driven drift, the stability statement is not established. No internal inconsistency was found in the derivations, but the missing γ-sensitivity evidence is the weakest link in the central argument.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the nonlinear merger dynamics of two vertically separated toroidal Bose-Einstein condensates carrying different initial vorticities, modeled by a three-dimensional weakly dissipative Gross-Pitaevskii equation with a gradually eliminated Gaussian barrier. For pancake-shaped traps with vorticities (m1,m2)=(1,0), the authors report a threshold initial population imbalance, Pcr in the interval 0.20<Pcr<0.22, below which the merged condensate relaxes to zero angular momentum and above which it relaxes to a single-charged persistent current; for (2,0) they report Pcr≈0.29 and no final m=1 state. For elongated (prolate) traps with m1=1, m2=-1, the simulations produce long-lived hybrid vortex structures with horizontal Josephson vortices, and the authors claim that a residual nonvanishing barrier pins these fluxons and makes the hybrids completely stable, with tunable angular momentum per particle.","tokens_in":12958,"tokens_out":2940,"duration_ms":34797,"significance":"If the reported effects hold, the paper contributes interesting and potentially testable predictions for persistent-current dynamics in atomtronic devices: a dominant superfluid ring can impose its vorticity on a less populated ring, and the final vorticity is controlled by the initial population imbalance and the trap aspect ratio. The model is standard and physically motivated, the explored regimes are clearly separated in parameter space, and the central claims are falsifiable by experiment and by direct simulation. The comparison with earlier 2D studies of Kelvin-Helmholtz instability is useful. However, the quantitative thresholds and the stability statements rest on numerical evidence that is not fully documented, and the asserted insensitivity to the dissipation parameter is not demonstrated.","major_comments":[{"comment":"The paper states that \"we have verified that results reported below do not essentially depend on a specific value of gamma << 1\" but presents no data, figure, or table supporting this claim. This is load-bearing because the relaxation mechanism itself—vortex drift, fluxon bending and splitting, antivortex annihilation, and the vertical drift of fluxons—is explicitly described in Sections III.A and III.B as driven by gamma. A spatially uniform gamma=0.03 is a particular choice, and the threshold Pcr and the claimed hybrid lifetimes are dynamical outcomes that can in principle shift with gamma. Please provide a scan over gamma (e.g., final Lp versus P for gamma=0.01, 0.03, 0.06) or otherwise quantitatively demonstrate the insensitivity asserted in the text.","section":"Section II, Eq. (1) and text after Eq. (1)"},{"comment":"The critical intervals 0.20<Pcr<0.22 for (m1,m2)=(1,0) and Pcr≈0.29 for (2,0) are determined from very sparse data: Fig. 4 shows only two evolution curves and its inset shows a handful of final points, while Fig. 5 shows four curves. No grid spacing, time step, domain size, or convergence tests are reported anywhere in the paper. Because the thresholds are central quantitative claims, the manuscript needs a statement of the numerical parameters and a demonstration that the threshold location is stable under resolution and time-step refinement, together with additional P values bracketing each threshold.","section":"Section III.A, Figs. 4 and 5"},{"comment":"The claim that a residual barrier makes the hybrid \"completely stable\" is supported only by finite-time simulations at a single value of gamma and without a quantitative stability analysis. The text states that hybrids are transient and will eventually transform into ordinary states, yet later calls the pinned configuration completely stable; the distinction is not established. Please provide either longer-time simulations with a stated duration, a scan over residual barrier height and gamma, or a linear stability analysis around the pinned state before using the word \"stable\" rather than \"long-lived on the simulated timescale.\"","section":"Section III.B, Fig. 8"}],"minor_comments":[{"comment":"There are several typos and inconsistent notations: \"superfl uid\" in the header, \"filed\" in the Introduction, and \"withm = 0\" in the caption of Fig. 3 should be \"with m2 = 0.\"","section":"General/Introduction"},{"comment":"The wave function is denoted psi in Eq. (1) and Eq. (5) but Psi in Eq. (8); please use a single symbol or explicitly define the change. Also, Eq. (10) has a stray period after the closing brace.","section":"Section II, Eqs. (1) and (8)"},{"comment":"The asymmetry parameter P is defined through N1 and N2, but the text discussing \"N1 > N2 for z0 > 0\" would benefit from an explicit statement of how z0 is related to the barrier center and the ring populations, since z0 is also used as a possible shift of the barrier position.","section":"Section II, text after Eq. (10)"},{"comment":"The sentences \"they remain robust in the course of long evolution, and do not decay even being strongly perturbed\" and \"the emerging hybrid states are transient ones, as they will eventually transform into usual ones\" are in tension; please clarify whether hybrids are finite-lifetime or effectively stable on experimental timescales.","section":"Section III.B, paragraph on hybrid decay"},{"comment":"Reference [9] is cited as an arXiv preprint; if it has been published, the published version should be cited. Additionally, the experimental parameters quoted from Refs. [26,27] would be easier to check if the corresponding trap geometries were summarized in one place.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely within scope for a cold-atoms/atomtronics journal and the central scenario is plausible, but the quantitative claims need stronger numerical support. The main concerns are not about the model choice itself but about the missing dissipation-sensitivity evidence and the lack of convergence tests for the reported thresholds. These are fixable within the manuscript's scope, so I recommend major revision rather than rejection. I would also encourage the editor to ask for the numerical resolution parameters and for a more careful use of the word \"stable\" in the elongated-trap section."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What you should know: this paper reports a genuinely new numerical result. It goes beyond the existing 2D counterflow and static fluxon studies by simulating the full 3D dissipative merger of two toroidal BECs with different vorticities. The two headline findings are (i) a threshold in initial population imbalance Pcr ≈ 0.21 for the (m1,m2)=(1,0) case (and ≈0.29 for (2,0)) that determines whether the final condensate has zero or unit circulation, and (ii) long-lived 3D hybrid vortex states in axially elongated traps, which can be made stationary by a residual barrier and whose angular momentum per particle is tunable between integer vorticities. Both are concrete predictions an experimental group could attempt to test.\n\nThe modeling is standard: a weakly dissipative Gross-Pitaevskii equation with uniform γ = 0.03 and a time-dependent barrier. The initial states are prepared cleanly by imaginary-time propagation with imprinted phases. The paper also gives a useful energy argument for why horizontally oriented fluxons are preferred in elongated traps, and it is honest that the hybrids are transients unless pinned.\n\nThe soft spots are real but not fatal. The biggest is the asserted insensitivity to γ. The words “we have verified that results reported below do not essentially depend on a specific value of γ ≪ 1” appear in Section II, but no scan, plot, or data is shown. That matters because the relaxation path—vortex drift, annihilation, fluxon bending—is driven by this dissipation. A spatially nonuniform γ, or a different value, could move Pcr outside the quoted 0.20–0.22 interval. The stress-test note has this right. The fix is simple: show Lp(t) for γ = 0.01, 0.05, 0.1, and preferably a spatially dependent γ test.\n\nSecond, the paper gives no numerical convergence details: grid spacing, time step, domain size, or resolution checks. The Pcr intervals are determined from sparse points on an inset plot. Third, “completely stable” hybrids are demonstrated only for finite times and a single barrier configuration; more evidence is needed to distinguish a long-lived transient from a genuinely stable state. The claim that Lp can be tuned to any value between m1 and m2 rests on a few data points in Fig. 8. The qubit suggestion is speculative, but clearly labeled as such.\n\nAll that said, the central scenario is plausible and internally consistent. The citation pattern is fine; the authors build on their own earlier fluxon work, which is appropriate. This is not a desk-reject. It deserves a serious referee, and the revision should be moderate: add the γ-sensitivity evidence, the numerical methods, and tighten the stability language. I would cite it if I worked on persistent currents in toroidal BECs.","headline":"Solid new 3D simulation results with a real gap: the threshold Pcr and stable hybrids are plausible, but the paper's asserted insensitivity to the dissipation rate γ needs to be shown before publication.","tokens_in":13443,"tokens_out":2714,"would_cite":true,"duration_ms":30570,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["03.75.Lm","67.85.De"],"model":"deepseek-v4-flash","headline":"The final state of two merging superfluid rings is decided by their population imbalance and trap shape.","keywords":["Josephson vortices","Bose-Einstein condensates","toroidal traps","persistent currents","population imbalance","Gross-Pitaevskii equation","Kelvin-Helmholtz instability","hybrid vortex structures"],"falsifier":"Prepare a double-ring condensate with vorticities $(m_1,m_2)=(1,0)$ and a pancake trap, ramp the barrier off, and measure the final angular momentum per particle for imbalances below and above $P=0.21$. Observing a gradual rather than steplike change of $L_p$ with $P$, or finding that the threshold shifts with the dissipation rate $\\gamma$, would contradict the claim. Conversely, a direct measurement of a stable hybrid in an elongated trap with a residual barrier, whose $L_p$ changes when the barrier position is moved, would confirm the central prediction.","tokens_in":12448,"feed_emoji":"🌀","tokens_out":7467,"duration_ms":77019,"temperature":0.7,"pith_summary":"This paper asks what happens when two parallel ring-shaped Bose-Einstein condensates carrying different amounts of circulation are encouraged to merge by lowering the barrier between them. It claims the final state is set by the initial population imbalance between the rings and by the shape of the three-dimensional trap. In flat (pancake) traps, a ring carrying a single-quantum persistent current drags an initially non-rotating ring into the same rotating state only if its population is large enough; below a threshold $P_{\\rm cr}\\approx 0.21$, the merged condensate relaxes to zero angular momentum. In axially elongated traps, the merger instead produces long-lived three-dimensional hybrid vortex structures, which become completely stable when a weak residual barrier pins the Josephson vortices, and their angular momentum per particle can be tuned continuously. Because the dynamics are driven by the splitting, bending, and drift of Josephson vortices, the result connects persistent currents, collective vortex motion, and the fate of superflows after merging.","feed_headline":"Merging superfluid rings: imbalance sets the final spin","feed_subtitle":"3D simulations show a threshold population gap decides whether the merged condensate ends up still or rotating by one quantum.","key_machinery":"The load-bearing object is the Josephson vortex, a rotational fluxon that appears in the low-density barrier between two superflows with different topological charges because the tunneling current has azimuthal periodicity. The argument runs on the three-dimensional weakly dissipative Gross-Pitaevskii equation with a single dissipation constant $\\gamma=0.03$, a toroidal trapping potential whose aspect ratio $A=\\omega_z/\\omega_r$ is either large (pancake) or small (elongated), and a time-dependent sheet barrier that vanishes linearly in time. The fluxon's fate—bending and splitting into vertically oriented vortex-antivortex pairs in pancake traps, or staying radially oriented in elongated traps—is what determines whether the final condensate acquires, keeps, or loses the angular momentum of the initially dominant ring. A geometric energy argument based on the ratio of Thomas-Fermi widths explains why the radial orientation is favored in elongated traps.","core_discovery":"The paper's central claim is that the merger of two toroidal condensates with different vorticities is controlled by the nonlinear three-dimensional dynamics of Josephson vortices (rotational fluxons) that form in the tunneling barrier between the rings. For a pancake-shaped trap with initial vorticities $(m_1,m_2)=(1,0)$, the angular momentum per particle after relaxation is a step function of the initial population imbalance $P$: zero for $P<P_{\\rm cr}$ and one quantum for $P>P_{\\rm cr}$, with $0.20<P_{\\rm cr}<0.22$. For $(m_1,m_2)=(2,0)$ the same logic selects a final state with either $m=0$ or $m=2$, never $m=1$, with $P_{\\rm cr}\\approx 0.29$. In an axially elongated trap with $m_1=1$, $m_2=-1$, the fluxons keep a horizontal (radial) orientation, the Kelvin-Helmholtz instability does not develop, and the long-lived hybrid state—two axially separated parts with different vorticities connected by radially oriented Josephson vortices—emerges. Keeping a nonzero residual barrier pins the fluxons and turns the hybrid into a completely stable stationary vortex complex whose angular momentum per particle can be tuned anywhere in $m_1<L_p<m_2$.","pith_inferences":["If the central claim is right, a direct measurement of the step in $L_p$ at $P\\approx 0.21$ in a pancake trap would be a clean experimental test, and the step location should be independent of ramp speed if the dissipation assumption is correct.","The fluxon-splitting picture suggests that similar merger thresholds should appear in atomtronic circuits whenever two current-carrying loops with different circulation are connected by a tunable junction.","If real dissipation is spatially nonuniform, the threshold $P_{\\rm cr}$ likely becomes geometry-dependent rather than a universal number; varying the temperature of the thermal cloud could expose this.","The completely stable hybrid with two axial vorticities and pinned fluxons is a natural platform for topologically protected information storage, since its two sectors carry distinct integer charges joined by Josephson vortices."],"forward_implications":["In pancake traps, the final angular momentum of a merged $(1,0)$ double ring is a sharply quantized switch: zero below $P_{\\rm cr}\\approx 0.21$, one quantum above it.","A $(2,0)$ input cannot settle at $m=1$; symmetry forces the final state to either $m=0$ or $m=2$, with threshold $P_{\\rm cr}\\approx 0.29$.","In elongated traps with $m_1=1,m_2=-1$, the 3D vortex lattice does not roll up into Kelvin-Helmholtz turbulence, even under noise of a few percent, so long-lived hybrids rather than turbulent decay are the expected outcome.","Retaining a residual barrier pins the Josephson vortices and yields completely stable 3D hybrid vortex complexes with tunable angular momentum per particle over $m_1<L_p<m_2$.","Tuning the barrier position after the hybrid forms changes $L_p$ dynamically, providing an in-situ control knob rather than only an initial-condition effect."],"supporting_citations":[{"why":"Establishes that coupled rings with different vorticities host Josephson vortices with zero net tunneling current, the starting configuration for the merger.","marker":"[9]"},{"why":"Provides the 2D counterflow simulations in which vortex streets undergo Kelvin-Helmholtz roll-up, the baseline the 3D elongated-trap results contrast with.","marker":"[15]"},{"why":"Supplies the weakly dissipative Gross-Pitaevskii framework and the value $\\gamma=0.03$ used for quantum-vortex dynamics.","marker":"[12]"},{"why":"Derives the weakly dissipative Gross-Pitaevskii equation from coupling to a thermal reservoir, the model's foundation.","marker":"[22]"},{"why":"Introduced hybrid states built from persistent currents with different topological charges, which the stabilized 3D hybrids resemble.","marker":"[10]"},{"why":"Provides the pancake-shaped toroidal trap parameters used in the simulations.","marker":"[26]"},{"why":"Reports the experimental condensate lifetime used to set the slow particle-decay time scale $t_0=10$ s.","marker":"[40]"},{"why":"Supports the energy argument that persistent currents of different topological charges are locally stable in toroidal condensates.","marker":"[43]"}],"fun_headline_variants":["Critical population imbalance sets spin after superfluid ring merger","Merging rings: a sharp threshold in matter determines final vortex","Population gap steers vortex takeover in colliding superfluid rings","Hybrid vortex structures from merging persistent-current superflows"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that one small, spatially uniform friction parameter (set to 0.03 in the simulations) faithfully describes how real vortex lines drift and annihilate after the rings merge; the paper asserts without showing evidence that the outcomes do not depend on the exact value.","fun_headline_variants_meta":{"raw":{"variants":["Critical population imbalance sets spin after superfluid ring merger","Merging rings: a sharp threshold in matter determines final vortex","Population gap steers vortex takeover in colliding superfluid rings","Hybrid vortex structures from merging persistent-current superflows"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000478,"raw_usage":{"total_tokens":2404,"prompt_tokens":1017,"completion_tokens":1387,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":633,"completion_tokens_details":{"reasoning_tokens":1319}},"tokens_in":633,"tokens_out":1387,"duration_ms":12470,"temperature":1.0,"reasoning_tokens":1319,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:42:31.481582+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Prepare a double-ring condensate with vorticities $(m_1,m_2)=(1,0)$ and a pancake trap, ramp the barrier off, and measure the final angular momentum per particle for imbalances below and above $P=0.21$. Observing a gradual rather than steplike change of $L_p$ with $P$, or finding that the threshold shifts with the dissipation rate $\\gamma$, would contradict the claim. Conversely, a direct measurement of a stable hybrid in an elongated trap with a residual barrier, whose $L_p$ changes when the barrier position is moved, would confirm the central prediction.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes that coupled rings with different vorticities host Josephson vortices with zero net tunneling current, the starting configuration for the merger."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the 2D counterflow simulations in which vortex streets undergo Kelvin-Helmholtz roll-up, the baseline the 3D elongated-trap results contrast with."},{"cited_title":"Tsubota, M","cited_arxiv_id":null,"evidence_quote":"Supplies the weakly dissipative Gross-Pitaevskii framework and the value $\\gamma=0.03$ used for quantum-vortex dynamics."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Derives the weakly dissipative Gross-Pitaevskii equation from coupling to a thermal reservoir, the model's foundation."},{"cited_title":"Driben, Y","cited_arxiv_id":null,"evidence_quote":"Introduced hybrid states built from persistent currents with different topological charges, which the stabilized 3D hybrids resemble."},{"cited_title":"Carretero-Gonzalez, N","cited_arxiv_id":null,"evidence_quote":"Provides the pancake-shaped toroidal trap parameters used in the simulations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the experimental condensate lifetime used to set the slow particle-decay time scale $t_0=10$ s."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supports the energy argument that persistent currents of different topological charges are locally stable in toroidal condensates."}],"review_version":1}