{"id":"f087179f-0d32-40c6-82c5-5555a25c4ef1","arxiv_id":"2603.23622","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":4,"one_line_summary":"Viscosity damps neutron-star radial modes on ms timescales, shifts frequencies by up to ~1% at ζ∼10^30 g/cm/s, produces overdamped modes above ∼10^31, and cannot stabilize unstable stars in Eckart or BDNK theory.","lead":"Viscosity damps radial oscillations of cold polytropic neutron stars on millisecond timescales and shifts their frequencies by up to about one percent at bulk viscosities around 10^30 g/cm/s. The work quantifies these shifts in both Eckart and causal BDNK hydrodynamics and shows viscosity cannot stop linear gravitational collapse.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"BDNK threshold-shift claim rests on sparse time-domain runs near criticality that may be contaminated by numerical viscosity and KO dissipation.","rationale":"The reader already isolates the precise soft spot (sparse BDNK time-domain evidence near criticality, idealized cold polytropes). The remainder of the central claim—ms damping, percent-level frequency shifts, overdamping at ζ ∼ 10^{31}, and Eckart’s inability to stabilize—is cross-validated by independent FD/TD methods, public code, two EOS, and three causal frames, and is consistent with the small-viscosity analytic results of Caballero & Yunes. No stronger internal inconsistency appears; the CONDITIONAL verdict is therefore left unchanged.","tokens_in":25848,"tokens_out":593,"duration_ms":17203,"concrete_test":"Re-run the BDNK constrained system (frame A) on a denser ϵ_c grid (step 0.0005 × 10^{15} g cm^{-3}) around 5.66–5.67 for both ˆζ=10^{-3} and ˆζ=0.01, extracting Im(ω) from exponential fits to |δu(t,R_*)| over T ≥ 50 ms at h = 5 m with KO coeff. halved; if the zero-crossing remains within 0.001 of the Eckart value 5.663 for both viscosities, the “slightly modifies” claim is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The strongest claim’s BDNK half—that viscosity “slightly modifies the threshold of collapse” yet still cannot stabilize an unstable inviscid star—rests on three time-domain evolutions of Eckart-eigenvector initial data at ˆζ=10^{-3} (Fig. 7, ϵ_c = 5.65, 5.66, 5.67 × 10^{15} g cm^{-3}). Unlike the Eckart sector, no frequency-domain eigenvalue problem is solved for BDNK; the constrained first-order system (App. C) requires KO dissipation (coeff. 0.2) and already exhibits a numerical-viscosity plateau in damping rates for ˆζ ≲ 10^{-3} (Fig. 3). A 0.001 shift in the critical density is therefore comparable to both the TOV termination tolerance (p = p_c × 10^{-6}) and residual numerical dissipation, so the reported modification could be an artifact rather than a genuine frame effect. The over-damping and percent-level shift results themselves remain robust (Eckart FD + TD cross-checks, two EOS, frame robustness for ζ ≲ 10^{30}).","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The paper studies linear radial oscillations of cold, polytropic, spherically symmetric neutron stars in two first-order viscous hydrodynamics frameworks: Eckart and BDNK. Using both frequency-domain eigenvalue methods (Eckart) and time-domain evolutions (Eckart and BDNK), it reports that bulk viscosity damps radial modes on millisecond timescales and produces fractional frequency shifts that grow with compactness and viscosity, reaching the percent level for the fundamental mode near ζ ∼ 10^30 g cm^{-1} s^{-1}. For ζ ≳ 10^31 g cm^{-1} s^{-1} the fundamental frequency vanishes (overdamped). Viscosity in Eckart theory leaves the linear collapse threshold unchanged (consistent with prior analytic work) while slowing the instability rate; numerical evidence is offered that BDNK viscosity likewise cannot stabilize an unstable inviscid star but slightly shifts the threshold. Public code and extensive cross-checks (perfect-fluid recovery of Kokkotas & Ruoff, frame robustness, independent residuals) are provided.","tokens_in":26201,"tokens_out":1287,"duration_ms":11852,"significance":"If the results hold, the work supplies concrete, quantitative targets for viscous asteroseismology with third-generation detectors: percent-level frequency shifts and millisecond damping at post-merger viscosities, plus the existence of arbitrarily low-frequency overdamped radial modes. The public NeutronStarOscillations.jl package, the Eckart frequency-domain formulation, and the demonstrated agreement between Eckart and BDNK at moderate viscosity are reusable assets. The collapse analysis extends recent analytic stability criteria into the large-viscosity regime and supplies the first numerical indication of a BDNK threshold shift, even if that indication remains provisional.","major_comments":[{"comment":"Sec. V.C and Fig. 7: the claim that BDNK viscosity “slightly modifies the threshold of collapse” rests on three time-domain runs of Eckart-eigenvector initial data at ζ̂ = 10^{-3} near ε_c^* ≈ 5.663 × 10^{15} g cm^{-3}. No BDNK frequency-domain eigenvalue problem is solved; the constrained system (App. C) employs KO dissipation (coeff. 0.2) and already shows a numerical-viscosity plateau for ζ̂ ≲ 10^{-3} (Fig. 3). A 0.001 shift in critical density is comparable to the TOV termination tolerance (p = p_c \times 10^{-6}) and residual numerical dissipation. Either a systematic BDNK eigenvalue scan (or a carefully controlled resolution study that isolates the threshold) is needed, or the claim should be rephrased as a tentative indication pending further work.","section":"Sec. V.C, Fig. 7"},{"comment":"Secs. II–III, Eq. (23): all quantitative results (percent-level shifts, overdamping at ζ ∼ 5 \times 10^{31}, collapse timescales) are obtained for two fixed cold polytropes with η = ζ/10 and zero heat conductivity. While the authors note this limitation, the abstract and conclusions present the numbers as generic for neutron-star viscosities. A short discussion of how the quoted thresholds and fractional shifts are expected to change under finite-temperature or tabulated EOS would strengthen the central claim.","section":"Secs. II–III, Eq. (23)"}],"minor_comments":[{"comment":"Fig. 3 caption and surrounding text: the faint plateau of the BDNK curves as ζ̂ \to 0 is correctly attributed to numerical viscosity, but the figure itself would benefit from an explicit annotation or a higher-resolution inset so that readers do not misread the plateau as a physical effect.","section":"Fig. 3"},{"comment":"Appendix B, Tables II–III: the perfect-fluid frequencies are stated to agree with Kokkotas & Ruoff (2001) to ≲ 1 %. Quoting the absolute differences (or a short comparison column) would make the validation more transparent.","section":"Appendix B"},{"comment":"Eq. (9) and the definition of L: the length scale that converts dimensionless transport coefficients into dimensionful viscosities is never given a concrete numerical value. Stating the choice used for the reported ζ_c values would aid reproducibility.","section":"Eq. (9)"},{"comment":"Sec. IV.B.2: the Crank–Nicholson scheme and KO coefficient 0.2 are mentioned, but the precise form of the KO operator (and whether it is applied to all variables) is left implicit. A one-sentence clarification would help readers re-implement the code.","section":"Sec. IV.B.2"},{"comment":"Typographical: “Einstein-NA VIER-STOKES” (Sec. II heading) and occasional missing spaces around “ζ∼” should be cleaned.","section":"Sec. II"}],"recommendation":"minor_revision","confidential_remarks":"The BDNK threshold claim is the only load-bearing point that is currently under-supported; once it is either strengthened or carefully caveated, the paper is a solid contribution. The public code and the Eckart–BDNK cross-checks are genuine strengths that raise the manuscript above a pure methods note."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The new numbers here are the ones that matter: viscosity damps radial modes on ~ms timescales and shifts the fundamental frequency by up to ~1% at ζ ~ 10^30 g/cm/s, with the mode going overdamped near 10^31. Eckart and BDNK agree closely in that regime, and Eckart cannot move the collapse threshold (matching the recent analytic small-viscosity result). That package is useful for anyone thinking about viscous asteroseismology with 3G detectors.\n\nWhat the paper does well is the cross-check. Frequency-domain and time-domain Eckart agree; perfect-fluid limits recover Kokkotas & Ruoff to ≲1%; three BDNK frames show frame robustness for the shifts and damping; independent residuals converge; code is public. The overdamping transition appears for both polytropes they tried, which is a clean numerical fact even if the microphysics is idealized. Tables of modes are there for reuse.\n\nSoft spots are real but limited. Everything is cold barotropic polytropes, linear order, spherical symmetry, η = ζ/10, zero heat conductivity. That is the stated scope, not a hidden flaw. The BDNK threshold shift is the weakest piece: it rests on a few time-domain runs of Eckart-eigenvector initial data near criticality, with KO dissipation and a known numerical-viscosity plateau at low ζ. A 0.001 shift in ε_c is comparable to their TOV termination and residual dissipation, so I would treat “slightly modifies” as suggestive rather than established. The paper itself phrases it as numerical evidence, which is fair. The main claims on damping, percent-level shifts, and Eckart stability do not depend on that point.\n\nThis is for people working on relativistic viscous hydro or NS mode extraction. It is careful enough and reproducible enough that a serious editor should send it to referees. I would cite the shift/damping/overdamping numbers and the Eckart threshold result; I would not lean hard on the BDNK threshold modification without more work. Engage with it.","headline":"Solid linear numerics on viscous radial modes with public code; the percent-level shifts and overdamping are robust, while the BDNK threshold claim is thinner than the rest.","tokens_in":26806,"tokens_out":545,"would_cite":true,"duration_ms":6936,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Viscosity damps neutron-star radial modes in milliseconds and can erase their oscillation frequency at high bulk viscosity, but does not stop gravitational collapse.","keywords":["neutron stars","radial oscillations","bulk viscosity","Eckart hydrodynamics","BDNK hydrodynamics","gravitational collapse","asteroseismology","polytropic stars"],"falsifier":"A full frequency-domain eigenvalue scan of BDNK stars across a dense grid of central densities and viscosities, or a nonlinear radial simulation that shows whether an inviscid-unstable configuration remains unstable once finite-amplitude and finite-temperature effects are restored.","tokens_in":26765,"feed_emoji":"🌊","tokens_out":657,"duration_ms":6069,"temperature":0.7,"pith_summary":"Neutron-star oscillation modes are a prime target for future gravitational-wave detectors because they carry information about the dense nuclear matter inside the star. This paper asks what bulk viscosity does to the purely radial modes of cold, polytropic stars. Working to linear order in two first-order relativistic hydrodynamics theories—one acausal (Eckart) and one causal (BDNK)—the authors show that viscosity damps the modes on millisecond timescales and shifts their frequencies by up to the percent level at bulk viscosities of order 10^30 g/cm/s. At still higher viscosity the fundamental frequency falls to zero and the mode becomes overdamped. Viscosity cannot stabilize a star that is already unstable to collapse, though it can slow the collapse rate dramatically and, in the causal theory, slightly moves the critical density. The results supply concrete numbers for viscous asteroseismology with next-generation detectors.","feed_headline":"Viscosity damps star quakes in ms, can kill their frequency","feed_subtitle":"Percent-level shifts at merger viscosities; overdamping above 10^31; collapse still wins","key_machinery":"Linearized radial master equations (or constrained wave-plus-constraint systems) for the Eckart and BDNK stress-energy tensors, solved both as frequency-domain eigenvalue problems and as time-domain evolutions of single-mode and Gaussian initial data on polytropic TOV backgrounds.","core_discovery":"Viscosity damps radial modes of cold polytropic neutron stars on millisecond timescales and produces fractional frequency shifts that grow with both compactness and viscosity, reaching the percent level for the fundamental mode near ζ ∼ 10^30 g/cm/s; for ζ ≳ 10^31 g/cm/s the frequency vanishes (overdamped). Viscosity in Eckart theory leaves the linear collapse threshold unchanged; numerical evidence indicates BDNK viscosity is likewise unable to prevent collapse while only slightly shifting the threshold.","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["Viscosity damps neutron-star radial modes in milliseconds","Percent-level frequency shifts from bulk viscosity in compact stars","Modes overdamp for ζ ≳ 10^31; collapse threshold barely shifts","Eckart and BDNK viscosity leave unstable stars free to collapse","Radial oscillations silenced by viscosity on ms timescales"],"cache_read_input_tokens":16512,"weakest_assumption_plain":"The stars are treated as cold barotropic polytropes with fixed shear-to-bulk ratio and zero heat conductivity, and all statements about stability are made only at linear order in spherical symmetry.","fun_headline_variants_meta":{"raw":{"variants":["Viscosity damps neutron-star radial modes in milliseconds","Percent-level frequency shifts from bulk viscosity in compact stars","Modes overdamp for ζ ≳ 10^31; collapse threshold barely shifts","Eckart and BDNK viscosity leave unstable stars free to collapse","Radial oscillations silenced by viscosity on ms timescales"]},"model":"grok-4.5","effort":"low","cost_usd":0.003648,"raw_usage":{"total_tokens":1254,"prompt_tokens":873,"num_sources_used":0,"completion_tokens":90,"cost_in_usd_ticks":36480000,"prompt_tokens_details":{"text_tokens":873,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":291,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":873,"tokens_out":90,"duration_ms":2748,"temperature":1.0,"reasoning_tokens":291,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-13T19:27:32.873188+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"A full frequency-domain eigenvalue scan of BDNK stars across a dense grid of central densities and viscosities, or a nonlinear radial simulation that shows whether an inviscid-unstable configuration remains unstable once finite-amplitude and finite-temperature effects are restored.","supporting_citations":[],"review_version":1}