Modellierung

   Prinzipiell sind bei der Modellierung des Hochfeldtransports in HL mit Mehrtalbandstruktur zwei Ansätze möglich: Die Einführung der Ladungsträgergrößen der einzelnen Täler (Konzentration, Energie etc.) als unabhängige Variablen in den Transportgleichungen. Dies impliziert Formulierungen der abhängigen physikalischen Parameter wie Beweglichkeit etc. für die talspezifischen Variablen als Funktionen der externer Größen [28]. Die Teilchenübergänge zwischen den Tälern werden so zu Generations-/Rekombinationstermen, für die zusätzliche Modelle gefunden werden müssen, die kaum oder nur schwer experimentell unterlegbar sind. Dieser Ansatz erhöht nicht nur die Anzahl der Unabhängigen und Differentialgleichungen sondern auch die Modellkomplexität enorm und wird daher in der Praxis kaum verfolgt. Der alternative, auf den ersten Blick vielleicht physikalisch weniger motivierte Weg besteht in der Beibehaltung des Bildes eines einzelnen Elektronengases unter Berücksichtigung der Mehrtalbandstruktur in der Modellierung durch effektive Größen des gesamten Ladungsträgerensembles. Diese Vorstellung wird auch in der Folge benutzt.

Die Modellierung des Feldeinflusses auf die Beweglichkeit im Fall des DD Transportmodells  geschieht allgemein durch eine Formulierung $\mu_{\mathrm{}}^{} = \mu_{\mathrm{}}^{}(\mu_{\mathrm{0}}^{},F)$.Das heißt, man nimmt eine Trennung vor in die Nullfeldbeweglichkeit $\mu_{\mathrm{0}}^{}$ inklusive aller Gitter-, Legierungs- und Dotierungsabhängigkeiten und eine explizite Abhängigkeit vom Betrag der treibenden elektrischen Kraft

$$\frac{\vec{\vec{F}_{\mathrm{ext}}}}{e_0} = \vec{F} - \frac{\nabla n\,k_{\mathrm{B}}\,T}{e_0\,n}\,,$$ (6.68)

die unter homogenen Bedingungen dem elektrischen Feld $\vec{F}$ gleich ist.

Meist wird für III-V HL der Ansatz [131]

$$\mu_{\mathrm{}}^{} = \frac{\mu_{\mathrm{0}}^{}+\displaystyl... ...} {1+\left(\displaystyle\frac{F}{F_{\mathrm{crit}}}\right)^m}$$ (6.69)

gebraucht, welcher durch die Dotierungsunabhängigkeit der Sättigungsgeschwindigkeit motiviert ist. (6.69) ist eine rein empirische Beschreibung der effektiven Beweglichkeit, die den Zwischentaltransfer implizit enthält. Die Parameter $F_{\mathrm{crit}}$ und m bestimmen die Lage des Maximums von $v_{\mathrm{d}}$ und die Steilheit des Übergangs vom linearen Verlauf $\mu_{\mathrm{0}}^{}\,F$ zu $v_{\mathrm{sat}}$. Eine konzeptuell gleiche Beschreibung mit zusätzlichen Abhängigkeiten der Modellparameter $v_{\mathrm{sat}}$ und m von $\mu_{\mathrm{0}}^{}$ wurde in [223] vorgeschlagen, um die Meßwerte genauer wiederzugeben. Eine andere Berechnung für die effektive Gesamtbeweglichkeit basierend auf der Beschreibung eines Zweitalsystems mit feldabhängiger Kopplung unter Einhaltung von Teilchen- und Energieerhaltung wurde in [131] vorgeschlagen.

