Revisiting Theoretical Modeling of Seebeck Coefficient of Semiconductors

Technische Samenvatting: Een herziening van de theoretische modellering van de Seebeck-coëfficiënt van halfgeleiders

Probleemstelling
Het artikel identificeert fundamentele thermodynamische gebreken in bestaande methodologieën die worden gebruikt om de Seebeck-coëfficiënt ($\alpha_S$) van halfgeleiders te modelleren, met name die welke steunen op de Boltzmann-transportvergelijking (BTE) en specifieke definities van warmteflux. De auteur wijst op drie primaire problemen in conventionele benaderingen:

1. Potentiële verwarring: Een frequente samensmelting van elektrochemische potentiaal (Fermi-niveau, $E_F$) en elektrostatische potentiaal ($\varphi$) in definities, wat leidt tot incorrecte afleidingen.
2. Willekeurige verwaarlozing: In BTE-afleidingen worden de ruimtelijke gradiënten van de bodem van de geleidingsband ($E_C$) en het Fermi-niveau ($E_F$) vaak willekeurig verwaarloosd bij het berekenen van de gradiënt van de distributiefunctie, ondanks het feit dat deze variabelen over een monster heen veranderen, zelfs onder open-circuit (OC) condities.
3. Onjuiste warmteflux: De definitie van lokale warmteflux ($J_Q = T J_S = J_E - \mu J_N$) berust op een sterke aanname van lokale quasi-statische (reversibele) processen met nul entropieproductie. Dit spreect fundamentele principes van niet-evenwichtsthermodynamica tegen, maar wordt breed toegepast in BTE-, Landauer-Büttiker- en Kubo-formule-toepassingen.

Methodologie
Om deze gebreken aan te pakken, stelt de auteur een alternatieve methodologie voor op basis van de drift-diffusievergelijking inclusief het Soret-effect onder open-circuit condities. Deze benadering vermijdt expliciete presentaties van warmteflux en de ambiguïteit tussen elektrochemische en elektrostatische potentialen.

De afleiding verloopt in drie stadia:

1. Nabij evenwicht (Bandterm): Door uit te gaan van lokaal evenwicht (LEA) en lokale elektrische neutraliteit, leidt de auteur een gesloten model af voor de bandterm van $\alpha_S$. Door de netto ladingflux gelijk te stellen aan nul en gebruik te maken van gegeneraliseerde Einstein-relaties, wordt een "reversibel model" teruggewonnen. Dit model elimineert expliciete elektronensnelheid en relaxatietijd-termen, en vertrouwt in plaats daarvan op de ruimtelijke gradiënt van het Fermi-niveau.
2. Fonon-drag term: De methodologie wordt uitgebreid om het fonon-drag effect te bevatten door een ladingflux-term ($J_d$) toe te voegen die voortkomt uit momentumuitwisseling met niet-evenwichtsfononen. De resulterende expressie voor de fonon-drag term wordt getoond identiek te zijn in vorm aan het conventionele BTE-gebaseerde model.
3. Ballistisch geval: Voor het ballistische regime (waar LEA ongeldig is), maakt de auteur gebruik van de Landauer-formulering om ladingfluxen van koude en warme terminals te berekenen, waarbij de conventionele ballistische Seebeck-coëfficiënt expressie wordt afgeleid.

De auteur verenigt deze modellen (BTE, Reversibel, en Ballistisch) in een enkele integrale vorm die afhankelijk is van een energie-exponent $r$, die varieert op basis van verstrooiingsmechanismen en de dimensionaliteit van het materiaal.

Belangrijkste Bijdragen

• Thermodynamische Correctie: Het artikel bekritiseert rigoureus de thermodynamische inconsistenties in standaard BTE-afleidingen en warmteflux-definities.
• Afleiding van het Reversibele Model: Het toont aan dat het reversibele model (voorheen gepresenteerd in het eerdere werk van de auteur) geldig is voor de bandterm van de Seebeck-coëfficiënt in nabij-evenwicht gevallen, ondanks het weglaten van snelheid en relaxatie-termen.
• Theoretische Grenzen: De studie stelt vast dat het reversibele model een theoretische ondergrens vormt voor de bandterm van de Seebeck-coëfficiënt (overeenkomend met $r=1/2$ voor 3D-materialen), terwijl het BTE-model met zwak gescreende geïoniseerde onzuiverheidsverstrooiing een bovengrens vormt ($r=3$).
• Verenigd Kader: Het presenteert een verenigd wiskundig kader om de reversibele, ballistische en BTE-modellen te vergelijken op basis van de energieafhankelijkheid van de toestandsdichtheid, snelheid en relaxatietijd.

Resultaten
Een casestudy werd uitgevoerd op sterk n-gedoteerd Silicium ($n_d > 1 \times 10^{19} \text{ cm}^{-3}$), waarbij de drie modellen werden vergeleken met experimentele gegevens van Geballe et al. en Sadhu et al.:

• Reversibel Model: Leverde "redelijk goede" voorspellingen met lage computationele kosten. Het diende consistent als de ondergrens voor de experimentele gegevens, met gemiddelde relatieve afwijkingen van 13%, 18% en 5% voor dopingconcentraties van respectievelijk $2,7 \times 10^{19}$, $6,0 \times 10^{19}$ en $7,0 \times 10^{19} \text{ cm}^{-3}$.
• Ballistisch Model: Toonde een gemengde prestatie, met afwijkingen van 3%, 6% en 16%. Het overschatte de gegevens voor het nanowire-geval ($7,0 \times 10^{19} \text{ cm}^{-3}$).
• BTE-model: Bij gebruik van zwak gescreende geïoniseerde onzuiverheidsverstrooiing (de bovengrens-aanname), overschatte het model de experimentele waarden aanzienlijk, met gemiddelde relatieve afwijkingen oplopend tot 60% tot 90%. Dit benadrukt de gevoeligheid van BTE-resultaten voor de keuze van de relaxatietijd-formulering.

Betekenis en Claims
Het artikel claimt dat het voorgestelde reversibele model een beknopte, theoretisch solide en computationeel efficiënte alternatief biedt voor het evalueren van thermoelektrische eigenschappen, met name in nabij-evenwicht en reversibele gevallen. Door de noodzaak van complexe verstrooiingstijd-berekeningen weg te nemen en thermodynamische paradoxen te vermijden, wordt het model gepresenteerd als een levensvatbaar instrument voor high-throughput berekeningen en materiaalscreening. De auteur benadrukt dat hoewel het model irreversibele mechanismen (verstrooiing) verwaarloost, het de ondergrens van de Seebeck-coëfficiënt accuraat vastlegt, waardoor het een robuuste baseline vormt voor de evaluatie van thermoelektrische eigenschappen.