Pith. sign in

REVIEW 3 major objections 6 minor 12 references

Spectral tuning of hyperbolic shear polaritons in monoclinic gallium oxide via isotopic substitution

T0 review · 3 major / 6 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Substituting 18O for 16O in beta-Ga2O3 redshifts its hyperbolic shear polariton band by about 40 cm^-1 while leaving the polariton geometry essentially unchanged.

desk verdict A convincing demonstration that isotope substitution shifts hyperbolic shear polaritons in beta-Ga2O3 by ~40 cm^-1, with three independent methods agreeing; the main soft spot is a sloppy skin-depth estimate in the SI, not the central claim. read the letter →

arxiv 2507.20896 v1 pith:RZABYULM submitted 2025-07-28 physics.optics cond-mat.mtrl-sci

Giulia Carini , Mohit Pradhan , Elena Gelzinyte , Andrea Ardenghi , Saurabh Dixit , Maximilian Obst , Aditha S. Senarath , Niclas S. Mueller
show 21 more authors
Gonzalo Alvarez-Perez Katja Diaz-Granados Ryan A. Kowalski Richarda Niemann Felix G. Kaps Jakob Wetzel Raghunandan Balasubramanyam Iyer Piero Mazzolini Mathias Schubert J. Michael Klopf Johannes T. Margraf Oliver Bierwagen Martin Wolf Karsten Reuter Lukas M. Eng Susanne Kehr Joshua D. Caldwell Christian Carbogno Thomas G. Folland Markus R. Wagner Alexander Paarmann
This is my paper · ORCID
classification physics.opticscond-mat.mtrl-sci
keywords hyperbolicshearpolaritonsbeta-Ga2O3isotopicsubstitutionnear-fieldopticalmicroscopyphononreststrahlenbands-SNOMspectraltuning
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper reports that substituting the oxygen isotope 18O for natural 16O in the monoclinic crystal $\beta$-Ga2O3 shifts the hyperbolic shear polariton band downward by roughly 40 cm$^{-1}$, about 5 percent of the frequency of the high-energy transverse-optical phonons that define the band. The shift is practically relevant because hyperbolic shear polaritons, which guide infrared light along nanoscale directional rays, exist only inside narrow reststrahlen bands fixed by the crystal's phonon frequencies; isotopic substitution moves those bands to frequencies the natural crystal cannot access. The authors show the shift can be measured directly from near-field images of the polariton rays, through the orientation of the optical axis and the opening angle of the hyperbolic wave fronts, without fitting an 18O dielectric tensor. Polarized infrared reflectance and density-functional-theory calculations agree with the near-field estimate.

What carries the argument

The central object is the hyperbolic shear polariton (HShP), a surface-bound infrared mode in monoclinic crystals whose in-plane isofrequency contour is a hyperbola with two asymmetric arms; the paper images these modes in real space after launching them from a 2-µm gold disc that acts as a high-momentum antenna. The argument is carried by two angle quantities extracted from the near-field images: the optical-axis dispersion angle $\gamma(\omega)$, defined from the eigenvector of $\mathrm{Re}(\varepsilon)$ as $\gamma(\omega) = \mathrm{arctan2}(m_x, m_y)$, and the hyperbola opening angle $\alpha(\omega) = 2\arctan\sqrt{-\mathrm{Re}(\varepsilon_{nn})/\mathrm{Re}(\varepsilon_{mm})}$. Isotopic substitution enters through the Lorentz-oscillator model of the permittivity, in which the transverse-optical frequencies shift with oxygen mass while the oscillator orientation vectors $S_i$ remain essentially unchanged, so the $\gamma$ and $\alpha$ values read from the ray directions already contain the isotope shift without any need to know the 18O dielectric tensor.

What would settle it

Thin the 18O epilayer to about half a micrometer, or excite the same film with a larger antenna that launches lower momenta, and check whether the apparent 40 cm$^{-1}$ shift splits into two bands or moves toward the 16O substrate response; the semi-infinite assumption predicts no substrate contribution, while a substrate-coupling alternative predicts a mixed or thickness-dependent shift.

