{"id":"a2b6445c-5b78-49b8-8f2c-e34a4de0555e","arxiv_id":"1909.01278","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"low","formal_verification":"none","parameter_count":2,"one_line_summary":"Nucleon-nucleon-axion bremsstrahlung cooling in binary neutron star merger remnants is too weak to leave a detectable imprint on gravitational waves or ejecta, so it cannot improve axion mass constraints.","lead":"Simulations of binary neutron star mergers that include axion cooling show the extra cooling changes the remnant's temperature, but the effect on gravitational waves and ejected mass is too small to measure. The result suggests this particular axion channel cannot improve axion mass limits with current or planned detectors.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'best case' argument ignores EOS sensitivity: one EOS (ALF2) may not maximize the T^6/T^3.5 axion cooling rate, so the null result is not yet general.","rationale":"The paper is a careful null result with a transparent cooling implementation and constraint checks (Appendix B). The reader is right that ignoring reabsorption and neutrinos biases the result toward a larger axion effect, so those simplifications support the null result. The remaining gap is the single-EOS restriction. The steep temperature dependence of the emissivity means the cooling rate is not a small perturbation around ALF2; it can vary by orders of magnitude across plausible EOS and parameter choices. Thus the claimed 'best case' is only a best case within the absent reabsorption/neutrino assumptions, not within the EOS dimension. The authors' own Sec. V caveat admits this possibility. A single rerun with a hotter EOS would either validate the general conclusion or force a scoped claim, so this is a decisive, concrete test. The verdict remains near ACCEPT; CONDITIONAL reflects that the abstract should either be qualified to the tested EOS or supported by one additional EOS run.","tokens_in":12950,"tokens_out":15470,"duration_ms":165127,"concrete_test":"Repeat the highest-mass simulation (Case-1) with the same cooling implementation for ma = 1 eV versus ma = 0, but with a second, softer EOS known to produce hotter remnants (e.g., LS220 or SLy). Compare the postmerger GW frequency shift and collapse-time difference to the estimated 3G measurement uncertainties and to the EOS-systematic scatter. If the shifts remain below those thresholds, the central claim survives; if they exceed them, the conclusion must be restricted to the tested EOS and configurations.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's central claim that nucleon-nucleon-axion bremsstrahlung cannot improve axion-mass constraints is defended as a 'best case' scenario (Sec. II A, Sec. V). That defense covers two omissions that indeed suppress the cooling effect: free-streaming with no reabsorption and no neutrino cooling. It does not cover the choice of a single EOS (ALF2, Sec. III, Table I). The emissivity scales as epsilon_D proportional to T^6 and epsilon_ND proportional to T^3.5 (Sec. II A), so the cooling rate is highly sensitive to the postmerger temperature. A different plausible EOS or mass ratio can produce a hotter remnant, increasing the cooling rate by orders of magnitude. The authors explicitly concede in Sec. V that 'other equations of state ... might lead to larger, more noticeable effects.' Therefore the abstract's unqualified statement that a more thorough investigation is unlikely to change the conclusions is not supported by the presented evidence; the null result is established only for ALF2 equal-mass irrotational configurations. Without EOS sensitivity testing, the 'best case' argument is incomplete, and the central claim as worded is stronger than the simulations justify.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents the first numerical relativity simulations of binary neutron star mergers that include a phenomenological description of nucleon-nucleon-axion bremsstrahlung as an additional cooling channel. The cooling is implemented through a covariant radiation four-force with the Brinkmann-Turner emissivity, and a set of equal-mass irrotational ALF2 configurations is evolved for axion masses of 0, 1e-3, 1e-2, 0.1, and 1 eV. The authors find that axion cooling noticeably changes the temperature profile and, for the most massive case, the remnant lifetime, but the imprints on the postmerger gravitational-wave spectrogram and the dynamical ejecta mass are too small to be detectable with current or future detectors. They conclude that this cooling channel is unlikely to yield new axion mass constraints and that their setup represents a best-case scenario for the axion effect.","tokens_in":13170,"tokens_out":10843,"duration_ms":102214,"significance":"If the null result is robust, it closes a specific and previously unexplored avenue for using BNS postmerger observations to constrain the QCD axion mass. The paper has several strengths: it uses established numerical relativity codes (SGRID and BAM), reports constraint-violation checks in Appendix B, and deliberately neglects neutrino cooling and axion reabsorption so that the simulated cooling effect is an upper bound. The numerical framework is clearly described and the code behavior is validated. The main weakness is that the ``best case'' interpretation is incomplete with respect to equation-of-state and mass-ratio variations, and the abstract overstates the generality of the conclusion. The work is a solid first step that is likely to be influential for future axion searches.","major_comments":[{"comment":"The ``best case'' argument does not account for the sensitivity of the result to the nuclear equation of state, and the abstract therefore overstates the generality of the null result. The simulations use only the ALF2 EOS (Sec. III, Table I), while the axion emissivities in Sec. II A scale as epsDot_D proportional to rho^(1/3) T^6 and epsDot_ND proportional to rho^2 T^3.5; a hotter remnant, which may be produced by another EOS or by a different mass ratio, could increase the cooling rate by orders of magnitude. The manuscript itself concedes in Sec. V that ``other equations of state ... might lead to larger, more noticeable effects,'' yet the abstract claims that ``a more thorough investigation is unlikely to change the conclusions.'' The evidence presented establishes the null result only for the ALF2 equal-mass irrotational configurations considered. The abstract and conclusion should be reworded to restrict the claim to those configurations, or a quantitative argument for why ALF2 represents a near-maximal temperature case must be added.","section":"Abstract and Sec. V"}],"minor_comments":[{"comment":"The last paragraph of the introduction refers to ``Ref. IV'' when describing the results; this should read ``Sec. IV.''","section":"Sec. I"},{"comment":"The main text states that the GW spectrograms are computed for an inclination of iota = pi/4, while the caption of Fig. 4 says iota = pi/2; please reconcile this inconsistency.","section":"Sec. IV B and Fig. 4 caption"},{"comment":"The text says ``the ejection ... of the two neutron''; this should be ``the two neutron stars.''","section":"Sec. IV C"},{"comment":"The conclusion that ``with an increasing mass of the axion, the lifetime of the remnant before collapsing to a BH decreases'' is contradicted by the non-monotonic behavior reported in Sec. IV A, where the ma = 1 eV case has a longer lifetime than the ma = 0.1 eV case; please qualify this statement as a general trend rather than a strict monotonic decrease.","section":"Sec. V"},{"comment":"The footnote mentioning that the most recent GW170817 analysis suggests radii slightly smaller than ALF2 predicts is relevant to the representativeness of the EOS choice and would be more visible if discussed in the main text.","section":"Sec. III and footnote [92]"},{"comment":"The sentence ``The central core has a low temperature for all studied systems'' is confusing because the temperature maps in Fig. 2 show high temperatures in much of the inner region; please clarify the meaning of ``central core'' in this context.","section":"Sec. IV A"}],"recommendation":"major_revision","confidential_remarks":"The numerical work is sound and the paper fits the journal's scope. The main issue is the wording of the abstract, which claims a generality that the single-EOS simulations cannot support. The reader's report suggested acceptance, but I believe a major revision is appropriate to soften the claims and explicitly state the restrictions of the study."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper does something genuine: it is the first numerical relativity study to add nucleon-nucleon-axion bremsstrahlung cooling to BNS merger simulations, and it does so carefully. The null result for ALF2 equal-mass binaries is well supported: they use established codes, check constraint violations, scan axion masses from 0 to 1 eV, and deliberately stack the deck in favor of the axion by neglecting reabsorption and neutrino cooling. Under those conditions, the GW postmerger frequency shifts by less than 50 Hz and the ejecta mass changes are within the usual simulation uncertainty, so the conclusion that this channel will not sharpen axion constraints for the configurations studied is solid. The temperature differences they do find are physically sensible, and the non-monotonic lifetime for the most massive case is discussed honestly rather than smoothed over. Credit is due for a clean, well-documented setup and a sensible interpretation of a null run. The soft spot is real but narrow. The abstract and Section V claim this is a \"best case\" scenario and therefore a more thorough investigation is unlikely to change the conclusions. That is only true for the simplifications they explicitly model: free-streaming axions and no neutrino cooling. It does not cover equation-of-state sensitivity. The emissivity scales as T^6 (degenerate) and T^3.5 (non-degenerate), so a hotter remnant from a different EOS or mass ratio could boost the cooling by orders of magnitude. The authors concede in Section V that other EOSs might lead to larger effects, which undercuts the unqualified wording. The central null result still holds for ALF2 and for the physical reason that axion cooling is a small perturbation on the dynamical timescale, but generalizing it to all plausible BNS remnants requires at least one hotter EOS test or a clear argument that ALF2 is not unusually cold. This is a proportionality issue, not a fatal one. The paper deserves a serious referee. The overreach can be fixed with softer language and a caveat in the abstract, and the simulation itself is a useful reference point for future BNS-axion work, especially for people who want a baseline before trying more exotic cooling channels. I would not cite it in my own papers, but I would send it to review and would not be surprised if it gets published with minor revisions.","headline":"A clean, honest null result on axion cooling in BNS mergers that slightly overgeneralizes from a single EOS.","tokens_in":13694,"tokens_out":1545,"would_cite":false,"duration_ms":15820,"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":"The paper argues that axion cooling of neutron-star merger remnants via nucleon-nucleon-axion bremsstrahlung, though it changes the remnant's temperature profile, imprints too weakly on gravitational-wave signals and ejecta masses to…","keywords":["axion","binary neutron star merger","nucleon-nucleon-axion bremsstrahlung","numerical relativity","postmerger gravitational waves","axion cooling","QCD axion","multimessenger astronomy"],"falsifier":"A detection of a postmerger gravitational-wave signal whose dominant frequency is shifted by more than about 50 Hz relative to axion-free simulations, in a direction consistent with the axion-cooling template, would overturn the paper's null conclusion for this channel.","tokens_in":12709,"feed_emoji":"⚛️","tokens_out":5822,"duration_ms":53669,"temperature":0.7,"pith_summary":"Binary neutron star mergers heat matter to temperatures where the hypothetical QCD axion, if it exists, should be produced copiously through nucleon-nucleon-axion bremsstrahlung, providing an extra cooling channel. The paper tests whether this cooling leaves a detectable signature in the postmerger remnant—in its collapse time, gravitational-wave emission, or ejected mass. Using numerical relativity simulations with a covariant phenomenological cooling term, it finds that while the remnant's temperature profile and shape do change, the resulting shifts in the gravitational-wave frequency and ejecta mass are far below what current or future planned detectors could measure. Because the simulations omit neutrino cooling and axion reabsorption, they represent a best-case scenario, so the paper concludes this particular channel cannot yield new axion-mass constraints from binary neutron star observations. If correct, supernova 1987A neutrino observations remain the strongest stellar-cooling bound on this axion coupling.","feed_headline":"Postmerger signals won't pin down the axion mass","feed_subtitle":"Numerical relativity shows this cooling channel's imprint on gravitational waves and ejecta stays below detectability.","key_machinery":"The central object is the nucleon-nucleon-axion bremsstrahlung emission rate from Brinkmann and Turner, applied as a volumetric cooling term Λ = min(ϵ̇_D, ϵ̇_ND), with degenerate and non-degenerate rates proportional to $ρ^{{1/3}}$ $T^{6}$ and $ρ^{2}$ $T^{{3.5}}$ respectively. This is coupled to the hydrodynamics through a covariant radiation four-force G^α = -u^α Λ added to the general-relativistic energy-momentum equations, the same prescription used for neutrino cooling in earlier work. The min() selection stitches together the two regimes, and the four-force converts the microphysical emission rate into a local cooling of the fluid; the axion mass enters only through the coupling g ∝ m_a, so larger masses mean more cooling until reabsorption would set in.","core_discovery":"The central claim is that nucleon-nucleon-axion bremsstrahlung, modeled with the Brinkmann-Turner emission rates and a free-streaming radiation four-force, does cool the postmerger remnant—reducing its temperature, making it more spherical, and for the most massive binary shortening the time to black-hole collapse—but every observable imprint is too small to detect. The shift in the main postmerger gravitational-wave frequency stays below 50 Hz, too little for existing or upcoming analysis techniques, and the differences in dynamical ejecta mass are comparable to or smaller than current modeling uncertainties. The authors deliberately choose the most optimistic setup: no neutrino cooling, no axion reabsorption, and an equation of state consistent with GW170817. Even under these favorable assumptions, the effect is insufficient to place new limits on the axion mass, and they argue a more complete treatment would only weaken the signal further.","pith_inferences":["Other axion production mechanisms, such as axion-photon conversion in the merger's magnetic field or orbital effects during the inspiral, may still offer multimessenger constraints; this paper's null result applies only to the bremsstrahlung cooling channel.","The non-monotonic remnant lifetime seen for m_a = 1 eV (longer than for 0.1 eV) hints that chaotic core oscillations can mask small cooling effects; population-level analyses of remnant lifetimes would need many events to average over such dynamics.","A longer-lived remnant, e.g., from a softer equation of state or a lower-mass binary, would amplify the integrated cooling effect; targeted simulations of such systems could test whether any configuration pushes the imprint above detectability.","The same covariant cooling scheme can be applied to other feebly interacting particles (e.g., ALPs with different couplings) to decide which particle models are testable with postmerger observations."],"forward_implications":["Postmerger gravitational-wave observations, even with third-generation detectors, will not constrain the QCD axion mass through this cooling channel.","The remnant's faster sphericalisation and earlier collapse under strong axion cooling are real dynamical effects but shift collapse times only by a few milliseconds, below observability.","Kilonova ejecta masses from a binary neutron star merger carry no usable axion signal: the differences are within current ejecta-mass uncertainties.","Since the simulations omit reabsorption and neutrinos, more realistic modeling would make the axion imprint even smaller, reinforcing the null result.","Supernova 1987A bounds on the axion-nucleon coupling remain the leading stellar-cooling constraints, with no independent check from merger observations via this channel."],"supporting_citations":[{"why":"Supplies the degenerate and non-degenerate nucleon-nucleon-axion bremsstrahlung emission rates that define the cooling term.","marker":"[19]"},{"why":"Provides the covariant radiation four-force formalism used to couple cooling to general-relativistic hydrodynamics.","marker":"[36]"},{"why":"Companion paper establishing the covariant cooling scheme for neutron star merger simulations.","marker":"[57]"},{"why":"Sets the supernova 1987A cooling bounds that motivate the axion mass range tested here.","marker":"[32]"},{"why":"Identifies the reabsorption threshold that defines the paper's best-case scenario.","marker":"[33]"},{"why":"Provides the hadronic ALF2 equation of state used in all simulations.","marker":"[58]"},{"why":"Establishes the postmerger gravitational-wave frequency measurement capability against which the predicted shift is judged undetectable.","marker":"[82]"}],"fun_headline_variants":["Axion cooling too weak to show in merger signals","Merger remnants can't pin down axion mass","Axion imprint on postmerger waves undetectable","Neutron star mergers won't constrain axion emission","Axion bremsstrahlung leaves no merger signature"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central conclusion depends on the assumption that the Brinkmann-Turner nucleon-nucleon-axion bremsstrahlung rates, applied with a free-streaming prescription and without neutrino cooling, represent an upper bound on how much axions can cool a merger remnant.","fun_headline_variants_meta":{"raw":{"variants":["Axion cooling too weak to show in merger signals","Merger remnants can't pin down axion mass","Axion imprint on postmerger waves undetectable","Neutron star mergers won't constrain axion emission","Axion bremsstrahlung leaves no merger signature"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000291,"raw_usage":{"total_tokens":1726,"prompt_tokens":997,"completion_tokens":729,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":613,"completion_tokens_details":{"reasoning_tokens":650}},"tokens_in":613,"tokens_out":729,"duration_ms":7040,"temperature":1.0,"reasoning_tokens":650,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:21:55.026055+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A detection of a postmerger gravitational-wave signal whose dominant frequency is shifted by more than about 50 Hz relative to axion-free simulations, in a direction consistent with the axion-cooling template, would overturn the paper's null conclusion for this channel.","supporting_citations":[{"cited_title":"Importance of cooling in triggering the collapse of hypermassive neutron stars","cited_arxiv_id":"1208.5487","evidence_quote":"Provides the covariant radiation four-force formalism used to couple cooling to general-relativistic hydrodynamics."},{"cited_title":"Burrows, M","cited_arxiv_id":null,"evidence_quote":"Sets the supernova 1987A cooling bounds that motivate the axion mass range tested here."}],"review_version":1}