Simulazione del diodo a giunzione PN in InP: deriva-diffusione, raddrizzamento e ricombinazione delle portanti di carica.Simulazione del diodo a giunzione PN in InP: deriva-diffusione, raddrizzamento e ricombinazione delle portanti di carica.

1. Introduzione: una giunzione raddrizzatrice in InP1. Introduzione: una giunzione raddrizzatrice in InP

Il fosfuro di indio (InP) è un semiconduttore di tipo III-V con banda proibita diretta, utilizzato come base per dispositivi elettronici ad alta velocità e tecnologie fotoniche integrate. Un diodo a giunzione PN è una delle sue strutture di dispositivo più semplici e più utili: una regione di tipo p incontra una regione di tipo n, creando un campo elettrico intrinseco e una barriera di potenziale che controlla il movimento di elettroni e lacune. Il diodo conduce in modo efficiente in polarizzazione diretta e sopprime l'iniezione di portatori in polarizzazione inversa, producendo una rettifica elettrica.Il fosfuro di indio (InP) è un semiconduttore di tipo III-V con banda proibita diretta, utilizzato come base per dispositivi elettronici ad alta velocità e tecnologie fotoniche integrate. Un diodo a giunzione PN è una delle sue strutture di dispositivo più semplici e più utili: una regione di tipo p incontra una regione di tipo n, creando un campo elettrico intrinseco e una barriera di potenziale che controlla il movimento di elettroni e lacune. Il diodo conduce in modo efficiente in polarizzazione diretta e sopprime l'iniezione di portatori in polarizzazione inversa, producendo una rettifica elettrica.

L’InP presenta una banda proibita a temperatura ambiente di circa 1,34 eV e una mobilità degli elettroni relativamente elevata. Grazie a queste proprietà, il trasporto nelle giunzioni di InP è rilevante per l’elettronica dei semiconduttori compositi e per la comprensione delle interfacce dei semiconduttori utilizzate nelle piattaforme di dispositivi basate su InP. L’esempio presentato in questo tutorial illustra il funzionamento elettrico di una giunzione omogiunzione p-InP/n-InP compatta con uno strato di contatto n⁺.L’InP presenta una banda proibita a temperatura ambiente di circa 1,34 eV e una mobilità degli elettroni relativamente elevata. Grazie a queste proprietà, il trasporto nelle giunzioni di InP è rilevante per l’elettronica dei semiconduttori compositi e per la comprensione delle interfacce dei semiconduttori utilizzate nelle piattaforme di dispositivi basate su InP. L’esempio presentato in questo tutorial illustra il funzionamento elettrico di una giunzione omogiunzione p-InP/n-InP compatta con uno strato di contatto n⁺.

Invece di adattare la curva J–V a una semplice equazione di diodo ideale, OghmaNano risolve le equazioni di Poisson e di continuità per la deriva-diffusione, considerate come un sistema accoppiato. I dati ottenuti rivelano sia la corrente ai terminali sia la distribuzione interna delle energie delle bande di conduzione e di valenza, le densità di elettroni e lacune, i livelli quasi-Fermi e i processi di ricombinazione. Questo ci permette di spiegare perché il diodo svolge la funzione di raddrizzatore, anziché limitarci a osservare che lo fa.Invece di adattare la curva J–V a una semplice equazione di diodo ideale, OghmaNano risolve le equazioni di Poisson e di continuità per la deriva-diffusione, considerate come un sistema accoppiato. I dati ottenuti rivelano sia la corrente ai terminali sia la distribuzione interna delle energie delle bande di conduzione e di valenza, le densità di elettroni e lacune, i livelli quasi-Fermi e i processi di ricombinazione. Questo ci permette di spiegare perché il diodo svolge la funzione di raddrizzatore, anziché limitarci a osservare che lo fa.

La finestra principale della simulazione, mostrata in ??, visualizza il dispositivo InP a strati verticali. Sebbene la geometria sia rappresentata in tre dimensioni, il problema elettrico risolto in questo caso è unidimensionale lungo la direzione di crescita.La finestra principale della simulazione, mostrata in ??, visualizza il dispositivo InP a strati verticali. Sebbene la geometria sia rappresentata in tre dimensioni, il problema elettrico risolto in questo caso è unidimensionale lungo la direzione di crescita.

2. Creazione di una nuova simulazione2. Creazione di una nuova simulazione

Nella finestra principale di OghmaNano, selezionare «Nuova simulazione». Si aprirà la finestra di selezione delle simulazioni, con le categorie di esempi disponibili (??). Aprire le dimostrazioni III–V/InP e scegliere l’esempio relativo al diodo InP PN (??). La posizione esatta dell’esempio potrebbe variare a seconda della versione di OghmaNano.Nella finestra principale di OghmaNano, selezionare «Nuova simulazione». Si aprirà la finestra di selezione delle simulazioni, con le categorie di esempi disponibili (??). Aprire le dimostrazioni III–V/InP e scegliere l’esempio relativo al diodo InP PN (??). La posizione esatta dell’esempio potrebbe variare a seconda della versione di OghmaNano.

Dopo aver aperto l’esempio, selezionate la scheda «Struttura del dispositivo» e aprite l’editor dei livelli. Verificate i nomi dei livelli, i loro spessori e i modelli dei materiali elettrici associati prima di avviare la simulazione. Le sezioni seguenti spiegano in che modo ciascuna regione contribuisce alla risposta complessiva del diodo.Dopo aver aperto l’esempio, selezionate la scheda «Struttura del dispositivo» e aprite l’editor dei livelli. Verificate i nomi dei livelli, i loro spessori e i modelli dei materiali elettrici associati prima di avviare la simulazione. Le sezioni seguenti spiegano in che modo ciascuna regione contribuisce alla risposta complessiva del diodo.

3. Analisi della struttura del dispositivo, dei contatti e del processo di drogaggio.3. Analisi della struttura del dispositivo, dei contatti e del processo di drogaggio.

3.1 Struttura dello strato semiconduttore3.1 Struttura dello strato semiconduttore