Die Formulierung der Beweglichkeit im Rahmen des hydrodynamischen (HD)  oder Energietransportmodells wird allgemein durch einen Ansatz $\mu_{\mathrm{}}^{} = \mu_{\mathrm{}}^{}(\mu_{\mathrm{0}}^{},w_{\mathrm{}}^{})$ vorgenommen. Die Abhängigkeit von $\mu_{\mathrm{}}^{}$ vom elektrischen Feld ist hier implizit in der Feldabhängigkeit der Ladungsträgerenergie $w_{\mathrm{}}^{}(F)$ enthalten. Dieser Zusammenhang ist nichtlokal, wenn $w_{\mathrm{}}^{}$ aus der Lösung der Energiebilanzgleichung in der HD Simulation resultiert. Zum Zweck der Modellierung wird Homogenität des Feldes angenommen und man hat somit wieder einen eindeutigen, üblicherweise monotonen Zusammenhang $w_{\mathrm{}}^{}(F)$.

Unter speziellen Bedingungen können HD und DD Formulierung von $\mu_{\mathrm{}}^{}$ eindeutig ineinander übergeführt werden. Nimmt man an, daß das Abklingen der mittleren Impulskomponente in Feldrichtung p zum Gleichgewichtswert Null mit der Impulsrelaxationszeit  $\tau_{\mathrm{m}}$ erfolgt,

$${\left\langle {\frac{d p}{d t}}\right\rangle}=-\frac{p}{\tau_{\mathrm{m}}}\,,$$ (6.70)

und berücksichtigt man die Definition der Beweglichkeit (6.16), so stellt die Impulsbilanz (6.23) unter stationären homogenen Bedingungen die Stromrelation

$$\vec{j} = e_0\,n\,\mu_{\mathrm{}}^{}\,\vec{F}$$ (6.71)

dar. Nimmt man in gleicher Weise eine Relaxation der mittleren Energie zum Gleichgewichtswert $w_{\mathrm{0}}^{}$ mit der Energierelaxationszeit  $\tau_{w}$ an,

$${\left\langle {\frac{d E}{d t}}\right\rangle}=-\frac{w_{\mathrm{}}^{}-w_{\mathrm{0}}^{}}{\tau_{w}}\,,$$ (6.72)

so lautet die stationäre Energiebilanz (6.22)

$$\mbox{${\vec{\vec{j}} \cdot \vec{F}}$} = n\,\frac{w_{\mathrm{}}^{}-w_{\mathrm{0}}^{}}{\tau_{w}}\,.$$ (6.73)

Die Größen $\tau_{\mathrm{m}}$ und $\tau_{w}$ sind im Prinzip von den im HD Modell berücksichtigten ersten vier Momenten der Distribution ($n,v_{\mathrm{d}},w_{\mathrm{}}^{},S$) abhängig, und können im Rahmen eines selbstkonsistenten Systems daher nicht frei gewählt werden. Hänsch hat gezeigt, daß man unter Verwendung eines DD Modells für $\tau_{\mathrm{m}}$ bzw. $\mu_{\mathrm{}}^{}(F)$, das den Sättigungsverlauf der Driftgeschwindigkeit in Si gut beschreibt, den einfachen Zusammenhang

$$\mu_{\mathrm{}}^{}(w_{\mathrm{}}^{}) = \frac{\mu_{\mathrm{0... ...qquad \eta = \frac{\mu_0}{e_0\,\tau_{w}\,v_{\mathrm{sat}}^2}$$ (6.74)

für die HD Beweglichkeit mit einer konstanten Energierelaxationszeit erhält [77].

Die Gleichungen (6.71) und (6.73) geben den expliziten $w_{\mathrm{}}^{}(F)$ Zusammenhang unter Verwendung der erwähnten Modelle für $\mu_{\mathrm{}}^{}(F)$ und $\tau_{w}$. Weitab vom Gleichgewicht ist $\mu_{\mathrm{}}^{}$ umgekehrt proportional zu $w_{\mathrm{}}^{}$, $w_{\mathrm{}}^{}$ steigt linear mit F, was einer Ladungsträgerbewegung mit konstanter Driftgeschwindigkeit $v_{\mathrm{sat}}$ entspricht.