Pith. sign in

REVIEW 3 major objections 4 minor 41 references

The Optical Extinction Law Depends on Magnetic Field Orientation: The $R_V$-$\psi$ Relation

T0 review · 3 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash

Pith's one-line read The extinction law's reddening parameter $R_V$ varies from 3.21 to 3.05 depending on the angle between the magnetic field and the line of sight.

desk verdict A model-independent R_V–ψ relation with a real zero-point uncertainty: the assumed <sin²ψ>=2/3 for the composite extinction curve is likely wrong, and the abstract overstates the comparison to observed scatter. read the letter →

arxiv 2411.16839 v1 pith:7JNL64IW submitted 2024-11-25 astro-ph.GA

classification astro-ph.GA
keywords interstellardustextinctionlawR_VparametergrainalignmentpolarizedSerkowskimagneticfieldorientationGalacticmedium
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 establishes that the Milky Way's optical extinction law is not a fixed property of dust but depends on the orientation of the interstellar magnetic field relative to the line of sight. Because aspherical aligned grains extinguish more strongly when the field points toward the observer, and because aligned and unaligned grains have different wavelength dependence, the ratio $R_V \equiv A_V/E(B-V)$ must shift with the field angle $\psi$. Using only observationally determined total and polarized extinction curves of the diffuse interstellar medium, the paper derives a model-independent relation predicting $R_V = 3.21$ at $\psi = 0$ and $R_V = 3.05$ at $\psi = 90^\circ$. This 0.16 swing matches the observed scatter $\sigma(R_V)\simeq 0.2$ in stellar-spectroscopy surveys, suggesting much reported variation in the extinction law is geometric rather than compositional.

What carries the argument

Equation (9) is the load-bearing identity: $\tau(\lambda,\psi)=\tau_{\rm ran}(\lambda)+p(\lambda,90^\circ)(2/3-\sin^2\psi)$. It follows from the modified picket fence approximation and the split of grains into an aligned fraction $f$ and a randomly oriented fraction $1-f$, converting magnetic orientation into a smooth modulation of the extinction curve. From it, Equation (10) yields $R_V(\psi)$ directly, so the relation is fully prescribed by measured quantities and does not require a particular dust composition model.

What would settle it

Measure $R_V$ toward many individual stars or small cloud cores where the magnetic field angle $\psi$ is known independently from starlight polarization angles or background dust polarization. If $R_V$ does not follow the predicted roughly 0.16 increase from $\psi=90^\circ$ to $\psi=0$ on such sightlines, the relation fails in the regime where it should be cleanest.

Watch

Extended reading notes

Core claim

The central claim is that for partially aligned aspherical grains the total extinction toward a sightline can always be written as $\tau(\lambda,\psi)=\tau_{\rm ran}(\lambda)+p(\lambda,90^\circ)(2/3-\sin^2\psi)$, where $\tau_{\rm ran}$ is the orientation-averaged extinction and $p(\lambda,90^\circ)$ is the polarized extinction for a field in the plane of the sky. Since the observed polarized extinction has a different wavelength shape than the total extinction, Equation (10) follows: $R_V$ becomes a quotient of such terms and must decrease as the magnetic field moves from the line of sight into the plane of the sky. With $R^{\rm ran}_V=3.1$, $p_B/p_V=0.96$ from the Serkowski law, and $(p_V/\tau_V)_{\rm max}\simeq 0.046$, the paper obtains $R_V(0)=3.21$ and $R_V(90^\circ)=3.05$. The same geometric term predicts that the polarization fraction of dust emission and the polarization angle dispersion both correlate with $R_V$.

Load-bearing premise

The argument stands on the assumption that the composite Milky Way extinction curve used to set $R^{\rm ran}_V=3.1$ is truly the orientation-averaged curve, meaning its sightlines sample $\sin^2\psi$ with mean $2/3$; if they are biased toward a different field orientation, the quoted 3.21 and 3.05 shift.

Editorial extensions

