pith. sign in

arxiv: 2510.21201 · v1 · submitted 2025-10-24 · ⚛️ physics.ao-ph · physics.optics

Light scattering by random convex polyhedron in geometric optics approximation

Quan Mu This is my paper

Pith reviewed 2026-05-18 05:20 UTC · model grok-4.3

classification ⚛️ physics.ao-ph physics.optics
keywords light scatteringice crystalsgeometric opticsconvex polyhedroncirrus cloudsscattering matrixpolarizationradiative transfer
0
0 comments X

The pith

A convex hull construction method generates models of any convex polyhedral ice crystal for which full light scattering matrices are computed in geometric optics.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a geometrical model of ice crystals based on the convex hull construction algorithm to study scattering properties relevant to cirrus cloud particles. Scattering matrices that include complete polarization information are calculated for randomly oriented large crystals having either random or prescribed convex polyhedron shapes. The construction method and associated computational scheme apply to any convex polyhedron in the geometric optics regime, demonstrated by calculations for three distinct shapes under one framework. The hexagonal column case produces results that agree with earlier work, and the approach targets radiative transfer simulations and remote sensing applications in planetary atmospheres.

Core claim

Based on the convex hull construction algorithm, a new geometrical model of ice crystals is proposed to investigate the scattering properties of cirrus clouds particles. Light scattering matrices involving complete polarization information are calculated in geometric optics approximation for randomly oriented large crystals with random and given convex polyhedron shape. The proposed model construction method and computational scheme of light scattering matrix works for any convex polyhedron within the scope of geometrical optics.

What carries the argument

Convex hull construction algorithm that produces the polyhedral geometry, followed by a unified geometric optics ray-tracing procedure to obtain the full scattering matrix.

If this is right

  • Scattering matrices for arbitrary convex shapes become computable under a single unified scheme.
  • Results for the classical hexagonal column match those reported by other authors.
  • The framework applies directly to radiative transfer simulations in terrestrial and planetary atmospheres.
  • Remote sensing data interpretation for cirrus clouds and similar particles can draw on the calculated matrices.

Where Pith is reading between the lines

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

  • The method could accommodate a broader range of irregular particle shapes observed in actual cloud samples.
  • Integration into existing atmospheric radiative transfer codes would allow testing against satellite observations of polarization signals.
  • For particles near the wavelength scale, adding diffraction corrections would be a direct extension of the current geometric optics treatment.

Load-bearing premise

The geometric optics approximation remains valid and sufficient for large randomly oriented convex polyhedral crystals without corrections for diffraction or absorption.

What would settle it

A measured or rigorously computed scattering matrix for a specific large convex polyhedron that deviates systematically from the model's prediction would show the geometric optics limit has been exceeded.

Figures

Figures reproduced from arXiv: 2510.21201 by Quan Mu.

