Pith. sign in

REVIEW 3 major objections 5 minor 1 cited by

Variability of the near-infrared extinction curve towards the Galactic centre

T0 review · 3 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

Pith's one-line read The near-infrared extinction curve toward the Galactic centre is not a single power law: the extinction index is 2.43 ± 0.03 for J–H but 2.23 ± 0.03 for H–Ks, Δα = 0.19 ± 0.05, with no line-of-sight or extinction variation.

desk verdict Plausible Δα, but the 4σ significance rests on an unquantified assumption about the red clump's intrinsic color–magnitude slope. read the letter →

arxiv 1909.02494 v1 pith:SZCX42KN submitted 2019-09-05 astro-ph.SR astro-ph.GA

classification astro-ph.SRastro-ph.GA
keywords near-infraredextinctionindexGalacticcentreredclumpstarsnuclearstellardiscinterstellardustGALACTICNUCLEUSsurvey
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper aims to establish that the near-infrared extinction curve toward the Galactic centre — the law that converts observed $JHK_s$ photometry into intrinsic starlight — cannot be represented by a single power-law index. Using roughly 62,000 red clump stars from the high-angular-resolution GALACTICNUCLEUS catalogue and two independent methods, it finds the extinction index between $J$ and $H$ is $\alpha_{JH} = 2.43 \pm 0.03$ while between $H$ and $K_s$ it is $\alpha_{HK_s} = 2.23 \pm 0.03$, a difference of $\Delta\alpha = 0.19 \pm 0.05$. The same analysis finds no significant variation of the index across the nuclear stellar disc or with the amount of extinction itself. The result matters because extinction corrections built on a single $\alpha$ bias distance estimates to the inner Galaxy by up to about a kiloparsec and can change the inferred stellar type of a dereddened star.

What carries the argument

The argument runs on the extinction-index definition $A_\lambda \propto \lambda^{-\alpha}$ and on red clump stars as a standard stellar population whose spread along the reddening vector traces differential extinction. Two independent estimators are combined. The grid method reddens synthetic red clump magnitudes, built from atmosphere models with adopted effective temperature, surface gravity, radius, metallicity and a Galactic-centre distance of $8.0 \pm 0.1$ kpc, over a fine grid of extinction $A_{1.61}$ and index $\alpha$, and minimises the corresponding $\chi^2$. The slope method instead measures the slope $m$ of the red clump feature in colour–magnitude diagrams and converts it to an extinction index through the identity $\alpha = -\log(1 + 1/m) / \log(\lambda_{\mathrm{eff},1}/\lambda_{\mathrm{eff},2})$, with $\lambda_{\mathrm{eff}}$ the effective wavelengths of the filters. The near-agreement of the two methods' $\Delta\alpha$, despite their different systematic sensitivities, is what lets the paper quote $\Delta\alpha = 0.19 \pm 0.05$ as its best estimate.

What would settle it

Observe the red clump feature in $K_s$ versus $H-K_s$ in a low-extinction field where differential extinction is negligible: a nonzero slope there would prove an intrinsic colour–magnitude correlation and break the slope method's premise. A complementary check is to compare spectroscopically measured extinctions of individual red clump stars with the values predicted by $\alpha_{JH} = 2.43$ and $\alpha_{HK_s} = 2.23$; a systematic offset that grows with colour would falsify the band-dependent indices.

Watch

Extended reading notes

Core claim

The central claim is that the NIR extinction index toward the Galactic centre depends on wavelength: $\alpha_{JH} = 2.43 \pm 0.03$ and $\alpha_{HK_s} = 2.23 \pm 0.03$, giving $\Delta\alpha = 0.19 \pm 0.05$ when the two methods are combined. The slope of the bright red clump feature alone yields $\Delta\alpha = 0.19 \pm 0.05$, a detection at roughly $4\sigma$ significance, and the fully independent grid method returns $\Delta\alpha = 0.22 \pm 0.13$; the authors take the agreement between these two systematically different estimators as the evidence for a real wavelength dependence. They also find the index constant, within uncertainties, across the studied region of the nuclear stellar disc and along the reddening vector, so the wavelength dependence is not a line-of-sight artefact. The paper further argues that this dependence reconciles the wide spread of $\alpha$ values in the literature (roughly 1.5 to 2.5), because a single-index fit over $JHK_s$ measures an average, here $\alpha_{JHK_s} = 2.32 \pm 0.09$, which sits between the two band-specific indices.