Watch

Extended reading notes

Core claim

On the paper's own terms, the central discovery is that 18O bGO supports the same hyperbolic shear polaritons as 16O bGO, only shifted about 40 cm$^{-1}$ lower in frequency: a near-field image at 682 cm$^{-1}$ on the 18O epilayer reproduces the pattern seen at 720 cm$^{-1}$ on 16O bGO. The frequency-dependent optical-axis angle and hyperbola half-opening angle extracted from real-space ray propagation match the curves computed from the FT-IR permittivity of 18O bGO, and both are the 16O curves displaced by the same amount. Table 1 quantifies the shift in the in-plane transverse-optical phonons: the high-frequency B$_u$ modes #1-#4 shift by 4.4-5.6% (FT-IR) and 4.9-5.5% (DFT), while the oscillator orientation vectors stay essentially fixed. Because modes #1 and #2 dictate the permittivity in the hyperbolic band, the ~5% shift of these modes yields the observed mid-band shift of about 40 cm$^{-1}$. The authors conclude that isotopic substitution tunes the spectral position of hyperbolic shear polaritons while preserving their propagation characteristics.

Load-bearing premise

The load-bearing assumption is that the 1.2-micrometer-thick 18O epilayer behaves as a semi-infinite crystal for the launched polaritons, so the near-field images carry no contribution from the underlying 16O substrate.

Editorial extensions

If this is right

  • The 40 cm-1 redshift places the 18O bGO hyperbolic band at frequencies where natural bGO is elliptical, effectively widening the spectral coverage of bGO-based nanophotonics.
  • Because the optical-axis dispersion and hyperbola opening angle are unchanged, the ray-like propagation design rules derived for natural bGO carry over to the shifted band.
  • The near-field extraction works on a 1.2 µm epitaxial film, where conventional far-field reflectance fits must model both film and substrate, making it a practical tool for thin-film isotope assessment.
  • The agreement between FT-IR and DFT relative shifts for the high-frequency modes supports using phonon masses to predict the tuning range before growth.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the semi-infinite assumption is tested by varying the epilayer thickness, the apparent shift should weaken below roughly half a micrometer; a persistent shift at smaller thickness would instead indicate the isotope effect is robust to substrate coupling.
  • The same near-field angle-extraction recipe could be applied to other monoclinic or triclinic polar crystals, turning real-space ray angles into a model-free probe of compositional or isotopic changes where the dielectric tensor is unknown.
  • Isotopically patterned heterostructures, with 16O and 18O regions on one chip, could confine the same type of mode to frequency-selective domains; the paper does not attempt such patterning.
  • A direct propagation-length measurement in the 18O film would clarify whether the larger fitted damping found for some low-frequency modes alters losses in the shifted band, since the near-field images here are not primarily loss measurements.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The paper reports spectral tuning of hyperbolic shear polaritons (HShPs) in monoclinic beta-Ga2O3 by substituting 18O for 16O. Near-field imaging of a 1.2 micrometer 18O epilayer homoepitaxially grown on a 16O substrate shows ray-like polariton patterns at 667-706 cm^-1 that closely resemble 16O patterns shifted by about 40 cm^-1, with the 16O image at 720 cm^-1 matching the 18O image at 682 cm^-1. Polarized FT-IR reflectance and DFT calculations yield TO phonon frequency redshifts of roughly 5% for the high-frequency Bu modes, while the optical-axis dispersion and hyperbola opening angle remain essentially unchanged. The authors argue that the near-field images alone provide a model-free estimate of the isotope-induced frequency shift, without requiring knowledge of the 18O dielectric tensor.

Significance. If the near-field interpretation is correct, the work demonstrates a practical route to spectral tuning of highly directional shear polaritons in a non-van-der-Waals, low-symmetry 3D crystal, with potential applications in nanophotonic devices operating in previously inaccessible frequency ranges. The study is strengthened by the combination of three independent approaches: real-space near-field imaging, azimuth-dependent FT-IR reflectance fitted with a multilayer model, and DFT calculations with isotope masses changed from first principles. The claim that the optical-axis orientation and hyperbola opening angle are nearly unchanged under isotopic substitution is supported by both the FT-IR-derived permittivity and the DFT results. The near-field analysis is genuinely model-free in the sense that the ray directions and the resulting gamma and alpha values are extracted without fitting a dielectric tensor to the near-field data, which is a useful methodological contribution for thin epitaxial layers and small samples.

