solve_init compute the initial solution of the 3D Poisson equation, when dealing with quasi-1D channels, displaced along the z-axis. In particular, it solves the Laplace equation over the whole domain, imposing Dirichlet boundary conditions in correspondence of the gate, and in correspondence of the atoms belonging to the doped reservoirs (defined by means of the dope_reservoirs command). In particular, the potential E* to be imposed in the reservoirs is computed assuming complete ionization and solving the equation
where N1D is the 1D density of states self-computed by the solve_init procedure, while the left hand term is the fixed charge, equal to the molar fraction of the reservoir.
If zmean is the mean z-coordinate of the structure, we also make the assumption that atoms belonging to a reservoir and with z-coordinate > zmean have an electrochemical potential equal to channel.mu2, and equal to channel.mu1 if z-coordinate < zmean.
Here an example of the solve_init command (the command is highlighted in red)
from NanoTCAD_ViDES import *
# I create the grid
# Now I define the gate regions
# I take care of the solid
# I dope the reservoirs
# Here I compute the initial solution
# Now that I have an initial solution, I can start computing the
# self consistent problem