Load-bearing premise

The headline value $\Delta\alpha = 0.19 \pm 0.05$ rests on the slope method's premise that the slope of the red clump feature in the colour–magnitude diagram is caused entirely by differential extinction, with no intrinsic colour–magnitude correlation among red clump stars and with correctly adopted effective wavelengths; if either assumption fails, both $\alpha_{JH}$ and $\alpha_{HK_s}$ shift and $\Delta\alpha$ changes.

Editorial extensions

If this is right

  • A single-index extinction correction toward the Galactic centre is systematically biased: fitting one power law across JHKs returns an average index (2.32 ± 0.09 in this paper) that lies between the true J–H and H–Ks indices.
  • Distance moduli to the nuclear disc derived from NIR photometry can be biased by up to about 1 kpc if the wavelength dependence is ignored, since a 10–15% change in the extinction index shifts the absolute extinction by roughly 0.3 mag.
  • Stellar classifications that use dereddened NIR colours can change type when a band-dependent extinction law is applied instead of a single index.
  • Within the studied region of the nuclear stellar disc, a single constant extinction law can safely be adopted for structural studies, as long as the band pair used is specified.
  • Published extinction indices toward the Galactic centre spanning roughly 1.5 to 2.5 need not contradict one another; they may reflect different effective band pairs convolved with a wavelength-dependent law.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the wavelength dependence is a property of the dust rather than of this particular sightline, a similar Δα between adjacent NIR bands should show up toward other heavily extincted inner-Galaxy fields, which could be checked by applying the slope method to existing catalogues.
  • A natural physical cause would be a change in the dust grain population; comparing Δα across sightlines with different extinction curves would turn the effect into a probe of grain size or composition that needs no absolute photometric calibration.
  • Because the slope method's output depends on the assumed effective wavelengths, calibrating intrinsic red clump colours spectroscopically in a low-extinction field would directly test the premise that the observed red clump slope is purely extinction-driven and could sharpen Δα below the current 0.05.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper uses the GALACTICNUCLEUS JHK_s survey of 14 fields covering ~90 pc x 20 pc around Sgr A* to study the near-infrared extinction index alpha toward the Galactic centre. Two red-clump-based methods are applied: a 'grid method' that reddens synthetic stellar atmosphere models over a grid of extinction and alpha, and a 'slope method' that measures the slope of the red-clump feature in colour-magnitude diagrams and converts it to alpha via Eq. (1). The authors report no significant line-of-sight or extinction dependence of alpha within the nuclear stellar disc, but a wavelength dependence Delta_alpha = alpha_JH - alpha_HKs = 0.19 +/- 0.05, with mean values alpha_JH = 2.43 +/- 0.03 and alpha_HKs = 2.23 +/- 0.03. They interpret the result as evidence that the NIR extinction curve towards the Galactic centre is not a single power law, and suggest this explains discrepancies in previously published alpha values.

Significance. If the result is robust, it is an important step for Galactic-centre studies: a wavelength-dependent NIR extinction index affects distance estimates and stellar classifications based on NIR photometry, and it offers a natural explanation for the spread of alpha values in the literature. The two methods agree on the sign and rough magnitude of Delta_alpha, and the high angular resolution of the GALACTICNUCLEUS survey is a genuine strength. However, the headline significance rests heavily on the slope method (Section 4.4), whose central assumption is not quantified, and the effective wavelengths enter the conversion without a propagated uncertainty. The paper is concise and readable, but the current form does not fully secure the central claim.

