A stochastic simulation of the dislocation-mediated etching of porous GaN distributed Bragg reflectors

Titolo

Una simulazione stocastica dell'incisione mediata da dislocazioni di riflettori di Bragg distribuiti (DBR) in GaN poroso.

1. Il Problema

L'incisione elettrochimica selettiva del nitruro di gallio (GaN) permette di creare strutture nanoscopiche porose, fondamentali per applicazioni come i riflettori di Bragg distribuiti (DBR) che migliorano l'estrazione della luce nei LED e nei laser.
Il processo tradizionale richiede strati alternati di GaN fortemente drogato (n-doped) e non intenzionalmente drogato (NID). Poiché l'incisione avviene solo nelle regioni drogate, è necessario un meccanismo per far arrivare l'agente chimico agli strati droganti sepolti sotto lo strato superficiale NID.

• Approccio precedente: Utilizzo di trincee profonde create litograficamente per esporre gli strati laterali. Questo metodo è costoso e complesso.
• Approccio alternativo (senza litografia): Sfruttamento delle dislocazioni a vite (Threading Dislocations - TD) presenti naturalmente nel GaN eteroepitassiale. Queste agiscono come "tubi" (nanopipes) che permettono all'agente chimico di attraversare gli strati NID.
• La sfida scientifica: Il modello iniziale ("modello kebab") ipotizzava che ogni TD formasse un canale verticale continuo attraverso tutti gli strati. Tuttavia, dati sperimentali recenti (tomografia FIB) hanno rivelato una morfologia più complessa chiamata "a cascata": i percorsi di incisione possono deviare lateralmente attraverso uno strato drogato, spegnere una TD e riattivare un'altra TD adiacente. La fisica esatta che governa la probabilità di questi eventi (quando una TD si attiva, si disattiva o si propaga lateralmente) non è ancora pienamente compresa, rendendo difficile prevedere la morfologia finale e la cinetica di incisione.

2. Metodologia

Gli autori hanno sviluppato una simulazione stocastica in Python (ottimizzata con la libreria Numba) per modellare il processo di incisione.

• Modello Computazionale:

  • La struttura DBR è rappresentata come un array di pixel (5 nm x 5 nm).
  • Gli elementi includono: strati NID (probabilità di incisione = 0), strati drogati (probabilità $P_{doped}$), dislocazioni (TD) e agente chimico.
  • Le TD sono modellate come linee verticali con una spaziatura di 500 nm.
  • Algoritmo: L'agente chimico si diffonde cercando i pixel adiacenti. L'incisione di un pixel è determinata da una probabilità assegnata.
  • Criterio di non fusione: È stato aggiunto un vincolo per impedire che pori adiacenti si fondano, mantenendo la separazione dei domini di porosità associati a ciascun nucleo di TD, come osservato sperimentalmente.
  • Metriche: Ogni pixel inciso genera un'unità di carica simulata, permettendo di generare curve di chronoamperometria (corrente vs tempo).
  • Classificazione morfologica: I risultati sono classificati in "Kebab" (KX, dove una TD attraversa X strati consecutivi) e "Cascata" (C, dove una TD si disattiva e si riattiva più in basso dopo un percorso laterale).

• Esperimenti di Confronto:

  • Sono stati cresciuti campioni DBR su wafer di Si tramite MOVPE.
  • Campioni:
    • Campion A: Strati drogati e NID di spessore uniforme.
    • Campioni B e C: Spessori degli strati drogati variabili (crescenti o decrescenti verso la superficie).
  • Incisione: I campioni sono stati incisi elettrochimicamente in acido ossalico a diverse tensioni (6V, 7V, 8V, 12V) senza litografia.
  • Caratterizzazione: Misurazioni di chronoamperometria e imaging SEM/FIB per analizzare la morfologia dei pori.

3. Contributi Chiave

