Pith. sign in

REVIEW 3 major objections 5 minor 39 references

Euclid imaging shows the near-infrared extinction curve in LDN 1641 flattens by up to 27% toward denser regions, consistent with grain growth inside one molecular cloud.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

In LDN 1641 the NIR extinction power-law index α drops by up to 27% from diffuse edges to dense cores, implying substantial grain growth at high AV.

T0 review reviewed 2026-07-11 challenge →

load-bearing objection Solid Euclid Q1 measurement of NIR extinction flattening inside LDN 1641; the grain-growth claim is plausible but rests on small core samples and empirical cuts. the 3 major comments →

arxiv 2607.04286 v1 pith:JCI56JRU submitted 2026-07-05 astro-ph.GA

Euclid Q1 reveals spatial variations of the extinction law in the dense cloud LDN 1641

Rui Chen , Shu Wang , Xiaodian Chen , Kun Wang This is my paper
classification astro-ph.GA
keywords dust extinctionmolecular cloudsnear-infrared extinctiongrain growthLDN 1641Euclid photometryreddening lawOrion A
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The reading

This paper uses deep Euclid Quick Data Release 1 VIS and NISP photometry of the Orion A dark cloud LDN 1641 to measure how starlight is extinguished from the optical through the near infrared along sightlines reaching AV about 30 mag. For the full high-extinction footprint the near-infrared law is a power law Aλ proportional to λ to the minus α with α = 1.57 ± 0.06, giving relative extinctions A_VIS/A_H = 4.23, A_Y/A_H = 2.13, and A_J/A_H = 1.47. When the same analysis is split by local extinction, α is systematically smaller in dense cores than on the cloud outskirts, with differences of about 17% between combined samples and up to 27% between individual subregions. The authors read that flattening as evidence that large grains become more important where the column is highest. A sympathetic reader cares because extinction laws set the accuracy of photometry and distances, and because a clear density-dependent change inside a single nearby cloud would mean grain growth is not only theoretical but observationally mapped from edge to core.

Core claim

From curvature-corrected Euclid (λ−H) versus (Y−H) color–color slopes the authors derive color-excess ratios, convert the near-infrared ratio into a power-law index α under Aλ ∝ λ−α, and obtain α = 1.57 ± 0.06 for the full LDN 1641 study region. High-extinction core subregions give smaller α (combined α ≈ 1.61; individual cores 1.69 and 1.41) than low-extinction reference regions (α ≈ 1.95) or a cleaned 2MASS peripheral sample (α ≈ 2.25). That spatial drop in α means a flatter extinction curve in denser gas and is presented as support for substantial grain growth from the diffuse outskirts to the dense core of one molecular cloud.

What carries the argument

Color-excess ratios E(λ−H)/E(Y−H) measured as the slopes of curvature-corrected (λ−H) versus (Y−H) diagrams; under the assumption Aλ ∝ λ−α the near-infrared ratio fixes α and the relative extinctions Aλ/A_H.

Load-bearing premise

The central claim treats the slope of the observed color–color diagram as the true color-excess ratio after only a model-based filter-wavelength correction, without star-by-star intrinsic colors or spectral types.

What would settle it

Spectroscopic classification of the same Euclid sources that shows the high-extinction cores have a different mix of stellar types or metallicities than the outskirts, such that correcting star-by-star for intrinsic colors erases the reported difference in α, would overturn the grain-growth interpretation.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • Reddening corrections inside dense molecular clouds must allow a spatially varying near-infrared power-law index rather than a single universal α.
  • Smaller α at high AV supports an enhanced large-grain population and therefore coagulation and ice-mantle growth in cloud cores.
  • Average Galactic NIR indices near α ≈ 2 better describe diffuse or cloud-edge sightlines than dense interiors.
  • Relative extinctions such as A_VIS/A_H ≈ 4.23 should be preferred over diffuse-ISM defaults when correcting photometry toward dense Orion cores.
  • Broader Euclid releases can map the same density-dependent flattening over larger Galactic-plane areas and greater depths.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If grain growth this strong is typical, synthetic photometry of young stellar objects may misestimate luminosities unless pipelines adopt environment-dependent extinction curves.
  • Applying the same Euclid color–color slope method to other dark clouds would test whether the 17–27% α swing is universal or tied to local star-formation activity.
  • A spectroscopic campaign that supplies spectral types for the Euclid LDN 1641 sample would cleanly separate residual stellar-population mixing from true grain-size change.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 5 minor