major comments (3)
  1. [Section 2; SI S9] The claim that the 1.2 micrometer 18O epilayer acts as a semi-infinite crystal is not supported by a consistent quantitative estimate. The main text states that 2 micrometer Au discs launch momenta around 10 k0 with a skin depth of about 0.5 micrometers into the film, while SI S9 defines d = 1/kappa with kappa = k0 sqrt(kx0^2 - epsilon), which for k0 = 0.44 micrometer^-1, kx0 = 10, and epsilon_zz about -4 gives d about 0.22 micrometer; the text then quotes 0.46 micrometer as a factor-of-two 'skin depth' and also gives 0.84 micrometer for a '4 micrometer thick' disc at kx0 about 5. Because the assertion that the 16O substrate does not contribute to the 18O near-field images is load-bearing for the model-free shift estimate, please replace this estimate with a quantitative layered calculation, for example a transfer-matrix or full-wave simulation of the disk-launched field at the measured frequencies, and state unambiguously whether field-amplitude or intensity decay is being quoted.
  2. [Methods; SI S1 and SI S4] The signal processing is not identical for the two isotopes: the 18O images in Fig. 2 are analyzed as the raw O2A amplitude, whereas the 16O images in Fig. S4.1 are processed as S2 = -O2A * cos(O2P) to compensate for self-homodyne mixing of amplitude and phase. Since the central comparison of gamma and alpha in Fig. 4 mixes these two processing routes, and the Methods section states that amplitude and phase cannot be separated in the self-homodyne scheme, the relative extracted angles could be biased. Please apply the same processing to both isotopes, or demonstrate explicitly (for example by re-extracting the ray angles from both channels for at least one frequency per isotope) that the phase correction does not change the extracted directions.
  3. [SI S5; Table 1] The FT-IR fit for 18O is initialized from DFT TO frequencies, and Table 1 is then presented as an FT-IR/DFT confirmation of the isotope shift. Please state explicitly which parameters were free in the final fit and provide a sensitivity test, such as refitting the 18O reflectance starting from the 16O TO frequencies, to show that the fitted 18O TO positions are determined by the reflectance data rather than by the DFT seed. This would also clarify how independent the far-field confirmation is from the ab initio prediction.
minor comments (6)
  1. [Methods, DFT calculations] The text says the calculations used 'a X x X x X k-grid', which appears to be a placeholder; please supply the actual k-grid and the supercell sizes used for the finite-displacement phonon and anharmonic calculations.
  2. [Section 4] There are typos in the sentence about the hyperbolic range: 'spannig' should be 'spanning' and 'bonunded' should be 'bounded'.
  3. [Fig. 3 caption] The caption says 'blue and orange lines, respectively' while the text says 'orange and blue curves'; please make the color-to-isotope assignment consistent in both places.
  4. [Eq. (1) and surrounding text] The text refers to 'the principal eigenvector of Re(epsilon_xy)', but the optical-axis direction should be determined from the full in-plane real permittivity tensor, not from the xy component alone; please correct or clarify this notation.
  5. [SI S9] In addition to the consistency issue raised above, 'a 4 micrometer thick Au disc' should presumably read 'a 4 micrometer diameter Au disc', and the factor-of-two distinction between intensity and field decay should be defined explicitly or removed.
  6. [References] Reference 13 appears incomplete, as it lists authors and an arXiv identifier but no title; please check whether this citation is appropriate for the claim it supports.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the ~40 cm^-1 redshift is directly imaged in near-field data and independently reproduced by DFT with changed oxygen masses.

full rationale