major comments (3)
  1. [Section 4.4, Eq. (1)] The slope method interprets the observed slope m of the red-clump feature as entirely due to differential extinction, i.e. it assumes a zero intrinsic colour-magnitude correlation for red-clump stars. The paper's best estimate, Delta_alpha = 0.19 +/- 0.05, is carried by this assumption, since the grid method alone yields only Delta_alpha = 0.22 +/- 0.13, a ~1.7 sigma signal. The manuscript's only caveat (Section 5) concerns completeness at the faint end and does not bound the intrinsic slope. A nonzero intrinsic slope, possibly different between the (J-H, H) and (H-Ks, Ks) planes, could shift alpha_JH and alpha_HKs in opposite directions and mimic a wavelength dependence even for a single true extinction index. Please quantify the intrinsic RC slope using synthetic populations or external estimates, or present the grid method as the primary result with the slope method as supporting evidence.
  2. [Section 4.4, Eq. (1)] The conversion from the measured slope to alpha uses adopted effective wavelengths lambda_eff, which enter logarithmically in the denominator of Eq. (1). The paper does not propagate uncertainties in lambda_eff or in the filter transmission curves into the quoted uncertainties on alpha_JH, alpha_HKs, or Delta_alpha. A systematic error in the ratio lambda_J/lambda_H relative to lambda_H/lambda_Ks would produce exactly the reported wavelength dependence. Please provide a quantitative estimate of this systematic term, or state explicitly why it is negligible compared with the quoted 0.05 uncertainty on Delta_alpha.
  3. [Sections 4.1-4.2 and Table 1] The grid method's systematic uncertainties are computed by varying parameters 'independently in their uncertainty ranges', but the adopted red-clump template parameters (Kurucz 1993 models, T_eff = 4750 K, log g = +2.5, R = 10 R_sun, twice-solar metallicity) are fixed inputs. The sensitivity of the derived Delta_alpha to the template choice, especially metallicity and effective temperature, is not shown. Since both methods use the same red-clump selection from the same survey data, their agreement is not a fully independent confirmation; please add a table or figure showing how Delta_alpha changes over the physically allowed parameter ranges.
minor comments (5)
  1. [Section 3] The text refers to 'blue dashed trapezoids' while Section 4.1 and the caption of Fig. 2 use 'blue dashed parallelograms'; please unify the terminology.
  2. [Section 4.4, Eq. (1)] Please state the sign convention for the slope m (e.g., m = d lambda_eff2 / d(lambda_eff1 - lambda_eff2)) and whether positive slopes correspond to features aligned with the reddening vector, so that readers can reproduce the calculation without consulting Nogueras-Lara et al. (2018b).
  3. [Table 3] Table 3 reports statistical and systematic uncertainties separately, but the combined values quoted in the text (alpha_HKs = 2.29 +/- 0.02, alpha_JH = 2.52 +/- 0.09, Delta_alpha = 0.23 +/- 0.09, and Delta_alpha = 0.19 +/- 0.05) do not state how the two components were combined; please make the propagation explicit.
  4. [Abstract and body] The abstract contains a 'Conclusions' heading with no following text, and the body has no separate conclusions section; either add a brief summary or remove the empty heading.
  5. [References] The references to 'Nogueras-Lara et al., submitted' should be updated to published or arXiv identifiers, and the two distinct submitted papers should be clearly distinguished.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the extinction-index measurements are direct inversions of photometric data, not fitted inputs renamed as predictions.

full rationale

The paper derives the NIR extinction indices from two independent methods rather than from a parameter fitted to the claimed conclusion. In the grid method (Sect. 4.1) the authors redden Kurucz synthetic stellar models over a grid of extinction and alpha values and minimize chi^2, so alpha is an output of an external stellar-atmosphere model and photometric data. In the slope method (Sect. 4.4) the observed slope m of the red-clump feature is converted to alpha via Eq. (1), alpha = -log(1 + 1/m) / log(lambda_eff1/lambda_eff2); this is a direct algebraic inversion of the assumed power-law extinction curve, not a fit that has been relabeled as a prediction. The two methods use the same survey data but different observables, and their agreement (Delta_alpha = 0.22 +/- 0.13 versus 0.19 +/- 0.05) is presented as a consistency check, not as a construction of the result. Self-citations to Nogueras-Lara et al. (2018a,b) provide the measurement techniques, but those techniques are restated in the present paper and are anchored by external inputs (Kurucz 1993 models, adopted stellar parameters, effective wavelengths); the cited prior work does not itself supply the wavelength-dependence claim as a premise. The assumption that the RC feature slope is entirely due to differential extinction, and possible systematic errors in the adopted effective wavelengths, are genuine astrophysical caveats that could bias alpha_JH and alpha_HKs, but they are not circularity: no equation in the paper defines the result in terms of itself, and no fitted parameter is subsequently presented as an independent prediction. Accordingly the derivation chain is self-contained with respect to the circularity criteria, and the appropriate score is 0.

