REVIEW 3 major objections 3 minor 21 references
Bisingular surface polaritons at the interface of two uniaxial media
T0 review · 3 major / 3 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read Two cubic equations set all bisingular polariton angles
desk verdict New closed-form analytics for bisingular surface polaritons, with one load-bearing algebra step left in the footnotes. 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 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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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$.
Reading between the lines
- 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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (3)
- [Section II / Appendix B] 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 III.C, Eq. (32)] 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 III.C, Eqs. (33)-(34)] 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.
minor comments (3)
- [Eq. (8)] 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 III.C] The symbol Δ0 is used in the argument about simultaneous roots of P(t*,t) and P(t*^{-1},t) but is never defined.
- [Appendix B] 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.
Circularity Check
No circularity found: the derivation is self-contained, with benchmark checks against independent published results.
full rationale
The paper's central claim is that bisingular surface polaritons at the interface of two uniaxial media exist only for specific angles between the optic axes, determined by cubic equations (11) and (13), with all polariton parameters expressible in terms of permittivities and that angle. The derivation chain starts from Maxwell's equations, takes the singular solution ansatz (2)-(6) derived in Appendix A from the characteristic equation (A1), imposes standard electromagnetic boundary conditions, and reduces the resulting 4x4 determinant to the closed-form equations (11) and (13). None of these steps is circular: the target result (existence angles and polariton parameters) is not used as an input; permittivities and the geometry are the only inputs, and no parameter is fitted to any data. The reduction from the determinant to (11) and (13) is asserted rather than displayed ('after some rather cumbersome calculations'), but an omitted algebraic step is a correctness/verifiability concern, not circularity, especially since the paper provides independent checks: the isotropic-medium limit reproduces the known result of Marchevskii et al. and Golenitskii (Eq. (17) compared with Refs. [11,12,15]), and Eq. (22) for identical media is independently re-derived from the Dyakonov dispersion equation of Ref. [16]. These are external benchmarks, not self-referential inputs. The only self-citations (Ref. [15] by one of the present authors and Ref. [16] by the other) are used as comparison points for limiting cases, and they do not carry the load of the general existence analysis, which stands on Eqs. (11), (13), and (32). The restriction that the optic axes be parallel to the interface is imported from the literature (Ref. [4]) but is stated as a configuration assumption, not as a consequence derived from the target result, and it does not make the derivation circular. The numerical study of (32) is used only to illustrate root counts, not to fit anything. Consequently, no step in the derivation reduces by construction to its own inputs, and no fitted parameter is renamed as a prediction. The paper merits a circularity score of 0; any concerns about the unshown determinant algebra belong to correctness risk rather than circularity.
Assumptions & free parameters
assumptions (3)
- domain assumption Maxwell's equations in the frequency domain describe the electromagnetic field in nonabsorbing anisotropic dielectrics, with the dielectric tensor diagonal in principal axes.
- domain assumption The interface is planar and ideal, and the two media are homogeneous, nonmagnetic, and nonabsorbing, with real permittivities.
- domain assumption Singular surface polaritons require the optic axis to be parallel to the interface, as concluded in ref. [4].
Cite this review
Pith. "Pith review of Bisingular surface polaritons at the interface of two uniaxial media." pith.science (2026). https://pith.science/paper/JHWVZ22D
@misc{pith2026241116923,
author = {Pith},
title = {Pith review of: Bisingular surface polaritons at the interface of two uniaxial media},
year = {2026},
howpublished = {\url{https://pith.science/paper/JHWVZ22D}},
note = {Machine review of arXiv:2411.16923}
}
abstract
In some anisotropic bulk media (for example, biaxial weakly absorbing crystals) there are special directions along which the plane wave field distribution has a singular profile of the form $\propto (\mathbf{n} \mathbf{r}) \exp(i q \mathbf{n} \mathbf{r})$. They are also known as Voigt waves. Similar singular profiles also arise in the theory of surface electromagnetic waves in anisotropic media. In this work we have considered surface polaritons at the interface of two, generally different, uniaxial media. Optic axis of both media is parallel to the interface. One of the specific solutions, called bisingular, greatly simplifies the dispersion equation for surface polaritons. In this case the analytical solution in closed form is found and existence conditions have been determined. It is shown that bisingular surface polariton exists only for certain angles between the optic axes, which are found from two cubic equations. All parameters of the bisingular surface polariton depend only on permittivities and this angle. If one medium is weakly anisotropic, or both media are almost the same then two angles exist. In the general case of two arbitrary media there can be from two to six such angles.
Figures
Reference graph
Works this paper leans on
-
[17]
N. S. Averkiev and M. I. Dyakonov, Opt. Spectrosc. 68, 653 (1990)
work page 1990
-
[15]
A. Lakhtakia, T. G. Mackay, and C. Zhou, Eur. J. Phys. 42, 015302 (2020)
work page 2020
- [1]
-
[2]
A. P. Khapalyuk, Opt. Spectrosc. (USSR) 12, 52 (1962)
work page 1962
-
[3]
F. I. Fedorov and A. M. Goncharenko, Opt. Spectrosc. (USSR) 14, 51 (1963)
work page 1963
-
[4]
F. I. Fedorov, Teoriya Girotropii (Nauka i tekhnika, Minsk, 1976)
work page 1976
- [5]
- [6]
Show all 21 references
-
[7]
T. G. Mackay and A. Lakhtakia, in Adventures in Contemporary Electromagnetic Theory , edited by T. G. Mackay and A. Lakhtakia (Springer International Pub- lishing, Cham, 2023) pp. 61–91
2023
-
[8]
Brenier, J
A. Brenier, J. Opt. 17, 075603 (2015)
2015
-
[9]
V. I. Alshits and V. N. Lyubimov, J. Exp. Theor. Phys. 98, 870 (2004)
2004
-
[10]
Turpin, Y
A. Turpin, Y. V. Loiko, T. K. Kalkandjiev, and J. Mom- part, Laser Photonics Rev. 10, 750 (2016)
2016
-
[11]
Weak anisotropy does not affect the existence condition, even for corner cases ( ε′ ⊥ ≈ ε⊥, or ε∥ ≈ ε⊥ + 2ε′ ⊥, or ε′ ∥ ≈ −ε′ ⊥)
for a singular surface polariton. Weak anisotropy does not affect the existence condition, even for corner cases ( ε′ ⊥ ≈ ε⊥, or ε∥ ≈ ε⊥ + 2ε′ ⊥, or ε′ ∥ ≈ −ε′ ⊥). Our preliminary estimates show that for large η (or |η′|) and ε′ ⊥ ≈ ε⊥ requires a more careful analysis. This is...
-
[12]
F. N. Marchevski ˇi, V. L. Strizhevski ˇi, and S. V. Strizhevskiˇi, Fiz. Tverd. Tela (Leningrad) 26, 1501 (1984)
1984
-
[13]
T. G. Mackay, C. Zhou, and A. Lakhtakia, Proc. R. Soc. A. 475, 20190317 (2019)
2019
-
[14]
C. Zhou, T. G. Mackay, and A. Lakhtakia, Phys. Rev. A 100, 033809 (2019)
2019
-
[16]
K. Yu. Golenitskii, Phys. Rev. B 110, 035301 (2024)
2024
-
[18]
Lakhtakia and T
A. Lakhtakia and T. G. Mackay, J. Opt. Soc. Am. B 37, 2444 (2020)
2020
-
[19]
Jahani and Z
S. Jahani and Z. Jacob, Nat. Nanotechnol. 11, 23 (2016)
2016
-
[20]
Korzeb, M
K. Korzeb, M. Gajc, and D. A. Pawlak, Opt. Express 23, 25406 (2015)
2015
-
[21]
Galiffi, G
E. Galiffi, G. Carini, X. Ni, G. ´Alvarez-P´ erez, S. Yves, E. M. Renzi, R. Nolen, S. Wasserroth, M. Wolf, P. Alonso-Gonzalez, A. Paarmann, and A. Al` u, Nat. Rev. Mater. 9, 9 (2023)
2023
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.