[Wien] Error: 'LOPW plane waves exhausted'
Laurence Marks
L-marks at northwestern.edu
Thu Jul 28 16:18:29 CEST 2011
This is tricky, particularly as something in your email triggered
gmail's spam filter.
More information than you want to know. A change was made in Wien2k 11
was made to improve the orthogonality of the orbitals, to get rid of
what I call "ghost forces" that showed up with MSR1a. These are
sudden, anomalously large valence forces that occurred. When this
change was made Peter was concerned that this might lead to a "LOPW
plane wave exhausted" problem that you are seeing for very high
symmetry structures.
Probably you can solve this by just increasing RKMAX to 8 or 9 from
the default of 7. This increases the number of plane waves so should
cure the problem.
On Tue, Jul 26, 2011 at 4:49 AM, pascal boulet
<pascal.boulet at univ-provence.fr> wrote:
>
> -----BEGIN PGP SIGNED MESSAGE-----
>
>
> Dear all,
>
> I am trying to set up calculations on a Gamma brass structure. The
> problem is that W2K (version 11) stops with the error 'LOPW plane waves
> exhausted' in lapw1. The topic has already been reported on the list
> (Feb. 2011) but the answers there are not very clear to me. Here is what
> I am doing.
>
> So, I am starting with Copper with a gamma brass structure (group number
> 217). The structure generated by W2K from a CIF file (4 inequivalent
> positions) looks correct. The local symmetry of the atomic sites agrees
> with that given on the Bilbao website.
>
> Now comes the generation of the atomic densities. According to the
> electronic structure of Cu which is:
>
> E-up(Ry) E-dn(Ry) Occupancy q/sphere core-state
> 1S -650.760163 -650.760146 1.00 1.00 1.0000 T
> 2S -78.505196 -78.461906 1.00 1.00 1.0000 T
> 2P* -68.813819 -68.782078 1.00 1.00 1.0000 T
> 2P -67.308434 -67.275562 2.00 2.00 1.0000 T
> 3S -8.856436 -8.771764 1.00 1.00 0.9999 F
> 3P* -5.848176 -5.764727 1.00 1.00 0.9996 F
> 3P -5.655968 -5.573289 2.00 2.00 0.9995 F
> 3D* -0.747458 -0.672780 2.00 2.00 0.9718 F
> 3D -0.725704 -0.651398 3.00 2.00 0.9708 F
> 4S -0.433190 -0.420827 1.00 1.00 0.4335 F
>
> I choose the core-valence energy separation to be -10 Ry because the
> default -6.0 Ry would separate the 3s and the 3p states which does not
> seem wise to me. I think this is what is done in the example for Zn in
> the manual. Is it correct to choose -10 Ry?
>
> Here is a sample of the file case.in1c obtained automatically using w2web:
>
> WFFIL EF= 0.50000 (WFFIL, WFPRI, ENFIL, SUPWF)
> 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)
> 0 0.30 0.000 CONT 1
> 0 -8.56 0.001 STOP 1
> 1 0.30 0.000 CONT 1
> 1 -5.36 0.001 STOP 1
> 2 0.30 0.005 CONT 1
> etc. for the other 3 atoms.
>
> which I changed to:
>
> WFFIL EF= 0.50000 (WFFIL, WFPRI, ENFIL, SUPWF)
> 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)
> 0 0.30 0.000 CONT 0
> 0 -8.56 0.001 STOP 0
> 1 0.30 0.000 CONT 0
> 1 -5.36 0.001 STOP 0
> 2 0.30 0.005 CONT 1
>
> in order to match the example for TiO2. Whether I change the file or not
> does not change the final result anyway, I mean the crash!
>
> In the mailing-list it is suggested to increase RKmax, which I did up to
> 9. No change.
>
> For information here is the cif file of the structure:
> data_generated_by_bilbao_crystallographic_server
>
> _cell_length_a 8.878
> _cell_length_b 8.878
> _cell_length_c 8.878
> _cell_angle_alpha 90
> _cell_angle_beta 90
> _cell_angle_gamma 90
> _symmetry_space_group_name_H-M 'I-43m'
> _symmetry_Int_Tables_number 217
>
> loop_
> _symmetry_equiv_pos_site_id
> _symmetry_equiv_pos_as_xyz
> 1 'x,y,z'
> 2 '-x,-y,z'
> 3 '-x,y,-z'
> 4 'x,-y,-z'
> 5 'z,x,y'
> 6 'z,-x,-y'
> 7 '-z,-x,y'
> 8 '-z,x,-y'
> 9 'y,z,x'
> 10 '-y,z,-x'
> 11 'y,-z,-x'
> 12 '-y,-z,x'
> 13 'y,x,z'
> 14 '-y,-x,z'
> 15 'y,-x,-z'
> 16 '-y,x,-z'
> 17 'x,z,y'
> 18 '-x,z,-y'
> 19 '-x,-z,y'
> 20 'x,-z,-y'
> 21 'z,y,x'
> 22 'z,-y,-x'
> 23 '-z,y,-x'
> 24 '-z,-y,x'
> 25 'x+1/2,y+1/2,z+1/2'
> 26 '-x+1/2,-y+1/2,z+1/2'
> 27 '-x+1/2,y+1/2,-z+1/2'
> 28 'x+1/2,-y+1/2,-z+1/2'
> 29 'z+1/2,x+1/2,y+1/2'
> 30 'z+1/2,-x+1/2,-y+1/2'
> 31 '-z+1/2,-x+1/2,y+1/2'
> 32 '-z+1/2,x+1/2,-y+1/2'
> 33 'y+1/2,z+1/2,x+1/2'
> 34 '-y+1/2,z+1/2,-x+1/2'
> 35 'y+1/2,-z+1/2,-x+1/2'
> 36 '-y+1/2,-z+1/2,x+1/2'
> 37 'y+1/2,x+1/2,z+1/2'
> 38 '-y+1/2,-x+1/2,z+1/2'
> 39 'y+1/2,-x+1/2,-z+1/2'
> 40 '-y+1/2,x+1/2,-z+1/2'
> 41 'x+1/2,z+1/2,y+1/2'
> 42 '-x+1/2,z+1/2,-y+1/2'
> 43 '-x+1/2,-z+1/2,y+1/2'
> 44 'x+1/2,-z+1/2,-y+1/2'
> 45 'z+1/2,y+1/2,x+1/2'
> 46 'z+1/2,-y+1/2,-x+1/2'
> 47 '-z+1/2,y+1/2,-x+1/2'
> 48 '-z+1/2,-y+1/2,x+1/2'
>
> loop_
> _atom_site_label
> _atom_site_type_symbol
> _atom_site_symmetry_multiplicity
> _atom_site_Wyckoff_symbol
> _atom_site_fract_x
> _atom_site_fract_y
> _atom_site_fract_z
> Cu1 Cu ? 8c 0.10890 0.10890 0.10890
> Cu2 Cu ? 8c 0.82800 0.82800 0.82800
> Cu3 Cu ? 12e 0.35580 0.00000 0.00000
> Cu4 Cu ? 24g 0.31280 0.31280 0.03660
>
> I can provide additional information from W2K output files if necessary.
>
> Thank you for your help
> Best regards
> Pascal
>
> * Anglais - détecté
> * Anglais
> * Français
>
> * Anglais
> * Français
>
> <javascript:void(0);>
> -----BEGIN PGP SIGNATURE-----
> Version: GnuPG v1.4.10 (GNU/Linux)
> Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org/
>
> iEYEARECAAYFAk4ujYwACgkQNLjdBN1V75kChQCfc+d/KY5CX3RuFj44LKCZJOA9
> sFAAnidVQa6Pih00ee9IJFVyHlTr8fEM
> =EuJr
> -----END PGP SIGNATURE-----
>
> _______________________________________________
> Wien mailing list
> Wien at zeus.theochem.tuwien.ac.at
> http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien
>
--
Laurence Marks
Department of Materials Science and Engineering
MSE Rm 2036 Cook Hall
2220 N Campus Drive
Northwestern University
Evanston, IL 60208, USA
Tel: (847) 491-3996 Fax: (847) 491-7820
email: L-marks at northwestern dot edu
Web: www.numis.northwestern.edu
Chair, Commission on Electron Crystallography of IUCR
www.numis.northwestern.edu/
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