- 4.1.8.1 Artificial Boundaries
- 4.1.8.2 Semiconductor/Metal

- 4.1.8.3 Insulator/Metal
- 4.1.8.4 Semiconductor/Insulator
- 4.1.8.5 Semiconductor/Semiconductor

For the solution of the semiconductor equations a closed domain is required. Boundary conditions for the unknowns in a certain segment must be specified at the boundaries of that segment. Depending on the number of segments (one or two) and the type of the quantity of the unknown (distributed or non-distributed) several model types can be separated [267]. However, in order to illustrate the basic boundary conditions such a differentiation is not needed and therefore not provided here.

The simulation domain usually includes only a single device. In order to separate it from the neighboring devices artificial boundaries, which have no physical origin, are introduced. Through the Neumann boundary conditions a self-contained domain is guaranteed. The fluxes across the boundary are set to zero:

is an outward-orientated vector normal to the boundary. (4.28) define the boundary conditions for a semiconductor segment, while for an insulator the first two conditions are sufficient.

Ohmic contacts are defined by Dirichlet boundary conditions: the contact potential , the carrier contact concentration and , and in the case of a HD simulation the carrier contact temperatures and are fixed. The metal quasi-Fermi level (which is specified by the contact potential ) is equal to the semiconductor quasi-Fermi level. The model assumes charge-neutrality on the boundary. The contact potential at the semiconductor boundary is:

where is the built-in potential:

is the net concentration of dopants and other charged defects at the contact boundary. The variables and are defined by [268]

(4.32) | |||

(4.33) |

where and are the effective density of states.

The carrier concentrations in the semiconductor are pinned to the
equilibrium carrier concentrations at the contact:

(4.34) | |||

(4.35) |

This model works properly, only if a high doping ( cm) in the semiconductor is guaranteed.

The carrier temperatures and are set equal to the lattice temperature :

If a contact temperature is specified, the lattice temperature is calculated using and a thermal resistance . The thermal heat flow density at the contact boundary reads:

If no thermal resistance is specified, an isothermal boundary condition is assumed and the lattice temperature is set equal to the contact temperature:

In case of a DD simulation with self-heating an additional thermal energy is accounted for. It is generated, when the carriers transcend the potential difference between the conduction or valence band and the metal quasi-Fermi level. The energy equation reads:

The expression denotes the surface divergence of the thermal heat flux at the boundary considered. In case of a HD simulation with self-heating the thermal heat flow across the boundary is accounted for self-consistently.

It is possible to include an electric line resistance of the contact using:

(4.40) |

with the applied terminal voltage and the current through the contact.

At the Schottky contact mixed boundary conditions apply. The contact potential , the carrier contact concentration and , and in the case of a HD simulation, the contact carrier temperatures and are fixed. The semiconductor contact potential is the difference between the metal quasi-Fermi level and the metal work function difference :

where

The difference between the conduction band energy and the metal work function energy gives the work function difference energy , which is the so called barrier height of the Schottky contact, also denoted as . The applied boundary conditions are:

(4.43) | |||

(4.44) |

where and are the recombination velocities. The carrier concentrations read:

(4.45) | |||

(4.46) |

where is the Fermi level in the semiconductor. If the values of and are set to zero, the current through the Schottky contact is suppressed. The Schottky contact boundary conditions for the carrier temperatures and and for the lattice temperature are similar to those which apply for the Ohmic contact: (4.36) and (4.37) or (4.38).

For the equilibrium situation, the quasi-equilibrium concentrations and can be rewritten:

(4.47) | |||

(4.48) |

The recombination velocities are calculated from:

(4.49) |

where is the Richardson constant.

Typical values for the Schottky barrier heights of common
semiconductor barriers are listed in Table 4.1. Values
measured by different methods (*I-V* and *C-V*) for
n-GaN/metal differ significantly. This is caused by defects in the
surface region, which enhance the tunneling effect and therefore have
an impact on the Richardson constant [269,270]. The exact
value of the barrier height depends also on orientation, stress, and
polarity of the GaN layer [271]. No data is provided for
contacts on InN as most metals show ohmic
behavior [272]. Schottky barrier height varies with
annealing temperature: generally it is reduced after annealing.

