pith. sign in

arxiv: 2510.11042 · v3 · submitted 2025-10-13 · ⚛️ physics.optics · physics.comp-ph

Lattice Boltzmann Method for Electromagnetic Wave Scattering

Mohd. Meraj Khan , Sumesh P. Thampi , Anubhab Roy This is my paper

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

classification ⚛️ physics.optics physics.comp-ph
keywords lattice boltzmann methodelectromagnetic scatteringmaxwell equationslorenz-mie solutiondielectric cylinderdielectric spherenumerical wave propagation
0
0 comments X

The pith

The lattice Boltzmann method solves Maxwell's equations for electromagnetic scattering and matches analytical solutions for dielectric cylinders, spheres and hexagons.

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

The paper tests whether the lattice Boltzmann method can serve as a practical time-domain solver for electromagnetic wave scattering problems. It applies the method to planar interfaces, two-dimensional circular cylinders, three-dimensional spheres, and non-circular hexagonal cylinders. In each case the computed scattering intensities agree closely with exact Lorenz-Mie solutions and other reference calculations across a range of size-to-wavelength ratios and dielectric contrasts. The work therefore presents LBM as an explicit, grid-based alternative that is naturally suited to parallel computation.

Core claim

When the lattice Boltzmann method is formulated for Maxwell's equations and run on structured grids, the resulting angular scattering intensities reproduce analytical Lorenz-Mie results for infinitely long circular dielectric cylinders and for dielectric spheres, and also match the Discretized-Mie Formalism for hexagonal cylinders, over a range of size parameters and material contrasts.

What carries the argument

Lattice Boltzmann discretization of Maxwell's equations using explicit streaming and collision steps with relaxation parameters chosen for electromagnetic propagation on a Cartesian grid.

If this is right

  • LBM scattering results remain accurate when the dielectric constant of the scatterer is varied.
  • The same framework extends directly to three-dimensional problems such as dielectric spheres.
  • Non-circular cross-sections such as hexagons produce results consistent with independent semi-analytical methods.
  • The explicit time-stepping structure allows straightforward parallel implementation on regular grids.

Where Pith is reading between the lines

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

  • The grid-based nature of the method may allow easy handling of irregular or time-varying scatterers without major reformulation.
  • Because LBM is already used for fluid flow, the same code base could be extended to coupled light-fluid problems.
  • Stability limits for very high dielectric contrasts or absorbing materials remain to be mapped systematically.

Load-bearing premise

The chosen lattice Boltzmann discretization and relaxation parameters stay stable and accurate for the dielectric contrasts and grid resolutions used in the tested scattering geometries.

What would settle it

A clear mismatch between lattice Boltzmann scattering patterns and the corresponding analytical Lorenz-Mie curves at a higher size-to-wavelength ratio or larger dielectric contrast than those already tested would show the method is not generally reliable.

Figures

Figures reproduced from arXiv: 2510.11042 by Anubhab Roy, Mohd. Meraj Khan, Sumesh P. Thampi.

