Spectrally smooth broadband response via autocorrelation-constrained inverse design

Technische Samenvatting: Spectraal Glad Brede Respons via Autocorrelatie-beperkte Inverse Design

Probleemstelling
Time-domain inverse design in de fotonica wordt algemeen erkend als efficiënt voor het maximaliseren van de optische prestaties van optische componenten over brede frequentiebereiken. Standaard benaderingen maken gebruik van tijd-geïntegreerde veldgroottes (bijv. energie, flux of gedissipeerd vermogen) afgeleid van de stelling van Poynting. Door de stelling van Parseval komt overeen dat deze tijd-geïntegreerde doelstellingen corresponderen met frequentie-geïntegreerde responsen. De auteurs identificeren echter een cruciale beperking: hoewel deze geïntegreerde doelstellingen computationeel efficiënt zijn, meten ze slechts de totale respons over een bandbreedte. Ze slagen er niet in om ongewenste in-band spectrale rimpelingen, smalle spectrale kenmerken (zoals dalen) of zijlobben te detecteren of te bestraffen. Bijgevolg kan een optimizer convergeren naar een lokaal maximum van de totale energie, terwijl de resulterende spectrale respons "ruig" is met scherpe defecten die voor de doelstelling effectief onzichtbaar zijn. Bestaande oplossingen, zoals adaptieve spectrale weging, vereisen de expliciete identificatie van onderpresterende regio's en de dynamische modificatie van de excitatiepuls tijdens de optimalisatie, wat de complexiteit verhoogt.

Methodologie
De auteurs stellen een time-domain metriek voor die ontworpen is om spectrale gladheid te kwantificeren zonder dat aanpassing van de puls nodig is. De kern van de methodologie rust op de stelling van Wiener–Khinchin, die vaststelt dat de autocorrelatiefunctie van een signaal en de energiedichtheid van het spectrum een Fourier-transformatiepaar vormen.

• De Metriek ($T_{AC}$): De auteurs definiëren een metriek gebaseerd op de gewogen lang-lag autocorrelatie-energie. Een brede, gladde spectrum correspondeert met een snel afnemende autocorrelatiefunctie. Daarentegen manifesteren smalle spectrale kenmerken (zoals dalen of rimpelingen) zich als langzaam afnemende, oscillerende staarten in de autocorrelatie.
• Formulering: De metriek $T_{AC}[s; s_0]$ wordt gedefinieerd als de ratio van de gewogen lang-lag autocorrelatie-energie van het respons-signaal $s(t)$ ten opzichte van die van een breedbandig referentie-signaal $s_0(t)$. De wegefunctie neemt kwadratisch toe met de lag ($\tau^2$) om de "ringing"-staarten in de late tijd te benadrukken.
  • $T_{AC} = 1$ impliceert een respons die vergelijkbaar is met het gladde referentiemodel.
  • $T_{AC} > 1$ duidt op de aanwezigheid van scherpe spectrale kenmerken.
• Optimalisatiekader: Deze metriek wordt geïntegreerd in een density-based topology optimization (TopOpt) kader met behulp van de Finite-Difference Time-Domain (FDTD) methode. De metriek wordt op twee manieren toegepast:
  1. Als een Constraint: Het maximaliseren van reflectie-efficiëntie ($\eta_r$) onder de voorwaarde $T_{AC} \leq 1$.
  2. Als een Strafterm (Penalization): Het maximaliseren van een gewogen doelstelling $J = (1-\alpha)\eta_r - \alpha T_{AC}$, waarbij $\alpha$ de balans bepaalt tussen efficiëntie en spectrale gladheid.

Belangrijkste Bijdragen