The central claim — that 18O substitution redshifts the HShP band by ~40 cm^-1 — is directly observed in the near-field images (Fig. 1d,e: 720 cm^-1 for 16O vs 682 cm^-1 for 18O) and through model-free extraction of the ray directions and angles from s-SNOM data (SI S1-S2), without using the dielectric tensor. The DFT calculations are first-principles: the oxygen mass is changed in the harmonic phonon and oscillator-strength equations (Methods, Eq. 4), and the resulting TO shifts agree with the measured ones. The only mild coupling is that the FT-IR fit for 18O was initialized with DFT TO frequencies (SI S5), but the final fitted values differ from the DFT estimates (e.g., 703.6 vs 696.76 cm^-1 for mode #1), so the FT-IR column is not forced to equal the DFT column by construction. The analytical gamma(omega) and alpha(omega) curves derived from FT-IR permittivities are used for comparison with, not as the source of, the near-field observations. The semi-infinite-film and skin-depth assumptions (SI S9) concern experimental validity rather than circular reasoning. Overall, the derivation chain is self-contained: observation, first-principles calculation, and far-field fitting corroborate one another without any load-bearing step reducing to its own input.

Assumptions & free parameters 4 free parameters · 4 assumptions · 0 invented entities

The central 40 cm^-1 shift is not a fitted constant: it appears in raw near-field images, in the FT-IR TO frequencies, and in DFT. The Lorentz model parameters are free parameters but are used only for the analytic comparison curves, not to generate the shift. One fragile ad hoc assumption is the semi-infinite treatment of the 18O film.

free parameters (4)
  • Eight Bu TO phonon frequencies (omega_TO,i) = SI Table S1, e.g. 703.6 cm^-1 for 18O mode 1
    Fitted to polarized FT-IR reflectance using a Lorentz oscillator model; used to compute permittivity and analytic gamma/alpha curves.
  • Eight oscillator strengths Si (magnitude and orientation angle) = SI Table S1
    Fitted to reflectance maps; needed for epsilon_nn and epsilon_mm in eq. 2.
  • Eight damping constants gamma_i = SI Table S1
    Fitted to reflectance; broaden the polariton bands and affect the width of the hyperbolic region.
  • Electronic permittivity components epsilon_infinity = SI Table S2 (3.69 to 4.04)
    Fitted or taken from literature for the FT-IR model; small effect in the polariton range.
assumptions (4)
  • domain assumption The Lorentz oscillator model (eq. 3) with eight Bu modes describes the in-plane dielectric tensor of bGO.
    Used to convert fitted phonon parameters into permittivity maps; standard for polar crystals.
  • domain assumption DFT/LDA and DFPT in the harmonic approximation give phonon frequencies and oscillator strengths accurate enough for isotope-shift comparison.
    The paper compares DFT to FT-IR and reports discrepancies in oscillator strengths and mode orientations; the frequency shifts still agree.
  • ad hoc to paper Evanescent fields launched by 2-micrometer disks decay within the 1.2-micrometer film and do not couple to the 16O substrate.
    This is load-bearing for assigning the observed modes to 18O bGO; justified in SI S9 by a rough isotropic skin-depth estimate.
  • domain assumption The near-field intensity dips and extracted ray directions map directly to the asymptotes of the in-plane hyperbolic isofrequency surface, with the optical axis midway between rays.
    Underpins the model-free gamma and alpha extraction; based on prior literature on hyperbolic shear polaritons.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Spectral tuning of hyperbolic shear polaritons in monoclinic gallium oxide via isotopic substitution." pith.science (2026). https://pith.science/paper/RZABYULM

@misc{pith2026250720896,
  author       = {Pith},
  title        = {Pith review of: Spectral tuning of hyperbolic shear polaritons in monoclinic gallium oxide via isotopic substitution},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/RZABYULM}},
  note         = {Machine review of arXiv:2507.20896}
}
abstract

