<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
Dear Konstantin,<br>
If I well understand you would like to draw the difference electron
density for the Ti d level for TiN.<br>
The Ti is +3 (d1 configuration), so you expect 1 electron into the d
level (you have one Titanium per primitive unit cell).<br>
You should first find the energy separation between the N-P states and
Ti-d level (in Ry) and you can use the DOS for this (in Ry) or look
into the scf1 what is the energy of the 9th level.<br>
I assume that you do non spin polarized calculation, so you will find
into the scf1<br>
1 level for the Ti 3s (2 electrons)<br>
3 levels for the Ti 3p (6electrons)<br>
1 level for the Ni 2s (2 electrons)<br>
3 levels for the Ni 3p (6 electrons)<br>
and then the 5 Ti d levels<br>
The Emin should be just below this last level.<br>
If you keep the same integration scheme (tetra), then the number of
electron is OK : 17. The number of electron is used to calculate the
Fermi level and you should check that the Fermi level is still the same
as in the SCF calculation.<br>
Emin is just use to decide from what level to what level you generate
the clmval file (be careful, your clmval file will not be correct to
restart an scf cycle).<br>
<br>
You can effectively do an only one k-point calculation. You do first
the scf with the dense k-mesh. Then you save the calculation and
generate a new k-list file (or klist_band). <br>
You do the calculation of the vector file with lapw1 for only this
k-point and then generate the clmval file from this vector using ROOT
option for efmod in th in2, keeping Emin and ne at the same values.<br>
<br>
regards<br>
Florent<br>
<br>
<br>
<br>
<br>
<br>
<br>
Konstantin Nefedev a écrit :
<blockquote type="cite">
<div>Dear All, please could you help me with calculation of
the difference electron density for TiN. This is end of the<br>
exercise 1 from wien2k workshop tasks (Singapore 6-9 july). I use old<br>
version of wien2k and I have not Xcrysden for visualization of
structure. <br>
<br>
The content of *.in2 file is<br>
<br>
TOT<br>
0.350 17.0 0.50 0.05 EMIN, NE, ESEPERMIN, ESEPER0<br>
TETRA 0.000<br>
0 0 4 0 4 4 6 0 6 4<br>
0 0 4 0 4 4 6 0 6 4<br>
12. GMAX<br>
NOFILE <br>
<br>
<br>
<br>
Emin I have took from *.scf1 file as a bottom<br>
<br>
<br>
<br>
LOCAL ORBITAL<br>
E( 2)= 0.3500 E(BOTTOM)= 0.350 E(TOP)= -200.000<br>
APW+lo<br>
<br>
<br>
The content of *.inst file <br>
<br>
<br>
Ti<br>
Ar 2<br>
3, 2,2.0 P<br>
3, 2,0.0 N<br>
4,-1,1.0 N<br>
4,-1,1.0 N<br>
N<br>
He 3<br>
2,-1,1.0 N<br>
2,-1,1.0 N<br>
2, 1,1.0 N<br>
2, 1,0.0 N<br>
2,-2,2.0 N<br>
2,-2,0.0 N<br>
****<br>
****<br>
<br>
I have use x lstart [-sigma] and for (100) direction have chose *.in5
in <br>
following form<br>
<br>
-1 -1 0 4<br>
-1 3 0 4<br>
3 -1 0 4<br>
3 2 3<br>
100 100<br>
DIFF SUB<br>
ANG VAL NODEBUG<br>
ORTHO<br>
</div>
<div>What is wrong?</div>
<div><br>
But I can not to obtain the right "occupied Ti-d band" plot, in form as
I can see in<br>
exercise. My value of density charge is 5, but I see it must be 0.08 in<br>
maximum.</div>
<div><br>
How it is possible to modify *.in2 for obtain the density in a certain<br>
E-interval? I have search in UG for case.in2 and have not found
information about<br>
such modification...<br>
<br>
Which way I need use for plotting of charge density only of
d-electrons? <br>
<br>
Other question concerns of band structure. Is it possible to leave only<br>
interesting direction in BZ for example Gamma L and X, and cut other?<br>
</div>
<div>Thank you for replay.</div>
<div><br>
Regards,<br>
Konstantin</div>
<br clear="all">
<pre wrap="">
<hr size="4" width="90%">
_______________________________________________
Wien mailing list
<a class="moz-txt-link-abbreviated" href="mailto:Wien@zeus.theochem.tuwien.ac.at">Wien@zeus.theochem.tuwien.ac.at</a>
<a class="moz-txt-link-freetext" href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a>
</pre>
</blockquote>
<br>
<br>
<pre class="moz-signature" cols="80">--
-------------------------------------------------------------------------
| Florent BOUCHER | |
| Institut des Matériaux Jean Rouxel | <a class="moz-txt-link-freetext" href="Mailto:Florent.Boucher@cnrs-imn.fr">Mailto:Florent.Boucher@cnrs-imn.fr</a> |
| 2, rue de la Houssinière | Phone: (33) 2 40 37 39 24 |
| BP 32229 | Fax: (33) 2 40 37 39 95 |
| 44322 NANTES CEDEX 3 (FRANCE) | <a class="moz-txt-link-freetext" href="http://www.cnrs-imn.fr">http://www.cnrs-imn.fr</a> |
-------------------------------------------------------------------------</pre>
</body>
</html>