Materials Science and Engineering/Doctoral review questions/Daily Discussion Topics/01132008

Carrier Model Equations

Density of States and Fermi Function

$g_c(E)=\frac{m_n*\sqrt{2m_n*(E-E_c)}}{{\pi }^2{\hslash }^3},E\ge E_c$

$g_v(E)=\frac{m_p*\sqrt{2m_p*(E_v-E)}}{{\pi }^2{\hslash }^3},E\le E_c$

 $f(E)=\frac{1}{1+e^{(E-E_F)/kT}}$

Carrier Concentration Relationships

$n=N_C\frac{2}{\sqrt{\pi }}{F}_{1/2}({\eta }_c)$

$p=N_V\frac{2}{\sqrt{\pi }}{F}_{1/2}({\eta }_c)$

$N_C=2{\left[\frac{m_n*kT}{2\pi {\hslash }^2}\right]}^{3/2}$

$N_V=2{\left[\frac{m_p*kT}{2\pi {\hslash }^2}\right]}^{3/2}$

$n=N_C{e}^{(E_F-E_C)/kT}$

$p=N_v{e}^{E_V-E_F)/kT}$

$n=n_i{e}^{(E_F-E_i)/kT}$

$p=n_i{e}^{(E_i-E_F)/kT}$

ni, np-Product, and Charge Neutrality

$n_i=\sqrt{N_c{N}_v}{e}^{-E_G/2kT}$

$np=n_i^2$

$p-n+N_D-N_A$

n,p, and Fermi Level Computational Relationships

$n=\frac{N_D-N_A}{2}+{[{\left(\frac{N_D-N_A}{2}\right)}^2+n_i^2]}^{1/2}$

$n\approx N_D$

$p\approx n_i^2/N_D$

$p\approx N_A$

$N_A\gg N_D,N_A\gg n_i$

$n\approx n_i^2/N_A$

 $E_i=\frac{E_c+E_v}{2}+\frac{3}{4}kT\ln \left(\frac{m_p*}{m_n*}\right)$

$E_F-E_i=kT\ln (n/n_i)=-kT\ln (p/n_i)$

$N_D\gg N_A,N_D\gg n_i$

$E_F-E_i=kT\ln (N_D/n_i)$

$N_A\gg N_D,N_A\gg n_i$

$E_i-E_F=kT\ln (N_A/n_i)$

Einstein Relationship

The Einstein relationship connects $D$ and $\mu$

The net carrier current needs to be zero in the case of equilibrium conditions

$J_{N|drift}+J_{N|diff}=q\mu nE+q{D}_N\frac{dn}{dx}=0$

$E=\frac{1}{q}\frac{d{E}_i}{dx}$

$n=n_i{e}^{(E_F-E_i)/kT}$

The change of the Fermi energy with respect to position is 0, $d{E}_F/dx=0$

$\frac{dn}{dx}=-\frac{n_i}{kT}{e}^{(E_F-E_i)/kT}\frac{d{E}_i}{dx}=-\frac{q}{kT}nE$

$(qnE){\mu }_n-(qnE)\frac{q}{kT}{D}_N=0$

The electric field, $E$, does not equal zero

 $\frac{D_N}{{\mu }_n}=\frac{kT}{q}$
 $\frac{D_P}{{\mu }_p}=\frac{kT}{q}$

The Built-in Potential (Vbi)

The electric field is related to the derivative of the potential with position

$E=-\frac{dV}{dx}$

$-\underset{-x_p}{\overset{x_n}{\int }}Edx=\underset{V(-x_p)}{\overset{V(x_n)}{\int }}dV=V(x_n)-V(-x_p)=V_{bi}$

The following is true under equilibrium

$J_N=q{\mu }_{n}nE+q{D}_N\frac{dn}{dx}$

Use the Einstein relationship

$E=-\frac{D_N}{{\mu }_n}\frac{dn/dx}{n}$

$E=-\frac{kT}{q}\frac{dn/dx}{n}$

$V_{bi}=-\underset{-x_p}{\overset{x_n}{\int }}Edx$

$V_{bi}=\frac{kT}{q}\underset{n(-x_p)}{\overset{n(x_n)}{\int }}\frac{dn}{n}$

$V_{bi}=\frac{kT}{q}\ln \left[\frac{n(x_n)}{n(-x_p)}\right]$

Consider nondegenerately doped step junction with $N_D$ and $N_A$ as the n- and p-side doping concentrations.

$n(x_n)=N_D$

$n(-x_p)=\frac{n_i^2}{N_A}$

 $V_{bi}=\frac{kT}{q}\ln \left(\frac{N_A{N}_D}{n_i^2}\right)$

LCAO and Tight-Binding

LCAO

A linear combination of atomic orbitals or LCAO is a quantum superposition of atomic orbitals and a technique for calculating molecular orbitals in quantum chemistry [1]. In quantum mechanics, electron configurations of atoms are described as wavefunctions. In mathematical sense, these wave functions are the basis set of functions, the basis functions, which describe the electrons of a given atom. In chemical reactions, orbital wavefunctions are modified, i.e. the electron cloud shape is changed, according to the type of atoms participating in the chemical bond.

It was introduced in 1929 by Sir John Lennard-Jones with the description of bonding in the diatomic molecules of the first main row of the periodic table, but had been used earlier by Pauling for H2+.

A mathematical description is

${\varphi }_i=c_1{\chi }_1+c_2{\chi }_2+c_3{\chi }_3+\cdots +c_n{\chi }_n$

or

