fixed width

8.1 Model driftu a difuze

Podrobnější popis

1        Základní rovnice pro polovodiče

Popis elementární podstaty fyzikálních jevů v strukturách z polovodičových materiálů je umožněn matematickými nástroji kvantové fyziky. Jejich použitím získáme model chování částic na atomární úrovni ve formě vlnových funkcí, jejichž kmitočet a vlnová délka je provázána s energií a momentem hybnosti částic. Chování elektronu v krystalu polovodiče popisujeme jako superpozici různých vlnových funkcí, tzv. vlnový paket (balík). Energie a hybnosti tohoto paketu získáme řešením Schrödingerovy rovnice, která je srdcem teorie kvantové fyziky. Zkoumání pohybu elektronu v struktuře polovodiče potom představuje problém jednoho vlnového paketu, který je vystaven vlivu vnějších polí a vlivu vzájemné interakce (kolizí) s krystalovou mřížkou.

Avšak modelování polovodičových prvků a součástek v praxi představuje problém transportu částic v makroskopickém měřítku. V takovém případě uvažujeme popis pohybu nosičů náboje (např. elektronu) mezi jednotlivými kolizemi Newtonovými zákony a samotné interakce popisujeme statisticky nástroji kvantové fyziky. Takový systém je možné popsat Boltzmanovou kinetickou rovnicí, jejímž řešením je distribuční funkce nosičů náboje v objemu polovodiče. Takové řešení však vyžaduje náročný matematický aparát, proto byla zavedena řada zjednodušení, které redukují takový systém v jednoduchou soustavu parciálních diferenciálních rovnic, které tvoří matematický popis polovodičů nazývaný model driftu a difůze. Tyto rovnice jsou také v některých publikacích nazývány jako Základní rovnice pro polovodiče. Jejich forma vychází z Maxwellových rovnic (8.1B.1), ze kterých se dají pomocí idealizujících předpokladů snadno odvodit, jak se přesvědčíme v následujících odstavcích.

×H=J+Dt, ( 8.1B.1a )
×E=Bt ( 8.1B.1b )
D=ρ, ( 8.1B.1c )
B=0. ( 8.1B.1d )

V uvedených rovnicích E značí elektrického pole a D elektrickou indukci, H intenzitu magnetického pole a B magnetickou indukci, J značí proudovou hustotu, ρ hustotu volného náboje a t čas.

Pro úplný popis je potřeba definovat tzv. materiálové vztahy

D=εE, ( 8.1B.2a )
B=μH, ( 8.1B.2b )

kde ε představuje permitivitu a μ permeabilitu prostředí. Obecně se jedná o tenzory druhého řádu (anizotropní charakter), avšak pro běžné materiály používané polovodičových prvcích mohou být považovány za skaláry (isotropní charakter). Hodnoty relativní permitivity pro nejběžnější materiály jsou uvedeny v tab. 8.1B.1. Vztahy (8.1B.2) jsou platné pouze v polovodičových materiálech s časově invariantní permitivitou a permeabilitou a při zanedbání polarizace vlivem mechanických sil. Z tohoto důvodu nejsou vztahy (8.1B.2) dostatečné pro popis piezoelektrických a feroelektrických jevů.

Tab. 8.1B.1Typické hodnoty relativní permitivity [44].
Materialεr []
Si11,7
SiO23,9
Si3N47,2
GaAs12,5
Ge16,1

1.1     Poissonova rovnice

Třetí Maxwellovu rovnici (8.1B.1c), která má formu Poissonovy parciální diferenciální rovnice, můžeme použitím vztahů (8.1B.2a) a (8.1B.3) rozepsat do tvaru pro elektrický potenciál Ψ (8.1B.4).

E=Ψ, ( 8.1B.3 )
(εΨ)=ρ ( 8.1B.4 )

Pro použití rovnice (8.1B.4) v analýze polovodičů definujeme hustotu volného náboje ρ jako součin elementárního náboje q a koncentrace nosičů náboje, tj, součtu koncentrace kladně nabitých děr p, koncentrace záporně nabitých elektronů n, které jsou kompenzovány koncentrací pevného náboje příměsí C (donorů či akceptorů) v atomové mřížce polovodiče.

ρ=q(pn+C), ( 8.1B.5 )