Summary. The paper uses Euclid Q1 VIS and NISP (Y, J, H) photometry of LDN 1641 to measure optical-to-NIR extinction along sightlines reaching high column density. Color-excess ratios are obtained from linear fits to curvature-corrected (λ−H) vs (Y−H) diagrams, converted under a power-law assumption A_λ ∝ λ^{-α} into α and relative extinctions. For the full high-extinction footprint the authors report α = 1.57 ± 0.06, A_VIS/A_H = 4.23 ± 0.24, A_Y/A_H = 2.13 ± 0.18 and A_J/A_H = 1.47 ± 0.11. Subregion analysis and a 2MASS comparison are then used to argue that α is systematically smaller (flatter curve) in denser cores than in low-extinction reference regions, with differences of ~17% between combined samples and up to ~27% between individual subregions, interpreted as grain growth.

Significance. If the spatial variation is real, the work supplies a rare continuous measurement of the NIR extinction law from the outskirts to the dense interior of a single nearby molecular cloud, using Euclid's depth and resolution to reach A_V ~ 30 mag. The full-sample CER pipeline (quality cuts, branch selection, PARSEC-based effective-wavelength curvature correction, linear fits, conversion via Eq. 4) is standard and yields a well-documented average law. The environmental dependence, if robust, would be a direct observational constraint on grain-growth models and a useful calibration for future Euclid extinction work. The 2MASS cross-check usefully anchors the low-extinction end.

