As discussed here:

You just need to set the material parameters, using a function like `SetSiliconParameters`

.

The scripts in `python_packages`

are not documented in the manual. You would have to look at the code. I will work on adding docstrings in the future.

https://github.com/devsim/devsim/issues/146

This function solves the Potential equation in a semiconductor region, assuming equilbrium.

This function sets the boundary condition at a contact for the potential equation in equilibrium.

Setups up the Potential, Electron, and Hole equations for drift-diffusion simulation.

Sets up an ideal ohmic contact for drift diffusion simulation.

Set the material parameters for your material. The existing equations should work for a single region device.

Please check the discussion here:

https://forum.devsim.org/t/simulation-of-simple-2d-sic-diode/328/2

One approach is to use the Step function. However, if you want to have an abrupt junction that is independent of meshing effects, you can have a P region and N region with an interface between them. Then the

`CreateSiliconSiliconInterface`

function would be useful.