Thank you Prof. Giuseppe,<span name="Paolo Giannozzi" class="gD"> and Prof. Paolo</span>, i will try from MT pseudopotential and try to built it myself. <br><br><div class="gmail_quote">On Mon, Sep 24, 2012 at 1:57 PM, Giuseppe Mattioli <span dir="ltr">&lt;<a href="mailto:giuseppe.mattioli@ism.cnr.it" target="_blank">giuseppe.mattioli@ism.cnr.it</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
Dear Brahma<br>
I&#39;m afraid I&#39;m not. I&#39;ve never employed Bachelet-Hamann-Schluter pseudopotentials in my calculations<br>
and I do not know if there are still codes that generate them. Why do not you try with Martins-<br>
Troullier pseudopotentials? There is plenty of them in the ps library, and you can build your own<br>
pseudopotential with ld1.x<br>
Yours<br>
Giuseppe<br>
<br>
On Friday 21 September 2012 19:43:28 Bramha Pandey wrote:<br>
&gt; Thank you very much Dear Mattioli Sir , it  works.<br>
<div class="im">&gt;<br>
&gt; Can please give a clue to obtained the Al.pz-bhs.UPF (Al LDA with bhs<br>
&gt; funtionals essentially NormConversing) pseudopotential as i am very poor to<br>
&gt; generate it myself and not obtained  it in pslibrary.<br>
</div>&gt; it&#39;s a very bad habit i know but i have no option except taking help at<br>
&gt; this stage for my research.<br>
&gt;<br>
&gt; thanks in advanced for your kind help.<br>
&gt;<br>
&gt; On Fri, Sep 21, 2012 at 7:45 PM, Giuseppe Mattioli &lt;<br>
&gt;<br>
&gt; <a href="mailto:giuseppe.mattioli@ism.cnr.it">giuseppe.mattioli@ism.cnr.it</a>&gt; wrote:<br>
&gt; &gt; Dear Brahma<br>
&gt; &gt;<br>
&gt; &gt; &gt;  &amp;ELECTRONS<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;                     conv_thr = 1.0d-11 ,<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;  /<br>
&gt; &gt;<br>
&gt; &gt; You require a very strict convergence threshold on wfcs...<br>
&gt; &gt; 1.0d-7 may be sufficient for scf calculations.<br>
&gt; &gt;<br>
&gt; &gt; &gt;             nqx1 = 1, nqx2 = 1, nqx3 = 1,<br>
&gt; &gt;<br>
&gt; &gt; You should check the convergence of such mesh. The above 1 1 1 could be<br>
&gt; &gt; quite poor.<br>
&gt; &gt;<br>
&gt; &gt; HTH<br>
&gt; &gt;<br>
&gt; &gt; Giuseppe<br>
&gt; &gt;<br>
&gt; &gt; On Friday 21 September 2012 15:42:03 Bramha Pandey wrote:<br>
&gt; &gt; &gt; Dear Developers and Users,<br>
&gt; &gt; &gt; I was calculating the band structure of ZB-InN using &#39;hse&#39; hybrid<br>
&gt; &gt; &gt; functional. It gives the reasonable band gap opening but the run was<br>
&gt; &gt; &gt; goes on and can&#39;t stop. Here i am giving some sort of it&#39;s output.<br>
&gt; &gt; &gt; Firstly sorry for giving long output data...but have no choice...<br>
&gt; &gt; &gt; Please give me a ray of hope to converges of hse calculation in pw.x.<br>
&gt; &gt; &gt; input....<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; &amp;control<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;     calculation = &#39;scf&#39;<br>
&gt; &gt; &gt;     restart_mode=&#39;restart&#39;,<br>
&gt; &gt; &gt;     pseudo_dir = &#39;/pseudo/&#39; ,<br>
&gt; &gt; &gt;     outdir=&#39;/tmp/&#39; ,<br>
&gt; &gt; &gt;     prefix = &#39;innexx&#39; ,<br>
&gt; &gt; &gt;     tprnfor = .true. ,<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;  /<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; &amp;SYSTEM<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;                        ibrav = 2,<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;                    celldm(1) = 9.90,<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;                          nat = 2,<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;                         ntyp = 2,<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;             nbnd = 8,<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;                      ecutwfc = 130,<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;         input_dft = &#39;hse&#39;,<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;             nqx1 = 1, nqx2 = 1, nqx3 = 1,<br>
&gt; &gt; &gt;             x_gamma_extrapolation = .TRUE.,<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;  /<br>
&gt; &gt; &gt;  &amp;ELECTRONS<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;                     conv_thr = 1.0d-11 ,<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;                  mixing_beta = 0.7 ,<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;  /<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; ATOMIC_SPECIES<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;    In  114.81400  In.pz-bhs.UPF<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;     N   14.01000  N.pz-vbc.UPF<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; ATOMIC_POSITIONS alat<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;    In      0.000000000    0.000000000    0.000000000<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;     N      0.250000000    0.250000000    0.250000000<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; K_POINTS<br>
&gt; &gt; &gt; 36<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;    0.5 0.5 0.5  1<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; ..<br>
&gt; &gt; &gt; ..<br>
&gt; &gt; &gt; ...<br>
&gt; &gt; &gt; Output, At k= 0 0 0 points only..<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;  k = 0.0000 0.0000 0.0000 (  5985 PWs)   bands (ev):<br>
&gt; &gt; &gt;    -10.0435   2.6536   2.8935   2.8935   2.8935  11.6834  11.6834<br>
&gt; &gt;<br>
&gt; &gt;  11.6834<br>
&gt; &gt;<br>
&gt; &gt; &gt;      highest occupied, lowest unoccupied level (ev):     2.8935<br>
&gt; &gt; &gt;      2.8935 EXX: now go back to refine exchange calculation<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;   -2.4962685539301832<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; highest occupied, lowest unoccupied level (ev):     1.4312    4.4034<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;   -2.4962685539301832       -2.4966288035394730<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; -2.4978367300949738<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;      est. exchange err (dexx)  =       0.00042384 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; !    total energy              =     -24.15032730 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; highest occupied, lowest unoccupied level (ev):     1.4134    4.4023<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;   -2.4978367300949738       -2.4980591338705209<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; -2.4984413942151975<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;      est. exchange err (dexx)  =       0.00007993 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; !    total energy              =     -24.15054868 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;  highest occupied, lowest unoccupied level (ev):     1.4066    4.4806<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;   -2.4995084229633417       -2.4995690173056424<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; -2.4997151620876128<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;      est. exchange err (dexx)  =       0.00004278 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; !    total energy              =     -24.15103886 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; ....<br>
&gt; &gt; &gt; ...<br>
&gt; &gt; &gt; ....<br>
&gt; &gt; &gt; ....<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;      est. exchange err (dexx)  =       0.00000009 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; !    total energy              =     -24.15145297 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;    highest occupied, lowest unoccupied level (ev):     1.4052    4.5809<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;   -2.5006821135594719       -2.5006822803381836<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; -2.5006825016241492<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;      est. exchange err (dexx)  =       0.00000003 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; !    total energy              =     -24.15145336 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;    highest occupied, lowest unoccupied level (ev):     1.4052    4.5809<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;   -2.5006830182691306       -2.5006823182635491<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; -2.5006816409825148<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;      est. exchange err (dexx)  =       0.00000001 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; !    total energy              =     -24.15145346 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;  highest occupied, lowest unoccupied level (ev):     1.4052    4.5809<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;   -2.5006827878685880       -2.5006825627249212<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; -2.5006823476519933<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;      est. exchange err (dexx)  =       0.00000001 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; !    total energy              =     -24.15145351 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;   highest occupied, lowest unoccupied level (ev):     1.4052    4.5809<br>
&gt; &gt; &gt;   -2.5006823476519933       -2.5006827558468370<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; -2.5006831719007772<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;      est. exchange err (dexx)  =       0.00000000 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; !    total energy              =     -24.15145351 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;   highest occupied, lowest unoccupied level (ev):     1.4052    4.5809<br>
&gt; &gt; &gt;   -2.5006831719007772       -2.5006824586209802<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; -2.5006817515742186<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;      est. exchange err (dexx)  =       0.00000000 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; !    total energy              =     -24.15145352 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;   highest occupied, lowest unoccupied level (ev):     1.4052    4.5809<br>
&gt; &gt; &gt;   -2.5006817515742186       -2.5006821655655500<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; -2.5006825845917415<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;      est. exchange err (dexx)  =       0.00000000 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; !    total energy              =     -24.15145353 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;   highest occupied, lowest unoccupied level (ev):     1.4052    4.5809<br>
&gt; &gt; &gt;   -2.5006827994750860       -2.5006825695782151<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; -2.5006823413433099<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;      est. exchange err (dexx)  =       0.00000000 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; !    total energy              =     -24.15145367 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;   highest occupied, lowest unoccupied level (ev):     1.4052    4.5809<br>
&gt; &gt; &gt;   -2.5006823413433099       -2.5006827420502393<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; -2.5006831444189488<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;      est. exchange err (dexx)  =       0.00000000 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; !    total energy              =     -24.15145368 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;   highest occupied, lowest unoccupied level (ev):     1.4052    4.5809<br>
&gt; &gt; &gt;   -2.5006831444189488       -2.5006824325837682<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; -2.5006817224109170<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;      est. exchange err (dexx)  =       0.00000000 Ry<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; !    total energy              =     -24.15145368 Ry<br>
&gt; &gt; &gt; ....<br>
&gt; &gt; &gt; .....<br>
&gt; &gt; &gt; .....<br>
&gt; &gt; &gt; goes on for atleast total time will take approximate 14 hrs but not<br>
&gt; &gt; &gt; stop the run.<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; Please any type of comment is appreciable.<br>
&gt; &gt;<br>
&gt; &gt; --<br>
&gt; &gt; ********************************************************<br>
&gt; &gt; - Article premier - Les hommes naissent et demeurent<br>
&gt; &gt; libres et ègaux en droits. Les distinctions sociales<br>
&gt; &gt; ne peuvent être fondèes que sur l&#39;utilitè commune<br>
&gt; &gt; - Article 2 - Le but de toute association politique<br>
&gt; &gt; est la conservation des droits naturels et<br>
&gt; &gt; imprescriptibles de l&#39;homme. Ces droits sont la libertè,<br>
&gt; &gt; la propriètè, la sùretè et la rèsistance à l&#39;oppression.<br>
&gt; &gt; ********************************************************<br>
&gt; &gt;<br>
&gt; &gt;    Giuseppe Mattioli<br>
&gt; &gt;    CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<br>
&gt; &gt;    v. Salaria Km 29,300 - C.P. 10<br>
&gt; &gt;    I 00015 - Monterotondo Stazione (RM)<br>
&gt; &gt;    Tel + 39 06 90672836 - Fax +39 06 90672316<br>
&gt; &gt;    E-mail: &lt;<a href="mailto:giuseppe.mattioli@ism.cnr.it">giuseppe.mattioli@ism.cnr.it</a>&gt;<br>
<div class="im">&gt; &gt;<br>
&gt; &gt; _______________________________________________<br>
&gt; &gt; Pw_forum mailing list<br>
&gt; &gt; <a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
&gt; &gt; <a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
<br>
</div>--<br>
********************************************************<br>
- Article premier - Les hommes naissent et demeurent<br>
libres et ègaux en droits. Les distinctions sociales<br>
ne peuvent être fondèes que sur l&#39;utilitè commune<br>
- Article 2 - Le but de toute association politique<br>
est la conservation des droits naturels et<br>
imprescriptibles de l&#39;homme. Ces droits sont la libertè,<br>
la propriètè, la sùretè et la rèsistance à l&#39;oppression.<br>
********************************************************<br>
<span class="HOEnZb"><font color="#888888"><br>
   Giuseppe Mattioli<br>
   CNR - ISTITUTO DI STRUTTURA DELLA MATERIA<br>
   v. Salaria Km 29,300 - C.P. 10<br>
   I 00015 - Monterotondo Stazione (RM)<br>
   Tel + 39 06 90672836 - Fax +39 06 90672316<br>
   E-mail: &lt;<a href="mailto:giuseppe.mattioli@ism.cnr.it">giuseppe.mattioli@ism.cnr.it</a>&gt;<br>
</font></span><div class="HOEnZb"><div class="h5">_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://www.democritos.it/mailman/listinfo/pw_forum" target="_blank">http://www.democritos.it/mailman/listinfo/pw_forum</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br>Thanks and Regards<br>Bramha Prasad Pandey<br>Ph.D Student Indian School of Mines(ISM)<br>Dhanbad, INDIA.<br><br>