kde C = ND - NA odpovídá koncentraci příměsí definované jako součet koncentrace ionizovaných donorů ND a ionizovaných akceptorů NA včetně znaménka kompenzovaného náboje.

Nyní můžeme sloučením (8.1B.4) a (8.1B.5) obdržet finální tvar Poissonovy rovnice (8.1B.6) vhodnou pro modelování elektrostatického pole polovodičových struktur.

(Ψ)=qε(npC) ( 8.1B.6 )

1.2     Rovnice kontinuity

Rovnici kontinuity můžeme odvodit přímo z první Maxwellovy rovnice (8.1B.1a). Aplikováním operátoru divergence na tuto rovnici a s uvážením toho, že divergence rotace vektoru je rovna nule, obdržíme

×H=J+ρt=0. ( 8.1B.7 )

Pokud vektor proudové hustoty J rozložíme na součet složky Jp způsobenou děrami a složky Jn způsobenou elektrony

J=Jp+Jn ( 8.1B.8 )

a uvážením toho, že pevný náboj v krystalické mřížce polovodiče je časově invariantní, tj.

Ct=0, ( 8.1B.9 )

dostáváme

(Jp+Jn)+qt(pn)=0. ( 8.1B.10 )

Tato rovnice vyjadřuje to, že zřídlové oblasti proudové hustoty jsou plně kompenzovány časovou změnou volného náboje. Abychom obdrželi rovnici kontinuity zvlášť pro elektrony a díry, musíme zavést společnou funkci R, která vyjadřuje četnost generace a rekombinace elektronů a děr. Potom dostáváme

(Jn)qnt=qR, ( 8.1B.11a )
(Jp)+qpt=qR. ( 8.1B.11b )

Vyjádření četnosti R vyžaduje dobrou znalost fyzikálních mechanismů generace a rekombinace v polovodičových materiálech a je jedním z klíčových parametrů pro získání relevantních výsledku simulací polovodičových struktur.

1.3     Kinetická rovnice nosičů náboje

V předchozí kapitole jsme uvedli vektor proudové hustoty rozdělený na dvě složky (8.1B.8). Proudovou hustotu nabitých částic můžeme obecně rozepsat jako součin náboje jedné částice, koncentrace částic a průměrné rychlosti (driftové rychlosti) částic. Potom pro proudovou hustotu elektronů a proudovou hustotu děr můžeme psát

Jn=qnvn, ( 8.1B.12a )
Jp=qpvp, ( 8.1B.12b )

kde vn představuje vektor driftové rychlosti elektronů a vp vektor driftové rychlosti děr. Hlavním problémem je určení vztahu driftové rychlosti nosičů vůči silám působícím na nosiče. Pomoci řady zjednodušení a předpokladů [44] je možné příčinu nenulové driftové rychlosti rozložit na působení Lorenzovy síly a difúzní síly. Potom můžeme psát

Jn=qnμnE+qDnn, ( 8.1B.13a )
Jp=qpμpEqDpp. ( 8.1B.13b )

U obou rovnic, první člen na pravé straně představuje složku proudové hustoty způsobenou Lorenzovou silou s uvážením vlivu pouze elektrického pole (vliv magnetické indukce je zanedbána) a zavádí efektivní pohyblivost elektronů μn a děr μp, které jsou funkcí průměrné četnosti kolize nosičů, a jejich hodnota je zpravidla definována empirickými vztahy. Pohyblivost nosičů náboje je veličina, která má velký vliv na vlastnosti polovodičových součástek a proto je její správná formulace velmi důležitá. Druhý člen na pravé straně rovnic (8.1B.13) představuje vliv difůze nosičů ve směru gradientu jejich koncentrace a zavádí Einsteinovi difúzní konstanty

Dn=μnkTq, ( 8.1B.14a )
Dp=μpkTq, ( 8.1B.14b )

kde k je Boltzmanova konstanta a T je teplota.

2        Fyzikální parametry

Uvedené rovnice (8.1B.6), (8.1B.11) a (8.1B.13) tvoří základní rovnice umožňující analýzu a simulaci většiny polovodičových prvků. Avšak mohou nastat případy, kdy nabývají na významu fyzikální jevy jež nejsou dostatečně v těmito rovnicemi podchyceny. Tak je tomu v případě tunelových proudů skrz velmi tenké vrstvy, které mají významný vliv na vlastnosti prvku. V praxi se ukázalo efektivnější aplikovat v těchto případech empirické korekce, než daný systém rovnic dále rozšiřovat.