Assumptions & free parameters 3 free parameters · 4 assumptions · 0 invented entities

The paper's central numbers rest on adopted model inputs (effective wavelengths, red clump parameters, distance) that come from prior literature, plus the standard power-law extinction assumption. None of these are tuned to force the wavelength dependence, but they set the scale of α and are the main source of systematic uncertainty.

free parameters (3)
  • Effective wavelengths λ_eff for J, H, Ks = Not given in the paper; taken from filter curves
    Used in Eq. 1 to convert the measured CMD slope to the extinction index; a small bias in the wavelength ratios changes both α values and Δα.
  • Red clump synthetic model parameters (T_eff, log g, radius, metallicity) = 4750 K, 2.5, 10.0 R_sun, 2x solar
    Used in the grid method (Sect. 4.1) to synthesize intrinsic RC magnitudes; the paper varies these to estimate systematics, but biases in their values would shift the fitted α.
  • Distance to the Galactic center = 8.0 kpc
    Adopted in the grid method from recent literature; affects the absolute magnitude scale and weakly the extinction index.
assumptions (4)
  • domain assumption The NIR extinction curve follows a power law A_λ ∝ λ^{-α} over JHKs
    The entire analysis and the definition of α assume this form; if the curve deviates, the index comparison loses meaning.
  • domain assumption The slope of the red clump feature in the CMD is due entirely to differential extinction
    Load-bearing for the slope method (Sect. 4.4); any intrinsic color-magnitude correlation of RC stars would bias α.
  • domain assumption The two red clump features can be separated with a two-Gaussian mixture and each traces a single stellar population
    Used in Sect. 4.4 to assign slopes to the bright and faint RC features.
  • domain assumption The GALACTICNUCLEUS photometric catalog is complete in the selected CMD region after masking
    The slope method removes faint-end bins and masked regions because completeness is known to drop; the correction assumes the remaining data are complete enough.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Variability of the near-infrared extinction curve towards the Galactic centre." pith.science (2026). https://pith.science/paper/SZCX42KN

@misc{pith2026190902494,
  author       = {Pith},
  title        = {Pith review of: Variability of the near-infrared extinction curve towards the Galactic centre},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/SZCX42KN}},
  note         = {Machine review of arXiv:1909.02494}
}
abstract

Due to the extreme extinction towards the Galactic centre ($A_{V} \sim 30$ mag), its stellar population is mainly studied in the near-infrared (NIR) regime. Therefore, a proper analysis of the NIR extinction curve is necessary to fully characterise the stellar structure and population of the inner part of the galaxy. We studied the dependence of the extinction index ($\alpha_\lambda$) in the NIR on the line of sight, wavelength, and extinction. We used the GALACTICNUCLEUS imaging survey, a high angular resolution catalogue ($0.2''$) for the inner part of the Galaxy in $JHK_s$, and studied the spatial variation in the extinction index. We also applied two independent methods based on red clump stars to compute the extinction index between different bands and its variation with wavelength. We did not detect any significant line-of-sight or extinction variation in $\alpha$ within the studied region in the nuclear stellar disc. The extinction index between $JH$ and $HK_s$ differs by $0.19 \pm 0.05$. We obtained mean values for the extinction indices $\alpha_{JH} = 2.43\pm0.03$ and $\alpha_{HK_s} = 2.23\pm0.03$. The dependence of the extinction index on the wavelength could explain the differences obtained for $\alpha_\lambda$ in the literature since it was assumed constant for the NIR regime.

