{"id":"c346d6ae-d817-4204-8d9f-e93efe0550f6","arxiv_id":"2506.10649","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In the gapless superfluid phase, the zero-temperature normal-fluid fraction of nucleon superfluids is finite and can approach unity at the critical superfluid velocity.","lead":"A neutron star's superfluid can, in a 'gapless' state, develop a normal-fluid component that carries mass and momentum even at absolute zero, not only at finite temperature. The authors derive an exact formula for this component and argue it could change how neutron-star glitches and oscillations are modeled.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (38) is algebraically sound, but the central claim assumes the homogeneous gapless superflow is a stable equilibrium; negative-energy quasiparticles make phase separation or LOFF a plausible lower-energy competitor, and no stability test is provided.","rationale":"I found no algebraic error in the derivation of Eq. (38); the T→0 reduction from Eq. (32) is internally consistent and reproduces the expected zero at the Landau velocity and unity at the critical velocity in the weak-coupling limit. The paper's own Section IV candidly describes the pinning-force and critical-lag uncertainties, and the reader's conditional verdict already rests on the unverified reachability of the gapless regime. My stress-test refines this: before worrying whether stellar vortices can reach V>VL, one must know whether the homogeneous gapless state itself is the stable solution at fixed density and current. Because the quasiparticle spectrum has a negative-energy shell, the uniform state is a candidate for phase separation or LOFF-type modulation, and no calculation in this manuscript rules those out. The negative proton normal-fluid density in Fig. 5 does not by itself invalidate the formalism—normal/superfluid densities are response functions and need not be positive—but it reinforces that the 'normal fluid mass' should be tied to a well-defined thermodynamic state. The concrete test proposed above would settle the stability issue. Since the paper is already CONDITIONAL and the concern lands within the stated condition, no verdict change is warranted.","tokens_in":16857,"tokens_out":13333,"duration_ms":167390,"concrete_test":"Fix a baryon density n and a total current density J whose associated effective velocity lies in the gapless interval VLq < Vq < Vcq. Solve Eqs. (43)–(44) for the uniform gapless solution and evaluate its T=0 grand potential Ω = E − μN − V·J (equivalently free energy at fixed J). Then evaluate Ω for a phase-separated configuration with a fraction f of the volume in the subcritical superfluid at V=VLq and the remainder in the normal phase, and for a LOFF/FF ansatz of the form Δ(r)=Δ_q e^{2iq·r} with the same average density and current, optimizing q. If either competitor has lower Ω, the uniform gapless state is thermodynamically unstable; Eq. (38) would then not be an equilibrium normal-fluid fraction, and the paper's central astrophysical claim would need to be withdrawn or reframed as a metastable-state result.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The transfer from the T→0 limit of Eq. (32) to Eq. (38) checks out; the weak-coupling benchmark (52) and the vanishing at the Landau velocity are consistent. The load-bearing premise is the one stated in Sec. II.A and taken from Ref. [24]: for VLq < Vq < Vcq the stationary homogeneous solution exists with quasiparticle energies E_{k}^{(q)} = ℏk·Vq + sqrt(ε_k^2 + Δ_q^2) negative on the momentum shell x− < x < x+ (Eq. (35)). At T=0 these negative-energy levels are occupied, and it is this occupation that produces the nonzero normal-fluid fraction quantified by Eq. (38). But a negative-energy pocket is exactly the situation in which the uniform state can lower its energy by phase separation into subcritical superfluid and normal regions, or by a LOFF/FF modulated condensate with a finite pairing wavevector. The manuscript does not compute the T=0 free energy of the gapless solution, does not examine the Hessian or the BdG spectrum including phase fluctuations, and does not compare with any competing ordered state; it inherits the existence of the state from the authors' earlier paper [24] without supplying the stability test. If the homogeneous gapless state is only a saddle point, Eq. (38) describes a transient configuration rather than the equilibrium normal-fluid fraction, and the neutron-star implications lose their target. The astrophysical pinning uncertainty discussed in Sec. IV is a second barrier, but the stability question is logically prior.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies hydrodynamic properties of homogeneous neutron-proton superfluid mixtures in the gapless superflow regime, building on earlier work within time-dependent nuclear energy-density functional theory. The main formal result is a closed-form expression, Eq. (38), for the quantity Y_q that controls the normal-fluid densities and the entrainment matrix at zero temperature when the effective superfluid velocity lies between Landau's velocity and the critical velocity. The authors also derive a weak-coupling universal formula, solve the self-consistent gap and density equations in the gapless phase, and apply the formalism to npeμ matter in neutron-star cores using the BSk24 functional. They conclude that a zero-temperature normal-fluid component can appear in the gapless phase and discuss whether such a phase could be realized in neutron stars.","tokens_in":17179,"tokens_out":14642,"duration_ms":186544,"significance":"If the homogeneous gapless superflow is a genuine equilibrium state, the paper provides a valuable analytic input for two-fluid hydrodynamics of neutron-star matter. The derivation is self-contained and reduces to known weak-coupling limits (Leggett, Vollhardt-Maki), which is a strength. The numerical application self-consistently solves the gap and density equations rather than fitting the normal-fluid fractions, and the predicted weak-coupling universality of Y_q is a clear, falsifiable statement. However, the central physical interpretation depends on the stability of the gapless state, which the manuscript does not establish. The result may describe a stationary but unstable configuration, in which case the astrophysical implications would need substantial revision. The paper is therefore significant in its formal content but currently conditional in its physical claim.","major_comments":[{"comment":"The central claim that a normal-fluid component exists at T=0 in the gapless phase presupposes that the homogeneous stationary solution with negative-energy quasiparticle states is the equilibrium state. The paper does not test this premise: it does not compute the T=0 free energy of the gapless solution, does not examine the Hessian or the Bogoliubov-de Gennes spectrum including phase fluctuations, and does not compare with phase-separated or LOFF/FF modulated states. The occupation of the negative-energy pocket x- < x < x+ is precisely the mechanism producing Y_q > 0 in Eq. (38), so an instability of the homogeneous state would invalidate the interpretation of Eq. (38) as the equilibrium normal-fluid fraction. This is load-bearing for the neutron-star implications in Sections III and IV; the pinning uncertainties discussed in Section IV are secondary to this question. Please add an explicit stability analysis, or at the very least state clearly that the results apply only to metastable configurations and explain what observable consequences would survive in that case.","section":"Sec. II.A, Eq. (35), and Sec. IV"}],"minor_comments":[{"comment":"As printed, the asymptotic expression for the case E_x > 2 V̄_q sqrt(x) appears dimensionally inconsistent and inconsistent with the subsequent derivation: the logarithm should behave as 2 V̄_q sqrt(x)/T̄_q, not as 2 V̄_q E_x/T̄_q. The derivation of Eqs. (43) and (44) is consistent with the former form, so this is likely a typesetting error, but it should be corrected.","section":"Eq. (42)"},{"comment":"The symbol Y_p is overloaded: in Section II, Y_p denotes the response function defined in Eq. (25), while at the beginning of Section III.A it is redefined as the proton fraction ρ_p/ρ. This makes Eqs. (59)-(63) hard to read, since the same letter appears with two different meanings in the same formulas. Please use a distinct symbol such as x_p for the proton fraction.","section":"Sec. III.A, Eqs. (59)-(63)"},{"comment":"The statement that no approximation has been made so far is strong but accurate only within the TDHFB model and in the thermodynamic limit. The continuum limit and the T→0 limit are mathematical limits, but the sentence could be made less ambiguous by specifying that the result is exact within the chosen theoretical framework.","section":"Sec. II.F"}],"recommendation":"major_revision","confidential_remarks":"The formal derivation appears sound, but the physical premise of a stable homogeneous gapless superflow is unverified. A stability analysis, even a local one or a simple LOFF/phase-separation comparison, would substantially strengthen the paper. Without it, the central astrophysical conclusion remains conditional. I would encourage the editor to seek such an addition before reconsidering the manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, if you care about two-fluid neutron-star models, read this. Equation (38) is an exact closed form for the T=0 normal-fluid fraction Y_q in the gapless phase, and it is derived properly: direct integration of Eq. (32), with the right limits to Leggett, Vollhardt–Maki, and Parmenter. That is a genuinely new result, not a repackaging. The numerical work solves the self-consistent gap and density equations rather than fitting the target, and the authors are candid that the astrophysical realization depends on pinning forces uncertain by orders of magnitude. Credit is due for that honesty and for the formal care.\n\nThe soft spot is the load-bearing premise. The gapless phase is assumed to exist as a stationary homogeneous solution for VL < V < Vc, inherited from their earlier paper. But the negative-energy quasiparticle pocket that produces the normal fluid is exactly the kind of configuration that can be unstable to phase separation or a LOFF-type modulated condensate. The paper does not compute the T=0 free energy, does not examine the BdG spectrum against phase fluctuations, and does not compare with any competing ordered state. If the homogeneous gapless solution is only a saddle point, Eq. (38) describes a transient, not an equilibrium normal-fluid fraction. That is logically prior to the pinning question, and the manuscript leaves it unaddressed. A secondary concern: the entrainment matrix from Refs. [18,19] is used at T=0 with negative-energy quasiparticles; a short justification or re-derivation would strengthen that step. The citation pattern is otherwise fine—the self-citations point to the framework this paper extends.\n\nWho gets value: specialists in neutron-star superfluidity and glitch modeling. If the gapless state is stable, this changes glitch reservoir estimates, mutual friction, and oscillation modes. If not, the formula is still a useful limiting calculation, but the astrophysical implications are moot. This paper deserves a serious referee. Send it to review, but require the authors to either demonstrate stability of the homogeneous gapless state or explicitly frame the results as conditional on an untested assumption.","headline":"The math in Eq. (38) is real and clean; the physical conclusion hangs on an untested stability assumption that the paper does not address.","tokens_in":17751,"tokens_out":3735,"would_cite":true,"duration_ms":47048,"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":"Inside cold neutron stars, a superfluid pushed past Landau's velocity can develop a normal-fluid component that carries mass even at zero temperature, so two-fluid models must include a zero-temperature normal part in the gapless regime.","keywords":["gapless superfluidity","neutron stars","normal-fluid fraction","two-fluid model","superfluid velocity","Landau critical velocity","nuclear energy-density functional","pulsar glitches"],"falsifier":"Compute the zero-temperature phase stability of a homogeneous superflow in the interval between Landau's and critical velocities: if a spatially modulated state has lower free energy, the predicted normal fraction does not exist as a homogeneous phase. Observationally, a firm upper bound on vortex unpinning lags below Landau's velocity would show that neutron-star vortices never enter the gapless regime.","tokens_in":16613,"feed_emoji":"🌌","tokens_out":6246,"duration_ms":70704,"temperature":0.7,"pith_summary":"Superfluid helium acquires a normal component only by heating; a cold neutron star's superfluid was thought to be entirely superfluid at zero temperature. This paper argues that in the gapless superfluid phase, which sets in when the local superfluid velocity lies between Landau's velocity and a higher critical velocity, quasiparticle excitations appear even at $T=0$. Those excitations carry mass and momentum but no entropy, so the normal-fluid mass density of the star is nonzero at absolute zero. The authors derive an exact formula for the normal-fluid fraction of an arbitrary neutron-proton superfluid mixture and give a universal weak-coupling expression. If realized in nature, the result means the hydrodynamics of cold neutron stars has an extra mass-carrying component that current models omit.","feed_headline":"Neutron-star superflow can carry normal mass even at absolute zero","feed_subtitle":"A gapless superfluid adds a mass-carrying normal component at T=0 and shrinks the reservoir behind pulsar glitches.","key_machinery":"The engine of the argument is the function $Y_q(T, \\mathbf{V}_q)$, defined as the ratio of the momentum carried by quasiparticle excitations of species $q$ to the total momentum of that species; when effective masses equal bare masses, $Y_q$ is exactly the normal-fluid fraction. In the gapless phase the paper reduces $Y_q$ at $T=0$ to a closed algebraic expression, Eq. (38), involving the reduced chemical potential $\\bar\\mu_q$, the reduced pairing field $\\bar\\Delta_q$, and the reduced effective superfluid velocity $\\bar V_q$. The interval boundaries are Landau's velocity $V_{Lq}$, where the quasiparticle energy gap closes while the order parameter is unchanged, and the critical velocity $V^{(0)}_{cq}$, where $\\Delta_q$ vanishes and superfluidity disappears. Through the entrainment matrix of the time-dependent Hartree-Fock-Bogoliubov theory, $Y_q$ feeds the normal-fluid densities $\\rho_N^{(n)}$ and $\\rho_N^{(p)}$; the weak-coupling form matches the function Vollhardt and Maki derived for superfluid $^3$He, which makes the neutron-star result a direct corollary of a known superfluid universality.","core_discovery":"The central claim is that within the gapless superfluid phase, defined by effective superfluid velocities $V_{Lq} < V_q < V^{(0)}_{cq}$ for species $q$, the normal-fluid mass density is nonzero even at zero temperature. For each nucleon species the normal-fluid fraction is controlled by the function $Y_q$, and the paper obtains the exact zero-temperature result $$Y_q = \\left( \\bar\\mu_q + \\bar $V_q^{2}$ - \\frac{1}{4} \\frac{\\bar\\$Delta_q^{2}$}{\\bar $V_q^{2}$} \\right)^{3/2}$$ in the gapless interval, with $Y_q$ rising from 0 at Landau's velocity to 1 at the critical velocity. In the weak-coupling limit this becomes a universal function of $V_q/V_{Lq}$, independent of composition, density, and the chosen nuclear energy-density functional. The authors show that the true velocity at which nucleons are transported differs from the superfluid velocity in this phase, and that in neutron-star matter the neutron normal-fluid fraction can be substantial while the proton one is negative, because normal-fluid densities are current-current response coefficients rather than literal particle densities. The paper concludes that the neutron superfluid reservoir in the outer core may be significantly reduced in the gapless phase, complicating the standard picture of pulsar glitches.","pith_inferences":["If the zero-temperature normal component exists, its mass-carrying but entropy-free character may produce a distinctive inertia for torsional oscillations and gravitational-wave modes that differs from either a fully superfluid or fully normal core; future asteroseismic modeling could look for that signature.","A direct stability calculation of the homogeneous zero-temperature gapless superflow against spatial modulation would settle a premise the present paper inherits: if a LOFF-like state has lower energy anywhere in the gapless interval, the predicted homogeneous normal fraction never appears.","The same $Y_q$ machinery could be transferred to ultracold Fermi gases with moving superflows, where the velocity interval between Landau's velocity and the critical velocity is experimentally accessible and the predicted normal fraction at $T\\sim 0$ could be measured directly.","The negative proton normal-fluid density underscores that 'normal-fluid fraction' is a response coefficient, not a literal particle count, so observational tests should target total mass currents rather than local particle densities."],"forward_implications":["Two-fluid hydrodynamic models of neutron-star cores must include a zero-temperature normal component whenever the local neutron superfluid velocity lies in the gapless interval.","The available superfluid neutron reservoir in the outer core is reduced relative to the standard picture, so the angular-momentum budget for pulsar glitches shrinks in the gapless phase.","Because $Y_q$ is universal in the weak-coupling limit, predictions for the normal-fluid fraction do not depend on the uncertain details of the nuclear functional, only on the ratio $V_q/V_{Lq}$.","In the gapless phase the velocity with which nucleons are actually transported is not the superfluid velocity, so simulations that identify the two will miscompute mass currents.","Protons co-moving with the normal fluid stay below Landau's velocity, so in the standard two-fluid star the gapless normal component is carried by neutrons."],"supporting_citations":[{"why":"Establishes the existence of the gapless nuclear superfluid phase and the closure of the quasiparticle energy gap at Landau's velocity.","marker":"[24]"},{"why":"Gives the exact entrainment matrix for arbitrary temperatures and effective superfluid velocities in TDHFB theory, from which the normal-fluid densities are obtained.","marker":"[18]"},{"why":"Provides the quasiparticle energies and the self-consistent gap and number equations, and defines the effective superfluid velocity used throughout.","marker":"[19]"},{"why":"Derives the temperature-velocity integral representation of $Y_q$ and the low-temperature asymptotic forms of the gap equations.","marker":"[20]"},{"why":"Supplies the two-fluid mass-current decompositions and the relations between normal densities and the entrainment matrix.","marker":"[33]"},{"why":"Gives Leggett's normal-fluid fraction for superfluid Fermi liquids, to which the present $Y_q$ reduces in the current-free limit.","marker":"[34]"},{"why":"Gives the Vollhardt-Maki function for superfluid $^3$He, which the weak-coupling form of $Y_q$ matches and which establishes the universality.","marker":"[35]"},{"why":"Provides the BSk24 energy-density functional used for the numerical neutron-star applications.","marker":"[37]"}],"fun_headline_variants":["Zero-temperature normal flow emerges in gapless neutron-star superfluids","Gapless superfluidity in neutron stars yields normal mass at T=0","Neutron-star gapless superfluid: normal fraction survives at absolute zero","Normal-fluid fraction persists at zero Kelvin in neutron-star gapless phase"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"For the zero-temperature normal component to appear in a real star, the local neutron superfluid velocity must actually be pushed past Landau's velocity before vortices unpin, and the gapless state must remain stable against decay into other phases.","fun_headline_variants_meta":{"raw":{"variants":["Zero-temperature normal flow emerges in gapless neutron-star superfluids","Gapless superfluidity in neutron stars yields normal mass at T=0","Neutron-star gapless superfluid: normal fraction survives at absolute zero","Normal-fluid fraction persists at zero Kelvin in neutron-star gapless phase"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00025,"raw_usage":{"total_tokens":1576,"prompt_tokens":993,"completion_tokens":583,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":609,"completion_tokens_details":{"reasoning_tokens":506}},"tokens_in":609,"tokens_out":583,"duration_ms":6928,"temperature":1.0,"reasoning_tokens":506,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T04:23:56.558463+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the zero-temperature phase stability of a homogeneous superflow in the interval between Landau's and critical velocities: if a spatially modulated state has lower free energy, the predicted normal fraction does not exist as a homogeneous phase. Observationally, a firm upper bound on vortex unpinning lags below Landau's velocity would show that neutron-star vortices never enter the gapless regime.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the two-fluid mass-current decompositions and the relations between normal densities and the entrainment matrix."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives Leggett's normal-fluid fraction for superfluid Fermi liquids, to which the present $Y_q$ reduces in the current-free limit."}],"review_version":1}