Figure 1
Figure 1. Figure 1: Demonstration of mesh optimization: before (left) and after (right) merging coplanar trian [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Examples of random irregular convex polyhedron generated by the program [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Examples of regular convex polyhedron generated by the program [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Schematic representation of the incident, reflected, and refracted rays, together with the [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Comparison of six Mueller matrix elements for randomly oriented hexagonal column obtained [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Mueller matrix elements for faceted ellipsoid obtained by the program [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Mueller matrix elements for randomly constructed convex hull obtained by the program [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗
read the original abstract

Based on the convex hull construction algorithm, a new geometrical model of ice crystals is proposed to investigate the scattering properties of cirrus clouds particles. Light scattering matrices involving complete polarization information are calculated in geometric optics approximation for randomly oriented large crystals with random and given convex polyhedron shape. The proposed model construction method and computational scheme of light scattering matrix works for any convex polyhedron within the scope of geometrical optics. To illustrate the broad applicability of the proposed ice crystal model, scattering matrices for three ice crystal examples with different geometrical shapes are calculated under a unified computational framework. Diffraction and absorption are not considered in this work. The calculated results for the classical hexagonal column model show overall agreement with those reported by other authors. The crystal model and scattering matrix computational framework developed in this study are applicable to radiative transfer simulations and remote sensing data interpretation in terrestrial and planetary atmospheres.

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 / 2 minor

Summary. The manuscript proposes a convex-hull-based geometrical model for ice crystals and a ray-tracing scheme to compute full scattering matrices (including polarization) in the geometric optics approximation for randomly oriented large convex polyhedra. The central claim is that the construction method and computational framework apply to any convex polyhedron; this is illustrated with three example shapes, with the hexagonal-column case showing overall agreement to prior literature. Diffraction and absorption are explicitly omitted.

Significance. If the generality of the scheme is established, the work supplies a flexible forward model for scattering by arbitrary convex ice-crystal shapes that could support radiative-transfer simulations and remote-sensing retrievals in cirrus clouds. The convex-hull approach for generating random shapes and the unified computational framework are strengths when the implementation proves robust across facet counts.

major comments (2)
  1. [Abstract] Abstract: the claim that the model construction method and computational scheme 'works for any convex polyhedron' is central yet supported only by three specific examples. Because orientation averaging and multi-bounce ray accounting both scale with facet number and arrangement, additional test cases (high-face-count or near-degenerate polyhedra) or pseudocode demonstrating algorithmic completeness are required to substantiate universality.
  2. [Numerical results] Numerical results section: the validation against prior hexagonal-column results is described only as 'overall agreement'; quantitative differences in individual scattering-matrix elements and error estimates from the orientation averaging are not reported, weakening the assessment of numerical accuracy for the new scheme.
minor comments (2)
  1. [Abstract] The abstract states that diffraction and absorption are not considered; this scope limitation should be restated explicitly in the conclusions to prevent readers from overgeneralizing the applicability to smaller particles.
  2. [Figures] Figure captions and axis labels for the scattering-matrix plots should include the wavelength and size parameter used, as these are essential for reproducibility in geometric-optics calculations.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments, which help clarify the scope and validation of our work. We address each major comment below.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the claim that the model construction method and computational scheme 'works for any convex polyhedron' is central yet supported only by three specific examples. Because orientation averaging and multi-bounce ray accounting both scale with facet number and arrangement, additional test cases (high-face-count or near-degenerate polyhedra) or pseudocode demonstrating algorithmic completeness are required to substantiate universality.

    Authors: The convex-hull construction is based on standard, general-purpose algorithms that generate a convex polyhedron from any finite set of vertices in 3D space, independent of the specific number or arrangement of facets. The subsequent ray-tracing scheme operates on the resulting facet list and connectivity without assuming particular symmetries or bounded facet counts. Nevertheless, to strengthen the demonstration of universality, we will add pseudocode for the core construction and ray-tracing steps in an appendix and include scattering results for one additional high-facet-count polyhedron in the revised manuscript. revision: yes

  2. Referee: [Numerical results] Numerical results section: the validation against prior hexagonal-column results is described only as 'overall agreement'; quantitative differences in individual scattering-matrix elements and error estimates from the orientation averaging are not reported, weakening the assessment of numerical accuracy for the new scheme.

    Authors: We agree that a purely qualitative statement limits the ability to assess accuracy. In the revised manuscript we will report quantitative comparisons, specifically the maximum relative difference for each non-zero element of the scattering matrix against the reference literature values, together with the estimated statistical uncertainty arising from the finite number of orientations used in the Monte Carlo averaging. revision: yes

Circularity Check

0 steps flagged

No significant circularity; forward geometric model independent of target data

full rationale

The paper constructs a forward scattering model from convex-hull geometry and standard geometric-optics ray tracing rules. Scattering matrices are computed directly from polyhedron facet definitions and orientation averaging; no parameters are fitted to the output scattering data itself. The hexagonal-column validation references external literature rather than a self-citation chain that would close the loop. The claim of applicability to arbitrary convex polyhedra is a statement of algorithmic generality, not a definitional reduction. No equations or steps reduce the reported results to inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

The central claim rests on the validity of the geometric-optics ray-tracing rules for large particles and on the assumption that convex-hull shapes adequately represent real ice-crystal geometry; no explicit free parameters are introduced in the abstract, but the random-point generation step implicitly contains choices for point distribution and number.

axioms (1)
  • domain assumption Geometric optics approximation is valid for the modeled large, randomly oriented crystals
    Abstract states diffraction and absorption are not considered and that the scheme operates within geometric optics.
invented entities (1)
  • Random convex polyhedron ice-crystal model via convex hull no independent evidence
    purpose: Generate arbitrary convex shapes for scattering calculations
    New geometrical model proposed to replace or extend classical hexagonal forms.

pith-pipeline@v0.9.0 · 5664 in / 1326 out tokens · 29577 ms · 2026-05-18T05:20:50.009926+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

28 extracted references · 28 canonical work pages

  1. [1]

    Light scattering by nonspher- ical particles: theory, measurements, and applications.Measurement Science and Technology, 11(12):1827–1827, 2000

    Michael I Mishchenko, Joop W Hovenier, and Larry D Travis. Light scattering by nonspher- ical particles: theory, measurements, and applications.Measurement Science and Technology, 11(12):1827–1827, 2000

    work page 2000

  2. [2]

    Cambridge University Press, 2016

    Kuo-Nan Liou and Ping Yang.Light scattering by ice crystals: fundamentals and applications. Cambridge University Press, 2016

    work page 2016

  3. [3]

    The dependence of phase functions of large transparent particles on their refractive index and shape.Journal of Physics D: Applied Physics, 31(11):1329, 1998

    AA Kokhanovsky and TY Nakajima. The dependence of phase functions of large transparent particles on their refractive index and shape.Journal of Physics D: Applied Physics, 31(11):1329, 1998

    work page 1998

  4. [4]

    Applicability of regular particle shapes in light scat- tering calculations for atmospheric ice particles.Applied optics, 35(21):4291–4296, 1996

    Andreas Macke and Michael I Mishchenko. Applicability of regular particle shapes in light scat- tering calculations for atmospheric ice particles.Applied optics, 35(21):4291–4296, 1996

    work page 1996

  5. [5]

    Radiative transfer in cirrus clouds

    Y Takano and KN Liou. Radiative transfer in cirrus clouds. part iii: Light scattering by irregular ice crystals.Journal of Atmospheric Sciences, 52(7):818–837, 1995

    work page 1995

  6. [6]

    Grynko and Yu

    Ye. Grynko and Yu. Shkuratov. Scattering matrix calculated in geometric optics approximation for semitransparent particles faceted with various shapes.Journal of Quantitative Spectroscopy and Radiative Transfer, 78(3):319–340, 2003

    work page 2003

  7. [7]

    Trigonal ice crystals in earth’s atmosphere.Bulletin of the American Meteorological Society, 96(9):1519–1531, 2015

    Benjamin J Murray, Christoph G Salzmann, Andrew J Heymsfield, Steven Dobbie, Ryan R Neely III, and Christopher J Cox. Trigonal ice crystals in earth’s atmosphere.Bulletin of the American Meteorological Society, 96(9):1519–1531, 2015

    work page 2015

  8. [8]

    A review of ice cloud optical property models for passive satellite remote sensing.Atmosphere, 9(12):499, 2018

    Ping Yang, Souichiro Hioki, Masanori Saito, Chia-Pang Kuo, Bryan A Baum, and Kuo-Nan Liou. A review of ice cloud optical property models for passive satellite remote sensing.Atmosphere, 9(12):499, 2018

    work page 2018

  9. [9]

    The effective equivalence of geometric irregular- ity and surface roughness in determining particle single-scattering properties.Optics express, 22(19):23620–23627, 2014

    Chao Liu, R Lee Panetta, and Ping Yang. The effective equivalence of geometric irregular- ity and surface roughness in determining particle single-scattering properties.Optics express, 22(19):23620–23627, 2014

    work page 2014

  10. [10]

    Single scattering properties of atmo- spheric ice crystals.Journal of Atmospheric Sciences, 53(19):2813–2825, 1996

    Andreas Macke, Johannes Mueller, and Ehrhard Raschke. Single scattering properties of atmo- spheric ice crystals.Journal of Atmospheric Sciences, 53(19):2813–2825, 1996

    work page 1996

  11. [11]

    The effects of surface roughness on the scattering properties of hexagonal columns with sizes from the rayleigh to the geometric optics regimes

    Chao Liu, R Lee Panetta, and Ping Yang. The effects of surface roughness on the scattering properties of hexagonal columns with sizes from the rayleigh to the geometric optics regimes. Journal of Quantitative Spectroscopy and Radiative Transfer, 129:169–185, 2013

    work page 2013

  12. [12]

    Modelling lidar-relevant optical properties of complex mineral dust aerosols.Tellus B: Chemical and Physical Meteorology, 63(4):725–741, 2011

    Josef Gasteiger, Matthias Wiegner, Silke Groß, Volker Freudenthaler, Carlos Toledano, Matthias Tesche, and Konrad Kandler. Modelling lidar-relevant optical properties of complex mineral dust aerosols.Tellus B: Chemical and Physical Meteorology, 63(4):725–741, 2011

    work page 2011

  13. [13]

    Coherent and incoherent backscattering by a single large particle of irregular shape.Optics Express, 27(23):32984–32993, 2019

    Victor Shishko, Alexander Konoshonkin, Natalia Kustova, Dmitry Timofeev, and Anatoli Borovoi. Coherent and incoherent backscattering by a single large particle of irregular shape.Optics Express, 27(23):32984–32993, 2019. 12

    work page 2019

  14. [14]

    Computer-aided construction of three-dimensional convex bodies of arbitrary shapes.Computational technologies, 27(2):54–61, 2022

    Q Mu, BA Kargin, and EG Kablukova. Computer-aided construction of three-dimensional convex bodies of arbitrary shapes.Computational technologies, 27(2):54–61, 2022

    work page 2022

  15. [15]

    Numerical stochastic simulation of optical radiation scattering by ice crystals of irregular random shapes.Computational Technologies, 27(2):4–18, 2022

    BA Kargin, EG Kablukova, and Q Mu. Numerical stochastic simulation of optical radiation scattering by ice crystals of irregular random shapes.Computational Technologies, 27(2):4–18, 2022

    work page 2022

  16. [16]

    Springer Science & Business Media, 2012

    Franco P Preparata and Michael I Shamos.Computational geometry: an introduction. Springer Science & Business Media, 2012

    work page 2012

  17. [17]

    Springer, 2008

    Mark De Berg, Otfried Cheong, Marc Van Kreveld, and Mark Overmars.Computational geometry: algorithms and applications. Springer, 2008

    work page 2008

  18. [18]

    Monte carlo method for numerical simu- lation of solar energy radiation transfer in crystal clouds.Numerical Analysis and Applications, 17(2):140–151, 2024

    BA Kargin, EG Kablukova, Q Mu, and SM Prigarin. Monte carlo method for numerical simu- lation of solar energy radiation transfer in crystal clouds.Numerical Analysis and Applications, 17(2):140–151, 2024

    work page 2024

  19. [19]

    Kattawar, Si-Chee Tsay, Bryan A

    Zhibo Zhang, Ping Yang, George W. Kattawar, Si-Chee Tsay, Bryan A. Baum, Yongxiang Hu, Andrew J. Heymsfield, and Jens Reichardt. Geometrical-optics solution to light scattering by droxtal ice crystals.Appl. Opt., 43(12):2490–2499, Apr 2004

    work page 2004

  20. [20]

    Light scattering and absorption by nonspherical ice crystals

    Ping Yang and Kuo-Nan Liou. Light scattering and absorption by nonspherical ice crystals. In Light Scattering Reviews: Single and Multiple Light Scattering, pages 31–71. Springer, 2006

    work page 2006

  21. [21]

    Principles of optics (60th anniversary of first edition, 20th anniversary of seventh edition), 2019

    M Born, E Wolf, and P Knight. Principles of optics (60th anniversary of first edition, 20th anniversary of seventh edition), 2019

    work page 2019

  22. [22]

    Courier Corporation, 1981

    Hendrik Christoffel Hulst and Hendrik C van de Hulst.Light scattering by small particles. Courier Corporation, 1981

    work page 1981

  23. [23]

    Springer Science & Business Media, 2014

    Joop W Hovenier, Cornelis VM Van der Mee, and Helmut Domke.Transfer of polarized light in planetary atmospheres: basic concepts and practical methods, volume 318. Springer Science & Business Media, 2014

    work page 2014

  24. [24]

    John Wiley & Sons, 2007

    Hiroyuki Fujiwara.Spectroscopic ellipsometry: principles and applications. John Wiley & Sons, 2007

    work page 2007

  25. [25]

    Beam splitting algorithm for the problem of light scattering by atmospheric ice crystals

    AV Konoshonkin, NV Kustova, and AG Borovoi. Beam splitting algorithm for the problem of light scattering by atmospheric ice crystals. part 1. theoretical foundations of the algorithm. Atmospheric and Oceanic Optics, 28(5):441–447, 2015

    work page 2015

  26. [26]

    Beam splitting algorithm for the problem of light scattering by atmospheric ice crystals

    AV Konoshonkin, NV Kustova, and AG Borovoi. Beam splitting algorithm for the problem of light scattering by atmospheric ice crystals. part 2. comparison with the ray tracing algorithm. Atmospheric and Oceanic Optics, 28:448–454, 2015

    work page 2015

  27. [27]

    Scattering of light by polyhedral ice crystals.Applied Optics, 32(15):2780–2788, 1993

    Andreas Macke. Scattering of light by polyhedral ice crystals.Applied Optics, 32(15):2780–2788, 1993

    work page 1993

  28. [28]

    A. Macke. rt-crystal, 2020. Available athttps://doi.org/10.5281/zenodo.3965488. 13

    work page doi:10.5281/zenodo.3965488 2020