A globally convergent algorithm for the solution of the steady‐state semiconductor device equations

An iterative method for solving the discretized steady‐state semiconductor device equations is presented. This method uses Gummel’s block iteration technique to decouple the nonlinear Poisson and electron‐hole current continuity equations. However, the main feature of this method is that it takes advantage of the diagonal nonlinearity of the discretized equations, and solves each equation iteratively by using the nonlinear Jacobi method. Using the fact that the diagonal nonlinearities are monotonically increasing functions, it is shown that this method has two important advantages. First, it has global convergence, i.e., convergence is guaranteed for any initial guess. Second, the solution of simultaneous algebraic equations is avoided by updating the value of the electrostatic and quasi‐Fermi potentials at each mesh point by means of explicit formulae. This allows the implementation of this method on computers with small random access memories, such as personal computers, and also makes it very attractiv...