Fyzikální parametry v těchto rovnicích definují geometrii prvku, druh polovodičového materiálu a technologii použitou při jeho výrobě. Jejich modelováním se budeme zabývat v následujících odstavcích.

2.1     Dotační profil

Distribuce příměsí v objemu polovodiče, mající za následek nevlastní vodivost typu N nebo P, určuje geometrii a funkci polovodičového prvku. Proto je to naprosto zásadní vstupní informace a tudíž je pochopitelné, že pro přesné výsledky simulací je zapotřebí věrné modelování používaných výrobních procesů, jako jsou iontová implantace, difůze, tepelná oxidace, epitaxní narůstání a jiné, které tuto distribuci (dotační profil) vytvářejí. Díky šíři tohoto tématu se omezíme jen na procesy iontové implantace a difůze atomů příměsi do objemu polovodičového substrátu. Výsledné profily těchto procesů jsou nejčastěji modelovány Gaussovou funkcí

G(x)=ae(xb)22c2, ( 8.1B.15 )

kde x představuje prostorovou souřadnici, a, b a c jsou reálné konstanty a e je eulerovo číslo. Gaussova funkce je sudá (souměrná podle osy y) a má tvar „zvonu“, kterému přisuzujeme charakter distribuce příměsí v objemu polovodiče jako důsledek lokálně aplikovaných uvedených technologických procesů. Potom konstanta a určuje špičkovou koncentraci příměsí (výšku „zvonu“), b pozici aplikace technologického procesu (souřadnici vrcholu „zvonu“) a c je úměrná difúzi příměsí do substrátu (šířce „zvonu“). V (8.1B.5) jsme zavedli koncentraci pevného náboje C, která odpovídá distribuci příměsí v polovodiči. Pokud uvažujeme, že ND a NA značí špičkovou koncentraci akceptorů a donorů (odpovídají koeficientu a), potom můžeme pro obecný P-N přechod v případě jedno rozměrného (1D) problému psát

C=NDi+NDG(xD,bD,cD)NAG(xA,bA,cA), ( 8.1B.16 )

kde NDi představuje koncentraci příměsí výchozího substrátu, který bývá zpravidla lehce dotován donory, tj. je typu N.

Uvažujme P-N diodu redukovanou na 1D problém o celkové délce L = 8 μm vyrobené ze substrátu křemíku s nízkou výchozí koncentrací donorů NDi = 1021 m-3 a vysoce dopovanými oblastmi na protilehlých koncích substrátu tvořené koncentracemi donorů ND = 1023 m-3 a akceptorů NA = 1023 m-3. Šířku P-N přechodu požadujeme LJ = 1 μm. Výsledný dotační profil, zobrazený na obr. 8.1B.1, potom vyjádříme následujícím předpisem

C=NDi+NDe((x0)22c2)NAe((xL)22c2), ( 8.1B.17 )

kde c vyjádříme pomocí žádané šířky přechodu jako

c=LJ2ln(NA/NDi). ( 8.1B.18 )
Obr. 8.1B.1
Obr. 8.1B.1Příklad dotačního profilu P-N diody

2.2     Pohyblivost nosičů

Pohyblivost μ je důležitá kvantita, která dává do souvislosti driftovou rychlost nosičů v a elektrické pole E, kterému jsou tyto nosiče vystaveny, vztahem

v=μE. ( 8.1B.19 )

Pohyb nosičů v polovodiči je narušen kolizemi s krystalovou mřížkou, příměsemi a defekty, tudíž dráha pohybu nosiče mezi dvěma body není přímka. Z důvodu snadnější představitelnosti tohoto jevu se zavedla právě veličina pohyblivosti, která se modeluje s ohledem na mechanismy, které ovlivňují kolize nosičů. Tyto mechanizmy mohou být způsobeny např.

  • atomovou mřížkou materiálu
  • ionizovanými nečistotami
  • kolizí nosič-nosič
  • neutrálními nečistotami
  • intenzitou elektrického pole