Figure 1
Figure 1. Figure 1: Schematic of scattering problems studied: (a) planar dielectric interface, (b) circular cylinder, (c) hexagonal cylinder, and (d) sphere. [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Schematic of the near-to-far-field transformation. (a) The fields are first recorded on a fictitious boundary enclosing the scatterer. (b) [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Normalized amplitudes of reflected (black) and transmitted (red) electric fields relative to the incident field at a planar interface under normal incidence. (a) Dependence on dielectric constant εr with µr = 1. (b) Dependence on relative permeability µr with εr = 1. Solid lines represent analytical solutions, and markers denote LBM results. For normal incidence, the amplitude ratios of the reflected and t… view at source ↗
Figure 4
Figure 4. Figure 4: Real part of the total electric field for circular cylinders under plane-wave incidence. The top row shows perfect electrically conducting [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Comparison of the normalized scattering width ( [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Comparison between analytical solutions (black solid lines) and LBM results (red dashed lines) for the normalized scattering width [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Comparison of normalized scattering width ( [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Scattering of a hexagonal dielectric cylinder with [PITH_FULL_IMAGE:figures/full_fig_p012_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Schematic of angular discretization in the DMF for a regular hexagon. The angular domain is sampled at [PITH_FULL_IMAGE:figures/full_fig_p013_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Spherical coordinate system used for scattering from a dielectric sphere. The polar angle [PITH_FULL_IMAGE:figures/full_fig_p014_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Absolute value of the total electric field (incident + scattered) around a dielectric sphere with [PITH_FULL_IMAGE:figures/full_fig_p015_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Comparison of radar cross section (RCS) of a dielectric sphere with [PITH_FULL_IMAGE:figures/full_fig_p016_12.png] view at source ↗
read the original abstract

In this work, the lattice Boltzmann method (LBM) is assessed as a time-domain numerical approach for electromagnetic wave scattering. Owing to its explicit formulation and suitability for parallel computation on structured grids, LBM provides an alternative framework for solving Maxwell's equations. The formulation is first validated using canonical benchmarks, including reflection and refraction at a planar dielectric interface and two-dimensional scattering from infinitely long circular cylinders, where the computed angular scattering intensities are compared with analytical Lorenz-Mie solutions. Additional comparisons are performed for circular cylinders with varying dielectric constants to examine performance across different material contrasts. The framework is then extended to three-dimensional scattering from dielectric spheres, representing the most computationally demanding case considered in this work, and the resulting angular scattering intensities are compared with exact Lorenz-Mie solutions. To further examine performance for non-circular geometries, scattering from an infinitely long hexagonal dielectric cylinder is investigated and benchmarked against results obtained using the Discretized-Mie Formalism. Across all cases, the LBM predictions show close agreement with analytical and semi-analytical reference solutions over a range of size-to-wavelength ratios.

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

1 major / 2 minor

Summary. The manuscript assesses the lattice Boltzmann method (LBM) as a time-domain numerical approach for electromagnetic wave scattering. It validates the formulation against analytical Lorenz-Mie solutions for reflection/refraction at planar dielectric interfaces and 2D scattering from infinitely long circular cylinders, extends comparisons to cylinders with varying dielectric constants, then to 3D scattering from dielectric spheres, and finally benchmarks scattering from an infinitely long hexagonal dielectric cylinder against the Discretized-Mie Formalism. Across all cases the LBM predictions are reported to show close agreement with the reference solutions over a range of size-to-wavelength ratios and material contrasts.

Significance. If the reported agreements hold under the stated conditions, the work is significant as a demonstration that an explicit, grid-based LBM formulation can serve as a stable and accurate alternative to conventional time-domain methods for canonical EM scattering problems. The direct, external validations against exact Lorenz-Mie and Discretized-Mie results for multiple geometries and contrasts provide a clear, falsifiable test of the chosen discretization and relaxation parameters.

major comments (1)
  1. The central claim rests on the stability and accuracy of the specific LBM discretization and relaxation parameters for the dielectric contrasts and grid resolutions employed. While the abstract states close agreement, the manuscript does not appear to tabulate or derive these parameters explicitly (e.g., relaxation time, equilibrium distribution, or collision operator for Maxwell's equations), making it difficult to assess generality beyond the presented benchmarks.
minor comments (2)
  1. Add quantitative error measures (e.g., integrated squared difference or maximum relative error in angular intensity) alongside the qualitative statements of agreement in the results sections for the 2D cylinder and 3D sphere cases.
  2. Clarify boundary-condition implementation at dielectric interfaces and the treatment of the scattered-field formulation, as these details are load-bearing for reproducibility of the reported planar-interface and cylinder results.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and the recommendation of minor revision. We address the single major comment below.

read point-by-point responses

  1. Referee: The central claim rests on the stability and accuracy of the specific LBM discretization and relaxation parameters for the dielectric contrasts and grid resolutions employed. While the abstract states close agreement, the manuscript does not appear to tabulate or derive these parameters explicitly (e.g., relaxation time, equilibrium distribution, or collision operator for Maxwell's equations), making it difficult to assess generality beyond the presented benchmarks.

    Authors: We thank the referee for this observation. Section 2 of the manuscript derives the LBM collision operator for Maxwell's equations, the equilibrium distribution functions, and the relation between the relaxation time and the dielectric permittivity (via the effective speed of light on the lattice). The specific relaxation times are selected to maintain numerical stability while satisfying the Courant-Friedrichs-Lewy condition for each refractive index; these choices are stated in the text accompanying each benchmark. To improve reproducibility and allow readers to assess generality more readily, we will add an explicit table in the revised manuscript that lists the relaxation time, grid resolution, time step, and dielectric contrast for every test case. We do not claim broad generality beyond the validated regimes; the reported agreements with exact Lorenz-Mie and Discretized-Mie solutions across multiple geometries and contrasts serve as the primary evidence that the chosen parameters are appropriate for the problems considered. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The manuscript formulates the lattice Boltzmann method for time-domain Maxwell equations and validates it exclusively against independent external references: exact Lorenz-Mie analytical solutions for planar interfaces, circular cylinders, and spheres, plus the Discretized-Mie Formalism for hexagonal cylinders. These benchmarks are not derived from the LBM discretization parameters, relaxation times, or any fitted quantities internal to the paper; they constitute separate, closed-form or semi-analytical results. No self-definitional steps, fitted-input predictions, or load-bearing self-citations appear in the derivation or validation chain. The central claim of agreement across dielectric contrasts and size-to-wavelength ratios therefore rests on external falsifiability rather than reduction to the method's own inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review; no explicit free parameters, axioms, or invented entities are stated. Typical LBM work implicitly relies on discretization choices and relaxation-time parameters whose values are not reported here.

pith-pipeline@v0.9.0 · 5720 in / 1074 out tokens · 29033 ms · 2026-05-18T08:07:01.984333+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

54 extracted references · 54 canonical work pages · 1 internal anchor

  1. [1]

    A. V . Osipov, S. A. Tretyakov, Modern electromagnetic scattering theory with applications, John Wiley & Sons, 2017

    work page 2017

  2. [2]

    Kerker, The scattering of light and other electromagnetic radiation., Elsevier, 1969

    M. Kerker, The scattering of light and other electromagnetic radiation., Elsevier, 1969

    work page 1969

  3. [3]

    K.-N. Liou, P. Yang, Light scattering by ice crystals: fundamentals and applications, Cambridge University Press, 2016

    work page 2016

  4. [4]

    J. Li, B. E. Carlson, Y . L. Yung, D. Lv, J. Hansen, J. E. Penner, H. Liao, V . Ramaswamy, R. A. Kahn, P. Zhang, Scattering and absorbing aerosols in the climate system, Nature Reviews Earth & Environment 3 (6) (2022) 363–379

    work page 2022

  5. [5]

    M. T. Falconi, F. S. Marzano, Weather Radar Data Processing and Atmospheric Applications: An overview of tools for monitoring clouds and detecting wind shear, IEEE Signal Processing Magazine 36 (4) (2019) 85–97

    work page 2019

  6. [6]

    Ochoa-Rodriguez, L.-P

    S. Ochoa-Rodriguez, L.-P. Wang, P. Willems, C. Onof, A review of radar-rain gauge data merging methods and their potential for urban hydrological applications, Water Resources Research 55 (8) (2019) 6356–6391

    work page 2019

  7. [7]

    Rashidi, S

    A. Rashidi, S. Domínguez-Medina, J. Yan, D. S. Efremenko, A. A. Vasilyeva, A. Doicu, T. Wriedt, C. L. Wirth, Developing scattering morphology resolved total internal reflection microscopy (SMR-TIRM) for orientation detection of colloidal ellipsoids, Langmuir 36 (43) (2020) 13041–13050

    work page 2020

  8. [8]

    W. Shao, T. McCollough, Advances in microwave near-field imaging: Prototypes, systems, and applications, IEEE microwave magazine 21 (5) (2020) 94–119

    work page 2020

  9. [9]

    H. C. van de Hulst, Light scattering by small particles, Courier Corporation, 1981

    work page 1981

  10. [10]

    C. A. Balanis, Advanced engineering electromagnetics, John Wiley & Sons, 2012

    work page 2012

  11. [11]

    J. A. Stratton, Electromagnetic theory, John Wiley & Sons, 2007

    work page 2007

  12. [12]

    F. M. Kahnert, Numerical methods in electromagnetic scattering theory, Journal of Quantitative Spectroscopy and Radiative Transfer 79 (2003) 775–824

    work page 2003

  13. [13]

    M. Kahnert, Numerical solutions of the macroscopic Maxwell equations for scattering by non-spherical particles: a tutorial review, Journal of Quantitative Spectroscopy and Radiative Transfer 178 (2016) 22–37. 18

    work page 2016

  14. [14]

    M. I. Mishchenko, Electromagnetic scattering by particles and particle groups: an introduction, Cambridge University Press, 2014

    work page 2014

  15. [15]

    Atlas, M

    D. Atlas, M. Kerker, W. Hitschfeld, Scattering and attenuation by non-spherical atmospheric particles, Journal of Atmospheric and Terrestrial Physics 3 (2) (1953) 108–119

    work page 1953

  16. [16]

    T. A. Seliga, V . Bringi, Potential use of radar differential reflectivity measurements at orthogonal polarizations for measuring precipitation, Journal of Applied Meteorology and Climatology 15 (1) (1976) 69–76

    work page 1976

  17. [17]

    Matrosov, Prospects for the measurement of ice cloud particle shape and orientation with elliptically polarized radar signals, Radio science 26 (4) (1991) 847–856

    S. Matrosov, Prospects for the measurement of ice cloud particle shape and orientation with elliptically polarized radar signals, Radio science 26 (4) (1991) 847–856

    work page 1991

  18. [18]

    S. Y . Matrosov, Theoretical study of radar polarization parameters obtained from cirrus clouds, Journal of At- mospheric Sciences 48 (8) (1991) 1062–1070

    work page 1991

  19. [19]

    S. Y . Matrosov, Radar reflectivity in snowfall, IEEE transactions on geoscience and remote sensing 30 (3) (1992) 454–461

    work page 1992

  20. [20]

    Barber, C

    P. Barber, C. Yeh, Scattering of electromagnetic waves by arbitrarily shaped dielectric bodies, Applied optics 14 (12) (1975) 2864–2872

    work page 1975

  21. [21]

    Seliga, V

    T. Seliga, V . Bringi, Differential reflectivity and differential phase shift: Applications in radar meteorology, Radio Science 13 (2) (1978) 271–275

    work page 1978

  22. [22]

    C. Yeh, R. Woo, A. Ishimaru, J. Armstrong, Scattering by single ice needles and plates at 30 ghz, Radio Science 17 (06) (1982) 1503–1510

    work page 1982

  23. [23]

    Vivekanandan, W

    J. Vivekanandan, W. Adams, V . Bringi, Rigorous approach to polarimetric radar modeling of hydrometeor ori- entation distributions, Journal of Applied Meteorology (1988-2005) (1991) 1053–1063

    work page 1988

  24. [24]

    Aydin, C

    K. Aydin, C. Tang, Millimeter wave polarimetric scattering from ice crystals, in: 26th International Conference on radar meteorology, Publ by American Meteorological Soc, 1993, pp. 123–125

    work page 1993

  25. [25]

    S. G. O’Brien, G. H. Goedecke, Scattering of millimeter waves by snow crystals and equivalent homogeneous symmetric particles, Applied optics 27 (12) (1988) 2439–2444

    work page 1988

  26. [26]

    Vivekanandan, W

    J. Vivekanandan, W. Adams, Theoretical investigation of multiparameter radar scattering characteristics of ice crystals, in: Preprints, 26th Int. Conf. on Radar Meteorology, Norman, OK, Amer. Meteor. Soc, V ol. 109, 1993, p. 111

    work page 1993

  27. [27]

    C. E. Dungey, C. F. Bohren, Backscattering by nonspherical hydrometeors as calculated by the coupled-dipole method: An application in radar meteorology, Journal of Atmospheric and Oceanic Technology 10 (4) (1993) 526–532

    work page 1993

  28. [28]

    T. L. Schneider, G. L. Stephens, Backscattering by nonspherical ice particles at millimeter wavelengths: a theo- retical study (1994)

    work page 1994

  29. [29]

    Yee, Numerical solution of initial boundary value problems involving maxwell’s equations in isotropic media, IEEE Transactions on antennas and propagation 14 (3) (1966) 302–307

    K. Yee, Numerical solution of initial boundary value problems involving maxwell’s equations in isotropic media, IEEE Transactions on antennas and propagation 14 (3) (1966) 302–307

    work page 1966

  30. [30]

    P. C. Waterman, Symmetry, unitarity, and geometry in electromagnetic scattering, Physical review D 3 (4) (1971) 825

    work page 1971

  31. [31]

    E. M. Purcell, C. R. Pennypacker, Scattering and absorption of light by nonspherical dielectric grains, Astro- physical Journal, V ol. 186, pp. 705-714 (1973) 186 (1973) 705–714

    work page 1973

  32. [32]

    P. P. Silvester, R. L. Ferrari, Finite elements for electrical engineers, Cambridge university press, 1996. 19

    work page 1996

  33. [33]

    J. L. V olakis, A. Chatterjee, L. C. Kempel, Finite element method for electromagnetics, IEEE Press, 1998

    work page 1998

  34. [34]

    Sheng, J.-M

    X.-Q. Sheng, J.-M. Jin, J. Song, C.-C. Lu, W. C. Chew, On the formulation of hybrid finite-element and boundary-integral methods for 3-d scattering, IEEE transactions on antennas and propagation 46 (3) (1998) 303–311

    work page 1998

  35. [35]

    Morgan, C

    M. Morgan, C. Chen, S. Hill, P. Barber, Finite element-boundary integral formulation for electromagnetic scat- tering, Wave Motion 6 (1) (1984) 91–103

    work page 1984

  36. [36]

    X. Yuan, D. R. Lynch, J. W. Strohbehn, Coupling of finite element and moment methods for electromagnetic scattering from inhomogeneous objects, IEEE Transactions on Antennas and Propagation 38 (3) (2002) 386– 393

    work page 2002

  37. [37]

    Zhang, G

    J. Zhang, G. Yan, X. Shi, Lattice boltzmann model for wave propagation, Physical Review E—Statistical, Non- linear, and Soft Matter Physics 80 (2) (2009) 026706

    work page 2009

  38. [38]

    Dellar, Electromagnetic waves in lattice Boltzmann magnetohydrodynamics, Europhysics Letters 90 (5) (2010) 50002

    P. Dellar, Electromagnetic waves in lattice Boltzmann magnetohydrodynamics, Europhysics Letters 90 (5) (2010) 50002

    work page 2010

  39. [39]

    Mendoza, J

    M. Mendoza, J. Munoz, Three-dimensional lattice Boltzmann model for electrodynamics, Physical Review E 82 (5) (2010) 056708

    work page 2010

  40. [40]

    Liu, A multi-energy-level lattice boltzmann model for maxwell’s equations without sources, Journal of Elec- trostatics 69 (6) (2011) 564–570

    Y . Liu, A multi-energy-level lattice boltzmann model for maxwell’s equations without sources, Journal of Elec- trostatics 69 (6) (2011) 564–570

    work page 2011

  41. [41]

    S. M. Hanasoge, S. Succi, S. A. Orszag, Lattice Boltzmann method for electromagnetic wave propagation, EuroPhysics Letters 96 (1) (2011) 14002

    work page 2011

  42. [42]

    Chen, C.-H

    P. Chen, C.-H. Wang, J.-R. Ho, A lattice boltzmann model for electromagnetic waves propagating in a one- dimensional dispersive medium, Computers & Mathematics with Applications 65 (6) (2013) 961–973

    work page 2013

  43. [43]

    Y . Liu, G. Yan, A lattice Boltzmann model for Maxwell’s equations, Applied Mathematical Modelling 38 (5-6) (2014) 1710–1728

    work page 2014

  44. [44]

    Hauser, J

    A. Hauser, J. Verhey, Stable lattice Boltzmann model for Maxwell equations in media, Physical Review E 96 (6) (2017) 063306

    work page 2017

  45. [45]

    Hauser, J

    A. Hauser, J. Verhey, Comparison of the lattice-Boltzmann model with the finite-difference time-domain method for electrodynamics, Physical Review E 99 (3) (2019) 033301

    work page 2019

  46. [46]

    Y . Liu, H. Wang, Simulations of the rectangular wave-guide pattern in the complex maxwell vorticity equations by lattice boltzmann method, Mathematics and Computers in Simulation 173 (2020) 1–15

    work page 2020

  47. [47]

    M. M. Khan, S. P. Thampi, A. Roy, Electromagnetic scattering by curved surfaces and calculation of radiation force: Lattice boltzmann simulations, Journal of Applied Physics 136 (19) (2024)

    work page 2024

  48. [48]

    Krüger, H

    T. Krüger, H. Kusumaatmaja, A. Kuzmin, O. Shardt, G. Silva, E. M. Viggen, The lattice Boltzmann Method, Springer, 2017

    work page 2017

  49. [49]

    Mohamad, Lattice Boltzmann Method, V ol

    A. Mohamad, Lattice Boltzmann Method, V ol. 70, Springer, 2011

    work page 2011

  50. [50]

    M. M. Khan, S. P. Thampi, A. Roy, Radiation forces and torques on janus cylinders, arXiv preprint arXiv:2509.22308 (2025)

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  51. [51]

    J. B. Schneider, Understanding the finite-difference time-domain method,www.eecs.wsu.edu/~schneidj/ ufdtd, accessed: March 20, 2025 (2010)

    work page 2025

  52. [52]

    Taflove, S

    A. Taflove, S. C. Hagness, M. Piket-May, Computational electromagnetics: the finite-difference time-domain method, The Electrical Engineering Handbook 3 (629-670) (2005) 15. 20

    work page 2005

  53. [53]

    D. J. Griffiths, Introduction to electrodynamics, Pearson, 2013

    work page 2013

  54. [54]

    Rother, K

    T. Rother, K. Schmidt, The discretized mie-formalism for electromagnetic scattering-summary, Journal of elec- tromagnetic waves and applications 11 (12) (1997) 1619–1625. 21

    work page 1997