REVIEW 3 major objections 4 minor 43 references
Combining color and polarization maps with dust modeling, this paper argues that comet 67P's 2021/22 coma was dominated by particles larger than 10 micrometers and that unusually high polarization points to changed dust properties since 201
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 →
T0 review · deepseek-v4-flash
2026-08-01 10:58 UTC pith:DPKDV4QE
load-bearing objection Solid new 67P color/polarization maps from the 2021/22 apparition; the observational work is careful and worth having, but the >10 μm dust-size conclusion is model-dependent and the paper itself gives the reason. the 3 major comments →
A comprehensive study of comet 67P/Churyumov-Gerasimenko in the 2021/2022 apparition. II. Colorimetry, polarimetry, modeling
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The paper establishes that the observed spatial gradients in color (red near the nucleus, bluer outward) and polarization (increasing with distance) are the signature of dynamic sorting of dust particles, not primarily of compositional zoning. Using ballistic particle-cluster aggregates made of equal parts silicates and processed organics, the model reproduces the observed maps only when particles larger than 10 micrometers dominate the scattering cross section at all distances; smaller grains in the 1-10 micrometer range are confined to small regions where polarization is locally enhanced. The same data show that the 2021/22 polarization exceeds the reference curve for high-polarization com
What carries the argument
The central object is the ballistic particle-cluster aggregate (BPCA), a fluffy dust grain built from sub-micron monomers, which reproduces the negative polarization branch. Light scattering is computed exactly with a fast superposition T-matrix method for small particles and with a Monte Carlo radiative-transfer solver for larger ones; a Monte Carlo dust-dynamics model adds gas drag, solar gravity, and radiation pressure. The dynamic-sorting mechanism - small grains accelerated faster and blown outward, large grains lingering near the nucleus and tail - is what turns a single size distribution into the observed spatial color and polarization pattern.
Load-bearing premise
The size conclusion depends on the assumed grain optics - specifically, a fairly absorbing mix of silicates and processed organics with a particular aggregate structure; if the grains are less absorbing, the polarization minimum shifts to much larger sizes and the 'small grains only in patches' reading loses its footing.
What would settle it
Simultaneous color and polarization imaging at an intermediate phase angle near 30 degrees would test the model's predicted size-color-polarization correlation; a mismatch in sign or magnitude would rule out the assumed grain optics. A mid-infrared spectrum checking for a strong 10-micrometer silicate feature would also test whether small silicate grains are actually abundant.
If this is right
- In 2021/22, the coma's optical properties were mostly set by grains larger than 10 micrometers, so remote-sensing interpretations for this apparition should not assume a small-grain-dominated coma.
- Local polarization enhancements are tracers of narrow regions with 1-10 micrometer grains or less-absorbing material, giving a way to map fine-grain enrichment from polarimetry alone.
- The unusually high positive polarization places 67P at the upper end of high-polarization comets, which bears on how processed organic material is distributed in its dust.
- Because the 2021/22 coma was smoother in color and polarization than the 2015/16 one, the change in dust properties between apparitions is a measurable, testable trend for future returns.
Where Pith is reading between the lines
- If future mid-infrared observations show a strong 10-micrometer silicate feature, the large-grain dominance would be hard to sustain; a weak feature would support the authors' interpretation. This is a direct, testable consequence of the paper's size claim.
- The same combined light-scattering-plus-dynamics approach could be applied to other Jupiter-family comets with smooth comae; the paper's method effectively inverts color-polarization maps into a size-and-composition field.
- The positive color-polarization correlation the authors see in the southern hemisphere suggests compositional (not purely size) variations; a mapping campaign at multiple phase angles could separate size from composition more cleanly than this single-return dataset allows.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper reports new imaging colorimetry and linear polarimetry of comet 67P/Churyumov–Gerasimenko obtained during the 2021/22 apparition, together with aperture polarimetry from two additional telescopes. The data show a red inner coma (g−r ≈ 0.8 mag) that becomes bluer with distance, and polarization of about 10–12% at α = 47.9° and about −1% to −2.5% at α = 10.5°, with a more homogeneous spatial distribution than in 2015/16. The authors model the observations with BPCA aggregates of equal-mass silicates and processed organics, using FaSTMM/SIRIS4 light scattering and a Monte Carlo dynamical dust model. They conclude that the coma is dominated by particles larger than 10 μm, that local polarization enhancements are due to 1–10 μm grains or less-absorbing material, and that the 2021/22 polarization exceeds the reference high-Pmax comet curve, suggesting a change in dust properties since 2015/16.
Significance. The observational dataset is valuable: the color and polarization maps are carefully reduced, gas contamination is estimated spectroscopically and found to be less than 10%, and the aperture polarimetry extends the phase-angle coverage of 67P across four apparitions. The modeling effort is sophisticated, combining asymptotically exact light-scattering computations with an explicit dynamical model, and the paper is transparent about most parameters. If the size inference held, it would be strong evidence for inter-apparition dust evolution. However, the central size conclusion rests on a specific assumed refractive index and aggregate morphology, and the paper's own Section 5.4 acknowledges an alternative that would remove the need for small grains. The model also fails to reproduce several observed features. The color/polarization maps themselves should be publishable; the size interpretation needs either sensitivity analysis or a more cautious statement.
major comments (3)
- [Section 5.4, Fig. 8] The conclusion that the coma is dominated by >10 μm particles and that local polarization enhancements require 1–10 μm grains is directly contingent on the assumed refractive index (m=1.754+i0.157 at 440 nm, m=1.778+i0.133 at 650 nm) and BPCA morphology. The text itself notes that less-absorbing particles shift the polarization minimum to ~10–100 μm, which would make the outer-coma negative polarization consistent with large grains. Since this refractive index was selected to reproduce the aperture-averaged polarization, color, and albedo (Section 5.2), the inference is partly circular. A sensitivity analysis over a plausible range of complex refractive indices and porosities should be provided, or the claim should be explicitly re-phrased as conditional on this composition/morphology.
- [Section 5.4, radial profiles] The same section reports that the model overestimates the radial increase of polarization and underestimates the radial decrease of color, and cannot reproduce the 6 February 2022 negative polarization or the small-scale jet structures. Given these discrepancies, the fitted parameters (r_min=0.25 μm, power-law index −3.3, r_max=2.6 mm, ν=1.1, γ=0.175) are not strongly constrained. The paper should include quantitative residuals of the model versus observed radial profiles and discuss whether the imperfect fit biases the inferred size threshold.
- [Section 5.3–5.4, activity history] The dynamical model assumes continuous activity began 12 months before the observation and adopts Qgas ∝ rh^-6.5 from Rosetta's 2015/16 apparition. Since the observed coma is a convolution of the emission history with particle velocities and radiation pressure, the conclusion that only large particles remain in the tail depends on this assumption. A sensitivity test with different activity onset times and gas production scalings (or at least an explicit justification) is needed to establish that the >10 μm result is not an artifact of the assumed activity history.
minor comments (4)
- [Section 5.2] Typo: 'the resulting refractive using the Maxwell-Garnett mixing rule' should read 'the resulting refractive index using the Maxwell-Garnett mixing rule'.
- [Section 4.1] The phrase 'the polarization decreases sharply from ~1% near the nucleus to ~4%' is ambiguous for negative polarization; suggest 'the absolute value of the negative polarization increases from about 1% near the nucleus to 4%'.
- [Section 5.1] The sentence 'for particles larger than up to 1 cm' is awkward; suggest 'for particles up to 1 cm'.
- [Data availability] The statement 'Data will be made available on request' is weaker than standard practice; consider depositing the reduced maps and polarimetric tables in a permanent repository.
Circularity Check
No circular derivation found; the >10 μm dust-size conclusion is model-dependent and explicitly acknowledged as such, not a reduction to the model's inputs.
full rationale
The derivation chain is: observed color/polarization maps -> light-scattering (FaSTMM/SIRIS4) and dynamical Monte Carlo models with assumed BPCA morphology, monomer size distribution, 50/50 silicate/organic composition, refractive index m=1.754+i0.157 (440 nm) and m=1.778+i0.133 (650 nm), and free dynamical and size-distribution parameters (nu, gamma, rmin, power-law index, rmax) fitted to the observations. The central claim that the coma is dominated by particles >10 um is read off the computed Fig. 8 saturation of color/polarization with size for the chosen model. This is an inverse-modeling inference, i.e., a fitted/assumption-dependent result, but it is not circular: no quantity is defined in terms of the target conclusion, and no fitted parameter is relabeled as an independent prediction. The paper itself flags the degeneracy in Section 5.4: 'One possible explanation is the presence of less absorbing particles, for which the polarization minimum shifts to larger sizes (~10-100 um), allowing the observations to be reconciled with the dynamical constraints.' That is an honest statement of model non-uniqueness (a correctness risk), not a circular reduction. Self-citations to the modeling codes and to prior apparition papers are methodological or data references and are not used as an unverified uniqueness theorem. Therefore no specific circular step can be exhibited, and the appropriate score is 0.
Axiom & Free-Parameter Ledger
free parameters (5)
- Effective refractive index of organic-silicate mixture =
m = 1.754 + i0.157 (440 nm); m = 1.778 + i0.133 (650 nm)
- Minimum particle size r_min =
0.25 μm
- Power-law index of size distribution =
-3.3
- Maximum particle size r_max =
2.6 mm
- Gas-drag velocity parameters ν, γ =
ν = 1.1, γ = 0.175
axioms (5)
- domain assumption Dust particles are moderately porous BPCA aggregates of spherical monomers with log-normal size distribution (mean 50 nm, σ = 30 nm)
- domain assumption Refractive index of the mixture follows Maxwell-Garnett mixing with 50/50 organics/silicates
- domain assumption Dust emission follows a cosine law from a spherical surface of radius 10 km with no night-side emission, and gas production ∝ r_h^−6.5
- ad hoc to paper Continuous activity began 12 months before the observation date
- domain assumption The approximate Monte Carlo radiative transfer solver SIRIS4 is valid for particles >5 μm
Cite this review
Pith. "Pith review of A comprehensive study of comet 67P/Churyumov-Gerasimenko in the 2021/2022 apparition. II. Colorimetry, polarimetry, modeling." pith.science (2026). https://pith.science/paper/DPKDV4QE
@misc{pith2026260720038,
author = {Pith},
title = {Pith review of: A comprehensive study of comet 67P/Churyumov-Gerasimenko in the 2021/2022 apparition. II. Colorimetry, polarimetry, modeling},
year = {2026},
howpublished = {\url{https://pith.science/paper/DPKDV4QE}},
note = {Machine review of arXiv:2607.20038}
}
read the original abstract
This paper presents the second part of our multi-technique observational study of comet 67P during its 2021-2022 apparition. While the first paper focused on photometry and spectroscopy, here we investigate the dust color and linear polarization of the coma. Imaging in the g-sdss and r-sdss filters was obtained with the 6-m BTA SAO Telescope, and aperture polarimetry was performed with the 2.6-m Shajn Telescope (CrAO) and the 2-m RCC Telescope. Observations at phase angles of 47.9 deg (6 October 2021) and 10.5 deg (6 February 2022) show negative polarization of about 0.5-2.5% at 10.5 deg and positive polarization of about 11% at 47.9 deg. The dust color was red, reaching about 0.8 mag near the nucleus and in the tail, and decreasing to about 0.2 mag in the coma and jet within 20,000 km of the nucleus. Spatial maps reveal a gradual decrease in color index and a corresponding increase in polarization with distance from the nucleus. Numerical modeling using ballistic particle-cluster aggregates composed of equal fractions of silicates and organics indicates that the coma is dominated by particles larger than 10 micrometers, while local polarization enhancements are consistent with regions containing smaller (1-10 micrometer) grains. Combining our aperture measurements with published data from four previous apparitions confirms that comet 67P belongs to the high-polarization class. The unusually high positive polarization observed during the 2021/22 apparition, together with the relatively isotropic coma morphology, suggests that the dust properties may have changed after the 2015/16 apparition.
Figures
Reference graph
Works this paper leans on
-
[11]
Imaging polarimetry of comet 67P/ Churyumov–Gerasimenko: homogeneous distribution of polarization and its implications. Mon. Not. R. Astron. Soc. 531, 1638–1652. https://doi.org/10.48550/ arXiv.2405.09297. Groussin, O., Attree, N., Brouet, Y., et al.,
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2405.09297
-
[18]
The global shape, density and rotation of comet 67P/Churyumov-Gerasimenko from preperihelion Rosetta/OSIRIS observations. Icarus 277, 257–278. https://doi.org/10.1016/j.icarus.2016.05.002. Keller, H.U., Mottola, S., Hviid, S.F., et al.,
-
[19]
Seasonal mass transfer on the nucleus of comet 67P/Chuyumov–Gerasimenko. Mon. Not. R. Astron. Soc. 469, S357–S371. https://doi.org/10.1093/mnras/stx1726. Kelley, M.S., Woodward, C.E., Harker, D.E., et al.,
-
[24]
Interaction of Electromagnetic Radiation with Cometary Dust
Interaction of electromagnetic radiation with cometary dust. In: Comets III. Arizona University Press, pp. 621–652. https://doi.org/10.48550/arXiv.2501.03117. Kwon, Y.G., Bagnulo, S., Markkanen, J., et al.,
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2501.03117
-
[25]
VLT spectropolarimetry of comet 67P: dust environment around the end of its intense southern summer. Astron. Astrophys. 657, A40. https://doi.org/10.1051/0004-6361/202141865. Larson, S.M., Sekanina, Z.,
-
[27]
The gas production of 14 species from comet 67P/Churyumov–Gerasimenko based on DFMS/COPS data from 2014 to
2014
-
[28]
Mon. Not. R. Astron. Soc. 498, 3995–4004. https://doi.org/10.1093/mnras/ staa2643. Levasseur-Regourd, A.C., Agarwal, J., Cottin, H., et al.,
-
[30]
Dust of comet 67P/Churyumov- Gerasimenko collected by Rosetta/ MIDAS: classification and extension to the nanometer scale. Astron. Astrophys. 630, A26. https://doi.org/10.1051/0004-6361/ 201834851. Markkanen, J., Yuffa, A.J.,
-
[31]
Fast superposition T-matrix solution for clusters with arbitrarily-shaped constituent particles. J. Quant. Spectrosc. Radiat. Transf. 189, 181–188. https://doi.org/10.1016/j.jqsrt.2016.11.004. Markkanen, J., Agarwal, J., V¨ais¨anen, T., Penttil¨a, A., Muinonen, K.,
-
[32]
Astrophys
Interpretation of the phase functions measured by the OSIRIS instrument for comet 67P/ Churyumov-Gerasimenko. Astrophys. J. Let. 868, L16. https://doi.org/10.3847/ 2041-8213/aaee10. Myers, R.V., Nordsieck, K.H.,
2041
-
[34]
Spatial variations of brightness, colour and polarization of dust in comet 67P/Churyumov–Gerasimenko. Mon. Not. R. Astron. Soc. 469, S475–S491. https://doi.org/10.1093/mnras/stx2003. Rosenbush, V., Ivanova, O., Kleshchonok, V., et al.,
-
[35]
Imaging photometry and long-slit spectroscopy
Comet 2P/Encke in apparitions of 2013 and 2017: I. Imaging photometry and long-slit spectroscopy. Icarus 348, 113767. https://doi.org/10.1016/j.icarus.2020.113767. Rosenbush, V., Kleshchonok, V., Ivanova, O., et al.,
arXiv 2013
-
[39]
On the nucleus structure and activity of comet 67P/Churyumov–Gerasimenko. Science 347, aaa1044. https://doi.org/ 10.1126/science.aaa1044. Stinson, A., Bagnulo, S., Tozzi, G.P., et al.,
-
[40]
Polarimetry of comets 67P/ Churyumov–Gerasimenko, 74P/Smirnova–Chernykh, and 152P/Helin–Lawrence. Astron. Astrophys. 594, A110. https://doi.org/10.1051/0004-6361/201527696. Thomas, N., Davidsson, B., Jorda, L., Kührt, E., Marschall, R., Snodgrass, C., Rodrigo, R. (Eds.),
-
[41]
Proceedings of SPIE Observatory Operations: Strategies, Processes, and Systems VII
SALT and SAAO strategy, focusing on the time-domain: process, plans, and challenges. Proceedings of SPIE Observatory Operations: Strategies, Processes, and Systems VII. 107040A. https://doi.org/ 10.1117/12.2313272. Willmer, C.N.A.,
-
[42]
The absolute magnitude of the sun in several filters. Astrophys. J. Suppl. Ser. 236, 47, 14 pp. https://doi.org/10.3847/1538-4365/aabfdf. Agarwal, J., Kim, Y., Kelley, M., Marschall, R.,
-
[43]
Dust emission and dynamics. In: Meech, K.J., et al. (Eds.), Comets III. University of Arizona Press, pp. 653–678. https://doi.org/10.2458/azu_uapress_9780816553631-ch020. Attree, N., Guti´errez, P., Schuckart, C., et al.,
-
[44]
Geometrical model of jets in cometary comae. Icarus 425, 116300. https://doi.org/10.1016/j.icarus.2024.116300. Kleshchonok, V., Sierks, H., Güttler, C.,
arXiv 2024
-
[347]
https://doi.org/10.1126/ science.aaa0276 id. aaa0276. Ivanova, O.V., Rosenbush, V.K., Kiselev, N.N., et al., 2017a. Post-perihelion observations of comet 67P/Churyumov–Gerasimenko at the 6 m BTA telescope: optical spectroscopy. Mon. Not. R. Astron. Soc. 469, S386–S395. https://doi.org/10.1093/ mnras/stx1725. Ivanova, O., Rosenbush, V., Afanasiev, V., et a...
-
[1910]
Astron. J. 89, 571–578. https://doi.org/10.1086/113551. L¨auter, M., Kramer, T., Rubin, M., Altwegg, K.,
-
[1982]
On standard polarized stars. Astrophys. J. 262, 732–738. https://doi.org/10.1086/160467. Chernova, G.P., Kiselev, N.N., Jockers, K.,
-
[1984]
Spectropolarimetry of comets Austin and Churyumov-Gerasimenko. Icarus 58, 431–439. https://doi.org/10.1016/0019-1035 (84)90088-5. Rosenbush, V.K., Ivanova, O.V., Kiselev, N.N., et al.,
-
[1992]
The Hubble space telescope northern- hemisphere grid of stellar polarimetric standards. Astron. J. 104, 1563–1567. https://doi.org/10.1086/116341. Schulz, R., Stüwe, J.A., Boehnhardt, H.,
-
[1993]
Polarimetric characteristics of dust particles as observed in 13 comets: comparisons with asteroids. Icarus 103, 144–158. https://doi.org/10.1006/icar.1993.1063. Dahlen, D., Kwon, Y.G., Masiero, J.R., et al.,
arXiv 1993
-
[2000]
The HB narrowband comet filters: standardizing calibration and analysis techniques. Icarus 147, 180–204. https://doi. org/10.1006/icar.2000.6420. Fornasier, S., Hasselmann, P.H., Barucci, M.A., et al.,
arXiv 2000
-
[2004]
Rosetta target comet 67P/Churyumov- Gerasimenko Postperihelion gas and dust production rates. Astron. Astrophys. 422, L19–L21. https://doi.org/10.1051/0004-6361:20040190. Serkowski, K.,
-
[2006]
A Spitzer study of comets 2P/ Encke, 67P/Churyumov-Gerasimenko, and C/2001 HT50 (LINEAR-NEAT). Astrophys. J. 651, 1256–1271. https://doi.org/10.48550/arXiv.astro-ph/0607416. Kiselev, N., Rosenbush, V., Levasseur-Regourd, A., Kolokolova, L.,
-
[2007]
Two different evolutionary types of comets proved by polarimetric and infrared properties of their dust. Astron. Astrophys. 463, 1189–1196. https://doi.org/10.1051/0004-6361:20065069. Kolokolova, L., Kelley, M.S., Kimura, H., Hoang, T.,
-
[2010]
Astron. Astrophys. 517, A86. https://doi.org/10.1051/0004-6361/201014167. Hadamcik, E., Levasseur-Regourd, A.C., Hines, D.C., et al.,
-
[2011]
SCORPIO on the 6m telescope: current state and perspectives for spectroscopy of galactic and extragalactic objects. Balt. Astron. 20, 363–370. https://doi.org/10.1515/astro-2017-0305. Afanasiev, V.L., Rosenbush, V.K., Kiselev, N.N.,
-
[2012]
Technique of polarimetric observations of faint objects at the 6-m BTA telescope. Astrophys. Bull. 67, 438–448. https://doi. org/10.1134/S1990341312040074. Afanasiev, V.L., Moiseev, A.V.,
-
[2014]
Polarimetry of major Uranian moons at the 6-m telescope. Astrophys. Bull. 69, 211–223. https://doi.org/10.1134/ S1990341314020096. V. Rosenbush et al. 14 activity model on comet 67P/Churyumov–Gerasimenko. Mon. Not. R. Astron. Soc. 541, 771–783. https://doi.org/10.1093/mnras/staf1040. Bardyn, A., Baklouti, D., Cottin, H., et al.,
-
[2015]
Spectrophotometric properties of the nucleus of comet 67P/Churyumov-Gerasimenko from the OSIRIS instrument onboard the ROSETTA spacecraft. Astron. Astrophys. 583, A30, 18 pp. https://doi.org/10.1051/0004-6361/201525901. Gray, Z., Bagnulo, S., Boehnhardt, H., et al.,
-
[2016]
Properties of dust particles in comets from photometric and polarimetric observations of 67P. Mon. Not. R. Astron. Soc. S507–S515. https://doi.org/10.1093/mnras/stx030. H¨assig, M., Altwegg, K., Balsiger, et al.,
-
[2017]
Carbon-rich dust in comet 67P/Chur - yumov-Gerasimenko measured by COSIMA/Rosetta. Mon. Not. R. Astron. Soc. 469, S712–S722. https://doi.org/10.1093/mnras/stx2640. Boehnhardt, H., Lara, L., Gray, Z., Bagnulo, S.,
-
[2018]
Cometary dust. Space Sci. Rev. 214, 64, 56 pp. https://doi.org/10.1007/s11214-018-0496-3. Li, A., Greenberg, J.M.,
-
[2019]
The thermal, mechanical, structural, and dielectric properties of cometary nuclei after Rosetta. Space Sci. Rev. 215, 29, 51 pp. https://doi.org/10.1007/s11214-019-0594-x. Hadamcik, E., Sen, A.K., Levasseur-Regourd, A.C., et al.,
-
[2021]
Observations of distant comet C/ 2011 KP36 (Spacewatch): photometry, spectroscopy, and polarimetry. Astron. Astrophys. 651, A29. https://doi.org/10.1051/0004-6361/202039668. Ivanova, O., Rosenbush, V., Luk'yanyk, I., et al.,
-
[2022]
New polarimetric data for the Galilean satellites: Europa observations and modeling. Planet. Sci. J. 3, 134, 13 pp. https://doi.org/10.3847/PSJ/ac6bef. Churyumov –Gerasimenko during increased activity at perihelion on 2015 august 22–23. Mon. Not. R. Astron. Soc. 542, 241–255. https://doi.org/10.1093/mnras/ staf1077. Kolokolova, L., Jockers, K., Gustafson,...
-
[2023]
Quasi-simultaneous photometric, polarimetric, and spectral observations of distant comet C/2014 B1 (Schwartz). Astron. Astrophys. 672, A76. https://doi.org/10.1051/0004-6361/202244686. Jenniskens, P.,
-
[2024]
Coma and tail of comet 67P/ Churyumov-Gerasimenko during the 2021-2022 apparition. Mon. Not. R. Astron. Soc. 531, 3912–3926. https://doi.org/10.1093/mnras/stae1412. Hsu, J.C., Breger, V.,
-
[2025]
Kete: predicting known minor bodies in images. arXiv:2509.04666 astro-ph.EP. https://doi.org/10.48550/ arXiv.2509.04666. Dorschner, J., Begemann, B., Henning, Th., et al.,
-
[2026]
A comprehensive study of comet 67P/Churyumov-Gerasimenko in the 2021/2022 apparition. I. Photometry, spectroscopy, morphology. Icarus 444, 116799. https://doi.org/10.1016/j. icarus.2025.116799. Schmidt, G.D., Elston, R., Lupie, O.L.,
arXiv 2021
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.