Il dispositivo elettrico contiene tre strati di InP, con uno spessore totale del materiale semiconduttore pari a 700 nm. Gli strati p-InP e n-InP formano la giunzione omogiunzione raddrizzante, mentre la regione n⁺-InP costituisce una regione catodica con una maggiore concentrazione di impurità. La disposizione degli strati è visibile nell'editor (??).Il dispositivo elettrico contiene tre strati di InP, con uno spessore totale del materiale semiconduttore pari a 700 nm. Gli strati p-InP e n-InP formano la giunzione omogiunzione raddrizzante, mentre la regione n⁺-InP costituisce una regione catodica con una maggiore concentrazione di impurità. La disposizione degli strati è visibile nell'editor (??).

La regione di aria di 200 nm visualizzata nell’editor di struttura fa parte della rappresentazione grafica e della geometria ottica circostante, e non dello strato elettrico del semiconduttore. Pertanto, il suo spessore non viene considerato nel calcolo dello spessore totale del dispositivo elettrico, che è di 700 nm.La regione di aria di 200 nm visualizzata nell’editor di struttura fa parte della rappresentazione grafica e della geometria ottica circostante, e non dello strato elettrico del semiconduttore. Pertanto, il suo spessore non viene considerato nel calcolo dello spessore totale del dispositivo elettrico, che è di 700 nm.

3.2 Drogaggio ed elettrostatica intrinseca3.2 Drogaggio ed elettrostatica intrinseca

Aprire la sezione «Parametri elettrici» → «Drogaggio/Ioni» per visualizzare le concentrazioni dei donatori e degli accettori (??). La variazione netta del drogaggio cambia segno all’interfaccia tra il materiale p-InP e il materiale n-InP, a circa 250 nm dalla parte superiore dello strato semiconduttore attivo. La seconda variazione, tra il materiale n-InP e il materiale n⁺-InP, modifica la densità degli elettroni in direzione del catodo, ma non crea una nuova giunzione PN.Aprire la sezione «Parametri elettrici» → «Drogaggio/Ioni» per visualizzare le concentrazioni dei donatori e degli accettori (??). La variazione netta del drogaggio cambia segno all’interfaccia tra il materiale p-InP e il materiale n-InP, a circa 250 nm dalla parte superiore dello strato semiconduttore attivo. La seconda variazione, tra il materiale n-InP e il materiale n⁺-InP, modifica la densità degli elettroni in direzione del catodo, ma non crea una nuova giunzione PN.

All’equilibrio, gli elettroni tendono a diffondersi dalla regione di tipo n alla regione di tipo p, mentre le lacune si diffondono nella direzione opposta. Gli atomi droganti ionizzati non compensati generano un campo che si oppone all’ulteriore diffusione. La regione di svuotamento e la barriera di potenziale interna che si formano costituiscono la base della risposta di raddrizzamento.All’equilibrio, gli elettroni tendono a diffondersi dalla regione di tipo n alla regione di tipo p, mentre le lacune si diffondono nella direzione opposta. Gli atomi droganti ionizzati non compensati generano un campo che si oppone all’ulteriore diffusione. La regione di svuotamento e la barriera di potenziale interna che si formano costituiscono la base della risposta di raddrizzamento.

In questa formula, \(\phi\) rappresenta il potenziale elettrostatico, \(\varepsilon\) la permittività dielettrica e \(\rho\) la densità di carica locale. Gli atomi droganti ionizzati e le cariche mobili determinano il campo elettrico locale, mentre la tensione applicata modifica la barriera effettiva per l’iniezione di cariche.In questa formula, \(\phi\) rappresenta il potenziale elettrostatico, \(\varepsilon\) la permittività dielettrica e \(\rho\) la densità di carica locale. Gli atomi droganti ionizzati e le cariche mobili determinano il campo elettrico locale, mentre la tensione applicata modifica la barriera effettiva per l’iniezione di cariche.

3.3 Contatti elettrici3.3 Contatti elettrici

Il modello prevede una barriera ohmica con prevalenza di lacune all’anodo e una barriera ohmica con prevalenza di elettroni al catodo. Queste condizioni al contorno sono adatte per studiare la risposta della giunzione senza introdurre una barriera di iniezione aggiuntiva nei punti di contatto. Nei dispositivi fabbricati, la metallizzazione e la lavorazione dei contatti determinano quanto i contatti si avvicinino al comportamento ohmico ideale.Il modello prevede una barriera ohmica con prevalenza di lacune all’anodo e una barriera ohmica con prevalenza di elettroni al catodo. Queste condizioni al contorno sono adatte per studiare la risposta della giunzione senza introdurre una barriera di iniezione aggiuntiva nei punti di contatto. Nei dispositivi fabbricati, la metallizzazione e la lavorazione dei contatti determinano quanto i contatti si avvicinino al comportamento ohmico ideale.

4. Modelli dei materiali e parametri di trasporto4. Modelli dei materiali e parametri di trasporto

4.1 Modelli dei materiali Lua: da dove provengono i parametri di InP.4.1 Modelli dei materiali Lua: da dove provengono i parametri di InP.

Una delle caratteristiche importanti di questo esempio di diodo InP è che i parametri elettrici del materiale sono definiti utilizzando modelli di materiale Lua. Ciò stabilisce una chiara connessione tra le proprietà fisiche di un semiconduttore e la soluzione relativa alla deriva e alla diffusione. Le regioni p-InP, n-InP e n⁺-InP fanno tutte riferimento alla stessa definizione di materiale, Inorganic/InP/InP.lua. Il loro diverso comportamento elettrico deriva principalmente dal profilo di drogaggio, piuttosto che dall’assegnazione di tre serie distinte di costanti del materiale.Una delle caratteristiche importanti di questo esempio di diodo InP è che i parametri elettrici del materiale sono definiti utilizzando modelli di materiale Lua. Ciò stabilisce una chiara connessione tra le proprietà fisiche di un semiconduttore e la soluzione relativa alla deriva e alla diffusione. Le regioni p-InP, n-InP e n⁺-InP fanno tutte riferimento alla stessa definizione di materiale, Inorganic/InP/InP.lua. Il loro diverso comportamento elettrico deriva principalmente dal profilo di drogaggio, piuttosto che dall’assegnazione di tre serie distinte di costanti del materiale.