Často si vystačíme pouze s tabulkovou hodnotou μn,p300, udávanou pro daný materiál při teplotě 300K, avšak při tvorbě věrného modelu polovodičové struktury je formulace pohyblivosti nosičů klíčová, protože podchycuje fyzikální jevy, které mají zásadní vliv na vlastnosti polovodičových prvků. V následujících odstavcích se soustředíme na nejběžnější modely pohyblivosti.

Závislost pohyblivosti nosičů na teplotě mřížky

Jedním ze základních mechanismu ovlivňující kolize jsou vibrace atomů v mřížce čistého krystalu polovodiče. Intenzita těchto vibrací je úměrná teplotě mřížky. Modelování tohoto jevu empirickými vztahy (8.1B.20) se opírá o experimentálně naměřené hodnoty pohyblivosti nosičů při různých teplotách, ze kterých byly určeny hodnoty koeficientů αn a αp.

μnL=μn300(T300K)αn, ( 8.1B.20a )
μpL=μp300(T300K)αp. ( 8.1B.20b )

Vynesené závislosti pohyblivosti nosičů na teplotě v krystalu Si a GaAs jsou uvedeny na obr. 8.1B.2. Typické hodnoty pohyblivostí a koeficientů uvádí tab. 8.1B.2.

Obr. 8.1B.2aObr. 8.1B.2b
a) krystal Sib) krystal GaAs
Obr. 8.1B.2Závislost pohyblivosti elektronů a děr na teplotě krystalu polovodiče
Tab. 8.1B.2Hodnoty pohyblivostí a koeficientů jejich tepelné závislosti [44].
SiGaAs
μn300 [m2V-1s-1]0,140,8
μp300 [ m2V-1s-1]0,40,4
αn []2,21
αp []2,22,1
Závislost pohyblivosti nosičů na koncentraci ionizovaných příměsí

Ionizované příměsi představují především příměsi donorů a akceptorů. Vztahy popisující vliv jejich koncentrace na pohyblivost nosičů jsou definovány pomocí tzv. celkové koncentrace ionizovaných příměsí CI. Tato je definována jako součet koncentrací pevného náboje všech m druhů ionizovaných příměsí násobený absolutní hodnotou jejich normovaného náboje Z

CI=i=1m|Zi|Ci. ( 8.1B.21 )

Pro obvyklé příměsi je |Z|=1. Avšak, v některých aplikacích, jako jsou solární články s příměsí zinku, mohou být akceptory dvojitě ionizované, což odpovídá hodnotě |Z|=2. Empirické vztahy (8.1B.22) popisují vliv koncentrace příměsí na pohyblivost nosičů v objemu Si a zároveň postihují vliv teploty mřížky krystalu. Na obr. 8.1B.3 jsou vykresleny jejich průběhy při teplotách 300K a 400K.

μnLI=0,0088(T300K)0,57+0,1252(T300K)2,331+CI1,4321023(T300K)2,546 ( 8.1B.22a )
μpLI=0,00543(T300K)0,57+0,0407(T300K)2,331+CI2,671023(T300K)2,546 ( 8.1B.22b )
Obr. 8.1B.3aObr. 8.1B.3b
a) pohyblivost elektronůb) pohyblivost děr
Obr. 8.1B.3 Pohyblivost nosičů náboje v Si v závislosti na koncentraci příměsí při teplotách 300K a 400K
Závislost pohyblivosti nosičů na intenzitě elektrického pole

Dalším významný jev jehož důsledkem je změna pohyblivostí nosičů je saturace jejich driftové rychlosti při vysokých intenzitách elektrického pole. Pro křemík je možné použít následující empirický vztah

μn,pLE=μn,pL(1+(En,pEn,pcrit)βn,p)1/βn,p, ( 8.1B.23 )

kde En,p značí velikost složky elektrického pole ve směru vektoru proudové hustoty. Mezní hodnota elektrického pole En,pcrit a koeficient βn,p jsou volené tak, aby výsledný průběh pohyblivosti odpovídal experimentálně zjištěným hodnotám. Tyto koeficienty a jejich závislost na teplotě vyjadřují následující vztahy.

Encrit=857103(T300K)1,55, ( 8.1B.24a )
Epcrit=1800103(T300K)1,68, ( 8.1B.24b )
βn=1,11(T300K)0,66, ( 8.1B.25a )
βp=1,21(T300K)0,17. ( 8.1B.25b )

