REVIEW 3 major objections 4 minor 9 references
Effects of off-diagonal permittivity terms on polarization singularities in anisotropic grating system
T0 review · 3 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read By rotating the optical axis of a uniaxial material inside a fixed grating, the paper shows that off-diagonal permittivity terms can shift bound states in the continuum, split them into pairs of circularly polarized C points, and even…
desk verdict This is a credible numerical study that usefully catalogs how different off-diagonal permittivity components affect BICs and C points in an anisotropic grating, but the headline TE attribution to ε_yz is not cleanly isolated from simultaneous diagonal changes. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The controlling object is the off-diagonal part of the anisotropic permittivity tensor, obtained by rotating the optic axis of a uniaxial material through polar angle α and azimuth β. In the wave equation inside the material, the eigensolutions show that TE modes (Ey ≠ 0, Ex = Ez = 0) depend on εyz and εzy plus εyy, while TM modes (Ey = 0) involve all three off-diagonal terms εxy, εxz, εyz and the diagonals εxx, εzz. Rotating in the x–y plane generates εxy only; rotating in the y–z plane generates εyz; rotating in the x–z plane (β = 90°, α ≠ 0) generates εxz. The paper tracks the polarization orientation angle φ in momentum space, whose winding number (Eq. 2) gives the topological charge, and uses ellipticity maps to identify C points; the creation and annihilation of C points are governed by conservation of total topological charge.
What would settle it
A full vectorial Bloch-mode calculation of the same grating at, say, β = 10° for the TE1 band, without the pure-TE ansatz, would show whether the eigenmode at the claimed C-point positions really has |Ex|, |Ez| ≪ |Ey|; if not, or if the BICs no longer split into exactly two C points per BIC once hybrid character is included, the attribution of splitting to εyz alone is not solid. Alternatively, an angle-resolved measurement of the far-field polarization for a fabricated sample at α = 0°, β = 10° should show six C points with the predicted positions and handedness; observing fewer or differently handed singularities would falsify the scenario.
Extended reading notes
Core claim
The central finding is a cause–effect map between individual off-diagonal permittivity components and the fate of polarization singularities in an anisotropic grating (period P = 900 nm, uniaxial material with ε1 = 9, ε2 = 9.2). The symmetry-protected BIC at Γ and two accidental BICs on the TE1 band all carry integer topological charges. Introducing εxy by rotating the optic axis in the x–y plane moves the accidental BICs along ky but never splits them. Introducing εyz by rotating in the y–z plane splits every BIC into a pair of C points of the same half-integer charge and opposite handedness; over a full rotation period these pairs annihilate, regenerate, and merge, and with a slightly different grating height the net effect is reversal of the Γ-point BIC's topological charge. For the TM2 band, εyz again splits BICs (along kx, with the opposite splitting direction from TE), while εxz—a term with no counterpart in the TE eigensolution—splits BICs along ky and can reverse the Γ-point charge in its own way. The paper therefore establishes that out-of-plane anisotropy terms split BICs while the in-plane term εxy acts only as a position shifter.
Load-bearing premise
The argument assumes the eigenmodes stay purely TE (Ey ≠ 0, Ex = Ez = 0) or purely TM (Ey = 0) when off-diagonal permittivity terms are present; Appendix B shows this requires conditions like kxky + ω²μ₀ε₀εxy = 0 that cannot hold for every Fourier component of a Bloch mode in a periodic grating, so the modes may actually be hybrid.
Editorial extensions
If this is right
- Rotating the optic axis, rather than fabricating a new lattice, can shift accidental BICs, split them into C points, and merge them back, offering dynamic tuning of polarization singularities in a fixed structure.
- Because the splitting direction depends on which off-diagonal term is active (εyz splits TE BICs along ky; εxz and εyz split TM BICs along ky and kx respectively), one can choose the emission direction of the half-integer-charge C points by the rotation plane.
- Charge reversal at Γ (from +1 to −1 for TE, from −1 to +1 for TM) follows purely from conservation of total topological charge during the C-point evolution, giving a route to flip the vortex charge of a BIC.
- The same mechanism should appear in any periodic slab whose material has a rotatable uniaxial axis, not just the specific grating studied here.
Reading between the lines
- A likely extension: the pure-TE/pure-TM eigensolution assumption used in Appendix B may hold only approximately in periodic Bloch modes; if the modes hybridize, the exact critical angles and positions of the C points could shift, though the qualitative split-and-merge scenario should survive weak coupling.
- A quantitative prediction implicit in the paper is that the BIC-to-C-point splitting distance grows with the magnitude of the active off-diagonal component (maximal at β = 45° for εyz), which could be tested by measuring far-field polarization maps at intermediate angles.
- The same off-diagonal knob could be combined with magneto-optical or nonlinear materials to make the singularities switchable by an external field or pump rather than by a fixed rotation.
- Because εxy only shifts BICs, it could serve as a fine position aligner, while εyz and εxz act as splitters; combining these rotations would allow independent control of BIC location and charge conversion in a single device.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper numerically studies how off-diagonal permittivity tensor elements affect polarization singularities (BICs and C points) in a one-dimensional anisotropic grating. By rotating the optical axis of a uniaxial material, the authors vary the permittivity tensor components and track the evolution of three BICs in both TE and TM bands. They report that in-plane anisotropy (ε_xy) only shifts the accidental BICs, while out-of-plane anisotropy (ε_yz for TE; ε_xz and ε_yz for TM) splits BICs into pairs of C points, with topological charge conservation governing the creation, annihilation, and merging of these singularities. They also observe charge reversal at the Γ point in certain parameter ranges.
Significance. If the central claim holds, the work identifies off-diagonal permittivity terms as a new degree of freedom for controlling polarization singularities without altering the grating geometry, which could be relevant for chiral light manipulation and dynamic modulation of singular optics. The numerical evidence is rich: the reported polarization maps, Q-factor distributions, ellipticity plots, and evolutionary trajectories are consistent with known physics and topological charge conservation. The paper also includes a useful comparative analysis of TE and TM modes. However, the causal attribution of splitting to specific off-diagonal components is not rigorously established, as detailed in the major comments. The manuscript is a computational study with no experimentally falsifiable predictions, but it does provide a clear, reproducible parameter sweep.
major comments (3)
- [§2.1, Fig. 2 and Eq. (1)] The attribution of BIC splitting to the off-diagonal term ε_yz is not isolated from simultaneous changes in the diagonal components. When rotating β with α=0, Eq. (1) shows that ε_yy = ε1 + Δ cos²β and ε_zz = ε1 + Δ sin²β change alongside ε_yz = -Δ sinβ cosβ. No control calculation is reported that keeps ε_yy and ε_zz fixed while independently varying ε_yz, so the observed splitting could in principle be driven by the diagonal-index change rather than by ε_yz. This ambiguity is load-bearing for the paper's central claim that off-diagonal terms are the new control knob; a control simulation with an artificially constructed tensor (e.g., fixed diagonals and nonzero ε_yz) would resolve it.
- [Appendix B, Eq. (A8)] The plane-wave derivation of pure TE modes is not applicable to the periodic grating structure. Equation (A8) requires kx ky + ω² μ0 ε0 ε_xy = 0 and ky kz + ω² μ0 ε0 ε_zy = 0 for a single plane wave, but a Bloch mode in a periodic grating contains many Fourier components with kx = k0x + nG. These conditions cannot hold simultaneously for all Fourier orders at finite ky, so the modes treated in the numerical simulations are at best quasi-TE. Consequently, the statement that the eigenfields depend on ε_yz (or ε_xz for TM) is not a rigorous derivation for the periodic structure and does not by itself justify the attribution of the splitting to specific off-diagonal components.
- [§2.2, Fig. 5(b)] The independent control of ε_xz, ε_xx, and ε_zz shown in Fig. 5(b) is not actually independent: varying α at fixed β=90° changes these components simultaneously, as Eq. (2) (permittivity tensor) makes explicit. The claim that 'nonzero ε_xz will result in the splitting of each BIC along ky direction' is inferred from curves at different α where the diagonal components also differ. Without a control that varies ε_xz alone (e.g., by artificially setting the off-diagonal term while keeping diagonals constant), the conclusion that ε_xz is the cause of the splitting is not fully supported.
minor comments (4)
- [Throughout] Many subscripts are corrupted or missing in the text, e.g., 'εy?', 'εx?', 'εz?' in Section 2 and Appendix B. These should be corrected to the proper tensor indices (ε_yz, ε_xz, ε_xy, ε_yy, ε_zz) to avoid ambiguity.
- [§2.1, Fig. 1(c)] The sign convention for topological charge q is not defined before Eq. (2); consider stating whether q is defined by the winding of the polarization orientation angle with a clockwise or counterclockwise path.
- [§2.1, Fig. 2 caption] The caption says 'blue and red solid dots represent left-hand (LH) and right-hand (RH) C point' but the handedness definition based on the polarization ellipses is not given in the text; adding a sentence defining handedness from the propagation direction would improve clarity.
- [§2.2, Fig. 4(b)] The ellipticity χ = sin(2δ) is introduced, but the relation between δ and the polarization ellipse axes is not fully specified; a brief definition would help readers reproduce the plots.
Circularity Check
No significant circularity: the BIC-to-C-point evolutions are emergent numerical results from Maxwell simulations, not fitted parameters renamed as predictions or conclusions forced by self-citations.
full rationale
The paper's central results are obtained by numerically solving Maxwell's equations for a fixed grating geometry and a specified permittivity tensor, and the observed shifts, splittings, annihilations, and mergings of BICs and C points are emergent features of the computed far-field polarization maps. The initial parameters (epsilon_1=9, epsilon_2=9.2, d=1.1P, rho=0.46P) are hand-picked so that three BICs exist at the starting point, but the subsequent evolution is not imposed by any fitting procedure or by an equation whose output is defined to match the singularity positions. The analytic statements in Appendix B are simple algebraic consequences of the wave equation under the pure-TE and pure-TM ansatz; they are used only to motivate which tensor components may matter and are not invoked as a uniqueness theorem or as a fitted predictor. Self-citations (e.g., Ref. [38]) appear only in the background discussion of prior work on anisotropic cavities and half vortices and are not load-bearing for the present claims. The nearest concern is that rotating the optical axis changes diagonal permittivity components simultaneously with the off-diagonal ones, so the causal attribution to epsilon_yz is not fully isolated by a control calculation; however, this is a support/confound issue rather than a circular reduction. No fitted parameter is renamed as a prediction, and no derived quantity is equivalent to an input by construction.
Assumptions & free parameters
free parameters (4)
- epsilon_1, epsilon_2 (initial permittivities) =
9, 9.2
- Grating height d =
1.1P for TE; 1.09P in Fig 3; 1.33P for TM
- Filling ratio rho =
0.46P
- Period P =
900 nm
assumptions (5)
- standard math Maxwell's equations in anisotropic media (curl E = i*omega*mu0*H, curl H = -i*omega*epsilon0*epsilon_A*E)
- standard math Bloch's theorem for the periodic grating
- standard math Topological charge of a polarization singularity is defined by winding of the polarization orientation angle (Eq. 2)
- domain assumption The uniaxial permittivity tensor for arbitrary OA orientation is obtained by coordinate rotation (Eq. 1)
- domain assumption The numerical simulations accurately resolve the band structure and far-field polarization
Cite this review
Pith. "Pith review of Effects of off-diagonal permittivity terms on polarization singularities in anisotropic grating system." pith.science (2026). https://pith.science/paper/TOMZFXXF
@misc{pith2026250518537,
author = {Pith},
title = {Pith review of: Effects of off-diagonal permittivity terms on polarization singularities in anisotropic grating system},
year = {2026},
howpublished = {\url{https://pith.science/paper/TOMZFXXF}},
note = {Machine review of arXiv:2505.18537}
}
read the original abstract
The evolutions of polarization singularities, including bound states in the continuum (BICs) and circularly polarized states (C points), are usually realized by tuning the geometric parameters of photonic crystal slabs. Here, we use the off-diagonal terms of permittivity tensor to manipulate polarization singularities without breaking the structural symmetry in an anisotropic grating system. By controlling the optical axis of anisotropic media, BICs can be shifted to different positions or split into C points, meanwhile, the creation and annihilation of multiple C points are also observed during the evolution process for both TE and TM modes, respectively. Remarkably, two different splitting directions of BICs can be achieved by tuning the off-diagonal terms of permittivity tensor for the two modes. This work illustrates the important role of off-diagonal terms on the far-field polarization singularities and provide an alternative way to precisely manipulate optical singularities
Reference graph
Works this paper leans on
-
[1]
V.N.Kolobanov,I.A.Kamenskikh,V.V.Mikhailin,etal.,Opticalandluminescentpropertiesofanisotropic tungstatecrystals,Nucl.Instrum.MethodsPhys.Res.A486496-503(2002)
work page 2002
-
[2]
N.C.Passler,X.Ni,G.Hu,etal.,Hyperbolicshearpolaritonsinlow-symmetrycrystals,Nature602595- 600(2022)
work page 2022
-
[3]
Z.Wang,Y.D.Chong,J.D.Joannopoulos,etal.,Reflection-freeone-wayedgemodesinagyromagneticphotonic crystal,Phys.Rev.Lett.100013905(2008)
work page 2008
-
[4]
J.Gomis-Bresco,D.ArtigasandL.Torner,TransitionfromDiracpointstoexceptionalpointsinanisotropic waveguides,Phys.Rev.Res.1(3)033010(2019)
work page 2019
-
[5]
P.Pujol-ClosaandD.Artigas,Diracpointsembeddedinthecontinuum,Phys.Rev.B108205106(2023)
work page 2023
-
[6]
E.M.Renzi,E.Galiffi,X.Ni,etal.,HyperbolicShearMetasurfaces,Phys.Rev.Lett.132(26)263803(2024)
work page 2024
-
[7]
G.Hu,W.Ma,D.Hu,etal.,Real-spacenanoimagingofhyperbolicshearpolaritonsinamonocliniccrystal,Nat. Nanotechnol.1864-70(2023)
work page 2023
-
[8]
W.Lv,H.Qin,Z.Su,etal.,RobustgenerationofintrinsicCpointswithmagneto-opticalboundstatesinthe continuum,Sci.Adv.10eads0157(2024)
work page 2024
Show all 9 references
-
[9]
[10]Q.A.Tu,H.Zhou,D.Zhao,etal.,Magneticallytunableboundstatesinthecontinuumwitharbitrarypolarization andintrinsicchirality,PhotonicsRes.12(12):2972(2024)
J.Li,Z.Wu,D.Zhang,etal.,Nonreciprocaltoroidaldipoleresonanceandone-wayquasi-boundstateinthe continuum,Opt.Lett.491313-1316(2024). [10]Q.A.Tu,H.Zhou,D.Zhao,etal.,Magneticallytunableboundstatesinthecontinuumwitharbitrarypolarization andintrinsicchirality,PhotonicsRes.12(12):2972...
2024
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.