Per esaminare il modello, aprire la sezione Struttura dispositivo → Parametri elettrici (l’editor dei materiali elettrici) e selezionare uno degli strati di InP. Nella scheda Editor codice viene visualizzato lo script Lua relativo al materiale associato a tale regione; nella scheda Parametri vengono visualizzate le proprietà del materiale calcolate. Queste due sezioni devono essere utilizzate insieme. Il codice definisce come vengono ottenute le proprietà, mentre i valori calcolati mostrano quali sono i dati effettivamente ricevuti dal risolutore elettrico.Per esaminare il modello, aprire la sezione Struttura dispositivo → Parametri elettrici (l’editor dei materiali elettrici) e selezionare uno degli strati di InP. Nella scheda Editor codice viene visualizzato lo script Lua relativo al materiale associato a tale regione; nella scheda Parametri vengono visualizzate le proprietà del materiale calcolate. Queste due sezioni devono essere utilizzate insieme. Il codice definisce come vengono ottenute le proprietà, mentre i valori calcolati mostrano quali sono i dati effettivamente ricevuti dal risolutore elettrico.

Le immagini in ?? e ?? mostrano come accedere all’editor di codice. La configurazione dei livelli della simulazione è la fonte di riferimento per il modello di materiale assegnato: in questo esempio, ogni regione attiva è associata al file Lua di InP. I nomi dei file immagine derivano dalla sequenza di acquisizione degli screenshot e non determinano quale semiconduttore viene simulato.Le immagini in ?? e ?? mostrano come accedere all’editor di codice. La configurazione dei livelli della simulazione è la fonte di riferimento per il modello di materiale assegnato: in questo esempio, ogni regione attiva è associata al file Lua di InP. I nomi dei file immagine derivano dalla sequenza di acquisizione degli screenshot e non determinano quale semiconduttore viene simulato.

Nella simulazione fornita, gli interruttori db_overrides sono impostati su false per i parametri chiave relativi al trasporto e alla banda. Di conseguenza, OghmaNano è configurato per ottenere i valori corrispondenti dal database dei materiali/modello Lua, anziché utilizzare i valori predefiniti memorizzati altrove nel file JSON. Questo aspetto è importante quando si esamina la configurazione: i valori effettivi visualizzati nella scheda Parametri forniscono informazioni più dettagliate rispetto ai valori predefiniti memorizzati in un campo di sovrascrittura inattivo.Nella simulazione fornita, gli interruttori db_overrides sono impostati su false per i parametri chiave relativi al trasporto e alla banda. Di conseguenza, OghmaNano è configurato per ottenere i valori corrispondenti dal database dei materiali/modello Lua, anziché utilizzare i valori predefiniti memorizzati altrove nel file JSON. Questo aspetto è importante quando si esamina la configurazione: i valori effettivi visualizzati nella scheda Parametri forniscono informazioni più dettagliate rispetto ai valori predefiniti memorizzati in un campo di sovrascrittura inattivo.

4.2 Struttura a bande, densità degli stati e mobilità4.2 Struttura a bande, densità degli stati e mobilità

Nella scheda «Parametri valutati», i primi valori da controllare sono il divario di energia \(E_g\) e l’affinità elettronica \(\chi\). Il divario di energia definisce la separazione tra la banda di conduzione e la banda di valenza, mentre l’affinità elettronica indica la posizione del bordo della banda di conduzione rispetto a un riferimento di energia del vuoto. Per l’InP, il divario di energia a temperatura ambiente è approssimativamente \(1.34\,\mathrm{eV}\). Un utile controllo preliminare di coerenza consiste nel verificare se la separazione tra \(E_c\) e \(E_v\) nel diagramma di banda risultante corrisponde al divario di energia valutato.Nella scheda «Parametri valutati», i primi valori da controllare sono il divario di energia \(E_g\) e l’affinità elettronica \(\chi\). Il divario di energia definisce la separazione tra la banda di conduzione e la banda di valenza, mentre l’affinità elettronica indica la posizione del bordo della banda di conduzione rispetto a un riferimento di energia del vuoto. Per l’InP, il divario di energia a temperatura ambiente è approssimativamente \(1.34\,\mathrm{eV}\). Un utile controllo preliminare di coerenza consiste nel verificare se la separazione tra \(E_c\) e \(E_v\) nel diagramma di banda risultante corrisponde al divario di energia valutato.

Le densità di stati effettive delle bande di conduzione e di valenza, \(N_c\) e \(N_v\), influenzano la relazione tra la concentrazione dei portatori di carica e la posizione del livello di quasi-Fermi. Di conseguenza, contribuiscono alle concentrazioni di elettroni e lacune all’equilibrio e alla risposta all’iniezione in condizioni di polarizzazione diretta. La permittività relativa \(\varepsilon_r\) entra nell’equazione di Poisson e influenza il potenziale elettrostatico e la larghezza della zona di svuotamento. Queste quantità devono essere considerate come un insieme coerente, e non come parametri numerici indipendenti da regolare.Le densità di stati effettive delle bande di conduzione e di valenza, \(N_c\) e \(N_v\), influenzano la relazione tra la concentrazione dei portatori di carica e la posizione del livello di quasi-Fermi. Di conseguenza, contribuiscono alle concentrazioni di elettroni e lacune all’equilibrio e alla risposta all’iniezione in condizioni di polarizzazione diretta. La permittività relativa \(\varepsilon_r\) entra nell’equazione di Poisson e influenza il potenziale elettrostatico e la larghezza della zona di svuotamento. Queste quantità devono essere considerate come un insieme coerente, e non come parametri numerici indipendenti da regolare.

Le mobilità degli elettroni e delle lacune, \(\mu_n\) e \(\mu_p\), determinano la risposta alla deriva dei portatori di carica in presenza di un campo elettrico; inoltre, la diffusione dipende dai gradienti di concentrazione dei portatori di carica. Nel modello unidimensionale, le densità di corrente sono espresse come segue:Le mobilità degli elettroni e delle lacune, \(\mu_n\) e \(\mu_p\), determinano la risposta alla deriva dei portatori di carica in presenza di un campo elettrico; inoltre, la diffusione dipende dai gradienti di concentrazione dei portatori di carica. Nel modello unidimensionale, le densità di corrente sono espresse come segue:

In questo caso, \(n\) e \(p\) rappresentano rispettivamente la densità locale degli elettroni e delle lacune, \(E\) è il campo elettrico e \(D_n\) e \(D_p\) sono i coefficienti di diffusione. Le due componenti della corrente possono variare notevolmente all’interno del diodo, anche quando la corrente totale in regime stazionario è approssimativamente indipendente dalla posizione. Per questo motivo, nelle sezioni successive vengono analizzati \(J_n\) e \(J_p\) insieme ai profili di densità delle cariche.In questo caso, \(n\) e \(p\) rappresentano rispettivamente la densità locale degli elettroni e delle lacune, \(E\) è il campo elettrico e \(D_n\) e \(D_p\) sono i coefficienti di diffusione. Le due componenti della corrente possono variare notevolmente all’interno del diodo, anche quando la corrente totale in regime stazionario è approssimativamente indipendente dalla posizione. Per questo motivo, nelle sezioni successive vengono analizzati \(J_n\) e \(J_p\) insieme ai profili di densità delle cariche.

4.3 Come controllare e confrontare le impostazioni dei materiali4.3 Come controllare e confrontare le impostazioni dei materiali

Nell’editor di dispositivi elettrici, selezionate il materiale p-InP ed esaminate il relativo codice Lua e i parametri calcolati. Quindi, passate al materiale n-InP e, infine, al materiale n⁺-InP. Confrontate il modello di materiale assegnato, il gap di banda, l’affinità elettronica, le mobilità, le densità di stati effettive, la permittività e i coefficienti di ricombinazione. Per questa giunzione omogiunzione, viene utilizzato lo stesso modello di base di InP in tutte e tre le regioni; la principale differenza intenzionale tra le regioni è la loro concentrazione di accettori o donatori.Nell’editor di dispositivi elettrici, selezionate il materiale p-InP ed esaminate il relativo codice Lua e i parametri calcolati. Quindi, passate al materiale n-InP e, infine, al materiale n⁺-InP. Confrontate il modello di materiale assegnato, il gap di banda, l’affinità elettronica, le mobilità, le densità di stati effettive, la permittività e i coefficienti di ricombinazione. Per questa giunzione omogiunzione, viene utilizzato lo stesso modello di base di InP in tutte e tre le regioni; la principale differenza intenzionale tra le regioni è la loro concentrazione di accettori o donatori.

È utile distinguere due concetti: il modello dei materiali descrive le proprietà dell’InP, mentre il profilo di drogaggio stabilisce quali portatori di carica prevalgono localmente e dove si trova la giunzione. Questa distinzione diventa particolarmente utile quando si passa da una giunzione omopolare PN a una giunzione eteropolare a base di InP, dove le differenze tra i modelli dei materiali possono a loro volta generare differenze di energia tra le bande di conduzione e di valenza. In questo diodo, la curvatura delle bande di conduzione e di valenza, visibile nelle immagini, riflette principalmente l’elettrostatica attraverso la transizione p/n.È utile distinguere due concetti: il modello dei materiali descrive le proprietà dell’InP, mentre il profilo di drogaggio stabilisce quali portatori di carica prevalgono localmente e dove si trova la giunzione. Questa distinzione diventa particolarmente utile quando si passa da una giunzione omopolare PN a una giunzione eteropolare a base di InP, dove le differenze tra i modelli dei materiali possono a loro volta generare differenze di energia tra le bande di conduzione e di valenza. In questo diodo, la curvatura delle bande di conduzione e di valenza, visibile nelle immagini, riflette principalmente l’elettrostatica attraverso la transizione p/n.

Per ottenere risultati riproducibili, mantieni fisso il modello dei materiali Lua quando modifichi le concentrazioni di drogaggio. Al contrario, se devi analizzare la sensibilità del risultato rispetto a un coefficiente di mobilità o di ricombinazione, modifica intenzionalmente tale parametro del materiale e confronta la curva J–V, le densità di portatori e i profili di corrente interna risultanti. L’editor dei parametri e le visualizzazioni del codice Lua forniscono il collegamento tra le ipotesi di input fisiche e l’output simulato.Per ottenere risultati riproducibili, mantieni fisso il modello dei materiali Lua quando modifichi le concentrazioni di drogaggio. Al contrario, se devi analizzare la sensibilità del risultato rispetto a un coefficiente di mobilità o di ricombinazione, modifica intenzionalmente tale parametro del materiale e confronta la curva J–V, le densità di portatori e i profili di corrente interna risultanti. L’editor dei parametri e le visualizzazioni del codice Lua forniscono il collegamento tra le ipotesi di input fisiche e l’output simulato.

4.4 Parametri di ricombinazione4.4 Parametri di ricombinazione

Le regioni di InP comprendono la ricombinazione di Shockley-Read-Hall (SRH) e un canale di ricombinazione di portatori liberi, come previsto dal modello dei materiali utilizzato. La ricombinazione SRH descrive la cattura di elettroni e lacune facilitata da difetti; la ricombinazione tra portatori liberi è particolarmente rilevante nell’InP a banda proibita diretta, poiché le transizioni elettrone-lacuna possono avvenire in modo radiativo. Nel presente esercizio, entrambi i processi sono utili in quanto contribuiscono alla perdita di portatori e modificano la corrente elettrica.Le regioni di InP comprendono la ricombinazione di Shockley-Read-Hall (SRH) e un canale di ricombinazione di portatori liberi, come previsto dal modello dei materiali utilizzato. La ricombinazione SRH descrive la cattura di elettroni e lacune facilitata da difetti; la ricombinazione tra portatori liberi è particolarmente rilevante nell’InP a banda proibita diretta, poiché le transizioni elettrone-lacuna possono avvenire in modo radiativo. Nel presente esercizio, entrambi i processi sono utili in quanto contribuiscono alla perdita di portatori e modificano la corrente elettrica.