Na obr. 8.1B.4 jsou znázorněny průběhy odpovídající výše uvedeným vztahům. Obr. 8.1B.5 znázorňuje efekt saturace driftové rychlosti nosičů vlivem velikosti intenzity elektrického pole, kdy pro driftovou rychlost platí

vn,p=μn,pLEEn,p. ( 8.1B.26 )

Obr. 8.1B.4aObr. 8.1B.4b
a) pohyblivost elektronůb) pohyblivost děr
Obr. 8.1B.4Pohyblivost nosičů náboje v Si v závislosti na velikosti intenzity elektrického pole při teplotách 300K a 400K
Obr. 8.1B.5aObr. 8.1B.5b
a) driftová rychlost elektronůb) driftová rychlost děr
Obr. 8.1B.5Driftová rychlost nosičů náboje v Si v závislosti na velikosti intenzity elektrického pole při teplotách 300K a 400K

Pro elektrony v GaAs je možno tento jev modelovat následujícím vztahem

μnLE=μnL+vsat(En)3(Encrit)41+(EnEncrit)4, ( 8.1B.27 )

kde En opět značí velikost složky elektrického pole ve směru vektoru proudové hustoty Encrit je jeho kritická hodnota a vsat je tepelně závislá saturační rychlost elektronů, která je modelovaná vztahem

vsat=vsat300(1An)+An(T300K), ( 8.1B.28 )

kde hodnota koeficientu An je určována experimentálně. Pro typické hodnoty koeficientů Encrit = 400×103 Vm-1, An = 0,56 a vsat300 = 72×103 ms-1, jsou na obr. 8.1B.6 vykresleny závislosti pohyblivosti a driftové rychlosti elektronů na velikosti intenzity elektrického pole. Když si uvědomíme, že podle (8.1B.12) je velikost driftové rychlosti přímo úměrná velikosti proudu, můžeme na obr. 8.1B.6b pozorovat oblast záporného diferenciálního odporu, který je důsledkem Gunnova jevu, jež je podstatou činnosti mnoha mikrovlnných aktivních prvků vyrobených z GaAs.

Obr. 8.1B.6aObr. 8.1B.6b
a) pohyblivost elektronůb) driftová rychlost elektronů
Obr. 8.1B.6Pohyblivost a driftová rychlost elektronů v GaAs v závislosti na velikosti intenzity elektrického pole při teplotách 300K a 400K

2.3     Generace a rekombinace nosičů

V předchozích odstavcích byla zmíněna funkce R, která popisuje četnost generace a rekombinace elektronů a děr. Uvažujme homogenně dotovaný polovodič, který je v tepelné rovnováze, tj. procesy generace a rekombinace vybuzené tepelnou energií jsou v dynamické rovnováze. Pro tento stav můžeme definovat rovnovážné koncentrace elektronů n0a děr p0, jež jsou svázány následujícím vztahem

n0p0=ni2, ( 8.1B.29 )

kde ni představuje intrinsickou koncentraci. Pokud polovodič vybudíme vnějším stimulem, rovnováha mezi procesy generace a rekombinace bude narušena a tím budou narušeny i koncentrace elektronů a děr z jejich rovnovážných hodnot. V rovnicích kontinuity (8.1B.11) je tohoto narušení je vyjádřeno funkcí R. Podle toho, který fyzikální jev je příčinou generace a rekombinace nosičů náboje, rozlišujeme

  • RSRH (Shockley-Read-Hall) přechody fononů
  • ROPT přechody fotonu
  • RSURF povrchové jevy
  • RAU Augerovy trojité přechody
  • RII nárazová ionizace

Jednotlivé formulace R kombinujeme pouhým sečtením R = RSRH + ROPT + RAU + RSURF + RII.

Dominantním mechanismem je Shockley-Read-Hall rekombinace a generace, která je modelována jako past mezi valenčním a vodivostním pásmem polovodiče (v zakázaném pásmu). Výsledná funkce má následující podobu

RSRH=npni2τp(n+n1)+τn(p+p1), ( 8.1B.30 )

