Quantized plasmon modes for metallic nanoparticles of arbitrary shape with a generic dielectric function

Definizione del Problema
Gli attuali approcci multiscala per la modellazione di nanoparticelle (NP) metalliche plasmoniche accoppiate con molecole si basano tipicamente sull'elettrodinamica classica per la NP e sulla chimica quantistica ab initio per la molecola. Sebbene efficaci per l'accoppiamento debole, questa descrizione classica decade nei regimi di accoppiamento forte plasmone-molecola o di forti campi di pilotaggio (driving fields), dove una descrizione quantizzata della risposta della NP è necessaria per catturare la corretta dinamica del sistema. Inoltre, i metodi di quantizzazione esistenti spesso si affidano a funzioni dielettriche semplificate di tipo Drude o Drude-Lorentz (DL). Questi modelli analitici non riescono a descrivere accuratamente metalli reali come l'argento e l'oro, dove molteplici transizioni di interbanda si sovrappongono alle risonanze plasmoniche degli elettroni liberi, portando a imprecisioni nella risposta dielettrica reale. È necessario un quadro di quantizzazione che possa ospitare NP di forma arbitraria e funzioni dielettriche realistiche e dipendenti dalla frequenza derivate da dati sperimentali.

Metodologia
Gli autori introducono un approccio efficace per quantizzare la risposta elettromagnetica di NP metalliche di forma arbitraria utilizzando un framework di Metodo degli Elementi di Confine (BEM) nel limite quasi-statico, specificamente il Modello del Continuo Polarizzabile - Nanoparticella (PCM-NP).

1. Fondazione Classica: Il metodo inizia con l'equazione della risposta classica PCM-NP, dove la densità di carica superficiale indotta da un campo esterno viene risolta numericamente tramite tessellazione superficiale. La risposta ottica del metallo è definita da una funzione dielettrica generica, $\epsilon(\omega)$, adattata ai dati sperimentali.
2. Fitting Dielettrico: La funzione dipendente dalla frequenza $f(\omega) = (\epsilon(\omega)-1)/(\epsilon(\omega)+1)$ viene adattata a una somma di $N$ poli di tipo Drude-Lorentz. Ciò consente di decomporre le cariche di polarizzazione in componenti dipendenti dai poli.
3. Mappatura Quantistica: Le equazioni classiche dell'oscillatore armonico smorzato e accoppiato derivate dalla decomposizione dei poli vengono mappate sull'equazione di risposta lineare di un sistema quantistico (analogamente alla teoria del funzionale della densità dipendente dal tempo). Separando i termini risonanti e anti-risonanti e utilizzando l'espansione degli autovettori delle matrici di Calderon dipendenti dalla geometria, gli autori derivano un insieme di oscillatori quantistici accoppiati e dissipativi.
4. Costruzione dell'Hamiltoniana: Questa mappatura produce un'Hamiltoniana della NP quantizzata ($\hat{H}_{NP}$) composta da modi plasmonici discreti con specifiche energie di eccitazione e tassi di decadimento. Gli elementi di transizione per gli operatori di carica superficiale vengono identificati, permettendo il calcolo del tensore di polarizzabilità dipolo-dipolo in forma spettrale.
5. Accoppiamento Plexitonico: I modi della NP quantizzata sono accoppiati a una descrizione chimica-quantistica molecolare tramite un'Hamiltoniana totale ($\hat{H}_{tot}$). Questa include gli effetti di polarizzazione dello stato fondamentale e dello stato eccitato e termini di accoppiamento di tipo Jaynes-Cummings tra le transizioni molecolari e i modi plasmonici.

Risultati Chiave

• Validazione della Quantizzazione: Il metodo è stato testato su NP di argento e oro ellissoidali. Si è dimostrato che la parte immaginaria della polarizzabilità dipolo-dipolo calcolata utilizzando i derivati modi quantistici (Q-PCM-NP) riproduce esattamente i risultati ottenuti dalle equazioni di Maxwell macroscopiche classiche (PCM-NP), validando la procedura di quantizzazione.
• Superiorità rispetto ai Modelli Drude-Lorentz: Nelle simulazioni di una NP di argento ellissoidale accoppiata a una molecola di chinolone, l'approccio della funzione dielettrica generica ha catturato con successo gli effetti di accoppiamento coinvolgenti molteplici risonanze plasmoniche. Al contrario, un modello DL a singolo oscillatore semplificato, tarato per corrispondere al picco plasmonico primario, non è riuscito a descrivere l'accoppiamento con le transizioni molecolari a energia più elevata (intorno a 4,5 eV). Il modello DL mostrava un mixing trascurabile, mentre l'approccio generico ha rivelato un accoppiamento e un mixing significativi, coerenti con la presenza di un secondo modo plasmonico.
• Robustezza: Il metodo è stato confermato come robusto attraverso diverse forme di NP (testate su sfere ed ellissoidi) e non richiede procedure di fitting specifiche per il sistema per quanto riguarda la quantizzazione stessa, poiché la quantizzazione deriva naturalmente dalla risposta dielettrica e dalla geometria.

Significatività e Rivendicazioni
L'articolo sostiene di fornire un framework efficace e generale per la quantizzazione delle risposte plasmoniche in nanostrutture di forma arbitraria utilizzando funzioni dielettriche realistiche e derivate sperimentalmente. Gli autori affermano che questa metodologia:

• Recupera correttamente la polarizzazione della risposta lineare classica, consentendo al contempo una descrizione quantistica necessaria per i regimi di accoppiamento forte e di forte pilotaggio.
• Supera i limiti dei modelli dielettrici analitici (come il semplice Drude) che non possono catturare le complesse transizioni di interbanda nei metalli reali.
• Apre la strada a una modellazione accurata di sistemi plexitonici in cui l'forte accoppiamento plasmone-molecola richiede un trattamento completamente quantizzato sia dell'emittente che dell'ambiente plasmonico.
• Offre una descrizione più robusta e trasferibile rispetto agli approcci basati su DL, che spesso richiedono una laboriosa taratura dei parametri specifica per il sistema per riprodurre specifici tratti spettrali.

Il lavoro stabilisce che la rappresentazione della risposta della NP come un insieme di modi discreti, accoppiati e dissipativi è una caratteristica generale delle nanostrutture all'interno del framework quasi-statico, facilitando l'integrazione di proprietà metalliche realistiche nelle simulazioni di elettrodinamica quantistica-chimica quantistica.