Il coefficiente \(B\), i tempi di ricombinazione nelle trappole \(\tau_n\) e \(\tau_p\) e le densità ausiliarie \(n_1\) e \(p_1\) determinano l’intensità di ciascun processo. L’aspetto fondamentale di questa guida è illustrare come la distribuzione spaziale della ricombinazione cambia quando una tensione di polarizzazione induce un maggior numero di elettroni e lacune nella stessa regione.Il coefficiente \(B\), i tempi di ricombinazione nelle trappole \(\tau_n\) e \(\tau_p\) e le densità ausiliarie \(n_1\) e \(p_1\) determinano l’intensità di ciascun processo. L’aspetto fondamentale di questa guida è illustrare come la distribuzione spaziale della ricombinazione cambia quando una tensione di polarizzazione induce un maggior numero di elettroni e lacune nella stessa regione.

4.5 Rete e soluzione numerica4.5 Rete e soluzione numerica

Il calcolo elettrico utilizza una griglia monodimensionale di 150 punti su una lunghezza di 700 nm, corrispondente a una spaziatura media di circa 4,7 nm. Le equazioni non lineari accoppiate vengono risolte mediante l’iterazione di Newton e il risolutore lineare sparso UMFPACK. Questa griglia rappresenta un utile punto di partenza per l’analisi della giunzione e delle regioni quasi neutre.Il calcolo elettrico utilizza una griglia monodimensionale di 150 punti su una lunghezza di 700 nm, corrispondente a una spaziatura media di circa 4,7 nm. Le equazioni non lineari accoppiate vengono risolte mediante l’iterazione di Newton e il risolutore lineare sparso UMFPACK. Questa griglia rappresenta un utile punto di partenza per l’analisi della giunzione e delle regioni quasi neutre.

L’output include informazioni sulla convergenza e profili dei dispositivi che dipendono dalla polarizzazione. Un buon primo controllo consiste nel verificare che la curva corrente-tensione (J–V) sia regolare e quindi confermare che le correnti di elettroni e lacune rispettino la legge di conservazione della corrente in condizioni di regime stazionario. L’editor della griglia è illustrato in ??.L’output include informazioni sulla convergenza e profili dei dispositivi che dipendono dalla polarizzazione. Un buon primo controllo consiste nel verificare che la curva corrente-tensione (J–V) sia regolare e quindi confermare che le correnti di elettroni e lacune rispettino la legge di conservazione della corrente in condizioni di regime stazionario. L’editor della griglia è illustrato in ??.

5. Esecuzione della simulazione e interpretazione della curva caratteristica corrente-tensione (J–V).5. Esecuzione della simulazione e interpretazione della curva caratteristica corrente-tensione (J–V).

Torna alla finestra principale e seleziona l’opzione «Esegui simulazione». OghmaNano esegue una scansione della tensione applicata e registra la soluzione convergente in ogni punto. Il file di esempio salvato utilizza una scansione da 0 a 3 V con incrementi di 25 mV. La finestra di output registra l’avanzamento del processo di risoluzione durante la scansione, mentre la scheda «Output» consente di accedere ai dati risultanti.Torna alla finestra principale e seleziona l’opzione «Esegui simulazione». OghmaNano esegue una scansione della tensione applicata e registra la soluzione convergente in ogni punto. Il file di esempio salvato utilizza una scansione da 0 a 3 V con incrementi di 25 mV. La finestra di output registra l’avanzamento del processo di risoluzione durante la scansione, mentre la scheda «Output» consente di accedere ai dati risultanti.

Una volta completata la simulazione, i risultati saranno disponibili nella scheda «Output» (??). Aprire jv.csv per visualizzare il grafico della densità di corrente in funzione della tensione applicata. Gli screenshot includono sia una rappresentazione lineare della curva J–V (??) sia una rappresentazione logaritmica della corrente (??).Una volta completata la simulazione, i risultati saranno disponibili nella scheda «Output» (??). Aprire jv.csv per visualizzare il grafico della densità di corrente in funzione della tensione applicata. Gli screenshot includono sia una rappresentazione lineare della curva J–V (??) sia una rappresentazione logaritmica della corrente (??).

5.1 Comprendere le caratteristiche5.1 Comprendere le caratteristiche

A basse tension di alimentazione, la barriera di giunzione intrinseca limita notevolmente l’iniezione di portatori di carica. Aumentando la polarizzazione diretta, tale barriera si riduce, consentendo agli elettroni e alle lacune di attraversare più facilmente il lato opposto, con conseguente aumento della corrente totale. Nel regime intermedio, la corrente è spesso approssimativamente proporzionale alla tensione; a correnti più elevate, i fenomeni di trasporto e le perdite resistive possono modificare la pendenza della curva.A basse tension di alimentazione, la barriera di giunzione intrinseca limita notevolmente l’iniezione di portatori di carica. Aumentando la polarizzazione diretta, tale barriera si riduce, consentendo agli elettroni e alle lacune di attraversare più facilmente il lato opposto, con conseguente aumento della corrente totale. Nel regime intermedio, la corrente è spesso approssimativamente proporzionale alla tensione; a correnti più elevate, i fenomeni di trasporto e le perdite resistive possono modificare la pendenza della curva.

Questa espressione che descrive il comportamento ideale di un diodo fornisce un utile punto di riferimento per l’analisi di un intervallo approssimativamente esponenziale, ma non è una condizione imposta dal modello di deriva-diffusione. Il fattore di idealità \(n\) e il prefattore \(J_0\) devono essere determinati solo all’interno di un intervallo di polarizzazione appropriato; possono variare in funzione della ricombinazione, della resistenza in serie e del comportamento del contatto.Questa espressione che descrive il comportamento ideale di un diodo fornisce un utile punto di riferimento per l’analisi di un intervallo approssimativamente esponenziale, ma non è una condizione imposta dal modello di deriva-diffusione. Il fattore di idealità \(n\) e il prefattore \(J_0\) devono essere determinati solo all’interno di un intervallo di polarizzazione appropriato; possono variare in funzione della ricombinazione, della resistenza in serie e del comportamento del contatto.