If this is right

  • If the relation is right, $R_V$ varies by 0.16 across the sky, so a major part of the observed $\sigma(R_V)\simeq 0.18$ can be attributed to magnetic field geometry rather than to real changes in dust properties.
  • High-latitude sightlines with a single dominant cloud should show the cleanest signal: $R_V$ should track magnetic field orientation and correlate with dust polarization fraction.
  • Because the effect is a lower limit (if $(p_V/E(B-V))_{\rm max}$ exceeds 13% mag$^{-1}$, the swing grows), the predicted $R_V$ range may widen as polarization-to-reddening measurements improve.
  • The relation predicts $R_V \propto S^{-1/2}$ at high latitudes, where $S$ is the polarization angle dispersion, giving a direct observational test from existing polarization maps.

Reading between the lines

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

  • A natural extension is to treat magnetic field orientation as a nuisance parameter in 3D extinction and reddening maps, since the relation converts $\psi$ into an $R_V$ correction with no new dust physics.
  • In galaxies with a coherent magnetic field, the extinction law should vary systematically across the disk with viewing geometry; current extragalactic extinction corrections rarely include this term.
  • Combining the predicted linear $R_V$ versus $\sin^2\psi$ relation with all-sky polarization angle dispersion maps could produce a predicted $R_V$ sky that stellar spectroscopy can then test.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 4 minor

Summary. The paper derives a model-independent relation between the optical extinction parameter R_V and the angle ψ between the line of sight and the interstellar magnetic field. Using the modified picket fence approximation and assuming a population of perfectly aligned grains plus a randomly oriented population, the author writes the total extinction as τ(λ, ψ) = τ_ran(λ) + p(λ, 90°) (2/3 − sin²ψ) (Eq. 9). With observed total and polarized extinction curves, this yields Eq. (10), an algebraic relation R_V(ψ). Inserting R_ran_V = 3.1, p_B/p_V = 0.96, and (p_V/τ_V)_max = 0.046, the paper finds R_V(0°) = 3.21 and R_V(90°) = 3.05, a difference of 0.16. It argues this can explain much of the observed sky variation in R_V, especially at high Galactic latitudes, and predicts correlations with polarized dust emission diagnostics.

Significance. If the result holds, it connects two widely studied observables—the extinction law and dust polarization—and has practical implications for 3D dust mapping and extinction corrections. The derivation is transparent and algebraic: the R_V–ψ relation is not fitted to the target R_V variation but follows from independent measurements of total and polarized extinction curves, and it yields falsifiable predictions (e.g., R_V ∝ S^{-1/2} at high latitude). The qualitative conclusion that R_V depends on magnetic field orientation is physically well motivated and likely robust. However, the quantitative calibration depends on an untested assumption about the orientation-averaged nature of the composite extinction curve, and the abstract's comparison with observed σ(R_V) uses a mismatched statistic. These issues do not invalidate the core relation but require revision of the stated precision and of the abstract's claims.

major comments (3)
  1. [Section 3, Eq. (10)] The identification of the composite Milky Way extinction curve (R_ran_V = 3.1, Fitzpatrick et al. 2019) with the orientation-averaged curve τ_ran(λ) in Eq. (9) requires that the contributing sightlines have ⟨sin²ψ⟩ = 2/3. The composite is built primarily from low-latitude UV-bright OB stars; for a Galactic magnetic field lying mostly in the plane, the expected average is ⟨sin²ψ⟩ ≈ 1/2, not 2/3, so the composite contains a polarized contamination term p(λ,90°)(2/3 − ⟨sin²ψ⟩) that biases the zero-point of the R_V–ψ relation. A concrete estimate: with ⟨sin²ψ⟩ = 1/2, the effective R_ran_V entering Eq. (10) is ≈3.08 rather than 3.1, shifting R_V(0°) from 3.21 to ≈3.18 and R_V(90°) from 3.05 to ≈3.02. Because this systematic is neither tested nor propagated, the three-significant-figure values in the abstract are unsupported, even though the qualitative conclusion that R_V varies with ψ by ~0.16 is robust.
  2. [Section 3, abstract] The comparison between the predicted effect and observed scatter uses mismatched statistics. Equation (10) with a uniform distribution of cos ψ yields σ(R_V) = 0.048, as the paper states, whereas the observed σ(R_V) = 0.18 from Schlafly et al. (2016). The paper instead compares the total range 0.16 to σ = 0.18 and concludes in the abstract that the effect 'could therefore account for much of the large-scale R_V variation.' A range-to-standard-deviation comparison is not appropriate: the predicted σ is roughly one quarter of the observed σ, corresponding to only ~7% of the variance. Please either present a statistically consistent comparison or temper the abstract's claim.
  3. [Section 3, Eq. (10) and Figure 1] The headline values R_V(0°) = 3.21 and R_V(90°) = 3.05 are quoted without propagated uncertainties, even though the inputs (p_V/E(B−V))_max, p_B/p_V, and R_ran_V all carry substantial uncertainties. For (p_V/E(B−V))_max, Panopoulou et al. (2019) span 0.13–0.182 mag⁻¹, and Figure 2 shows a corresponding shaded band for the extinction curves, but no corresponding uncertainty is propagated to the R_V(ψ) values. Please provide an error estimate for R_V(ψ) or explicitly state that the quoted values are nominal and subject to the systematics discussed above.
