$Withionmove

When the system has ion to perform as drift-diffusion equations, we solve the time dependent drift-diffusion for slow ion move simulation. Since the ion moves may not be governed fermi-level concept. We simply treat is a tradiational drift-diffusion equations.

$\frac{\partial M_{ion}}{\partial t} = \nabla \left( (q_{sign}) e\mu M_{ion} \vec{E} - q D_{M} \nabla M_{ion} \right)$

Ideally, the ion density is given by initial setting. The total ion number should be fixed. The program is aim to model multi-ions drift-diffusion. The command is as following.

$Withionmove
$N_sweep$ $N_output$
$T_{1,stop}$  $dt_{1}$ $T_{2,stop}$  $dt_{2}$ ... $T_{N_{sweep},stop}$  $dT_{M_{sweep}}$
$N_{ions} ~~ q_{sign,1} ~~ q_{sign,2} ~~ q_{sign,3}  ....~~ q_{sign,N_{ions}}$
$Parameter_{type,1} ~~ M_{ions,1} \mu p3 p4 ....$   Parameters of the 1 layer
$Parameter_{type,2} ~~ M_{ions,2} \mu p3 p4 ....$   Parameters of the 2 layer
$Parameter_{type,3} ~~ M_{ions,3} \mu p3 p4 ....$   Parameters of the 3 layer
 ....
$Parameter_{type,N} ~~ M_{ions,N} \mu p3 p4 ....$   Parameters of the $totalregion layer