A dressed polarizability framework for interface-coupled meta-atoms and large-scale metasurfaces

Sintesi Tecnica: Un Framework di Polarizzabilità "Vestita" per Meta-Atomi Accoppiati a Interfacce e Metasuperfici su Larga Scala

Problema
La modellazione accurata di grandi collezioni di meta-atomi interagenti rimane una sfida centrale nella fotonica, in particolare quando queste strutture coinvolgono geometrie complesse, scattering multiplo e ambienti eterogenei. Sebbene le metasuperfici ordinate e disordinate condividano la stessa base fisica di scatterer interagenti, esse vengono tipicamente modellate utilizzando framework numerici disparati. I sistemi ordinati spesso si affidano a tecniche di discretizzazione full-wave (ad es., FDTD, RCWA), che soffrono di elevati costi computazionali, mentre i sistemi disordinati utilizzano formalismi di scattering multiplo come il metodo della T-matrix o la Discrete Dipole Approximation (DDA).

I metodi di scattering multiplo esistenti affrontano limitazioni significative quando applicati alle moderne architetture di metasuperfici:

1. Limitazioni della T-matrix: Per i pilastri dielettrici ad alto rapporto d'aspetto comuni nelle metasuperfici, la sfera circoscritta richiesta è spesso significativamente più grande dello scatterer fisico. Ciò porta a rappresentazioni di campo inefficienti e a prestazioni degradate, specialmente quando i pilastri sono strettamente spaziati, causando la sovrapposizione delle sfere.
2. Limitazioni della DDA: Rappresentare accuratamente i risonatori dielettrici ad alto indice richiede un gran numero di dipoli per ogni meta-atomo, aumentando la dimensione del problema di scattering multiplo e annullando i vantaggi computazionali derivanti dall'evitare la discretizzazione volumetrica.
3. Vincoli di interfaccia: Il metodo originale della Global Polarizability Matrix (GPM), che offre una rappresentazione compatta basata su operatori, è limitato a meta-atomi circondati da un mezzo omogeneo. Non può descrivere accuratamente meta-atomi che toccano, attraversano o sono immersi all'interno di interfacce dielettriche — una configurazione onnipresente nella nanofotonica.

Metodologia: La Dressed Global Polarizability Matrix (dGPM)
Per affrontare queste limitazioni, gli autori introducono la Dressed Global Polarizability Matrix (dGPM). Questo framework estende il formalismo GPM originale agli ambienti fotonici eterogenei incorporando esplicitamente l'influenza delle interfacce vicine nello operatore di scattering.

