[Wien] C1 space group
Laurence Marks
L-marks at northwestern.edu
Fri Mar 14 14:19:42 CET 2014
An alternative approach (not well documented in UG, perhaps not
documented at all) with lstart is to use a number between 0.97 and 1.0
rather than an energy, e.g.
x lstart
SELECT XCPOT:
recommended: 13: PBE-GGA (Perdew-Burke-Ernzerhof 96)
5: LSDA
11: WC-GGA (Wu-Cohen 2006)
19: PBEsol-GGA (Perdew etal. 2008)
19
SELECT ENERGY to separate core and valence states:
recommended: -6.0 Ry (check how much core charge leaks out of MT-sphere)
ALTERNATIVELY: specify charge localization
(between 0.97 and 1.0) to select core state
.99
LSTART ENDS
I use this sometimes, then look in case.outputst, case.in1_st to see
which states have gone into the core.
On Fri, Mar 14, 2014 at 8:07 AM, Jinjan Ren <ren at uni-muenster.de> wrote:
> Dear Peter Blaha:
> Thank you very much for your email! What I have done is first xnn,
> x sgroup, select yes to accept the proposed C-centered monoclinic
> structure; and then x symmetry, copy struct_st, and finally x
> lstart .
> When I set default -6, yes it show the same as you wrote. But when
> I set the Ry value from -10 to -12, it always show the information
> below. there is still CORE electrons leak out of MT-sphere.
> Best,
> Jinjun Ren
>
>
> SELECT XCPOT:
> recommended: 13: PBE-GGA (Perdew-Burke-Ernzerhof 96)
> 5: LSDA
> 11: WC-GGA (Wu-Cohen 2006)
> 19: PBEsol-GGA (Perdew etal. 2008)
> SELECT ENERGY to separate core and valence states:
> recommended: -6.0 Ry (check how much core charge leaks out of MT-sphere)
> ALTERNATIVELY: specify charge localization
> (between 0.97 and 1.0) to select core state
>
> :WARNING: 0.117 P CORE electrons leak out of MT-sphere !!!!
> :WARNING: touch .lcore and run scf-cycle with core density superposition
> :WARNING: Or: rerun lstart with lower E-core separation energy
>
>
>
>
>
> Peter Blaha wrote on 2014-03-14:
>> I hope that after
>
>> x cif2struct
>
>> you run
>
>> setrmt case; copied the new struct file and
>
>> x sgroup
>
>> accepting the proposed C-centered monoclinic structure.
>
>> instgen
>
>> Then try: x lstart
>
>> when you specify -13.5 Ry it does not continue, but repeats the
>> question again, because lstart thinks this is an input error.
>> So you can't run lstart with -13.5 Ry.
>
>> When you run with the default -6 it tells you (among others ...):
>
>> :WARNING: ORBITAL: 3S -10.873 -10.872 for Ga
>> :WARNING: ORBITAL: 3P* -7.433 -7.432
>> :WARNING: ORBITAL: 3P -7.167 -7.165
>> and
>> :WARNING: ORBITAL: 2P* -9.172 -9.141 for P
>> :WARNING: ORBITAL: 2P -9.104 -9.074
>
>> So "natural" energies are either -8.0 (only Ga 3p), -10. (also P 2p)
>> or
>> -11. (also Ga 3s).
>
>> And -10. or -11. does not give any leakage.
>
>
>> On 03/14/2014 10:40 AM, Jinjan Ren wrote:
>> >Dear Wien2k user:
>> > When I set the Ry value from -6 to -13.5 in the wien2k
>> > calculation of Ga(PO3)3, there is always CORE electrons leak out
>> > of MT-sphere. The .in1_st file is below. And the structure file
>> > is enclosed. since in the space group of wien2k, there is no
>> >C1c1(9), I
>> > select Bb(9), and exchange all the values of a and c axis
>> >(including the a,
>> > c and the angle alfa, gamma in cell and the atomic parameters),
>> >. Thanks a
>> > lot!
>> >Best,
>> >Jinjun
>
>
>
>
>
>> > 7.00 10 4 (R-MT*K-MAX; MAX L IN WF, V-NMT
>> > 0.30 5 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 1 0.30 0.000 CONT 1
>> > 1 -6.87 0.001 STOP 1
>> > 2 -1.09 0.002 CONT 1
>> > 2 0.30 0.000 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 0.30 5 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 1 0.30 0.000 CONT 1
>> > 1 -6.87 0.001 STOP 1
>> > 2 -1.09 0.002 CONT 1
>> > 2 0.30 0.000 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 0.30 5 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 1 0.30 0.000 CONT 1
>> > 1 -6.87 0.001 STOP 1
>> > 2 -1.09 0.002 CONT 1
>> > 2 0.30 0.000 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -0.79 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -0.79 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -0.79 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -0.79 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -0.79 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -0.79 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -0.79 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -0.79 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -0.79 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> > 0.30 3 0 (GLOBAL E-PARAMETER WITH n OTHER CHOICES, global
>> >APW/LAPW)
>> > 0 -1.55 0.002 CONT 1
>> > 0 0.30 0.000 CONT 1
>> > 1 0.30 0.000 CONT 1
>> >K-VECTORS FROM UNIT:4 -11.0 1.5 1061 emin / de (emax=Ef+de) / nband
>
>
>
>> >_______________________________________________
>> >Wien mailing list
>> >Wien at zeus.theochem.tuwien.ac.at
>> >http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
>> >SEARCH the MAILING-LIST at:
>> >http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.ht
>> >ml
>
>
>> --
>
>> P.Blaha
>> ---------------------------------------------------------------------
>> -----
>> Peter BLAHA, Inst.f. Materials Chemistry, TU Vienna, A-1060 Vienna
>> Phone: +43-1-58801-165300 FAX: +43-1-58801-165982
>> Email: blaha at theochem.tuwien.ac.at WWW:
>> http://info.tuwien.ac.at/theochem/
>> ---------------------------------------------------------------------
>> -----
>> _______________________________________________
>> Wien mailing list
>> Wien at zeus.theochem.tuwien.ac.at
>> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
>> SEARCH the MAILING-LIST at:
>> http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.htm
>> l
> _______________________________________________
> Wien mailing list
> Wien at zeus.theochem.tuwien.ac.at
> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
> SEARCH the MAILING-LIST at: http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html
--
Professor Laurence Marks
Department of Materials Science and Engineering
Northwestern University
www.numis.northwestern.edu 1-847-491-3996
"Research is to see what everybody else has seen, and to think what
nobody else has thought"
Albert Szent-Gyorgi
More information about the Wien
mailing list