PLaSK User Manual

Page Contents


<electrical solver="ShockleyCyl">

Corresponding Python class: electrical.shockley.ShockleyCyl.

Two-dimensional phenomenological solver in cylindrical geometry based on Shockley equation and using finite-element method.

  • name (required) – Solver name.

Geometry for use by this solver.

  • ref (required) – Name of a Cylindrical geometry defined in the <geometry> section.

Rectangular2D mesh used by this solver.

  • ref (required) – Name of a Rectangular2D mesh defined in the <grids> section. (mesh)
  • include-empty – Should empty regions (e.g. air) be included into computation domain? (bool, default is ‘no’)

Voltage boundary conditions. See subsection Boundary conditions.


Configuration of the self-consistent loop.

  • maxerr – Maximum allowed current density error. (float [%], default 0.05 %)

Configuration of the matrix solver.

  • algorithm – Algorithm used for solving set of linear positive-definite equations. (‘cholesky’, ‘gauss’, or ‘iterative’, default is ‘cholesky’)
  • itererr – Maximum allowed residual error for the iterative algorithm. (float, default 1e-08)
  • iterlim – Maximum number of iterations for the iterative algorithm. (int, default 10000)
  • logfreq – Number of iterations after which the progress is logged. (int, default 500)

Configuration of the effective model of p-n junction.

  • beta# – Junction coefficients. This is an inverse of the junction thermal voltage. (float [1/V])
  • js# – Reverse bias current densities. (float [A/m2])
  • pnjcond – Initial vertical conductivity of the junctions. (float [S/m], default 5.0 S/m)

Properties of the contact layers.

  • pcond – p-contact conductivity. (float [S/m], default 5.0 S/m)
  • ncond – n-contact conductivity. (float [S/m], default 50.0 S/m)