Simulazione di filtri ottici mediante il metodo della matrice di trasferimentoSimulazione di filtri ottici mediante il metodo della matrice di trasferimento
In questo tutorial utilizzeremo il risolutore della matrice di trasferimento ottica di OghmaNano per simulare la riflessione, la trasmissione e l’interferenza ottica in filtri ottici multistrato a film sottile. Tali strutture sono ampiamente utilizzate nei riflettori di Bragg, negli specchi dielettrici, nei rivestimenti antiriflesso, nei display OLED e nelle celle solari a film sottile, dove la risposta ottica è controllata dall’indice di rifrazione e dallo spessore di ciascun strato.In questo tutorial utilizzeremo il risolutore della matrice di trasferimento ottica di OghmaNano per simulare la riflessione, la trasmissione e l’interferenza ottica in filtri ottici multistrato a film sottile. Tali strutture sono ampiamente utilizzate nei riflettori di Bragg, negli specchi dielettrici, nei rivestimenti antiriflesso, nei display OLED e nelle celle solari a film sottile, dove la risposta ottica è controllata dall’indice di rifrazione e dallo spessore di ciascun strato.
1. Premessa1. Premessa
Il metodo della matrice di trasferimento (TMM) è ampiamente utilizzato per simulare strutture ottiche multistrato, come le microcavità OLED, le celle solari a film sottile, gli specchi dielettrici, i riflettori di Bragg distribuiti (DBR), i rivestimenti antiriflesso e i filtri ottici selettivi in lunghezza d’onda. In queste strutture, quando la luce incide normalmente, viene parzialmente riflessa e trasmessa ad ogni interfaccia tra materiali con indici di rifrazione diversi, generando interferenza ottica tra le onde elettromagnetiche che si propagano in avanti e all’indietro. Il metodo della matrice di trasferimento modella questa interferenza in modo rigoroso, consentendo di calcolare con precisione la riflessione, la trasmissione e l’assorbimento ottico dell’intera struttura multistrato. La risposta ottica delle strutture multistrato è determinata dall’interferenza tra le onde elettromagnetiche riflesse e trasmesse all’interno della struttura. Una rappresentazione schematica della struttura multistrato e dei campi ottici che si propagano è mostrata nella figura ??.Il metodo della matrice di trasferimento (TMM) è ampiamente utilizzato per simulare strutture ottiche multistrato, come le microcavità OLED, le celle solari a film sottile, gli specchi dielettrici, i riflettori di Bragg distribuiti (DBR), i rivestimenti antiriflesso e i filtri ottici selettivi in lunghezza d’onda. In queste strutture, quando la luce incide normalmente, viene parzialmente riflessa e trasmessa ad ogni interfaccia tra materiali con indici di rifrazione diversi, generando interferenza ottica tra le onde elettromagnetiche che si propagano in avanti e all’indietro. Il metodo della matrice di trasferimento modella questa interferenza in modo rigoroso, consentendo di calcolare con precisione la riflessione, la trasmissione e l’assorbimento ottico dell’intera struttura multistrato. La risposta ottica delle strutture multistrato è determinata dall’interferenza tra le onde elettromagnetiche riflesse e trasmesse all’interno della struttura. Una rappresentazione schematica della struttura multistrato e dei campi ottici che si propagano è mostrata nella figura ??.
Nell’ambito della formalizzazione della matrice di trasferimento, il campo elettrico in ciascun strato viene espresso come la somma di un’onda che si propaga in avanti e di un’onda che si propaga all’indietro.Nell’ambito della formalizzazione della matrice di trasferimento, il campo elettrico in ciascun strato viene espresso come la somma di un’onda che si propaga in avanti e di un’onda che si propaga all’indietro.
\[ E_i(z)=E_i^{+}e^{-\zeta_i z}+E_i^{-}e^{\zeta_i z} \]dove \(E_i^{+}\) e \(E_i^{-}\) rappresentano rispettivamente l’ampiezza del campo che si propaga in avanti e quella che si propaga all’indietro, edove \(E_i^{+}\) e \(E_i^{-}\) rappresentano rispettivamente l’ampiezza del campo che si propaga in avanti e quella che si propaga all’indietro, e
\[ \zeta_i=\frac{2\pi}{\lambda}\bar{n}_i \]dove \(\bar{n}=n+j\kappa\) rappresenta l’indice di rifrazione complesso del materiale. Applicando le equazioni di Maxwell e le condizioni al contorno elettromagnetiche a ciascuna interfaccia, si ottengono i coefficienti di riflessione e trasmissione di Fresnel.dove \(\bar{n}=n+j\kappa\) rappresenta l’indice di rifrazione complesso del materiale. Applicando le equazioni di Maxwell e le condizioni al contorno elettromagnetiche a ciascuna interfaccia, si ottengono i coefficienti di riflessione e trasmissione di Fresnel.
\[ r_{ij}=\frac{n_i-n_j}{n_i+n_j}, \qquad t_{ij}=\frac{2n_i}{n_i+n_j} \]consentendo di mettere in relazione i campi ottici su entrambi i lati di un’interfaccia tramite una matrice di trasferimento.consentendo di mettere in relazione i campi ottici su entrambi i lati di un’interfaccia tramite una matrice di trasferimento.
\[ \begin{pmatrix} E_i^{+} \\ E_i^{-} \end{pmatrix} = \frac{1}{t_{ij}} \begin{pmatrix} 1 & r_{ij} \\ r_{ij} & 1 \end{pmatrix} \begin{pmatrix} E_j^{+} \\ E_j^{-} \end{pmatrix} \]mentre la propagazione attraverso uno strato di spessore \(d_i\) può essere espressa come:mentre la propagazione attraverso uno strato di spessore \(d_i\) può essere espressa come:
\[ \begin{pmatrix} E_i^{+}(z+d_i) \\ E_i^{-}(z+d_i) \end{pmatrix} = \begin{pmatrix} e^{\zeta_i d_i} & 0 \\ 0 & e^{-\zeta_i d_i} \end{pmatrix} \begin{pmatrix} E_i^{+}(z) \\ E_i^{-}(z) \end{pmatrix} \]Pertanto, la struttura ottica multistrato completa si ottiene moltiplicando tra loro le matrici di interfaccia e di propagazione di ogni singolo strato.Pertanto, la struttura ottica multistrato completa si ottiene moltiplicando tra loro le matrici di interfaccia e di propagazione di ogni singolo strato.
\[ M=M_1M_2M_3 \cdots M_n \]a partire dalla quale è possibile calcolare direttamente la potenza ottica riflessa, trasmessa e assorbita. Questo approccio è efficiente dal punto di vista computazionale, numericamente stabile e richiede solo semplici operazioni di moltiplicazione matriciale, il che lo rende particolarmente adatto per la simulazione di filtri ottici multistrato e di strutture a film sottili che producono fenomeni di interferenza.a partire dalla quale è possibile calcolare direttamente la potenza ottica riflessa, trasmessa e assorbita. Questo approccio è efficiente dal punto di vista computazionale, numericamente stabile e richiede solo semplici operazioni di moltiplicazione matriciale, il che lo rende particolarmente adatto per la simulazione di filtri ottici multistrato e di strutture a film sottili che producono fenomeni di interferenza.
Modificando lo spessore e gli indici di rifrazione dei diversi strati, è possibile progettare filtri ottici con proprietà spettrali personalizzate. Ad esempio, un singolo strato dielettrico con spessore pari a un quarto della lunghezza d’onda può ridurre la riflessione superficiale a una determinata lunghezza d’onda, mentre l’alternanza di strati con indice di rifrazione elevato e basso crea un riflettore di Bragg con una banda di blocco ottico molto marcata.Modificando lo spessore e gli indici di rifrazione dei diversi strati, è possibile progettare filtri ottici con proprietà spettrali personalizzate. Ad esempio, un singolo strato dielettrico con spessore pari a un quarto della lunghezza d’onda può ridurre la riflessione superficiale a una determinata lunghezza d’onda, mentre l’alternanza di strati con indice di rifrazione elevato e basso crea un riflettore di Bragg con una banda di blocco ottico molto marcata.
2. Primi passi2. Primi passi
Per avviare la prima simulazione di calcolo di un filtro ottico, apri la finestra «Nuova simulazione» dalla scheda «File» nel menu principale. Fai doppio clic sull’esempio di «Filtro ottico» (vedi ??) e salva la simulazione in una cartella sul tuo disco rigido. A questo punto, vedrai la finestra principale (vedi ??) con una struttura multistrato composta da circa dieci strati alternati. Clicca sul pulsante «Esegui simulazione» (play) per calcolare gli spettri; una volta completata, i risultati appariranno nella scheda «Output» (vedi ??).Per avviare la prima simulazione di calcolo di un filtro ottico, apri la finestra «Nuova simulazione» dalla scheda «File» nel menu principale. Fai doppio clic sull’esempio di «Filtro ottico» (vedi ??) e salva la simulazione in una cartella sul tuo disco rigido. A questo punto, vedrai la finestra principale (vedi ??) con una struttura multistrato composta da circa dieci strati alternati. Clicca sul pulsante «Esegui simulazione» (play) per calcolare gli spettri; una volta completata, i risultati appariranno nella scheda «Output» (vedi ??).
3. Analisi dei risultati.3. Analisi dei risultati.
Dopo aver eseguito la simulazione, fare doppio clic su «Uscita ottica» nel menu ??. Si aprirà l’editor di simulazione ottica. L’editor contiene diverse schede. La prima scheda, «Distribuzione dei fotoni», viene visualizzata automaticamente (vedere ??). In questa scheda viene visualizzata la densità dei fotoni all’interno della cavità ed è chiaramente visibile la struttura a strati del filtro, rappresentata da bande verticali. La seconda scheda, «Distribuzione dei fotoni assorbiti» (vedere ??), mostra dove vengono assorbiti i fotoni. In questo esempio, l’assorbimento è debole ma non nullo, perché uno dei materiali è stato impostato con un piccolo coefficiente di assorbimento \(\alpha\). Infine, la scheda «Luce riflessa» (vedere ??) visualizza lo spettro di riflettanza. Il risultato dimostra una forte riflessione tra circa 500 nm e 800 nm, mentre la luce al di fuori di questa banda viene trasmessa in misura variabile, caratteristica di un filtro di tipo Bragg.Dopo aver eseguito la simulazione, fare doppio clic su «Uscita ottica» nel menu ??. Si aprirà l’editor di simulazione ottica. L’editor contiene diverse schede. La prima scheda, «Distribuzione dei fotoni», viene visualizzata automaticamente (vedere ??). In questa scheda viene visualizzata la densità dei fotoni all’interno della cavità ed è chiaramente visibile la struttura a strati del filtro, rappresentata da bande verticali. La seconda scheda, «Distribuzione dei fotoni assorbiti» (vedere ??), mostra dove vengono assorbiti i fotoni. In questo esempio, l’assorbimento è debole ma non nullo, perché uno dei materiali è stato impostato con un piccolo coefficiente di assorbimento \(\alpha\). Infine, la scheda «Luce riflessa» (vedere ??) visualizza lo spettro di riflettanza. Il risultato dimostra una forte riflessione tra circa 500 nm e 800 nm, mentre la luce al di fuori di questa banda viene trasmessa in misura variabile, caratteristica di un filtro di tipo Bragg.
Lo spettro della luce trasmessa conferma il comportamento di blocco di banda del filtro. La luce con lunghezza d’onda compresa tra circa 300 e 500 nm viene trasmessa in modo efficace, mentre le lunghezze d’onda nell’intervallo 500–800 nm vengono bloccate in modo significativo. A lunghezze d’onda maggiori, superiori a 800 nm, si osserva una parziale ripresa della trasmissione, il che dimostra la natura multi-banda della risposta del filtro.Lo spettro della luce trasmessa conferma il comportamento di blocco di banda del filtro. La luce con lunghezza d’onda compresa tra circa 300 e 500 nm viene trasmessa in modo efficace, mentre le lunghezze d’onda nell’intervallo 500–800 nm vengono bloccate in modo significativo. A lunghezze d’onda maggiori, superiori a 800 nm, si osserva una parziale ripresa della trasmissione, il che dimostra la natura multi-banda della risposta del filtro.
4. Modifica del filtro4. Modifica del filtro
Per esaminare o modificare lo stack, aprire la scheda «Struttura dispositivo» nella finestra principale e fare clic su «Editor livelli» (vedere ??). L’editor elenca ogni livello del dispositivo, indicandone lo spessore, il materiale ottico e le impostazioni. È possibile modificare direttamente lo spessore dei livelli nella tabella, cambiare i materiali, aggiungere o rimuovere livelli e riordinarli secondo le esigenze del progetto del filtro.Per esaminare o modificare lo stack, aprire la scheda «Struttura dispositivo» nella finestra principale e fare clic su «Editor livelli» (vedere ??). L’editor elenca ogni livello del dispositivo, indicandone lo spessore, il materiale ottico e le impostazioni. È possibile modificare direttamente lo spessore dei livelli nella tabella, cambiare i materiali, aggiungere o rimuovere livelli e riordinarli secondo le esigenze del progetto del filtro.
5. Misurazione dello spessore del filtro ottico mediante scansione.5. Misurazione dello spessore del filtro ottico mediante scansione.
Successivamente, analizzeremo la possibilità di eseguire automaticamente la scansione delle larghezze dei filtri. Questo vi permetterà di progettare filtri, eseguirne la scansione automatica e studiare come la struttura del filtro influisce sia sulla trasmissione che sulla riflessione. Eseguendo scansioni sistematiche, potrete generare rapidamente una panoramica di diversi progetti di filtri e scegliere quello più adatto alla vostra applicazione.Successivamente, analizzeremo la possibilità di eseguire automaticamente la scansione delle larghezze dei filtri. Questo vi permetterà di progettare filtri, eseguirne la scansione automatica e studiare come la struttura del filtro influisce sia sulla trasmissione che sulla riflessione. Eseguendo scansioni sistematiche, potrete generare rapidamente una panoramica di diversi progetti di filtri e scegliere quello più adatto alla vostra applicazione.
Nella finestra principale, individuate la scheda «Automazione» e fate clic su «Analisi parametri» (vedere ??). Si aprirà la finestra «Analisi parametri», che dovrebbe già contenere una nuova voce denominata «nuova». Fate doppio clic su questa voce e si aprirà la finestra di configurazione dell’analisi (vedere ??). Infine, fate clic su «Avvia analisi» per iniziare il calcolo. L’analisi potrebbe richiedere un po’ di tempo per essere completata, a seconda del numero di valori dei parametri selezionati.Nella finestra principale, individuate la scheda «Automazione» e fate clic su «Analisi parametri» (vedere ??). Si aprirà la finestra «Analisi parametri», che dovrebbe già contenere una nuova voce denominata «nuova». Fate doppio clic su questa voce e si aprirà la finestra di configurazione dell’analisi (vedere ??). Infine, fate clic su «Avvia analisi» per iniziare il calcolo. L’analisi potrebbe richiedere un po’ di tempo per essere completata, a seconda del numero di valori dei parametri selezionati.
Facendo clic su «Esegui scansione», si avvia il piccolo programma mostrato in ??. Spiegheremo questo programma in seguito, ma in sostanza esso modifica lo spessore dello strato con l’indice di rifrazione più elevato, facendo variare tale valore all’interno di un intervallo predefinito. Aprire la scheda «Risultati» nella finestra «Scansione parametri» (vedere ??): si vedranno quattro cartelle, ciascuna delle quali corrisponde a uno degli spessori scansionati. Ogni cartella contiene una simulazione completa con i file abituali, differenziandosi solo per lo spessore dello strato in questione.Facendo clic su «Esegui scansione», si avvia il piccolo programma mostrato in ??. Spiegheremo questo programma in seguito, ma in sostanza esso modifica lo spessore dello strato con l’indice di rifrazione più elevato, facendo variare tale valore all’interno di un intervallo predefinito. Aprire la scheda «Risultati» nella finestra «Scansione parametri» (vedere ??): si vedranno quattro cartelle, ciascuna delle quali corrisponde a uno degli spessori scansionati. Ogni cartella contiene una simulazione completa con i file abituali, differenziandosi solo per lo spessore dello strato in questione.
Nella cartella principale della scansione troverete anche icone speciali «multi-curva» (file CSV con un simbolo multilinea). Queste raccolgono le curve corrispondenti di tutte le sottosimulazioni. Fate doppio clic su «optical_output» per aprirlo (vedere ??); vedrete reflect.csv e transmit.csv. Aprendo reflect.csv, verrà visualizzato un grafico che mostra la riflettanza di tutti gli spessori scansionati, come illustrato in ??.Nella cartella principale della scansione troverete anche icone speciali «multi-curva» (file CSV con un simbolo multilinea). Queste raccolgono le curve corrispondenti di tutte le sottosimulazioni. Fate doppio clic su «optical_output» per aprirlo (vedere ??); vedrete reflect.csv e transmit.csv. Aprendo reflect.csv, verrà visualizzato un grafico che mostra la riflettanza di tutti gli spessori scansionati, come illustrato in ??.
reflect.csv e transmit.csv.All’interno di una sottocartella di scansione: dati aggregati reflect.csv e transmit.csv.
6. Comprendere il funzionamento del programma di scansione dei parametri6. Comprendere il funzionamento del programma di scansione dei parametri
Se si osserva nuovamente ??, si noterà che sono elencate cinque righe. Ogni riga specifica un parametro da analizzare. In questo esempio, il parametro è lo spessore (dy) di specifici strati nello strato epitassiale. Le voci sono epitaxy.layer0.dy, epitaxy.layer2.dy, epitaxy.layer4.dy, epitaxy.layer6.dy e epitaxy.layer8.dy. Queste righe corrispondono agli strati con un elevato indice di rifrazione nella struttura del dispositivo.Se si osserva nuovamente ??, si noterà che sono elencate cinque righe. Ogni riga specifica un parametro da analizzare. In questo esempio, il parametro è lo spessore (dy) di specifici strati nello strato epitassiale. Le voci sono epitaxy.layer0.dy, epitaxy.layer2.dy, epitaxy.layer4.dy, epitaxy.layer6.dy e epitaxy.layer8.dy. Queste righe corrispondono agli strati con un elevato indice di rifrazione nella struttura del dispositivo.
Nella prima riga è possibile visualizzare un elenco di valori di spessore: 1.0e-7, 2e-7, 4e-7, 8e-7. Questi sono i valori che la scansione assegnerà a layer0. L’operazione per questa riga è impostata su «scansione», il che significa che il programma varierà sistematicamente lo spessore di layer0 in base a questi valori specificati.Nella prima riga è possibile visualizzare un elenco di valori di spessore: 1.0e-7, 2e-7, 4e-7, 8e-7. Questi sono i valori che la scansione assegnerà a layer0. L’operazione per questa riga è impostata su «scansione», il che significa che il programma varierà sistematicamente lo spessore di layer0 in base a questi valori specificati.
Per gli altri strati (layer2, layer4, layer6 e layer8), l’operazione è collegata a epitaxy.layer0.dy. Nella colonna dei valori, questo appare come duplicato. Questo comando indica al programma di copiare il valore attualmente impostato per layer0 e di applicarlo agli altri strati con indice elevato. In pratica, ogni volta che layer0 cambia, gli altri strati con indice elevato vengono aggiornati automaticamente alla stessa dimensione.Per gli altri strati (layer2, layer4, layer6 e layer8), l’operazione è collegata a epitaxy.layer0.dy. Nella colonna dei valori, questo appare come duplicato. Questo comando indica al programma di copiare il valore attualmente impostato per layer0 e di applicarlo agli altri strati con indice elevato. In pratica, ogni volta che layer0 cambia, gli altri strati con indice elevato vengono aggiornati automaticamente alla stessa dimensione.
In sintesi, questo programma analizza una serie di valori di spessore per uno strato ad alto indice di rifrazione e quindi replica tale spessore in tutti gli altri strati ad alto indice di rifrazione presenti nello stack. Il risolutore viene eseguito per ciascun caso, consentendo di analizzare come le prestazioni del filtro variano in funzione dello spessore degli strati.In sintesi, questo programma analizza una serie di valori di spessore per uno strato ad alto indice di rifrazione e quindi replica tale spessore in tutti gli altri strati ad alto indice di rifrazione presenti nello stack. Il risolutore viene eseguito per ciascun caso, consentendo di analizzare come le prestazioni del filtro variano in funzione dello spessore degli strati.
7. Riepilogo7. Riepilogo
In questa parte del tutorial, hai imparato come automatizzare la progettazione di filtri ottici utilizzando la scansione dei parametri. Hai modificato sistematicamente lo spessore degli strati, hai replicato i parametri su più strati e hai analizzato come queste modifiche influiscono sulla trasmissione e sulla riflessione. Grazie a questi strumenti, ora puoi valutare rapidamente numerosi progetti di filtri e identificare le strutture più adatte alla tua applicazione.In questa parte del tutorial, hai imparato come automatizzare la progettazione di filtri ottici utilizzando la scansione dei parametri. Hai modificato sistematicamente lo spessore degli strati, hai replicato i parametri su più strati e hai analizzato come queste modifiche influiscono sulla trasmissione e sulla riflessione. Grazie a questi strumenti, ora puoi valutare rapidamente numerosi progetti di filtri e identificare le strutture più adatte alla tua applicazione.