Il grafico in scala logaritmica aiuta a individuare l’inizio di una corrente significativa e le variazioni della pendenza. Il grafico lineare facilita l’interpretazione dell’ampiezza e della forma della curva di corrente elevata. Per ottenere una rappresentazione più chiara del funzionamento del diodo, è opportuno esaminare entrambi i grafici insieme ai profili della banda interna e della densità di portatori.Il grafico in scala logaritmica aiuta a individuare l’inizio di una corrente significativa e le variazioni della pendenza. Il grafico lineare facilita l’interpretazione dell’ampiezza e della forma della curva di corrente elevata. Per ottenere una rappresentazione più chiara del funzionamento del diodo, è opportuno esaminare entrambi i grafici insieme ai profili della banda interna e della densità di portatori.

6. Funzionamento interno: come funziona il diodo.6. Funzionamento interno: come funziona il diodo.

La directory «snapshots» contiene i profili spaziali salvati nei punti di polarizzazione selezionati durante la scansione J–V. Le immagini fornite mostrano quattro posizioni lungo tale scansione e illustrano una progressione dall’iniezione bassa verso una polarizzazione diretta più intensa. Le quattro immagini all’interno di ciascun gruppo corrispondono alle stesse posizioni del cursore selezionate, consentendo di confrontare direttamente le diverse quantità fisiche.La directory «snapshots» contiene i profili spaziali salvati nei punti di polarizzazione selezionati durante la scansione J–V. Le immagini fornite mostrano quattro posizioni lungo tale scansione e illustrano una progressione dall’iniezione bassa verso una polarizzazione diretta più intensa. Le quattro immagini all’interno di ciascun gruppo corrispondono alle stesse posizioni del cursore selezionate, consentendo di confrontare direttamente le diverse quantità fisiche.

6.1 Bande di conduzione e di valenza e livelli di quasi-Fermi6.1 Bande di conduzione e di valenza e livelli di quasi-Fermi

Aprire il visualizzatore delle immagini e selezionare Ec.csv, Ev.csv, Fn.csv e Fp.csv. Questi rappresentano rispettivamente il limite della banda di conduzione, il limite della banda di valenza, il livello quasi-Fermi degli elettroni e il livello quasi-Fermi delle lacune. La loro separazione spaziale e le loro pendenze rivelano la barriera di potenziale, la risposta alla polarizzazione applicata e la deviazione dall’equilibrio.Aprire il visualizzatore delle immagini e selezionare Ec.csv, Ev.csv, Fn.csv e Fp.csv. Questi rappresentano rispettivamente il limite della banda di conduzione, il limite della banda di valenza, il livello quasi-Fermi degli elettroni e il livello quasi-Fermi delle lacune. La loro separazione spaziale e le loro pendenze rivelano la barriera di potenziale, la risposta alla polarizzazione applicata e la deviazione dall’equilibrio.

A basse correnti di iniezione, i bordi della banda di conduzione si incurvano attraverso la giunzione a causa del campo elettrostatico intrinseco. Con l’aumentare della polarizzazione diretta, la barriera potenziale effettiva diminuisce e i livelli quasi-Fermi degli elettroni e delle lacune si separano man mano che i portatori di carica attraversano il dispositivo. L’aspetto fondamentale è che l’InP rimane lo stesso semiconduttore su entrambi i lati della giunzione: le variazioni dei profili riflettono le sue caratteristiche elettrostatiche, determinate dalla concentrazione di impurità, e l’iniezione di portatori di carica.A basse correnti di iniezione, i bordi della banda di conduzione si incurvano attraverso la giunzione a causa del campo elettrostatico intrinseco. Con l’aumentare della polarizzazione diretta, la barriera potenziale effettiva diminuisce e i livelli quasi-Fermi degli elettroni e delle lacune si separano man mano che i portatori di carica attraversano il dispositivo. L’aspetto fondamentale è che l’InP rimane lo stesso semiconduttore su entrambi i lati della giunzione: le variazioni dei profili riflettono le sue caratteristiche elettrostatiche, determinate dalla concentrazione di impurità, e l’iniezione di portatori di carica.

Il confronto tra ?? e ?? illustra come la polarizzazione dei terminali influenzi la struttura delle bande interne. In particolare, i livelli quasi-Fermi forniscono una chiara indicazione dello stato di non equilibrio che si stabilisce durante la conduzione in condizioni stabili.Il confronto tra ?? e ?? illustra come la polarizzazione dei terminali influenzi la struttura delle bande interne. In particolare, i livelli quasi-Fermi forniscono una chiara indicazione dello stato di non equilibrio che si stabilisce durante la conduzione in condizioni stabili.

6.2 Densità di elettroni e lacune6.2 Densità di elettroni e lacune

Successivamente, si analizzano i grafici Q_nfree.csv e Q_pfree.csv, che mostrano la distribuzione spaziale degli elettroni liberi e delle lacune. A basse tensioni, gli elettroni sono predominanti sul lato di tipo n, mentre le lacune lo sono sul lato di tipo p. Le loro concentrazioni diminuiscono notevolmente in prossimità della giunzione, riflettendo la zona di svuotamento.Successivamente, si analizzano i grafici Q_nfree.csv e Q_pfree.csv, che mostrano la distribuzione spaziale degli elettroni liberi e delle lacune. A basse tensioni, gli elettroni sono predominanti sul lato di tipo n, mentre le lacune lo sono sul lato di tipo p. Le loro concentrazioni diminuiscono notevolmente in prossimità della giunzione, riflettendo la zona di svuotamento.

All’aumentare della polarizzazione diretta, le cariche minoritarie vengono iniettate nei lati opposti. Ciò aumenta la sovrapposizione delle popolazioni di elettroni e lacune e modifica la forma di entrambi i profili di densità. A livelli di iniezione più elevati, le popolazioni di cariche libere ai due lati della giunzione differiscono notevolmente dai loro valori a bassa polarizzazione.All’aumentare della polarizzazione diretta, le cariche minoritarie vengono iniettate nei lati opposti. Ciò aumenta la sovrapposizione delle popolazioni di elettroni e lacune e modifica la forma di entrambi i profili di densità. A livelli di iniezione più elevati, le popolazioni di cariche libere ai due lati della giunzione differiscono notevolmente dai loro valori a bassa polarizzazione.