major comments (3)
  1. The central claim of spatial variation rests on very small high-core samples (Fig. 4: N = 40 and 53 for cores 1 and 2; N = 93 combined) and on the assumption that the fitted slope of (J−H) vs (Y−H) equals the true CER after only a single PARSEC-based curvature correction (§3.1, Eqs. 1–2). The paper itself notes residual foreground/low-extinction contamination and sparse-N instability as dominant systematics (end of §4). Because the low-extinction reference already has a shallower slope (k = 0.400), incomplete removal of that population from the cores systematically raises high-core k and therefore lowers α (Eq. 4). The 2MASS comparison only reconfirms the low-extinction end and does not protect the high-core measurement. The reported 17–27% flattening could therefore be an artifact of selection purity rather than grain growth. A quantitative contamination test (e.g., controlled injection
  2. Table 1 and §4: the high-core α values (1.69 ± 0.12, 1.41 ± 0.14, combined 1.61 ± 0.09) and the low-extinction value (1.95 ± 0.21) overlap within ~1–2σ once the large low-extinction uncertainty is taken into account. The abstract and §4 language of 'significant spatial variations' and 'systematically smaller α' therefore overstates the statistical weight of the difference. Either a formal significance test that folds in the acknowledged systematics, or a more cautious wording that presents the trend as suggestive, is required.
  3. §3.1 and the free-parameter list: the branch-selection cut (J−H) < 0.42(Y−H)+0.01, the J−H > 0.4 mag foreground cut, the mean (Y−H) thresholds that define high/low subregions, and the fixed PARSEC red-giant parameters used for λ_eff are all free choices that directly affect the measured slopes. Sensitivity of α to reasonable variations of these cuts should be quantified and reported; without that, the load-bearing conversion from slope to α remains under-tested.
minor comments (5)
  1. Figure 2 caption and text: the low-extinction polygons are described as black in the text and white/black in the caption; unify the description.
  2. Equation (3) and surrounding text: state explicitly that A_Y/A_H is obtained from the power-law solution of Eq. (4) before being inserted into Eq. (3), so the logical order is unambiguous.
  3. Table 1: the 2MASS row lists only α; adding the measured E(H−K_S)/E(J−H) = 0.467 ± 0.010 in the table body (already in the notes) would make the comparison self-contained.
  4. Abstract and §1: 'A_V ~ 30 mag' is stated without a direct conversion from the Euclid colors; a brief note on how that estimate is obtained would help the reader.
  5. Several references appear twice or with slightly inconsistent formatting (e.g., Euclid Collaboration entries); a single consistent citation style would improve polish.

Circularity Check

0 steps flagged

Empirical CER slopes converted under a standard power-law ansatz; no load-bearing circular reduction.

full rationale

The derivation chain is: (1) fit slopes of (λ−H) vs (Y−H) color–color diagrams to obtain CERs k_λ = E(λ−H)/E(Y−H); (2) assume the standard NIR form A_λ ∝ λ^{−α} and invert Eq. 4 for α from the measured E(J−H)/E(Y−H); (3) convert CERs to A_λ/A_H via Eq. 3. None of these steps is self-definitional: α is solved from an independent photometric slope, not fitted and then re-predicted. Spatial variation is obtained by repeating the same fit on high- vs low-(Y−H) subregions and on an independent 2MASS (H−K_S)/(J−H) sample; the reported α differences are therefore data-driven comparisons, not forced by construction. Citations to Wang & Chen (2019, 2024) supply only the MW average α for external comparison and methodological precedent for slope-as-CER; they do not enter the LDN 1641 fit or uniqueness argument. The PARSEC-based effective-wavelength curvature correction (Eqs. 1–2) is a model prior on filter response, not a circular reuse of the target α. No uniqueness theorem, ansatz smuggling, or renaming of a known result is load-bearing. Residual selection purity and small core N are correctness/systematics concerns, not circularity.

Axiom & Free-Parameter Ledger

6 free parameters · 4 axioms · 0 invented entities

The central claim rests on standard photometric assumptions plus a handful of empirical selection choices and one model-dependent curvature correction. No new physical entities are postulated; free parameters are the fitted CER slopes themselves and the hand-chosen sample boundaries that define which stars enter those fits.

free parameters (6)
  • E(J−H)/E(Y−H) full-sample slope = 0.422 ± 0.004
    Fitted linear slope that is converted into α; the numerical value of α is completely determined by this fit once wavelengths are fixed.
  • E(VIS−H)/E(Y−H) full-sample slope = 2.855 ± 0.028
    Fitted slope used with AY/AH to obtain AVIS/AH.
  • Branch-selection line (J−H) < 0.42(Y−H)+0.01 = 0.42, +0.01
    Empirical cut that discards the upper branch; changes the sample from 20 456 to 14 210 sources and therefore affects all subsequent slopes.
  • J−H > 0.4 mag foreground cut = 0.4 mag
    Hand-chosen threshold to suppress low-extinction/foreground stars.
  • Mean (Y−H) thresholds for high/low subregions = >2.35 / <1.3 mag
    High cores defined by mean (Y−H)>2.35, low references by mean (Y−H)<1.3; these cuts define the samples that produce the claimed 17–27% α difference.
  • PARSEC red-giant parameters for λ eff = log g=2.0, log Teff=3.64, [M/H]=-0.1
    log g=2.0, log Teff=3.64, [M/H]=-0.1 chosen to compute unreddened effective wavelengths used in the curvature correction and in Eq. 4.
axioms (4)
  • domain assumption NIR extinction follows a pure power law A_λ ∝ λ^{-α} between the Y, J and H bands.
    Used in Eq. 4 to convert the measured CER into α; standard in the field but not re-derived here.
  • domain assumption After the model-based effective-wavelength correction, the slope of the observed color-color diagram equals the true color-excess ratio (intrinsic-color scatter only increases dispersion, not slope).
    Stated in §3.1; underpins every CER measurement.
  • ad hoc to paper The lower branch in the (J−H)–(Y−H) diagram is the relatively metal-rich giant population appropriate for Orion A.
    Justified by comparison with PARSEC tracks but remains an empirical choice that defines the working sample.
  • domain assumption Filter-averaged extinction can be corrected by convolving a single representative red-giant spectrum with the Euclid transmission curves.
    Eqs. 1–2 and the curvature-correction step in §3.1.

reviewed 2026-07-11 · how reviews work

0 comments
Cite this review

Pith. "Pith review of Euclid Q1 reveals spatial variations of the extinction law in the dense cloud LDN 1641." pith.science (2026). https://pith.science/paper/JCI56JRU

@misc{pith2026260704286,
  author       = {Pith},
  title        = {Pith review of: Euclid Q1 reveals spatial variations of the extinction law in the dense cloud LDN 1641},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/JCI56JRU}},
  note         = {Machine review of arXiv:2607.04286}
}
Share X Bluesky LinkedIn Reddit HN
read the original abstract

