[Wien] case.inhf

[ali ghafari]

Dear Prof. Tran

Thank you very much for your replay. I'm using cluster for this calculation.
the :BAN in case.scf is in the following, as you said the last column is occupancy but it is less than 11?is it right?
with Best Regards
Ali

        Bandranges (emin - emax) and occupancy:
:BAN00070:  70    0.428082    0.509292  1.00000000
:BAN00071:  71    0.440837    0.511967  1.00000000
:BAN00072:  72    0.460794    0.520858  1.00000000
:BAN00073:  73    0.471917    0.522531  1.00000000
:BAN00074:  74    0.471972    0.529510  1.00000000
:BAN00075:  75    0.495224    0.541786  1.00000000
:BAN00076:  76    0.496657    0.544179  1.00000000
:BAN00077:  77    0.509514    0.760574  0.51995576
:BAN00078:  78    0.509898    0.765628  0.40435460
:BAN00079:  79    0.604807    0.829075  0.03620490
:BAN00080:  80    0.615796    0.867992  0.01107466
:BAN00081:  81    0.801158    0.905668  0.00000000
:BAN00082:  82    0.802653    0.905668  0.00000000
:BAN00083:  83    0.846712    0.913216  0.00000000
:BAN00084:  84    0.866694    0.925919  0.00000000
:BAN00085:  85    0.920677    1.009571  0.00000000  




________________________________
 From: "f.tran at pci.uzh.ch" <f.tran at pci.uzh.ch>
To: A Mailing list for WIEN2k users <wien at zeus.theochem.tuwien.ac.at> 
Sent: Monday, January 7, 2013 1:29 PM
Subject: Re: [Wien] case.inhf
 

11 can not be the number of occupied bands. Search for :BAN in case.scf
and then you can see how many bands are occupied (the last column is the occupancy).

Furthermore, your system (Bi2Sr2CaCu2O8) looks maybe too big for hybrid calculations.
I don't know what computer ressources you have, but I doubt that you can do a hybrid calculations
on this system (hybrid calculations are VERY expensive).

F. Tran

-----wien-bounces at zeus.theochem.tuwien.ac.at wrote: -----
To: A Mailing list for WIEN2k users <wien at zeus.theochem.tuwien.ac.at>
From: ali ghafari 
Sent by: wien-bounces at zeus.theochem.tuwien.ac.at
Date: 01/07/2013 12:51PM
Subject: Re: [Wien] case.inhf


Dear Prof. Tran
regarding your suggestion about ' nband' in case.inhf.
In the Bi2Sr2CaCu2O8 compound I would like to use B3LYP, while the occupation number is about 11 at the case.scf.
The question is that when I put any numbers less than 100 the error appears due to the less value of 'nband' and calculation is stoped.   
What is your advice?
 Wwith Best Regards
Ali




________________________________
 From: "tran at theochem.tuwien.ac.at" <tran at theochem.tuwien.ac.at>
To: A Mailing list for WIEN2k users <wien at zeus.theochem.tuwien.ac.at> 
Sent: Thursday, September 20, 2012 1:40 PM
Subject: Re: [Wien] case.inhf
 
Search for :BAN in case.scf. The last column is the occupation number.

On Thu, 20 Sep 2012, ali ghafari wrote:

> Dear Prof. Tran
> Thank you very much for your replay.
> But  Could you please explain to me how can I find the number of partially occupied band in case.scf?
> I can not find any explanation in the UG.
> 
> Best Regards
> Ali
> 
> 
> 
> ________________________________
>  From: "tran at theochem.tuwien.ac.at" <tran at theochem.tuwien.ac.at>
> To: A Mailing list for WIEN2k users <wien at zeus.theochem.tuwien.ac.at> 
> Sent: Wednesday, September 19, 2012
 1:06 PM
> Subject: Re: [Wien] case.inhf
>  
> Why not, but I think that it is not really more simple than just looking
> at the number n_occ of (partially) occupied bands in case.scf and choosing
> nband = n_occ + a few more bands. But remember that nband in case.inhf is
> a parameter (similar to R*K_max) which has to be tested.
> 
> One last thing: the purpose of nband in case.in1 is not at all the same
> as nband in case.inhf.
> 
> F. Tran
> 
> On Wed, 19 Sep 2012, ali ghafari wrote:
> 
> > Dear Prof. Blaha
> > 
> > For hybrid functionals case.inhf is necessary as discussed in the UG on pages 49 and 97. But the value of "nband' at case.inhf is confusing. on the page 97 of UG has mentioned "  .... nband should be at least equal to the number of (partially) occupied bands plus one". While at the line 5 of case.in1 the value of nband
 is automatically determined ( UG page 102) by "   nband = ne ∗ 2.0 + 5 " which is number of eigenvalues. Can we consider inband of case.inhf is equals: (eigenvalues -5)/2+1 which is extracted from case.in1?
> >    
> > 
> > Best Regards
> > Ali
> _______________________________________________
> Wien mailing list
> Wien at zeus.theochem.tuwien.ac.at
> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
_______________________________________________
Wien mailing list
Wien at zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien



_______________________________________________
Wien mailing list
Wien at zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
 
_______________________________________________
Wien mailing list
Wien at zeus.theochem.tuwien.ac.at
http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20130107/036b78fa/attachment.htm>