kde τn a τp jsou doby života nosičů, n1 a p1 jsou koncentrace definované energetickou hladinou pasti. Pokud je hladina pasti právě uprostřed zakázaného pásma, jsou tyto koncentrace rovné intrinsické, tj. p1 = n1 = ni.

3        Analýza základních rovnic a okrajové podmínky

První sekci byl uveden soubor základních rovnic (8.1B.6), (8.1B.11) a (8.1B.13). Je důležité zmínit, že ve vztazích proudových hustot (8.1B.13) byly opomenuty složky proudu způsobené zúžením zakázaného pásma a teplotním gradientem, protože jsou vlivy těchto efektů považovány za zanedbatelné. Avšak je nutné si uvědomit, že vztahy pro proudy se tímto mohou stát nekorektní, pokud některý ze zmíněných efektů zásadně nabude na významu.

Použitím (8.1B.6) a substitucí rovnic proudové hustoty (8.1B.13) do rovnic kontinuity (8.1B.11) dostáváme systém tří parciálních diferenciálních rovnic (8.1B.31) s proměnnými Ψ, n a p.

(ψ)qε(npC)=0, ( 8.1B.31a )
(Dnnnμnψ)R=nt, ( 8.1B.31b )
(Dpp+pμpψ)R=pt. ( 8.1B.31c )

Za předpokladu, že parametry Dn, Dp, μn a μp jsou vždy kladné a pokud funkce R neobsahuje derivace druhého (nebo vyššího) řádu podle prostorových souřadnic a první (nebo vyšší) derivace koncentarací n a p podle času, potom jsou (8.1B.31b) a (8.1B.31c) parabolické parciální diferenciální rovnice. Rovnice (8.1B.31a) představuje eliptickou parciální diferenciální rovnici.

Pro matematickou analýzu potřebujeme znát počáteční odhad proměnných Ψ, n a p v doméně, ve které jsou tyto rovnice platné a okrajové podmínky. Tato ohraničená doména D je obecně trojrozměrná, tak jako prakticky všechny polovodičové struktury. Avšak v mnoha případech můžeme tuto doménu považovat za dvojrozměrnou, či dokonce za jednorozměrnou, čímž se analyzovaný problém zásadně zjednoduší.

Nechť ∂D značí dílčí hranice domény D. Tato může být z principu rozdělena na dvě části

D=DpDa, ( 8.1B.32 )

kde ∂Dp značí ty části hranice, které představují skutečné, fyzické, hranice, jako jsou kontakty a rozhraní. Naopak ∂Da značí umělé rozhraní, které zavádíme, například, za účelem vyčlenění dílčí struktury z rozměrného substrátu nebo oddělení sousedních prvků na jednom substrátu, tudíž nepředstavují fyzickou hranici.

Obr. 8.1B.7
Obr. 8.1B.72D geometrie planárního tranzistoru MOSFET [44]

Zmíněné rozdělení hranic ilustrativně popíšeme na idealizované geometrii 2D planárního MOS tranzistoru, viz. obr. 8.1B.7. Celá doména je reprezentována mnohoúhelníkem A-B-C-D-E-F-G-H-B. Rovnice (8.1B.31) platí pouze v subdoméně A-B-E-F-G-H-B. Izolant, ohraničený body B-C-D-E-B, můžeme charakterizovat buď Laplaceovou rovnicí pro elektrostatický potenciál (8.1B.33) nebo předpokladem absence nosičů náboje (8.1B.34).

(Ψ)=0 ( 8.1B.33 )
n=p=C=0 ( 8.1B.34 )

Avšak použití těchto rovnic znemožňuje určení proudu procházející bází a vliv náboje izolační vrstvy oxidu.

Hranice A-B, E-F, C-D a B-E můžeme pokládat za fyzické hranice představující tři kontakty a jedno rozhraní mezi polovodičem a izolantem. Hranice A-H, B-C, D-E, F-G a G-H pokládáme za umělé hranice. Tyto sice v reálné struktuře neexistují, ale usnadňují simulaci tím, že strukturu separují od okolního prostředí.

Fyzické hranice mohou být rozděleny do tří kategorii.

DP=DODSDI, ( 8.1B.35 )

kde ∂DO značí ohmický kontakt, ∂DS značí Schottkyho kontakt a ∂DI značí rozhraní s izolantem.

3.1     Ohmický kontakt