Hyperbolic phonon polaritons - hybridized modes arising from the ultrastrong coupling of infrared light to strongly anisotropic lattice vibrations in uniaxial or biaxial polar crystals - enable to confine light to the nanoscale with low losses and high directionality. In even lower symmetry materials, such as monoclinic $\beta$-Ga$_2$O$_3$ (bGO), hyperbolic shear polaritons (HShPs) further enhance the directionality. Yet, HShPs are intrinsically supported only within narrow frequency ranges defined by the phonon frequencies of the host material. Here, we report spectral tuning of HShPs in bGO by isotopic substitution. Employing near-field optical microscopy to image HShPs in $^{18}$O bGO films homo-epitaxially grown on a $^{16}$O bGO substrate, we demonstrate a spectral redshift of $\sim~40~$cm$^{-1}$ for the $^{18}$O bGO, compared to $^{16}$O bGO. The technique allows for direct observation and a model-free estimation of the spectral shift driven by isotopic substitution without the need for knowledge of the dielectric tensor. Complementary far-field measurements and ab initio calculations - in good agreement with the near-field data - confirm the effectiveness of this estimation. This multifaceted study demonstrates a significant isotopic substitution induced spectral tuning of HShPs into a previously inaccessible frequency range, creating new avenues for technological applications of such highly directional polaritons.

Figures

Figures reproduced from arXiv: 2507.20896 by the authors.

Figure 1
Figure 1. Isotopic substitution-induced spectral shift of the hyperbolic shear polaritons in β-Ga2O3. (a) 3D sketch of the Free Electron Laser(FEL)-coupled s-SNOM setup, showing the launching of HShPs by a Au disc on the surface of a 1.2 µm thick 18O bGO film homo-epitaxially grown on a 16O bGO substrate. The inset displays a schematic stick-and-ball representation of the conventional unit cell of bGO. The lattice constants w… view at source ↗
Figure 2
Figure 2. Near-field propagation of low-symmetry hyperbolic polaritons in isotopically sub￾stituted 18O bGO film on 16O bGO substrate. (a-h) Near-field microscopic images of HShPs in a 1.2 µm thick isotopically substituted bGO film homoepitaxially grown on a 16O bGO substrate. The excitation of the modes occurs though a Au disk with a 2 µm diameter. The panels show images recorded at different excitation frequencies. For all … view at source ↗
Figure 3
Figure 3. Experimental and theoretical permittivity data for the 16O and 18O bGO. (a-f) Real parts of εxx, εyy, and εxy for 16O and 18O bGO, blue and orange lines, respectively, as extracted from azimuth-dependent polarized reflection measurements (a-c) and ab initio calculations (d-f). The solid vertical lines indicate the positions of the theoretically derived TO phonon frequencies for the 16O bGO isotope. The numbers refer… view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: Experimental verification of the spectral tuning of the optical axis dispersion γ and the hyperbola opening angle α by isotopic substitution. (a) Illustrative sketch of the an￾gles γ and α. While γ is defined between the crystallographic a-axis and the major polarizabi…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