${\varphi }_i=\sum _r{c}_{ri}{\chi }_r$

where ${\varphi }_i$ is a molecular orbital represented as the sum of n atomic orbitals ${\chi }_r$ (chi), each multiplied by a corresponding coefficient $c_r$. The coefficients are the weights of the contributions of the n atomic orbitals to the molecular orbital.

The orbitals are thus expressed as linear combinations of basis functions, and the basis functions are one-electron functions centered on nuclei of the component atoms of the molecule. The atomic orbitals used are typically those of hydrogen-like atoms since these are known analytically i.e. Slater-type orbitals but other choices are possible like Gaussian functions from standard basis sets.

Tight-Binding

In the tight binding model, it is assumed that the full Hamiltonian $H$ of the system may be approximated by the Hamiltonian of an isolated atom centred at each lattice point. The atomic orbitals ${\psi }_n$, which are eigenfunctions of the single atom Hamiltonian $H_{at}$, are assumed to be very small at distances exceeding the lattice constant. This is what is meant by tight-binding. It is further assumed that any corrections to the atomic potential $\mathrm{\Delta }U$, which are required to obtain the full Hamiltonian $H$ of the system, are appreciable only when the atomic orbitals are small. A solution to the time-independent single electron Schrödinger equation $\mathrm{\Phi }$ is then assumed to be a linear combination of atomic orbitals ${\psi }_n$

$\mathrm{\Phi }(\overrightarrow{r})=\sum _{n,\overrightarrow{R}}{b}_{n,\overrightarrow{R}}{\psi }_n(\overrightarrow{r}-\overrightarrow{R})$

where n refers to the n-th atomic energy level and ${\overrightarrow{R}}^{}$ is an atomic site in the crystal lattice.

Using this approximate form for the wavefunction, and assuming only the m-th atomic energy level is important for the m-th energy band, the Bloch energies ${\epsilon }_m$ are of the form

${\epsilon }_m(\overrightarrow{k})=E_m-\frac{{\beta }_m+\sum _{\overrightarrow{R}\ne 0}{\gamma }_m(\overrightarrow{R})e^{i\overrightarrow{k}\cdot \overrightarrow{R}}}{b_{m,\overrightarrow{R}}+\sum _{\overrightarrow{R}\ne 0}{\alpha }_m(\overrightarrow{R})e^{i\overrightarrow{k}\cdot \overrightarrow{R}}}$,

where $E_m$ is the energy of the $m$th atomic level,

${\beta }_m=-\int {\psi }_m^{*}(\overrightarrow{r})\mathrm{\Delta }U(\overrightarrow{r})\mathrm{\Phi }(\overrightarrow{r})d\overrightarrow{r}$,

${\alpha }_m(\overrightarrow{R})=\int {\psi }_m^{*}(\overrightarrow{r})\mathrm{\Phi }(\overrightarrow{r}-\overrightarrow{R})d\overrightarrow{r}$,

and

${\gamma }_m(\overrightarrow{R})=-\int {\psi }_m^{*}(\overrightarrow{r})\mathrm{\Delta }U(\overrightarrow{r})\mathrm{\Phi }(\overrightarrow{r}-\overrightarrow{R})d\overrightarrow{r}$

are the overlap integrals.

The tight-binding model is typically used for calculations of electronic band structure and energy gaps in the static regime. However, in combination with other methods such as the random phase approximation (RPA) model, the dynamic response of systems may also be studied.

Phase Diagram

The Iron–Iron Carbide (Fe–Fe₃C) phase diagram. The percentage of carbon present and the temperature define the phase of the iron carbon alloy and therefore its physical characteristics and mechanical properties. The percentage of carbon determines the type of the ferrous alloy: iron, steel or cast iron

A phase diagram for a binary system displaying a eutectic point.

Complex types of phase diagrams can be constructed, particularly when more than one pure component is present. In that case concentration becomes an important variable. Phase diagrams with more than two dimensions can be constructed that show the effect of more than two variables on the phase of a substance. Phase diagrams can use other variables in addition to or in place of temperature and pressure and composition, for example the strength of an applied electrical or magnetic field and they can also involve substances that take on more than just three states of matter.

One type of phase diagram plots temperature against the relative concentrations of two substances in a binary mixture called a binary phase diagram, as shown at right. Such a mixture can be either a solid solution, eutectic or peritectic, among others. These two types of mixtures result in very different graphs. A textbook example of a eutectic phase diagram is that of the olivine (forsterite and fayalite) system.

Another type of binary phase diagram is a boiling point diagram for a mixture of two components, i. e. chemical compounds. For two particular volatile components at a certain pressure such as atmospheric pressure, a boiling point diagram shows what vapor (gas) compositions are in equilibrium with given liquid compositions depending on temperature. In a typical binary boiling point diagram, temperature is plotted on a vertical axis and mixture composition on a horizontal axis.

A simple example diagram with hypothetical components 1 and 2 in a non-azeotropic mixture is shown at right. The fact that there are two separate curved lines joining the boiling points of the pure components means that the vapor composition is usually not the same as the liquid composition the vapor is in equilibrium with. See Vapor-Liquid Equilibrium for a fuller discussion.

In addition to the above mentioned types of phase diagrams, there are thousands of other possible combinations. Some of the major features of phase diagrams include congruent points, where a solid phase transforms directly into a liquid. There is also the peritectoid, a point where two solid phases combine into one solid phase during heating. The inverse of this, when one solid phase transforms into two solid phases during heating, is called the eutectoid.

External Link

• A collection of about 150 alloy phase diagrams and some PT diagrams