Dust extinction laws are essential for precision photometry and provide a direct probe of grain properties, but their behaviour in dense molecular clouds remains poorly constrained at high extinction. Using Euclid Quick Data Release 1 (Q1) imaging of the Orion A dark cloud Lynds Dark Nebula 1641 (LDN 1641), we measured the extinction law from the broad Visible Instrument (VIS) band and the Near-Infrared Spectrometer and Photometer (NISP) $Y$, $J$, and $H$ bands along sightlines reaching $A_V\sim 30$ mag towards the cloud core. We derived colour-excess ratios $E(\lambda-H)/E(Y-H)$ from linear fits to colour--colour diagrams of $(\lambda-H)$ versus $(Y-H)$ and converted them into relative extinctions, $A_\lambda/A_H$. The near-infrared extinction in LDN 1641 is well described by a power law, $A_\lambda\propto \lambda^{-\alpha}$, with $\alpha=1.57 \pm 0.06$, corresponding to $A_{\rm VIS}/A_H=4.23 \pm 0.24$, $A_Y/A_H=2.13 \pm 0.18$, and $A_J/A_H=1.47 \pm 0.11$. We further find significant spatial variations: $\alpha$ changes by up to $27\%$, with systematically smaller values and therefore flatter extinction curves in higher-extinction regions. This flattening is consistent with an enhanced large-grain population and supports substantial grain growth from the diffuse outskirts to the dense core of a single molecular cloud.

Figures

Figures reproduced from arXiv: 2607.04286 by Kun Wang, Rui Chen, Shu Wang, Xiaodian Chen.

Figure 1
Figure 1. Figure 1: Euclid (J − H) versus (Y − H) color–color diagrams for sources in the LDN 1641 field. Left: Scatter plot for N = 20456 sources after photometric-quality cuts, showing the bifurcated source distribution. For guidance, we overplot illustrative giant and dwarf evolutionary tracks for [M/H]= −0.1 and −1.5, reddened assuming a power-law extinction curve with α = 1.57 and AV = 7 mag. The black arrow shows an ill… view at source ↗
Figure 2
Figure 2. Figure 2: Spatial distributions of Herschel dust emission and the Euclid NIR binned mean color (Y − H) in LDN 1641. Left: Herschel/SPIRE 500 µm surface brightness map (I500; MJy sr−1 ), used as an empirical proxy sensitive to dust column density and temperature. Middle: Euclid pseudo￾color composite image of the region. Right: Interpolated map of the 20 × 20 binned mean (Y − H) color, serving as a reddening/extincti… view at source ↗
Figure 3
Figure 3. Figure 3: Optical and NIR color–color diagrams and extinction-law fits for the LDN 1641 region. Upper panels: Fits to the observed data without curvature correction. The left and right panels show (VIS − H) versus (Y − H) and (J − H) versus (Y − H), respectively. The black dashed, red dotted, and green solid lines represent the best-fitting linear, second-order, and third-order polynomials, respectively, illustratin… view at source ↗
Figure 4
Figure 4. Figure 4: Color excess E(J − H) versus E(Y − H) diagrams for subregions with different extinction levels in LDN 1641. Top row: Individual high￾extinction core subregions 1 (left) and 2 (right), selected with (Y − H) > 2.35 mag. Bottom row: Combined high-extinction core subregions 1 and 2 (left) and the combined low-extinction reference subregions 1–5 (right, mean (Y − H) < 1.3 mag). In each panel, red circles repres… view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

