High-Accuracy Semi-Analytical Method for Solving the Problem of Electromagnetic Wave Scattering by Arbitrary Ensembles of Parallel Circular Cylinders

Sintesi Tecnica: Metodo Semi-Analitico ad Alta Precisione per lo Scattering di Onde Elettromagnetiche da Parte di Insiemi Arbitrari di Cilindri Circolari Paralleli

Definizione del Problema
Il documento affronta il problema bidimensionale dello scattering di onde elettromagnetiche da parte di un cluster di un numero arbitrario ($P$) di cilindri circolari destri, paralleli, infinitamente lunghi, omogenei e non sovrapposti. I cilindri possono possedere raggi arbitrari e permittività complesse, e i loro assi possono occupare posizioni arbitrarie nel piano trasversale. Sebbene esistano soluzioni esatte per singoli cilindri tramite la teoria di Mie, i cluster di due o più elementi generalmente mancano di soluzioni analitiche semplici. Gli esistenti metodi numerici, come gli schemi a differenze finite o a elementi finiti, spesso faticano a mantenere l'accuratezza nelle regioni di near-field caratterizzate da variazioni spaziali brusche (ad esempio, gap stretti tra i cilindri) e richiedono una discretizzazione estremamente fine, portando a elevati costi computazionali e problemi di convergenza.

Metodologia
Gli autori propongono un approccio semi-analitico basato sull'espansione del campo elettromagnetico in armoniche cilindriche attorno all'asse di ciascun cilindro. Il nucleo del metodo prevede:

1. Formulazione Analitica:

   • Il problema viene ridotto a un'equazione d'onda scalare per la componente del campo parallela agli assi dei cilindri ($E_z$ o $H_z$).
   • Il campo totale all'esterno di ogni cilindro è espresso come una sovrapposizione del campo incidente e dei campi diffusi dagli altri cilindri, espansi in funzioni di Hankel. Il campo all'interno di ogni cilindro è espanso in funzioni di Bessel.
   • Il Teorema di Addizione di Graf è impiegato per traslare i campi diffusi dal sistema di coordinate locale di un cilindro al sistema locale di un altro. Ciò consente di accoppiare tutti i cilindri in un unico sistema autosufficiente.
   • Le condizioni al contorno (continuità dei campi tangenziali) vengono applicate per derivare un sistema infinito di equazioni lineari per i coefficienti di espansione ignoti (coefficienti multipolari).