minor comments (4)
  1. [Title and Abstract] The title has a typo: 'TheRV' should be 'The R_V'. In the abstract, 'induces' should be 'induce' (subject–verb agreement). In the Introduction, 'aspehrical' should be 'aspherical'.
  2. [Section 2, Eq. (2)] The expression for C_align_ext is algebraically correct but the notation is initially confusing because of the nested 1/2 factors. A brief parenthetical explaining that the second term is (1/2)C1 plus (1/4)(C2+C3) would improve readability.
  3. [Section 4, Eq. (16)] The statement that the observed polarization fraction is 'approximately linearly correlated with R_V' is ambiguous about sign. Since R_V decreases with sin²ψ while \tilde{p}/(1+\tilde{p}) increases, the correlation is negative; please state this explicitly.
  4. [Figures 1 and 2] Figure 1 shows R_V(ψ) without any uncertainty band. Since Figure 2 already displays a shaded band for the (p_V/E(B−V))_max range, propagating that range through Eq. (10) would give a useful visual representation of the calibration uncertainty.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: the R_V-psi relation is derived algebraically from observed total and polarized extinction curves rather than fitted to the target R_V variation.

full rationale

Equation (10) is obtained by substituting Equations (5) and (7) into the definition R_V = tau_V/(tau_B - tau_V); the quoted endpoints R_V(0)=3.21 and R_V(90)=3.05 are outputs of that substitution, not quantities used to set any parameter. The inputs are observationally determined: R_ran_V=3.1 (Fitzpatrick et al. 2019), a Serkowski-law polarized extinction curve (Whittet 2003), and (p_V/E(B-V))_max=0.13 mag^-1 (Panopoulou et al. 2019; Planck Collaboration XII 2020). The author's prior work is invoked for the MPFA (Draine & Hensley 2021), for compiled tau_ran and p curves (Hensley & Draine 2021), and for the sin^2 psi proportional to S^-1/2 relation (Hensley et al. 2019), but none of these citations already contains the R_V-psi relation or the specific functional form of Equation (10). The MPFA is an approximation validated by numerical tests, not a result equivalent to the paper's conclusion. The main caveat, flagged in Section 3, is the identification of the composite Milky Way extinction curve with tau_ran(lambda), which requires <sin^2 psi>=2/3 over the contributing sightlines; this is an empirical assumption that affects the zero-point of the relation, but it is not circular because R_ran_V is not defined in terms of the predicted R_V variation. Overall, the central derivation is self-contained, and the self-citations are not load-bearing in a circular sense.

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

The derivation itself introduces no free parameters or invented entities. The result rests on empirical inputs (mean R_V, maximum polarization efficiency, Serkowski parameters) and on the MPFA alignment model inherited from prior work, the main assumptions of which are listed.