Confronta la variazione della densità con la posizione della transizione di drogaggio e, successivamente, con le zone in cui si verifica una forte curvatura della banda di energia. In questo modo si mette in relazione la concentrazione di portatori forniti dal contatto con la barriera elettrostatica prodotta dalla giunzione PN.Confronta la variazione della densità con la posizione della transizione di drogaggio e, successivamente, con le zone in cui si verifica una forte curvatura della banda di energia. In questo modo si mette in relazione la concentrazione di portatori forniti dal contatto con la barriera elettrostatica prodotta dalla giunzione PN.

6.3 Densità di corrente degli elettroni e delle lacune.6.3 Densità di corrente degli elettroni e delle lacune.

Rappresentare graficamente Jn.csv e Jp.csv. Questi grafici mostrano il contributo degli elettroni e delle lacune alla corrente elettrica che attraversa il dispositivo. I singoli componenti possono variare a seconda della posizione, poiché il processo di ricombinazione trasferisce le cariche tra le popolazioni di elettroni e lacune, ma la corrente totale in condizioni stabili dovrebbe rimanere approssimativamente costante quando non si accumula una carica netta.Rappresentare graficamente Jn.csv e Jp.csv. Questi grafici mostrano il contributo degli elettroni e delle lacune alla corrente elettrica che attraversa il dispositivo. I singoli componenti possono variare a seconda della posizione, poiché il processo di ricombinazione trasferisce le cariche tra le popolazioni di elettroni e lacune, ma la corrente totale in condizioni stabili dovrebbe rimanere approssimativamente costante quando non si accumula una carica netta.

A correnti molto basse, l’annullamento reciproco tra i termini di deriva e di diffusione può far sì che i singoli profili di corrente appaiano irregolari o presentino picchi marcati. Le immagini acquisite a bassa polarizzazione mostrano questo comportamento in modo più evidente rispetto alle immagini acquisite a correnti più elevate. È utile esaminare i profili su una scala appropriata e confrontare la corrente totale con il valore J–V al terminale prima di interpretare le caratteristiche numeriche più dettagliate.A correnti molto basse, l’annullamento reciproco tra i termini di deriva e di diffusione può far sì che i singoli profili di corrente appaiano irregolari o presentino picchi marcati. Le immagini acquisite a bassa polarizzazione mostrano questo comportamento in modo più evidente rispetto alle immagini acquisite a correnti più elevate. È utile esaminare i profili su una scala appropriata e confrontare la corrente totale con il valore J–V al terminale prima di interpretare le caratteristiche numeriche più dettagliate.

6.4 Profilo di ricombinazione6.4 Profilo di ricombinazione

La ricombinazione determina la posizione in cui gli elettroni e le lacune iniettate vengono rimossi dalle rispettive popolazioni di portatori di carica. Per confrontare i canali di ricombinazione diretta e quelli assistiti da difetti, aprire la finestra di visualizzazione delle immagini con R_nfree_to_pfree.csv e R_ss_srh.csv. Entrambi i processi dipendono dalla disponibilità di elettroni e lacune, pertanto la loro distribuzione spaziale si modifica in base alle variazioni delle popolazioni di portatori di carica iniettate.La ricombinazione determina la posizione in cui gli elettroni e le lacune iniettate vengono rimossi dalle rispettive popolazioni di portatori di carica. Per confrontare i canali di ricombinazione diretta e quelli assistiti da difetti, aprire la finestra di visualizzazione delle immagini con R_nfree_to_pfree.csv e R_ss_srh.csv. Entrambi i processi dipendono dalla disponibilità di elettroni e lacune, pertanto la loro distribuzione spaziale si modifica in base alle variazioni delle popolazioni di portatori di carica iniettate.

All’inizio della conduzione, il processo di ricombinazione può essere concentrato in prossimità della giunzione, dove le popolazioni di elettroni e lacune si sovrappongono. In caso di iniezione diretta più intensa, la sovrapposizione si estende a una porzione maggiore del semiconduttore e i profili di ricombinazione possono diventare più ampi. Confrontare la posizione e l’ampiezza dei picchi di ricombinazione con le densità di portatori di carica alle stesse condizioni di polarizzazione.All’inizio della conduzione, il processo di ricombinazione può essere concentrato in prossimità della giunzione, dove le popolazioni di elettroni e lacune si sovrappongono. In caso di iniezione diretta più intensa, la sovrapposizione si estende a una porzione maggiore del semiconduttore e i profili di ricombinazione possono diventare più ampi. Confrontare la posizione e l’ampiezza dei picchi di ricombinazione con le densità di portatori di carica alle stesse condizioni di polarizzazione.

Quando si confrontano i due canali, è necessario considerare sia il valore di picco sia l’integrale spaziale. Un picco locale marcato non deve necessariamente prevalere sull’integrale della ricombinazione se un altro canale è attivo su un volume molto più ampio. Questo è particolarmente importante quando si cerca di mettere in relazione i profili di ricombinazione con la corrente elettrica totale.Quando si confrontano i due canali, è necessario considerare sia il valore di picco sia l’integrale spaziale. Un picco locale marcato non deve necessariamente prevalere sull’integrale della ricombinazione se un altro canale è attivo su un volume molto più ampio. Questo è particolarmente importante quando si cerca di mettere in relazione i profili di ricombinazione con la corrente elettrica totale.

6.5 Profili della velocità di generazione netta di elettroni e lacune6.5 Profili della velocità di generazione netta di elettroni e lacune

Le immagini contrassegnate con R_n.csv e R_p.csv forniscono ulteriori informazioni sui valori utilizzati nelle equazioni di continuità per gli elettroni e le lacune. Queste immagini possono evidenziare i punti in cui l’equilibrio tra i termini di generazione, cattura e ricombinazione locali è più marcato. Interpretate i loro segni in base alle convenzioni di rappresentazione grafica utilizzate da OghmaNano.Le immagini contrassegnate con R_n.csv e R_p.csv forniscono ulteriori informazioni sui valori utilizzati nelle equazioni di continuità per gli elettroni e le lacune. Queste immagini possono evidenziare i punti in cui l’equilibrio tra i termini di generazione, cattura e ricombinazione locali è più marcato. Interpretate i loro segni in base alle convenzioni di rappresentazione grafica utilizzate da OghmaNano.