2. Strategia di Soluzione Numerica:

   • Troncamento: Le serie infinite sono troncate a un ordine finito $N$, determinato dai criteri di convergenza e dal parametro di dimensione ($k_0 a$).
   • Sistema Matriciale: Il sistema troncato forma un'equazione di matrice a blocchi di dimensione $P \times M$ (dove $M = 2N+1$), che coinvolge matrici diagonali di coefficienti di scattering di singolo cilindro e matrici di traslazione off-diagonal derivate dal teorema di Graf.
   • Solver Adattivo: Per gestire la severa ill-condizione che sorge in configurazioni risonanti o densamente impacchettate (gap sub-lunghezza d'onda), gli autori implementano un solver adattivo a quattro livelli:
     • Livello 1: Decomposizione LU standard in doppia precisione (float64).
     • Livello 2: Scalatura della matrice (equilibratura) e raffinamento iterativo.
     • Livello 3: Aritmetica ad precisione arbitraria utilizzando la libreria mpmath (ad esempio, 120 cifre decimali) per sopprimere gli errori di arrotondamento durante la fattorizzazione.
     • Livello 4: Eliminazione esatta sul campo dei numeri razionali gaussiani ($\mathbb{Q}(i)$) utilizzando sympy.DomainMatrix. Questo elimina interamente gli errori di arrotondamento durante l'algebra lineare trattando gli elementi della matrice come frazioni esatte, fungendo da cross-validazione per il Livello 3.
   • Riduzione della Simmetria: Per geometrie simmetriche, il sistema viene ridotto sfruttando la simmetria speculare, diminuendo il numero di variabili indipendenti e migliorando il numero di condizionamento.

3. Verifica:

   • Il metodo è rigorosamente verificato attraverso molteplici controlli indipendenti: convergenza rispetto all'ordine di troncamento $N$, analisi dei residui delle condizioni al contorno, confronto degli indicatrici di scattering calcolate tramite diversi metodi (ampiezza di scattering vs integrazione del vettore di Poynting), validazione del teorema ottico e controlli di bilancio energetico.

Caso di Studio: Trimero di Alluminio
Il metodo è dimostrato su un trimero simmetrico di tre identici nanocilindri di alluminio ($a=10$ nm) disposti in un triangolo equilatero con un gap di $g=5$ nm ($g/a=0.5$), illuminato da un'onda piana normalmente incidente, polarizzata linearmente ($H_z$) a $\lambda=116$ nm (regime plasmonico).

• Fisica del Near-Field: Il campo del vettore di Poynting rivela una complessa ridistribuzione dell'energia. Una separatrice corre lungo l'asse di simmetria, impedendo lo scambio diretto di energia tra i due cilindri alla base, mentre si verifica un forte accoppiamento tra i cilindri alla base e il cilindro apicale. I cilindri alla base schermano il cilindro apicale, riducendo l'intensità del campo in quel punto, sebbene i flussi di energia provenienti dalla base compensino parzialmente questo effetto.
• Caratteristiche del Far-Field: Nonostante la complessa topologia del near-field, l'indicatrice di scattering nel far-field è dominata dal modo dipolo, somigliando strettamente al pattern di un singolo cilindro e di un dipolo elettrico.
• Sezioni Trasversali: Le sezioni trasversali di scattering, assorbimento ed estinzione per il trimero sono calcolate. Le efficienze del trimero sono inferiori a quelle di un singolo cilindro a causa degli effetti di schermatura, sebbene l'assorbimento sia distribuito in modo non uniforme tra i cilindri.
• Stabilità Numerica: Lo studio evidenzia come, al diminuire del gap $g$, il numero di condizionamento della matrice di sistema cresca esponenzialmente (raggiungendo $\sim 10^{40}$ per $g/a=0.5$). Il solver proposto gestisce con successo questi regimi in cui i metodi standard a doppia precisione falliscono, mantenendo un'alta accuratezza attraverso l'estensione della precisione e l'aritmetica esatta.

Significatività e Rivendicazioni
Il documento sostiene di fornire una soluzione ad alta accuratezza e verificata numericamente per una specifica ma importante classe di problemi di scattering. La sua significatività risiede nel fatto che:

1. Gestione dell'Ill-Condizionamento: Supera l'instabilità numerica inerente alla modellazione di strutture risonanti e densamente impacchettate a scala sub-lunghezza d'onda, un limite noto dei precedenti metodi semi-analitici.
2. Capacità di Benchmarking: Gli autori dichiarano esplicitamente che il metodo non è inteso a sostituire i solver universali (come FEM o FDTD) per geometrie arbitrarie. Al contrario, funge da potente strumento complementare per:
   • Studi parametrici di ensemble di cilindri.
   • Generazione di soluzioni di benchmark per validare altri schemi numerici universali.
   • Design inverso e screening preliminare nella nanofotonica (ad esempio, SERS, potenziamento della fluorescenza).
3. Controllo Analitico: Il metodo soddisfa esattamente la condizione di radiazione di Sommerfeld, evita errori di discretizzazione spaziale (staircasing) e calcola i campi vicini e lontani dallo stesso set di coefficienti senza interpolazione.

Il lavoro sottolinea che, sebbene l'ambito di applicabilità sia limitato ai cilindri circolari paralleli, l'approccio offre un livello di precisione e di intuizione fisica che è difficile da raggiungere con i metodi numerici basati su griglia, particolarmente nelle critiche regioni di near-field delle nanostrutture plasmoniche.