{"id":"433ba699-2f5d-4201-8db9-7b545b9c3025","arxiv_id":"2411.17047","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"M33 X-7, NGC 300 X-1 and IC 10 X-1 show spectral evidence for a hot slim accretion disk plus a cool standard disk, with black hole masses around 9-30 solar masses.","lead":"This paper analyzes all X-ray observations of three massive black hole binaries outside our Galaxy and finds their hot accretion disks may be 'slim', puffed-up disks rather than the standard thin disks seen in most such systems. The work also gives new mass ranges for the three black holes and measures the sizes of the stars that eclipse them.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The slim-disc inference is not statistically secured: Model-1 and Model-2 are often statistically equivalent, and the diskpbb index p pegs at its lower boundary, so the hot slim disc and derived masses are not uniquely supported.","rationale":"The reader's weakest_assumption identifies the correct load-bearing concern: the data do not uniquely force the replacement of a Comptonized power-law with a hot slim disc. The paper's own language in Section 3.1 concedes that Model-1 is statistically satisfactory and that the difference 'may not be significant,' which undercuts the central claim. I add a concrete technical detail: in many Model-2 fits the parameter p—the very parameter that distinguishes a slim disc (p ≈ 0.5) from a standard disc (p = 0.75)—pegs at its lower boundary, so the claimed radial temperature profile is not measured but imposed. All subsequent results, including the SPL classification, the inferred vertically extended inner region, and the continuum-fitting masses from slimbh, inherit this model degeneracy. That said, the paper is transparent about its limitations, uses public data and standard XSPEC tools, and provides a sanity check for M33 X-7: the slimbh mass (14.86 M☉) is consistent with the dynamical mass (15.65 M☉). This is real supporting evidence, so the appropriate outcome is not rejection but a conditional acceptance requiring a proper model comparison and a treatment of the p boundary. The reader's verdict of CONDITIONAL is therefore appropriate and should be kept.","tokens_in":32118,"tokens_out":5181,"duration_ms":49183,"concrete_test":"Refit the highest-S/N non-eclipse spectra (e.g., M33 X-7 MN18+MX19, NGC 300 X-1 NX7/NX8, IC 10 X-1 IX2) with Model-1 and Model-2 using identical energy ranges, binning, and background files, and compute a non-nested model comparison such as AIC/BIC or a simulation-calibrated Δχ² distribution. Report the per-epoch ΔBIC; if the preference for Model-2 is ΔBIC < 5 in most epochs, or if the profile likelihood in p includes p = 0.75, the slim-disc claim is not statistically established and the paper should remain conditional.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3.1 explicitly states that Model-1 (Tbabs*Tbabs(diskbb+apec+powerlaw)) 'provides statistically satisfactory fits' and that the difference relative to Model-2 (Tbabs*Tbabs(diskbb+apec+diskpbb)) 'may not be significant'; Model-2 is chosen primarily because Model-1 is deemed physically inadequate. This is a physical preference, not a statistical detection. No non-nested model comparison (AIC/BIC, Bayesian evidence, or simulation-calibrated Δχ²) is presented for Model-1 versus Model-2; the F-test mentioned in Section 3.1 is used only to justify adding a soft diskbb component to diskpbb, not to justify replacing powerlaw with diskpbb. In many epochs the χ² values are nearly identical (e.g., M33 X-7 MX1: Model-1 χ²/dof = 235.46/221 versus Model-2 236.15/222). Furthermore, Table 3 shows that the diskpbb index p is pegged at the hard lower limit 0.50 in a large fraction of Model-2 fits, so the claimed p = 0.5–0.66 profile is not actually constrained below the boundary in those epochs. All downstream conclusions—the hot slim disc geometry, the p < 0.75 advection signature, the SPL-state classification, and the slimbh mass ranges in Section 3.2—presuppose that Model-2 is the correct continuum. If a Comptonized power-law describes the same data equally well, the central claim does not follow.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a comprehensive X-ray spectral and timing study of three eclipsing extragalactic black hole X-ray binaries: M33 X-7, NGC 300 X-1, and IC 10 X-1, using all available XMM-Newton and NuSTAR observations. The authors first fit the non-eclipse spectra with a diskbb+powerlaw continuum (Model-1), finding steep photon indices and sub-Eddington luminosities that they associate with the steep power-law (SPL) state. They then replace the powerlaw with a second thermal component, diskpbb, yielding Model-2 (diskbb+apec+diskpbb), which they interpret as a hot slim disc with temperature profile T(r) ∝ r^{-p}, p = 0.5–0.66, alongside a cooler standard disc at 0.1–0.2 keV. Using the relativistic slim-disc model slimbh, they estimate black hole masses of 8.9–14.9 M_sun, 8.7–28 M_sun, and 10.2–30 M_sun for M33 X-7, NGC 300 X-1, and IC 10 X-1. Eclipse light-curve modelling yields companion-star radii of ~10, ~10, and ~18 R_sun, and eclipse spectral modelling is used to infer that the soft disc component is completely obscured while a hard component remains, supporting a vertically extended inner accretion region.","tokens_in":32513,"tokens_out":4938,"duration_ms":44511,"significance":"If the central spectral interpretation is correct, the paper would provide evidence for slim-disc accretion in sub-Eddington wind-fed extragalactic binaries, potentially linking their SPL-like spectra to the ultraluminous states seen in ULXs, and would add new mass estimates for three high-mass black hole X-ray binaries. The paper makes good use of a large, multi-epoch dataset, includes simultaneous broad-band fits where available, and attempts a consistency check of the slimbh continuum-fitting method on M33 X-7, whose dynamically measured mass is reproduced. The eclipse analysis is also a useful addition. However, the load-bearing inference—that the non-eclipse spectra require a hot slim disc rather than a Comptonized powerlaw—rests on a model preference that the authors themselves describe as statistically marginal, and several of the quoted slim-disc parameters are at hard parameter boundaries. The mass and geometry conclusions are conditional on this model choice and on additional assumptions about spin, inclination, and the mass of IC 10 X-1.","major_comments":[{"comment":"The preference for Model-2 over Model-1 is not established statistically, yet the entire slim-disc interpretation depends on it. The text states that Model-1 'provides statistically satisfactory fits' and that 'the difference may not be significant,' and the F-test reported is used to justify adding a soft diskbb to diskpbb, not to justify replacing the powerlaw by diskpbb. In several epochs the chi-square values are nearly identical, e.g., M33 X-7 MX1: Model-1 chi2/dof = 235.46/221 versus Model-2 = 236.15/222, and MX6: 234.80/219 versus 234.90/220. A formal non-nested comparison (AIC/BIC, Bayesian evidence, or simulation-calibrated Delta chi2) is needed to support the claim that the hot slim disc is actually required; without it, the paper should explicitly present the slim-disc scenario as conditional on a physically motivated but statistically indistinguishable alternative.","section":"Section 3.1"},{"comment":"The quoted radial temperature index p = 0.5–0.66 is not actually constrained in a large fraction of the fits because the parameter is pegged at its hard lower limit. Many entries in Table 3 carry the '‡' marker indicating p = 0.50 at the hard limit, including most M33 X-7 epochs and several NGC 300 X-1 epochs. In these cases the data only provide an upper limit, not a measurement of p < 0.75, so the discussion's inference that low p indicates advection is supported only by the few epochs where p is free and significantly above 0.5 (notably IC 10 X-1). The paper should state how many epochs have p pegged and weaken the advection claim accordingly.","section":"Table 3 and Section 3.1"},{"comment":"The abstract and discussion claim that the sources consist of a hot slim disc 'along with a cooler standard disc,' but Table 3 shows that most M33 X-7 epochs are fit by diskpbb alone, with no diskbb component listed (e.g., MX1–MX3, MX5–MX13, MX15–MX17). Only a few epochs (MX4, MX14, MN18+MX19) require the soft diskbb. The two-component geometry presented in Figure 8 therefore applies only to a subset of the data, and the discussion should be qualified to reflect that the cool standard disc is not generally detected.","section":"Section 3.1 and Table 3"},{"comment":"The eclipse spectral analysis that is used to infer the vertically extended hard emitter is based on Model-1 (powerlaw), not the preferred Model-2. The text reports that during eclipse 'neither of the disk models (diskbb or diskpbb) adequately fit the spectrum' and that the eclipse spectrum could only be fitted with an absorbed powerlaw. Using the very component rejected in the non-eclipse analysis to draw conclusions about the accretion geometry is internally inconsistent, and it weakens the claim that the eclipse results independently confirm the Model-2 picture. The paper should either reconcile this inconsistency or clearly present the eclipse interpretation as tentative.","section":"Section 4.2"},{"comment":"Several derived quantities depend on assumptions that are not independently constrained: the Eddington-scaled luminosity of IC 10 X-1 assumes a mass of 17 M_sun copied from NGC 300 X-1, and the companion-star radii in Table 4 use the average slimbh mass derived in Section 3.2. Additionally, the slimbh mass ranges for NGC 300 X-1 and IC 10 X-1 assume high spin (a > 0.8) and fixed inclination angles, with Figure 4 showing that the mass is strongly degenerate with spin. These internal assumptions should be stated prominently wherever the mass ranges and luminosities are quoted, and the inferred masses should be presented as conditional on the adopted model and prior constraints.","section":"Section 3.2 and Table 3"}],"minor_comments":[{"comment":"In the paragraph describing Table 3, the text says the p parameter is '0.5−5.56 in NGC 300 X-1,' which appears to be a typo for 0.5–0.56; please correct it.","section":"Section 3.1"},{"comment":"Several entries in Table 1 are incomplete or run together, such as the NX7 exposure columns ('20 10 9-') and the IX1 row; the table should be reformatted so eclipse timings are unambiguous.","section":"Table 1"},{"comment":"The caption of Figure 4 should state more explicitly that the spin range for M33 X-7 is restricted to 0.79–0.89 while NGC 300 X-1 and IC 10 X-1 use 0.8–0.98; the current text is slightly ambiguous.","section":"Figure 4"},{"comment":"The statement that the diskbb normalization gives a disc radius of 1000–6500 km (<130 R_g) should specify the black hole mass used to convert kilometres to gravitational radii, since the conversion is mass-dependent.","section":"Section 5.3"},{"comment":"Equation (1) defines R_in using correction factors kappa and xi, but the text does not discuss how uncertainties in these factors propagate into the mass ranges; a brief comment would be helpful.","section":"Section 3.2"}],"recommendation":"major_revision","confidential_remarks":"The paper has interesting and potentially valuable data products, but the central claim is not yet statistically secured: the authors themselves acknowledge that Model-1 fits are statistically satisfactory and that the difference from Model-2 may not be significant, and the diskpbb index is pegged at its boundary in many epochs. I would encourage the editor to request a strengthened model-comparison section or a reframing of the results as conditional rather than definitive. The eclipse analysis using Model-1 while the spectral analysis prefers Model-2 is an internal inconsistency that needs to be resolved. The manuscript is otherwise within scope for a journal like MNRAS, and the dataset is substantial."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper delivers a genuinely useful thing: the first systematic multi-epoch XMM-Newton plus NuSTAR spectral census of M33 X-7, NGC 300 X-1, and IC 10 X-1, together with first slimbh continuum-fitting mass estimates for these sources. The eclipse timing and spectral analysis is careful and adds something new: the complete obscuration of the soft component while hard flux persists is a real observational result, and the derived companion radii are consistent with WR-star expectations. The M33 X-7 mass from slimbh (about 15 solar masses) agrees with the dynamical value, which is a good sanity check. I give credit for using public data in a standard way and for disclosing many limitations in the text.\n\nThe soft spot is exactly what the stress-test note identifies. The central claim is that a hot slim disc (diskpbb) plus a cool standard disc (diskbb) describes these spectra better than a standard disc plus power law. But the authors themselves state in Section 3.1 that Model-1 provides statistically satisfactory fits and that the difference \"may not be significant.\" The choice of Model-2 is largely justified by physical preference, not by statistical necessity. No non-nested model comparison (AIC, BIC, or simulation-calibrated delta-chi-squared) is presented for Model-1 versus Model-2. In several epochs the chi-squared values are nearly identical. Also, the diskpbb index p is pegged at the hard lower limit 0.50 in a large fraction of the fits, so the quoted p = 0.5–0.66 range overstates what is actually constrained. The downstream inference of vertical extension, advection, and SPL classification all rest on this model choice.\n\nA related but more minor concern is the chain of assumptions in the mass estimates: spin, inclination, and distance are fixed or scanned, and for IC 10 X-1 the Eddington-scaled luminosity assumes its mass equals that of NGC 300 X-1. These are disclosed, and the mass range is stated as conditional, but it means the numbers should be read as illustrative rather than firm.\n\nNone of this sinks the paper. The observational work is solid and the limitations are transparent. But the abstract and conclusions are more assertive than the statistics justify. With a proper model comparison and a clearer separation of fitted parameters from assumed priors, this could be a solid contribution.\n\nRecommendation: send to peer review. A good referee can push for the missing model comparison and a re-presentation of the p values that respects the boundary. The data analysis and eclipse results are worth publishing regardless of whether the slim-disc interpretation survives.","headline":"A thorough, useful X-ray census of three extragalactic BH-XRBs whose headline slim-disc claim is physically motivated but not statistically forced; the authors are candid about that, so the paper deserves a serious referee with a request for a proper model comparison.","tokens_in":33063,"tokens_out":1551,"would_cite":true,"duration_ms":16411,"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":"Three extragalactic black hole binaries are better described by hot slim discs plus cool outer discs than by a disc plus corona.","keywords":["X-ray binaries","accretion discs","slim disc model","steep power-law state","black hole masses","M33 X-7","NGC 300 X-1","IC 10 X-1"],"falsifier":"Take the $3$--$79$ keV spectrum of M33 X-7 with NuSTAR to high signal to noise: a true slim disc should show a Wien-like rollover at a few times the inner temperature and no unresolved power-law tail, whereas a corona would produce a distinct hard tail and a Compton hump above $\\sim 10$ keV; the paper's NuSTAR data only reach $10$ keV, so the test remains open. Alternatively, a measured spin for NGC 300 X-1 or IC 10 X-1 that falls far below $0.8$ would break the mass estimates.","tokens_in":31846,"feed_emoji":"🕳️","tokens_out":13026,"duration_ms":100324,"temperature":0.7,"pith_summary":"This paper argues that the persistent X-ray emission of three massive, eclipsing extragalactic black hole X-ray binaries, M33 X-7, NGC 300 X-1, and IC 10 X-1, is not the usual standard thin disc plus Comptonizing corona. Instead, the non-eclipse spectra are best fitted with a hot slim disc at $1$--$2$ keV (temperature profile $T(r)\\propto r^{-p}$ with $p=0.5$--$0.66$) together with a cooler standard disc at $0.1$--$0.2$ keV. If this interpretation is right, the sources are in a sub-Eddington steep power-law state with a vertically extended inner accretion region, and continuum fitting with a relativistic slim-disc model yields black hole masses of roughly $9$--$15$, $9$--$28$, and $10$--$30$ solar masses. The eclipse analysis supports the geometry: the soft disc component is completely hidden during eclipse while a hard component from an extended region persists, implying a puffed-up inner flow.","feed_headline":"X-ray spectra reveal hot slim discs around three massive black holes","feed_subtitle":"New modelling places these eclipsing binaries in a steep power-law state with masses of 9–30 solar masses.","key_machinery":"The central object is the slim-disc model, represented by the phenomenological diskpbb component and the relativistic slimbh model. diskpbb is a multicolour disc model with a free radial temperature index $p$ in $T(r)\\propto r^{-p}$; the fits give $p=0.5$--$0.66$ and $kT_{in}=1$--$2$ keV for the hot inner component, whereas a standard thin disc has $p=0.75$. This index is the diagnostic of vertical extension: the paper reads $p<0.75$ as an advection-dominated, optically thick flow with scale height $H/R\\sim 1$. slimbh then turns the same continuum into a black hole mass via the relation between inner radius and mass, using assumed spin and inclination. The corroborating mechanism is the eclipse spectral analysis, where the in-eclipse spectra of IC 10 X-1 and NGC 300 X-1 contain only a power-law component with high covering fraction ($\\sim 0.9$), indicating that the soft disc is entirely obscured while the extended hard-emitting region is not.","core_discovery":"The central claim is that the $0.3$--$8$ keV spectra of M33 X-7, NGC 300 X-1, and IC 10 X-1, assembled from all XMM-Newton and NuSTAR observations, are better reproduced by two thermal components than by the conventional disc-blackbody-plus-power-law model. The hot component, a diskpbb slim disc with inner temperature $kT_{in}=1$--$2$ keV and radial index $p=0.5$--$0.66$, replaces the power law; the lower value of $p$ compared with the standard $0.75$ is interpreted as evidence of advection and a vertically extended inner flow. A cool diskbb component at $0.1$--$0.2$ keV accounts for the soft excess, which the paper attributes to a faint outer thin disc rather than a wind. The authors argue that the power-law component in the conventional model is an extrapolation artifact from the soft band and that Comptonization models leave unconstrained parameters. They then apply the relativistic slim-disc model slimbh in place of diskpbb, obtaining black hole masses of $8.9$--$14.9$, $8.7$--$28$, and $10.2$--$30$ $M_\\odot$ for the three sources, consistent with dynamical estimates where available.","pith_inferences":["If the slim-disc interpretation is correct, the same two-thermal-component description could apply to other wind-fed HMXBs and to Galactic sources in the steep power-law state; this is testable with broadband spectra that resolve the soft excess.","The paper's mass ranges for NGC 300 X-1 and IC 10 X-1 assume spins above $0.8$; a measurement of spin via reflection modelling or with future X-ray polarimetry would either tighten or invalidate the claimed masses.","A dedicated search for X-ray polarization of the $1$--$2$ keV component could distinguish a true slim disc, which should show a specific polarization-angle swing, from a hot corona, which would not, providing a clean test of the model choice."],"forward_implications":["The three sources are interpreted as being caught in the steep power-law state, a sub-Eddington state previously identified in Galactic transients but not firmly in these extragalactic wind-fed systems.","The derived black hole masses place M33 X-7, NGC 300 X-1, and IC 10 X-1 among the most massive stellar-mass black holes known, making them plausible progenitors of merging binary black holes.","The two-zone accretion geometry — hot slim disc inside, cool standard disc outside — would describe wind-fed systems at luminosities $0.1$--$0.7\\,L_{\\rm Edd}$, extending the regime where slim-disc behaviour appears.","Because the soft component is fully eclipsed while hard emission persists, any complete model of these systems must put the hard emitter on a larger scale height than the compact soft emitter, a constraint independent of the spectral decomposition."],"supporting_citations":[{"why":"Introduces the slim-disc accretion model that motivates replacing the power-law with a hot thermal component.","marker":"Abramowicz et al. 1988"},{"why":"Defines the standard thin-disc spectrum with $T\\propto r^{-0.75}$, the contrast used to argue that $p<0.75$ implies advection.","marker":"Shakura & Sunyaev 1973"},{"why":"Provides the precedent of describing sub-Eddington black hole spectra with diskpbb, justifying the same move here.","marker":"Kubota & Done 2004"},{"why":"Supplies the relativistic slim-disc model slimbh used for the continuum-fitting mass estimates.","marker":"Sadowski 2011"},{"why":"Gives the geometric factor $\\xi=0.353$ used to convert the disc normalization into an inner radius and hence a mass.","marker":"Vierdayanti et al. 2008"},{"why":"Provides the dynamical mass of M33 X-7 used to validate the slimbh continuum-fitting method.","marker":"Orosz et al. 2007"},{"why":"Provides the eclipse-duration formula used to derive the companion star radii from the measured eclipse periods.","marker":"Laycock et al. 2015"},{"why":"Supplies the inclination range and binary parameters of NGC 300 X-1 used as inputs in the slimbh fits.","marker":"Binder et al. 2021"}],"fun_headline_variants":["Hot slim discs found around three massive black holes","Spectra show hot slim discs in three extragalactic black holes","New model reveals hot slim discs and masses of three black holes","Three black holes host hot slim discs according to X-ray data","X-ray data expose hot slim discs around three massive black holes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the hard continuum is thermal emission from a hot slim disc rather than Comptonized radiation from a corona; the paper itself concedes that the disc-plus-power-law model yields statistically acceptable fits, so the switch to the slim-disc model rests on physical interpretation and not on fit quality alone.","fun_headline_variants_meta":{"raw":{"variants":["Hot slim discs found around three massive black holes","Spectra show hot slim discs in three extragalactic black holes","New model reveals hot slim discs and masses of three black holes","Three black holes host hot slim discs according to X-ray data","X-ray data expose hot slim discs around three massive black holes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000943,"raw_usage":{"total_tokens":4191,"prompt_tokens":1273,"completion_tokens":2918,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":889,"completion_tokens_details":{"reasoning_tokens":2833}},"tokens_in":889,"tokens_out":2918,"duration_ms":20071,"temperature":1.0,"reasoning_tokens":2833,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:36:03.190584+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the $3$--$79$ keV spectrum of M33 X-7 with NuSTAR to high signal to noise: a true slim disc should show a Wien-like rollover at a few times the inner temperature and no unresolved power-law tail, whereas a corona would produce a distinct hard tail and a Compton hump above $\\sim 10$ keV; the paper's NuSTAR data only reach $10$ keV, so the test remains open. Alternatively, a measured spin for NGC 300 X-1 or IC 10 X-1 that falls far below $0.8$ would break the mass estimates.","supporting_citations":[{"cited_title":"A., et al., 2007, @doi [ ] 10.1038/nature06218 , https://ui.adsabs.harvard.edu/abs/2007Natur.449..872O 449, 872","cited_arxiv_id":null,"evidence_quote":"Provides the dynamical mass of M33 X-7 used to validate the slimbh continuum-fitting method."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the eclipse-duration formula used to derive the companion star radii from the measured eclipse periods."}],"review_version":1}