Recognition: no theorem link
Millimeter dust continuum and polarization in protoplanetary disks with scattering: A slab model
Pith reviewed 2026-05-15 03:05 UTC · model grok-4.3
The pith
Common analytic approximations underestimate the millimeter continuum emission from protoplanetary disks by 10 to 15 percent.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We numerically solve the radiative transfer equation in an isothermal, constant-density plane-parallel slab including dust absorption, emission, and self-scattering with full Stokes parameters. Commonly used analytic approximations for the continuum emission are systematically about 10 to 15% lower than our numerical solutions. We provide empirical fitting formulae that reproduce our numerical results for the continuum emission and polarization fraction.
What carries the argument
Numerical solution of the radiative transfer equation with full Stokes parameters in an isothermal constant-density plane-parallel slab.
If this is right
- SED analyses using old approximations overestimate optical depth and thus disk mass.
- Old approximations lead to overestimated dust temperature estimates.
- Albedo is underestimated, altering constraints on grain size from polarization data.
- New empirical formulae enable more accurate and efficient analysis of (sub)millimeter observations.
Where Pith is reading between the lines
- The corrections may revise dust grain size estimates in many protoplanetary disk observations.
- More realistic vertical structures in disks could amplify or reduce the reported discrepancies.
- The new polarization formulae can be directly compared to resolved ALMA polarization maps.
Load-bearing premise
The isothermal, constant-density plane-parallel slab geometry adequately represents the emission and polarization properties of real protoplanetary disks at millimeter wavelengths.
What would settle it
Comparing the new fitting formulae to full 3D radiative transfer calculations in stratified disk models with varying density and temperature would show if the 10-15 percent difference holds.
Figures
read the original abstract
Millimeter continuum emission and self-scattering polarization from protoplanetary disks are widely used to constrain dust properties. Interpreting these observations requires practical prescriptions for the disk emission. However, only approximate formulae are available for the continuum emission, and no widely applicable formula has yet been established for the polarized emission. We aim (i) to assess the validity of commonly used analytic approximations for the (sub)millimeter continuum emission from protoplanetary disks, and (ii) to derive realistic prescriptions for the disk emission for both the continuum and the polarization. We numerically solve the radiative transfer equation in an isothermal, constant-density plane-parallel slab, including dust absorption, emission, and self-scattering with full Stokes parameters. We find that commonly used analytic approximations for the continuum emission are systematically about 10 to 15% lower than our numerical solutions. Consequently, SED analyses of (sub)millimeter observations that adopt these formulae are likely to overestimate the optical depth (and thus the disk mass) and the dust temperature, and underestimate the albedo (and thus altering the inferred constraints on grain size). We also provide empirical fitting formulae that reproduce our numerical results for the continuum emission and polarization fraction. These formulae will enable observational data analyses to be carried out more accurately and efficiently than with the conventional approaches. For the analysis of (sub)millimeter observations, we recommend using our new empirical formulae or interpolation of our numerical results, rather than commonly used approximations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript numerically solves the radiative transfer equation in an isothermal, constant-density plane-parallel slab including dust absorption, emission, and self-scattering with full Stokes parameters. It reports that commonly used analytic approximations for the millimeter continuum emission are systematically 10-15% lower than the numerical solutions, leading to likely overestimates of optical depth, disk mass, and dust temperature (and underestimates of albedo) in SED analyses; empirical fitting formulae are provided for both total intensity and polarization fraction as practical replacements.
Significance. If the slab results generalize, the quantified 10-15% offset would refine mass and grain-size constraints from (sub)mm observations of protoplanetary disks, and the new fitting formulae would offer a reproducible, efficient improvement over existing approximations. The direct numerical integration of the RTE and post-processed empirical matches to those integrations constitute a controlled, standard-method benchmark that strengthens the comparison.
major comments (2)
- [Abstract] Abstract: the recommendation that 'SED analyses ... are likely to overestimate the optical depth (and thus the disk mass)' is load-bearing for the paper's applied claim, yet rests entirely on the isothermal constant-density slab; no demonstration is given that the 10-15% offset persists (or even retains sign) once vertical temperature gradients and exponential density stratification are introduced.
- [Methods/Results] Methods/Results: the empirical fitting formulae are derived exclusively from the slab solutions; without a sensitivity test to realistic T(z) profiles or radial optical-depth gradients, the assertion that these formulae should replace conventional approaches in observational data analyses remains unvalidated for the target application.
minor comments (1)
- Tabulate the validity ranges (optical depth, albedo, wavelength) for each empirical fitting formula and state the maximum residual relative to the numerical solutions.
Simulated Author's Rebuttal
We thank the referee for the constructive comments and recommendation. Our study is explicitly limited to the isothermal constant-density slab as a controlled benchmark to quantify biases in common analytic approximations. We agree the claims require qualification for realistic stratified disks and will revise the abstract, add a limitations discussion, and soften recommendations accordingly.
read point-by-point responses
-
Referee: [Abstract] Abstract: the recommendation that 'SED analyses ... are likely to overestimate the optical depth (and thus the disk mass)' is load-bearing for the paper's applied claim, yet rests entirely on the isothermal constant-density slab; no demonstration is given that the 10-15% offset persists (or even retains sign) once vertical temperature gradients and exponential density stratification are introduced.
Authors: We acknowledge that the 10-15% offset and resulting bias direction are demonstrated only for the isothermal constant-density slab. Without additional calculations for stratified T(z) and density profiles, we cannot confirm persistence or sign of the offset. In revision we will update the abstract to state that the overestimation applies to slab-based SED analyses and add a discussion paragraph noting that the bias may vary with realistic vertical structure, recommending the formulae be used cautiously until further tests are available. revision: partial
-
Referee: [Methods/Results] Methods/Results: the empirical fitting formulae are derived exclusively from the slab solutions; without a sensitivity test to realistic T(z) profiles or radial optical-depth gradients, the assertion that these formulae should replace conventional approaches in observational data analyses remains unvalidated for the target application.
Authors: The fitting formulae and numerical solutions are derived exclusively from the slab, consistent with the paper's scope as a benchmark study. We agree sensitivity tests to T(z) and radial gradients are needed for full validation in complex disks. In the revision we will qualify the recommendation in the abstract and conclusions to specify that the formulae replace conventional approximations for slab-like models, provide the numerical grid for interpolation, and note that users should verify applicability for stratified geometries. revision: partial
- Demonstration that the 10-15% offset persists (or retains sign) with vertical temperature gradients and exponential density stratification
Circularity Check
Numerical RTE solutions independent; fits are post-hoc matches
full rationale
The paper's core results follow from direct numerical integration of the full-Stokes radiative transfer equation in an isothermal slab. The reported 10-15% systematic offset is measured against external analytic approximations, not against any quantity defined inside the paper. The new empirical fitting formulae are explicitly constructed as post-processing matches to those numerical outputs and do not redefine or presuppose the inputs. No self-citation chain, uniqueness theorem, or definitional loop appears in the derivation; the slab geometry is stated as an assumption rather than derived from the results themselves.
Axiom & Free-Parameter Ledger
free parameters (1)
- coefficients in the empirical fitting formulae
axioms (1)
- domain assumption Isothermal, constant-density plane-parallel slab
Reference graph
Works this paper leans on
-
[1]
Effects of Particle Size Distribution on Opacity Curves of Protoplanetary Disks around T Tauri Stars. , keywords =. doi:10.1006/icar.1993.1156 , adsurl =
-
[2]
Dust properties of protoplanetary disks in the Taurus-Auriga star forming region from millimeter wavelengths. , keywords =. doi:10.1051/0004-6361/200913403 , archivePrefix =. 0912.3356 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1051/0004-6361/200913403
-
[3]
The Radial Distribution of Dust Particles in the HL Tau Disk from ALMA and VLA Observations. , keywords =. doi:10.3847/1538-4357/ab3d33 , archivePrefix =. 1908.07140 , primaryClass =
-
[4]
Characterizing the dust content of disk substructures in TW Hydrae. , keywords =. doi:10.1051/0004-6361/202039812 , archivePrefix =. 2102.04648 , primaryClass =
-
[5]
Scattering-induced Intensity Reduction: Large Mass Content with Small Grains in the Inner Region of the TW Hya disk. , keywords =. doi:10.3847/1538-4357/ab8223 , archivePrefix =. 2003.09353 , primaryClass =
-
[6]
Impact of Differential Dust Settling on the SED and Polarization: Application to the Inner Region of the HL Tau Disk. , keywords =. doi:10.3847/1538-4357/abf7b8 , archivePrefix =. 2104.05927 , primaryClass =
-
[7]
Massive Compact Dust Disk with a Gap around CW Tau Revealed by ALMA Multiband Observations. , keywords =. doi:10.3847/1538-4357/ac634d , archivePrefix =. 2203.16236 , primaryClass =
-
[8]
Molecules with ALMA at Planet-forming Scales (MAPS). XIV. Revealing Disk Substructures in Multiwavelength Continuum Emission. , keywords =. doi:10.3847/1538-4365/ac1431 , archivePrefix =. 2109.06433 , primaryClass =
-
[9]
Into the thick of it: ALMA 0.45 mm observations of HL Tau at a resolution of 2 au. , keywords =. doi:10.1051/0004-6361/202349046 , archivePrefix =. 2404.04164 , primaryClass =
-
[10]
Radiative processes in astrophysics
-
[11]
An analytical model of radial dust trapping in protoplanetary disks
An Analytical Model of Radial Dust Trapping in Protoplanetary Disks. , keywords =. doi:10.3847/1538-4357/ab1265 , archivePrefix =. 1903.08769 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-4357/ab1265 1903
-
[12]
One Solution to the Mass Budget Problem for Planet Formation: Optically Thick Disks with Dust Scattering. , keywords =. doi:10.3847/2041-8213/ab1f8c , archivePrefix =. 1904.02127 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/2041-8213/ab1f8c 2041
-
[13]
Opacity of fluffy dust aggregates
Opacity of fluffy dust aggregates. , keywords =. doi:10.1051/0004-6361/201323199 , archivePrefix =. 1312.1459 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1051/0004-6361/201323199
-
[14]
Millimeter-wave polarization of protoplanetary disks due to dust scattering
Millimeter-wave Polarization of Protoplanetary Disks due to Dust Scattering. , keywords =. doi:10.1088/0004-637X/809/1/78 , archivePrefix =. 1504.04812 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/0004-637x/809/1/78
-
[15]
Grain size constraints on HL Tau with polarization signature
Grain Size Constraints on HL Tau with Polarization Signature. , keywords =. doi:10.3847/0004-637X/820/1/54 , archivePrefix =. 1507.08902 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/0004-637x/820/1/54
-
[16]
The Evidence of Radio Polarization Induced by the Radiative Grain Alignment and Self-scattering of Dust Grains in a Protoplanetary Disk. , keywords =. doi:10.3847/2041-8213/aa7e33 , archivePrefix =. 1707.01612 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/2041-8213/aa7e33 2041
-
[17]
Inclination-induced polarization of scattered millimetre radiation from protoplanetary discs: the case of HL Tau. , keywords =. doi:10.1093/mnras/stv2633 , archivePrefix =. 1507.08353 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1093/mnras/stv2633
-
[18]
Scattering-produced (sub)millimetre polarization in inclined discs: optical depth effects, near-far side asymmetry and dust settling. , keywords =. doi:10.1093/mnras/stx1951 , archivePrefix =. 1705.05432 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1093/mnras/stx1951
-
[19]
Unveiling Dust Aggregate Structure in Protoplanetary Disks by Millimeter-wave Scattering Polarization. , keywords =. doi:10.3847/1538-4357/ab45f0 , archivePrefix =. 1907.00189 , primaryClass =
-
[20]
Porous Dust Particles in Protoplanetary Disks: Application to the HL Tau Disk. , keywords =. doi:10.3847/1538-4357/acdb4e , archivePrefix =. 2306.00158 , primaryClass =
-
[21]
Support for fragile porous dust in a gravitationally self-regulated disk around IM Lup. Nature Astronomy , keywords =. doi:10.1038/s41550-024-02308-6 , archivePrefix =. 2406.07427 , primaryClass =
-
[22]
The Disk Substructures at High Angular Resolution Project (DSHARP). V. Interpreting ALMA Maps of Protoplanetary Disks in Terms of a Dust Model. , keywords =. doi:10.3847/2041-8213/aaf743 , archivePrefix =. 1812.04043 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/2041-8213/aaf743 2041
- [23]
-
[24]
Thermal emission and scattering by aligned grains: Plane-parallel model and application to multiwavelength polarization of the HL Tau disc. , keywords =. doi:10.1093/mnras/stac753 , archivePrefix =. 2112.10998 , primaryClass =
-
[25]
, year = 1952, month = jan, volume =
The alignment of galactic dust. , year = 1952, month = jan, volume =. doi:10.1093/mnras/112.2.215 , adsurl =
-
[26]
Subsonic Mechanical Alignment of Irregular Grains
Subsonic Mechanical Alignment of Irregular Grains. , keywords =. doi:10.1086/523849 , archivePrefix =. 0707.3805 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1086/523849
-
[27]
Badminton birdie-like aerodynamic alignment of drifting dust grains by subsonic gaseous flows in protoplanetary discs. , keywords =. doi:10.1093/mnras/stae2248 , archivePrefix =. 2407.10025 , primaryClass =
-
[28]
Theory of microwave remote sensing by Tsang , year = 1985, pages =
Theory of microwave remote sensing. Theory of microwave remote sensing by Tsang , year = 1985, pages =
work page 1985
-
[29]
Absorption and scattering of light by small particles
-
[30]
Volume 1 - Single scattering and transport theory
Wave propagation and scattering in random media. Volume 1 - Single scattering and transport theory. doi:10.1016/B978-0-12-374701-3.X5001-7 , adsurl =
-
[31]
Multiwavelength Dust Characterization of the HL Tau Disk and Implications for Planet Formation. , keywords =. doi:10.3847/1538-4357/adf214 , archivePrefix =. 2507.14443 , primaryClass =
- [32]
-
[33]
Journal of Geophysical Research (Atmospheres) , keywords =
Optical constants of ice from the ultraviolet to the microwave: A revised compilation. Journal of Geophysical Research (Atmospheres) , keywords =. doi:10.1029/2007JD009744 , adsurl =
-
[34]
Disentangling the independently controllable factors of variation by interacting with the world
Interstellar Dust Grains. , keywords =. doi:10.1146/annurev.astro.41.011802.094840 , archivePrefix =. astro-ph/0304489 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1146/annurev.astro.41.011802.094840
-
[35]
Dust opacities for protoplanetary accretion disks: influence of dust aggregates. , keywords =
-
[36]
The size distribution of interstellar grains. , keywords =. doi:10.1086/155591 , adsurl =
-
[37]
Particle Emissivity in Circumstellar Disks. , keywords =. doi:10.1086/170646 , adsurl =
-
[38]
Protoplanetary Disks and Their Evolution
Protoplanetary Disks and Their Evolution. , keywords =. doi:10.1146/annurev-astro-081710-102548 , archivePrefix =. 1103.0556 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1146/annurev-astro-081710-102548
-
[39]
Observations of Protoplanetary Disk Structures. , keywords =. doi:10.1146/annurev-astro-031220-010302 , archivePrefix =. 2001.05007 , primaryClass =
-
[40]
Protoplanetary Disk Chemistry. , keywords =. doi:10.1146/annurev-astro-022823-040820 , archivePrefix =. 2309.05685 , primaryClass =
-
[41]
, year = 1983, month = sep, volume =
The determination of cloud masses and dust characteristics from submillimetre thermal emission. , year = 1983, month = sep, volume =
work page 1983
-
[42]
A Survey for Circumstellar Disks around Young Stellar Objects. , keywords =. doi:10.1086/115385 , adsurl =
-
[43]
Circumstellar Dust Disks in Taurus-Auriga: The Submillimeter Perspective
Circumstellar Dust Disks in Taurus-Auriga: The Submillimeter Perspective. , keywords =. doi:10.1086/432712 , archivePrefix =. astro-ph/0506187 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1086/432712
-
[44]
Progress of Theoretical Physics , year = 1976, month = dec, volume =
The gas drag effect on the elliptical motion of a solid body in the primordial solar nebula. Progress of Theoretical Physics , year = 1976, month = dec, volume =. doi:10.1143/PTP.56.1756 , adsurl =
-
[45]
Aerodynamics of solid bodies in the solar nebula. , keywords =. doi:10.1093/mnras/180.2.57 , adsurl =
-
[46]
Dust Coagulation. , keywords =. doi:10.1086/172562 , adsurl =
-
[47]
Coagulation, fragmentation and radial motion of solid particles in protoplanetary disks
Coagulation, fragmentation and radial motion of solid particles in protoplanetary disks. , keywords =. doi:10.1051/0004-6361:20077759 , archivePrefix =. 0711.2192 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1051/0004-6361:20077759
-
[48]
Panchromatic (Sub)millimeter polarization observations of HL Tau unveil aligned scattering grains. , keywords =. doi:10.1093/mnras/stae040 , archivePrefix =. 2309.10055 , primaryClass =
-
[49]
Scattering by Interstellar Dust Grains. I. Optical and Ultraviolet. , keywords =. doi:10.1086/379118 , archivePrefix =. astro-ph/0304060 , primaryClass =
work page internal anchor Pith review Pith/arXiv arXiv doi:10.1086/379118
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.