Figures

Figures reproduced from arXiv: 1909.02494 by the authors.

Figure 1
Figure 1. Image of the studied region produced combining the Ks , H, and J bands in red, green, and blue, respectively. Sagittarius A* and the Arches and Quintuplet clusters are indicated by arrows. The black rectangle near the Arches cluster corresponds to a field with incomplete data. The white dashed contours indicate regions dominated by dark clouds. The cross-shaped region corresponds to a low completeness region due to … view at source ↗
Figure 3
Figure 3. Extinction-index maps: a) JH-map; b) HKs-map; c) JKs-map. Cross-shaped pixels indicate that there are not enough stars for a reliable estimate. 2.3 2.4 2.5 ®JH 0 30 60 90 # pixels ® = 2.39 ¾ = 0.04 2.1 2.2 2.3 2.4 ®HKs ® = 2.21 ¾ = 0.06 2.2 2.3 2.4 ®JKs ® = 2.30 ¾ = 0.03 [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 2
Figure 2. Colour-magnitude diagrams H vs. J − H (upper panel) and Ks vs. H − Ks (lower panel). The RC is marked by the blue dashed paral￾lelograms. The two blue arrows show a double feature in the RC. The black arrow depicts the reddening vector with an extinction AKs = 0.5 mag (computed using α = 2.30 ± 0.08, Nogueras-Lara et al. 2018a). The insets show the RC region with the two features obtained applying GMM and their unce… view at source ↗
Figures from the paper (1 more)
Figure 5
Figure 5. Figure 5: Upper panels: Histograms obtained for αJH (left panel) and αHKs (right panel) using the method presented in Sect. 4.1. Lower panels: Histograms obtained for A1.61, associated with the calculations using JH and HKs . The red lines show the Gaussian fits to the data. The…

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Unveiling the Milky Way with a Gaia DR3 census of OB-type stars within 2 kpc. I. Tracing local Galactic structure, massive star-forming regions and core-collapse supernova progenitors

    astro-ph.GA 2026-07 conditional novelty 4.0 of 10

    A Gaia DR3-based census of 105,971 OB stars within 2 kpc maps local Galactic structure and identifies over 4,200 core-collapse supernova or black hole progenitor candidates.

Reference graph

Works this paper leans on

30 extracted references · 16 canonical work pages · cited by 1 Pith paper

  1. [1]

    1974, Automatic Control, IEEE Transactions on, 19, 716 Alonso-García, J., Minniti, D., Catelan, M., et al

    Akaike, H. 1974, Automatic Control, IEEE Transactions on, 19, 716 Alonso-García, J., Minniti, D., Catelan, M., et al. 2017, ApJ, 849, L13

  2. [2]

    L., Rix, H.-W., et al

    Bovy, J., Nidever, D. L., Rix, H.-W., et al. 2014, ApJ, 790, 127

  3. [3]

    Chaplin, W. J. & Miglio, A. 2013, ARA&A, 51, 353

  4. [4]

    2000, A&AS, 147, 335

    Diolaiti, E., Bendinelli, O., Bonaccini, D., et al. 2000, A&AS, 147, 335

  5. [5]

    2019, arXiv e-prints, arXiv:1907.10731

    Do, T., Hees, A., Ghez, A., et al. 2019, arXiv e-prints, arXiv:1907.10731

  6. [6]

    2018, ApJ, 855, L5

    Do, T., Kerzendorf, W., Konopacky, Q., et al. 2018, ApJ, 855, L5

  7. [7]

    2015, ApJ, 809, 143

    Do, T., Kerzendorf, W., Winsor, N., et al. 2015, ApJ, 809, 143

  8. [8]

    Draine, B. T. 1989, in ESA Special Publication, V ol. 290, Infrared Spectroscopy in Astronomy, ed. E. Böhm-Vitense, 93–98

