{"id":"cb977c6f-8999-4767-a625-c6fc5f27fa40","arxiv_id":"2601.21144","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A trap-shaping and phase-imprinting protocol can prepare high-fidelity superpositions of clockwise and counter-clockwise persistent currents in a toroidal Bose-Einstein condensate, confirmed numerically.","lead":"The authors propose a method to create quantum superposition states of rotating currents in a ring-shaped Bose-Einstein condensate by shaping the trap with repulsive barriers and then imprinting phases with a light pulse. Numerical simulations show the method reaches high fidelity (>90%) and the states remain stable for seconds, which could be useful for compact rotation sensors.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Idealized instantaneous barrier removal and phase mask could inflate fidelity; no sensitivity analysis is provided.","rationale":"The reader's weakest assumption—instantaneous, perfectly aligned barrier removal and phase imprint—is indeed the most load-bearing physical-realizability concern. I agree with the reader's assessment and do not find a more severe internal inconsistency. The finite-operation effect is structurally distinct from the numerical idealization and directly determines whether the experimental claim holds. While estimates suggest the required timescales are favorable (trap period ~0.22 s, mode dynamics ~0.6 s), the paper provides no quantitative analysis, so a conditional verdict is appropriate. The analytical two-state model amplitude discrepancy (~30%) is not load-bearing because the central claim of high fidelity and stability is based on the GPE numerics, not the model. The paper has independent support: it reports a grid-convergence check (Sec. III) and a physical explanation of the radial oscillations (Sec. V), which strengthen confidence. Therefore, I recommend keeping the verdict CONDITIONAL rather than upgrading to ACCEPT, pending the proposed sensitivity test.","tokens_in":13840,"tokens_out":28979,"duration_ms":299097,"concrete_test":"Use the same GPE simulation and parameters as Sec. III, but replace the instantaneous barrier removal with a Gaussian ramp of duration τ_off, and replace the ideal phase mask with a smooth-edged mask: Φ(ϕ)=exp[iπ s_ε(ϕ)] where s_ε is an error-function step of width w. Add a lateral offset δ to the mask relative to the barrier minima. Compute the fidelity F(τ_off, w, δ) and the subsequent autocorrelation R(t) for m=3 and m=9, N=10^3 and 10^4. Scan τ_off up to 10 ms, w up to 3.3 μm, and δ up to 2 μm. If F≥0.90 for all experimentally plausible choices, the concern is resolved; if F falls below 0.90 for any such choice, the central claim needs to be qualified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim depends on the Sec. II assumption that the 2m barriers are removed 'suddenly' and simultaneously with an ideal π-phase mask on alternate sectors, with perfect spatial alignment. The numerical fidelity F>90% is computed for this ideal operation. Real experiments will have finite switching times (τ_off for the barrier ramp, τ_pulse for the phase imprint), a smooth phase-mask edge of width w (optical resolution), and possible lateral misalignment δ between the mask and the barrier positions. Finite τ_pulse introduces a mean-field phase evolution g2D|ψ|²τ/ħ that is spatially modulated by the ground-state density, adding an unintended phase pattern. Even for large V0, this term is not automatically negligible. A smooth or misaligned phase step places the π-phase jump where the condensate density is finite, creating a momentum kick ℏ/w and exciting modes; the reported F already includes the ideal sharp-step effect, but the paper does not analyze the additional degradation from finite w and δ. Since the claim is that the method 'can be realized experimentally with existing light sculpting techniques,' the lack of any timescale or robustness analysis leaves a gap between the numerical result and the physical protocol. The relevant timescales (trap period 0.22 s, mode spacing ~0.6 s) suggest fast switching may suffice, but this is not demonstrated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a protocol for engineering persistent-current superpositions in a toroidal BEC: prepare the ground state of a ring with 2m Gaussian barriers, then suddenly remove the barriers while imprinting a π phase on alternate sectors. The resulting |ENG⟩ state is compared via fidelity F=|⟨OAM|ENG⟩|² to the ideal counter-rotating superposition |OAM⟩=(|m⟩+|−m⟩)/√2. 2D GPE simulations for ⁸⁷Rb parameters give F>90% for m=3 and 9 with N=0, 10³, 10⁴ after optimization over the barrier height. The autocorrelation R(t) remains above 0.90 for several seconds. An analytical two-state model reproduces the dominant higher-mode dynamics and derives the selection rule m₂=3m₁.","tokens_in":14207,"tokens_out":14555,"duration_ms":155923,"significance":"If the protocol is experimentally robust, it offers a simple and high-efficiency alternative to Raman/Laguerre-Gauss methods for creating persistent-current superpositions, with potential applications in atomtronics and rotation sensing. The numerical parameters are clearly specified (r₀=50 μm, a_ho=r₀/10, FWHM=3.3 μm, grid 401 points, time step 5×10⁻⁶ s), and grid-convergence is checked. The analytical model in Appendix D is a genuine strength: it is derived from the GPE nonlinearity without fitted parameters and yields quantitative predictions (selection rule, oscillation period, amplitude) that agree with the numerics. The central numerical result is plausible and the paper is within the scope of the journal.","major_comments":[{"comment":"The central numerical claim F>90% is obtained under the idealization that the barriers are removed instantaneously and the π-phase mask is applied simultaneously with perfect spatial alignment (§II, Fig. 2). The abstract's statement that the method 'can be realized experimentally' therefore goes beyond what is demonstrated. Finite switching times, a finite phase-mask edge width, and lateral misalignment will introduce non-adiabatic excitations and phase errors. Please add a sensitivity study (e.g., ramp times τ_off/τ_pulse, edge width w, displacement δ) or explicitly restrict the realizability claim. Given the trap period of 0.22 s and the ~0.6 s mode-coupling period, a short analysis may suffice, but it is presently absent.","section":"§II and Abstract"},{"comment":"For m=3, N=10⁴, the manuscript states that the optimal barrier height 'is not reached in the investigated range'; the reported maximum therefore lies at the upper edge h_barrier/ω_trap=50. Since this case is one of the headline examples (F>90%), the optimization and the choice of barrier height used in §V are not fully established. Please either extend the scan while checking the coherence/fragmentation bound, or provide a quantitative criterion for the maximal admissible barrier height and show that the chosen value is the physically meaningful optimum.","section":"§IV, Fig. 3"},{"comment":"The two-state linearized model is derived under |c₂|² ≪ |c₁|², and its quantitative comparison with numerics is shown only for N=10³ (Figs. 7 and 9). In §V and Appendix C the model is also invoked to explain the N=10⁴ results, especially the weak oscillations for m=9, N=10⁴. The manuscript should state the range of g₂D (or N) over which the quantitative predictions for A and the period are expected to hold, and ideally show a decomposition or direct comparison for N=10⁴ as well.","section":"Appendix D and Figs. 7, 9"}],"minor_comments":[{"comment":"The linear-trap demonstration is qualitative; a fidelity value would make the claimed generality to arbitrary motional states more concrete.","section":"§II, Fig. 1"},{"comment":"Reference [36] cites only software documentation. Please provide a version/DOI or a more archival methods reference for the Trotter-Suzuki package.","section":"§III"},{"comment":"The notation h_barrier/ω_trap mixes an energy with a frequency; define h_barrier explicitly as an energy, or write h_barrier/(ℏω_trap).","section":"§IV"},{"comment":"Equation (9) and the following text contain '(2π/0.362)T' with a missing parenthesis; the dimensionless-to-physical conversion factor is clear but should be typeset cleanly.","section":"Appendix D"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope and the central numerical result is credible. The main gap is the disconnect between the idealized simulated operations and the claimed experimental realizability; a robustness analysis or a careful softening of the claim is needed before publication. No concerns about citation practice or novelty."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The core idea is a natural extension of phase-imprint and trap-shaping techniques, but the specific application to persistent-current superpositions is new, and the analytical two-state model with the m2=3m1 selection rule is a real contribution. The GPE simulations are clearly specified (87Rb parameters, ring radius 50 µm, barrier FWHM 3.3 µm close to DMD resolution), and the reported fidelities above 90% for m=3 and 9 across non-interacting and interacting cases are credible given the simulation setup. The paper is also honest about some of its own limitations: it explicitly notes that the optimal barrier height for m=3, N=10^4 is not reached in the scanned range, so that maximum sits at the edge of the scan. That is a minor caveat, not a fatal one.\n\nTwo soft spots deserve attention. First, the analytical amplitude estimate gives A≈0.088 against a numerical value of 0.068; the paper calls this good agreement, which is generous. The oscillation period and selection rule match well, but a ~30% discrepancy in the amplitude should be acknowledged more carefully. Second, the claim that node positions remain constant is supported only by 'analysis not shown here,' which is a weak evidence trail. The bigger physical concern is the idealized operation: barriers are removed suddenly and the phase mask is applied with perfect alignment and a sharp edge. No sensitivity analysis is provided for finite switching times, optical resolution, or lateral misalignment. The stability section shows fast timescales (trap period 0.22 s, mode oscillations ~0.6 s), which suggest the idealization may be harmless, but the paper does not demonstrate that. Since the abstract claims experimental feasibility with existing light sculpting, this gap is real but likely addressable with a short robustness check.\n\nOverall, the central numerical result is plausible, the analytical model is a genuine contribution, and the limitations are mostly acknowledged. This deserves a serious referee. I would want the referee to ask for a sensitivity analysis on the switching and phase-mask parameters, and a tighter treatment of the amplitude discrepancy.","headline":"Solid numerical proposal for persistent-current superpositions via barrier shaping plus phase imprint, with a genuine two-state selection-rule result; the main caveat is idealized switching and a fidelity maximum sitting at the edge of the scan.","tokens_in":14632,"tokens_out":2320,"would_cite":true,"duration_ms":25534,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A simple barrier-plus-phase-flip sequence can place a Bose-Einstein condensate into a superposition of counter-rotating persistent currents, with numerical fidelity above 90%.","keywords":["Bose-Einstein condensate","persistent current","superposition","phase imprint","toroidal trap","optical potential","Gross-Pitaevskii equation","atomtronics"],"falsifier":"Measure the fidelity as the barrier-removal time increases: if the fidelity drops below 90% when the removal time is a small fraction of the trap period (~0.22 s), the sudden-switch assumption fails. Alternatively, measure the population of the |9±⟩ mode for an interacting m=3 condensate with N=10^3; the two-state model predicts an oscillation amplitude of about 0.088 and a period of about 0.59 s, and a clear deviation would disprove the selection rule.","tokens_in":13773,"feed_emoji":"🌀","tokens_out":3728,"duration_ms":42082,"temperature":0.7,"pith_summary":"The paper proposes that the motional state of a Bose-Einstein condensate can be engineered by separately controlling density (with trapping barriers) and phase (with a sudden phase imprint). Applied to a toroidal trap, this yields a superposition of clockwise and counter-clockwise persistent currents, a state not yet created experimentally. Gross-Pitaevskii simulations show the engineered state matches the ideal cosine state with fidelity greater than 90% for m=3 and m=9, both without interactions and with up to 10^4 atoms, and stays stable with autocorrelation above 0.9 for several seconds. A two-state analytical model captures the dominant nonlinear coupling and predicts a selection rule in which the main higher mode has angular momentum 3m.","feed_headline":"Barrier trick builds superposed ring currents","feed_subtitle":"Numerical method reaches >90% fidelity and stays stable for seconds, enabling compact atom interferometers.","key_machinery":"The load-bearing mechanism is the combination of amplitude shaping by thin Gaussian barriers (which pre-impose the node structure of the target state) and a sudden π phase imprint that flips alternate sectors, converting the density modulation into a phase winding. The key analytic object is a two-state linearized model of the nonlinear Gross-Pitaevskii equation, which yields a selection rule: the dominant coupled higher mode has angular momentum 3m, with an oscillation period set by the chemical-potential difference (scaling as the square of angular momentum) and an amplitude given by U/(4πΔμ).","core_discovery":"The central claim is that a condensate wave function can be treated as amplitude and phase, each controllable independently: barriers shape the density to approximate the target's nodes, and a π phase flip on alternate sectors converts that density into the target phase structure. For a ring trap, starting from the ground state with 2m repulsive barriers, suddenly removing the barriers while imprinting the phase produces a state very close to cos(mφ), i.e., an equal superposition of |m⟩ and |−m⟩ persistent currents. The fidelity is optimized by the barrier height; interactions and higher m lower the achievable fidelity but keep it above 90% in the studied range. The engineered state is stabl","pith_inferences":["The protocol's sudden-switch assumption is the main practical risk: finite barrier-removal time or imperfect alignment will introduce non-adiabatic excitations and phase errors, so the actual fidelity in an experiment may fall below the numerical 90%.","The selection rule m2 = 3m1 arises from the quadratic nonlinearity of the GPE, suggesting that any weakly interacting ring condensate will exhibit this mode-coupling structure; the two-state model could be adapted to predict interaction-strength-dependent dephasing in other ring geometries.","If the method is extended to imbalanced superpositions, it could produce arbitrary persistent-current superpositions, opening a path to magnetic-field-gradient sensing or qubit encodings in the motional state.","The numerical claim of R(t) > 0.90 for several seconds assumes a pure mean-field condensate; finite-temperature effects, atom loss, or trap anharmonicity are not modeled, so robustness to these effects is a natural next test."],"forward_implications":["If correct, this method offers a high-efficiency route to persistent-current superpositions, avoiding the roughly 50% transfer-efficiency limit of two-photon Laguerre-Gauss methods.","The same amplitude-and-phase control applies to linear traps, so arbitrary excited states or superpositions of motional states could be engineered.","The stability of the node positions suggests the engineered states are usable for Sagnac rotation sensing, where the precession of the nodes measures rotation.","The two-state model predicts a specific oscillation period (~0.59 s) and population for the dominant higher mode, providing a clear experimental signature.","The protocol's generality means it could be extended to imbalanced superpositions, enabling richer interferometric and sensing schemes."],"fun_headline_variants":["Light barriers shape BEC into superposed ring currents","Optical potentials craft persistent-current superpositions","High-fidelity BEC superposition via dynamic light barriers","Barriers plus phase flip yield stable ring currents","Sculpted light superposes BEC persistent currents"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The protocol assumes the barriers can be removed and the phase imprinted suddenly and with perfect spatial alignment, so that the post-operation state is exactly the initial ground state multiplied by the phase mask.","fun_headline_variants_meta":{"raw":{"variants":["Light barriers shape BEC into superposed ring currents","Optical potentials craft persistent-current superpositions","High-fidelity BEC superposition via dynamic light barriers","Barriers plus phase flip yield stable ring currents","Sculpted light superposes BEC persistent currents"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.0005,"raw_usage":{"total_tokens":2226,"prompt_tokens":632,"completion_tokens":1594,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":376,"completion_tokens_details":{"reasoning_tokens":1520}},"tokens_in":376,"tokens_out":1594,"duration_ms":12702,"temperature":1.0,"reasoning_tokens":1520,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T07:04:25.478862+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the fidelity as the barrier-removal time increases: if the fidelity drops below 90% when the removal time is a small fraction of the trap period (~0.22 s), the sudden-switch assumption fails. Alternatively, measure the population of the |9±⟩ mode for an interacting m=3 condensate with N=10^3; the two-state model predicts an oscillation amplitude of about 0.088 and a period of about 0.59 s, and a clear deviation would disprove the selection rule.","supporting_citations":[],"review_version":1}