{"id":"1d134587-4b95-4924-a382-76d5665fbf0e","arxiv_id":"2501.18354","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"ATR spectroscopy reveals distinct spectral footprints of ghost and leaky hyperbolic polaritons in bulk crystal quartz, with signatures controlled by the crystal's anisotropy orientation.","lead":"This paper uses infrared attenuated total reflection (ATR) spectroscopy, backed by simulations, to show how ghost and leaky hyperbolic polaritons in crystal quartz appear in far-field spectra. Because ATR is a tabletop technique, this could make these exotic nanoscale light modes accessible to more laboratories and useful for direction-dependent optical devices.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"ATR dips are assigned to GHPs/LHPs by shape comparison to prior s-SNOM work only; the TMM minima are never verified at the mode level (eigenvector, decay, Poynting flux), so the central far-field-footprint claim rests on an unvalidated identification.","rationale":"The reader's weakest assumption is the one I would also flag: the simulated dips are matched to GHPs and LHPs by qualitative resemblance to near-field images, not by an independent test of the mode's physical nature. My stress-test sharpens this into a concrete, checkable gap. Section 2 and Fig. 4 present a hyperbola-shaped reflectance drop and call it a GHP because it 'mirrors' ref 8; Section 3 and Fig. 6 present a lenticular drop and call it an LHP because it resembles ref 11. The TMM described in Appendix B computes only the reflectance coefficients (B47-B54); it does not output the excited mode's Poynting vector or decay constants inside the crystal. Since the paper's novelty is precisely the far-field 'optical footprint' of these modes, the reflectance minimum itself is not enough: a dip in R_p can also come from Otto-prism coupling to ordinary SPhPs or from the continuum of bulk hyperbolic waves when the air gap is finite. The experimental data in Figs. 1(c) and 3 validate only the elliptical SPhP at kx/k0=1.66; no experiment is shown for the GHP (d=0.1 micron, epsilon_p=50) or LHP (d=10 micron, epsilon_p=2.2 or 5.5) configurations. Appendix F's calcite GHP calculation is another simulation, so it does not independently establish the identification. I therefore agree with the CONDITIONAL verdict: the physics is plausible and the transfer-matrix machinery is standard, but the central claim requires either a field-profile analysis or an experimental ATR measurement of the GHP/LHP configurations before it can be accepted as demonstrated. I would not move the verdict to ACCEPT or REJECT on the current evidence.","tokens_in":20597,"tokens_out":6369,"duration_ms":63998,"concrete_test":"At the frequency/angle of each claimed GHP/LHP dip (Figs. 4(b,d), 6(b,d)), reconstruct the 4x4 TMM field eigenvectors inside the quartz layer and compute the Poynting vector and complex k_z for the excited mode. For a GHP, the mode must show Re(S_z)=0 (or negligible) with decay away from the surface and in-plane S_parallel; for an LHP, S_z must be nonzero into the crystal (leakage). If the reflectance-minimum mode instead has bulk-like Poynting flux or an ordinary evanescent SPhP signature, the 'ghost/leaky' identification fails and the claim should be downgraded. Alternatively, compare the extracted pole loci of r_p(kx,ky) with the analytical GHP/LHP dispersion of refs 8/11.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that the simulated ATR spectra reveal the optical footprint of GHPs and LHPs—depends on the assumption that the reflectance minima in Figs. 4(b,d) and 6(b,d) actually correspond to those modes. The paper only compares the angular shape of the dips (hyperbola; lenticule) with s-SNOM images from refs 8 and 11; it never computes the mode profile, decay constant, or Poynting vector inside quartz at the dip. In an Otto configuration, a TMM reflectance minimum can arise from ordinary surface-phonon-polariton coupling, from frustrated TIR into bulk hyperbolic modes, or from improper (leaky) branch contributions. Without extracting the associated pole/eigenvector of r_p and checking the defining characteristics (for GHPs: wavefronts slanted away from the surface and Poynting vector parallel to the surface; for LHPs: energy leakage into the bulk), the identification is qualitative. This matters because no independent far-field measurement of these specific modes is reported; the only experimental validation is for the ordinary elliptical SPhP (Figs. 1(c), 3), and Appendix F's calcite check is itself a simulation. Therefore the headline result is conditional on a mode-identity verification that is not yet present.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the attenuated total reflection (ATR) response of crystal quartz in the Otto configuration, with the anisotropy axis tilted by angles φ and β relative to the surface and incidence plane. Using a 4×4 transfer-matrix method, the authors compute reflectance maps over frequency and in-plane wavevector for three polariton classes: ordinary elliptical surface phonon polaritons, ghost hyperbolic polaritons (GHPs), and leaky hyperbolic polaritons (LHPs). They report that GHPs produce a hyperbola-shaped reflectance dip in the Type II hyperbolic region and LHPs produce a lenticular dip in the Type I/ENZ region, with the tilt of the anisotropy axis controlling the shape and the cross-polarization conversion. The paper includes experimental ATR spectra for the ordinary surface-polariton case and compares them to simulations with a fitted air gap, while the GHP and LHP results are purely simulated and are identified by qualitative comparison to prior s-SNOM images.","tokens_in":20919,"tokens_out":5264,"duration_ms":53067,"significance":"If the identification of the simulated reflectance dips as GHPs and LHPs is correct, the paper offers a practical far-field spectroscopic route to probing these recently discovered near-field modes, which could be valuable for device design. The manuscript has several strengths: the 4×4 transfer-matrix formulation is given in full in Appendices A and B; the parameter dependence on air gap, prism permittivity, and anisotropy orientation is explored systematically in Appendices C–H; and the experimental validation for the ordinary elliptical surface polariton provides a baseline check of the numerical method. The claims of direction-dependent cross-polarization conversion are also potentially interesting. However, the central claim is currently supported only by simulated reflectance maps whose mode identity is asserted by shape resemblance to earlier near-field work, not by a mode-level analysis or by far-field experiments on GHPs and LHPs.","major_comments":[{"comment":"The hyperbola-shaped reflectance dips in Figs. 4(b,d) and 5 are identified as ghost hyperbolic polaritons solely because their angular shape mirrors the s-SNOM images of Ref. [8]. The transfer-matrix computation returns only the reflected intensity R_p; it does not by itself establish that the minimum corresponds to a ghost mode. In an Otto configuration, a reflectance minimum can also arise from frustrated total internal reflection into bulk hyperbolic modes, from coupling to an ordinary surface phonon polariton with modified dispersion, or from an improper/leaky branch contribution. To support the central claim, the authors should extract the complex pole of r_p (or the relevant eigenvector/eigenvalue of the 4×4 transfer matrix) at the dip and verify the defining GHP characteristics: phase wavefronts slanted away from the surface, Poynting vector parallel to the surface, and evanescent decay away from the surface. A comparison with the analytic dispersion relation of GHPs would also strengthen the identification.","section":"Section 2, Figs. 4(b,d) and 5"},{"comment":"The same identification issue applies to leaky hyperbolic polaritons. The lenticular reflectance minimum in Figs. 6(b,d) and 7 is linked to the lenticular isofrequency contours of Ref. [11], but the simulation does not demonstrate that the mode at the minimum has the characteristic LHP properties: a complex propagation constant, energy leakage into the bulk, and exponentially growing field behavior in the air gap on the improper branch. Without this mode-level check, the dip could be a bulk-wave artifact or a frustrated-TIR feature. The manuscript should also clarify how a rather thick air gap (d = 10 µm) still permits evanescent coupling to the leaky mode, since the text states that leaky waves grow in air but does not quantify the coupling condition.","section":"Section 3, Figs. 6(b,d) and 7"},{"comment":"The abstract states that 'Our findings reveal that the ATR spectra of GHPs exhibit a distinct hyperbolic behaviour' and that the authors can 'discern the effects of large asymmetry due to cross-polarisation conversion.' However, no experimental ATR data for GHPs or LHPs are presented; the only experimental comparison (Figs. 1(c) and 3) is for ordinary elliptical surface polaritons, and even there the air gap is a fitted parameter (d = 2 µm, with variation 1.5–3 µm). The GHP and LHP conclusions rest entirely on simulations. The authors should either present far-field measurements for these modes or explicitly reframe the claims as theoretical predictions, rather than empirical 'findings.' This distinction is load-bearing for the paper's stated contribution.","section":"Abstract and Section 5"},{"comment":"The GHP simulations use a prism permittivity of ε_p = 50 (corresponding to k_x/k_0 up to about 7), while the LHP simulations use ε_p = 2.2 and d = 10 µm. These parameter choices are motivated only by the need to reach the relevant wavevector range and are not discussed in terms of physical feasibility or sensitivity. In particular, a prism with ε_p = 50 is far outside the range of common infrared prism materials (e.g., Ge, Si, KRS-5), and the conclusions about direction-dependent behavior depend on this choice. The authors should justify that such parameters are achievable or show that the qualitative conclusions are robust over a range of ε_p and d.","section":"Section 2, 'In order to obtain the results above' and Fig. 4"}],"minor_comments":[{"comment":"The caption for Fig. 4(b) states 'the anisotropy orientation is unchanged (φ = 0◦)' but the text in Section 2 describes the GHP case with φ = 90◦; this is a typographical inconsistency that should be corrected.","section":"Fig. 4 caption"},{"comment":"The panel labels in the text and in the captions of Figs. D1 and G1 are inconsistent. For example, in Appendix D the text refers to tilting the anisotropy to φ = 60◦ in Fig. D1(b), while the caption for D1(b) says φ = 90◦, and the air-gap/no-air-gap assignments appear swapped between panels (b) and (c). The same issue affects Fig. G1. These mismatches make the appendix difficult to follow.","section":"Appendices D and G"},{"comment":"The text and figure caption contain an unresolved placeholder: the incident angle is given as 'θ = XX◦', which should be replaced with the actual numerical value.","section":"Fig. 7 and Section 3"},{"comment":"Equation (3) uses ε_xx for the in-plane dielectric component, but the surrounding text defines the in-plane components as ε⊥. The notation should be harmonized to avoid ambiguity, especially since later sections use ε_xx, ε_yy, and ε_zz for the rotated tensor.","section":"Eq. (3)"},{"comment":"The experimental section states that silica spacers of 1.5 µm were used but that the best-fitting theoretical air gap is 2 µm, with variation between 1.5 and 3 µm. The procedure for estimating this uncertainty and its effect on the fitted reflectance should be described more precisely, and the fitted value should be reported with an error bar.","section":"Section 5, experimental methods"},{"comment":"There are several typographical errors that should be corrected, including 'between between 450 and 480 cm-1' in the Introduction, 'S-NOM' for 's-SNOM' in Section 2, 'at with a constant' in the captions of Figs. F1 and H1, and the reference to 'Estevam et al.' in the text versus 'Estevâm da Silva' in Ref. [25].","section":"Various"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the transfer-matrix calculations are transparent and reproducible in principle, but the central claim is currently supported by simulated spectra whose mode identity is asserted by shape comparison to prior near-field work. The authors should be encouraged to add a mode-level analysis (poles of r_p, field profiles, Poynting vector) and to either present experimental ATR data for GHPs/LHPs or clearly label the results as theoretical predictions. The manuscript also contains several unresolved placeholders and caption inconsistencies that suggest it was not fully proofread; these should be fixed in revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: this is a solid transfer-matrix simulation of ATR spectra for quartz, and the new qualitative result—the predicted far-field footprints of ghost and leaky hyperbolic polaritons, including the cross-polarization asymmetry when the anisotropy is tilted—is worth taking seriously. But the central identification of those dips as GHP/LHP relies on shape-resemblance to prior s-SNOM images. The authors never verify at the mode level that the minima correspond to modes with the defining characteristics (Poynting vector parallel to surface for GHPs, leakage into bulk for LHPs). That is the load-bearing weakness, and the stress-test note is on point.\n\nWhat the paper does well: the 4x4 Berreman formulation is clearly laid out, the code appears reproducible from the equations, and the experimental validation for ordinary surface phonon polaritons (Figs. 1c, 3) gives confidence the TMM is correctly implemented. The appendices on air-gap dependence and the calcite check are genuinely useful. The asymmetry in Rps under tilted anisotropy is a new prediction, not just a rehash of earlier work. The authors are transparent about the fitted air gap and the spacer movement, which is honest.\n\nWhere it's soft: first, no experiment for the GHP/LHP footprints. That alone wouldn't be fatal if the mode identity were established theoretically, but it isn't—the paper identifies dips by eye against refs. 8 and 11. Could those dips be frustrated TIR into bulk modes or ordinary SPhPs? Possibly not, but the authors should check by computing the reflection pole or the field profile at the dip. Second, there are artifacts that suggest haste: \"θ = XX◦\" in Fig. 7, mislabeled figure references in Appendix G (D1 instead of G1), and an inconsistent φ=60° vs φ=70° for the leaky case. These are minor but should be cleaned.\n\nWho this is for: nanophotonics people using ATR or far-field methods to probe hyperbolic polaritons; it will be a useful citation as a predicted spectral fingerprint. I would not yet treat the fingerprint as established, but the paper is a legitimate step.\n\nRecommendation: send it to peer review. The method is sound, the claims are clear, and a referee can ask for the mode-identity check or a softened claim. With that addition, the paper would be substantially stronger.","headline":"Plausible, well-executed TMM study of ATR fingerprints of ghost and leaky hyperbolic polaritons, but the core identification relies on shape-matching to prior near-field work and needs mode-level verification to be fully convincing.","tokens_in":21406,"tokens_out":2886,"would_cite":true,"duration_ms":26621,"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":"This paper argues that ghost and leaky hyperbolic polaritons, modes previously observed only with near-field optical microscopy, have identifiable far-field signatures in attenuated total reflection spectra of crystal quartz, and that…","keywords":["attenuated total reflection","ghost hyperbolic polaritons","leaky hyperbolic polaritons","crystal quartz","anisotropy orientation","cross-polarisation conversion","4x4 transfer matrix method","hyperbolic dispersion"],"falsifier":"Measure the ATR reflectance of a polished quartz slab with the prism/air-gap parameters of the GHP calculation (high-index prism, $d=0.1$ µm, $\\varphi=90^\\circ$, 460 cm$^{-1}$): if the predicted hyperbola-shaped dip does not appear at the quoted $k_x/k_0$ and $\\beta$ values, or if it persists without the air gap, the ghost-polariton assignment fails. A computational falsifier is to compute the complex poles of the reflection coefficient and show that the dips coincide with modes whose Poynting vectors are parallel to the surface (GHP) or tilted into the bulk (LHP).","tokens_in":20414,"feed_emoji":"","tokens_out":11276,"duration_ms":97975,"temperature":0.7,"pith_summary":"Ghost and leaky hyperbolic polaritons, exotic surface modes recently imaged with near-field microscopy, leave distinct far-field footprints in attenuated total reflection (ATR) spectra of crystal quartz. The paper argues that a hyperbola-shaped reflectivity dip in the Type II hyperbolic band is the optical footprint of ghost hyperbolic polaritons (GHPs), and that a lenticular dip just beyond the critical angle in the Type I/epsilon-near-zero region is the footprint of leaky hyperbolic polaritons (LHPs). Tilting the crystal's anisotropy axis away from the surface weakens the GHP dip and, for LHPs, produces a large asymmetric cross-polarisation conversion between p- and s-polarised reflected light. The authors conclude that ATR spectroscopy can serve as a far-field probe for these modes and that crystal orientation is a control parameter for direction-dependent infrared optics. The experimental evidence covers ordinary surface phonon polaritons, while the GHP and LHP spectra are simulated with the transfer-matrix method.","feed_headline":"Far-field dips reveal ghost and leaky polaritons","feed_subtitle":"A tabletop infrared measurement could replace s-SNOM for spotting these modes.","key_machinery":"The engine of the calculation is a 4×4 transfer-matrix method for stratified anisotropic media, applied to the prism/air-gap/crystal stack that couples evanescent waves to surface modes. The dielectric tensor of quartz, uniaxial with components $\\varepsilon_\\parallel$ and $\\varepsilon_\\perp$, is rotated by two angles: $\\varphi$ tilts the anisotropy axis away from the surface normal, and $\\beta$ rotates the incidence plane around that axis; the rotated tensor acquires off-diagonal elements. The method returns the full set of reflection coefficients $r_{pp}$, $r_{ps}$, $r_{sp}$, $r_{ss}$, and the paper tracks the total p-polarised reflectance $R_p = R_{pp}+R_{ps}$ as a function of azimuthal angle $\\beta$ and reduced in-plane momentum $k_x/k_0$. The spectral 'footprints' are the shapes of the dips in this reflectance: a hyperbola for GHPs and a lens for LHPs, with their position and intensity controlled by $\\varphi$, $\\beta$, the air gap distance $d$, and the prism index.","core_discovery":"The central claim is that the ghost and leaky hyperbolic polaritons discovered in bulk anisotropic crystals through near-field imaging can be recognised in the far field by their ATR spectra. For GHPs in the Type II hyperbolic region of quartz at 460 cm$^{-1}$, the calculated reflectance, plotted against azimuthal angle and in-plane momentum, shows a hyperbola-shaped dip outside the bulk propagation bands; this is described as mirroring the angular dispersion of the original near-field ghost-polariton observation. For LHPs at 545 cm$^{-1}$, a lenticular dip appears just beyond the critical angle when a 10 µm air gap is introduced, matching the lenticular isofrequency contours seen for leaky polaritons in calcite. When the anisotropy is tilted away from the surface ($\\varphi = 60^\\circ$ for GHPs, $\\varphi = 70^\\circ$ for LHPs), the dips shift and weaken, and cross-polarisation conversion becomes strongly asymmetric, with $R_{ps}$ values reaching roughly 0.7 on one side and 0.1 on the other for LHPs. The paper presents this orientation-controlled response as the basis for direction-dependent optical devices.","pith_inferences":["The mode assignment could be tested independently by comparing dip positions with a direct complex-wavevector dispersion calculation of GHPs and LHPs in the same quartz geometry; the paper stops short of such an eigenvalue analysis.","If the assignment holds, the opening angle of the LHP lenticular dip as the frequency crosses the epsilon-near-zero region could become a contact-free spectroscopic estimator of $\\varepsilon_\\parallel(\\omega)$ in the reststrahlen band.","Because the GHP and LHP spectra are purely simulated, a high-index-prism ATR experiment with controlled air gap and crystal cut is the natural next step; it would also test the predicted $R_{ps}$ asymmetry of about 0.7 versus 0.1.","The reciprocal asymmetric cross-polarisation response suggests that a single tilted anisotropic crystal might act as a passive polarization-dependent beam splitter in the infrared, although absorption losses will set the practical efficiency."],"forward_implications":["ATR spectroscopy becomes a far-field, tabletop screening tool for ghost and leaky hyperbolic polaritons in bulk crystals, complementing near-field s-SNOM imaging.","Rotating the crystal orientation ($\\varphi$ and $\\beta$) gives a practical way to shift, deepen, or suppress the polariton dips, enabling direction-dependent infrared filters and couplers.","The large asymmetric cross-polarisation conversion of tilted LHPs could be used for direction-selective or one-way infrared signal routing.","GHPs respond at high in-plane momenta, so they could serve as wide-field-of-view receivers, while LHPs operate in a narrow momentum window near the light line, suited to directional emitters.","The same transfer-matrix approach can be applied to shear hyperbolic polaritons and twisted-layer polariton systems to search for analogous far-field fingerprints."],"supporting_citations":[{"why":"Supplies the near-field s-SNOM discovery of ghost hyperbolic polaritons in calcite, the angular dispersion that the GHP ATR dip is said to mirror.","marker":"[8]"},{"why":"Supplies the near-field observation of directional leaky polaritons and the lenticular isofrequency contours that the LHP ATR dip is said to reproduce.","marker":"[11]"},{"why":"Provides the quartz dielectric function parameters used in every transfer-matrix simulation and in the experimental fits.","marker":"[25]"},{"why":"Documents the asymmetric cross-polarisation conversion in tilted hyperbolic media that the paper invokes to explain the $R_{ps}$ asymmetry.","marker":"[34]"},{"why":"Provides the anisotropic-multilayer resonance treatment used to compute evanescent polariton coupling in the ATR configuration.","marker":"[37]"},{"why":"Establishes the attenuated-total-reflection method on quartz for surface polaritons, the measurement geometry extended here.","marker":"[56]"},{"why":"Introduces the 4x4 transfer-matrix formalism that underlies all calculated ATR spectra in the paper.","marker":"[68]"},{"why":"Shows how rotating the anisotropy axis rotates the hyperbolic dispersion and creates off-diagonal tensor components, a central control mechanism.","marker":"[5]"},{"why":"Provides the temperature-dependent phonon-mode model on which the quartz dielectric tensor is built.","marker":"[73]"}],"fun_headline_variants":["Far-field spectra reveal ghost and leaky polaritons","Anisotropy tilts polariton optical footprints","ATR spots ghost and leaky hyperbolic polaritons","Orientation controls polariton spectral dips","Ghost and leaky polaritons leave far-field traces"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the simulated ATR dips really are ghost and leaky hyperbolic polaritons; the paper identifies them by qualitative resemblance to near-field s-SNOM images from earlier experiments, not by a far-field measurement or an independent mode analysis of its own.","fun_headline_variants_meta":{"raw":{"variants":["Far-field spectra reveal ghost and leaky polaritons","Anisotropy tilts polariton optical footprints","ATR spots ghost and leaky hyperbolic polaritons","Orientation controls polariton spectral dips","Ghost and leaky polaritons leave far-field traces"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000156,"raw_usage":{"total_tokens":1222,"prompt_tokens":950,"completion_tokens":272,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":566,"completion_tokens_details":{"reasoning_tokens":197}},"tokens_in":566,"tokens_out":272,"duration_ms":3177,"temperature":1.0,"reasoning_tokens":197,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T23:49:02.755193+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the ATR reflectance of a polished quartz slab with the prism/air-gap parameters of the GHP calculation (high-index prism, $d=0.1$ µm, $\\varphi=90^\\circ$, 460 cm$^{-1}$): if the predicted hyperbola-shaped dip does not appear at the quoted $k_x/k_0$ and $\\beta$ values, or if it persists without the air gap, the ghost-polariton assignment fails. A computational falsifier is to compute the complex poles of the reflection coefficient and show that the dips coincide with modes whose Poynting vectors are parallel to the surface (GHP) or tilted into the bulk (LHP).","supporting_citations":[{"cited_title":"Ghost hyperbolic surface polaritons in bulk anisotropic crystals,","cited_arxiv_id":null,"evidence_quote":"Supplies the near-field s-SNOM discovery of ghost hyperbolic polaritons in calcite, the angular dispersion that the GHP ATR dip is said to mirror."},{"cited_title":"Observation of directional leaky polaritons at anisotropic crystal interfaces,","cited_arxiv_id":null,"evidence_quote":"Supplies the near-field observation of directional leaky polaritons and the lenticular isofrequency contours that the LHP ATR dip is said to reproduce."},{"cited_title":"Far-infrared slab lensing and subwavelength imaging in crystal quartz,","cited_arxiv_id":null,"evidence_quote":"Provides the quartz dielectric function parameters used in every transfer-matrix simulation and in the experimental fits."},{"cited_title":"Asymmetric Reflection Induced in Reciprocal Hyperbolic Materials,","cited_arxiv_id":null,"evidence_quote":"Documents the asymmetric cross-polarisation conversion in tilted hyperbolic media that the paper invokes to explain the $R_{ps}$ asymmetry."},{"cited_title":"Layer-resolved resonance intensity of evanescent polariton modes in anisotropic multilayers,","cited_arxiv_id":null,"evidence_quote":"Provides the anisotropic-multilayer resonance treatment used to compute evanescent polariton coupling in the ATR configuration."},{"cited_title":"Optics in Stratified and Anisotropic Media: 4×4-Matrix Formulation,","cited_arxiv_id":null,"evidence_quote":"Introduces the 4x4 transfer-matrix formalism that underlies all calculated ATR spectra in the paper."},{"cited_title":"Oriented Asymmetric Wave Propagation and Refraction Bending in Hyperbolic Media,","cited_arxiv_id":null,"evidence_quote":"Shows how rotating the anisotropy axis rotates the hyperbolic dispersion and creates off-diagonal tensor components, a central control mechanism."},{"cited_title":"Temperature dependence of transverse and longitudinal optic modes in the $\\ensuremath{\\alpha}$ and $\\ensuremath{\\beta}$ phases of quartz,","cited_arxiv_id":null,"evidence_quote":"Provides the temperature-dependent phonon-mode model on which the quartz dielectric tensor is built."}],"review_version":1}