{"id":"ff03f5ad-5769-4903-b733-c0813c349977","arxiv_id":"2501.05301","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A rotating magnetic polarization protocol creates persistent azimuthal circulation in a three-well ring of dipolar bosons, with optimal frequencies predictable by a (N-1)U/J scaling relation.","lead":"This paper proposes a magnetostirring protocol that rotates the polarization direction of dipolar bosons in a three-well ring to create persistent circulating currents. The authors use exact numerical simulations to show the protocol works for intermediate rotation speeds, and they propose a scaling relation to predict optimal parameters for larger systems.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The core protocol is credible for N<=40, but the advertised scalability rests on an unvalidated (N-1)U/J frequency-scaling whose mean-field motivation fails in the circulation lobe.","rationale":"The reader's conditional verdict is the right one. The protocol itself has a clean mechanism: the final U=U_d configuration turns the interaction part into a constant, so the current operator is conserved after t_f, and the numerical evidence for N<=40 is consistent with that picture. I do not see an internal inconsistency in the central construction. The advertised scalability, however, is exactly the 'tentative' part of the paper. It is not enough to show that fixed small-N optimal frequencies produce growing circulation at constant (N-1)U/J; one must show that the location of the optimum in frequency space is independent of N. The mean-field argument that would supply that independence is discredited by Appendix B in the region where circulation is generated. Since a three-site Bose-Hubbard model has polynomial Hilbert-space dimension in N, this missing check is directly executable and should be a condition for the large-N claim. I therefore keep the verdict conditional rather than accepting the abstract's unqualified scalability statement.","tokens_in":14040,"tokens_out":15236,"duration_ms":150810,"concrete_test":"Perform exact time-dependent diagonalization at fixed g=(N-1)U/J=35 for N=60 and N=80 (Fock dimensions 1891 and 3321), scan (omega_xy, omega_z) on the same grid as Fig. 3, locate the circulation lobe, and compare the optimal frequencies and peak circulation with the small-N values and the Fig. 8 power-law extrapolation. If the large-N optima move outside the fit's confidence interval, or if using the small-N optimal frequencies at N=80 gives a circulation substantially below the re-optimized value, the scaling claim fails. As a complementary check, rerun Fig. 7 with the optimum chosen at N0=40 instead of N0=15; if the constant-g curve stops increasing, finite-N corrections to the scaling are not negligible.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The primary protocol claim is well supported for the simulated sizes: with U=U_d the final Hamiltonian is free up to a constant, so the generated average circulation is conserved for t>t_f, and exact diagonalization through N=40 shows circulation creation. The load-bearing weakness is the extrapolation to large N. Section 4.5 and Appendix A assume that the optimal frequencies (omega_xy, omega_z) depend on the exact many-body dynamics only through g=(N-1)U/J, so optima found on a small system can be reused for a larger system. The mean-field equations of B.1 motivate this scaling, but B.3 explicitly shows that the mean-field description has condensed fraction below 0.75 in most of the circulation lobe, i.e. it is invalid in much of the region where the protocol works. The exact evidence is Fig. 7, which fixes frequencies at the N0=15 optimum and shows that circulation grows with N at constant g; that does not test whether the optimal frequency pair itself is N-independent. Fig. 8's extrapolating fit is labelled 'tentative' and is never validated at large g. Moreover, for three sites the Fock dimension is only (N+1)(N+2)/2, so exact evolution well beyond N=40 is computationally feasible; the absence of such a check leaves the abstract's claim that the method 'overcomes computational limitations' unsupported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a magnetostirring protocol for a three-well ring of dipolar bosons, in which the dipole orientation follows a spherical spiral from the y-axis to the z-axis while the on-site interaction is tuned to the dipolar strength (U=Ud). At the end of the protocol, the interaction term cancels and the final Hamiltonian is free, so any generated azimuthal circulation is conserved. Using exact diagonalization for up to N=40 bosons and fourth-order Runge-Kutta time evolution, the authors show that the protocol creates a substantial average circulation for intermediate driving frequencies, that the final condensed fraction is often large, and that the protocol is robust against small deviations of the on-site interaction. They also propose a scaling rule based on keeping (N−1)U/J constant and a power-law fit to predict optimal frequencies for larger systems. The main numerical results are clearly presented and the protocol is physically plausible.","tokens_in":14339,"tokens_out":5288,"duration_ms":50252,"significance":"If the central protocol claim holds, the paper offers a new and experimentally accessible way to create persistent currents in atomtronic circuits, complementing existing stirring, phase-imprinting, and gauge-field methods. The paper is careful to use exact numerics for the simulated sizes and to test robustness, and the final-Hamiltonian argument that the generated circulation is conserved is clean and load-bearing. The condensed-fraction analysis and the mean-field comparison are also useful. However, the advertised scalability to large numbers of bosons rests on an empirical scaling assumption that is not validated by direct many-body calculations, and the mean-field approximation that motivates the scaling is shown to fail in most of the circulation-producing regime. The significance of the paper therefore depends on closing this gap; with that, it would be a solid contribution to the atomtronics literature.","major_comments":[{"comment":"The claim that optimal protocol frequencies can be transferred between different particle numbers at constant (N−1)U/J is only tested by fixing the N0=15 optimal frequencies and then varying N (Fig. 7). That test shows that the circulation fraction increases with N along this path, but it does not show that the optimal frequency pair itself is independent of N at constant g. The power-law fit in Fig. 8 is explicitly labeled 'tentative' and is never validated at large (N−1)U/J. Because the three-site Fock dimension is only (N+1)(N+2)/2, exact diagonalization for N=100 or more is computationally inexpensive; the paper should either perform such a check or substantially soften the abstract's and Sec. 5's claims about overcoming computational limitations.","section":"Sec. 4.5 and Appendix A"},{"comment":"The authors demonstrate in Appendix B.3 that the final condensed fraction is below the 0.75 threshold in most of the circulation lobe (Figs. 10 and 11), i.e., the mean-field model is invalid in much of the parameter region where the protocol is effective. Since the (N−1)U/J scaling in Sec. 4.5 is motivated by the structure of the mean-field equations (Appendix B.1), the mean-field failure removes the main theoretical support for the scaling. The scaling should be either derived from the exact many-body dynamics or supported by direct exact calculations at larger N.","section":"Appendix B"},{"comment":"The statement that the protocol 'overcomes computational limitations' and 'enables application to systems with large numbers of bosons' is not justified by the presented evidence. For three sites, the Hilbert-space dimension grows polynomially with N, and the authors already reach N=40; exact runs at N=100 or larger would be straightforward with standard exact diagonalization. The large-N extrapolation is therefore not necessary for the system at hand, and the claim as written overstates the gap being filled.","section":"Abstract and Sec. 5"}],"minor_comments":[{"comment":"Please specify the base of the logarithmic color scale used in panel (a), since the caption says 'logarithmic scale' but does not indicate base 10 or natural logarithm.","section":"Fig. 2"},{"comment":"It would be helpful to state the grid resolution used in the parameter scans of Figs. 3, 9, and 10; the text mentions that the grid is too coarse in the slow-frequency region without giving the actual step sizes in ω_xy and ω_z.","section":"Sec. 4.2 and Fig. 3"},{"comment":"The formula N2−1/N1−1 = U1/J1 / U2/J2 is correct but reads ambiguously because the left-hand side is a ratio of dimensionless numbers while the right-hand side involves ratios of interaction and hopping parameters; a brief derivation or a concrete example would improve clarity.","section":"Appendix A"},{"comment":"The definition of the one-body density matrix uses ρ_{j,k} = (1/N)⟨Ψ|a†_j a_k|Ψ⟩; please ensure the index ordering is consistent with the standard convention used in the natural-orbital analysis, since the ordering of j and k affects the interpretation of the off-diagonal elements.","section":"Eq. (5)"},{"comment":"The sentence 'to later rotate the direction of polarization, which induces a rotation in the condensate' could be smoothed to avoid the repetition of 'rotation'; this is a style point only.","section":"Introduction"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal, and the exact-diagonalization computations for N≤40 appear sound. The main concern is the unvalidated (N−1)U/J scaling that underpins the scalability claim; because exact simulations for larger N in this three-site geometry are computationally cheap, the authors can likely address this point directly. I would encourage the editor to ask for such a validation or for a clear rephrasing of the scalability claim before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The protocol works for the sizes tested; the scalability claim is the soft part.\n\nThe genuinely new piece is applying magnetostirring to a discrete three-well ring, with a spherical-spiral sweep of the polarization from in-plane to out-of-plane. The design is clean: set U=Ud so the final Hamiltonian is free, which makes the circulation conserved after the protocol. The exact-diagonalization results through N=40 look sound, the parameter-space lobe is physical (fast quench, adiabatic, and antidipolar limits all behave as expected), and the robustness scan against detuning U≠Ud is a useful experimental check. The authors also deserve credit for being upfront about where the mean-field model fails; Appendix B is honest and helpful.\n\nThe weak point is the large-N prediction. Section 4.5 and Appendix A assume the optimal frequencies depend on N only through (N−1)U/J. What Fig. 7 demonstrates is that if you fix the frequencies at the optimum for N0=15 and keep (N−1)U/J constant, the final circulation grows with N. That does not show the optimal frequency pair is unchanged with N. Fig. 8 fits a power law to the same optimal frequencies and calls it 'tentative'; extrapolating that fit to bigger (N−1)U/J is a guess, not a result. The mean-field equations that motivate the scaling are, by the authors' own analysis, invalid in most of the circulation lobe, so the motivation is shaky. And for three sites the Fock dimension is only (N+1)(N+2)/2; exact evolution up to N=100 or beyond is straightforward. Not doing that check leaves the abstract's 'overcomes computational limitations' unsupported. I would ask for a direct test at larger N and a toned-down abstract. I would also ask them to release code or at least data behind the main figures; it is a numerical paper and nothing is provided.\n\nThe citation pattern is fine: the prior magnetostirring and triple-well work is credited appropriately, with no citation inflation.\n\nThis paper will be of use to people building atomtronic circuits with dipolar species; as a protocol proposal it is plausible and worth taking seriously. But the selling point about scalability needs evidence. I would send it to peer review—the core numerical demonstration deserves refereeing—with the expectation of major or at least moderate revision on the scaling claims.","headline":"Credible protocol for N≤40, but the advertised large-N scaling is a tentative fit without a direct check.","tokens_in":14892,"tokens_out":2795,"would_cite":false,"duration_ms":27187,"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":"Rotating the dipole polarization drives dipolar bosons in a three-well ring into an excited state with a persistent azimuthal circulation.","keywords":["persistent currents","atomtronics","dipolar Bose gases","magnetostirring","Bose-Hubbard model","three-well ring","exact diagonalization","azimuthal circulation"],"falsifier":"Compute the optimal frequencies by exact diagonalization at an intermediate size, for instance $N=20$, with $(N-1)U/J$ fixed to the value used for the $N=8$ optimum, and check whether the $N=20$ circulation at the $N=8$-optimal frequencies matches the predicted growth; a mismatch larger than the fit's confidence interval would falsify the scaling rule. A complementary experimental falsifier is the absence of a persistent current signal in a three-well dipolar ring after the spherical-spiral rotation is applied at any intermediate frequency.","tokens_in":13825,"feed_emoji":"🌀","tokens_out":16085,"duration_ms":119847,"temperature":0.7,"pith_summary":"The paper proposes a protocol for creating persistent currents in a three-well atomtronic ring without stirring the trap or imprinting phase. The protocol starts with dipolar bosons polarized in the plane of the ring, then rotates the dipole orientation along a spherical spiral until the dipoles point perpendicular to the ring. The rotation transfers population into the third well and excites the system into a state with large average azimuthal circulation, and because the on-site and dipolar interactions are tuned to cancel at the end ($U=U_d$, dipoles along $z$), the final Hamiltonian is interaction-free and the circulation persists. Exact diagonalization of the Bose-Hubbard model for up to 40 bosons shows the circulation grows faster than the maximum single-atom value, meaning the protocol is more efficient for larger systems. The authors also show the optimal stirring frequencies for a small particle number can be reused for large systems if the combination $(N-1)U/J$ is kept fixed, making the scheme scalable beyond numerically tractable sizes.","feed_headline":"Rotating dipoles create a persistent current in a three-well ring","feed_subtitle":"A rotating magnetic field alone puts dipolar bosons into a persistent circulating state.","key_machinery":"The central object is the extended Bose-Hubbard Hamiltonian for a three-site ring with an anisotropic, long-range dipolar interaction, together with a time-dependent dipole orientation $\\boldsymbol{\\mu}(t)$ that traces a spherical spiral: $\\boldsymbol{\\mu}\\cdot\\mathbf{e}_x=\\sin(\\omega_{xy}t)\\cos(\\omega_z t)$, $\\boldsymbol{\\mu}\\cdot\\mathbf{e}_y=\\cos(\\omega_{xy}t)\\cos(\\omega_z t)$, $\\boldsymbol{\\mu}\\cdot\\mathbf{e}_z=\\sin(\\omega_z t)$. The anisotropic interaction breaks the reflection symmetry of the ring and, when the dipole turns in-plane, transfers bosons from the initially populated wells (sites 2 and 3) into the empty well (site 1). The protocol ends at $t_f=\\pi/(2\\omega_z)$ with $\\boldsymbol{\\mu}=\\mathbf{e}_z$; at that orientation, and with $U=U_d$, the inter-site interaction exactly cancels the on-site term in the shifted Hamiltonian, leaving $H'=-J\\sum_j(\\hat{a}^\\dagger_{j+1}\\hat{a}_j+\\mathrm{h.c.})$, a free-particle model. The azimuthal circulation is read from the operator $\\hat{L}_z$ that counts the imbalance between clockwise and counterclockwise tunneling, whose eigenstates are the Bloch waves of the triangle; the final state is a superposition of these circulation eigenstates. The two control frequencies $\\omega_{xy}$ and $\\omega_z$ are the only free protocol parameters, and the paper identifies the intermediate-frequency 'lobe' in which circulation is robustly created.","core_discovery":"The central claim is that a time-dependent rotation of the dipole polarization—the magnetostirring protocol—acts as an effective stirring force that drives an interacting dipolar Bose gas in a fully connected three-well ring from its ground state to an excited many-body state with a high average azimuthal circulation, and that this circulation is preserved after the driving stops. The preservation is guaranteed by a specific design: the on-site interaction strength is set equal to the dipolar coupling $U=U_d$, and the final dipole orientation is perpendicular to the ring plane, so the interaction term in the extended Bose-Hubbard Hamiltonian vanishes and the final dynamics is that of free bosons on a triangle. The paper supports this by exact time evolution (fourth-order Runge-Kutta on a Fock basis) for $N$ up to 40, showing that the final state overlaps with highly excited eigenstates, that the site occupations become imbalanced, and that the expectation value of the azimuthal circulation operator rises to a plateau. A further claim is that the optimal protocol frequencies ($\\omega_{xy}$ and $\\omega_z$) can be identified from small systems and transferred to larger ones using the scaling variable $(N-1)U/J$, so the protocol could be applied where direct simulation is impossible.","pith_inferences":["An experiment could use the lifetime of the final circulation as a direct probe of the $U=U_d$ cancellation: scanning a Feshbach resonance across the cancellation point should show the current surviving longest exactly at $U=U_d$ and decaying on either side; this diagnostic is implicit in the paper's robustness analysis.","The $(N-1)U/J$ scaling rule, if it holds, means the optimal frequencies are governed by the interaction energy per particle relative to tunneling; a tractable next computation is to find the lobe at an intermediate size (say $N=20$) and compare the optimum with the power-law extrapolation from $N=8$, which would test the rule where it is most needed.","Applying the same spherical-spiral stirring to a continuous toroidal condensate is a natural extension, but there the in-plane rotation rate must exceed the critical velocity for vortex nucleation; whether the discrete-ring result carries over is an open question the paper leaves implicit."],"forward_implications":["Persistent currents in atomtronic circuits can be created by a global, contactless manipulation—rotating the magnetic field that polarizes the dipoles—rather than by moving a stirrer or imprinting phase.","The final current survives because the interactions are engineered to vanish at the end of the protocol, so the persistent circulation is tied to a precise cancellation of $U$ and $U_d$.","For a fixed interaction strength and a fixed choice of frequencies within the effective lobe, the circulation per boson grows with particle number $N$, so the protocol is not only scalable but more efficient at larger sizes.","The optimal stirring frequencies for a large system can be predicted by simulating a small system with the same $(N-1)U/J$, reducing the computational cost from exponential in $N$ to a fixed small size.","The final state retains a high condensed fraction even though it is excited, because the interaction cancellation suppresses fragmentation at the end of the protocol."],"supporting_citations":[{"why":"First experimental observation of vortices in a dipolar condensate, establishing magnetostirring as a working technique that the protocol adapts.","marker":"[35]"},{"why":"Shows theoretically that rotating the polarization direction of a dipolar condensate creates vortex lattices, the mechanism the protocol exploits.","marker":"[36]"},{"why":"Identifies the triangular Bose-Hubbard trimer as the minimal model of a superfluid circuit and supplies the single-particle eigenbasis used to define the circulation measurement.","marker":"[42]"},{"why":"Demonstrates that choosing U = U_d with dipoles along z cancels the interaction term, leaving a free Hamiltonian; this is why the final circulation persists.","marker":"[46]"},{"why":"Provides the ground-state properties and anisotropy behavior of dipolar bosons in triple wells that set the initial condition and drive the population transfer.","marker":"[34]"},{"why":"Reviews persistent currents in ultracold gases and gives the azimuthal circulation operator underlying the protocol's target quantity.","marker":"[7]"}],"fun_headline_variants":["Magnetostirring turns on persistent currents in three-well atomtronic ring","Rotating dipoles induce lasting circulation in atomtronic ring","Dipolar bosons get stirred into persistent currents by rotating field","New magnetostirring protocol creates persistent flow in atomtronic loop","Stirring dipoles with magnets yields persistent current in ring circuit"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that optimal stirring frequencies depend on the system only through the combination $(N-1)U/J$—the effective interaction strength per particle relative to tunneling—so that frequencies found on a small, numerically tractable system remain optimal for a much larger one; the paper offers only a few numerical points and a tentative fit in support of this premise.","fun_headline_variants_meta":{"raw":{"variants":["Magnetostirring turns on persistent currents in three-well atomtronic ring","Rotating dipoles induce lasting circulation in atomtronic ring","Dipolar bosons get stirred into persistent currents by rotating field","New magnetostirring protocol creates persistent flow in atomtronic loop","Stirring dipoles with magnets yields persistent current in ring circuit"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00027,"raw_usage":{"total_tokens":1603,"prompt_tokens":903,"completion_tokens":700,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":519,"completion_tokens_details":{"reasoning_tokens":607}},"tokens_in":519,"tokens_out":700,"duration_ms":7337,"temperature":1.0,"reasoning_tokens":607,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:14:38.782342+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the optimal frequencies by exact diagonalization at an intermediate size, for instance $N=20$, with $(N-1)U/J$ fixed to the value used for the $N=8$ optimum, and check whether the $N=20$ circulation at the $N=8$-optimal frequencies matches the predicted growth; a mismatch larger than the fit's confidence interval would falsify the scaling rule. A complementary experimental falsifier is the absence of a persistent current signal in a three-well dipolar ring after the spherical-spiral rotation is applied at any intermediate frequency.","supporting_citations":[{"cited_title":"Klaus, T","cited_arxiv_id":null,"evidence_quote":"First experimental observation of vortices in a dipolar condensate, establishing magnetostirring as a working technique that the protocol adapts."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows theoretically that rotating the polarization direction of a dipolar condensate creates vortex lattices, the mechanism the protocol exploits."},{"cited_title":"Arwas, A","cited_arxiv_id":null,"evidence_quote":"Identifies the triangular Bose-Hubbard trimer as the minimal model of a superfluid circuit and supplies the single-particle eigenbasis used to define the circulation measurement."},{"cited_title":"Lahaye, T","cited_arxiv_id":null,"evidence_quote":"Demonstrates that choosing U = U_d with dipoles along z cancels the interaction term, leaving a free Hamiltonian; this is why the final circulation persists."},{"cited_title":"Rovirola, H","cited_arxiv_id":null,"evidence_quote":"Provides the ground-state properties and anisotropy behavior of dipolar bosons in triple wells that set the initial condition and drive the population transfer."}],"review_version":1}