The Fourier modal method for gratings with bi-anisotropic materials
Pith reviewed 2026-05-18 08:51 UTC · model grok-4.3
The pith
The Fourier modal method extended to bi-anisotropic materials converges faster with factorization rules even for large chirality.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The authors formulate the Fourier modal method for gratings with bi-anisotropic materials by deriving Fourier tensors of the material parameters for arbitrary 3x3 magneto-electric coupling tensors, and demonstrate that applying generalized factorization rules improves convergence without sacrificing the reduction to standard cases when coupling vanishes.
What carries the argument
Generalized Fourier factorization rules for the Fourier tensors of the macroscopic material parameters carrying the magneto-electric coupling.
If this is right
- The method enables rigorous simulations of periodic structures containing chiral materials.
- Explicit expressions allow direct implementation in existing Fourier modal method codes.
- Improved convergence holds for large macroscopic chirality coefficients.
- Non-reciprocal materials can be treated within the same framework.
Where Pith is reading between the lines
- This formulation might extend to other numerical methods for Maxwell equations in periodic media with similar tensorial responses.
- Applications could include designing optical isolators or polarization converters using bi-anisotropic gratings.
- Further work could test the method on structures with spatially varying chirality.
Load-bearing premise
The numerical stability of the standard Fourier modal expansion continues to hold when the material response is generalized to arbitrary three-by-three magneto-electric tensors.
What would settle it
Computing the convergence rate for a simple bi-anisotropic grating with large chirality coefficient using both schemes and finding no improvement or worse performance with the factorization rules would disprove the central result.
Figures
read the original abstract
We report an advanced formulation of the Fourier modal method developed for two-dimensionally periodic multilayered structures containing materials with non-zero macroscopic magneto-electric coefficients (also known as coefficients of chirality and bi-anisotropy) represented as arbitrary 3 by 3 tensors. We consider two numerical schemes for this formulation: with and without generalized Fourier factorization rules. For both schemes, we provide explicit expressions for the Fourier tensors of macroscopic material parameters and demonstrate that, in the absence of magneto-electric coupling, they reduce to the conventional factorization rules. We show that the scheme employing factorization rules facilitates improved convergence, even when the macroscopic chirality coefficient is large. The described formulation represents a fast and rigorous technique for theoretical studies of periodic structures with chiral, bi-anisotropic, or non-reciprocal materials in the widely used framework of the Fourier modal method.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports an advanced formulation of the Fourier modal method for two-dimensionally periodic multilayered structures with bi-anisotropic materials characterized by arbitrary 3x3 magneto-electric tensors. It describes two schemes (with and without generalized Fourier factorization rules), provides explicit expressions for the Fourier tensors of the material parameters, demonstrates reduction to conventional rules when magneto-electric coupling is absent, and claims that the factorization scheme leads to improved convergence even for large chirality coefficients. The work positions this as a fast and rigorous technique for theoretical studies of periodic structures involving chiral, bi-anisotropic, or non-reciprocal materials.
Significance. Should the formulation prove valid and the convergence improvement hold under numerical testing, this extension of the Fourier modal method would be significant for the optics community. It allows modeling of a wider range of materials, including those with chirality and non-reciprocity, within an established and computationally efficient framework. The explicit tensor expressions and the reduction to standard cases when coupling vanishes are particularly useful for implementation and validation.
major comments (2)
- [Formulation] The manuscript provides explicit Fourier expressions for the material tensors but does not re-derive the modal eigenvalue problem or the interface conditions from the coupled Maxwell equations for the generalized 3x3 magneto-electric tensors. This is load-bearing for the central claim, as it leaves unresolved whether off-diagonal or non-reciprocal terms introduce additional singularities, alter the Toeplitz structure, or affect the numerical stability of the standard FMM assumptions.
- [Numerical results / Convergence analysis] The claim that the scheme employing factorization rules facilitates improved convergence even when the macroscopic chirality coefficient is large is asserted without accompanying numerical benchmarks, error analysis, or supporting data. This weakens the evidence for the practical advantage of the factorization approach in the bi-anisotropic case.
minor comments (2)
- [Abstract] The abstract mentions 'we show that' the factorization facilitates improved convergence, but without reference to specific figures or tables demonstrating this.
- [References] Ensure that foundational papers on the Fourier modal method and prior extensions to anisotropic materials are cited for context.
Simulated Author's Rebuttal
We thank the referee for their thorough review and insightful comments on our manuscript. We address each of the major comments below and outline the revisions we will make to strengthen the paper.
read point-by-point responses
-
Referee: [Formulation] The manuscript provides explicit Fourier expressions for the material tensors but does not re-derive the modal eigenvalue problem or the interface conditions from the coupled Maxwell equations for the generalized 3x3 magneto-electric tensors. This is load-bearing for the central claim, as it leaves unresolved whether off-diagonal or non-reciprocal terms introduce additional singularities, alter the Toeplitz structure, or affect the numerical stability of the standard FMM assumptions.
Authors: We appreciate this observation. The formulation builds upon the standard Fourier modal method by generalizing the material tensors to include arbitrary 3x3 magneto-electric coupling. The modal eigenvalue problem is obtained by substituting the Fourier expansions of the fields and the material parameters into Maxwell's equations, leading to a system where the constitutive relations are incorporated via the Fourier coefficients of the tensors. We will revise the manuscript to include an explicit derivation of the eigenvalue problem and interface conditions for the bi-anisotropic case. This will demonstrate that the Toeplitz structure is preserved (as the Fourier coefficients form block matrices in the same manner), and that no additional singularities are introduced beyond those mitigated by the factorization rules. The numerical stability assumptions of the standard FMM remain valid as the approach is a direct extension. revision: yes
-
Referee: [Numerical results / Convergence analysis] The claim that the scheme employing factorization rules facilitates improved convergence even when the macroscopic chirality coefficient is large is asserted without accompanying numerical benchmarks, error analysis, or supporting data. This weakens the evidence for the practical advantage of the factorization approach in the bi-anisotropic case.
Authors: We agree that providing numerical evidence is crucial to substantiate the convergence improvement claim. In the revised version, we will add a dedicated section with numerical benchmarks. These will include convergence studies for structures with significant chirality parameters, comparing the two schemes, along with error metrics relative to known analytical or highly converged solutions. This will provide concrete data supporting the practical advantages of the generalized factorization rules. revision: yes
Circularity Check
No significant circularity; explicit tensor expressions extend standard FMM with consistency reduction to known rules
full rationale
The paper derives explicit Fourier expressions for arbitrary 3x3 magneto-electric tensors within the established Fourier modal method framework and verifies that these expressions reduce to conventional factorization rules when magneto-electric coupling is absent. This reduction serves as a consistency check rather than a definitional tautology or fitted-input prediction. Convergence benefits for large chirality are presented as numerical outcomes, not forced by construction. No load-bearing self-citations, uniqueness theorems from prior author work, or ansatzes smuggled via citation are indicated in the derivation chain. The central formulation adds new explicit tensor handling to a pre-existing, externally validated method without reducing its claims to its own inputs.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Fourier modal expansion remains valid and numerically stable for media with arbitrary 3x3 magneto-electric tensors.
Reference graph
Works this paper leans on
-
[1]
M. Moharam and T. Gaylord, Rigorous coupled-wave analysis of planar-grating diffraction., JOSA, 71(7), 811-
-
[2]
S. Tikhodeev, A. Y. E.A., M. N. Gippius, and T. Ishihara, Quasiguided modes and optical properties of photonic crystal slabs., Physical Review B, 66(4), 045102. (2002)
work page 2002
-
[3]
I. M. Fradkin, S. A. Dyakov, and N. A. Gippius, Nanopar- ticle lattices with bases: Fourier modal method and dipole approximation, Physical Review B102, 045432 (2020)
work page 2020
-
[4]
N. S. Salakhova, I. M. Fradkin, S. A. Dyakov, and N. A. Gippius, Fourier modal method for moir´ e lattices, Phys- ical Review B104, 085424 (2021)
work page 2021
-
[5]
W. C. Johnson, Secondary structure of proteins through circular dichroism spectroscopy, Annual Review of Bio- physics and Biophysical Chemistry17, 145 (1988)
work page 1988
- [6]
-
[7]
Y. Tang and A. E. Cohen, Enhanced enantioselectivity in excitation of chiral molecules by superchiral light, Science 332, 333 (2011)
work page 2011
- [8]
-
[9]
I. Hodgkinson and Q. h. Wu, Inorganic chiral optical ma- terials, Advanced materials13, 889 (2001)
work page 2001
- [10]
-
[11]
I. Lindell, A. Sihvola, S. Tretyakov, and A. J. Viitanen, Electromagnetic waves in chiral and bi-isotropic media (Artech House, 1994)
work page 1994
-
[12]
S. Almousa, T. Weiss, and E. Muljarov, Employing quasidegenerate optical modes for chiral sensing, Physi- cal Review B109, L041410 (2024)
work page 2024
-
[13]
F. Nutskii, E. Barredo-Alamilla, and M. A. Gorlach, Complex-valued tellegen response, Applied Physics Let- ters126(2025)
work page 2025
-
[14]
S. A. Dyakov, N. S. Salakhova, A. V. Ignatov, I. M. Frad- kin, V. P. Panov, J.-K. Song, and N. A. Gippius, Chiral light in twisted fabry–p´ erot cavities, Advanced Optical Materials12, 2302502 (2024)
work page 2024
-
[15]
L. Li, Use of fourier series in the analysis of discontinu- ous periodic structures, Journal of the Optical Society of America A, 13(9), 1870. (1996)
work page 1996
-
[16]
L. . Li, Reformulation of the fourier modal method for surface-relief gratings made with anisotropic materials, Journal of Modern Optics, 45(7), 1313–1334. (1998)
work page 1998
-
[17]
L. Li, Fourier modal method for crossed anisotropic grat- ings with arbitrary permittivity and permeability ten- sors, Journal of Optics A: Pure and Applied Optics, 5(4), 345–355. (2003)
work page 2003
-
[18]
R. Redheffer, Inequalities for a matrix riccati equation, Journal of Mathematics and Mechanics , 349 (1959)
work page 1959
-
[19]
R. C. Rumpf, Improved formulation of scattering ma- trices for semi-analytical methods that is consistent with convention, Progress In Electromagnetics Research B35, 241 (2011)
work page 2011
-
[20]
B. E. Saleh and M. C. Teich,Fundamentals of photonics, 2 volume set(john Wiley & sons, 2019)
work page 2019
-
[21]
Simovski,Composite Media with Weak Spatial Disper- sion(Jenny Stanford Publishing, 2018)
C. Simovski,Composite Media with Weak Spatial Disper- sion(Jenny Stanford Publishing, 2018)
work page 2018
-
[22]
G.-G. Liu, S. Mandal, X. Xi, Q. Wang, C. Devescovi, A. Morales-P´ erez, Z. Wang, L. Yang, R. Banerjee, Y. Long,et al., Photonic axion insulator, Science387, 162 (2025)
work page 2025
-
[23]
Q. Yang, X. Wen, Z. Li, O. You, and S. Zhang, Gigantic tellegen responses in metamaterials, Nature Communica- tions16, 151 (2025)
work page 2025
-
[24]
S. Safaei Jazi, I. Faniayeu, R. Cichelero, D. C. Tzarouchis, M. M. Asgari, A. Dmitriev, S. Fan, and V. Asadchy, Optical tellegen metamaterial with sponta- neous magnetization, Nature Communications15, 1293 (2024)
work page 2024
-
[25]
L. Shaposhnikov, M. Mazanov, D. A. Bobylev, F. Wilczek, and M. A. Gorlach, Emergent axion response in multilayered metamaterials, Physical Review B108, 115101 (2023)
work page 2023
-
[26]
T. Z. Seidov and M. A. Gorlach, Unbounded tellegen re- sponse in media with multiple resonances, Physical Re- view A111, 033521 (2025)
work page 2025
-
[27]
W. Brown Jr, R. Hornreich, and S. Shtrikman, Upper bound on the magnetoelectric susceptibility, Physical Re- view168, 574 (1968)
work page 1968
-
[28]
J. Lekner, Optical properties of isotropic chiral media, Pure and Applied Optics: Journal of the European Op- tical Society Part A5, 417 (1996)
work page 1996
-
[29]
Y. Tang and A. E. Cohen, Optical chirality and its inter- action with matter, Physical Review Letters104, 163901 (2010)
work page 2010
- [30]
- [31]
- [32]
-
[33]
J. Petersen, J. Volz, and A. Rauschenbeutel, Chiral nanophotonic waveguide interface based on spin-orbit in- teraction of light, Science346, 67 (2014)
work page 2014
-
[34]
D. G. Baranov, C. Sch¨ afer, and M. V. Gorkunov, To- ward molecular chiral polaritons, ACS Photonics10, 2440 (2023)
work page 2023
-
[35]
W. Chen, Z. Wang, and M. Gorkunov, Uncovering max- imum chirality in resonant nanostructures, Nano Lett., 24, 9643-9649 (2024)
work page 2024
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.