QuantumATK Forum
QuantumATK => General Questions and Answers => Topic started by: mbjavan on May 30, 2012, 10:01
-
please help me to plot optical spectrums, calculated and saved in some *.nc files at a single figure.
I need a script :-\
-
Please check here:
http://www.quantumwise.com/documents/tutorials/latest/SiliconOptical/index.html/chap.Si.Band.html#sect1.running.si
-
Dear ATK team,
I'm trying to get the absorption for my system. I've already calculated the dielectric constants and it is OK. My problem is when I apply the script of : http://www.quantumwise.com/documents/tutorials/latest/SiliconOptical/index.html/chap.Si.Band.html#sect1.running.si, I have RuntimeWarning: divide by zero encountered in divide.
I guess I have to change values in this script and this values must be extracted from optical spectrum. Could you please tell me what is the meaning
of the values in red, I mean what the corresponding values in the dielectric constant spectrum
ax.axis([180,1000,2.2,6.4])
ax.set_ylabel(r"$n$", size=16)
ax.tick_params(axis='x', labelbottom=False, labeltop=True)
ax = pylab.subplot(212)
ax.plot(l,alpha,'r')
ax.axis([180,1000,0,0.24])
Thx
-
Where is the warning? Probably from the line dividing by energy. If so, it's not so critical, it just means you have a zero energy.
180 and 1000 are wavelengths, computed from the energies on the line commented as "# Calculate the wavelength" and 2.2 and 6.4 and the other ones are needed because axis() take a 4-element list (xmin, xmax, ymin, ymax).
-
Dear Anders,
Actually it is critical because I can't have the absorption ! how we can avoid the problem. all the optical properties are extracted from dielectric constant and this later is correct. I mean how to avoid this problem.
Best regards,
-
Dear Ander,
Now it's working, please avoid my last post.
best