1. Nieuwe Time-Domain Metriek: Introductie van een metriek die de detectie van smalle spectrale kenmerken vertaalt naar een time-domain grootheid (lang-lag autocorrelatie-energie), waardoor deze toegankelijk wordt voor standaard gradiënt-gebaseerde optimizers in het tijdsdomein.
2. Demonstratie op Bragg-roosters: Toepassing van de metriek op de inverse design van eendimensionale diëlektrische Bragg-spiegels. De studie toont aan dat het maximaliseren van de totale energie alleen onvoldoende is om spectrale defecten te voorkomen, terwijl de autocorrelatie-constraint deze defecten effectief onderdrukt.
3. Fysisch Inzicht: De methodologie laat zien dat bekende fysische ontwerpstrategieën (geapodiseerde en gechirpte roosters) puur door optimalisatie kunnen worden gereproduceerd, wat de capaciteit van de metriek valideert om de solver naar fysisch betekenisvolle, gladde breedbandige oplossingen te leiden.

Resultaten
De auteurs evalueerden de aanpak met behulp van 1D diëlektrische roosters met een zwak brekingsindexcontrast ($\Delta n \approx 0.22$) over een bandbreedte van 105 THz gecentreerd rond 300 THz.

• Onbeperkt versus Beperkt (Unconstrained vs. Constrained): In onbeperkte optimalisaties vertoonden de resulterende reflectiespectra vaak prominente smalle spectrale kenmerken (vergelijkbaar met $\pi$-verschoven roosterresonanties) met hoge $T_{AC}$-waarden (bijv. 8.6). Het afdwingen van de $T_{AC} \leq 1$ constraint verwijderde deze kenmerken consistent, wat resulteerde in gladde, flat-top reflectiespectra.
• Statistische Robuustheid: Een statistische evaluatie over 100 willekeurige initiële dichtheden toonde aan dat de constraint de mediane reflectie-efficiëntie ($\eta_r$) niet significant verminderde. De mediane $\eta_r$ bleef vrijwel identiek (0.8096 onbeperkt versus 0.8093 beperkt), wat aangeeft dat de constraint fungeert als een regularisator in plaats van een prestatiebeperker.
• Pareto-afwegingen: Door de weegparameter $\alpha$ in de gestrafte doelstelling te variëren, brachten de auteurs de afweging tussen efficiëntie en gladheid in kaart. Een zwakke straf ($\alpha = 10^{-2}$) verbeterde zelfs het lokale optimum voor efficiëntie terwijl een $T_{AC}$ nabij eenheid werd bereikt. Sterkere straffen leidden tot geapodiseerde profielen, waarbij een deel van de piekreflectie werd opgeofferd voor extreme gladheid.
• Ontdekking van Gechirpte Structuren: In een langer ontwerpgebied ($L = 17.2 \mu m$) convergeerde de optimalisatie natuurlijk naar een gechirpt roosterstructuur. Een parametrische analyse van het resulterende binaire ontwerp onthulde een lineaire variatie in de lokale roosterperiode, wat overeenkomt met de Bragg-vectoren van de beoogde spectrale grenzen. Een gereconstrueerd gechirpt rooster gebaseerd op dit patroon bereikte een vergelijkbare prestatie als het ruwe TopOpt-ontwerp, wat de fysische validiteit van de oplossing bevestigt.

Betekenis en Claims
Het artikel claimt dat de voorgestelde metriek een centrale beperking van geïntegreerde time-domain doelstellingen aanpakt door directe gevoeligheid te bieden voor gelokaliseerde spectrale defecten. De auteurs stellen dat:

• De metriek dient als een efficiënte time-domain maatstaf voor het beoordelen en onderdrukken van breedbandige spectrale fluctuaties zonder de noodzaak om de excitatiepuls aan te passen.
• Het kan worden gebruikt als een gewogen doelstelling, een constraint of een tool voor post-evaluatie.
• In tegenstelling tot methoden die een specifieke temporele enveloppe voorschrijven, legt het straffen van $T_{AC}$ een zwakkere beperking op, waardoor de optimizer de vrijheid behoudt om de gedetailleerde temporele vorm te bepalen terwijl de spectrale uniformiteit gewaarborgd blijft.
• De aanpak toepasbaar is buiten eenvoudige roosters voor andere uitdagende inverse-design problemen, zoals breedbandige near-field koppeling, achromatische focussering en breedbandige absorptie.

Dit werk wordt ondersteund door VILLUM FONDEN, en de onderliggende code is publiekelijk beschikbaar gemaakt.