{"id":"8129bbde-458e-47bc-b811-aa4ca6be2a90","arxiv_id":"2411.17586","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"In an SOAM-coupled spin-1 condensate with 4 angular momentum transfer, the zero-to-annular-stripe transition is preceded by the softening of a symmetric double roton mode at lq = ±4.","lead":"This paper calculates the wobble modes of an ultracold atomic gas that arranges itself into a stripey ring, a supersolid phase. It shows a specific low-energy mode softens as the system turns into the stripe phase, giving experiments a clear signature to look for.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central roton-softening claim is not verifiable because the BdG spectra are computed in an unstated, untested harmonic-oscillator truncation; the critical q≈-0.016 and gap closure could shift with basis size.","rationale":"The central claim would be true only if the numerical BdG spectrum is converged and correctly implements the linearized GP equations. The paper provides neither the basis size nor any convergence test, so the reported zero of the roton gap at q≈-0.016 cannot be checked. This is load-bearing because the phase-transition order and the precursor claim are both inferred from that zero. A small basis would most affect the high-angular-momentum modes, exactly where the roton minimum sits (lq=±4). The additional sign inconsistency in Eq. (11) reinforces the need for a reproducible method statement, but the primary obstacle is the missing truncation data. I am not claiming the physics is wrong; the concern is that the evidence as published is insufficient to distinguish a true softening from a numerical artifact. If the requested convergence test confirms the gap closure, the CONDITIONAL verdict can be upgraded; if not, the central claim should be revised.","tokens_in":14649,"tokens_out":19371,"duration_ms":180171,"concrete_test":"Run the circularly-symmetric BdG solver of Ref. [42] for Ω0=2, q∈[-0.04,0] with nmax_x=nmax_y = 10, 20, 30, 40 (Nb=121,441,961,1681) and record the lowest lq=±4 eigenvalue. Verify that at q=-0.016 the gap changes by less than 1% between the two largest bases and that it crosses zero within 0.0005 of -0.016; otherwise the claimed gap-closing point is a truncation artifact. Separately, re-derive the third row of Eq. (11) from Eq. (7); if the minus signs are typos, re-diagonalize the corrected matrix and recompute the AS and VN spectra in Figs. 8 and 9 to confirm the transition assignments.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim, that a double symmetric roton mode at lq=±4 closes at the direct ZAM-AS transition, is a numerical statement: the gap is a BdG eigenvalue computed on a truncated harmonic-oscillator basis. The Appendix (Eqs. (8)-(11)) defines the method but never states the value of Nb or nmax_x=nmax_y used, and no convergence test appears anywhere in the paper. The BdG matrix is 6Nb x 6Nb and its elements are integrals of products of ground-state densities and Raman terms; too small a basis artificially raises high-momentum/radial energies and can move the roton minimum or produce a spurious zero. The reported q≈-0.016 and the coexistence with a first-order energy jump in Fig. 6(a) therefore rest on an unstated numerical parameter. A second, independent issue is that the published third row of the block matrix in Eq. (11) has overall minus signs relative to Eq. (7); if the code follows the printed matrix, the AS and VN spectra in Figs. 8 and 9 are not solving the stated problem. Neither issue is resolved by the text.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a quasi-2D spin-1 Bose-Einstein condensate with spin-orbital-angular-momentum coupling produced by Laguerre-Gaussian beams with orbital angular momentum transfer l=4. Using imaginary-time Gross-Pitaevskii evolution for 23Na parameters, it maps the ground-state phase diagram in the Raman-coupling/quadratic-Zeeman plane and identifies three phases: annular stripe (AS), vortex necklace (VN), and zero angular momentum (ZAM). It then solves the Bogoliubov-de Gennes equations in a harmonic-oscillator basis and reports that at low Raman coupling the ZAM-to-AS transition is first order and is preceded by the softening of a double symmetric roton mode at lq=±4, which closes the roton gap at the transition. The paper also identifies low-lying dipole, breathing, spin-dipole, and spin-breathing modes in the symmetry-broken phases.","tokens_in":14851,"tokens_out":7581,"duration_ms":114232,"significance":"If the central numerical result holds, the paper gives a concrete and falsifiable prediction: in a SOAM-coupled spin-1 BEC with l=4, the direct ZAM-AS transition at low Raman coupling is first order and is signalled by a symmetric double roton at lq=±4, closely analogous to the stripe-supersolid transition in linearly SO-coupled BECs. The choice l=4 is experimentally motivated by improved stripe contrast, and the double roton is not imposed but emerges from the BdG calculation on the GP ground states, with no fitted parameters. The phase diagram is cross-checked against energy derivatives and the Goldstone-mode count is consistent with symmetry breaking. The main weaknesses are numerical reproducibility: the basis truncation is never stated, no convergence tests are shown, and the Appendix contains a sign inconsistency in the projected BdG matrix that calls into question whether the printed equations match the solved problem.","major_comments":[{"comment":"The block matrix in Eq. (11) is inconsistent with the BdG equation (7) for the third row. In Eq. (7), the third row must be (P1)_{3j} u_j + (P2)_{3j} v_j = ω u_{-1}, with the matrix elements of P1 as defined in the Appendix. However, the printed Mkl_31, Mkl_32, and Mkl_33 all carry an overall minus sign, and Mkl_34 through Mkl_36 are also written with minus signs. If the diagonalization used the printed matrix, the spectra in Figs. 8 and 9, and the ZAM modes in Fig. 7, do not solve the stated BdG problem. The authors must either correct the signs or explicitly state that the printed matrix is a typographical error and that the calculations used the correct signs.","section":"Appendix, Eq. (11) and Eq. (7)"},{"comment":"The numerical spectra are not reproducible because the basis truncation is not specified. The Appendix defines Nb = (nmax_x+1)^2 but never states the value of nmax_x or Nb actually used, and no convergence checks are reported. The central quantitative results—the closing of the double roton gap at q ≈ -0.016 in Fig. 7(a) and the critical couplings Ω0 ≈ 4.6 and 7.4 in Fig. 9—are BdG eigenvalues computed in this truncated basis. A too-small basis can shift the roton minimum and the gap-closing point, so the authors should state the basis size and show that the lq = ±4 gap and the phase boundaries are converged with respect to it.","section":"Appendix, Eqs. (8)-(11)"},{"comment":"The claimed phase boundaries, especially the direct ZAM-AS transition at Ω0 = 2, rely on imaginary-time propagation reaching the true ground state rather than a metastable state selected by the initial guess. The text states that random guesses were used in addition to single-particle-inspired guesses, but no systematic comparison of energies from different initial conditions is presented. This should be documented, for example by reporting the lowest energies obtained from several initial conditions in the vicinity of the phase boundaries.","section":"Section III, imaginary-time propagation"}],"minor_comments":[{"comment":"The text states that the 23Na condensate has antiferromagnetic interactions with c2 < 0, but Eq. (5) and the chosen scattering lengths (a2 = 55.01 aB, a0 = 50 aB) give c2 = +1.33. In the usual spin-1 convention, antiferromagnetic interactions correspond to c2 > 0; the sign statement should be corrected.","section":"Section III and Eq. (5)"},{"comment":"The symbol hcc (and h*_cc) used in Mkl_12, Mkl_21, Mkl_23, and Mkl_32 is never defined. Presumably it denotes HΩ or H*_Ω from the definition of P1, but it should be defined explicitly to make the projection unambiguous.","section":"Appendix, definitions after Eq. (11)"},{"comment":"The caption says the transition occurs for q ⪆ -0.017, while the text gives q ≈ -0.016; these values should be reconciled. The caption also says 'there is a phase transition from the AS to the circularly symmetric lz = 0 (ZAM) phases,' which is confusing because the sweep in Fig. 7(a) is described in the text as a ZAM-to-AS transition as q decreases.","section":"Fig. 7(a) caption and Sec. IV"},{"comment":"The label 'AS-ZAM transition' in Fig. 6(a) is inconsistent with the text's description of a ZAM-to-AS transition; the terminology should be unified across the caption and the body.","section":"Fig. 6(a)"}],"recommendation":"major_revision","confidential_remarks":"I recommend major revision rather than rejection. The central physics—a roton-softening precursor at a symmetry-breaking transition—is plausible and the overall methodology is standard, but the paper currently lacks the numerical details needed to verify the specific quantitative predictions. The sign inconsistency in the Appendix must be resolved; if the code followed the printed matrix, the authors need to recheck the spectra. The paper is within the scope of the journal and, with the requested convergence data and corrected presentation, could become acceptable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a competent numerical extension of the standard GP+BdG approach to a new slice of parameter space for SOAM-coupled spin-1 BECs, and the central roton-softening story is plausible. That said, the paper is not yet reproducible: the basis truncation used for the BdG spectra is never stated, no convergence tests appear anywhere, and the printed block matrix in Eq. (11) has a sign inconsistency relative to Eq. (7) that the authors need to fix. A referee should see this, but the manuscript needs major revision first.\n\nWhat is new and worth credit: the paper is the first to compute collective excitations for a spin-1 SOAM condensate with l=4 and a quadratic Zeeman field, covering the annular stripe and vortex necklace phases. The phase diagram in the Raman-coupling vs. Zeeman-field plane is clearly mapped, with a direct ZAM-to-AS transition at low coupling and a ZAM-VN-AS sequence at higher coupling. The identification of the double roton mode at lq=±4 as the precursor to the supersolid stripe is consistent with earlier work on SO-coupled BECs and gives experimentalists a concrete signature to look for.\n\nSoft spots, in order of severity. First, the numerical basis size is unspecified. The appendix defines Nb=(nmax_x+1)^2 but never gives nmax_x, and there is no convergence check anywhere. Since the roton gap and the critical q≈-0.016 are eigenvalues of a truncated matrix, those quantitative results could shift with basis size. This is a load-bearing omission for the paper's central quantitative claim. Second, the third row of the block matrix in Eq. (11) carries extra minus signs compared to the P1 matrix in the appendix; if the code follows the printed matrix, the spectra in Figs. 8 and 9 are not solving the stated BdG problem. The authors should clarify which is correct. Third, minor: the text says \"antiferromagnetic interactions (c2 < 0)\" but then sets c2 = 1.33; for 23Na, with a2 > a0, c2 is positive, so the parenthesis is a likely typo.\n\nNone of these issues necessarily invalidates the qualitative picture — roton softening at a symmetry-breaking transition is exactly what one expects, and the paper engages honestly with the existing literature. But the missing numerical details and the sign inconsistency make the current version untrustworthy for quantitative use.\n\nThis paper is for cold-atom theorists and experimentalists working on SOAM coupling and supersolid detection. It deserves a serious referee, not a desk rejection, but the referee should insist on the basis size, convergence tests, and a corrected sign convention. I would not cite it in its current form.","headline":"Plausible roton-softening story for a new parameter regime, but missing numerical details and a sign error in the BdG matrix keep it from being citable as is.","tokens_in":15406,"tokens_out":5590,"would_cite":false,"duration_ms":46644,"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":"This paper claims that a double symmetric roton mode at $l_q=\\pm4$ closes exactly at the direct first-order transition from the zero-angular-momentum phase to the annular stripe supersolid in a spin-orbital-angular-momentum-coupled spin-1…","keywords":["supersolid","annular stripe phase","spin-orbital-angular-momentum coupling","spin-1 Bose-Einstein condensate","Bogoliubov-de Gennes","roton mode","quadratic Zeeman effect","vortex necklace"],"falsifier":"Take the same $c_0=42.57$, $c_2=1.33$, $l=4$ parameters and recompute the BdG spectrum at $\\Omega_0=2$ with increasing basis size (for instance $N_b$ from 36 to 100 and beyond); if the double roton gap at $l_q=\\pm4$ does not descend to zero at $q\\approx-0.016$, or if the low-lying frequencies shift beyond numerical tolerance as $N_b$ grows, the claimed direct first-order ZAM-AS transition is not established. An experiment could also look for the roton dip at $l_q=\\pm4$ in the density response of a sodium-23 condensate with $4\\hbar$ Raman transfer.","tokens_in":14394,"feed_emoji":"🔄","tokens_out":10371,"duration_ms":91315,"temperature":0.7,"pith_summary":"This paper studies the collective excitations of a quasi-two-dimensional spin-1 Bose-Einstein condensate with spin-orbital-angular-momentum coupling, using a $4\\hbar$ orbital angular momentum transfer so that an annular stripe supersolid phase can appear. It establishes that for realistic antiferromagnetic interactions (sodium-23 parameters) the ground-state phase diagram contains three phases, namely the annular stripe, vortex necklace, and zero angular momentum phases, and that the transitions between them are accompanied by the softening of a double symmetric roton mode. The central result is that at low Raman coupling there is a direct, first-order transition from the zero angular momentum phase to the annular stripe phase, signaled by the roton gap at $l_q=\\pm4$ closing at $q\\approx-0.016$. A sympathetic reader would care because a closing roton gap is the standard precursor to supersolidity, and the paper identifies the parameter regime and the observable mode that would certify the annular stripe phase in experiments using larger orbital angular momentum transfer.","feed_headline":"Closing roton gap reveals direct route to annular supersolid","feed_subtitle":"The softening double mode is the fingerprint of the gas turning into a striped supersolid.","key_machinery":"The central object is the double symmetric roton mode: in the circularly symmetric ZAM phase, excitations carry a magnetic quantum number $l_q$, and the pair $l_q=\\pm4$ are degenerate because the Bogoliubov-de Gennes equations are invariant under $l_q\\to-l_q$ together with interchange of the $m=+1$ and $m=-1$ spin components. The paper tracks this pair's frequency as a function of $q$ and $\\Omega_0$; its softening to zero marks the onset of the AS phase. The computation linearizes the Gross-Pitaevskii equation around the mean-field ground state and solves the resulting BdG problem by expanding quasiparticle amplitudes in a truncated basis of two-dimensional harmonic oscillator eigenstates, then diagonalizing the resulting $6N_b\\times 6N_b$ matrix with a sparse eigensolver. The Raman coupling profile $\\Omega(r)=\\Omega_0\\,e^{(l/2)(r/w)^l}e^{-2r^2/w^2}$ with $l=4$ selects the $l_z=\\pm4$ single-particle minima that the AS phase occupies.","core_discovery":"On its own terms, the paper's central discovery is that the zero-angular-momentum (ZAM) phase of a SOAM-coupled spin-1 condensate with $l=4$ becomes unstable to the annular stripe (AS) phase through a double roton mode rather than through a competing intermediate phase. For $\\Omega_0=2$, as the quadratic Zeeman field $q$ is lowered, the two degenerate excitation branches with $l_q=\\pm4$ soften symmetrically and reach zero frequency at $q\\approx-0.016$, exactly where the ground state switches from ZAM to AS; the discontinuity in $\\partial E_0/\\partial q$ and in the dipole and breathing mode frequencies marks the transition as first order. At higher Raman coupling the same $l_q=\\pm4$ roton gap closes at the continuous ZAM-to-vortex-necklace boundary, showing that the double roton softening accompanies rotational symmetry breaking generally. The AS phase, which condenses in a superposition of $l_z=+4$ and $l_z=-4$ single-particle states, breaks both $U(1)$ and rotational symmetry and consequently supports two Goldstone modes; the low-lying dipole, breathing, spin-dipole, and spin-breathing modes are identified in each phase.","pith_inferences":["The paper implies, without stating it, that the roton minima sit at $l_q=\\pm l$ for any integer transfer $l$, so for $l=4$ the azimuthal stripe period should be $\\pi/2$; measuring that period would test the mechanism directly.","Because the same double roton softening appears at both the first-order ZAM-AS and the continuous ZAM-VN transitions, a natural extension is to look for an additional soft mode (for example a quadrupole branch) at the tricritical point $(\\Omega_0\\approx 3.5,\\;q\\approx -0.06)$ where the three phases meet.","A beyond-mean-field or finite-temperature calculation could decide whether the mean-field first-order character of the ZAM-AS transition survives fluctuations, since roton softening often signals an instability that thermal or quantum fluctuations can preempt.","The absence of a single dominant dipole/breathing frequency in the vortex necklace phase suggests that experimental identification of its collective modes would require angular-momentum-resolved probes rather than the usual trap-modulation spectroscopy."],"forward_implications":["At low Raman coupling ($\\Omega_0\\approx2$), the direct ZAM-to-AS transition is first order, so the roton gap at $l_q=\\pm4$ closes exactly at the phase boundary and the dipole and breathing mode frequencies jump discontinuously.","The same double roton mode at $l_q=\\pm4$ softens at both the direct ZAM-AS transition and the continuous ZAM-VN transition, making it a consistent precursor of the loss of rotational symmetry.","Because the AS phase breaks both $U(1)$ gauge and rotational symmetry, its spectrum contains two zero-energy Goldstone modes, in contrast to the single Goldstone mode of the ZAM phase.","For $\\Omega_0>3.5$ the ZAM phase is separated from the AS phase by the intermediate vortex necklace phase, so the experimentally cleanest route to the annular stripe supersolid is at low Raman coupling where the transition is direct.","The choice $l=4$, with the associated $\\Omega(r)$ profile, is what makes the annular stripe phase energetically accessible; the same model with $l=1$ does not produce the AS phase in the computed phase diagram."],"supporting_citations":[{"why":"Provides the single-particle SOAM Hamiltonian and the degenerate $l_z=\\pm4$ minima that seed the annular stripe phase.","marker":"[30]"},{"why":"Shows that larger orbital angular momentum transfer improves stripe contrast and supplies the Gaussian-profile $\\Omega(r)$ used in the model.","marker":"[39]"},{"why":"Earlier spin-1 SOAM BdG study that gives the harmonic-oscillator basis expansion and the $l_q\\to -l_q$ symmetry used here.","marker":"[42]"},{"why":"Analogous direct transition with double roton softening in a linear-Raman spin-1 BEC, the qualitative template for the ZAM-to-AS transition.","marker":"[25]"},{"why":"Establishes roton-mode softening as a precursor to stripe supersolidity in spin-orbit-coupled condensates.","marker":"[20]"},{"why":"Theoretical framework for roton modes and supersolid stripes used to interpret the $l_q=\\pm4$ branch.","marker":"[21]"},{"why":"Experimental realization of spin-1 SOAM coupling with $l=1$, setting the coupling geometry and spinor form the model adopts.","marker":"[34]"},{"why":"Companion spin-1 SOAM experiment supplying the $S_x,S_y,S_z$ formalism and the experimentally accessible regime.","marker":"[35]"},{"why":"First SOAM experiment with Laguerre-Gaussian beams, the physical basis of the laser geometry considered.","marker":"[33]"}],"fun_headline_variants":["Double roton softening reveals direct supersolid path","Roton gap collapse flags annular stripe phase","Two roton modes herald supersolid stripe transition","Annular supersolid emerges via double roton mode","Softening roton pair marks supersolid stripe onset"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The numerical excitation spectra rest on the truncated harmonic-oscillator basis being large enough to converge the Bogoliubov-de Gennes modes, but the paper does not state the basis size $N_b$ or $n_{\\max,x}$ used and gives no convergence tests; if the basis is too small, the roton branches and the reported critical values $q\\approx-0.016$, $\\Omega_0\\approx 4.6$, and $\\Omega_0\\approx 7.4$ could shift or be numerical artifacts.","fun_headline_variants_meta":{"raw":{"variants":["Double roton softening reveals direct supersolid path","Roton gap collapse flags annular stripe phase","Two roton modes herald supersolid stripe transition","Annular supersolid emerges via double roton mode","Softening roton pair marks supersolid stripe onset"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000202,"raw_usage":{"total_tokens":1432,"prompt_tokens":1044,"completion_tokens":388,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":660,"completion_tokens_details":{"reasoning_tokens":315}},"tokens_in":660,"tokens_out":388,"duration_ms":4600,"temperature":1.0,"reasoning_tokens":315,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:56:42.378915+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the same $c_0=42.57$, $c_2=1.33$, $l=4$ parameters and recompute the BdG spectrum at $\\Omega_0=2$ with increasing basis size (for instance $N_b$ from 36 to 100 and beyond); if the double roton gap at $l_q=\\pm4$ does not descend to zero at $q\\approx-0.016$, or if the low-lying frequencies shift beyond numerical tolerance as $N_b$ grows, the claimed direct first-order ZAM-AS transition is not established. An experiment could also look for the roton dip at $l_q=\\pm4$ in the density response of a sodium-23 condensate with $4\\hbar$ Raman transfer.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the single-particle SOAM Hamiltonian and the degenerate $l_z=\\pm4$ minima that seed the annular stripe phase."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows that larger orbital angular momentum transfer improves stripe contrast and supplies the Gaussian-profile $\\Omega(r)$ used in the model."},{"cited_title":"Banger, Rajat, A","cited_arxiv_id":null,"evidence_quote":"Earlier spin-1 SOAM BdG study that gives the harmonic-oscillator basis expansion and the $l_q\\to -l_q$ symmetry used here."},{"cited_title":"Yu, Phys","cited_arxiv_id":null,"evidence_quote":"Establishes roton-mode softening as a precursor to stripe supersolidity in spin-orbit-coupled condensates."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Theoretical framework for roton modes and supersolid stripes used to interpret the $l_q=\\pm4$ branch."},{"cited_title":"Chen, K.-Y","cited_arxiv_id":null,"evidence_quote":"Experimental realization of spin-1 SOAM coupling with $l=1$, setting the coupling geometry and spinor form the model adopts."},{"cited_title":"Chen, L.-R","cited_arxiv_id":null,"evidence_quote":"Companion spin-1 SOAM experiment supplying the $S_x,S_y,S_z$ formalism and the experimentally accessible regime."},{"cited_title":"Zhang, T","cited_arxiv_id":null,"evidence_quote":"First SOAM experiment with Laguerre-Gaussian beams, the physical basis of the laser geometry considered."}],"review_version":1}