Full-Wave Harmonic Balance Framework for Dispersive Time-Varying Photonic Structures

Sintesi Tecnica: Framework di Bilancio Armonico Full-Wave per Strutture Fotoniche Dispersive Tempo-Varianti

Definizione del Problema
Le strutture fotoniche tempo-varianti offrono capacità uniche come la conversione di frequenza, la non-reciprocità e il controllo dinamico del fronte d'onda, rompendo la simmetria di traslazione temporale. Tuttavia, modellare accuratamente questi sistemi rimane una sfida significativa quando i materiali coinvolti sono fortemente dispersivi (ad esempio, ossidi conduttori trasparenti, grafene, materiali epsilon-near-zero). In tali piattaforme, la modulazione temporale non può essere trattata come una semplice perturbazione istantanea di un parametro indipendente dalla frequenza. Al contrario, la causalità impone che la risposta del materiale dipenda dalla storia dell'eccitazione, il che significa che ogni armonica di Floquet generata sperimenta una risposta materiale distinta pur rimanendo accoppiata a tutte le altre armoniche attraverso la modulazione. I metodi computazionali esistenti affrontano compromessi: i metodi nel dominio del tempo (come FDTD) spesso richiedono simulazioni lunghe per raggiungere lo stato stazionario e faticano con aggiornamenti dispersivi causali, mentre i metodi analitici o di accoppiamento d'onda specializzati spesso mancano della flessibilità geometrica necessaria per gestire strutture 3D complesse, mezzi in omogenei o spettri di modulazione arbitrari.

Metodologia
Gli autori propongono un framework di bilancio armonico full-wave che risolve direttamente la risposta elettromagnetica nello stato stazionario nel dominio della frequenza, evitando l'avanzamento temporale. Il nucleo della metodologia consiste nel ricondurre il problema tempo-variante a un insieme accoppiato di sottoproblemi nel dominio della frequenza all'interno di uno spazio di Floquet esteso.

Elementi chiave della formulazione:

• Rappresentazione di Sorgente Secondaria: Le componenti tempo-varianti non sono rappresentate come condizioni al contorno tempo-dipendenti, ma come sorgenti secondarie indotte. Densità di polarizzazione volumetrica ($P$) sono utilizzate per i media bulk modulati, e densità di corrente superficiale ($J_s$) sono utilizzate per fogli conduttivi modulati.
• Espansione di Floquet: Tutti i campi e le sorgenti indotte sono espansi in una serie di Floquet troncata ($F(r, t) = \sum F_n(r)e^{j\omega_n t}$), dove $\omega_n = \omega_0 + n\Omega$.
• Formulazione del Sistema Accoppiato: Il problema viene risolto come un sistema lineare globale $(L - SC)X = Y_{ext}$.
  • $L$ rappresenta l'operatore di Maxwell nel dominio della frequenza per lo sfondo statico, valutato alla frequenza fisica di ciascuna armonica ($\omega_n$) per tenere conto della dispersione del materiale.
  • $S$ rappresenta gli operatori di iniezione di sorgente.
  • $C$ rappresenta gli operatori di generazione di sorgente derivati dall'espansione di Fourier temporale dei parametri del materiale. Questi operatori introducono un accoppiamento off-diagonal nello spazio armonico, rappresentando il meccanismo di conversione di frequenza.
• Implementazione: Il framework è implementato utilizzando il Metodo degli Elementi Finiti (FEM). Il vettore delle incognite standard del FEM viene elevato a un vettore nel dominio armonico, permettendo a ogni grado di libertà spaziale di trasportare molteplici coefficienti armonici. Ciò preserva la flessibilità geometrica della modellazione full-wave pur catturando l'accoppiamento inter-armonico.

Contributi Chiave

1. Framework Unificato: Il metodo fornisce una singola struttura algebrica per trattare sia i media volumetrici modulati che i fogli conduttivi modulati, supportando spettri di modulazione temporale arbitrari (monotono singolo, multi-tono o non sinusoidale).
2. Gestione della Dispersione: A differenza dei metodi che assumono una dispersione debole o utilizzano approssimazioni d'inviluppo, questo framework valuta la dispersione del materiale e gli operatori di radiazione alla specifica frequenza di ciascuna armonica generata.
3. Osservabili Risolti Armonicamente: La soluzione fornisce direttamente ampiezze, fasi e profili di campo in regime stazionario per tutte le armoniche di Floquet trattenute, consentendo il calcolo di scattering, assorbimento, flusso di potenza e distribuzioni di campo vicino specifici per ogni armonica senza necessità di post-elaborazione dei dati nel dominio del tempo.

Risultati
Il framework è stato validato e dimostrato attraverso due esempi rappresentativi:

1. Cilindro di Grafene Arrotolato a Pompa Parametrica:

   • Configurazione: Un cilindro dielettrico avvolto in un foglio di grafene con un livello di Fermi temporaneamente modulato ($\Omega = 2\omega_0$), operante in un regime di amplificazione parametrica degenerata.
   • Validazione: I risultati numerici per lo scattering e l'assorbimento risolti armonicamente hanno mostrato un eccellente accordo con una soluzione analitica di scattering di Floquet derivata per la geometria cilindrica.
   • Risultati: La simulazione ha catturato con successo la ridistribuzione dell'energia dalla portante all'armonica di conversione verso il basso ($n=-1$), confermando la capacità del framework di modellare la forte dispersione e la modulazione rapida in sistemi di plasmoni superficiali risonanti.

2. Metasuperficie Spazio-Temporale Ottica (basata su ITO):

   • Configurazione: Una metasuperficie riflettente composta da un reticolo di $TiO_2$ binario su uno strato di ITO (Indium Tin Oxide) epsilon-near-zero (ENZ) tempo-modulato, supportato da oro. La frequenza di plasma dell'ITO è stata modulata per guidare cambiamenti ultrafast nella permittività.
   • Applicazione: Il framework è stato utilizzato per progettare e ottimizzare la metasuperficie per la deflessione del fascio a conversione di frequenza. L'obiettivo era convertire un'onda incidente in bande laterali e reindirizzarle in un ordine di diffrazione spaziale specifico ($m=1$).
   • Risultati: Il design ottimizzato ha reindirizzato con successo le bande laterali generate ($\omega_{\pm 1}$) nel primo ordine di diffrazione spaziale. I risultati hanno dimostrato l'interazione tra la dispersione ENZ, l'accoppiamento Drude tempo-dipendente e la diffrazione spaziale, validando l'utilità del metodo per la progettazione di dispositivi ottici complessi e realistici.

Significatività
Il documento sostiene che questo framework stabilisce una piattaforma generale per modellare e progettare strutture fotoniche attive dispersive. Combinando la trasparenza fisica degli approcci di bilancio armonico con la flessibilità geometrica del FEM full-wave, il metodo supera i limiti degli esistenti solver per sistemi tempo-varianti. Esso permette ai ricercatori di quantificare direttamente come la modulazione ridistribuisce l'energia elettromagnetica tra i canali spostati in frequenza, fornendo osservabili essenziali (potenza, assorbimento e diffrazione risolti per armonica) per la progettazione di sistemi fotonici spazio-temporali di prossima generazione basati su grafene, ossidi conduttori trasparenti e metasuperfici risonanti. Il lavoro posiziona il framework non solo come uno strumento di analisi, ma anche come uno strumento per il design inverso di dispositivi che richiedono conversione di frequenza, amplificazione parametrica e controllo dinamico del fronte d'onda.