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 →
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
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.
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
- 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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)
- [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.
- [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.
- [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)
- [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'.
- [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.
- [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.
- [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
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
free parameters (3)
- R_ran_V =
3.1
- (p_V/E(B-V))_max =
0.13 mag^-1
- p_B/p_V =
0.96
assumptions (5)
- domain assumption Modified picket fence approximation (MPFA) accurately describes extinction cross sections of arbitrarily oriented grains.
- 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.
- domain assumption Magnetic dipole absorption is negligible, so extinction depends only on which principal axis the electric field is along.
- 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.
- 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.
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
Reference graph
Works this paper leans on
-
[1]
Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245, doi: 10.1086/167900
doi:10.1086/167900 1989
-
[2]
Clark, S. E. 2018, ApJL, 857, L10, doi: 10.3847/2041-8213/aabb54
-
[3]
Clayton, G. C., Anderson, C. M., Magalhaes, A. M., et al. 1992, ApJL, 385, L53, doi: 10.1086/186276
-
[4]
Dolginov, A. Z., & Mitrofanov, I. G. 1976, Ap&SS, 43, 291, doi: 10.1007/BF00640010
-
[5]
Draine, B. T., & Hensley, B. S. 2021, ApJ, 919, 65, doi: 10.3847/1538-4357/ac0050
-
[6]
Draine, B. T., & Weingartner, J. C. 1996, ApJ, 470, 551, doi: 10.1086/177887 —. 1997, ApJ, 480, 633, doi: 10.1086/304008
doi:10.1086/177887 1996
-
[7]
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]
Clayton, G. C. 2019, ApJ, 886, 108, doi: 10.3847/1538-4357/ab4c3a
Show all 41 references
-
[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
2009 doi
-
[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
1960 doi
-
[11]
M., & Chlewicki, G
Greenberg, J. M., & Chlewicki, G. 1987, QJRAS, 28, 321
1987
-
[12]
M., & Meltzer, A
Greenberg, J. M., & Meltzer, A. S. 1960, ApJ, 132, 667, doi: 10.1086/146970
1960 doi
-
[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
2018 doi
-
[14]
Hall, J. S. 1949, Science, 109, 166, doi: 10.1126/science.109.2825.166
1949 doi
-
[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
2020 doi
-
[16]
S., Zhang, C., & Bock, J
Hensley, B. S., Zhang, C., & Bock, J. J. 2019, ApJ, 887, 159, doi: 10.3847/1538-4357/ab5183
2019 doi
-
[17]
Hiltner, W. A. 1949, Science, 109, 165, doi: 10.1126/science.109.2825.165
1949 doi
-
[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
2016 doi
-
[19]
Hunter, J. D. 2007, Computing in Science and Engineering, 9, 90, doi: 10.1109/MCSE.2007.55
2007 doi
-
[20]
S., Babler, B
Indebetouw, R., Mathis, J. S., Babler, B. L., et al. 2005, ApJ, 619, 931, doi: 10.1086/426679
2005 doi
- [21]
-
[22]
Kim, S.-H., & Martin, P. G. 1995, ApJ, 444, 293, doi: 10.1086/175604
1995 doi
- [23]
-
[24]
Li, A., & Draine, B. T. 2001, ApJ, 554, 778, doi: 10.1086/323147
2001 doi
- [25]
-
[26]
G., Clayton, G
Martin, P. G., Clayton, G. C., & Wolff, M. J. 1999, ApJ, 510, 905, doi: 10.1086/306613
1999 doi
-
[27]
G., Adamson, A
Martin, P. G., Adamson, A. J., Whittet, D. C. B., et al. 1992, ApJ, 392, 691, doi: 10.1086/171470
1992 doi
-
[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
2019 doi
-
[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...
2020 doi
-
[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
1975 doi
-
[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
2016 doi
-
[32]
J., Finkbeiner, D
Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525, doi: 10.1086/305772
1998 doi
-
[33]
S., & Ford, V
Serkowski, K., Mathewson, D. S., & Ford, V. L. 1975, ApJ, 196, 261, doi: 10.1086/153410
1975 doi
-
[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
2023 doi
-
[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
1930 doi
-
[36]
Voshchinnikov, N. V. 1989, Astronomische Nachrichten, 310, 265, doi: 10.1002/asna.2113100405
1989 doi
-
[37]
C., & Draine, B
Weingartner, J. C., & Draine, B. T. 2001, ApJ, 548, 296, doi: 10.1086/318651
2001 doi
-
[38]
Whittet, D. C. B., ed. 2003, Dust in the galactic environment
2003
-
[39]
Anderson, C. M. 1997, ApJ, 478, 395, doi: 10.1086/303789
1997 doi
-
[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
2024 doi
- [41]
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.