A robust and efficient method to calculate electromagnetic modes on a cylindrical step-index nanofibre

Titel: Een robuuste en efficiënte methode voor het berekenen van elektromagnetische modi op een cilindrische stap-index nanofiber

Auteurs: Sebastian Golat en Francisco J. Rodríguez-Fortuño (King's College London)
Datum: 17 februari 2026

---

1. Het Probleem

De nauwkeurige berekening van geleide elektromagnetische modi in optische nanofibers is essentieel voor toepassingen in nanofotonica, variërend van kwantums interfaces tot vectoriële lichtsensoren. De huidige standaardmethode, zoals beschreven in leerboeken, vereist het oplossen van een $4 \times 4$ matrix-eigenwaardeprobleem om de modale velden te vinden.

Hoewel deze methode wijdverspreid is, heeft deze een fundamenteel en significant nadeel:

• Numerieke instabiliteit: De definitieve bepaling van de veldamplitudes vereist het vinden van de numerieke nulruimte (null space) van een theoretisch singuliere matrix.
• Fouten in longitudinale componenten: Dit is een slecht gesteld (ill-conditioned) probleem. Omdat de longitudinale veldcomponenten ($E_z$ en $H_z$) vaak veel kleiner zijn dan de transversale componenten, leiden kleine numerieke fouten in de nulruimteberekening tot grote relatieve fouten in deze cruciale longitudinale componenten.
• Impact: Deze fouten zijn kritiek voor toepassingen die afhankelijk zijn van chirale effecten en vectoriële licht-materie interacties, waar de polarisatie en de longitudinale veldsterkte van doorslaggevend belang zijn.

2. Methodologie

De auteurs introduceren een fundamenteel robuustere en efficiëntere semi-analytische methode die de standaard $4 \times 4$-benadering vervangt. De kern van de methode bestaat uit:

• Symmetrie-uitbuiting: De methode start vanuit de fundamentele symmetrieën van de cilindrische golfgeleider (translatie langs de as en rotatie rond de as).
• Modale expansie: De velden worden uitgedrukt als een expansie in eigenvectoren die de symmetrie-operatoren diagonaliseren, gebruikmakend van Bessel- en Hankelfuncties.
• Slimme normalisatie: Door een specifieke normalisatie van de veldamplitudes toe te passen (waarbij de longitudinale componenten aan de grens direct worden gekoppeld aan twee amplitudes $A$ en $B$), kan het oorspronkelijke stelsel van vier vergelijkingen analytisch worden gereduceerd.
• Reductie naar $2 \times 2$: De auteurs tonen aan dat het probleem equivalent is aan een veel eenvoudiger $2 \times 2$-systeem. De eerste twee rijen van de oorspronkelijke matrix geven direct $A=C$ en $B=D$, waardoor de onafhankelijke variabelen worden gereduceerd.

3. Belangrijkste Bijdragen

Deze herformulering biedt twee doorslaggevende voordelen:

1. Vereenvoudigde dispersierelatie: De dispersierelatie wordt verkregen door de determinant van de $2 \times 2$-matrix gelijk te stellen aan nul. Dit resulteert in een eenvoudige, goed-gedragde transcendentale vergelijking die numeriek zeer robuust op te lossen is (geen singulariteiten bij nul van de Bessel-functies).
2. Analytische bepaling van amplitudes: Zodra de propagatieconstante ($k_z$) en frequentie ($\omega$) bekend zijn, worden de veldamplitudes ($A$ en $B$) analytisch bepaald uit de overige twee vergelijkingen.
   • Dit elimineert volledig de noodzaak voor de numerieke berekening van de nulruimte.
   • Hierdoor blijft de nauwkeurigheid van de kleine longitudinale componenten behouden, wat essentieel is voor de volledige vectoriële veldstructuur.

De methode is geïmplementeerd in het open-source Python-pakket Anafibre.

4. Resultaten en Implementatie

• Dispersieberekening: De auteurs definiëren een genormaliseerde frequentie $V$ en een genormaliseerde propagatieconstante $b$. De wortels van de functie $F_\ell(b, V)$ leveren de effectieve brekingsindex $n_e$ voor elke modus.
• Modusclassificatie: De methode behandelt zowel niet-gedegenereerde modi ($\ell=0$, TE en TM) als gedegenereerde hybride modi ($\ell \neq 0$, HE en EH). Voor hybride modi wordt een nieuwe, meer algemene definitie voor de classificatie (HE vs. EH) voorgesteld, gebaseerd op de relatieve fase tussen $E_z$ en $H_z$ in plaats van alleen op de grootteverhouding.
• Veldreconstructie: De volledige vectoriële veldverdeling (elektrisch en magnetisch) kan exact worden gereconstrueerd binnen en buiten de fiber zonder numerieke integratie voor de normalisatie, dankzij analytische uitdrukkingen voor de vermogensintegralen.
• Generaliteit: De methode is geldig voor willekeurige materiaaleigenschappen, inclusief dispersieve en verliezende media (complex $\epsilon$ en $\mu$).

5. Significance (Betekenis)

Deze studie biedt een krachtig en betrouwbaar gereedschap voor het ontwerp en de analyse van nanofiber-gebaseerde apparaten.

• Kwaliteitsverbetering: Het elimineert de numerieke fouten die traditioneel de nauwkeurigheid van longitudinale veldcomponenten beperken.
• Toepassingsgebied: De methode is van cruciaal belang voor toepassingen waar precieze kennis van veldpolarisatie en longitudinale componenten vereist is, zoals:
  • Chirale kwantumoptica.
  • Vectoriële lichtsensoren.
  • Optische valsystemen voor atomen en deeltjes.
  • Interacties tussen licht en materie in de nabije veldregime.

De herformulering combineert conceptuele helderheid met praktische numerieke voordelen, waardoor het een nieuwe standaard wordt voor de modellering van cilindrische stap-index golfgeleiders.