REVIEW 5 major objections 12 minor 32 references
The X-ray Variability of the Ultraluminous X-ray Sources in the NGC 4631 galaxy
T0 review · 5 major / 12 minor · reviewed 2026-07-08 · glm-5.2
Pith's one-line read ULX variability mirrors AGN patterns across mass scales
desk verdict Competent 24-year variability census for five ULXs in NGC 4631; the headline AGN-analogy claim rests on three data points and is fragile. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The argument rests on two variability diagnostics: (1) normalized excess variance σ²_rms = [1/(N-1)x̄²] Σ(xᵢ-x̄)² - [1/Nx̄²] Σσ²_err,i, which subtracts measurement noise from total observed variance to isolate intrinsic source variability, normalized by the square of mean luminosity; and (2) the structure function SF²(Δt) = ⟨[log F(t+Δt)/F(t)]²⟩, fitted as a power law SF² = A²·(Δt)^γ, where A is the short-timescale variability amplitude and γ is the slope describing how variability accumulates with time. Together these two statistics, computed from 24 years of multi-instrument X-ray monitoring, allow the authors to place ULXs on the same variability diagnostic plane previously occupied only
What would settle it
If a larger sample of ULXs shows no anti-correlation between σ²_rms and luminosity, or if structure functions do not increase with timescale in the majority of sources, the parallel to AGN variability would dissolve. More immediately, if revised absorption corrections for X4 or recalibrated Swift count rates for X3 change either the σ²_rms or mean luminosity values, the three-point trend could vanish or reverse.
Extended reading notes
Core claim
The central finding is that three ULXs in NGC 4631 exhibit variability patterns — an anti-correlation between normalized excess variance and luminosity, and a growth of structure function with timescale — that qualitatively mirror established trends in AGN X-ray variability. The key statistical objects carrying this argument are σ²_rms (a dimensionless measure of intrinsic flux variability amplitude) and SF² (a measure of how flux differences scale with the time interval between observations). The anti-correlation between σ²_rms and luminosity means that among these three sources, the ones shining brighter show proportionally smaller fluctuations. The increase of SF² with timescale for X3 (γ
Load-bearing premise
The claim that ULX variability parallels AGN variability rests on three sources, each contributing a single σ²_rms data point, making the anti-correlation a trend across three points rather than a fitted relation with a significance test. The authors acknowledge this limitation, but the abstract and summary present the result more confidently. If even one source's variability or luminosity were revised due to absorption corrections or calibration differences, the qualitative
Editorial extensions
If this is right
- If the σ²_rms–luminosity anti-correlation is universal across ULXs and AGNs, structure function slopes could serve as a tool to estimate black hole masses in ULXs by analogy with AGN scaling relations, providing mass estimates where dynamical measurements are impossible.
- The ~311-day possible periodicity in X5, if confirmed by future monitoring, would imply a binary orbital or super-orbital timescale that could constrain the mass ratio and geometry of the system.
- The distinct structure function behaviors — X3 with steep slope (γ≈0.61) but low short-term amplitude, X4 with high short-term amplitude but flat slope (γ≈0.09) — suggest different accretion geometries or compact object types among ULXs, motivating classification schemes based on variability timescale rather than spectral shape alone.
- X1's identification as a white dwarf nuclear-burning source, if correct, extends the ULX phenomenon to accretors well below stellar-mass black holes, broadening the mass range over which super-Eddington accretion physics applies.
- The luminosity–hardness correlation in X5 (harder when brighter) supports super-Eddington funnel beaming models where increased accretion rate narrows the wind funnel and enhances the hard Comptonized component.
Reading between the lines
- The three-point σ²_rms–luminosity trend could be a selection effect rather than a physical relation: with only one data point per source, the trend conflates inter-source differences (different black hole masses, accretion rates, viewing angles) with the intra-source variability–luminosity relation seen in AGNs. A robust test would require measuring σ²_rms at multiple luminosity states within indi
- If ULXs and AGNs share variability physics, the structure function slope γ should systematically depend on black hole mass, with more massive objects showing shallower slopes (slower variability growth). Testing this would require a heterogeneous sample spanning stellar-mass through intermediate-mass to supermassive black holes with consistent SF analysis.
- The X4 rapid variability (hours timescale, flat SF slope) resembles behavior seen in neutron star ULXs with accretion column instabilities, raising the possibility that SF slope could discriminate between black hole and neutron star accretors in ULXs.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript presents a systematic long-term X-ray variability study of five ULXs (X1–X5) in NGC 4631, using 37 archival Chandra, XMM-Newton, and Swift observations spanning 24 years. The authors construct light curves, perform spectral modeling, compute normalized excess variance (σ²_rms), and apply structure function (SF) analysis to characterize variability on timescales from days to decades. They report that σ²_rms is anti-correlated with mean luminosity for three sources (X3, X4, X5) and that SF² increases with timescale for X3 and X5, drawing parallels to AGN variability behavior. The data reduction and spectral fitting follow standard procedures (CIAO, SAS, XSPEC), and the methods are correctly implemented per cited references. The central variability claims, however, rest on a very small number of data points and the manuscript does not adequately address cross-instrument systematics or test the sensitivity of its conclusions to individual observations.
Significance. The paper provides a useful compilation of 24 years of multi-mission X-ray data for the NGC 4631 ULXs and applies standard variability tools (σ²_rms, SF) to a sample that has not been systematically studied in this way before. The individual source descriptions (§5.1) are informative, particularly the discussion of X1's supersoft transient behavior and X4's rapid flux changes. The comparison to AGN scaling relations, while speculative given the sample size, is a legitimate motivating framework. The spectral fitting results are reproducible from the tabulated parameters. However, the significance of the central statistical claims is limited by the three-source sample and the absence of formal correlation tests or sensitivity analyses.
major comments (5)
- §4.3, Eq. (2), Table 4: The σ²_rms anti-correlation with luminosity is the paper's central quantitative claim, but it rests on only three data points (one per source: X3, X4, X5), with no formal correlation test (e.g., Pearson or Spearman coefficient, p-value). The claim is presented in the abstract and summary as an established result ('confirming that higher-luminosity ULXs exhibit smaller amplitude of flux fluctuations'), which overstates the evidence. A three-point trend without a significance test does not constitute confirmation. The authors should either downgrade the language to 'suggestive trend' or provide a statistical assessment, and the abstract/summary should be revised accordingly.
- §4.3, Table 4 (X3, Swift ObsID 00082263005): X3's σ²_rms is likely dominated by a single Swift observation with only 26 net full-band counts. In this observation, Γ_eff = 1.18±0.52 is derived from a hardness ratio with HR = 0.07±0.21, an uncertainty spanning nearly the full physically plausible range. The resulting luminosity (L_X = 2.5×10³⁹ erg/s) is a factor of 5–8× higher than every Chandra and XMM-Newton measurement of X3 (0.3–0.8×10³⁹; Tables 4, 9). This single point inflates both X3's mean luminosity and its variance. The measurement-error term σ²_err,i in Eq. (2) captures only statistical uncertainty on the luminosity, not the systematic uncertainty from the poorly constrained Γ_eff-to-flux conversion. The paper does not test the sensitivity of σ²_rms to this point—for instance, by fixing Γ_eff to the spectral-fitting value (Γ ≈ 1.0–1.3; Table 9) and recomputing. If this point is,
- §4.3, Table 4 (X3, Swift ObsID 00082263005), continued: revised downward, X3's σ²_rms could change substantially, potentially breaking the monotonic three-point anti-correlation. The authors should perform a sensitivity test by excluding or revising this single Swift data point and reporting the effect on the σ²_rms–luminosity trend. This is distinct from a generic small-sample concern: it identifies the specific observation and mechanism through which the central claim is most vulnerable.
- §4.4, Eqs. (4)–(6): The structure function is computed by combining luminosities from three instruments (Chandra, XMM-Newton, Swift) with different sensitivities, PSFs, and energy responses. Cross-instrument calibration systematics could introduce artificial scatter or bias the SF²–Δt relation, particularly for X3 and X5 where Swift data with poorly constrained Γ_eff contribute. The paper does not discuss this potential systematic or test whether the SF results are robust when restricted to a single instrument (e.g., Chandra-only). A check on the Chandra-only subset would strengthen the claim that the SF² increase with Δt is intrinsic rather than instrument-dependent.
- §4.4, Figures 6–8: The SF power-law slopes are reported as γ ≈ 0.61±0.24 (X3), γ ≈ 0.09±0.05 (X4), and a value for X5 that is not clearly stated in the text. The fit uses an unweighted least-squares method on binned log(SF²) values, but the binning criteria (minimum 10 points per bin, 0.1 dex minimum width) may result in very few bins for sources with sparse data. The number of bins used for each fit, the reduced χ² values, and the residuals are not tabulated. Without this information, the goodness of fit cannot be assessed. Please report the number of bins, χ²_red, and residual statistics for each SF fit.
minor comments (12)
- Abstract: 'The structure function values increases' should be 'increase' (subject-verb agreement).
- §1, paragraph 3: 'three additional ULXs (X6, X7, and X8) have been identified in NGC 4631 with new XMM-Newton data taken in 2025 (Ducci et al. 2025; Allak et al. 2026)' — the reference to Allak et al. 2026 appears to be a future-dated arXiv preprint; please verify the citation format and availability.
- §2: There is a stray period before 'Circular regions centered on the known source coordinates' (beginning of the paragraph describing extraction regions).
- Table 3 caption: 'The table lists, for each observation, the mission name, observation ID...' — this caption is repeated for Tables 4–6 via 'The columns have the same definitions as those in Table 3,' but Tables 4–6 include XMM-Newton data with different instrument configurations (MOS1, MOS2, PN) that are not mentioned in the Table 3 caption. Consider adding a note about the XMM-Newton detector identifiers.
- Table 9: The flux units are stated as 10⁻¹⁴ erg cm⁻² s⁻¹, whereas Tables 8 and 11 use 10⁻¹³. This inconsistency should be noted or corrected for clarity.
- §4.1: 'Source X3 reached the ULX luminosity threshold (10³⁹ erg s⁻¹) in a single Swift observation (ObsID = 00082263005)' — this observation is from 2013, not 2018 as stated later in §5.1.3 ('one Swift observation in 2018'). Please correct the date.
- §4.2: The Pearson correlation analysis for X5 excludes the Cycle 1 outlier (ObsID 797) to improve the correlation. The justification (ACIS-S quantum efficiency degradation) is reasonable, but the practice of post-hoc outlier removal before correlation testing should be noted explicitly as a caveat, and the result without removal should also be mentioned.
- §5.1.1: 'the thermonuclear burning on the surface of the white dwarf causes the photosphere to expand dramatically to a radius of ~10⁹ cm' — this radius estimate is derived from the blackbody assumption, but the text does not explicitly state the luminosity and temperature used. Please provide the specific values used in the calculation.
- §5.2: 'If a fraction of the ULX population are indeed scaled-down analogs of AGNs' — the grammar should be 'is' rather than 'are' (fraction is singular).
- Figures 1–8: The figure labels and annotations appear to have rendering issues (unicode escape sequences visible in the text). Please ensure figures are properly rendered in the final version.
- §3.1: 'We used PIMMS v4.15 to construct the HR–Γ_eff relation' — the acronym PIMMS should be expanded (Portable, Interactive, Multi-Mission Simulator) at first use, or referenced more formally.
- References: Several arXiv preprint references (e.g., Allak et al. 2026, Salvaggio et al. 2022) should be checked for published versions.
Simulated Author's Rebuttal
We thank the referee for a thorough and constructive report. The referee raises four major points: (1) the σ²_rms–luminosity anti-correlation rests on three data points with no formal correlation test and the language overstates the evidence; (2) X3's σ²_rms may be dominated by a single low-count Swift observation with a poorly constrained Γ_eff, and a sensitivity test is needed; (3) the structure function combines data from three instruments without addressing cross-instrument systematics or testing robustness with a single-instrument subset; and (4) SF fit diagnostics (number of bins, χ²_red, residuals) are not reported. We agree with all four points and will revise the manuscript accordingly. Specifically, we will downgrade the language from 'confirming' to 'suggestive,' add Spearman/Pearson coefficients where applicable, perform the requested sensitivity test on the X3 Swift data point, conduct a Chandra-only SF check, and tabulate all SF fit diagnostics. One partial limitation: for some sources the Chandra-only SF subset is too sparse to yield a meaningful power-law fit, and we will report this transparently.
read point-by-point responses
-
Referee: §4.3, Eq. (2), Table 4: The σ²_rms anti-correlation with luminosity rests on only three data points (X3, X4, X5) with no formal correlation test. The language in the abstract and summary ('confirming') overstates the evidence.
Authors: The referee is correct on both counts. With only three data points, no correlation test can yield a statistically meaningful p-value, and the word 'confirming' in the abstract and summary overstates the evidence. We will revise the manuscript as follows. First, we will replace 'confirming' with 'suggesting' or 'indicating a suggestive trend' throughout the abstract, Section 4.3, and the Summary (Section 6). Second, we will add an explicit caveat that N=3 precludes a robust statistical assessment. Third, although we recognize that a formal Spearman or Pearson test on three points carries limited inferential power, we will report the coefficients and p-values for completeness and transparency, so the reader can judge the evidence directly. We will also note in the text that this trend is consistent with the larger-sample results of González-Martín et al. (2011), which provides external support but does not substitute for a larger ULX sample. revision: yes
-
Referee: §4.3, Table 4 (X3, Swift ObsID 00082263005): X3's σ²_rms is likely dominated by a single Swift observation with only 26 net counts, Γ_eff = 1.18±0.52 from HR = 0.07±0.21, and L_X = 2.5×10³⁹ erg/s, a factor of 5–8× higher than all Chandra/XMM measurements. The measurement-error term captures only statistical uncertainty, not systematic uncertainty from the Γ_eff-to-flux conversion. A sensitivity test excluding or revising this point is needed.
Authors: This is a well-identified and specific concern. We agree that the Swift observation 00082263005 is the single most influential data point for X3's σ²_rms and that the Γ_eff uncertainty introduces a systematic error not captured by the σ²_err,i term in Equation (2). We will perform the requested sensitivity test in the revised manuscript. Specifically, we will recompute X3's σ²_rms and mean luminosity under two alternative assumptions: (a) excluding the Swift data point entirely, and (b) fixing Γ_eff to the spectral-fitting value (Γ ≈ 1.0–1.3 from Table 9) and recomputing the Swift luminosity. We will report the resulting σ²_rms values and show whether the three-point anti-correlation with X4 and X5 is preserved or broken. We will also add a discussion of the systematic uncertainty from the Γ_eff-to-flux conversion for low-count Swift observations. We note that the referee's comment was truncated in transmission, but the concern about the sensitivity of the σ²_rms–luminosity trend to this single point is clear, and we will address it fully. revision: yes
-
Referee: §4.4, Eqs. (4)–(6): The structure function combines luminosities from three instruments with different sensitivities, PSFs, and energy responses. Cross-instrument calibration systematics could introduce artificial scatter or bias the SF²–Δt relation. A Chandra-only subset check is needed.
Authors: We agree that cross-instrument systematics are a legitimate concern when combining Chandra, XMM-Newton, and Swift data in the structure function. We will add a discussion of this potential systematic in Section 4.4. We will also perform the requested Chandra-only SF analysis for X3, X4, and X5. We note that X3 has 11 Chandra detections, X4 has 10, and X5 has 12, which may yield enough pairs for a limited SF analysis, though the time baseline is shorter than the full multi-instrument dataset (Chandra data span 2000–2023, but the 2000 observation is a single epoch and the bulk of Chandra monitoring is from 2022–2023). We will report the Chandra-only SF results where feasible and state explicitly where the data are too sparse for a meaningful power-law fit. For X5, which has the most Chandra detections and the widest Chandra luminosity dynamic range, we expect the most informative single-instrument test. If the Chandra-only SF² still increases with Δt, this will strengthen the claim that the trend is intrinsic. revision: yes
-
Referee: §4.4, Figures 6–8: SF power-law slopes are reported but the number of bins, reduced χ² values, and residuals are not tabulated. The goodness of fit cannot be assessed. Please report the number of bins, χ²_red, and residual statistics for each SF fit.
Authors: The referee is correct that these fit diagnostics are essential for assessing goodness of fit and are currently missing from the manuscript. We will add a table (or expand the figure captions) reporting, for each source's SF fit: the number of bins used in the fit, the number of data pairs per bin, the reduced χ², and the residual statistics. We will also state the X5 slope value explicitly in the text, which is currently only shown in the figure. Additionally, we will verify that the binning criteria (minimum 10 points per bin, 0.1 dex minimum width) are clearly described and will report how many bins resulted from the merging procedure for each source. If any fit has very few bins (e.g., 2–3), we will note that the power-law slope is poorly constrained and interpret it accordingly. revision: yes
Circularity Check
No circularity: σ²_rms and SF² are computed from observed luminosities via standard formulas; AGN comparison uses external literature.
full rationale
The paper computes σ²_rms (Eq. 2) and SF² (Eqs. 4–6) directly from observed X-ray luminosities using standard formulas from Lanzuisi et al. (2014) and Prokhorenko et al. (2024). No parameter is fitted to a subset of data and then 'predicted' for a closely related quantity. The three-point anti-correlation between σ²_rms and mean luminosity is a direct observation, not a fitted relation repackaged as a prediction. The AGN comparison invokes external literature (Nandra et al. 1997; González-Martín et al. 2011; Prokhorenko et al. 2024) with no author overlap. The spectral model selection follows Soria & Ghosh (2009), but this is a methodological choice re-tested on the expanded dataset, not a self-citation chain that forces the variability results. The derivation chain is self-contained against external benchmarks.
Assumptions & free parameters
free parameters (4)
- Galactic NH (spectral fitting) =
1.2–1.3 × 10^20 cm^-2
- Power-law index Γ_eff (photometry) =
varies per observation, ~0.4–2.4
- SF normalization A =
X3: ~0.0975, X4: ~0.369, X5: not stated
- SF slope γ =
X3: 0.61±0.24, X4: 0.09±0.05, X5: not clearly stated
assumptions (4)
- domain assumption The standard thin-disk spectral model (diskbb) is an adequate approximation for the ULX spectra when photon statistics are sufficient for χ² fitting.
- domain assumption The power-law spectral model with Galactic absorption adequately describes the X5 spectrum across all epochs.
- domain assumption The structure function follows a power-law relation with time interval (Eq. 5).
- domain assumption Cross-instrument calibration between Chandra, XMM-Newton, and Swift is sufficiently consistent to combine luminosities in a single light curve.
Cite this review
Pith. "Pith review of The X-ray Variability of the Ultraluminous X-ray Sources in the NGC 4631 galaxy." pith.science (2026). https://pith.science/paper/DDZ2UICZ
@misc{pith2026260706156,
author = {Pith},
title = {Pith review of: The X-ray Variability of the Ultraluminous X-ray Sources in the NGC 4631 galaxy},
year = {2026},
howpublished = {\url{https://pith.science/paper/DDZ2UICZ}},
note = {Machine review of arXiv:2607.06156}
}
read the original abstract
We perform a systematic study on the long-term X-ray variability for the five ultraluminous X-ray sources (ULXs) in the NGC 4631 galaxy (X1-X5), using Chandra, XMM-Newton, and Swift observations covering a 24-year span. Light curves for the five ULXs are presented, while X-ray spectra were modeled for observations with sufficient counts. The normalized excess variance and structure function are utilized to study the X-ray variability behavior of the ULXs. The normalized excess variance is anti-correlated with average X-ray luminosity for three ULXs, indicating that objects with higher average luminosity tend to exhibit relatively lower amplitude of variability. The structure function values increases with time interval in two sources, showing that flux variations become more significant for longer timescales. These trends are also found in the X-ray variability of active galactic nuclei (AGNs). The similarity between ULXs and AGNs, if confirmed for a larger sample of sources, possibly indicates similar underlying physical mechanisms for their X-ray variability.
Figures
Figures from the paper (5 more)
Reference graph
Works this paper leans on
-
[1]
2026, arXiv e-prints, arXiv:2601.08047 2
Allak, S., Akyuz, A., Aladag, Y ., Ducci, L., & Santangelo, A. 2026, arXiv e-prints, arXiv:2601.08047 2
-
[2]
1989, Geochimica et Cosmochimica acta, 53, 197 7
Anders, E., & Grevesse, N. 1989, Geochimica et Cosmochimica acta, 53, 197 7
work page 1989
-
[3]
1996, in Astronomical data analysis software and systems v, V ol
Arnaud, K. 1996, in Astronomical data analysis software and systems v, V ol. 101, 17 5
work page 1996
- [4]
-
[5]
Ducci, L., Mereghetti, S., Pintore, F., et al. 2025, ApJ, 994, L38 2
work page 2025
-
[6]
Farrell, S. A., Webb, N. A., Barret, D., Godet, O., & Rodrigues, J. M. 2009, Nature, 460, 73 1
work page 2009
- [7]
-
[8]
Fruscione, A., McDowell, J. C., Allen, G. E., et al. 2006, in Observatory Operations: Strategies, Processes, and Systems, V ol. 6270, SPIE, 586 3
work page 2006
Show all 32 references
-
[9]
2004, in Astronomical Data Analysis Software and Systems (ADASS) XIII, V ol
Gabriel, C., Denby, M., Fyfe, D., et al. 2004, in Astronomical Data Analysis Software and Systems (ADASS) XIII, V ol. 314, 759 3
2004
-
[10]
e., et al
Gehrels, N., Chincarini, G., Giommi, P. e., et al. 2004, The Astrophysical Journal, 611, 1005 3 Gierli´nski, M., & Done, C. 2004, MNRAS, 347, 885 19
2004
-
[11]
C., Roberts, T
Gladstone, J. C., Roberts, T. P., & Done, C. 2009, MNRAS, 397, 1836 1, 20 Gonz´alez-Mart´ın, O., Papadakis, I., Reig, P., & Zezas, A. 2011, A&A, 526, A132 2, 14, 20
2009
-
[12]
2023, ApJ, 946, 72 2, 18
Guo, J., Wu, J., Feng, H., et al. 2023, ApJ, 946, 72 2, 18
2023
-
[13]
2001, Astronomy & Astrophysics, 365, L1 3
Jansen, F., Lumb, D., Altieri, B., et al. 2001, Astronomy & Astrophysics, 365, L1 3
2001
-
[14]
Kaaret, P., Feng, H., & Roberts, T. P. 2017, ARA&A, 55, 303 1
2017
-
[15]
2002, Monthly Notices of the Royal Astronomical Society, 337, L11 20
Kubota, A., Done, C., & Makishima, K. 2002, Monthly Notices of the Royal Astronomical Society, 337, L11 20
2002
-
[16]
2001, The Astrophysical Journal Letters, 560, L147 20
Kubota, A., Makishima, K., & Ebisawa, K. 2001, The Astrophysical Journal Letters, 560, L147 20
2001
-
[17]
2014, ApJ, 781, 105 13
Lanzuisi, G., Ponti, G., Salvato, M., et al. 2014, ApJ, 781, 105 13
2014
-
[18]
P., & Done, C
Luangtip, W., Roberts, T. P., & Done, C. 2016, Monthly Notices of the Royal Astronomical Society, 460, 4417 20
2016
-
[19]
2007, in IAU Symposium, V ol
Makishima, K. 2007, in IAU Symposium, V ol. 238, Black Holes from Stars to Galaxies – Across the Range of Masses, ed. V . Karas & G. Matt, 209 9
2007
-
[20]
1983, Astrophysical Journal, Part 1 (ISSN 0004-637X), vol
Morrison, R., & McCammon, D. 1983, Astrophysical Journal, Part 1 (ISSN 0004-637X), vol. 270, July 1, 1983, p. 119-122., 270, 119 7
1983
-
[21]
M., Mushotzky, R
Nandra, K., George, I. M., Mushotzky, R. F., Turner, T. J., & Yaqoob, T. 1997, ApJ, 488, L91 2, 20
1997
-
[22]
A., Sazonov, S
Prokhorenko, S. A., Sazonov, S. Y ., Gilfanov, M. R., et al. 2024, MNRAS, 528, 5972 2, 14, 15, 20
2024
-
[23]
N., et al
Pu, X., Luo, B., Brandt, W. N., et al. 2020, ApJ, 900, 141 3, 4
2020
-
[24]
J., et al
Robba, A., Pinto, C., Walton, D. J., et al. 2021, A&A, 652, A118 19
2021
-
[25]
2022, MNRAS, 512, 1814 1
Salvaggio, C., Wolter, A., Pintore, F., et al. 2022, MNRAS, 512, 1814 1
2022
-
[26]
E., Narayan, R., et al
Shafee, R., McClintock, J. E., Narayan, R., et al. 2006, ApJ, 636, L113 19
2006
-
[27]
Soria, R., & Ghosh, K. K. 2009, ApJ, 696, 287 2, 6, 7, 9, 17, 18, 19
2009
-
[28]
B., Courtois, H
Tully, R. B., Courtois, H. M., Dolphin, A. E., et al. 2013, AJ, 146, 86 2
2013
-
[29]
J., George, I
Turner, T. J., George, I. M., Nandra, K., & Turcan, D. 1999, ApJ, 524, 667 2, 20
1999
-
[30]
J., Harrison, F
Walton, D. J., Harrison, F. A., Grefenstette, B. W., et al. 2014, ApJ, 793, 21 1
2014
-
[31]
2012, Science, 337, 554 1
Webb, N., Cseh, D., Lenc, E., et al. 2012, Science, 337, 554 1
2012
-
[32]
2002, Publications of the Astronomical Society of the Pacific, 114, 1 3
Weisskopf, M., Brinkman, B., Canizares, C., et al. 2002, Publications of the Astronomical Society of the Pacific, 114, 1 3
2002
Reviewed July 8, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.