{"id":"c54c02e0-2480-4331-b4a4-a7b395b652f8","arxiv_id":"2412.19285","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Finite-temperature Hartree-Fock-Bogoliubov calculations show that thermal fluctuations melt the supersolid stripe phase of a homogeneous spin-orbit-coupled spin-1 BEC, while quantum fluctuations enlarge it.","lead":"This paper computes how temperature reshapes the phase diagram of a spin-orbit-coupled spin-1 Bose-Einstein condensate, showing that the supersolid stripe phase melts as temperature increases. It matters because it provides finite-temperature phase boundaries that ultracold-atom experiments could test.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The ST-PW phase boundary is mapped from the PW roton gap alone; because the stripe phase is never solved at finite T and no free-energy comparison is made, the reported boundary is likely a spinodal rather than the equilibrium coexistence curve.","rationale":"After reading the manuscript in good faith, I find the reader's weakest_assumption is the correct central risk. The paper's contribution is a finite-temperature phase diagram obtained by HFB-Popov theory, and every boundary in Figs. 3, 5(b), 6(c) is extracted from a single-phase excitation gap (roton gap for ST-PW/ST-ZM, condensate momentum for PW-ZM). The ST phase is never an input; it is inferred as the region where the PW or ZM phase becomes dynamically unstable. This matters because the equilibrium phase boundary at finite T is defined by equality of the grand potential between the two phases, and a roton-gap closure marks the point where the PW phase loses local stability. For a second-order quantum transition these coincide only at T=0; at finite T the transition can be first-order, and then the gap-closure point is the spinodal, lying inside the coexistence region. The paper provides no evidence that the transition remains second-order, and no calculation of the stripe-phase free energy. The 'quantum fluctuations amplify supersolidity' claim (Sec. 3.1, Fig. 4c) is similarly inferred from the shift of the PW roton gap, so it inherits the same assumption. I do not argue the qualitative direction (thermal fluctuations usually do melt supersolids) is impossible, but the quantitative phase diagrams and the asymmetry between quantum and thermal effects are not secured by the calculation as presented. Other concerns (unreported box size L, no error bars on extrapolated fits, possible grid resolution effects on the roton minimum) are real but secondary; they affect accuracy, whereas the roton-gap-vs-free-energy issue affects the definition of the object that is plotted. A direct free-energy comparison for one representative temperature slice would settle the point and is computationally feasible within the same HFB framework, though it requires a Bloch-periodic BdG treatment of the stripe phase. In the meantime, the conditional verdict is appropriate: the qualitative melting direction is plausible, but the reported boundaries should be treated as indicative rather than definitive.","tokens_in":12887,"tokens_out":6280,"duration_ms":58473,"concrete_test":"At c0n=1, c2n=0.1, ϵ=-1, T=0.4Tc, add the stripe condensate ansatz ψ(x)=√nc(α_+ e^{ikx}, α_0, α_- e^{-ikx}) to the HFB-Popov scheme (Bloch-periodic BdG equations), and compute the grand potential Ω[ψ] for both the stripe and PW solutions as Ω_Raman is varied across the suspected boundary. Determine where the two free energies cross. Compare that crossing with the roton-gap-closure value Ω_c1 from Fig. 1(b). If the crossing is shifted by more than ~5% in Ω_Raman (or if one solution never becomes globally stable), the reported ST-PW boundary is a spinodal, not the coexistence line, and the phase diagrams and melting claim should be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central finite-temperature phase diagrams (Figs. 3, 5(b), 6(c)) locate the ST-PW and ST-ZM boundaries by solving the HFB-Popov equations only for the homogeneous PW/ZM phases and detecting where the roton gap in that phase closes (Sec. 3.1, see also Secs. 3.2-3.3). The stripe phase is never constructed at finite T; its existence is inferred retroactively from the instability of the PW/ZM phase. The gap closure in a one-phase calculation locates the spinodal (limit of metastability) of that phase, which coincides with the equilibrium phase boundary only if the transition is continuous and the second phase's free energy crosses exactly at that point. The paper offers no evidence that the finite-T ST-PW transition remains second-order, and thermal fluctuations generically favor first-order transitions in supersolid systems. Consequently the reported 'melting' of the ST phase and the claimed opposite effects of quantum vs thermal fluctuations on the boundary (Sec. 4) are quantitative statements about a conjectured equilibrium curve, not one that has been computed from a free-energy comparison. The quadratic-polynomial extrapolation of the roton gap to zero further compounds this: the boundary is a fitting artifact unless verified by direct calculation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies a homogeneous three-dimensional Raman-induced spin-orbit-coupled spin-1 Bose-Einstein condensate with repulsive interactions, using Hartree-Fock-Bogoliubov theory with the Popov approximation. The authors solve the generalized Gross-Pitaevskii equations and Bogoliubov-de Gennes equations self-consistently at finite temperature, compute roton gaps in the plane-wave and zero-momentum phases, and extract phase boundaries between the stripe, plane-wave, and zero-momentum phases in the T-Ω and T-ε planes. They report that the supersolid stripe phase melts with increasing temperature, and that quantum fluctuations enlarge the stripe region while thermal fluctuations shrink it.","tokens_in":13195,"tokens_out":4617,"duration_ms":44169,"significance":"The paper addresses a genuine gap in the literature: finite-temperature phase diagrams of spin-orbit-coupled spin-1 condensates, including the quadratic Zeeman field as a tunable parameter. The explicit HFB-Popov equations, the self-consistent iteration scheme, and the observation of thermal roton-gap opening are useful and reproducible in structure. If the extracted boundaries were equilibrium phase boundaries, the contrasting roles of quantum and thermal fluctuations would be an interesting and nontrivial result. However, as detailed below, the central quantitative claim is currently supported only by a spinodal construction, so the significance is conditional on additional calculation or careful reframing.","major_comments":[{"comment":"The ST-PW and ST-ZM boundaries are obtained exclusively from the roton gap of the homogeneous PW or ZM phase; the stripe phase is never solved at finite temperature and no free-energy comparison is made. A roton-gap closing in a one-phase calculation locates the spinodal (limit of metastability) of that phase, which coincides with the equilibrium coexistence boundary only if the transition is continuous and the second phase's free energy crosses exactly at that point. Since thermal fluctuations can make the transition first order, the reported 'melting' curves are not demonstrated to be equilibrium phase boundaries. The central claims in Sec. 4 should be either supported by a finite-temperature stripe-phase calculation with free-energy comparison, or explicitly reframed as spinodal lines.","section":"Secs. 3.1-3.3, Figs. 3, 5(b), 6(c)"},{"comment":"The critical points are extracted by quadratic polynomial fits and linear extrapolations of numerical data with no error bars, no fit-window sensitivity analysis, and no convergence checks with respect to the BdG grid size or box length. In particular, the PW-ZM boundary is obtained by linear extrapolation after the self-consistent iteration fails to converge near Ωc2. The reported boundary shifts (e.g., Ωc1 decreasing and Ωc2 increasing with temperature) should be accompanied by estimates of the extrapolation uncertainty; otherwise the reader cannot judge whether the shifts are physical or artifacts of the fitting procedure.","section":"Sec. 3.1, Figs. 1(b), 2(a); Sec. 3.2 Fig. 5(a); Sec. 3.3 Fig. 6(b)"},{"comment":"The claim that quantum fluctuations 'amplify supersolidity' is based on comparing the roton gap at T=0 with and without quantum fluctuations at a single point, Ω=2.0. This is suggestive but does not establish that the boundary shift is robust across the phase diagram, especially because the self-consistent calculation breaks down near the boundary. The authors should either provide the Ω-dependence of the gap over the full PW region or qualify the conclusion accordingly.","section":"Sec. 3.1, 'Fluctuations and supersolidity' and Fig. 4"}],"minor_comments":[{"comment":"The caption states 'ϵ=1' while the text and all other figures use ϵ=-1; please correct the mismatch.","section":"Fig. 3"},{"comment":"There is a typo in 'the the PW-ZM phase boundary'; please fix.","section":"Sec. 3.1"},{"comment":"The notation c0n and c2n is ambiguous: it should be clarified that these are dimensionless products c0 n and c2 n, not products of two symbols where n appears both as a subscript and as density.","section":"Sec. 2, Eqs. (6)"},{"comment":"The box size L used for the BdG momentum grid is not stated; since the grid spacing is Δq=2π/L, L must be specified for reproducibility.","section":"Sec. 3.1"},{"comment":"The caption of the inset says 'Ω=1.96 at T=0 and T=0.4Tc'; please clarify which curves correspond to which temperature and whether Ω=1.96 lies in the PW phase at T=0.","section":"Fig. 1"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal, and the HFB-Popov machinery is standard. The main concern is that the phase boundaries are spinodal curves, not demonstrated equilibrium coexistence lines. The authors should be asked to either add a finite-temperature stripe-phase calculation with free-energy comparison or explicitly revise the claims to refer to spinodal lines. I see no issue with novelty of the method, but the quantitative claims in Sec. 4 need to be brought in line with what the calculation actually shows."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take. The paper does something genuinely new: it takes the HFB-Popov machinery that was already applied to pseudospin-1/2 Raman BECs and to spin-1 BECs without SO coupling, and runs it for a Raman-induced spin-1 SO-coupled BEC. The T-Omega and T-epsilon phase diagrams, including the double roton structure at finite T, are new. The qualitative conclusion—the stripe region shrinks as T rises, and quantum fluctuations push the ST-PW boundary the other way—is consistent with what I'd expect and is supported by the roton-gap behavior. That's a real contribution.\n\nThe main soft spot is structural. The ST-PW and ST-ZM boundaries are located by solving the homogeneous PW or ZM phase at finite T and finding where its roton gap closes. The stripe phase itself is never constructed at finite T, and there is no free-energy comparison between the phases. So the extracted boundary is the spinodal of the homogeneous phase, not the equilibrium coexistence curve, unless the transition remains continuous and the second phase's free energy crosses exactly at that point. The stress-test note makes this point and I think it lands. The paper cites refs [28,39] for the same approach, which is fair, but those papers have the same limitation. For a quantitative phase diagram, this matters: thermal fluctuations can turn a continuous transition first-order in supersolid-like systems. I'd want the authors to either check the nature of the transition at at least a few representative temperatures or label the curves as instability lines rather than equilibrium boundaries.\n\nThe other issues are smaller but real. The box size L is never given; with N=500 grid points and periodic boundary conditions, the momentum step is 2pi/L, and the results (especially near the roton minimum) depend on that. That's an omission that blocks reproduction. The critical points come from quadratic polynomial fits and linear extrapolations with no error bars and no convergence tests on the grid size or the under-relaxation parameter. On the positive side, the equations are written out, the self-consistent iteration is described, and the paper is honest about the fitting procedure.\n\nWho is this for? Someone working on finite-temperature properties of SO-coupled spinor BECs, and experimental groups studying stripe phases in spin-1 gases. I'd cite it for the qualitative T-Omega and T-epsilon maps, with a caveat about the boundary definition. It deserves a serious referee; the missing L and the spinodal-vs-coexistence question are addressable in revision, and the paper's core claim is plausible.","headline":"First finite-T HFB-Popov phase diagram for spin-1 Raman SO-coupled BECs; the qualitative melting story looks right, but the ST-PW boundary is computed from one-phase roton gaps and is likely a spinodal, not an equilibrium coexistence curve.","tokens_in":13716,"tokens_out":3418,"would_cite":true,"duration_ms":28911,"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":"Thermal fluctuations melt the supersolid stripe phase; quantum fluctuations amplify it.","keywords":["spin-orbit coupling","spin-1 Bose-Einstein condensate","supersolid stripe phase","Hartree-Fock-Bogoliubov theory","Popov approximation","roton gap","finite-temperature phase diagram","quantum and thermal fluctuations"],"falsifier":"A direct finite-temperature free-energy comparison of the stripe and plane-wave states would settle the point: if the coexistence value of $\\Omega$ or $\\epsilon$ differs from the roton-gap closing value, or if the transition shows hysteresis, the roton-closing criterion used to draw the phase boundaries is incomplete.","tokens_in":12695,"feed_emoji":"🌡️","tokens_out":10014,"duration_ms":82286,"temperature":0.7,"pith_summary":"This paper asks what happens to the supersolid stripe phase of a spin-orbit-coupled spin-1 Bose-Einstein condensate when the gas is heated. Using Hartree-Fock-Bogoliubov theory with the Popov approximation, the authors compute the finite-temperature phase diagram of a homogeneous three-dimensional gas and locate the phase boundaries by tracking the roton gap and the condensate momentum. They find that the stripe phase melts: as temperature rises, the stripe-to-plane-wave boundary moves to lower Raman coupling and the plane-wave-to-zero-momentum boundary moves to higher coupling, so the plane-wave phase widens at the expense of the stripe and zero-momentum phases. They also find that zero-point (quantum) fluctuations narrow the roton gap and thereby favour the stripe, while thermal fluctuations open the gap and destabilize it. The result matters because it gives a concrete prediction for how the supersolid phase of an experimentally accessible spinor gas disappears with temperature, with the quadratic Zeeman field as an extra tuning knob.","feed_headline":"Heat melts the supersolid stripe phase in a spin-1 Bose gas","feed_subtitle":"Finite-temperature maps show the stripe phase shrinks as thermal fluctuations open the roton gap.","key_machinery":"The central object is the roton gap, the minimum of the lowest Bogoliubov excitation branch at a finite wavevector, evaluated self-consistently in the plane-wave or zero-momentum phase. A vanishing roton gap marks the onset of a density-wave (stripe) instability, so the gap as a function of $\\Omega$ or $\\epsilon$ locates the stripe-phase boundary, while the condensate momentum $k$ locates the plane-wave-to-zero-momentum boundary. The numerical machinery is the Hartree-Fock-Bogoliubov framework with the Popov approximation, which keeps normal thermal and quantum fluctuations but drops anomalous densities, and the coupled generalized Gross-Pitaevskii and Bogoliubov-de-Gennes equations are solved self-consistently on a momentum grid.","core_discovery":"For a homogeneous antiferromagnetic spin-1 condensate with $c_0 n = 1$ and $c_2 n = 0.1$ at a fixed quadratic Zeeman field $\\epsilon = -1$, the paper reports that the critical Raman coupling $\\Omega_{c1}$ for the stripe-to-plane-wave transition falls from about $1.96$ at $T = 0$ as the temperature rises, while the plane-wave-to-zero-momentum critical coupling $\\Omega_{c2}$ rises from about $2.78$. The supersolid stripe phase therefore occupies a smaller region of the $T{-}\\Omega$ plane at finite temperature: it melts into the plane-wave phase, and the zero-momentum phase also loses territory to the plane-wave phase. The same melting appears in the $T{-}\\epsilon$ plane, both for the stripe-to-plane-wave boundary at $\\Omega = 1.8$ and for the direct stripe-to-zero-momentum boundary at $\\Omega = 0.2$. The paper attributes the melting to thermal fluctuations, which open the roton gap that closes at the instability; quantum fluctuations act in the opposite direction, closing the gap further and enlarging the stripe region relative to the $T=0$ mean-field prediction.","pith_inferences":["A direct finite-temperature free-energy comparison between the stripe and plane-wave states would test whether the roton-closing criterion coincides with the true coexistence line; the extrapolated boundaries would shift if the transition becomes first order.","The opposing quantum and thermal shifts imply that the stripe-to-plane-wave boundary may be non-monotonic in $T$ at very low temperature, with a slight initial strengthening of the stripe before thermal melting dominates; the paper does not resolve this regime.","Since the Popov approximation neglects anomalous densities, including them could move the quantitative boundary locations, especially near the roton minimum where pairing fluctuations are largest; this is a testable extension rather than a claim of the paper."],"forward_implications":["At any fixed temperature below $T_c$, the stripe supersolid survives only for Raman couplings below a threshold that decreases with temperature, so a sample that is a stripe at $T = 0$ can become a plane-wave superfluid when heated.","Heating widens the plane-wave phase from both sides: it consumes part of the stripe phase on the low-coupling side and part of the zero-momentum phase on the high-coupling side of the $T{-}\\Omega$ diagram.","Varying the quadratic Zeeman field drives the same physics: at fixed $\\Omega$, the stripe-to-plane-wave critical Zeeman field shifts to more negative values as $T$ rises, and at small $\\Omega$ the direct stripe-to-zero-momentum boundary also shifts downward.","Quantum and thermal fluctuations push the stripe-to-plane-wave boundary in opposite directions, so zero-point motion strengthens supersolidity while heat destroys it."],"supporting_citations":[{"why":"Supplies the spin-1 spin-orbit-coupled Hamiltonian and the experimental Raman-coupling and quadratic-Zeeman tuning knobs used throughout.","marker":"[4]"},{"why":"Introduces the roton-gap-closing method for locating the stripe-to-plane-wave boundary from the plane-wave excitation spectrum.","marker":"[28]"},{"why":"Provides the zero-temperature phase diagram and the roton and double-roton excitation structure for spin-1 SOC BECs that this paper extends to finite temperature.","marker":"[29]"},{"why":"Gives the BdG ground-state analysis and spinor forms for interacting spin-orbit-coupled spin-1 condensates used in the numerical scheme.","marker":"[30]"},{"why":"Supplies zero-temperature elementary-excitation results and critical couplings for spin-1 SOC BECs that serve as the $T=0$ baselines.","marker":"[34]"},{"why":"Reports the experimentally measured finite-temperature phase diagram of a spin-orbit-coupled BEC, the melting phenomenon this paper models theoretically.","marker":"[37]"},{"why":"Applies HFB-Popov theory to the finite-temperature phase diagram of a Raman spin-orbit-coupled Bose gas, the approach directly extended here to spin-1.","marker":"[39]"},{"why":"Provides the HFB-Popov treatment of quantum and thermal fluctuations for a spin-1 condensate without spin-orbit coupling, generalized here to the Raman-coupled case.","marker":"[43]"}],"fun_headline_variants":["Heat shrinks supersolid stripe in spin-1 condensate","Thermal fluctuations melt supersolid stripe in spin-1 BEC","Roton gap opens, stripe phase melts at finite T","Finite-T maps show stripe phase melting in spin-1 gas","Stripe phase loses ground to plane wave as T rises"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The phase boundaries are read off from where the roton gap closes in the plane-wave (or zero-momentum) phase, which assumes that this softening, not a first-order jump or stripe-phase fluctuations, is what actually marks the transition at finite temperature.","fun_headline_variants_meta":{"raw":{"variants":["Heat shrinks supersolid stripe in spin-1 condensate","Thermal fluctuations melt supersolid stripe in spin-1 BEC","Roton gap opens, stripe phase melts at finite T","Finite-T maps show stripe phase melting in spin-1 gas","Stripe phase loses ground to plane wave as T rises"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000192,"raw_usage":{"total_tokens":1349,"prompt_tokens":949,"completion_tokens":400,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":565,"completion_tokens_details":{"reasoning_tokens":314}},"tokens_in":565,"tokens_out":400,"duration_ms":4587,"temperature":1.0,"reasoning_tokens":314,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T00:44:46.262856+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct finite-temperature free-energy comparison of the stripe and plane-wave states would settle the point: if the coexistence value of $\\Omega$ or $\\epsilon$ differs from the roton-gap closing value, or if the transition shows hysteresis, the roton-closing criterion used to draw the phase boundaries is incomplete.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the spin-1 spin-orbit-coupled Hamiltonian and the experimental Raman-coupling and quadratic-Zeeman tuning knobs used throughout."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the roton-gap-closing method for locating the stripe-to-plane-wave boundary from the plane-wave excitation spectrum."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the zero-temperature phase diagram and the roton and double-roton excitation structure for spin-1 SOC BECs that this paper extends to finite temperature."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the BdG ground-state analysis and spinor forms for interacting spin-orbit-coupled spin-1 condensates used in the numerical scheme."},{"cited_title":"Phys.24 073041","cited_arxiv_id":null,"evidence_quote":"Supplies zero-temperature elementary-excitation results and critical couplings for spin-1 SOC BECs that serve as the $T=0$ baselines."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports the experimentally measured finite-temperature phase diagram of a spin-orbit-coupled BEC, the melting phenomenon this paper models theoretically."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Applies HFB-Popov theory to the finite-temperature phase diagram of a Raman spin-orbit-coupled Bose gas, the approach directly extended here to spin-1."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the HFB-Popov treatment of quantum and thermal fluctuations for a spin-1 condensate without spin-orbit coupling, generalized here to the Raman-coupled case."}],"review_version":1}