12 extracted references · 12 canonical work pages

  1. [1]

    https://www.jawoollam.com/ellipsometry-software/wvase

    WV ASE - ellipsometric analysis program. https://www.jawoollam.com/ellipsometry-software/wvase. Ac- cessed: 2025-05-16

    work page 2025

  2. [2]

    Tresguerres-Mata, C Lanza, E Terán-García, LF Álvarez-Tomillo, K Diaz-Granados, et al

    J Álvarez-Cuervo, M Obst, S Dixit, AI F. Tresguerres-Mata, C Lanza, E Terán-García, LF Álvarez-Tomillo, K Diaz-Granados, et al. Unidirectional ray polaritons in twisted asymmetric stacks. Nature communications, 15(1):9042, 2024

    work page 2024

  3. [3]

    Analytical approximations for the dispersion of electromagnetic modes in slabs of biaxial crystals

    Gonzalo Álvarez-Pérez, Kirill V Voronin, Valentyn S Volkov, Pablo Alonso-González, and Alexey Y Nikitin. Analytical approximations for the dispersion of electromagnetic modes in slabs of biaxial crystals. Physical Review B , 100(23):235408, 2019

    work page 2019

  4. [4]

    Modern scattering-type scanning near-field optical microscopy for advanced material research

    Xinzhong Chen, Debo Hu, Ryan Mescall, Guanjun You, DN Basov, Qing Dai, and Mengkun Liu. Modern scattering-type scanning near-field optical microscopy for advanced material research. Advanced Materials , 31(24):1804774, 2019

    work page 2019

  5. [5]

    Controlling the propagation asymmetry of hyperbolic shear polaritons in beta-gallium oxide

    Joseph Matson, Sören Wasserroth, Xiang Ni, Maximilian Obst, Katja Diaz-Granados, Giulia Carini, En- rico Maria Renzi, Emanuele Galiffi, Thomas G Folland, Lukas M Eng, et al. Controlling the propagation asymmetry of hyperbolic shear polaritons in beta-gallium oxide. Nature communications, 14(1):5240, 2023

    work page 2023

  6. [6]

    Towards smooth (010) β-ga2o3 films homoepitaxially grown by plasma assisted molecular beam epitaxy: the impact of substrate offcut and metal-to-oxygen flux ratio

    P Mazzolini and O Bierwagen. Towards smooth (010) β-ga2o3 films homoepitaxially grown by plasma assisted molecular beam epitaxy: the impact of substrate offcut and metal-to-oxygen flux ratio. Journal of Physics D: Applied Physics , 53(35):354003, jun 2020

    work page 2020

  7. [7]

    Determining the optical and polaritonic properties of isotopically pure hbn using cryogenic ftir micro- spectroscopy

    Siddharth Nandanwar, Aditya Desai, S Esfidani, Tristan McMillan, Eli Janzen, James H Edgar, and Thomas G Folland. Determining the optical and polaritonic properties of isotopically pure hbn using cryogenic ftir micro- spectroscopy. Applied Physics Letters , 126(1), 2025

    work page 2025

  8. [8]

    Terahertz twistoptics– engineering canalized phonon polaritons

    Maximilian Obst, Tobias Nörenberg, Gonzalo Álvarez-Pérez, Thales V AG de Oliveira, Javier Taboada-Gutiérrez, Flávio H Feres, Felix G Kaps, Osama Hatem, Andrei Luferau, Alexey Y Nikitin, et al. Terahertz twistoptics– engineering canalized phonon polaritons. ACS nano , 17(19):19313–19322, 2023

    work page 2023

Show all 12 references
  1. [9]

    Hyperbolic shear polaritons in low-symmetry crystals

    Nikolai C Passler, Xiang Ni, Guangwei Hu, Joseph R Matson, Giulia Carini, Martin Wolf, Mathias Schubert, Andrea Alù, Joshua D Caldwell, Thomas G Folland, et al. Hyperbolic shear polaritons in low-symmetry crystals. Nature, 602(7898):595–600, 2022. 17

    2022

  2. [10]

    Anisotropy, phonon modes, and free charge car- rier parameters in monoclinic β-gallium oxide single crystals

    Mathias Schubert, Rafal Korlacki, Sean Knight, Tino Hofmann, Stefan Schöche, Vanya Darakchieva, Erik Janzén, Bo Monemar, Daniela Gogova, Q-T Thieu, et al. Anisotropy, phonon modes, and free charge car- rier parameters in monoclinic β-gallium oxide single crystals. Physical Rev...

    2016

  3. [11]

    Phonon order and reststrahlen bands of polar vibrations in crystals with monoclinic symmetry

    Mathias Schubert, Alyssa Mock, Rafał Korlacki, and Vanya Darakchieva. Phonon order and reststrahlen bands of polar vibrations in crystals with monoclinic symmetry. Physical Review B , 99(4):041201, 2019

    2019

  4. [12]

    Si doping of β-ga2o3 by disilane via hybrid plasma-assisted molecular beam epitaxy

    Zhuoqun Wen, Kamruzzaman Khan, Xin Zhai, and Elaheh Ahmadi. Si doping of β-ga2o3 by disilane via hybrid plasma-assisted molecular beam epitaxy. Applied Physics Letters , 122(8):082101, 02 2023. 18

    2023

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.