{"id":"00c0eec8-6960-4c8c-9c7b-ff19f88cdab7","arxiv_id":"2505.00539","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"The quark mean-field model predicts larger crust clusters and a slower direct-Urca cooling phase than the relativistic mean-field model, while both reproduce the observed crustal cooling of the transient KS 1731-260.","lead":"This paper builds a family of neutron star equations of state using a quark-based nuclear model and compares them against the standard hadron-based model at three values of the symmetry energy slope. It then predicts how those interiors change neutron star radii, cooling histories, and the onset of r-mode spin instability.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Longer QMF dUrca relaxation time rests on volume fraction alone; local emissivity from higher effective mass and shared pairing gaps could reverse the ordering.","rationale":"I read the paper in good faith: it constructs genuinely unified QMF EoSs, directly computes crust composition, M-R relations, cooling curves, and r-mode damping, and it flags its own main limitations. The strongest-claim dynamics, however, rest on Sec. III B's inference from dUrca core volume fractions. The reader identified shared pairing gaps as the weakest assumption; I agree that this is a real limitation, and the paper itself acknowledges it. My stress-test adds a second, arguably more direct route to the same vulnerability: the paper compares f_dUrca_V but never integrates the dUrca emissivity over the two cores. Since QMF has higher effective masses in the high-density region, local dUrca emissivity is model-dependent in the opposite direction from volume, and the total cooling power could be comparable. This means the longer-relaxation claim is not robustly established even before pairing gaps are varied. I do not see an internal inconsistency or a reason to reject the paper; the static crust and M-R results are computed rather than assumed, and the limitation is disclosed. The quantitative dynamic prediction should be read as conditional pending a volume-integrated emissivity check and ideally self-consistent pairing gaps. My recommendation therefore preserves the reader's CONDITIONAL verdict, so the verdict_should_be field is UNCHANGED.","tokens_in":25841,"tokens_out":5546,"duration_ms":59833,"concrete_test":"Recompute the volume-integrated direct Urca neutrino luminosity Q_dUrca = integral of q_dUrca(rho, M*_n, M*_p, Y_p) dV for the QMF and RMF L0=80 EoSs at M = 1.4 and 1.8 Msun, using the same pairing-gap inputs as Sec. II E and the standard Lattimer et al. (1991) emissivity. Then rerun the NSCool cooling curves with this luminosity rather than using f_dUrca_V as a proxy. If QMF's total dUrca luminosity is not smaller than RMF's, or if the luminosity-drop epoch is not later, the longer-relaxation claim is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Sec. III B's headline dynamical result — that QMF predicts a longer dUrca thermal relaxation time than RMF — is inferred from dUrca core volume fractions (f_V = 2.6% vs 7.8% for 1.4 Msun, and 24.7% vs 33.7% for 1.8 Msun), with the statement that smaller dUrca cores give longer relaxation times (Refs. [79,80]). The paper never computes the total dUrca luminosity. Direct Urca neutrino emissivity scales with effective masses (roughly M*_n M*_p; see Ref. [52]) and with the dUrca-allowed volume. At high density QMF has M*/M_N ≈ 0.5 while RMF drops monotonically (Fig. 2a); a factor ~1.4 in M* alone gives a factor ~2 in local dUrca emissivity, partially offsetting the factor ~3 smaller volume. The paper only bounds the mUrca+bremsstrahlung emissivity ratio (<2) and PBF ratio (<1.15), not the dUrca emissivity ratio across the dUrca-allowed region. It also adopts identical neutron 3P2 and proton 1S0 critical temperatures in both cores (Sec. II E and III B), an approximation explicitly stated in Sec. IV to lie beyond the present scope; if QMF gaps differ, specific heat and dUrca suppression change, and the relaxation-time ordering could be reduced or reversed. Thus the longer-relaxation claim is not yet established by the presented evidence.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs unified equations of state (EoSs) within the quark mean-field (QMF) framework for three values of the symmetry energy slope L0 (40, 60, 80 MeV), with the crust and core described by the same effective Lagrangian and self-consistent compositions. The QMF unified EoSs are compared with analogous relativistic mean-field (RMF) EoSs through neutron star mass-radius relations, crust compositions (including pasta phases), thermal cooling of isolated neutron stars and of the X-ray transient KS 1731-260, and r-mode instability windows. The main reported results are that the QMF model predicts heavier nuclear clusters and larger Wigner-Seitz cells in the crust, a higher maximum mass, and, when direct Urca cooling is active (L0 ≳ 80 MeV), a longer thermal relaxation time than the RMF model, while both models can reproduce the observed crustal cooling curve of KS 1731-260. The paper also provides power-law fits for the L0 dependence of bulk and shear viscosities at saturation density.","tokens_in":26248,"tokens_out":3449,"duration_ms":36417,"significance":"If the main claims hold, the paper offers a useful, internally consistent set of unified EoSs that link quark-level input to neutron-star observables, including static structure, cooling, and r-mode stability. The systematic comparison between QMF and RMF at fixed saturation properties isolates the effect of the underlying many-body framework, which is valuable for the dense-matter community. Concrete strengths are the fully self-consistent crust-core construction (except for the quark-level caveat discussed below), the explicit comparison with an approximate unified-EoS method, the Markov-chain Monte Carlo fits to KS 1731-260, and the quantitative bounds on emissivity ratios for non-direct-Urca processes. However, the headline dynamical claim about longer direct-Urca relaxation time rests on volume fractions alone rather than on integrated emissivities, and the adopted approximation of identical pairing gaps in QMF and RMF cores weakens that claim. The paper is therefore a solid contribution that needs additional quantitative support for its central dynamical conclusion.","major_comments":[{"comment":"The cooling simulations assume identical neutron 3P2 and proton 1S0 critical temperatures in the QMF and RMF cores, even though these pairing gaps depend on effective mass and proton fraction, both of which differ between the models. The paper acknowledges in Sec. IV that a fully self-consistent determination lies beyond the present scope. Because the dUrca relaxation time is sensitive to superfluidity (through suppression of neutrino emissivity and modification of specific heat), the predicted longer QMF relaxation time could be reduced or reversed if the QMF gaps differ from the RMF gaps. At minimum, please provide a sensitivity estimate (e.g., vary the adopted critical temperatures within a reasonable range and show how the cooling curves and relaxation times change) or quantify the expected difference in the core pairing gaps from the differences in effective mass and proton fraction.","section":"Sec. II E and Sec. III B"}],"minor_comments":[{"comment":"The text says 'the QMF model predicts a longer thermal relaxation time than the RMF model' but does not define the quantitative measure of relaxation time (e.g., the time when the luminosity drop is steepest). Please specify the definition used in the context of the cooling curves.","section":"Sec. III B, Fig. 5"}],"recommendation":"major_revision","confidential_remarks":"The paper is a competent and useful contribution, but the 'unified from the quark level' claim is not supported by the crust calculation, and the direct-Urca relaxation-time ordering needs a full emissivity calculation. Both issues are fixable without changing the scope of the work. The reviewer is not concerned about the model-calibration workflow per se; the L0 dependence of the fits is legitimate. Please ensure the authors do not merely enclose the caveat in the conclusions but also revise the abstract and title-level claims accordingly."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is a solid, workmanlike paper that builds unified QMF equations of state with three values of L0 and compares them against RMF across structure, crust composition, cooling, and r-modes. That combination is genuinely new, and the paper does it in a disciplined way. The crust differences (heavier clusters, larger Wigner-Seitz cells, similar free-neutron density) are interesting; the M-R relations and the KS 1731-260 crust cooling fits are competently done; the viscosity power-law fits, though simple, will be convenient for people doing r-mode and Bayesian studies.\n\nThe soft spot is the one the stress-test note puts its finger on. The longer QMF thermal relaxation time in the dUrca case is inferred entirely from dUrca core volume fractions (2.6% vs 7.8% for 1.4 Msun). The paper never computes the total dUrca luminosity, and it explicitly bounds only the mUrca+bremsstrahlung and PBF emissivity ratios. Since QMF has a higher effective nucleon mass at high density, the local dUrca emissivity should be larger—roughly a factor of two from the effective masses alone—which partially cancels the factor of three in volume. That does not necessarily overturn the ordering, but it is not established by the evidence presented. The shared pairing-gap assumption in the core adds another layer of uncertainty; the authors acknowledge this and say a self-consistent treatment is beyond scope, which is honest but means the dynamical headline should be read as tentative.\n\nThe word \"unified\" also deserves a flag: the crust is derived from the same Lagrangian but not from the quark level, as the authors themselves admit in Sec. II B. That is a reasonable limitation, but it should be kept in mind when citing this as a fully unified QMF EoS.\n\nThese are not fatal flaws. The static structure, crust physics, and cooling comparisons hold up well; the dynamic dUrca result needs either a direct computation of the dUrca emissivity or a softened interpretation. I would bring this to a reading group and I would probably cite it for the crust and EoS comparisons. It deserves serious peer review—send it out, and ask the authors to address the dUrca luminosity point directly.","headline":"A careful, systematic QMF-vs-RMF EoS comparison with real new results, but the headline claim of longer QMF dUrca relaxation is plausible rather than proven because local dUrca emissivity and pairing gaps are not computed consistently.","tokens_in":26757,"tokens_out":1680,"would_cite":true,"duration_ms":20122,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A quark-level equation of state predicts heavier crust clusters, a higher maximum mass, and a longer direct-Urca cooling delay than the equivalent hadronic model.","keywords":["neutron star equation of state","quark mean-field model","unified EoS","neutron star cooling","direct Urca process","r-mode instability","symmetry energy slope","Wigner-Seitz cell"],"falsifier":"Recompute the proton $^1S_0$ and neutron $^3P_2$ critical temperatures from the Landau effective masses and proton fractions of the QMF and RMF unified EoSs and rerun the thermal evolution; if the QMF star no longer shows a longer direct-Urca relaxation time than the RMF star at the same $L_0=80$ MeV, the paper's central cooling distinction fails.","tokens_in":25680,"feed_emoji":"⭐","tokens_out":6785,"duration_ms":62279,"temperature":0.7,"pith_summary":"This paper constructs unified neutron-star equations of state from the quark mean-field (QMF) model, calibrated to three values of the symmetry-energy slope ($L_0 = 40, 60, 80$ MeV), and compares them with hadronic relativistic mean-field (RMF) models fitted to the same saturation properties. It seeks to show that the quark-level description changes the coherent picture of neutron-star structure: the QMF crust contains heavier nuclear clusters and larger Wigner-Seitz cells, the maximum mass is higher by about $0.13\\,M_\\odot$, and when the direct Urca process is open the star cools through a smaller rapid-cooling core with a longer thermal relaxation time. The paper also claims that both QMF and RMF reproduce the observed crustal cooling of the transient KS 1731-260, and that a larger $L_0$ widens the r-mode instability window while a larger stellar mass narrows it. The value of this claim is that static and dynamical neutron-star observables are traced back to one internally consistent microscopic input.","feed_headline":"Quark mean-field EoS delays rapid neutron-star cooling","feed_subtitle":"One microscopic model now links crust composition, maximum mass, cooling, and r-mode stability in one picture.","key_machinery":"The central object is the in-medium nucleon mass $M_N^*$ of the quark mean-field model, obtained by solving the Dirac equation for constituent quarks confined by a harmonic-oscillator potential in the presence of $\\sigma$, $\\omega$, and $\\rho$ meson fields, then adding center-of-mass, pion-cloud, and one-gluon-exchange corrections. This density-dependent mass controls the proton fraction, the direct Urca threshold, neutrino emissivities, and the viscous damping that sets the r-mode window. The crust is built in the same framework by solving coupled Klein-Gordon and Poisson equations in Wigner-Seitz cells under the Thomas-Fermi approximation, with five pasta geometries, and the core is joined to it at the crust-core transition density. The paper also supplies power-law parameterizations $\\log_{10}\\xi = 12.9285\\,L_0^{0.0767}$ and $\\log_{10}\\eta = 17.5192\\,L_0^{-0.0028}$ for the bulk and shear viscosities at saturation density.","core_discovery":"The central claim is that the QMF model, because its in-medium nucleon effective mass $M_N^*$ approaches roughly $0.5M_N$ at high density instead of declining monotonically, suppresses the proton fraction in the dense core relative to the RMF prediction. That suppression pushes the onset density of the direct Urca process upward, so for $L_0 = 80$ MeV the QMF star has a smaller dUrca-allowed core (volume fractions $2.6\\%$ versus $7.8\\%$ at $1.4\\,M_\\odot$, and $24.7\\%$ versus $33.7\\%$ at $1.8\\,M_\\odot$) and a correspondingly longer thermal relaxation time in the cooling curve. In the crust, the stronger $\\sigma$-meson field of QMF produces heavier nuclear clusters, larger Wigner-Seitz cells, and a lower outer-inner crust transition density, while the free neutron gas density remains nearly unchanged because the larger cells compensate. The paper argues that these differences leave cooling without dUrca almost model-independent, that both models fit the KS 1731-260 crustal cooling data with slightly different core temperatures and impurity parameters, and that the r-mode instability window widens with $L_0$ and narrows with stellar mass.","pith_inferences":["Because the cooling comparison assumes identical nucleon pairing gaps in QMF and RMF cores, a self-consistent gap calculation could shorten or reverse the predicted longer direct-Urca relaxation time; that is the natural next test.","The QMF suppression of the high-density proton fraction would also affect other proton-fraction-sensitive processes the paper does not model, such as ambipolar diffusion in magnetar fields or vortex pinning.","The saturation-density viscosity power laws are immediately usable in population-synthesis studies of r-modes, without recomputing the full unified EoS.","The near-degeneracy of QMF and RMF cooling curves when direct Urca is absent suggests that distinguishing the two models observationally requires either a star with $L_0$ large enough to open dUrca or very precise crust-cooling measurements."],"forward_implications":["If the QMF unified EoS is correct, neutron stars of the same mass are slightly larger and their maximum mass is about $0.13\\,M_\\odot$ higher than the RMF prediction with the same saturation properties.","For $L_0 \\gtrsim 80$ MeV, the QMF model predicts that rapid cooling via direct Urca is delayed to later times; an observed late luminosity drop in a massive isolated neutron star would support the quark-level effective-mass behavior.","Crustal cooling transients such as KS 1731-260 do not cleanly discriminate between QMF and RMF; both reproduce the data with modest shifts in fitted core temperature and impurity parameter.","Larger symmetry-energy slope widens the r-mode instability window, while larger stellar mass narrows it, with the QMF model giving slightly lower critical frequencies and temperatures at high mass.","The approximate unified-crust treatment reproduces the full unified EoS for masses above about $0.5\\,M_\\odot$, so cheaper approximations can be used in parameter scans."],"supporting_citations":[{"why":"It provides the earlier unified-crust calculation that this paper extends to the QMF model.","marker":"[45]"},{"why":"It supplies the in-medium nucleon mass formula and its quark-level corrections used in Eq. (3).","marker":"[21]"},{"why":"It gives the approximate crust EoS method used for the comparison in Sec. II D.","marker":"[49]"},{"why":"It defines the direct Urca threshold and emissivity that drive the rapid-cooling differences.","marker":"[52]"},{"why":"It supplies the envelope Te-Tb relation and the neutron 3P2 critical temperature model adopted in the cooling simulations.","marker":"[59]"},{"why":"It provides the KS 1731-260 crustal cooling observations used for the Markov-chain fits.","marker":"[67]"},{"why":"It supplies the bulk viscosity expressions and the r-mode instability comparison baseline.","marker":"[70]"},{"why":"It provides the isolated neutron star luminosity data used to compare cooling curves.","marker":"[77]"}],"fun_headline_variants":["Quark mean-field EoS pushes neutron-star dUrca onset deeper","Quark EoS slows neutron-star cooling by suppressing protons","Unified quark model links neutron-star cooling and r-modes","Neutron-star cooling delayed by quark model's proton suppression"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the same nucleon pairing gaps (neutron $^3P_2$ and proton $^1S_0$) apply in both QMF and RMF cores, although these gaps depend on effective mass and proton fraction; the paper states that self-consistent critical temperatures are beyond its scope.","fun_headline_variants_meta":{"raw":{"variants":["Quark mean-field EoS pushes neutron-star dUrca onset deeper","Quark EoS slows neutron-star cooling by suppressing protons","Unified quark model links neutron-star cooling and r-modes","Neutron-star cooling delayed by quark model's proton suppression"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000273,"raw_usage":{"total_tokens":1743,"prompt_tokens":1163,"completion_tokens":580,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":779,"completion_tokens_details":{"reasoning_tokens":508}},"tokens_in":779,"tokens_out":580,"duration_ms":6374,"temperature":1.0,"reasoning_tokens":508,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T04:39:45.723663+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the proton $^1S_0$ and neutron $^3P_2$ critical temperatures from the Landau effective masses and proton fractions of the QMF and RMF unified EoSs and rerun the thermal evolution; if the QMF star no longer shows a longer direct-Urca relaxation time than the RMF star at the same $L_0=80$ MeV, the paper's central cooling distinction fails.","supporting_citations":[{"cited_title":"Exploring nuclear force with pulsar glitch observation","cited_arxiv_id":"2412.09219","evidence_quote":"It provides the earlier unified-crust calculation that this paper extends to the QMF model."},{"cited_title":"Neutron star equation of state: QMF modeling and applications","cited_arxiv_id":"2007.05116","evidence_quote":"It gives the approximate crust EoS method used for the comparison in Sec. II D."},{"cited_title":"Neutron star properties and the equation of state for its core","cited_arxiv_id":"1611.01357","evidence_quote":"It defines the direct Urca threshold and emissivity that drive the rapid-cooling differences."},{"cited_title":"Flowers, M","cited_arxiv_id":null,"evidence_quote":"It supplies the envelope Te-Tb relation and the neutron 3P2 critical temperature model adopted in the cooling simulations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the bulk viscosity expressions and the r-mode instability comparison baseline."}],"review_version":1}