Materials | [eV] | Ref. | Materials | [eV] | Ref. |

n-GaN/Au | 0.87-1.1 | [273,274] | InAlN/Pd | 1.56 | [275] |

n-GaN/Ni | 0.95-1.13 | [273] | In Al N/Ni | 0.75 | [276] |

p-GaN/Ni | 2.68-2.87 | [277] | p-In Ga N/Ni | 0.39 | [278] |

Al Ga N/Ni | 0.94-1.24 | [279] | In Ga N/Pt | 0.62 | [270] |

Al Ga N/Ni | 1.26 | [280] | In Ga N/Ni | 1.39 | [270] |

Al Ga N/Ni | 1.02-1.30 | [279] | Al In GaN/Ni | 0.98-0.93 | [281] |

A model similar to the Schottky contact model is used to calculate the insulator contact potential. The semiconductor contact potential is the difference of the metal quasi-Fermi level and the metal work function difference potential similar to (4.41) and (4.42). The lattice temperature is set equal to the contact temperature (4.38).

In the absence of surface charges the normal component of the dielectric displacement and the potential are continuous:

(4.50) |

In the presence of surface charges along the interface, the dielectric displacement obeys the law of Gauß:

(4.51) |

where is the interface charge density. This distributed surface charge is introduced to describe the Fermi level pining at the surface due to a high trap density of states.

At the semiconductor/insulator interface the current densities and heat fluxes normal to the interface vanish.

(4.52) |

The lattice temperature at the interface is continuous.

The calculation of the electrostatic potential at the interfaces between two semiconductor segments is similar to that for semiconductor/insulator interfaces:

(4.53) |

The interface charge can take any value including zero. The way to calculate the exact value depending on the materials is discussed in Section 4.5. The subscripts are used to distinguish between the two semiconductor segments on both sides of of the interface.

To calculate the carrier concentrations and the carrier temperatures at the interface, three different approaches are considered:

- a model with continuous quasi-Fermi level across the interface (CQFL),
- a thermionic emission model,
- a thermionic field emission model.

Each model can be specified for electrons and holes for each semiconductor/semiconductor interface.

In the following denotes the normal to the interface component of the current density , the energy flux density component, and the difference in the conduction band edges . The effective electron mass is denoted by . The subscripts denote the semiconductor segment . Only the equations for electrons are given, those for holes can be deduced accordingly.

Dirichlet boundary conditions are applied, when the CQFL is used. The
carrier concentrations are determined so that the quasi-Fermi levels
remains continuous across the interface [134].

(4.54) | |||

(4.55) | |||

(4.56) | |||

(4.57) |

The model is suitable for use at homojunctions. In Chapter 5 it is applied at the interface between segments of the same material but different characteristics (e.g. material quality). However, it is not suitable for heterojunctions, as it ignores the band gap alignment.

To consider the band gap alignment as typical in a heterostructure, the
thermionic emission or thermionic field emission interface models must
be used.

(4.58) guarantees the particle conservation at the interface. (4.60) relates the energy and current fluxes on both sides of the interface. The second term corresponds to the conversion of potential energy into kinetic energy, when the electrons descend the step . (4.60) describes the thermionic emission across the energy barrier by relating the interface current to the quasi-Fermi level on both sides of the interface [134,282]. The exponential factor in (4.59) and (4.61) shows that only the fraction of the total electron concentration on Side 1, which is energetically above the barrier, is able to cross the potential step. In both (4.59) and (4.61) the same emission velocities occur:

(4.62) |

Tunneling can be included in this model by assuming that the tunnel effect lowers the barrier :

(4.63) |

The barrier lowering depends on the electric field orthogonal to the interface and the effective tunneling length . The latter is material dependent: for Al Ga As/In Ga As it is estimated to 7 nm [283] and for Al Ga N/GaN to 3 nm [138].

S. Vitanov: Simulation of High Electron Mobility Transistors