{"id":"6804106b-bc6c-490f-b484-999b932e6e57","arxiv_id":"2505.18537","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"By rotating the optical axis of a uniaxial material, the authors shift or split BICs into C points, with different splitting directions for TE and TM modes.","lead":"This paper shows that rotating the optical axis of an anisotropic material inside a grating can shift or split polarization singularities, such as bound states in the continuum and circularly polarized states, without changing the grating geometry. It offers a new material-based tuning knob for controlling chiral light and vortex beams in photonic devices.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Attribution of BIC splitting to ε_yz is not isolated from simultaneous diagonal-permittivity changes and relies on pure-TE conditions that periodic Bloch modes cannot satisfy.","rationale":"The reader identified the pure-TE/TM assumption as the weakest point. I agree that this is a real limitation, but the more consequential version is that the paper's component-specific causal claim is not isolated. Rotating the optical axis changes diagonal and off-diagonal tensor components together; the reported numerical experiments do not decompose these effects for the TE ε_yz case. The Appendix B derivation only applies to a homogeneous plane wave, not to a Bloch mode in a periodic grating, so it cannot justify the attribution either. The central observations may still be correct, and the TM case includes some independent component sweeps in Fig. 5(b), which is partial support. But the TE splitting claim needs a control simulation to be conclusive. Since the reader already returned CONDITIONAL, my concern reinforces that verdict rather than changing it.","tokens_in":10975,"tokens_out":7602,"duration_ms":75069,"concrete_test":"Run the same RCWA/mode solver for the TE1 band at α=0°, β=10° with two modified permittivity tensors for material A: (a) diag(ε_xx=9, ε_yy=ε1+Δ cos²β, ε_zz=ε1+Δ sin²β) with ε_yz set to zero, and (b) the baseline diagonal values ε_yy=9.2, ε_zz=9 kept fixed while ε_yz is set to its β=10° value. Track the three BIC positions and Q factors in each case. If case (a) also splits the BICs, the attribution to ε_yz fails; if case (b) does not split them, diagonal changes are responsible. Additionally record max(|Ex|,|Ez|)/|Ey| inside material A at the off-Γ BIC; if this ratio is not much smaller than 1, the Appendix B pure-TE conditions (A8) are violated for the actual Bloch mode.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that breaking out-of-plane anisotropy via ε_yz (TE) or ε_xz/ε_yz (TM) splits BICs into C points, while in-plane ε_xy only shifts them. The load-bearing step is the causal attribution to specific off-diagonal components. For the TE case, rotating β with α=0 in Eq. (1) changes ε_yz, but it also changes the diagonal components ε_yy=ε1+Δ cos²β and ε_zz=ε1+Δ sin²β. No control calculation is reported that isolates ε_yz while holding ε_yy and ε_zz fixed, so the observed BIC splitting could in principle be driven by the diagonal-index change rather than by ε_yz. The analytic support in Appendix B does not remove this ambiguity. 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, so these conditions cannot hold simultaneously for all n at finite ky. Thus the claimed 'pure TE mode' is at most a quasi-TE mode, and the statement that the eigensolutions depend on ε_yz is not a rigorous derivation for the periodic structure. If the diagonal components were the actual cause, the paper's central conclusion about off-diagonal control of BIC splitting would not be supported.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":11284,"tokens_out":3561,"duration_ms":30239,"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":[{"comment":"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.","section":"§2.1, Fig. 2 and Eq. (1)"},{"comment":"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.","section":"Appendix B, Eq. (A8)"},{"comment":"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.","section":"§2.2, Fig. 5(b)"}],"minor_comments":[{"comment":"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.","section":"Throughout"},{"comment":"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.","section":"§2.1, Fig. 1(c)"},{"comment":"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.","section":"§2.1, Fig. 2 caption"},{"comment":"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.","section":"§2.2, Fig. 4(b)"}],"recommendation":"major_revision","confidential_remarks":"The paper is a purely numerical study with a valuable qualitative finding, but the absence of control simulations isolating off-diagonal terms from diagonal changes is a serious gap in the causal argument. The plane-wave Appendix B does not support the periodic-structure claims. The authors should be asked to perform control calculations (e.g., fixed diagonals with artificially varied ε_yz or ε_xz) and to provide a symmetry or Bloch-mode argument for why out-of-plane anisotropy splits BICs while in-plane anisotropy does not. If these can be supplied, the paper would be a solid contribution; as it stands, the central attribution is not rigorously established."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a solid simulation paper with a genuinely useful new map: which off-diagonal permittivity component does what to polarization singularities. For TE modes, ε_xy only shifts accidental BICs; ε_yz splits them into C points. For TM modes, ε_yz and ε_xz both split, and the splitting directions differ from TE. They also show charge reversal at Γ and C-point creation/annihilation over a full rotation period. The numerical evidence—Q factors, polarization maps, ellipticity distributions—is consistent and qualitatively convincing.\n\nThe main soft spot is the attribution of the TE-mode splitting to ε_yz. Rotating β with α=0 changes ε_yy and ε_zz as well as ε_yz, and the paper says “combined with the effects of ε_yy” without showing the control. I think the conclusion is probably right, because the diagonal terms preserve the mirror symmetry that protects the BICs, so they would shift rather than split. But that argument should be made explicitly, not left to the reader. A simple control run with only ε_yy varied would close the gap.\n\nThe Appendix B wave-equation analysis is weaker than it looks. It assumes a single plane wave with constant k-vector components, but a Bloch mode in a periodic grating has many Fourier components, so the “pure TE” conditions k_x k_y + ω²μ₀ε₀ε_xy = 0 and k_y k_z + ... = 0 cannot hold for all harmonics. The paper should be honest that the supported statement is quasi-TE in a periodic structure, not rigorous pure-TE. This does not undermine the numerics, but it does mean the analytic narrative is overreach.\n\nMinor reproducibility issue: no code or data deposited; “available on request” is a real limitation for a simulation-only paper. Also, the abstract's “without breaking structural symmetry” is fine—they are breaking material symmetry, not lattice symmetry—so no overclaim there.\n\nWho benefits: researchers working on BIC-to-C-point transitions, topological polarization singularities, chiral emission, or anisotropic photonic crystal slabs. A serious referee would get useful value from this paper. I would send it to review, asking for the ε_yy control calculation and a tone-down of the analytic claims in Appendix B.","headline":"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.","tokens_in":11806,"tokens_out":3367,"would_cite":true,"duration_ms":30967,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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…","keywords":["polarization singularities","bound states in the continuum","C points","anisotropic grating","off-diagonal permittivity tensor","optical axis rotation","topological charge","momentum-space polarization"],"falsifier":"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.","tokens_in":10796,"feed_emoji":"🌀","tokens_out":7878,"duration_ms":52166,"temperature":0.7,"pith_summary":"The paper claims that off-diagonal permittivity terms—introduced simply by rotating the optical axis of a uniaxial material inside a grating—can control polarization singularities without touching the grating geometry. For TE modes, rotating in the x–y plane (εxy ≠ 0) only shifts accidental bound states in the continuum (BICs), while rotating out of that plane (εyz ≠ 0) splits each BIC into a pair of C points with half-integer topological charges. For TM modes, εxz and εyz cause splitting in different directions, and both TE and TM evolutions show creation, annihilation, and merging of C points with conserved total topological charge. If correct, this offers a post-fabrication tuning knob for chiral emission, vortex beams, and other singular-optics applications.","feed_headline":"Rotating the optical axis splits BICs into C points","feed_subtitle":"A fixed grating can move, split, and merge polarization singularities just by rotating a uniaxial crystal's axis.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Defines the topological charge of polarization singularities via the winding of the polarization orientation angle, the diagnostic used throughout the paper.","marker":"[12]"},{"why":"Shows that BICs spawn circularly polarized states (C points) upon symmetry breaking, the phenomenon this paper extends to off-diagonal permittivity terms.","marker":"[25]"},{"why":"Provides the generation-and-annihilation framework for topologically protected BICs and C points under symmetry breaking, used to interpret the observed splitting and merging.","marker":"[26]"},{"why":"Documents the dynamics of topological polarization singularities in momentum space, supplying the evolution scenario the paper reproduces with a different control knob.","marker":"[27]"},{"why":"Demonstrates anisotropy-induced bound states in the continuum in waveguides, establishing that anisotropy alone can create BICs.","marker":"[35]"},{"why":"Shows that rotating the optical axis of anisotropic media tunes topological half vortices in photonic crystal slabs, a direct predecessor for using optical-axis rotation to control singularities.","marker":"[38]"}],"fun_headline_variants":["Off-diagonal permittivity splits BICs into C points","Rotating the optical axis moves and splits singularities","Out-of-plane terms split BICs; in-plane shifts them","Anisotropic grating: off-diagonal terms steer singularities","Tune permittivity off-diagonals to split or shift BICs"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Off-diagonal permittivity splits BICs into C points","Rotating the optical axis moves and splits singularities","Out-of-plane terms split BICs; in-plane shifts them","Anisotropic grating: off-diagonal terms steer singularities","Tune permittivity off-diagonals to split or shift BICs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000767,"raw_usage":{"total_tokens":3407,"prompt_tokens":959,"completion_tokens":2448,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":575,"completion_tokens_details":{"reasoning_tokens":2361}},"tokens_in":575,"tokens_out":2448,"duration_ms":18614,"temperature":1.0,"reasoning_tokens":2361,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T14:29:10.676122+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[],"review_version":1}