{"id":"1b6b1083-f234-4c04-b667-b4e777517b1d","arxiv_id":"1908.02351","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"For few bosons on a ring, a rotating barrier reduces the superfluid fraction without significantly changing the condensate fraction, confirming that the two are independent.","lead":"This paper uses a beyond-mean-field quantum simulation method to compute how much a small cloud of bosons on a ring can still flow without friction when a moving barrier is added. The result shows the barrier reduces the superfluid fraction while leaving the condensate fraction almost unchanged, distinguishing two phenomena often conflated.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unverified MCTDHB convergence (M, grid, time step unreported) is the load-bearing gap; the abstract's 'exclusively' also overstates Table I.","rationale":"The paper's only evidence for the central claim is its own MCTDHB code. No machine-checked proof, no released code or data, no external comparison beyond a mention of MCTDH-X examples. The physical picture is plausible—a weak link impedes superflow while weak interactions keep a large condensate—so I do not recommend rejection. But because the quantitative curves in Fig. 5 and Table I are the basis of the claimed independence, the absence of any convergence report makes the reader's CONDITIONAL verdict appropriate. The concern is concrete: the hardest parameter point (gamma=30, lambda=10^4) is exactly where density features are narrowest and correlations strongest, and the Table I caption admits code limitations at gamma=30. A systematic M- and grid-convergence study at representative points would settle it. The abstract's 'exclusively' is internally inconsistent with Table I and should be softened to 'mostly' regardless of the convergence outcome. The reader identified the same primary weakness; I agree.","tokens_in":10623,"tokens_out":17452,"duration_ms":195159,"concrete_test":"Recompute the two extreme data points (N=5, gamma=30, lambda=10^4) and (N=11, gamma=10, lambda=10^3) with MCTDHB using M=2,3,4,5,6 orbitals, grid spacings Delta x = pi R/128, pi R/256, pi R/512, and an imaginary-time step small enough that the energy change per step is below machine precision; then compare <rho_s>_0 from Eq. (15) and the largest eigenvalue of n^(1)(x,x') with the paper's Fig. 5 and Table I values. If either observable changes by more than 1% under refinement, the reported curves are not converged and the independence claim is unverified; if they are stable or the change is below 1%, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is a numerical statement: barrier height drives <rho_s>_0 from 1 to 0 while the condensate fraction remains essentially fixed. For that statement to be supported, the MCTDHB results must be converged in the variational subspace (number of orbitals M, hence Fock configurations), the spatial grid, and the imaginary-time propagation. The manuscript reports none of these. The only convergence assertion is in the Conclusions: 'This method enable us to check convergence for the observables presented here, enlarging the basis...', with no data shown. The Table I caption is more candid: 'For gamma=30 we were able to perform the calculations only for 5 particles due to our code limitations.' This limitation matters precisely at the strong-interaction, high-barrier corner of Fig. 5 (lambda up to 10^4, density at the barrier vanishing), where a narrow barrier and strong correlations require many high-momentum orbitals. If M is too small, the largest eigenvalue of the 1-RDM is biased upward and the superfluid fraction can be quantitatively wrong. A secondary inconsistency: the abstract says the condensation fraction 'depends exclusively on the interaction strength,' but Table I shows it varies with lambda by up to about 7% (N=5, gamma=30: 0.75 to 0.70) and also depends on N. This does not by itself destroy the independence claim, but it shows the headline statement is stronger than the reported data.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies few bosons (N=5, 8, 11) confined to a ring and subject to a rotating barrier (weak link), using the multiconfigurational time-dependent Hartree method for bosons (MCTDHB). The authors derive the mass current from the continuity equation, express the superfluid fraction through the energy derivative with respect to the rotation frequency, and compute it numerically from ground states obtained by imaginary-time propagation. They report that the superfluid fraction at rest decreases from 1 toward 0 as the barrier height increases, while the condensate fraction (largest eigenvalue of the one-body density matrix) is only weakly affected, and they interpret the physics via the one-body correlation function g^(1). The central claim is that barrier-induced suppression of superfluidity is independent of condensation, which would mean a rotating weak link can impede dissipationless flow without destroying the condensate in a few-body system.","tokens_in":10890,"tokens_out":4971,"duration_ms":53661,"significance":"If the numerical results are quantitatively reliable, the paper provides a beyond-mean-field demonstration that in a few-boson ring the superfluid fraction and the condensate fraction respond very differently to a moving barrier, which is directly relevant to persistent-flow and atomtronic experiments. The derivation of the superfluid fraction from the energy derivative is standard and internally consistent, and the authors explicitly check their finite-difference procedure against the momentum formula with a claimed sub-1% deviation. The use of MCTDHB goes beyond the Gross-Pitaevskii level and allows access to correlation functions that are exactly unity in mean-field theory, yielding falsifiable predictions for g^(1). The main weakness is that the variational convergence of the MCTDHB calculations is not documented, and the abstract overstates what Table I actually shows.","major_comments":[{"comment":"The MCTDHB calculations are not documented with the number of single-particle orbitals M, the spatial grid size, or the imaginary-time step, and no convergence data are shown anywhere. The Conclusion states that the method 'enable us to check convergence ... enlarging the basis', but no such check is reported. This is a load-bearing gap because the high-barrier, strong-interaction corner (λ up to 10^4, γ=30) is precisely where a truncated variational basis and a finite grid are most likely to bias both the superfluid fraction and the largest eigenvalue of the one-body density matrix. Please provide convergence tests with increasing M and grid resolution for representative parameter points, or explicitly restrict the quantitative claims to the converged regime.","section":"Section IV, Fig. 5 and Table I"},{"comment":"The statement that the condensation fraction 'depends exclusively on the interaction strength' is stronger than the data in Table I support. For N=5, γ=30 the maximum/minimum condensation fraction changes from 0.75 to 0.70 as λ varies, and the values also depend on N (e.g., 0.9962 for N=5 versus 0.9936 for N=11 at γ=1). The claim should be softened to state that the condensation fraction depends weakly on the barrier height and predominantly on the interaction strength.","section":"Abstract and Conclusions"},{"comment":"The approximation ⟨ρ_s⟩0 ≈ ⟨ρ_s⟩(0.01), obtained by a finite difference between Ω=0 and Ω=0.02, is justified only by the assertion of a less than 1% deviation from Eq. (12); no data, table, or plot is provided to support this check. Please show the comparison between Eqs. (13) and (12) for the parameter points used, and ideally display ⟨ρ_s⟩(Ω) in the small-Ω region for several barrier heights so the extrapolation to Ω=0 can be assessed.","section":"Section IV, paragraph after Eq. (15)"},{"comment":"Because results for γ=30 are available only for N=5 (as acknowledged in the Table I caption), the conclusion that 'the number of particles and strength of interactions have a small impact in the form of the curves of ⟨ρ_s⟩0 as a function of λ' is not supported for strong interactions. Either extend the calculations to N=8 and N=11 at γ=30, or restrict the claim to the parameter range actually computed.","section":"Fig. 5 and Table I"}],"minor_comments":[{"comment":"The manuscript contains numerous typographical errors, including 'tunnable' (Introduction), 'eingenstates' and 'eingenvalue' (Section III), 'in unis of' (Section II), and 'the strenght' (Section IV). A careful proofread is needed.","section":"Throughout"},{"comment":"The crossing points of the single-particle parabolas are denoted Ω̃_j = (j + 1/2) in Section III, but in the Conclusions the same points are called Ω_j = (j + 1/2) and described as 'peaks'; the notation should be made consistent.","section":"Section III and Conclusions"},{"comment":"The y-axis of panel (B) is labeled ⟨ρ_s⟩(Ω), while the caption calls it 'current fraction' and the text refers to 'current fraction' in the discussion after Eq. (15). Please align the terminology used in the main text, equations, and figures.","section":"Fig. 3"},{"comment":"The period τ = 2πR/v is introduced only inside Eq. (11); defining it in the text before the equation would make the expression clearer.","section":"Eq. (11)"},{"comment":"The variational equations resulting from the MCTDHB ansatz are cited to Refs. [30,38] but are not written out; stating the coupled equations for C_α and the orbitals would make the numerical implementation more transparent.","section":"Section II"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope for a quantum-gas journal and the physical question is timely. The main barrier to publication is the missing convergence documentation for MCTDHB, which is essential for the quantitative claims in Fig. 5 and Table I. The abstract's 'exclusively' claim should also be revised. I do not see a fundamental flaw in the derivation of the superfluid fraction; the issues are fixable with additional numerical evidence and a more careful wording of the claims."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper is a straightforward beyond-mean-field numerical study of a few bosons in a ring with a rotating weak link. The new thing is the specific computation: using MCTDHB to extract the superfluid fraction from the energy derivative and to show that a moving barrier can suppress superfluidity while the condensate fraction stays roughly flat. That is a useful quantitative point for atomtronics and mesoscopic ring studies, and the paper does it cleanly enough. The current derivation is standard, the finite-difference check against the momentum formula at Omega=0.01 (reported deviation under 1%) is a good sanity check, and the one-body correlation maps in Fig. 7 are genuinely beyond mean field.\n\nThe main soft spot is convergence. The paper never reports M (number of orbitals), grid size, or imaginary-time step. The conclusions say the method \"enable us to check convergence... enlarging the basis,\" but no convergence data are shown. The Table I caption is more candid: for gamma=30 the calculations were done only for 5 particles \"due to our code limitations.\" That limitation matters most precisely in the strong-interaction, high-barrier corner of Fig. 5, where a narrow barrier and strong correlations demand many high-momentum orbitals. If M is too small, the largest eigenvalue of the one-body density matrix is biased upward and the superfluid fraction can be quantitatively off. This is not a fatal flaw, but it is the load-bearing gap, because the central claim is a numerical statement and the numerical convergence is unverified.\n\nA second, smaller issue: the abstract says the condensation fraction \"depends exclusively on the interaction strength,\" but Table I shows it varies with lambda by up to about 7% (N=5, gamma=30) and also depends on N. The qualitative independence claim survives, but the word \"exclusively\" is stronger than the reported data.\n\nThe citation pattern looks fine: the ring-superfluidity literature (including GP-based studies) and the MCTDHB method are properly cited. The paper ships no code or data, which is a missed opportunity given the convergence questions.\n\nThis paper is for readers working on few-body ring superfluidity and MCTDHB applications. It deserves a serious referee; the missing convergence details are fixable in revision, and the physics is coherent. I would ask the authors to report M and a convergence table, and to soften the abstract.","headline":"Useful MCTDHB study of few-boson ring superfluidity, but the central numerical claim rests on convergence details the paper never reports, and the abstract overstates the condensate-fraction independence.","tokens_in":11416,"tokens_out":1474,"would_cite":false,"duration_ms":16508,"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":"A rotating barrier suppresses superfluidity but leaves condensation intact","keywords":["superfluidity fraction","few bosons","ring geometry","rotating weak link","MCTDHB","Bose-Einstein condensation","one-body correlation function","persistent current"],"falsifier":"Repeat the same parameter sets with an exact or much larger-basis few-boson calculation, increasing the number of orbitals until all observables stop changing, and check whether the superfluid fraction at fixed barrier height shifts by more than the paper's stated accuracy; alternatively, measure the condensation fraction while sweeping barrier height at fixed interaction strength in a cold-atom ring and look for a change larger than the few-percent variation reported in the table.","tokens_in":10392,"feed_emoji":"🌀","tokens_out":6112,"duration_ms":65528,"temperature":0.7,"pith_summary":"This paper asks whether a moving potential barrier (a weak link) can degrade superfluid flow in a ring of just a few interacting bosons, and whether that degradation is tied to the loss of Bose-Einstein condensation. Using a beyond-mean-field variational method, the authors compute ground states for up to 11 particles and extract the superfluid fraction from the energy's response to rotation. They find that raising the barrier height steadily drives the superfluid fraction from 1 toward 0, while the condensation fraction, the largest natural-orbital occupation, depends almost exclusively on interaction strength. The conclusion is that, in this few-body system, the suppression of dissipationless flow and the loss of condensation are independent effects. This matters because it separates two concepts often conflated in cold-atom rings and yields experimentally relevant predictions for weak-link atomtronic devices.","feed_headline":"Rotating weak link kills superfluidity, leaves condensate intact","feed_subtitle":"Few-boson ring calculation separates flow suppression from condensate depletion, with experimental predictions for weak-link devices.","key_machinery":"The load-bearing machinery is the MCTDHB ansatz, in which the many-body wavefunction is a superposition of all Fock configurations distributing $N$ bosons over $M$ self-consistently optimized single-particle orbitals, with the Gross-Pitaevskii equation as the single-orbital special case. From the resulting ground state, the authors construct the reduced one-body density matrix, decompose the particle current into natural-orbital contributions, and extract the superfluid fraction from the derivative of the ground-state energy with respect to rotation, equivalently from the ratio of the moment of inertia to the rigid-body value. The normalized first-order correlation function $g^{(1)}(x,x')$, interpreted as the tunneling amplitude between two points weighted by the local densities, is the diagnostic that connects barrier height to the loss of superflow.","core_discovery":"The central claim is that the superfluid fraction $\\langle \\rho_s\\rangle_0$, defined as the zero-rotation limit of $\\langle\\rho_s\\rangle(\\Omega) = (2\\pi^2 N\\Omega)^{-1}\\partial E/\\partial\\Omega$, decreases monotonically as the rotating Gaussian barrier height $\\lambda$ is increased, while the condensation fraction, the largest eigenvalue of the reduced one-body density matrix, remains essentially constant across the same range of $\\lambda$. The authors show this for $N = 5, 8, 11$ bosons with contact interactions of strength $\\gamma$ from $0.1$ to $30$, using the multiconfigurational time-dependent Hartree method for bosons (MCTDHB). They also show that the ground-state energy is periodic in the dimensionless rotation frequency $\\Omega$, that strong barriers drag a finite fraction of particles even at infinitesimal rotation, and that the normalized one-body correlation function $|g^{(1)}(x,x')|^2$ develops an abrupt four-block suppression across the barrier that predicts the flow behavior.","pith_inferences":["If the independence of condensation and superfluid fraction survives at larger particle numbers, it would strengthen the SQUID analogy: the weak link controls dissipationless transport while the condensate, the analogue of the superconducting order parameter, remains intact.","A natural extension would be to vary the barrier width and shape at fixed height; the authors' argument suggests the controlling variable is the tunneling amplitude through the barrier, so different shapes with equal tunneling suppression should yield equal superfluid fractions.","The same energy-derivative definition could be applied to time-dependent protocols in which the barrier is suddenly raised, predicting metastable current decay rates for a ring prepared in a persistent-flow state."],"forward_implications":["If the central claim is right, a weak-link barrier can tune a few-boson ring continuously from a perfect superfluid to a near-rigid rotor without depleting the condensate.","Measuring the ground-state energy versus rotation frequency already determines the superfluid fraction, since the slope at $\\Omega=0$ gives $\\langle\\rho_s\\rangle_0$.","The correlation-function maps imply that particles are dragged by the barrier whenever the normalized tunneling amplitude across the barrier is suppressed, providing a static, measurable proxy for a dynamical property.","At very high barriers the density vanishes at the barrier peak while the angular-momentum distribution broadens only slightly, so loss of superflow is not accompanied by a large momentum change at rest."],"supporting_citations":[{"why":"Supplies the MCTDHB variational method used to compute the beyond-mean-field ground states.","marker":"[30]"},{"why":"Supplies the superfluid fraction definition via the energy derivative and moment-of-inertia ratio.","marker":"[43]"},{"why":"Earlier calculation of the superfluid fraction in a ring geometry that the present work extends beyond mean field.","marker":"[23]"},{"why":"Experimental realization of a ring Bose-Einstein condensate with a tunable weak link, motivating the barrier model.","marker":"[10]"},{"why":"Provides the reduced one-body density matrix and normalized correlation function used for the tunneling-amplitude diagnostic.","marker":"[27]"},{"why":"Describes the energy parabolas and swallow-tail structure for a rotating ring that underlies the periodicity argument.","marker":"[16]"},{"why":"Supports the use of the first-order correlation function in bosonic many-body systems for the numerical correlation maps.","marker":"[44]"}],"fun_headline_variants":["Weak link kills superflow, spares condensate","Barrier height drains superflow, leaves condensate","Superfluid fraction falls, condensate stays","Few-boson ring: barrier blocks flow, not order"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The numerical results assume that the truncated variational basis, the number of single-particle orbitals and Fock configurations in the MCTDHB calculation, is large enough to be converged for every interaction strength and barrier height reported, yet the paper never states the basis size and for the strongest interaction it reports only five particles.","fun_headline_variants_meta":{"raw":{"variants":["Weak link kills superflow, spares condensate","Barrier height drains superflow, leaves condensate","Superfluid fraction falls, condensate stays","Few-boson ring: barrier blocks flow, not order"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000556,"raw_usage":{"total_tokens":2604,"prompt_tokens":858,"completion_tokens":1746,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":474,"completion_tokens_details":{"reasoning_tokens":1683}},"tokens_in":474,"tokens_out":1746,"duration_ms":13404,"temperature":1.0,"reasoning_tokens":1683,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:46:33.448460+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Repeat the same parameter sets with an exact or much larger-basis few-boson calculation, increasing the number of orbitals until all observables stop changing, and check whether the superfluid fraction at fixed barrier height shifts by more than the paper's stated accuracy; alternatively, measure the condensation fraction while sweeping barrier height at fixed interaction strength in a cold-atom ring and look for a change larger than the few-percent variation reported in the table.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the MCTDHB variational method used to compute the beyond-mean-field ground states."},{"cited_title":"Muñoz Mateo, A","cited_arxiv_id":null,"evidence_quote":"Supplies the superfluid fraction definition via the energy derivative and moment-of-inertia ratio."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Experimental realization of a ring Bose-Einstein condensate with a tunable weak link, motivating the barrier model."},{"cited_title":"Navon, A","cited_arxiv_id":null,"evidence_quote":"Provides the reduced one-body density matrix and normalized correlation function used for the tunneling-amplitude diagnostic."},{"cited_title":"Eckel, J","cited_arxiv_id":null,"evidence_quote":"Describes the energy parabolas and swallow-tail structure for a rotating ring that underlies the periodicity argument."}],"review_version":1}