{"id":"efd6d8fd-b125-4b22-9677-1ae529cffc4e","arxiv_id":"1908.05815","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Radiative transfer simulations show that kilonova peak brightness and time in optical and near-infrared bands can distinguish prompt-collapse, magnetar-accelerated, and black hole-neutron star merger ejecta.","lead":"This paper simulates how kilonova light curves differ depending on what remains after two neutron stars or a neutron star and a black hole merge. It shows that the brightness and timing of the peak in optical and infrared bands could reveal which object powered the explosion.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Peak-magnitude classification rests on a handful of hand-picked ejecta models; parameter uncertainty is never propagated, so the claimed separability in Figs. 17/18 is not established.","rationale":"The reader's weakest assumption — that the ejecta profiles are fixed to a small set of NR-informed choices — is the same root issue I identify, and it is load-bearing. The central claim is not 'these light curves are diverse' (which is well supported) but that the peak location can be used to infer the central-engine type. That inference requires the scenario-to-ejecta mapping to be sufficiently narrow and the peak observables to be sufficiently separated, neither of which is quantified. The paper's own Section 7 explicitly defers the systematic parameter study needed to establish distinguishability, and Section 6.4 admits BH-NS events can resemble prompt-collapse or HMNS events for different ejecta masses. A Monte Carlo or grid test over the parameter ranges would settle whether the clusters in Figures 17/18 survive realistic scatter. The paper has real strengths: the updated r-process line list, the multi-component radiative transfer treatment, the explicit demonstration of viewing-angle effects, and the honest discussion of heating-rate and LTE limitations. These support the exploratory predictions but not the stronger classification claim. Since the reader already returned CONDITIONAL, my concern reinforces that condition rather than changing the verdict. I therefore recommend UNCHANGED.","tokens_in":53171,"tokens_out":5918,"duration_ms":63672,"concrete_test":"Run the same radiative-transfer code over a Monte Carlo grid of ejecta parameters drawn from the ranges cited in the paper, e.g., prompt collapse with Mpm = 0.001–0.01 Msun and Md = 0.0003–0.003 Msun, BH-NS with Md = 0.001–0.05 Msun, Mpm = 0–0.05 Msun, and Ye = 0.05–0.11, and HMNS/SMNS with corresponding mass, velocity, and Ye ranges. Compute the iJK peak magnitudes and peak times for each draw and build a confusion matrix between the four scenario regions in Figures 17/18. If a substantial fraction of prompt-collapse or BH-NS draws falls inside the HMNS cluster, the Section 6.4 inference claim should be weakened to a conditional prediction.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central inferential claim — that iJK peak brightness and peak time can reveal the merger/remnant type — requires each scenario to occupy a distinguishable region of ejecta-parameter space. But Table 1 provides only two or three hand-picked parameter sets per scenario, and Figures 17/18 draw cluster ellipses around those points without sampling the known ranges of ejecta mass, velocity, or Ye. The peak scalings in Eqs. (1)–(2) are steep: changing M by a factor of 3–10, as NR simulations do within a single scenario, shifts tpeak by sqrt(M) and peak magnitude by about 1 mag, comparable to or larger than several inter-cluster separations. The paper itself notes (Section 6.4 and Section 7) that BH-NS ejecta mass has a large variety and can mimic prompt-collapse or HMNS light curves, and that a quantitative distinguishability study is 'beyond the scope of this paper'. Thus the headline inference is not actually demonstrated; it is an illustrated conjecture whose false-positive rate is unquantified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents wavelength-dependent radiative transfer simulations of kilonova light curves for several post-merger scenarios: prompt collapse to a black hole, hypermassive neutron star with a long-lived remnant, supermassive neutron star with accelerated ejecta, and black hole-neutron star mergers. The ejecta density, velocity, and composition profiles are taken from representative numerical-relativity simulations, and a new r-process atomic line list is used. The authors identify three multi-component radiative transfer effects (polar diffusion, equatorial blocking, and heating of dynamical ejecta by post-merger ejecta), study the dependence of the light curves on ejecta mass, velocity, and electron fraction, and construct two models that reproduce the peak brightness of GW170817. The final section compares the peak magnitude and time of peak in the i, J, and K bands across the different scenarios and argues that these observables may allow inference of the central engine type.","tokens_in":53390,"tokens_out":3774,"duration_ms":38538,"significance":"If the inference from peak brightness and peak time to merger remnant type is robust, the paper offers a practical way to connect kilonova observations to the post-merger evolution of neutron star mergers, complementing gravitational-wave information. The work has several genuine strengths: it uses ejecta profiles motivated by numerical-relativity simulations, employs a complete atomic line list for r-process elements, explicitly models the non-spherical geometry of the dynamical ejecta, and includes a clear statement of the limitations of the radiative transfer treatment (early- and late-phase reliability, LTE assumption, heating-rate uncertainty). The scenario predictions are largely computed from first principles rather than fitted to the target light curves, and the only explicit fitting is the GW170817 model selection in Section 5. The central claim is falsifiable and clearly stated. However, as detailed below, the quantitative support for the central-engine inference is not yet established.","major_comments":[{"comment":"The headline claim that 'we may be able to infer the type of the central engine' from iJK peak magnitudes and times is not fully supported by the presented models. Table 1 provides only two or three hand-picked ejecta parameter sets per scenario, and the known ranges of ejecta mass, velocity, and Ye from numerical-relativity simulations are not sampled. Equations (1) and (2) imply that a factor of 3 to 10 variation in ejecta mass shifts the peak magnitude by roughly 1 mag and the peak time by a factor of sqrt(M), which is comparable to or larger than several inter-cluster separations in Figures 17 and 18. The authors themselves note in Section 6.4 that BH-NS ejecta with small mass could mimic prompt-collapse or HMNS light curves, and in Section 7 that quantitative distinguishability is 'beyond the scope of this paper.' This means the proposed classification is an illustrated conjecture whose false-positive rate is unquantified. I recommend either adding a systematic parameter scan with uncertainty propagation or explicitly reframing the conclusion as a qualitative illustration rather than a demonstrated inference.","section":"Section 6.4, Figures 17 and 18, and Section 7"},{"comment":"The peak magnitude is defined as the brightest magnitude for t ≥ 1 day, but for the SMNS (accelerated ejecta) models the light curves in Figure 15 appear to peak before 1 day, with the riz bands brighter than GW170817 at t ≲ 1 day and declining rapidly afterward. Since the early phase t < 1 day is excluded as unreliable, the SMNS points in Figures 17 and 18 may not represent the true peak magnitudes or peak times, which could bias the apparent separation between the SMNS cluster and other scenarios. The authors should quantify how the t ≥ 1 day restriction affects the SMNS peak values, or provide an approximate treatment of the early phase (for example, a blackbody estimate) to test whether the qualitative separation persists.","section":"Section 3.3 and Section 6.4 (peak definition)"}],"minor_comments":[{"comment":"There is a typo: 'BN-NS mergers' should read 'BH-NS mergers' in the paragraph describing electron fraction and abundances for black hole-neutron star ejecta.","section":"Section 3.2"},{"comment":"The word 'fssion' in the figure captions should be 'fission'.","section":"Captions of Figures 22 and 23"},{"comment":"The figure contains multiple repeated panels that are difficult to distinguish; the authors should reorganize it so that each panel is unique and clearly labeled.","section":"Figure 21"},{"comment":"The table is dense and the average velocities are given in parentheses without a clear explanation of the notation in the caption; a footnote defining v_ave and its relation to the kinetic energy would improve readability.","section":"Table 1"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid and useful parameter study, and the radiative transfer methodology is sound. The main issue is that the central inference about distinguishing merger scenarios from peak light curves is presented more strongly than the sparse sampling of the ejecta parameter space warrants. If the authors either add a systematic parameter exploration or clearly limit the claim to a qualitative illustration, the paper would be suitable for publication. The scope of the journal is appropriate for this type of multi-messenger follow-up study."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper deserves a serious referee, but the headline claim should not be oversold. What is new: the authors take an established radiative transfer scheme and run a systematic multi-band suite for prompt collapse, magnetar-accelerated (SMNS) ejecta, and BH-NS mergers, using an updated r-process line list and axisymmetric two-component ejecta profiles. The iJK peak-magnitude versus peak-time diagrams (Figs. 17/18) are genuinely useful organizing tools for observers planning follow-up, and the paper is careful about its own limitations: LTE breakdown, unreliable early and late phases, heating-rate and fission-fragment uncertainty, and the possibility that BH-NS ejecta mass varies enough to mimic other scenarios.\n\nThe soft spots are real but not fatal. The central inference — that central-engine type can be read off from peak brightness and peak time — is not actually demonstrated. Table 1 gives one to three hand-picked parameter sets per scenario, and the figures draw cluster ellipses around those points without sampling the known spread in ejecta mass, velocity, or Ye from numerical-relativity simulations. Since the analytic scalings (Eqs. 1–2) shift peak brightness by roughly a magnitude when M changes by a factor of 3–10, and NR simulations show such ranges within single scenarios, the apparent separation in the iJK peak plots could shrink or shift. The authors themselves write that a quantitative distinguishability study is beyond the scope of the paper. That is the right disclaimer, but it also means this should not be cited as proof of separability, only as an illustration.\n\nThere is no public code or atomic-data artifact, and no error bars on any predicted magnitude. The GW170817-matching models in Section 5 are chosen to reproduce the observed peak, so that section is model fitting rather than prediction, and the paper says so. The fission-heating sensitivity check for representative models is good practice, and the finding that BH-NS light curves depend on it more strongly is itself useful.\n\nBottom line: a workmanlike, useful simulation paper. The multi-component radiative-transfer effects (polar enhancement, blocking, heating) and the scenario diversity are worth knowing, and the limitations are stated honestly. What it does not do is establish a robust classification diagnostic. Send it to referee; the right referee will push for uncertainty quantification or at least a clearer statement that the clusters are illustrative, not predictive.\n\nMy recommendation: accept after minor-to-moderate revision.","headline":"Useful, honest survey of kilonova model diversity; the peak-magnitude classification is illustrative, not proven, but the paper deserves peer review.","tokens_in":53903,"tokens_out":1831,"would_cite":true,"duration_ms":19503,"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":"Kilonova light-curve peaks carry a fingerprint of the merger's central engine, letting observers distinguish prompt collapse, magnetar-accelerated ejecta, and black hole-neutron star mergers.","keywords":["kilonova","neutron star merger","radiative transfer","r-process nucleosynthesis","black hole-neutron star merger","magnetar","light curves","GW170817"],"falsifier":"Measure the iJK peak magnitudes and peak times of a gravitational-wave counterpart of known distance and inclination. A prompt-collapse candidate that is optically within about 1 mag of GW170817, or a magnetar-accelerated candidate that stays brighter than GW170817 after five days, would contradict the predicted separation; conversely, the detection of a population of faint-optical but bright-infrared kilonovae would confirm it.","tokens_in":52970,"feed_emoji":"🔭","tokens_out":4050,"duration_ms":38040,"temperature":0.7,"pith_summary":"This paper tries to establish that the peak brightness and time of peak of a kilonova's multi-band light curve encode which kind of merger produced it. Using ejecta profiles from numerical-relativity simulations, the authors compute light curves for neutron-star mergers that promptly collapse to a black hole, mergers with a long-lived remnant that accelerates its ejecta, and black hole-neutron star mergers. They find that these cases separate in the i, J, and K bands: prompt collapse is optically faint but infrared-bright, accelerated ejecta is briefly bright then declines fast, and BH-NS mergers are infrared-bright by 1-2 magnitudes. If correct, a handful of photometric points around peak could classify the central engine of a gravitational-wave counterpart.","feed_headline":"Kilonova peaks reveal the merger's hidden fate","feed_subtitle":"Brightness and timing in infrared bands separate prompt collapse, magnetar-boosted ejecta, and black hole-neutron star mergers.","key_machinery":"The load-bearing tool is a wavelength-dependent Monte Carlo radiative transfer code that treats two ejecta components together—a spherical post-merger ejecta and a non-spherical dynamical ejecta—so that photons diffusing preferentially toward the pole, being blocked near the equator, and reprocessed to heat the dynamical ejecta are captured self-consistently. It uses a new line list from atomic structure calculations for all r-process elements ($Z=26$-$92$) and r-process heating rates with thermalization. The diagnostic output is the peak magnitude and time of peak in the i, J, and K bands (Figures 17 and 18), where each merger scenario clusters.","core_discovery":"The central claim is that differences in ejecta properties are imprinted in the peak brightness and time of peak, so that observing the peak in multiple bands allows one to infer the type of central engine. Concretely: (i) optical emission from prompt-collapse mergers is fainter by $\\gtrsim 1$-$2\\,{\\rm mag}$ than GW170817 while the infrared stays as bright if post-merger ejecta is about $0.01\\,M_\\odot$; (ii) magnetar-accelerated ejecta outshines GW170817 by 1-2 mag for the first few days but fades below it within days; (iii) black hole-neutron star mergers with $\\gtrsim 0.02\\,M_\\odot$ of ejecta can be optically as bright as GW170817 and infrared-brighter by 1-2 mag. The paper argues that the resulting clustering in peak magnitude-peak time space is a practical diagnostic for the merger evolution.","pith_inferences":["The same peak-diagnostic logic could be applied to archival short gamma-ray burst afterglows with kilonova candidates, testing whether those events cluster by engine type.","If gravitational-wave measurements provide the inclination, the predicted viewing-angle dependence of the optical suppression could sharpen or falsify the classification before a large sample exists.","Extending the calculation to non-LTE and better late-time heating would directly test whether the fast-declining light curves are real or an artifact of the LTE assumption.","The 20-30 percent uncertainty from the capped electron fraction for BH-NS dynamical ejecta could be resolved by rerunning with $Y_e \\approx 0.05$ tables when available."],"forward_implications":["Prompt-collapse kilonovae may be missed in optical surveys but remain detectable in near-infrared for about a week.","A kilonova that peaks more than 1-2 mag brighter than GW170817 and fades within days points to magnetar-accelerated ejecta.","A BH-NS merger can masquerade as a normal optical kilonova but should stand out as infrared-bright.","Multi-band peak observations suffice to break degeneracies that single near-infrared bands have with respect to ejecta mass.","Ejecta mass estimates from optical brightness alone can be wrong by a factor of about two if polar diffusion enhancement is ignored."],"supporting_citations":[{"why":"Supplies the GW170817 optical and infrared light-curve data used as the reference for all brightness comparisons.","marker":"Villar et al. (2017)"},{"why":"Provides numerical-relativity results for the radial density profile of dynamical ejecta from NS mergers.","marker":"Kiuchi et al. (2017)"},{"why":"Together with Kiuchi et al., anchors the broken-power-law radial profile for the dynamical ejecta.","marker":"Hotokezaka et al. (2018)"},{"why":"Supplies the r-process nucleosynthesis element abundances and heating rates used in the radiative transfer.","marker":"Wanajo et al. (2014)"},{"why":"Provides the new atomic line list for all r-process elements that determines the opacities.","marker":"Tanaka et al. (2019)"},{"why":"Supplies the BH-NS ejecta density and electron fraction profiles used for the BH-NS merger models.","marker":"Foucart et al. (2017)"},{"why":"Supplies additional BH-NS ejecta profiles and the typical dynamical ejecta velocity for tidal disruption cases.","marker":"Kyutoku et al. (2018)"},{"why":"Establishes the prior multi-component radiative transfer treatment of non-spherical ejecta that this work extends.","marker":"Kawaguchi et al. (2018)"},{"why":"Identifies the blocking of optical emission by dynamical ejecta, a key effect quantified here.","marker":"Kasen et al. (2015)"}],"fun_headline_variants":["Kilonova peaks fingerprint the merger's engine","Peak brightness and timing decode kilonova origins","Kilonova peak times reveal the central engine type","Light-curve peaks expose the merger's hidden machinery","Infrared and optical peaks unmask kilonova engines"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The predictions rest on the ejecta mass, velocity, angular shape, and electron fraction being fixed to the representative profiles taken from a small set of numerical relativity simulations; if real ejecta differ, the claimed 1-2 magnitude separations and peak-time clusters could shrink or shift.","fun_headline_variants_meta":{"raw":{"variants":["Kilonova peaks fingerprint the merger's engine","Peak brightness and timing decode kilonova origins","Kilonova peak times reveal the central engine type","Light-curve peaks expose the merger's hidden machinery","Infrared and optical peaks unmask kilonova engines"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000224,"raw_usage":{"total_tokens":1519,"prompt_tokens":1065,"completion_tokens":454,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":681,"completion_tokens_details":{"reasoning_tokens":378}},"tokens_in":681,"tokens_out":454,"duration_ms":4695,"temperature":1.0,"reasoning_tokens":378,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:04:03.076068+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the iJK peak magnitudes and peak times of a gravitational-wave counterpart of known distance and inclination. A prompt-collapse candidate that is optically within about 1 mag of GW170817, or a magnetar-accelerated candidate that stays brighter than GW170817 after five days, would contradict the predicted separation; conversely, the detection of a population of faint-optical but bright-infrared kilonovae would confirm it.","supporting_citations":[{"cited_title":"A., et al","cited_arxiv_id":null,"evidence_quote":"Supplies the GW170817 optical and infrared light-curve data used as the reference for all brightness comparisons."},{"cited_title":"2017, Phys","cited_arxiv_id":null,"evidence_quote":"Provides numerical-relativity results for the radial density profile of dynamical ejecta from NS mergers."},{"cited_title":"2018, Astrophys","cited_arxiv_id":null,"evidence_quote":"Together with Kiuchi et al., anchors the broken-power-law radial profile for the dynamical ejecta."},{"cited_title":"2014, Astrophys","cited_arxiv_id":null,"evidence_quote":"Supplies the r-process nucleosynthesis element abundances and heating rates used in the radiative transfer."},{"cited_title":"2017, Class","cited_arxiv_id":null,"evidence_quote":"Supplies the BH-NS ejecta density and electron fraction profiles used for the BH-NS merger models."},{"cited_title":"2018, Phys","cited_arxiv_id":null,"evidence_quote":"Supplies additional BH-NS ejecta profiles and the typical dynamical ejecta velocity for tidal disruption cases."},{"cited_title":"2018, Astrophys","cited_arxiv_id":null,"evidence_quote":"Establishes the prior multi-component radiative transfer treatment of non-spherical ejecta that this work extends."},{"cited_title":"2015, Mon","cited_arxiv_id":null,"evidence_quote":"Identifies the blocking of optical emission by dynamical ejecta, a key effect quantified here."}],"review_version":1}