Considerando tutti questi elementi insieme, la struttura a bande, la densità dei portatori di carica, le correnti e i profili di variazione forniscono una spiegazione completa delle caratteristiche della giunzione a semiconduttore: la barriera interna regola l’iniezione, la tensione applicata modifica l’elettrostatica interna e il trasporto e la ricombinazione determinano la corrente di stato stazionario risultante.Considerando tutti questi elementi insieme, la struttura a bande, la densità dei portatori di carica, le correnti e i profili di variazione forniscono una spiegazione completa delle caratteristiche della giunzione a semiconduttore: la barriera interna regola l’iniezione, la tensione applicata modifica l’elettrostatica interna e il trasporto e la ricombinazione determinano la corrente di stato stazionario risultante.

7. Esercizi e interpretazione7. Esercizi e interpretazione

8. Interpretazione del modello e ulteriori analisi8. Interpretazione del modello e ulteriori analisi

Questo esempio utilizza un modello di deriva-diffusione InP a stato stazionario monodimensionale, una struttura semiconduttrice p/n/n⁺ compatta e terminali ohmici idealizzati. È particolarmente adatto a illustrare come la rettificazione elettrica emerga dalla barriera di giunzione e come tale comportamento sia collegato alla distribuzione spaziale delle cariche. Studi più approfonditi sui dispositivi possono includere le caratteristiche di contatto misurate, parametri materiali più precisi, variazioni di temperatura e dettagli geometrici aggiuntivi.Questo esempio utilizza un modello di deriva-diffusione InP a stato stazionario monodimensionale, una struttura semiconduttrice p/n/n⁺ compatta e terminali ohmici idealizzati. È particolarmente adatto a illustrare come la rettificazione elettrica emerga dalla barriera di giunzione e come tale comportamento sia collegato alla distribuzione spaziale delle cariche. Studi più approfonditi sui dispositivi possono includere le caratteristiche di contatto misurate, parametri materiali più precisi, variazioni di temperatura e dettagli geometrici aggiuntivi.

Le proprietà del materiale InP vengono ottenute dal modello Lua selezionato. Per studi quantitativi, verificare i valori calcolati per il gap di banda, l’affinità elettronica, la mobilità e i coefficienti di ricombinazione nell’editor dei materiali o nei dati di output. Controllare la conservazione della corrente e i residui del risolutore, e ripetere alcuni punti di polarizzazione rappresentativi utilizzando una griglia più fine per stabilire la stabilità numerica.Le proprietà del materiale InP vengono ottenute dal modello Lua selezionato. Per studi quantitativi, verificare i valori calcolati per il gap di banda, l’affinità elettronica, la mobilità e i coefficienti di ricombinazione nell’editor dei materiali o nei dati di output. Controllare la conservazione della corrente e i residui del risolutore, e ripetere alcuni punti di polarizzazione rappresentativi utilizzando una griglia più fine per stabilire la stabilità numerica.

Nell’esempio fornito, l’intervallo di variazione standard va da 0 a 3 V. Un esperimento utile da condurre è quello di includere un intervallo di polarizzazione inversa moderato, in modo da poter confrontare direttamente la conduzione in polarizzazione diretta con la corrente di dispersione in polarizzazione inversa. L’intervallo di polarizzazione inversa e le caratteristiche di dispersione rilevanti dal punto di vista fisico dipendono dal dispositivo e dai meccanismi inclusi nella simulazione.Nell’esempio fornito, l’intervallo di variazione standard va da 0 a 3 V. Un esperimento utile da condurre è quello di includere un intervallo di polarizzazione inversa moderato, in modo da poter confrontare direttamente la conduzione in polarizzazione diretta con la corrente di dispersione in polarizzazione inversa. L’intervallo di polarizzazione inversa e le caratteristiche di dispersione rilevanti dal punto di vista fisico dipendono dal dispositivo e dai meccanismi inclusi nella simulazione.

Per ottenere ulteriori informazioni, è possibile variare le concentrazioni di drogaggio una alla volta. In generale, l’aumento della concentrazione nominale di drogaggio restringe la regione di svuotamento e modifica l’equilibrio elettrostatico; la modifica dei coefficienti di ricombinazione influisce sulla corrente, che dipende dall’iniezione. Si tratta di confronti utili perché i materiali e la geometria dello strato possono rimanere invariati.Per ottenere ulteriori informazioni, è possibile variare le concentrazioni di drogaggio una alla volta. In generale, l’aumento della concentrazione nominale di drogaggio restringe la regione di svuotamento e modifica l’equilibrio elettrostatico; la modifica dei coefficienti di ricombinazione influisce sulla corrente, che dipende dall’iniezione. Si tratta di confronti utili perché i materiali e la geometria dello strato possono rimanere invariati.

9. Prossimi passi9. Prossimi passi

Il diodo PN in InP rappresenta un punto di partenza utile per la realizzazione di strutture semiconduttrici più avanzate basate su InP. Un’evoluzione naturale consiste nell’introdurre una regione leggermente drogata e confrontare il comportamento elettrostatico e di trasporto di una giunzione PN e di una giunzione PIN. Un’altra possibile direzione è quella di studiare eterostrutture che coinvolgono InGaAs/InP, dove le differenze nella composizione introducono ulteriori differenze di banda e nuovi comportamenti nell’iniezione di portatori.Il diodo PN in InP rappresenta un punto di partenza utile per la realizzazione di strutture semiconduttrici più avanzate basate su InP. Un’evoluzione naturale consiste nell’introdurre una regione leggermente drogata e confrontare il comportamento elettrostatico e di trasporto di una giunzione PN e di una giunzione PIN. Un’altra possibile direzione è quella di studiare eterostrutture che coinvolgono InGaAs/InP, dove le differenze nella composizione introducono ulteriori differenze di banda e nuovi comportamenti nell’iniezione di portatori.