free parameters (3)
  • R_ran_V = 3.1
    Mean Milky Way R_V from Fitzpatrick et al. (2019), adopted as the orientation-averaged value; an observational input, not tuned to the target.
  • (p_V/E(B-V))_max = 0.13 mag^-1
    Maximum observed polarization per reddening (Panopoulou et al. 2019; Planck Collaboration XII 2020), converted to (p_V/tau_V)_max = 0.046. The paper notes this is likely a lower limit and explores up to 0.182 in Figure 2.
  • p_B/p_V = 0.96
    Ratio from the Serkowski law with K = 0.87 and lambda_max = 5500 Angstrom (Whittet 2003); an observational approximation.
assumptions (5)
  • domain assumption Modified picket fence approximation (MPFA) accurately describes extinction cross sections of arbitrarily oriented grains.
    Invoked in Section 2; the authors state accuracy to about 10% for optical extinction and polarization (Draine & Hensley 2021).
  • domain assumption Perfect internal alignment, with grains rotating about their short axis, and two-state alignment: a fraction f are perfectly aligned and the rest are randomly oriented.
    Section 2 approximations 1 and 3; deviations would change the coefficient 2/3 in Equation (9).
  • domain assumption Magnetic dipole absorption is negligible, so extinction depends only on which principal axis the electric field is along.
    Section 2 approximation 2.
  • domain assumption The composite extinction curve is the orientation-averaged curve, i.e., the average of sin^2 psi over the sample is 2/3.
    Section 3, before Equation (10); if the sightline sample is biased toward plane-of-sky fields, R_ran_V is not the true orientation-averaged value.
  • domain assumption For polarized emission, all aligned grains have a single temperature and size-independent mass density, and there is one emitting region along the line of sight.
    Section 4, used only for the testable predictions involving polarization fraction and S.

how reviews work

0 comments
Cite this review

Pith. "Pith review of The Optical Extinction Law Depends on Magnetic Field Orientation: The $R_V$-$\psi$ Relation." pith.science (2026). https://pith.science/paper/7JNL64IW

@misc{pith2026241116839,
  author       = {Pith},
  title        = {Pith review of: The Optical Extinction Law Depends on Magnetic Field Orientation: The $R_V$-$\psi$ Relation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/7JNL64IW}},
  note         = {Machine review of arXiv:2411.16839}
}
abstract

For aspherical interstellar dust grains aligned with their short axes preferentially parallel to the local magnetic field, the amount of extinction per grain is larger when the magnetic field is along the line of sight and smaller when in the plane of the sky. To the extent that optical extinction arises from both aligned and unaligned grain populations with different extinction properties, changes in the magnetic field orientation induces changes in its wavelength dependence, parameterized by $R_V \equiv A_V/E(B-V)$. We demonstrate that the measured total and polarized extinction curves of the diffuse Galactic interstellar medium imply $R_V$ varies from 3.21 when the magnetic field is along the line of sight ($\psi = 0$) to $R_V = 3.05$ when in the plane of the sky ($\psi = 90^\circ$). This effect could therefore account for much of the large-scale $R_V$ variation observed across the sky ($\sigma(R_V) \simeq 0.2$), particularly at high Galactic latitudes.

Figures

Figures reproduced from arXiv: 2411.16839 by the authors.

