[Wien] case.inhf
f.tran at pci.uzh.ch
f.tran at pci.uzh.ch
Mon Jan 7 15:22:26 CET 2013
No. The number of occupied (fully or partially) bands is 80. The occupied bands from 1 to 69 are not printed in case.scf.
So, don't use less than 81 for nband (for safety I would not use less than, let's say, 83).
But again, I really doubt that your computational ressources (and patience) are enough for performing hybrid calculations on
your system.
-----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 02:29PM
Subject: Re: [Wien] case.inhf
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
_______________________________________________
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/3e9142d4/attachment.htm>
More information about the Wien
mailing list