Dimensional reduction for optical beams with thermal nonlocal nonlinearity

Technische Samenvatting: Dimensionale Reductie voor Optische Bundels met Thermische Nietlokale Nonlineariteit

Probleemstelling
Optische nonlineariteiten voortvloend uit het thermorefractieve effect zorgen voor een lang reikend materieel respons dat wordt beheerst door warmtediffusie en absorptie. In gradiëntindex (GRIN) media verandert deze nietlokaliteit de modale interacties fundamenteel, zoals het stabiliseren van vortexbundels en ringsolitonen die in lokale kubische media anders onderhev aan azimutale instabiliteit zouden zijn. Hoewel de volledige ruimtelijke nietlineaire Schrödinger-vergelijking (NLSE) gekoppeld aan een warmtediffusievergelijking nauwkeurig is, is het computationeel veeleisend om de dynamica hiervan te modelleren. Hoewel er gereduceerde dimensionale modellen bestaan voor lokale Kerr-nonlineariteiten (zoals de niet-polynomiale Schrödinger-vergelijking, of NPSE), schieten deze tekort in het verklaren van de nietlokale aard van thermische responsen. De auteurs adresseren de noodzaak voor een hanteerbare, gereduceerde beschrijving van bundeldynamica in thermische nietlokale media die de afhankelijkheid van nietlineaire interacties van orde van de modus en de graad van nietlokaliteit vastlegt.

Methodologie
De auteurs breiden het NPSE-raamwerk uit om thermisch gemedieerde nietlokale nonlineariteiten te incorporeren. De aanpak omvat de volgende stappen:

1. Beheersende Vergelijkingen: Het systeem wordt gemodelleerd door een paraxiale ruimtelijke NLSE voor de optische veldenvelop $\Phi$ te koppelen aan een stationaire warmtevergelijking voor de temperatuur-geïnduceerde brekingsindexverandering. Onder de aanname van zwakke absorptie en specifieke randvoorwaarden, wordt de warmtevergelijking gereduceerd tot een 2D gescreende Poisson-vergelijking. De resulterende optische voortplantingsvergelijking bevat een nietlokale operator gedefinieerd door een convolutie met een gemodificeerde Bessel-functie kern ($K_0$).
2. Variational Ansatz: Om dimensionale reductie uit te voeren, maken de auteurs gebruik van een variationele benadering gebaseerd op een Lagrangiaanse formulering van het verliesvrije systeem. Zij nemen aan dat het transversale veldprofiel een enkele Laguerre-Gauss (LG) modus is met een willekeurige azimutale orde $m$ (waarbij $S=|m|$) en een radiale index van nul ($p=0$), hoewel het formalisme algemene $p$-waarden toestaat. Cruciaal is dat de ansatz een radiale fase-chirp term bevat, $b(z)r^2/2$, die de auteurs als noodzakelijk identificeren om de diffractie in nietlokale media correct te beschrijven.
3. Afleiding van het Effectieve Model: Door de Lagrangiaanse dichtheid over het transversale vlak te integreren, leiden de auteurs een effectieve eendimensionale actie af. Dit proces levert expliciete analytische expressies op voor de nietlokale bijdrage, die afhankelijk is van de mode-indices en de variationele transversale breedte $\sigma(z)$. De nietlokale term wordt geëvalueerd met behulp van de Hankel-transformatie van het intensiteitsprofiel, wat leidt tot een gesloten vorm voor de fundamentele Gaussische modus ($p=0, m=0$) en een algemene integraalvorm voor hogere orden.
4. Bewegingsvergelijkingen: De Euler-Lagrange vergelijkingen afgeleid van deze effectieve actie resulteren in een gekoppeld systeem:
   • Een evolutievergelijking voor de axiale amplitude $A(z)$, die een $\sigma$-afhankelijke nietlineaire term bevat.
   • Een algebraïsche vergelijking die de transversale breedte $\sigma(z)$ bepaalt op basis van de balans tussen lokale diffractie, het GRIN-potentiaal en de nietlokale thermische druk.
   • Een vergelijking voor de radiale fase-chirp $b(z)$.

Belangrijkste Resultaten

• Analytische Expressies: Het artikel leidt expliciete expressies af voor de nietlokale integraal $I_{pS}(\sigma)$, die de sterkte van de thermische interactie kwantificeert. Voor de fundamentele Gaussische modus wordt dit uitgedrukt in termen van de exponentiële integraalfunctie $Ei$.
• Rol van de Radiale Chirp: De studie toont aan dat de standaard NPSE-ansatz (zonder een radiale fase-term) onvoldoende is voor nietlokale media. De inclusie van de chirp-parameter $b(z)$ is essentieel om de correcte bundeldynamica te vangen, met name de ademende oscillaties van de bundelbreedte.
• Stationaire Toestanden: De auteurs analyseren het effectieve potentiaal $V_{eff}(\sigma)$ dat de bundelbreedte beheerst. Zij tonen aan dat in het defocusing regime, het vergroten van de nietlokaliteitslengte $L_{eff}$ het potentiaalminimum naar grotere breedtes verschuift, wat resulteert in bredere stationaire bundels. In het focusing regime creëert nietlokaliteit een dieper potentiaalput bij kleine breedtes.
• Numerieke Validatie: Het gereduceerde model wordt gevalideerd tegen volledige 2D split-step Fourier simulaties van de paraxiale diffractievergelijking.
  • In regimes van kort bereik nietlokaliteit ($L_{eff} = 0.5$), vertoont de variationele ansatz een uitstekende overeenkomst met exacte stationaire radiale intensiteitsprofielen voor diverse vortex-ordes.
  • In regimes van sterke nietlokaliteit ($L_{eff} = 5.0$), wijken de exacte bundelprofielen af van de aangenomen LG functionele vorm. Bijgevolg vertoont het variationele model discrepanties in de gedetailleerde vorm van het stationaire intensiteitsprofiel.
  • Ondanks deze statische profieldiscrepanties, legt het gereduceerde model de dynamische evolutie van de bundel accuraat vast, inclusclusief het verval van vermogen door absorptie en de gedempte ademende oscillaties van de transversale breedte, zelfs voor grote nietlokaliteitslengtes ($L_{eff} = 40$).

Betekenis en Claims
Het artikel claimt een gegeneraliseerd raamwerk te bieden voor thermisch gemedieerde interacties in GRIN-vezels. De primaire bijdrage is de afleiding van een gereduceerd 1D-model dat expliciet de afhankelijkheid van nietlineaire interacties van de mode-orde en transversale breedte incorporeert. De auteurs benadrukken dat, hoewel het model niet perfect het exacte transversale profiel kan reproduceren in hoog-nietlokale regimes (waar de bundelvorm afwijkt van de LG-ansatz), het toch een krachtig instrument blijft voor het begrijpen van de algehele bundeldynamica, zoals stabiliteit en voortplantingsevolutie. Het werk breidt de toepasbaarheid van variationele dimensionale reductietechnieken uit naar systemen met gescreende Poisson-type nietlokaliteiten, die relevant zijn niet alleen voor thermische media, maar ook voor vloeibare kristallen, elektrostrictieve effecten en plasma-interacties. De auteurs merken op dat het huidige model beperkt is tot continue golven en conservatieve systemen (waarbij absorptie fenomeenologisch is toegevoegd), en dat toekomstige generalisaties multimode velden, verzadigde responsen en tijdsafhankelijke thermische relaxatie voor gepulste systemen kunnen adresseren.