Note: This discussion is about an older version of the COMSOL Multiphysics® software. The information provided may be out of date.

Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

Very simple semiconductor simulation

Kristian Storm
Hi! I'm trying to learn Comsol, but it seems have run into a wall. I've looked thoroughly through the semiconductor_diode model as well as the MOSFET model (in the showroom). However, as I implement these models myself, following the instructions point by point, the model diverges.

In general, I've noticed that I can get good convergence in a single Poisson simulation of a piece of semiconductor, with a given doping profile. As soon as I add transport of dilute species coupled with Poisson, it diverges.

In order to find my errors, I've created a very simple semiconductor model. It's basically a 1umX10um piece of silicon with a homogeneous doping concentration (but engineered such that it can be changed later). I provide an analytic initial value to an electrostatics module, which works fine and it solves the poisson equation. I then feed this solution to an electrostatic module coupled with two 'transport of dilute species' module, one for electrons and one for holes. As soon as this runs however, the model diverges.

I've attached the model (Comsol v. 4.1) and would be very grateful if someone with experience could point me in the right direction. I must be missing something fairly simple. If 'step 1' is solved, you see the potential for the first poisson solver (constant due to constant doping density). Step 2 then tries to solve poisson coupled with diffusion using the first solution as initial value, but diverges.

The essential variables are:
n0,p0,phi_init: analytic initial values fed to the first poisson solver.
phi0: potential solved by first poisson
phi,cn,cp: potential, electron and hole concentrations for second solver step.

I've been staring at this for a few days now and if I don't find the problem soon, I will have to give up using Comsol.

Thanks!


5 Replies Last Post 19/08/2013, 18:44 GMT-04:00
Kristian Storm
Posted: 6 years ago 11/05/2011, 04:05 GMT-04:00
I think I solved the problem. For future reference:
In the model builder, there is a menu if you click on the down arrow (v). There you can enable 'show more options'. In transport of diluted species I had to use the following settings:
- disable everything in 'consistent stabilization
- disable isotropic diffusion (this was already disabled)
- change convective term to conservative
- Change concentration discretization to quadratic

I think I solved the problem. For future reference: In the model builder, there is a menu if you click on the down arrow (v). There you can enable 'show more options'. In transport of diluted species I had to use the following settings: - disable everything in 'consistent stabilization - disable isotropic diffusion (this was already disabled) - change convective term to conservative - Change concentration discretization to quadratic

Posted: 6 years ago 06/06/2011, 05:02 GMT-04:00
Hello,
I am very interested in your problem. Did you find a solution?
My problem is mainly when I put a potential difference between the sides of the silicon bar.

If you are successful, thank you for your explanations.

Sincerely,

Gilles
Hello, I am very interested in your problem. Did you find a solution? My problem is mainly when I put a potential difference between the sides of the silicon bar. If you are successful, thank you for your explanations. Sincerely, Gilles

Carsten Madler
Posted: 5 years ago 06/12/2012, 17:40 GMT-05:00
I have a question about solving for the transport of holes and electrons using the Transport of Diluted Species Model. I am using Comsol 4.3 and in that model the species concentration has units of mol/cm^3 but the hole and electron concentration is obviously 1/m^3. Is there a way around that or was it different in previous versions of Comsol?

I have a question about solving for the transport of holes and electrons using the Transport of Diluted Species Model. I am using Comsol 4.3 and in that model the species concentration has units of mol/cm^3 but the hole and electron concentration is obviously 1/m^3. Is there a way around that or was it different in previous versions of Comsol?

Samuele Lilliu
Posted: 5 years ago 10/12/2012, 04:03 GMT-05:00
Hi,

Try to redefine nc as nc/(1 [mol]).

By the way I have similar convergence problems. I am trying to build a model for a pn solar cell. When I change the R and G terms in Transport of Diluted Species the simulation does not converge anymore. Switching from linear to quadratic does not help.

Please have a look at the attached file.

Cheers
Hi, Try to redefine nc as nc/(1 [mol]). By the way I have similar convergence problems. I am trying to build a model for a pn solar cell. When I change the R and G terms in Transport of Diluted Species the simulation does not converge anymore. Switching from linear to quadratic does not help. Please have a look at the attached file. Cheers


Ali Khalatpour
Posted: 4 years ago 19/08/2013, 18:44 GMT-04:00
Hi,
I believe I have the same problem in the semiconductor module. In simulation of Gunn effect we need to define field dependent mobility for electrons. This will bring serious convergence issues. I hope some one have done this before or have good suggestion for me.


Regards,
Al
Hi, I believe I have the same problem in the semiconductor module. In simulation of Gunn effect we need to define field dependent mobility for electrons. This will bring serious convergence issues. I hope some one have done this before or have good suggestion for me. Regards, Al

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.