pith. sign in

arxiv: 2606.21063 · v1 · pith:4M7VJ5XVnew · submitted 2026-06-19 · ⚛️ physics.optics · physics.ao-ph· physics.comp-ph

Boundary-conformal integration for the invariant-imbedding T-matrix method: high-order convergence for faceted particles

Zihua Wu , Yu Xiong This is my paper

Pith reviewed 2026-06-26 13:58 UTC · model grok-4.3

classification ⚛️ physics.optics physics.ao-phphysics.comp-ph
keywords invariant-imbedding T-matrixfaceted particlesboundary-conformal integrationlight scatteringhigh-order convergencehexagonal prismice crystalsstaircasing
0
0 comments X

The pith

A boundary-conformal integration scheme removes geometric non-smoothness and restores spectral plus fourth-order convergence in the invariant-imbedding T-matrix method for faceted particles.

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

The invariant-imbedding T-matrix method loses accuracy on faceted particles because standard quadrature encounters abrupt changes when the integration sphere crosses the particle faces and edges. These changes produce jumps, kinks, and half-integer branches whose origin is purely geometric. The paper introduces a boundary-conformal scheme that replaces the affected integrals with closed-form azimuthal coefficients, splits panels exactly at the analytically known tangency loci, and applies a square-root substitution to cancel the branch points. On a hexagonal prism this makes the azimuthal integral exact while the zenithal and radial integrals recover spectral and fourth-order convergence respectively. The same construction applies to any convex polyhedron because it depends only on the local contact geometry between sphere and particle.

Core claim

The non-smoothness of the dielectric-contrast integral has a single geometric origin in the tangencies of the integration sphere to the faces and edges of the particle. These tangencies generate jumps, kinks, and half-integer branches in all three coordinate directions. A boundary-conformal scheme eliminates them through closed-form azimuthal coefficients, panel splitting at the known tangency loci, and the substitution x maps to x_c plus t squared that absorbs the half-integer branches. For the hexagonal prism the azimuthal integration becomes exact, the zenithal direction recovers spectral convergence, and the radial direction recovers fourth-order convergence; the same scheme extends dire

What carries the argument

Boundary-conformal integration scheme that employs closed-form azimuthal coefficients, panel splitting at analytically known tangency loci, and square-root substitution x maps to x_c plus t squared.

If this is right

  • Azimuthal integration becomes exact for a hexagonal prism.
  • Zenithal integration recovers spectral convergence.
  • Radial integration recovers fourth-order convergence.
  • The scheme extends to any convex polyhedron because it uses only contact geometry.
  • Convergence orders are fixed by local contact geometry and remain independent of size up to kr_max equals 20.

Where Pith is reading between the lines

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

  • The same tangency-based splitting and substitution could be inserted into other volume-integral scattering codes that integrate over faceted domains.
  • Because the required resolution grows with particle size while the order stays fixed, the method supplies a clear scaling rule for atmospheric ice-crystal calculations at visible and infrared wavelengths.
  • The analytic location of tangency loci suggests that the scheme can be combined with adaptive quadrature that automatically refines only near edges and vertices.

Load-bearing premise

The non-smoothness originates solely from the tangencies of the integration sphere to the particle faces and edges.

What would settle it

A numerical test on a hexagonal prism that applies the full scheme yet fails to reach spectral convergence in the zenithal direction or fourth-order convergence in the radial direction.

Figures

Figures reproduced from arXiv: 2606.21063 by Yu Xiong, Zihua Wu.

