{"id":"7e543df9-431d-4b6e-a045-7573c20a7537","arxiv_id":"2411.16923","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Bisingular surface polaritons at the interface of two uniaxial crystals exist only when the angle between the optic axes satisfies one of two cubic equations.","lead":"This paper finds exact conditions for a special surface light wave, called a bisingular polariton, to travel along the boundary between two anisotropic crystals. It shows the wave exists only for particular angles between the crystals' optic axes and gives formulas for those angles.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The reduction of the 4x4 boundary-condition determinant to the closed-form cubics (11) and (13) is asserted but not shown; every existence and angle-count result depends on this unverified algebra.","rationale":"The reader identified the determinant-to-cubic reduction as the weakest assumption, and I agree. Every derived result, including the special cases, the unified cubic P(tau,t)=0 in Eq. (32), the sufficient existence condition (34), and the two-to-six angle count, flows from Eqs. (11) and (13). Since the paper explicitly marks their derivation as omitted, this is a genuine load-bearing gap rather than a stylistic quibble. The feasible symbolic check would settle it, and the special-case match with [11,12,15] is real but not sufficient to certify the general factorization. No independent evidence, such as machine-checked algebra or reproducible code, is provided. I do not see a critical internal inconsistency or a replacement for the central construction, so the appropriate verdict remains CONDITIONAL, unchanged from the reader's assessment.","tokens_in":12774,"tokens_out":8840,"duration_ms":79873,"concrete_test":"Use a computer algebra system to construct the 4x4 matrix M from the boundary conditions in Appendix B with fields (2)-(6), substitute q, lambda and the phase-matching relations (8)-(10), and symbolically simplify det M for alpha'<0 and alpha'>0. Verify that det M factors identically to Eq. (11) and Eq. (13), respectively, up to a nonzero prefactor. As a cross-check, evaluate det M at the roots of the cubics for several random permittivity sets; if the determinant is not zero to machine precision, the reduction is wrong.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section II says the dispersion equation reduces to (11) and (13) after 'straightforward algebraic transformations', and Appendix B says only 'after some rather cumbersome calculations'. This is the load-bearing step: (11)-(14) are the basis for the cubic reduction (32), the existence conditions (33)-(34), the special-case roots in Sections III.A-III.B, and the claim of two to six allowed angles. If the determinant algebra is wrong, misses a sign, or has a hidden factor, the central claim fails. The paper supplies no intermediate determinant expression, no factorization, and no independent check of the cubics against direct root-finding on the 4x4 boundary-condition matrix. The agreement with the isotropic limit [11,12,15] provides corroboration in a special case, but it cannot certify the general reduction, since spurious factors can vanish in that limit. This is an addressable, finite algebraic check, not an internal inconsistency.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript treats surface polaritons at the planar interface between two generally different uniaxial media whose optic axes are parallel to the interface. It constructs a singular (Voigt-type) ansatz for the electromagnetic fields in both media, writes the boundary conditions as a 4x4 linear system, and claims that the determinant condition reduces to two closed-form cubic equations, Eqs. (11) and (13), for the angle φ between the optic axes. The authors then reduce both equations to a common cubic, Eq. (32), derive existence conditions, and analyze the number of allowed angles, giving explicit or approximate solutions in the special cases of an isotropic/weakly anisotropic medium and of similar media.","tokens_in":111,"tokens_out":30405,"duration_ms":411838,"significance":"If the dispersion equations and the subsequent root-counting analysis were correct, the paper would provide a valuable closed-form characterization of bisingular surface polaritons, including exact existence conditions and a prediction of up to six allowed mutual orientations. The special-case results can be checked analytically, and the comparisons with prior work on singular surface polaritons (refs. [15,16]) are a strength. However, the central algebraic reduction is not shown and, as detailed below, appears to contain a concrete internal inconsistency; the general claims are therefore not yet supported in their present form.","major_comments":[{"comment":"The reduction of the 4x4 determinant of the boundary-condition system to Eqs. (11) and (13) is the load-bearing step of the paper, but it is only asserted via 'straightforward algebraic transformations' (Section II) and 'rather cumbersome calculations' (Appendix B). No intermediate expression for the determinant or its factorization is provided, and no numerical check against direct root-finding on the 4x4 system is given. Without this derivation or an equivalent verification, the special-case results in Sections III.A-III.B and the whole of Section III.C rest on an unverified algebraic claim. The authors should supply the full reduction or a reproducible symbolic/numerical verification.","section":"Section II / Appendix B"},{"comment":"The claimed reduction of (11) and (13) to the common cubic P(τ,t)=0 appears to be algebraically incorrect. Consider the admissible parameter set ε⊥=1, ε′⊥=9, ε∥=15+√106, ε′∥=√106−5. Then t*=1/2, Δ1=2, Δ2=−1/2. Equation (11) has an admissible root (satisfying the constraint (12)) with sin²(φ/2)≈0.421; the corresponding t=cos²(φ*/2)/cos²(φ/2)≈1.151. Substituting this t into the right-hand side of (32) gives approximately 0.139, not zero. In polynomial form, Eq. (11) reduces in t to 6t³−8t²+3t−2=0, whereas Eq. (32) with these parameters gives 6t³−8t²−(3/2)t+4=0. These two cubics are not equivalent, so the root-counting and existence results derived from (32) are not supported as they stand.","section":"Section III.C, Eq. (32)"},{"comment":"Even if Eq. (32) were accepted, the step from P(τ,1)<0 to the condition 1−Δ1−Δ2<0 is not justified. From the definition of P, P(τ,1)=τ/(1+τ)[(1−Δ1)(1−Δ2²)−Δ1Δ2], and this expression is not controlled in sign by 1−Δ1−Δ2. For example, with Δ1=0 and Δ2=−2, one obtains P(τ,1)<0 while 1−Δ1−Δ2=3>0. Thus the sufficient condition (34) is not established by the argument given. Once the correct cubic reduction is obtained, the existence conditions and the subsequent 'two to six angles' analysis must be redone.","section":"Section III.C, Eqs. (33)-(34)"}],"minor_comments":[{"comment":"Equation (8) appears to have a typesetting defect: the square-root symbols over ε′⊥ and over the parenthesis are not visible in the printed formula, making the expression dimensionally inconsistent with Eq. (4).","section":"Eq. (8)"},{"comment":"The symbol Δ0 is used in the argument about simultaneous roots of P(t*,t) and P(t*^{-1},t) but is never defined.","section":"Section III.C"},{"comment":"In the expression for B′/B, the notation 'sp' should be written as 's·p' or defined explicitly; as printed it looks like a single variable.","section":"Appendix B"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses an interesting problem and the special-case results may be valuable, but the central algebraic reduction is not reliable as written. I have verified the counterexample in major comment 2 by hand; the authors should be asked to reconcile it, and to provide the full determinant reduction or a reproducible numerical verification of Eqs. (11), (13), and (32). Until that is done, the general existence and root-counting claims cannot be accepted. I recommend major revision rather than rejection because the issue is an addressable algebraic verification/correction within the scope of the manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper solves analytically a problem that had only been treated numerically before: bisingular (Dyakonov–Voigt) surface polaritons at the interface of two uniaxial media. The authors derive closed-form dispersion equations, existence conditions, and a count of possible angular orientations. That is a real step beyond the prior biaxial numerical work [17] and the isotropic/uniaxial limit [15]. The derivation starts from Maxwell's equations in an appendix, and the boundary-condition system is spelled out. The limiting behavior checks out: the isotropic-medium limit recovers known results for singular surface polaritons, and the identical-media limit gives a simple condition, cos^2(φ/2)=2/|η|, that is also obtainable from the Dyakonov wave dispersion equation. That is a good sign.\n\nThe main soft spot is the reduction of the 4x4 determinant to the two cubic equations (11) and (13). The paper does not show the algebra; it says 'straightforward algebraic transformations' and 'rather cumbersome calculations'. Since every subsequent result—the cubic reduction, the existence conditions, the two-to-six angle count—hangs on that reduction, this is a load-bearing step. The limiting checks are reassuring but not a full certificate: a spurious factor could vanish in those limits. This is not a fatal flaw, because the algebra is finite and checkable, and nothing in the paper suggests the result is wrong. But a referee should ask to see the factorization, or at least a symbolic/numeric check of (11) and (13) against the determinant for random parameters.\n\nA minor second point: the claimed upper bound of six allowed angles is not illustrated with a concrete permittivity set. The paper gives an example with extreme anisotropy where three roots of one cubic satisfy t>1, but not a complete physical realization with actual material parameters. That weakens the 'six angles' claim a bit, though it does not undermine the method.\n\nOverall, the paper is honest, self-contained up to the hidden algebra, and likely correct. The citation pattern is reasonable, with due credit given to prior work. The result is worth having even if the algebra is filled in later.\n\nWho is it for: anyone working on surface polaritons in anisotropic media, Dyakonov waves, or singular optics. It deserves a serious referee, not a desk reject. The main request should be to show the determinant reduction, not to reject the work.\n\nRegards.","headline":"New closed-form analytics for bisingular surface polaritons, with one load-bearing algebra step left in the footnotes.","tokens_in":13459,"tokens_out":2401,"would_cite":true,"duration_ms":21805,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Two cubic equations set all bisingular polariton angles","keywords":["bisingular surface polariton","surface polaritons","uniaxial media","singular surface waves","Voigt waves","dispersion equation","Dyakonov surface waves","optic axes"],"falsifier":"Take a specific pair of uniaxial permittivities, solve the full $4\\times 4$ boundary-condition system from Appendix B numerically without using the reduction, and check whether a singular solution with equal localization constants appears exactly at the angles predicted by the roots of $P(\\tau,t)=0$ with $t>1$; a mismatch at a single parameter set would disprove the reduction.","tokens_in":1797,"feed_emoji":"🔬","tokens_out":2331,"duration_ms":64305,"temperature":0.7,"pith_summary":"This paper establishes that a bisingular surface polariton—a surface electromagnetic wave whose field profile grows linearly with distance from the interface in both adjoining media—can exist at the plane between two nonabsorbing uniaxial crystals only when the optic axes lie parallel to the interface and take one of a discrete set of mutual angles. Those angles are the roots of two cubic equations, which the authors reduce to a single cubic in a rescaled variable, so all properties of the polariton (wavevector, propagation direction, localization constants, and field distribution) are fixed by the permittivities and the chosen angle. The result matters because it turns a generally numerical surface-wave problem into a closed-form existence condition, and it predicts that the number of allowed orientations is always even, ranging from two to six for arbitrary media. It also shows that weakly anisotropic or nearly identical media admit exactly two such angles, with simple approximate formulas.","feed_headline":"Two cubic equations set all bisingular polariton angles","feed_subtitle":"For two uniaxial crystals, the allowed mutual orientations are closed-form angles—two normally, up to six for extreme anisotropy.","key_machinery":"The central object is the bisingular surface polariton: a surface wave whose electric field in each crystal has the form $(f_0+f_1 x)\\exp(\\mp\\lambda x)$, arising when the ordinary and extraordinary localization constants coincide. The argument is carried by the reduction of the $4\\times 4$ boundary-condition determinant to two trigonometric equations, each cubic in $\\cos^2(\\varphi/2)$ or $\\sin^2(\\varphi/2)$. Through the substitution $t=\\cos^2(\\varphi_\\star/2)/\\cos^2(\\varphi/2)$ for Eq. (11) and $t=\\sin^2(\\varphi_\\star/2)/\\sin^2(\\varphi/2)$ for Eq. (13), both become the same cubic $P(\\tau,t)=0$ with $t>1$; the roots for $\\tau=t_\\star$ and $\\tau=t_\\star^{-1}$ enumerate all allowed mutual orientations, and the constraints fix the sign of $\\alpha'$. The parameters $\\Delta_1$ and $\\Delta_2$ encode the anisotropy differences between the media, and their signs control existence.","core_discovery":"For two nonabsorbing uniaxial media with optic axes parallel to the interface, a surface polariton that is singular in both media exists only for angles $\\varphi$ between the optic axes that satisfy Eq. (11) when the upper axis is tilted negatively or Eq. (13) when it is tilted positively. After a change of variables, both equations become the same cubic $P(\\tau,t)=0$ with $t>1$, so the number of allowed orientations equals the number of roots $t>1$ of the two cubics; this number is always even, usually two, and up to six for extreme anisotropy. When such an angle exists, every parameter of the polariton—wavevector $q$, direction relative to the axes, localization constants, and the relative amplitudes $A,B,A',B'$—is expressed in closed form in terms of $\\varphi$ and the permittivities. The authors further derive simple necessary conditions: if $(\\varepsilon_\\parallel-\\varepsilon_\\perp-2\\varepsilon'_\\perp)(\\varepsilon'_\\parallel+\\varepsilon'_\\perp)>0$ then at least one solution exists, while if both $\\varepsilon'_\\parallel>\\varepsilon_\\parallel$ and $(\\varepsilon_\\parallel-\\varepsilon_\\perp)(\\varepsilon'_\\parallel-\\varepsilon'_\\perp)<0$, no solution exists.","pith_inferences":["A natural extension not pursued in the paper is the absorbing case: because Voigt-type bulk waves rely on complex wavevectors, bisingular surface polaritons may persist or even multiply when weak absorption is introduced, with the closed-form roots serving as low-loss starting points.","The reduction to a cubic suggests that the same technique could enumerate double-exceptional surface waves at interfaces involving biaxial crystals, where the singular-axis geometry is richer and no closed-form count is currently available.","The paper's condition that the optic axis must be parallel to the interface is borrowed from its reference [4]; testing whether small tilt angles destroy the bisingular solution or merely shift it would determine how robust the prediction is for real fabricated crystals.","One could search for bisingular polaritons experimentally using angle-resolved near-field or scattering measurements near a uniaxial/uniaxial interface, looking for the single propagation direction predicted by the cubic roots."],"forward_implications":["For any pair of nonabsorbing uniaxial media with optic axes parallel to the boundary, the existence question reduces to counting roots $t>1$ of a single cubic; no numerical root-finding for the full surface-wave problem is required.","If one medium is isotropic or weakly anisotropic, the bisingular polariton exists in exactly two configurations, and the propagation direction approaches the known singular surface polariton at the isotropic/uniaxial boundary.","For identical or nearly identical media, the bisingular polariton propagates along the bisector of the angle between the optic axes when $\\varepsilon_\\parallel>3\\varepsilon_\\perp$ or $\\varepsilon_\\parallel<-\\varepsilon_\\perp$, with existence condition $|\\eta|>2$.","The number of allowed mutual orientations is always even and lies between two and six; six orientations require very strong anisotropy and are likely unobservable.","Bisingular polaritons do not exist when $\\varepsilon'_\\parallel>\\varepsilon_\\parallel$ and $(\\varepsilon_\\parallel-\\varepsilon_\\perp)(\\varepsilon'_\\parallel-\\varepsilon'_\\perp)<0$."],"supporting_citations":[{"why":"Supplies the structural fact that singular surface polaritons require the optic axis parallel to the boundary; the paper restricts to this configuration.","marker":"[4]"},{"why":"Gives the earlier singular surface polariton at the boundary of an isotropic and a uniaxial medium, which this paper extends to two uniaxial media.","marker":"[11]"},{"why":"Introduces the exceptional and doubly-exceptional surface-wave terminology and analyzes related singular surface waves.","marker":"[12]"},{"why":"Provides an earlier analytical singular surface polariton result for a single uniaxial medium that the present work generalizes.","marker":"[15]"},{"why":"Supplies the Dyakonov surface wave dispersion equation for identical uniaxial media; the bisingular solution (22) is obtained by equating localization constants in that equation.","marker":"[16]"},{"why":"First analyzes bisingular surface waves at the interface of two biaxial crystals but solves the dispersion equation numerically; this paper provides the closed-form uniaxial version.","marker":"[17]"}],"fun_headline_variants":["Cubic equations set bisingular polariton angles","Closed-form angles for bisingular surface polaritons","Bisingular polaritons exist only at specific axis angles","Up to six allowed angles for bisingular polaritons"],"cache_read_input_tokens":15744,"weakest_assumption_plain":"The whole result depends on the unshown algebraic step that turns the four boundary conditions at the interface into the two compact angle equations (11) and (13); if that reduction is wrong, every derived angle and existence condition falls with it.","fun_headline_variants_meta":{"raw":{"variants":["Cubic equations set bisingular polariton angles","Closed-form angles for bisingular surface polaritons","Bisingular polaritons exist only at specific axis angles","Up to six allowed angles for bisingular polaritons"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000219,"raw_usage":{"total_tokens":1485,"prompt_tokens":1028,"completion_tokens":457,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":644,"completion_tokens_details":{"reasoning_tokens":389}},"tokens_in":644,"tokens_out":457,"duration_ms":4575,"temperature":1.0,"reasoning_tokens":389,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:44:10.470135+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a specific pair of uniaxial permittivities, solve the full $4\\times 4$ boundary-condition system from Appendix B numerically without using the reduction, and check whether a singular solution with equal localization constants appears exactly at the angles predicted by the roots of $P(\\tau,t)=0$ with $t>1$; a mismatch at a single parameter set would disprove the reduction.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the structural fact that singular surface polaritons require the optic axis parallel to the boundary; the paper restricts to this configuration."},{"cited_title":"Weak anisotropy does not affect the existence condition, even for corner cases ( ε′ ⊥ ≈ ε⊥, or ε∥ ≈ ε⊥ + 2ε′ ⊥, or ε′ ∥ ≈ −ε′ ⊥)","cited_arxiv_id":null,"evidence_quote":"Gives the earlier singular surface polariton at the boundary of an isotropic and a uniaxial medium, which this paper extends to two uniaxial media."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the exceptional and doubly-exceptional surface-wave terminology and analyzes related singular surface waves."},{"cited_title":"Lakhtakia, T","cited_arxiv_id":null,"evidence_quote":"Provides an earlier analytical singular surface polariton result for a single uniaxial medium that the present work generalizes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Dyakonov surface wave dispersion equation for identical uniaxial media; the bisingular solution (22) is obtained by equating localization constants in that equation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"First analyzes bisingular surface waves at the interface of two biaxial crystals but solves the dispersion equation numerically; this paper provides the closed-form uniaxial version."}],"review_version":1}