pith. machine review for the scientific record. sign in

arxiv: 2605.15014 · v1 · submitted 2026-05-14 · 🌌 astro-ph.EP

Recognition: no theorem link

Millimeter dust continuum and polarization in protoplanetary disks with scattering: A slab model

Authors on Pith no claims yet

Pith reviewed 2026-05-15 03:05 UTC · model grok-4.3

classification 🌌 astro-ph.EP
keywords protoplanetary disksmillimeter continuumself-scattering polarizationradiative transferdust propertiesempirical fitting formulaeslab modeloptical depth
0
0 comments X

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.

The paper tests how well simple analytic formulas work for calculating the millimeter light coming from protoplanetary disks, where dust both emits and scatters the light. By solving the full radiative transfer equation numerically in a simple slab geometry, the authors show that the usual approximations give intensities that are too low by 10 to 15 percent. This bias means that analyses of disk observations would overestimate the amount of dust and its temperature while underestimating how much light the grains reflect. The authors supply new fitting formulas that match their numerical results for both total emission and the degree of polarization, allowing more accurate interpretation of telescope data.

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

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

  • 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

Figures reproduced from arXiv: 2605.15014 by Akimasa Kataoka, Naoya Kitade.

Figure 1
Figure 1. Figure 1: Schematic representation of the plane-parallel slab and [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Schematic defining einc and eperp. eperp is defined as the direction perpendicular to both the z-axis and the line of sight, and einc as the direction perpendicular to both the line of sight and eperp. the mass extinction opacity. Throughout this paper, optical depth refers to the extinction optical depth. I, Q, U, V are the Stokes parameters, which have the same units as the specific intensity. We omit th… view at source ↗
Figure 4
Figure 4. Figure 4: Stokes I, Stokes Q, and polarization fraction (≡ Q/I) of the emergent intensity as functions of disk inclination for τmax = 10, obtained by numerically solving Eq. (1). Stokes I and Q are normalized by the Planck function, B. Here, the polar￾ization fraction is defined as 100× Q/I (in percent), without tak￾ing the absolute value. Each curve represents a different albedo value, ω. While the numerical result… view at source ↗
Figure 6
Figure 6. Figure 6: Peak emergent polarization fraction as a function of disk [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Stokes I as a function of z and the optical depth from the surface τ, obtained by numerically solving Eq. (1). The slab is inclined at 45°; results are shown for ω = 0.5 and 0.6. Stokes I is normalized by the Planck function, B. The lower axis shows ρκextz, while the upper axis shows τ. τ = 0 corresponds to the slab surface facing the observer. The “surface-layer effect,” which is the attenuation of Stokes… view at source ↗
Figure 8
Figure 8. Figure 8: Stokes I as a function of z and the optical depth from the surface τ, obtained by numerically solving Eq. (1). Results are shown i = 45° and 75°, with ω = 0.5. Stokes I is normalized by the Planck function, B. The lower axis shows ρκextz, while the upper axis shows τ. 10 3 10 2 10 1 10 0 10 1 ½ extz 10 3 10 2 10 1 I= B z-dependence of Stokes I (i = 45 ± ; ! = 0:5) ¿max = 0:1 ¿max = 1:0 ¿max = 10:0 [PITH_F… view at source ↗
Figure 9
Figure 9. Figure 9: Emergent Stokes I as a function of z, obtained by numer￾ically solving Eq. (1). The slab is inclined at 45° with ω = 0.5. Results are shown for τmax = 0.1, 1.0, and 10.0. Stokes I is nor￾malized by the Planck function, B. 3.2.2. Scattering of the polarized intensity Scattering of the incoming polarized component can make a non￾negligible contribution to the emergent polarization. This is no￾table because t… view at source ↗
Figure 11
Figure 11. Figure 11: Comparison of the scattering-angle dependence of the [PITH_FULL_IMAGE:figures/full_fig_p008_11.png] view at source ↗
Figure 10
Figure 10. Figure 10: Comparison of the emergent Stokes I and polarization fraction as functions of total optical depth (τmax) for a slab in￾clined at 45° with ω = 0.9, computed both including (blue curves, reproducing [PITH_FULL_IMAGE:figures/full_fig_p008_10.png] view at source ↗
Figure 12
Figure 12. Figure 12: Relative error in the emergent Stokes I between calcula￾tions neglecting and including the polarization of the incoming intensity, defined as (Iunpol − Ipol)/Ipol. By contrast, the polarized component of the incoming inten￾sity has a negligible impact on the emergent Stokes I. In the upper panel of [PITH_FULL_IMAGE:figures/full_fig_p008_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: The upper panels compare the approximate formulae for the emergent Stokes [PITH_FULL_IMAGE:figures/full_fig_p010_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Upper panel: Comparison between our numerical re [PITH_FULL_IMAGE:figures/full_fig_p011_14.png] view at source ↗
Figure 16
Figure 16. Figure 16: ω-dependence of τpeak for i = 45°, where τpeak is the total optical depth at which the polarization fraction of the emergent intensity peaks. τPF) for i = 15°, 30°, 45°, 60°, 75° are provided in Appendix F. These fitting parameters are available on the website1 . From Eq. (18), we can derive τpeak, the total optical depth at which the polarization fraction of the emergent intensity peaks. Specifically, τp… view at source ↗
Figure 17
Figure 17. Figure 17: Upper panels: Comparison, for the four dust models, between the numerically computed emergent Stokes [PITH_FULL_IMAGE:figures/full_fig_p013_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: Angular dependence of the scattering matrix elements [PITH_FULL_IMAGE:figures/full_fig_p014_18.png] view at source ↗
Figure 19
Figure 19. Figure 19: Comparison, for the four dust models, between the numerically computed emergent polarization including Mie scattering [PITH_FULL_IMAGE:figures/full_fig_p014_19.png] view at source ↗
Figure 20
Figure 20. Figure 20: Comparison of the emergent Stokes I and polarization fraction from a plane-parallel slab inclined at 45° between RADMC￾3D (Dullemond et al. 2012) simulations and our numerical solutions. Two cases are considered in the simulations: a Rayleigh scattering matrix with ω = 0.7, and a Mie scattering matrix corresponding to the case with amax f = 200 µm, f = 1, and a wavelength of 870 µm adopted in Section 6. T… view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

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)
  1. [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.
  2. [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)
  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

2 responses · 1 unresolved

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
  1. 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

  2. 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

standing simulated objections not resolved
  • Demonstration that the 10-15% offset persists (or retains sign) with vertical temperature gradients and exponential density stratification

Circularity Check

0 steps flagged

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

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the numerical RTE solution under a simplified slab geometry plus subsequent empirical fitting; no new physical entities are introduced.

free parameters (1)
  • coefficients in the empirical fitting formulae
    Parameters are tuned to reproduce the numerical RTE output for intensity and polarization.
axioms (1)
  • domain assumption Isothermal, constant-density plane-parallel slab
    Assumed to reduce the RTE to a one-dimensional numerical problem solvable without full 3D disk structure.

pith-pipeline@v0.9.0 · 5563 in / 1214 out tokens · 52121 ms · 2026-05-15T03:05:08.033363+00:00 · methodology

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

49 extracted references · 49 canonical work pages · 16 internal anchors

  1. [1]

    , keywords =

    Effects of Particle Size Distribution on Opacity Curves of Protoplanetary Disks around T Tauri Stars. , keywords =. doi:10.1006/icar.1993.1156 , adsurl =

  2. [2]

    Dust Properties of Protoplanetary Disks in the Taurus-Auriga Star Forming Region from Millimeter Wavelengths

    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 =

  3. [3]

    , keywords =

    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. [4]

    , keywords =

    Characterizing the dust content of disk substructures in TW Hydrae. , keywords =. doi:10.1051/0004-6361/202039812 , archivePrefix =. 2102.04648 , primaryClass =

  5. [5]

    , keywords =

    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. [6]

    , keywords =

    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. [7]

    , keywords =

    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. [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. [9]

    , keywords =

    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. [10]

    Radiative processes in astrophysics

  11. [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 =

  12. [12]

    One Solution to the Mass Budget Problem for Planet Formation: Optically Thick Disks with Dust Scattering

    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 =

  13. [13]

    Opacity of fluffy dust aggregates

    Opacity of fluffy dust aggregates. , keywords =. doi:10.1051/0004-6361/201323199 , archivePrefix =. 1312.1459 , primaryClass =

  14. [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 =

  15. [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 =

  16. [16]

    The evidence of radio polarization induced by the radiative grain alignment and self-scattering of dust grains in a protoplanetary disk

    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 =

  17. [17]

    Inclination-Induced Polarization of Scattered Millimeter Radiation from Protoplanetary Disks: The Case of HL Tau

    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 =

  18. [18]

    Scattering-Produced (Sub)millimeter Polarization in Inclined Disks: Optical Depth Effects, Near-Far Side Asymmetry, and Dust Settling

    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 =

  19. [19]

    , keywords =

    Unveiling Dust Aggregate Structure in Protoplanetary Disks by Millimeter-wave Scattering Polarization. , keywords =. doi:10.3847/1538-4357/ab45f0 , archivePrefix =. 1907.00189 , primaryClass =

  20. [20]

    , keywords =

    Porous Dust Particles in Protoplanetary Disks: Application to the HL Tau Disk. , keywords =. doi:10.3847/1538-4357/acdb4e , archivePrefix =. 2306.00158 , primaryClass =

  21. [21]

    Nature Astronomy , keywords =

    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. [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 =

  23. [23]

    2012 , note =

    RADMC-3D: A multi-purpose radiative transfer tool. 2012 , note =

  24. [24]

    , keywords =

    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. [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. [26]

    Subsonic Mechanical Alignment of Irregular Grains

    Subsonic Mechanical Alignment of Irregular Grains. , keywords =. doi:10.1086/523849 , archivePrefix =. 0707.3805 , primaryClass =

  27. [27]

    , keywords =

    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. [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 =

  29. [29]

    Absorption and scattering of light by small particles

  30. [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. [31]

    , keywords =

    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. [32]

    12.3, ver

    ALMA Technical Handbook, ALMA Doc. 12.3, ver. 1.5

  33. [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. [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 =

  35. [35]

    , keywords =

    Dust opacities for protoplanetary accretion disks: influence of dust aggregates. , keywords =

  36. [36]

    S., Rumpl, W., & Nordsieck, K

    The size distribution of interstellar grains. , keywords =. doi:10.1086/155591 , adsurl =

  37. [37]

    , keywords =

    Particle Emissivity in Circumstellar Disks. , keywords =. doi:10.1086/170646 , adsurl =

  38. [38]

    Protoplanetary Disks and Their Evolution

    Protoplanetary Disks and Their Evolution. , keywords =. doi:10.1146/annurev-astro-081710-102548 , archivePrefix =. 1103.0556 , primaryClass =

  39. [39]

    , keywords =

    Observations of Protoplanetary Disk Structures. , keywords =. doi:10.1146/annurev-astro-031220-010302 , archivePrefix =. 2001.05007 , primaryClass =

  40. [40]

    , keywords =

    Protoplanetary Disk Chemistry. , keywords =. doi:10.1146/annurev-astro-022823-040820 , archivePrefix =. 2309.05685 , primaryClass =

  41. [41]

    , year = 1983, month = sep, volume =

    The determination of cloud masses and dust characteristics from submillimetre thermal emission. , year = 1983, month = sep, volume =

  42. [42]

    , keywords =

    A Survey for Circumstellar Disks around Young Stellar Objects. , keywords =. doi:10.1086/115385 , adsurl =

  43. [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 =

  44. [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. [45]

    , keywords =

    Aerodynamics of solid bodies in the solar nebula. , keywords =. doi:10.1093/mnras/180.2.57 , adsurl =

  46. [46]

    , keywords =

    Dust Coagulation. , keywords =. doi:10.1086/172562 , adsurl =

  47. [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 =

  48. [48]

    , keywords =

    Panchromatic (Sub)millimeter polarization observations of HL Tau unveil aligned scattering grains. , keywords =. doi:10.1093/mnras/stae040 , archivePrefix =. 2309.10055 , primaryClass =

  49. [49]

    Scattering by Interstellar Dust Grains. I. Optical and Ultraviolet. , keywords =. doi:10.1086/379118 , archivePrefix =. astro-ph/0304060 , primaryClass =