Revisiting Theoretical Modeling of Seebeck Coefficient of Semiconductors

Sintesi Tecnica: Revisione della Modellazione Teorica del Coefficiente di Seebeck nei Semiconduttori

Problematica
Il documento identifica difetti termodinamici fondamentali nelle metodologie esistenti utilizzate per modellare il coefficiente di Seebeck ($\alpha_S$) nei semiconduttori, in particolare quelle basate sull'Equazione di Trasporto di Boltzmann (BTE) e su specifiche definizioni di flusso di calore. L'autore evidenzia tre problemi primari negli approcci convenzionali:

1. Potenziale Confusione: Una frequente confusione tra potenziale elettrochimico (livello di Fermi, $E_F$) e potenziale elettrostatico ($\varphi$) nelle definizioni, che porta a derivazioni errate.
2. Negligenza Arbitraria: Nelle derivazioni BTE, i gradienti spaziali del fondo della banda di conduzione ($E_C$) e del livello di Fermi ($E_F$) vengono spesso arbitrariamente trascurati nel calcolo del gradiente della funzione di distribuzione, nonostante queste variabili cambino attraverso il campione anche in condizioni di circuito aperto (OC).
3. Flusso di Calore Improprio: La definizione di flusso di calore locale ($J_Q = T J_S = J_E - \mu J_N$) si basa su un forte presupposto di processi locali quasi-statici (reversibili) con produzione di entropia nulla. Ciò contraddice i principi fondamentali della termodinamica fuori equilibrio, pur essendo ampiamente utilizzato nelle applicazioni BTE, Landauer-Büttiker e della formula di Kubo.

Metodologia
Per affrontare questi difetti, l'autore propone una metodologia alternativa basata sull'equazione di drift-diffusion che incorpora l'effetto Soret in condizioni di circuito aperto. Questo approccio evita presentazioni esplicite del flusso di calore e l'ambiguità tra potenziali elettrochimici ed elettrostatici.

La derivazione procede in tre fasi:

1. Vicino all'Equilibrio (Termine di Banda): Assumendo l'Equilibrio Locale (LEA) e la neutralità elettrica locale, l'autore deriva un modello a forma chiusa per il termine di banda di $\alpha_S$. Ponendo il flusso di carica netto uguale a zero e utilizzando relazioni di Einstein generalizzate, viene recuperato un "modello reversibile". Questo modello elimina i termini espliciti di velocità e tempo di rilassamento, basandosi invece sul gradiente spaziale del livello di Fermi.
2. Termine di Trascinamento Fononico (Phonon-Drag): La metodologia viene estesa per includere l'effetto di trascinamento fononico aggiungendo un termine di flusso di carica ($J_d$) derivante dallo scambio di quantità di moto con fononi fuori equilibrio. Si dimostra che l'espressione risultante per il termine di trascinamento fononico è identica nella forma al modello convenzionale basato su BTE.
3. Caso Ballistico: Per il regime ballistico (dove la LEA non è valida), l'autore utilizza la formulazione di Landauer per calcolare i flussi di carica dalle estremità calda e fredda, derivando l'espressione convenzionale del coefficiente di Seebeck ballistico.

Il documento unifica questi modelli (BTE, Reversibile e Ballistico) in un'unica forma integrale dipendente da un esponente energetico $r$, che varia in base ai meccanismi di scattering e alla dimensionalità del materiale.

Contributi Chiave

• Correzione Termodinamica: Il documento critica rigorosamente le incoerenze termodinamiche nelle derivazioni BTE standard e nelle definizioni del flusso di calore.
• Derivazione del Modello Reversibile: Dimostra che il modello reversibile (presentato nel lavoro precedente dell'autore) è valido per il termine di banda del coefficiente di Seebeck nei casi di quasi-equilibrio, nonostante l'omissione dei termini di velocità e rilassamento.
• Limiti Teorici: Lo studio stabilisce che il modello reversibile fornisce un limite inferiore teorico per il termine di banda del coefficiente di Seebeck (corrispondente a $r=1/2$ per materiali 3D), mentre il modello BTE con scattering di impurità ionizzate debolmente schermate fornisce un limite superiore ($r=3$).
• Quadro Unificato: Presenta un quadro matematico unificato per confrontare i modelli reversibile, ballistico e BTE basandosi sulla dipendenza energetica della densità degli stati, della velocità e del tempo di rilassamento.

Risultati
È stato eseguito un caso di studio sul Silicio altamente drogato con tipo n ($n_d > 1 \times 10^{19} \text{ cm}^{-3}$), confrontando i tre modelli con i dati sperimentali di Geballe et al. e Sadhu et al.:

• Modello Reversibile: Ha fornito previsioni "abbastanza buone" con basso costo computazionale. Ha agito costantemente come limite inferiore per i dati sperimentali, con scostamenti relativi medi del 13%, 18% e 5% per concentrazioni di drogaggio di $2,7 \times 10^{19}$, $6,0 \times 10^{19}$ e $7,0 \times 10^{19} \text{ cm}^{-3}$, rispettivamente.
• Modello Ballistico: Ha mostrato prestazioni miste, con scostamenti del 3%, 6% e 16%. Ha sovrastimato i dati per il caso del nanowire ($7,0 \times 10^{19} \text{ cm}^{-3}$).
• Modello BTE: Utilizzando l'assunzione di scattering di impurità ionizzate debolmente schermate (il limite superiore), il modello ha sovrastimato significativamente i valori sperimentali, con scostamenti relativi medi che hanno raggiunto il 60% - 90%. Ciò evidenzia la sensibilità dei risultati BTE alla scelta della formulazione del tempo di rilassamento.

Significatività e Rivendicazioni
Il documento sostiene che il modello reversibile proposto offre un'alternativa concisa, teoricamente solida ed efficiente dal punto di vista computazionale per valutare le proprietà termoelettriche, particolarmente nei casi di quasi-equilibrio e reversibili. Eliminando la necessità di complessi calcoli del tempo di scattering e evitando i paradossi termodinamici, il modello viene presentato come uno strumento valido per calcoli ad alto rendimento (high-throughput) e screening dei materiali. L'autore sottolinea che, sebbene il modello trascuri i meccanismi irreversibili (scattering), esso cattura accuratamente il limite inferiore del coefficiente di Seebeck, rendendolo una linea di base robusta per la valutazione delle proprietà termoelettriche.