39 extracted references · 1 linked inside Pith

  1. [1]

    2017, ApJ, 849, L13

    Alonso-García, J., Minniti, D., Catelan, M., et al. 2017, ApJ, 849, L13

    2017

  2. [2]

    2012, MNRAS, 427, 127

    Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127

    2012

  3. [3]

    2024, AJ, 168, 256

    Cao, Z., Jiang, B., Wang, S., & Li, J. 2024, AJ, 168, 256

    2024

  4. [4]

    A., Clayton, G

    Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245

    1989

  5. [5]

    2025, Euclid Quick Data Release (Q1) – Data release overview

    Collaboration, E., Aussel, H., Tereno, I., et al. 2025, Euclid Quick Data Release (Q1) – Data release overview

    2025

  6. [6]

    Draine, B. T. 2003, ARA&A, 41, 241 Euclid Collaboration, Aussel, H., Tereno, I., et al. 2025, arXiv e-prints, arXiv:2503.15302

    Pith/arXiv arXiv 2003

  7. [7]

    G., Hora, J

    Fazio, G. G., Hora, J. L., Allen, L. E., et al. 2004, ApJS, 154, 10

    2004

  8. [8]

    Fitzpatrick, E. L. 1999, PASP, 111, 63

    1999

  9. [9]

    M., Pipher, J

    Flaherty, K. M., Pipher, J. L., Megeath, S. T., et al. 2007, ApJ, 663, 1069

    2007

  10. [10]

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

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

    2011

  11. [11]

    W., & Li, A

    Gao, J., Jiang, B. W., & Li, A. 2009, ApJ, 707, 89

    2009

  12. [12]

    J., Bandyopadhyay, R

    Gosling, A. J., Bandyopadhyay, R. M., & Blundell, K. M. 2009, MNRAS, 394, 2247

    2009

  13. [13]

    S., Babler, B

    Indebetouw, R., Mathis, J. S., Babler, B. L., et al. 2005, ApJ, 619, 931

    2005

  14. [14]

    J., et al

    Juvela, M., Ristorcelli, I., Marshall, D. J., et al. 2015, A&A, 584, A93

    2015

  15. [15]

    1983, A&A, 128, 84

    Koornneef, J. 1983, A&A, 128, 84

    1983

  16. [16]

    Larson, K. A. & Whittet, D. C. B. 2005, ApJ, 623, 897

    2005

  17. [17]

    2024, The Flattest Infrared Extinction Curve in Four Isolated Dense Molecular Cloud Cores

    Li, J., Chen, B., Jiang, B., et al. 2024, The Flattest Infrared Extinction Curve in Four Isolated Dense Molecular Cloud Cores

    2024

  18. [18]

    2023, ApJ, 956, 26

    Li, L., Wang, S., Chen, X., & Jiang, Q. 2023, ApJ, 956, 26

    2023

  19. [19]

    Lombardi, M., Alves, J., & Lada, C. J. 2011, A&A, 527, A60 Maíz Apellániz, J., Pantaleoni González, M., Barbá, R. H., García-Lario, P., &

    2011

  20. [20]

    2020, MNRAS, 496, 4951

    Nogueras-Lara, F. 2020, MNRAS, 496, 4951

    2020

  21. [21]

    Martin, P. G. & Whittet, D. C. B. 1990, ApJ, 357, 113

    1990

  22. [22]

    L., & Gordon, K

    Massa, D., Fitzpatrick, E. L., & Gordon, K. D. 2020, ApJ, 891, 67

    2020

  23. [23]

    Mathis, J. S. 1990, ARA&A, 28, 37

    1990

  24. [24]

    2018, A&A, 614, A65

    Meingast, S., Alves, J., & Lombardi, M. 2018, A&A, 614, A65

    2018

  25. [25]

    J., Menten, K

    Messineo, M., Habing, H. J., Menten, K. M., et al. 2005, A&A, 435, 575

    2005

  26. [26]

    2006, ApJ, 638, 839

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

    2006

  27. [27]

    2019, A&A, 630, L3

    Nogueras-Lara, F., Schödel, R., Najarro, F., et al. 2019, A&A, 630, L3

    2019

  28. [28]

    W., Min, M., Tielens, A

    Ormel, C. W., Min, M., Tielens, A. G. G. M., Dominik, C., & Paszun, D. 2011, A&A, 532, A43

    2011

  29. [29]

    & Henning, T

    Ossenkopf, V . & Henning, T. 1994, A&A, 291, 943

    1994

  30. [30]

    Rieke, G. H. & Lebofsky, M. J. 1985, ApJ, 288, 618 Schödel, R., Najarro, F., Muzic, K., & Eckart, A. 2010, A&A, 511, A18

    1985

  31. [31]

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

    2009

  32. [32]

    2021, ApJ, 915, 74

    Uehara, H., Dobashi, K., Nishiura, S., Shimoikura, T., & Naoi, T. 2021, ApJ, 915, 74

    2021

  33. [33]

    & Chen, X

    Wang, S. & Chen, X. 2019, ApJ, 877, 116

    2019

  34. [34]

    & Chen, X

    Wang, S. & Chen, X. 2024, ApJ, 964, L3

    2024

  35. [35]

    W., Li, A., & Chen, Y

    Wang, S., Gao, J., Jiang, B. W., Li, A., & Chen, Y . 2013, ApJ, 773, 30

    2013

  36. [36]

    & Jiang, B

    Wang, S. & Jiang, B. W. 2014, ApJ, 788, L12

    2014

  37. [37]

    W., Roellig, T

    Werner, M. W., Roellig, T. L., Low, F. J., et al. 2004, ApJS, 154, 1

    2004

  38. [38]

    L., Eisenhardt, P

    Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868

    2010

  39. [39]

    R., Indebetouw, R., et al

    Zasowski, G., Majewski, S. R., Indebetouw, R., et al. 2009, ApJ, 707, 510 Article number, page 7

    2009

This paper was first reviewed by grok-4.5 on July 11, 2026.