Nejjednodušší okrajová podmínka pro ohmický kontakt je Dirichletova okrajová podmínka pro elektrostatický potenciál

ψ(t)ψbψD(t)=0, ( 8.1B.36 )

kde ΨD(t) značí vnější přiložené napětí a Ψb značí vestavěný potenciál, který přímo souvisí s příměsi v objemu polovodiče a může být vyjádřen následujícím vztahem

ψb=kTqarsinh(C2ni). ( 8.1B.37 )

Pokud je polovodič dotován pouze donory nebo akceptory, můžeme použít zjednodušené vztahy

ND>>NAψbkTqln(ND2ni), ( 8.1B.38a )
NA>>NDψbkTqln(NA2ni). ( 8.1B.38b )

Dále se budeme zabývat okrajovými podmínkami pro koncentraci nosičů. Pokud budeme uvažovat tepelnou rovnováhu (8.1B.29) a absenci prostorového náboje na ohmickém kontaktu, platí následující vztahy.

npni2=0, ( 8.1B.39 )
npC=0. ( 8.1B.40 )

Tyto dvě rovnice mohou být vyjádřeny Dirichletovou okrajovou podmínkou pro elektrony (8.1B.41a) a díry (8.1B.41b).

n=C2+4ni2+C2, ( 8.1B.41a )
p=C2+4ni2C2. ( 8.1B.41b )

3.2     Schottkyho kontakt

Vzhledem k tomu, že fyzikální popis Shottkyho kontaktu není triviální, používá se pro účely simulace velice zjednodušený model. Pro elektrostatický potenciál můžeme psát Dirichletovu okrajovou podmínku

ψ(t)ψbψsψD(t)=0, ( 8.1B.42 )

kde ΨD(t) opět značí vnější přiložené napětí a Ψs je velikost Schottkyho bariery, což je charakteristický parametr daný druhem kovu a polovodiče Schottkyho kontaktu. Ψb opět značí vestavěný potenciál. Okrajová podmínka pro koncentraci nosičů vychází z proudové hustoty protékající kontaktem. Použitím teorie tepelné emise a difůze, lze odvodit následující okrajové podmínky pro závěrně polarizovaný kontakt

Jnn=qun(nC2+4ni2+C2), ( 8.1B.43a )
Jpn=qup(pC2+4ni2C2), ( 8.1B.43b )

kde n značí normálový vektor kolmý na ∂D, un,p značí rychlosti tepelné rekombinace.

3.3     Rozhraní polovodič izolant

Další kategorii tvoří rozhraní mezi polovodičem a izolačním materiálem. Na takovém rozhraní musí platit Gaußův zákon,

εsemψn|semεinsψn|ins=Qint, ( 8.1B.44 )

kde εsem a εins značí permitivitu polovodiče a izolantu. Qint je náboj na rozhraní. Pro tenké izolanty jako jsou vrstvy oxidu, můžeme složku elektrického pole v izolantu kolmou na rozhraní zanedbat, čímž se nám vztah zjednodušuje na

εsemψn|sem=Qint. ( 8.1B.45 )

Pokus dále zanedbáme náboj na rozhraní Qint je předchozí podmínka redukována na Neumannovu okrajovou podmínku pro elektrostatický potenciál

ψn=0. ( 8.1B.46 )

Pro nosiče platí, že složky proudu kolmé k rozhraní se musí rovnat rychlosti povrchové rekombinace RSURF, tj. musí platit

Jnn=qRSURF ( 8.1B.47a )
Jpn=qRSURF. ( 8.1B.47b )

Velmi často se povrchová rekombinace zanedbává, a tak získáváme jednodušší podmínky

Jnn=0, ( 8.1B.48a )
Jpn=0. ( 8.1B.48b )

3.4     Umělé rozhraní

Posledním případem jsou umělé rozhraní. Zde předpokládáme buď okrajové podmínky (8.1B.49a), (8.1B.49b) a (8.1B.49c), které zaručí, že doména je uzavřená nebo specifické hodnoty pro elektrostatický potenciál a koncentraci nosičů dané vnějšími vlivy.

ψn=0, ( 8.1B.49a )
nn=0, ( 8.1B.49b )
pn=0. ( 8.1B.49c )

Copyright © 2010 FEEC VUT Brno All rights reserved.