• Base Concettuale: La dGPM rappresenta un meta-atomo complesso tramite un insieme compatto di dipoli elettrici e magnetici numerici. A differenza della DDA, dove i dipoli rappresentano correnti di polarizzazione fisiche con relazioni costitutive locali, i dipoli della dGPM sono costrutti numerici. La loro matrice di polarizzabilità associata è "appresa" da dati di scattering di riferimento, catturando sia i blocchi diagonali che quelli off-diagonali che codificano l'intera risposta elettromagnetica del meta-atomo, inclusi i accoppiamenti interni.
• Il Meccanismo di "Vestizione" (Dressing): Nella GPM originale, la superficie di apprendimento deve rimanere interamente all'interno di un mezzo omogeneo. La dGPM rimuove questa restrizione. L'operatore di scattering è appreso direttamente in presenza del mezzo stratificato (ad es., un'interfaccia). Di conseguenza, le funzioni di Green utilizzate nella fase di apprendimento includono i contributi diretti, riflessi, trasmessi ed evanescenti associati alle interfacce planari. Questo "veste" la matrice di polarizzabilità, assorbendo efficacemente le auto-interazioni mediate dal substrato (riflessioni di Fresnel e feedback del campo vicino) all'interno dell'operatore stesso.
• Implementazione:
  • Fase di Apprendimento: Simulazioni full-wave (utilizzando COMSOL) vengono eseguite per un set di sorgenti di dipolo indipendenti posizionate su una superficie di sorgente virtuale che circonda il meta-atomo. I campi diffusi vengono registrati su una superficie di apprendimento conforme. Una procedura di inversione dei minimi quadrati recupera l'operatore dGPM ($\tilde{\mathcal{A}}$) che mappa i campi incidenti nei momenti indotti.
  • Fase di Generalizzazione: La capacità predittiva dell'operatore viene validata rispetto a nuovi set di eccitazione non utilizzati durante l'addestramento.
  • Scattering Multiplo: Il framework combina naturalmente scatterer accoppiati all'interfaccia (dGPM) e nello spazio libero (GPM) all'interno di un unico formalismo di scattering multiplo utilizzando tensori di Green. La dGPM è limitata alle traslazioni parallele all'interfaccia e alle rotazioni attorno all'asse normale, mentre la standard GPM mantiene la piena libertà traslazionale e rotazionale 3D.

Risultati Chiave e Validazione
Il documento valida la dGPM attraverso diversi casi di test, dimostrando la sua capacità di gestire geometrie che vanno oltre l'applicabilità pratica degli approcci convenzionali T-matrix:

1. Meta-Atomi che Attraversano l'Interfaccia: Un micropilastro di GaAs con cappuccio in ZnS che interseca un'interfaccia dielettrica è stato modellato. La dGPM ha raggiunto un errore relativo medio di circa il 2,5% nella densità di energia e ha ricostruito accuratamente sia i campi prossimali che i pattern di radiazione in campo lontano, incluse le caratteristiche nette all'angolo critico. Questa geometria è intrattabile per i metodi T-matrix a causa della sovrapposizione delle sfere circoscritte e dell'attraversamento dell'interfaccia.
2. Dimeri Fortemente Accoppiati: È stato analizzato un sistema di due pilastri interagenti. La dGPM è rimasta accurata per separazioni fino a 0,6 µm, mentre i metodi T-matrix falliscono per separazioni inferiori a ~5,9 µm a causa della sovrapposizione delle sfere. La dGPM ha previsto con alta fedeltà i potenziamenti del campo vicino nel gap tra i pilastri e i pattern di interferenza in campo lontano.
   • Prestazioni: Rispetto alle simulazioni full-wave COMSOL, l'approccio dGPM ha offerto un'accelerazione di 700 volte nella velocità e una riduzione di 250 volte nell'uso della memoria per il caso del dimero.
3. Sistemi Ibridi: Il framework ha modellato con successo un sistema eterogeneo contenente un meta-atomo sepolto (modellato con la GPM standard) e pilastri superficiali (modellati con dGPM). L'approccio combinato ha predetto accuratamente le sezioni d'urto di assorbimento anche quando le superfici di apprendimento si sovrapponevano leggermente ($g = -0,2$ µm), dimostrando robustezza oltre i limiti di validità nominali.
4. Scalabilità: Il framework è stato applicato a ensemble fino a 1.000 meta-atomi (sia ordinati che disordinati). La scalabilità computazionale segue il previsto $O(N^2)$ per l'assemblaggio del tensore di Green e $O(N^3)$ per la risoluzione lineare densa. Gli autori osservano che, sebbene le attuali implementazioni in MATLAB non siano ottimizzate, il formalismo evita i problemi di dispersione numerica dei metodi full-wave all'aumentare delle dimensioni del sistema.

Significato e Rivendicazioni
Il documento afferma che la dGPM fornisce un framework elettromagnetico unificato a ordine ridotto che colma il divario tra la modellazione di metasuperfici ordinate e disordinate. Il suo significato primario risiede nel:

• Estensione della Modellazione Basata su Operatori: Estende con successo le rappresentazioni compatte basate su scatterer a ambienti eterogenei dove i meta-atomi interagiscono con le interfacce, uno scenario che precedentemente richiedeva la discretizzazione full-wave o metodi separati meno efficienti.
• Gestione di Geometrie Complesse: Supera i vincoli geometrici del metodo T-matrix, consentendo la simulazione di strutture ad alto rapporto d'aspetto e array densamente impacchettati dove le sfere circoscritte altrimenti si sovrapporrebbero.
• Efficienza e Accuratezza: Precomputando gli operatori specifici per l'ambiente, il metodo consente la simulazione efficiente e il design inverso di metasuperfici su larga scala con un'accuratezza comparabile alle simulazioni full-wave ma a una frazione del costo computazionale.
• Formalismo Unificato: Permette la combinazione fluida di scatterer nello spazio libero e accoppiati all'interfaccia, facilitando la modellazione di sistemi fotonici eterogenei complessi contenenti componenti sepolti, sospesi e accoppiati all'interfaccia.

Gli autori concludono che, sebbene rimangano sfide riguardanti il posizionamento ottimale dei dipoli numerici e la necessità di dataset di riferimento, il framework dGPM offre un prezioso complemento ai metodi esistenti per l'analisi e l'ottimizzazione di metasuperfici su larga scala e media fotoniche disordinate.