1. Sviluppo di un modello stocastico: Creazione di un codice di simulazione capace di riprodurre sia la cinetica di corrente (chronoamperometria) che la morfologia complessa dei pori (cascate vs kebab) senza parametri fisici microscopici complessi, ma basandosi su probabilità di transizione.
2. Correlazione Tensione-Probabilità: Dimostrazione che variare la tensione di incisione sperimentale può essere mappato nella simulazione variando la probabilità di incisione degli strati drogati ($P_{doped}$).
3. Spiegazione della morfologia "a cascata": Il modello conferma che la morfologia a cascata emerge naturalmente quando la probabilità di incisione laterale è significativa rispetto alla propagazione verticale, spiegando perché alcune TD si "spengono" e altre si "riattivano".
4. Generalizzabilità: Applicazione della simulazione a strutture non uniformi (spessori variabili) e a geometrie alternative (incisione da trincee laterali), mostrando che il meccanismo mediato da TD è rilevante anche in contesti diversi dal DBR standard.

4. Risultati

• Corrispondenza Chronoamperometrica:
  • La simulazione riproduce le oscillazioni di corrente osservate sperimentalmente, che corrispondono all'incisione strato per strato.
  • A tensioni più alte (e quindi $P_{doped}$ più elevate nella simulazione), le oscillazioni sono più marcate e l'incisione è più "strato per strato".
  • A tensioni più basse, le oscillazioni si smorzano rapidamente; la simulazione mostra che ciò è dovuto a una propagazione non uniforme del fronte di incisione, dove alcune TD penetrano più velocemente di altre, mescolando l'incisione di strati diversi.
• Morfologia dei Pori:
  • La simulazione riproduce la transizione da strutture "Kebab" (dominanti ad alte tensioni/$P_{doped}$) a strutture "Cascata" (dominanti a basse tensioni/$P_{doped}$).
  • I dati sperimentali mostrano che a 12V i pori sono più grandi e le TD tendono a formare percorsi più continui (K5), mentre a 6V prevale la morfologia a cascata. La simulazione replica questa tendenza variando $P_{doped}$.
• Limiti del Modello:
  • La simulazione non riesce a riprodurre la formazione di pori conici negli strati NID vicino alle TD (osservati sperimentalmente ad alta tensione), poiché il modello assume probabilità zero di incisione per gli strati NID (tranne che nel percorso della TD).
  • La simulazione tende a sottostimare l'incisione incompleta degli strati più profondi osservata sperimentalmente a basse tensioni.
• Applicazione a Geometrie Alternative:
  • Simulando l'incisione da una trincea laterale (senza contatto dalla superficie superiore), il modello predice che, man mano che ci si allontana dalla trincea, l'attivazione delle TD diventa più probabile, portando a una perdita di allineamento dei pori e alla scomparsa della birifrangenza osservata vicino alla trincea.

5. Significato e Implicazioni

Questo lavoro fornisce un ponte cruciale tra i dati sperimentali macroscopici (corrente) e la microstruttura nanoscopica complessa dei DBR in GaN poroso.

• Comprensione del Meccanismo: Suggerisce che l'effetto della tensione applicata non è uniforme su tutto il processo: aumenta la probabilità di incisione degli strati drogati più di quanto non faccia quella delle TD negli strati NID, indicando meccanismi fisici distinti (es. tunneling Zener).
• Strumento di Progettazione: La simulazione offre uno strumento predittivo per ottimizzare i parametri di incisione (tensione, tempo) per ottenere specifiche morfologie di pori e curve di corrente, utile per la fabbricazione di dispositivi fotonici ad alte prestazioni.
• Versatilità: La capacità del modello di adattarsi a strutture con spessori variabili e diverse geometrie di accesso all'agente chimico suggerisce che l'incisione mediata da dislocazioni è un fenomeno ubiquitario che influenza le proprietà dei materiali anche in configurazioni non convenzionali, aprendo nuove strade per la progettazione di dispositivi optoelettronici su larga scala.