Figure 1
Figure 1. Figure 1: Geometry of the three IITM quadratures for a regular hexagonal prism (apothem b, circumradius Rc = b/ cos(π/N), half-height h/2). (a) As the integration sphere of radius r sweeps outward it becomes tangent to the particle at three feature types (the vertical faces (r = b), the vertical edges (r = Rc), and the top/bottom faces (r = h/2)), and each tangency type seeds one or more of the non-smooth branches a… view at source ↗
Figure 2
Figure 2. Figure 2: Size dependence of the substituted (t 2 ) zenithal quadrature (hexagonal prism). The convergence is spectral at every size (all three curves reach the same ∼10−10 floor), but the resolution at which it enters the spectral regime recedes with krmax: the collapse moves from ≈ 12 Gauss points per panel at krmax = 2.5 (nmax = 10), through ≈ 16 at krmax = 10 (nmax = 21), to ≈ 24 at krmax = 20 (nmax = 33), as th… view at source ↗
Figure 3
Figure 3. Figure 3: Angular quadrature convergence for the hexagonal prism (radial resolution fixed); the two panels are separate studies with independent axes. Azimuth (abscissa Nφ; ordinate the relative error of the Fourier contrast coefficient): the equidistant-FFT coefficients converge as O(1/Nφ) (the staircasing of a boxcar, Proposition 1), while the closed-form arc coefficients (Lemma 1) are exact (dashed floor at ∼10−1… view at source ↗
Figure 4
Figure 4. Figure 4: Radial convergence (all panels: relative error in Csca versus the radial resolution Nr at fixed, machine-precision angular quadrature). All schemes are compared at equal radial work (Nr degrees of freedom); the Möbius curves labeled “no sub.” are the linear-lift integrator without the t 2 substitution (not a different scheme). The three particles isolate the effect: faceted non-axisymmetric, smooth axisymm… view at source ↗
Figure 5
Figure 5. Figure 5: The substitution locus across shapes and aspect ratios. Each panel is a meridian cross-section with the nested integration-sphere radii at which the shell becomes tangent to a feature. The √ r − rc branch absorbed by the t 2 substitution (red) sits at the inner radius rmin (prolate spheroid, tall cylinder), the outer radius rmax (oblate spheroid), or an interior radius (flat cylinder), and at the vertical-… view at source ↗
Figure 6
Figure 6. Figure 6: Boundary-conformal convergence on the solid hexagonal bullet (a convex ice habit with tilted pyramid faces; a = 1, Lc = 2, hp = 1.5, m = 1.311, krmax = 2.5, nmax = 10), relative error in Csca against the extrapolated reference C ⋆ sca. Zenith (Gauss points per panel): panel splitting alone leaves the tilted pyramid-face square-root branches unabsorbed and converges only algebraically (O(N −3 ), Proposition… view at source ↗
read the original abstract

The invariant-imbedding T-matrix method (IITM) is a standard tool for light scattering by large, sharply faceted, non-axisymmetric particles (atmospheric ice crystals and mineral dust) where the surface-based extended boundary condition method loses accuracy. Its accuracy is limited by "staircasing": the dielectric contrast of a faceted particle is integrated across boundaries that cut the quadrature grid, so standard quadrature converges at low algebraic order. We show that this non-smoothness has a single geometric origin, the tangencies of the integration sphere to the faces and edges of the particle, which produce jumps, kinks, and half-integer branches according to the tangency type, in all three integration directions. A boundary-conformal scheme removes them using closed-form azimuthal coefficients, panel splitting at the analytically known tangency loci, and a square-root substitution $x \mapsto x_c + t^2$ that absorbs the half-integer branches. For a hexagonal prism the azimuthal integration becomes exact and the zenithal and radial directions recover spectral and fourth-order convergence; because the construction depends only on the contact geometry, it extends to any convex polyhedron, demonstrated on the solid hexagonal bullet (a faceted ice habit with tilted faces). The zenithal crossing is a square-root branch rather than a kink, so the established interval-splitting alone gives only $\mathcal{O}(N^{-3})$, while the radial step removes the half-integer edge branch that caps the Riccati recurrence on faceted particles. The convergence orders are fixed by the local contact geometry and verified size-independent up to $k\,r_{\max} = 20$; what grows with size is the resolution needed to reach each asymptotic regime, not the order.

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

0 major / 3 minor

Summary. The paper claims that non-smoothness in the invariant-imbedding T-matrix method for faceted particles originates solely from tangencies of the integration sphere with faces and edges, producing jumps, kinks, and half-integer branches. A boundary-conformal scheme removes these via closed-form azimuthal coefficients, analytic panel splitting at known tangency loci, and the substitution x ↦ x_c + t². For a hexagonal prism this yields exact azimuthal integration, spectral zenithal convergence, and fourth-order radial convergence; the construction extends to any convex polyhedron and is demonstrated on the solid hexagonal bullet, with orders verified size-independent up to kr_max=20.

Significance. If the geometric analysis and quadrature recoveries hold, the work provides a parameter-free route to high-order convergence for scattering by sharply faceted non-axisymmetric particles, directly relevant to atmospheric ice crystals and mineral dust. The size-independent orders and explicit extension to general convex polyhedra are concrete strengths that would improve both accuracy and computational cost relative to staircased quadrature.

minor comments (3)
  1. [Abstract] Abstract and §1: the statement that 'the zenithal crossing is a square-root branch rather than a kink' is load-bearing for the O(N^{-3}) claim under interval-splitting alone; a short explicit contrast (perhaps via the local expansion of the integrand) would strengthen the exposition.
  2. Figure captions and §4: the radial and zenithal convergence plots would benefit from explicit annotation of the asymptotic slopes (e.g., 'slope 4' and 'spectral') to make the order recovery immediately visible without consulting the text.
  3. Notation: the definition of the tangency loci (used for panel splitting) appears in the geometric analysis; a compact table listing the three tangency types, their singularity class, and the corresponding remedy would improve readability for readers outside the immediate subfield.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive summary of our work, the recognition of its relevance to scattering by faceted particles, and the recommendation for minor revision. No specific major comments were provided in the report.

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper derives its boundary-conformal quadrature scheme directly from explicit geometric analysis of sphere-particle tangency loci (producing jumps, kinks, and square-root branches) together with standard substitutions (closed-form azimuthal integrals, analytic panel splits, and the x ↦ x_c + t² change of variable). These steps are constructed from the contact geometry itself and do not rely on any fitted parameter, self-referential definition, or load-bearing self-citation chain; the claimed spectral/4th-order recovery for the hexagonal prism and bullet is presented as a consequence of removing those specific singularities, with orders verified numerically rather than forced by construction. No step reduces the central result to an input by definition.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The paper rests on the geometric assumption that all non-smoothness originates from sphere-particle tangencies and on standard properties of convex polyhedra; no free parameters or new entities are introduced.

axioms (1)
  • domain assumption The non-smoothness has a single geometric origin, the tangencies of the integration sphere to the faces and edges of the particle, which produce jumps, kinks, and half-integer branches according to the tangency type, in all three integration directions.
    This premise directly enables the closed-form coefficients, panel splitting, and substitution; it is stated in the abstract as the origin of the staircasing problem.

pith-pipeline@v0.9.1-grok · 5857 in / 1456 out tokens · 32554 ms · 2026-06-26T13:58:47.649510+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

33 extracted references · 26 canonical work pages

  1. [1]

    P. C. Waterman, Symmetry, unitarity, and geometry in electromagnetic scattering, Physical Review D 3 (1971) 825–839. doi: 10.1103/PhysRe vD.3.825

    work page doi:10.1103/physre 1971

  2. [2]

    M. I. Mishchenko, L. D. Travis, A. A. Lacis, Scattering, Absorption, and Emission of Light by Small Particles, Cambridge University Press, Cambridge, 2002

    2002

  3. [3]

    Ganesh, S

    M. Ganesh, S. C. Hawkins, T-matrix computations for light scattering by penetrable particles with large aspect ratios, Journal of Quantitative Spectroscopy and Radiative Transfer 334 (2025) 109346. doi: 10.1016/ j.jqsrt.2025.109346

    arXiv 2025

  4. [4]

    B. R. Johnson, Invariant imbedding T matrix approach to electromag- netic scattering, Applied Optics 27 (1988) 4861–4873. doi: 10.1364/AO .27.004861

    work page doi:10.1364/ao 1988

  5. [5]

    L. Bi, P. Yang, G. W. Kattawar, M. I. Mishchenko, Efficient implemen- tation of the invariant imbedding T-matrix method and the separation of variables method applied to large nonspherical inhomogeneous parti- cles, Journal of Quantitative Spectroscopy and Radiative Transfer 116 (2013) 169–183. doi: 10.1016/j.jqsrt.2012.11.014

    work page doi:10.1016/j.jqsrt.2012.11.014 2013

  6. [6]

    L. Bi, P. Yang, Accurate simulation of the optical properties of at- mospheric ice crystals with the invariant imbedding T-matrix method, Journal of Quantitative Spectroscopy and Radiative Transfer 138 (2014) 17–35. doi: 10.1016/j.jqsrt.2014.01.013

    work page doi:10.1016/j.jqsrt.2014.01.013 2014

  7. [7]

    B. Sun, L. Bi, P. Yang, M. Kahnert, G. W. Kattawar, Invariant Imbed- ding T-matrix Method for Light Scattering by Nonspherical and Inho- mogeneous Particles, Elsevier, Amsterdam, 2020. doi: 10.1016/C201 8-0-02999-0

    work page doi:10.1016/c201 2020

  8. [8]

    P. Yang, L. Bi, B. A. Baum, K.-N. Liou, G. W. Kattawar, M. I. Mishchenko, B. Cole, Spectrally consistent scattering, absorption, and polarization properties of atmospheric ice crystals at wavelengths from 0.2 to 100 µm, Journal of the Atmospheric Sciences 70 (2013) 330–347. doi:10.1175/JAS-D-12-039.1

    work page doi:10.1175/jas-d-12-039.1 2013

  9. [9]

    Xiong, Z

    Y. Xiong, Z. Wu, TransitionMatrices.jl: a high-efficiency open-source Julia framework for electromagnetic scattering and polarimetric remote 29 sensing of nonspherical atmospheric particles, 2025. SSRN preprint, ht tps://doi.org/10.2139/ssrn.5335549 ; code: https://github.com /JuliaRemoteSensing/TransitionMatrices.jl

    work page doi:10.2139/ssrn.5335549 2025

  10. [10]

    S. Zhai, R. L. Panetta, P. Yang, Improvements in the computational efficiency and convergence of the invariant imbedding T-matrix method for spheroids and hexagonal prisms, Optics Express 27 (2019) A1441– A1457. doi: 10.1364/OE.27.0A1441

    work page doi:10.1364/oe.27.0a1441 2019

  11. [11]

    Doicu, T

    A. Doicu, T. Wriedt, N. Khebbache, An overview of the methods for deriving recurrence relations for T-matrix calculation, Journal of Quantitative Spectroscopy and Radiative Transfer 224 (2019) 289–302. doi:10.1016/j.jqsrt.2018.11.029

    work page doi:10.1016/j.jqsrt.2018.11.029 2019

  12. [12]

    S. Hu, L. Liu, Q. Zeng, T. Gao, F. Zhang, An investigation of the symmetrical properties in the invariant imbedding T-matrix method for the nonspherical particles with symmetrical geometry, Journal of Quantitative Spectroscopy and Radiative Transfer 259 (2021) 107401. doi:10.1016/j.jqsrt.2020.107401

    work page doi:10.1016/j.jqsrt.2020.107401 2021

  13. [13]

    J. Zhao, S. Hu, X. Liu, S. Li, The computational optimization of the invariant imbedding T matrix method for the particles with N-fold symmetry, Remote Sensing 14 (2022) 4061. doi: 10.3390/rs14164061

    work page doi:10.3390/rs14164061 2022

  14. [14]

    S. Hu, S. Li, Q. Zeng, L. Liu, Dimension-variable invariant imbedding (DVIIM) T-matrix computational method for the light scattering simu- lation of atmospheric nonspherical particles, Optics Express 31 (2023) 10052–10069. doi: 10.1364/OE.472809

    work page doi:10.1364/oe.472809 2023

  15. [15]

    B. Sun, C. Gao, L. Bi, R. Spurr, Analytical Jacobians of single scatter- ing optical properties using the invariant imbedding T-matrix method, Optics Express 29 (2021) 9635–9669. doi: 10.1364/OE.421886

    work page doi:10.1364/oe.421886 2021

  16. [16]

    B. Sun, C. Gao, D. Liang, Z. Liu, J. Liu, Capability and convergence of linearized invariant-imbedding T-matrix and physical-geometric optics methods for light scattering, Optics Express 30 (2022) 37769–37785. doi:10.1364/OE.473075

    work page doi:10.1364/oe.473075 2022

  17. [17]

    S. Hu, L. Liu, T. Gao, Q. Zeng, An analysis of the factors influencing the modeling accuracy of the invariant imbedding T-matrix method and the optimal design of the parameter settings for particles with different geometrical and optical properties, Journal of Quantitative 30 Spectroscopy and Radiative Transfer 256 (2020) 107306. doi: 10.1016/ j.jqsrt.2020.107306

    arXiv 2020

  18. [18]

    Zhang, J

    Y. Zhang, J. Ding, P. Yang, R. L. Panetta, Vector spherical wave func- tion truncation in the invariant imbedding T-matrix method, Optics Express 30 (2022) 30020–30037. doi: 10.1364/OE.459648

    work page doi:10.1364/oe.459648 2022

  19. [19]

    Z. Wang, L. Bi, S. Kong, Flexible implementation of the particle shape and internal inhomogeneity in the invariant imbedding T-matrix method, Optics Express 31 (2023) 29427–29444. doi: 10.1364/OE.498 190

    work page doi:10.1364/oe.498 2023

  20. [20]

    M. I. Mishchenko, Comprehensive thematic T-matrix reference database: a 2017–2019 update, Journal of Quantitative Spectroscopy and Radiative Transfer 242 (2020) 106692. doi: 10.1016/j.jqsrt.2019 .106692

    work page doi:10.1016/j.jqsrt.2019 2017

  21. [21]

    M. G. Duffy, Quadrature over a pyramid or cube of integrands with a singularity at a vertex, SIAM Journal on Numerical Analysis 19 (1982) 1260–1262. doi: 10.1137/0719090

    work page doi:10.1137/0719090 1982

  22. [22]

    J. C. F. Telles, A self-adaptive co-ordinate transformation for efficient numerical evaluation of general boundary element integrals, Interna- tional Journal for Numerical Methods in Engineering 24 (1987) 959–973. doi:10.1002/nme.1620240509

    work page doi:10.1002/nme.1620240509 1987

  23. [23]

    P. J. Davis, P. Rabinowitz, Methods of Numerical Integration, 2nd ed., Academic Press, Orlando, 1984

    1984

  24. [24]

    L. N. Trefethen, Approximation Theory and Approximation Practice, SIAM, Philadelphia, 2013

    2013

  25. [25]

    Dieci, T

    L. Dieci, T. Eirola, Preserving monotonicity in the numerical solution of Riccati differential equations, Numerische Mathematik 74 (1996) 35–47. doi: 10.1007/s002110050206

    work page doi:10.1007/s002110050206 1996

  26. [26]

    Schiff, S

    J. Schiff, S. Shnider, A natural approach to the numerical integration of Riccati differential equations, SIAM Journal on Numerical Analysis 36 (1999) 1392–1413. doi: 10.1137/S0036142996307946

    work page doi:10.1137/s0036142996307946 1999

  27. [27]

    Hairer, C

    E. Hairer, C. Lubich, G. Wanner, Geometric Numerical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations, volume 31 of Springer Series in Computational Mathematics , 2nd ed., Springer, Berlin, 2006. doi: 10.1007/3-540-30666-8 . 31

    work page doi:10.1007/3-540-30666-8 2006

  28. [28]

    S. G. Warren, R. E. Brandt, Optical constants of ice from the ultraviolet to the microwave: A revised compilation, Journal of Geophysical Re- search: Atmospheres 113 (2008) D14220. doi: 10.1029/2007JD009744

    work page doi:10.1029/2007jd009744 2008

  29. [29]

    C. G. Schmitt, A. J. Heymsfield, On the occurrence of hollow bullet rosette- and column-shaped ice crystals in midlatitude cirrus, Journal of the Atmospheric Sciences 64 (2007) 4515–4520. doi: 10.1175/2007JA S2317.1

    work page doi:10.1175/2007ja 2007

  30. [30]

    F. M. Kahnert, J. J. Stamnes, K. Stamnes, Application of the extended boundary condition method to particles with sharp edges: a comparison of two surface integration approaches, Applied Optics 40 (2001) 3101–

    2001

  31. [31]

    doi: 10.1364/AO.40.003101

    work page doi:10.1364/ao.40.003101

  32. [32]

    Calogero, Variable Phase Approach to Potential Scattering, vol- ume 35 of Mathematics in Science and Engineering , Academic Press, New York, 1967

    F. Calogero, Variable Phase Approach to Potential Scattering, vol- ume 35 of Mathematics in Science and Engineering , Academic Press, New York, 1967

    1967

  33. [33]

    Smooth sets of fields: A pedagogical introduction

    E. Hairer, S. P. Nørsett, G. Wanner, Solving Ordinary Differential Equa- tions I: Nonstiff Problems, volume 8 of Springer Series in Computa- tional Mathematics , 2nd ed., Springer, Berlin, 1993. doi: 10.1007/97 8-3-540-78862-1 . 32

    work page doi:10.1007/97 1993