Statistica di Boltzmann in un contatto GaAs fortemente drogato
Un gruppo di ricerca modella con un codice a statistica di Boltzmann il gas di elettroni nel contatto altamente drogato di un dispositivo in GaAs. A 300 K la densità elettronica nella regione attiva è n = 5·10²⁴ m⁻³ e la massa efficace è m* = 0,067 m_e. Il codice calcola il livello di Fermi rispetto al bordo di banda e, dal confronto con misure capacità-tensione, il gruppo deduce un offset di banda di 0,06 eV. Il revisore osserva che nel fit è stato usato un valore di densità efficace di banda N_c ≈ 4,4·10²³ m⁻³ e chiede di verificare la validità dell'approssimazione classica prima di pubblicare. Il gruppo risponde che il materiale è un semiconduttore e non un metallo, quindi la coda di Boltzmann sarebbe comunque rappresentativa, e che l'eventuale correzione sarebbe «un dettaglio del secondo ordine».
1. Quantificare la degenerazione del gas di elettroni e stimare l'errore commesso usando la statistica classica al posto di quella di Fermi-Dirac.
Read the worked answer
Lunghezza d'onda termica. λ = h/√(2π m* k T). Con m* = 0,067·9,109·10⁻³¹ = 6,10·10⁻³² kg e kT = 4,14·10⁻²¹ J (0,0259 eV): 2π m* kT = 1,59·10⁻⁵¹ → √ = 3,99·10⁻²⁶ kg·m/s → λ = 6,626·10⁻³⁴/3,99·10⁻²⁶ = 1,66·10⁻⁸ m = 16,6 nm; λ³ = 4,6·10⁻²⁴ m³.
Parametro di degenerazione. nλ³ = 5·10²⁴ × 4,6·10⁻²⁴ = 23; con il fattore di spin 2, n/N_c = nλ³/2 = 11,4. Il valore N_c = 2/λ³ = 4,4·10²³ m⁻³ coincide con quello usato dal gruppo: il conto è coerente.
Errore del modello classico. Boltzmann dà E_F − E_c = kT·ln(n/N_c) = 0,0259 × ln(11,4) = 0,0259 × 2,43 = 0,063 eV. Il valore degenere è E_F = (ħ²/2m*)(3π²n)^(2/3) = 9,11·10⁻³⁸ × 2,80·10¹⁷ = 2,55·10⁻²⁰ J = 0,159 eV (la correzione di Sommerfeld a T finita è dello 0,3%). L'errore è quindi ≈ 0,096 eV, quasi 4 kT.
Conseguenza. La condizione classica nλ³ ≪ 1 è violata di oltre un ordine di grandezza: la coda di Boltzmann è inapplicabile. L'offset di 0,06 eV «misurato» è dello stesso ordine dell'errore del modello: è un artefatto, non una proprietà del materiale.
2. Se il codice classico non è utilizzabile, che cosa si deve usare e quali altre approssimazioni vanno verificate prima di attribuire un offset di banda al dato sperimentale?
Read the worked answer
Statistica corretta. Integrale di Fermi-Dirac per il legame n–E_F, oppure l'approssimazione di Joyce-Dixon, che resta accurata fino a n/N_c dell'ordine di 10 (qui siamo al limite superiore). Per una stima rapida, Sommerfeld a T finita: E_F(T) ≈ E_F0[1 − (π²/12)(kT/E_F0)²], correzione dello 0,3%.
Altre approssimazioni da controllare. (i) Banda parabolica: con E_F a 0,16 eV sopra il bordo, la non parabolicità del GaAs corregge la densità degli stati di circa il 10%. (ii) Rinormalizzazione del gap e allargamento delle bande a drogaggio alto, con possibile formazione di banda di impurezza: a 5·10²⁴ m⁻³ si è vicini al limite di Mott, dove il modello a bande rigide inizia a cedere. (iii) Massa efficace dipendente dalla densità. (iv) Temperatura del reticolo: la T del gas non è necessariamente 300 K sotto iniezione o polarizzazione. (v) Il modello della misura C-V (giunzione brusca o graduale) e la sua incertezza, da confrontare con lo scarto che si vuole attribuire all'offset.
Conclusione. Il confronto teoria-misura va rifatto con statistica degenere e con un modello di banda esplicitato. Nella revisione conviene dichiarare esplicitamente che l'offset di 0,06 eV non è confermato, invece di presentarlo come una proprietà del materiale.
3. Esiste una condizione sperimentale o un intervallo in cui il codice classico tornerebbe valido per questa regione attiva? Giustificare quantitativamente.
Read the worked answer
Criterio. Serve nλ³ ≲ 1, cioè n/N_c ≲ 0,5; per un errore entro il 10% sulla densità a parità di livello di Fermi, in pratica n/N_c < 0,1.
Con n fissato a 5·10²⁴ m⁻³. Servirebbe N_c ≥ 5·10²⁵ m⁻³, cioè 115 volte il valore attuale. Poiché N_c ∝ T^(3/2), T = 300 × 115^(2/3) ≈ 300 × 23,6 ≈ 7·10³ K: temperatura alla quale il GaAs non esiste come cristallo.
Con T fissata a 300 K. Servirebbe n ≲ 5·10²² m⁻³, due ordini di grandezza sotto la densità reale del contatto: non è la condizione di funzionamento del dispositivo.
Conclusione operativa. Il modello non si può «salvare» aggiustando le condizioni al contorno: va cambiata la fisica del codice. La risposta al revisore deve contenere il parametro di degenerazione, l'errore stimato in E_F, la versione corretta del calcolo e la nuova interpretazione del confronto con la misura; è utile aggiungere un test di validazione del codice corretto nel limite nλ³ ≪ 1, dove entrambi i modelli devono convergere.