Figure 1
Figure 1. Variation of RV with the angle ψ between the interstellar magnetic field and the line of sight. The relation is fully prescribed by observationally determined quantities as per Equation (10). 0.1 0.2 0.5 1.0 2.0 4.0 λ [µm] 0.90 0.95 1.00 1.05 1.10 τ/τ ran ψ = 0 ψ = 90◦ [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Extinction curve for ψ = 0 (dotted) and ψ = 90◦ (dashed) relative to the orientation-averaged ex￾tinction curve (sin2 ψ = 2/3, solid). The curves are derived from Equation (9) using τ ran (λ) and p (λ, 90◦ ) of the dif￾fuse Galactic ISM as compiled by Hensley & Draine (2021). The shaded region corresponds to (pV /E(B − V ))max be￾tween 0.13 and 0.182 mag−1 . Given that the polarized extinction p (λ) is observed to h… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

41 extracted references · 9 canonical work pages

  1. [1]

    A., Clayton, G

    Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245, doi: 10.1086/167900

  2. [2]

    Clark, S. E. 2018, ApJL, 857, L10, doi: 10.3847/2041-8213/aabb54

  3. [3]

    C., Anderson, C

    Clayton, G. C., Anderson, C. M., Magalhaes, A. M., et al. 1992, ApJL, 385, L53, doi: 10.1086/186276

  4. [4]

    Z., & Mitrofanov, I

    Dolginov, A. Z., & Mitrofanov, I. G. 1976, Ap&SS, 43, 291, doi: 10.1007/BF00640010

  5. [5]

    T., & Hensley, B

    Draine, B. T., & Hensley, B. S. 2021, ApJ, 919, 65, doi: 10.3847/1538-4357/ac0050

  6. [6]

    T., & Weingartner, J

    Draine, B. T., & Weingartner, J. C. 1996, ApJ, 470, 551, doi: 10.1086/177887 —. 1997, ApJ, 480, 633, doi: 10.1086/304008

  7. [7]

    M., Ade, P

    Fissel, L. M., Ade, P. A. R., Angil` e, F. E., et al. 2016, ApJ, 824, 134, doi: 10.3847/0004-637X/824/2/134

  8. [8]

    Clayton, G. C. 2019, ApJ, 886, 108, doi: 10.3847/1538-4357/ab4c3a

Show all 41 references
  1. [9]

    D., Cartledge, S., & Clayton, G

    Gordon, K. D., Cartledge, S., & Clayton, G. C. 2009, ApJ, 705, 1320, doi: 10.1088/0004-637X/705/2/1320

  2. [10]

    Greenberg, J. M. 1960, Journal of Applied Physics, 31, 82, doi: 10.1063/1.1735423 —. 1968, in Nebulae and Interstellar Matter, ed. B. M. Middlehurst & L. H. Aller, 221

  3. [11]

    M., & Chlewicki, G

    Greenberg, J. M., & Chlewicki, G. 1987, QJRAS, 28, 321

  4. [12]

    M., & Meltzer, A

    Greenberg, J. M., & Meltzer, A. S. 1960, ApJ, 132, 667, doi: 10.1086/146970

  5. [13]

    2018, A&A, 610, A16, doi: 10.1051/0004-6361/201630271

    Guillet, V., Fanciullo, L., Verstraete, L., et al. 2018, A&A, 610, A16, doi: 10.1051/0004-6361/201630271

  6. [14]

    Hall, J. S. 1949, Science, 109, 166, doi: 10.1126/science.109.2825.166

  7. [15]

    S., & Draine, B

    Hensley, B. S., & Draine, B. T. 2020, ApJ, 895, 38, doi: 10.3847/1538-4357/ab8cc3 —. 2021, ApJ, 906, 73, doi: 10.3847/1538-4357/abc8f1 —. 2023, ApJ, 948, 55, doi: 10.3847/1538-4357/acc4c2 The RV –ψ Relation 7

  8. [16]

    S., Zhang, C., & Bock, J

    Hensley, B. S., Zhang, C., & Bock, J. J. 2019, ApJ, 887, 159, doi: 10.3847/1538-4357/ab5183

  9. [17]

    Hiltner, W. A. 1949, Science, 109, 165, doi: 10.1126/science.109.2825.165

  10. [18]

    2016, ApJ, 831, 159, doi: 10.3847/0004-637X/831/2/159

    Hoang, T., & Lazarian, A. 2016, ApJ, 831, 159, doi: 10.3847/0004-637X/831/2/159

  11. [19]

    Hunter, J. D. 2007, Computing in Science and Engineering, 9, 90, doi: 10.1109/MCSE.2007.55

  12. [20]

    S., Babler, B

    Indebetouw, R., Mathis, J. S., Babler, B. L., et al. 2005, ApJ, 619, 931, doi: 10.1086/426679

  13. [21]

    2024, arXiv e-prints, arXiv:2408.15909, doi: 10.48550/arXiv.2408.15909

    Karim, T., Singh, S., Rezaie, M., et al. 2024, arXiv e-prints, arXiv:2408.15909, doi: 10.48550/arXiv.2408.15909

  14. [22]

    Kim, S.-H., & Martin, P. G. 1995, ApJ, 444, 293, doi: 10.1086/175604

  15. [23]

    M., & Draine, B

    Lee, H. M., & Draine, B. T. 1985, ApJ, 290, 211, doi: 10.1086/162974

  16. [24]

    Li, A., & Draine, B. T. 2001, ApJ, 554, 778, doi: 10.1086/323147

  17. [25]

    F., & Witt, A

    Lillie, C. F., & Witt, A. N. 1976, ApJ, 208, 64, doi: 10.1086/154582

  18. [26]

    G., Clayton, G

    Martin, P. G., Clayton, G. C., & Wolff, M. J. 1999, ApJ, 510, 905, doi: 10.1086/306613

  19. [27]

    G., Adamson, A

    Martin, P. G., Adamson, A. J., Whittet, D. C. B., et al. 1992, ApJ, 392, 691, doi: 10.1086/171470

  20. [28]

    V., Hensley, B

    Panopoulou, G. V., Hensley, B. S., Skalidis, R., Blinov, D., & Tassis, K. 2019, A&A, 624, L8, doi: 10.1051/0004-6361/201935266

  21. [29]

    V., & Lenz, D

    Panopoulou, G. V., & Lenz, D. 2020, ApJ, 902, 120, doi: 10.3847/1538-4357/abb6f5 Planck Collaboration X. 2016, A&A, 594, A10, doi: 10.1051/0004-6361/201525967 Planck Collaboration XII. 2020, A&A, 641, A12, doi: 10.1051/0004-6361/201833885 Planck Collaboration Int. XIX. 2015, A...

  22. [30]

    Purcell, E. M. 1975, in The Dusty Universe, ed. G. B. Field & A. G. W. Cameron, 155–167 —. 1979, ApJ, 231, 404, doi: 10.1086/157204

  23. [31]

    F., Meisner, A

    Schlafly, E. F., Meisner, A. M., Stutz, A. M., et al. 2016, ApJ, 821, 78, doi: 10.3847/0004-637X/821/2/78

  24. [32]

    J., Finkbeiner, D

    Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525, doi: 10.1086/305772

  25. [33]

    S., & Ford, V

    Serkowski, K., Mathewson, D. S., & Ford, V. L. 1975, ApJ, 196, 261, doi: 10.1086/153410

  26. [34]

    2023, A&A, 670, A115, doi: 10.1051/0004-6361/202243860

    Siebenmorgen, R. 2023, A&A, 670, A115, doi: 10.1051/0004-6361/202243860

  27. [35]

    Trumpler, R. J. 1930, Lick Observatory Bulletin, 420, 154, doi: 10.5479/ADS/bib/1930LicOB.14.154T van der Walt, S., Colbert, S. C., & Varoquaux, G. 2011, Computing in Science and Engineering, 13, 22, doi: 10.1109/MCSE.2011.37

  28. [36]

    Voshchinnikov, N. V. 1989, Astronomische Nachrichten, 310, 265, doi: 10.1002/asna.2113100405

  29. [37]

    C., & Draine, B

    Weingartner, J. C., & Draine, B. T. 2001, ApJ, 548, 296, doi: 10.1086/318651

  30. [38]

    Whittet, D. C. B., ed. 2003, Dust in the galactic environment

  31. [39]

    Anderson, C. M. 1997, ApJ, 478, 395, doi: 10.1086/303789

  32. [40]

    P., Guillet, V., et al

    Ysard, N., Jones, A. P., Guillet, V., et al. 2024, A&A, 684, A34, doi: 10.1051/0004-6361/202348391

  33. [41]

    2024, arXiv e-prints, arXiv:2407.14594, doi: 10.48550/arXiv.2407.14594

    Zhang, X., & Green, G. 2024, arXiv e-prints, arXiv:2407.14594, doi: 10.48550/arXiv.2407.14594

Pith tools

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