Show all 30 references
  1. [9]

    2017, MNRAS, 464, 194

    Feldmeier-Krause, A., Kerzendorf, W., Neumayer, N., et al. 2017, MNRAS, 464, 194

  2. [10]

    K., Gillessen, S., Dodds-Eden, K., et al

    Fritz, T. K., Gillessen, S., Dodds-Eden, K., et al. 2011, ApJ, 737, 73

  3. [11]

    2016, ARA&A, 54, 95

    Girardi, L. 2016, ARA&A, 54, 95

  4. [12]

    J., Bandyopadhyay, R

    Gosling, A. J., Bandyopadhyay, R. M., & Blundell, K. M. 2009, MNRAS, 394, 2247 Gravity Collaboration, Abuter, R., Amorim, A., et al. 2018, A&A, 615, L15

  5. [13]

    R., Anderson, J., et al

    Hosek, Matthew W., J., Lu, J. R., Anderson, J., et al. 2018, ApJ, 855, 13

  6. [14]

    2008, A&A, 491, 941

    Kissler-Patig, M., Pirard, J.-F., Casali, M., et al. 2008, A&A, 491, 941

  7. [15]

    Kurucz, R. L. 1993, VizieR Online Data Catalog, 6039, 0

  8. [16]

    W., Bono, G., et al

    Matsunaga, N., Feast, M. W., Bono, G., et al. 2016, MNRAS, 462, 414

  9. [17]

    2003, in Proc

    Nagayama, T., Nagashima, C., Nakajima, Y ., et al. 2003, in Proc. SPIE, V ol. 4841, Instrument Design and Performance for Optical/Infrared Ground-based Telescopes, ed. M. Iye & A. F. M. Moorwood, 459–464

  10. [18]

    2018, MNRAS, 478, 4374

    Nandakumar, G., Ryde, N., Schultheis, M., et al. 2018, MNRAS, 478, 4374

  11. [19]

    M., Udalski, A., Gould, A., & Pinsonneault, M

    Nataf, D. M., Udalski, A., Gould, A., & Pinsonneault, M. H. 2011, ApJ, 730, 118

  12. [20]

    2006, ApJ, 638, 839

    Nishiyama, S., Nagata, T., Kusakabe, N., et al. 2006, ApJ, 638, 839

  13. [21]

    2008, ApJ, 680, 1174

    Nishiyama, S., Nagata, T., Tamura, M., et al. 2008, ApJ, 680, 1174

  14. [22]

    T., et al

    Nogueras-Lara, F., Schödel, R., Gallego-Calvente, A. T., et al. 2019, arXiv e- prints, arXiv:1908.10366

  15. [23]

    2011, Journal of Machine Learning Research, 12, 2825

    Pedregosa, F., Varoquaux, G., Gramfort, A., et al. 2011, Journal of Machine Learning Research, 12, 2825

  16. [24]

    Rieke, G. H. & Lebofsky, M. J. 1985, ApJ, 288, 618

  17. [25]

    Z., Hosek, Matthew W., J., Lu, J

    Rui, N. Z., Hosek, Matthew W., J., Lu, J. R., et al. 2019, arXiv e-prints, arXiv:1904.02395 Schödel, R., Najarro, F., Muzic, K., & Eckart, A. 2010, A&A, 511, A18+ Schödel, R., Yelda, S., Ghez, A., et al. 2013, MNRAS, 429, 1367

  18. [26]

    M., Origlia, L., et al

    Schultheis, M., Rich, R. M., Origlia, L., et al. 2019, A&A, 627, A152

  19. [27]

    1978, The Annals of Statistics, 6, 461

    Schwarz, G. 1978, The Annals of Statistics, 6, 461

  20. [28]

    Z., Stolovy, S

    Scoville, N. Z., Stolovy, S. R., Rieke, M., Christopher, M., & Yusef-Zadeh, F. 2003, ApJ, 594, 294

  21. [29]

    Stead, J. J. & Hoare, M. G. 2009, MNRAS, 400, 731

  22. [30]

    & Gerhard, O

    Wegg, C. & Gerhard, O. 2013, MNRAS, 435, 1874 Article number, page 5 of 5

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.