[Wien] Promoting Ge 3s shell to valence

Priyanka Seth priyanka.seth at polytechnique.edu
Thu Nov 19 11:13:13 CET 2015


Dear all,

I am trying to do a calculation on a Ge compound where I would like to 
promote the core 3s states of Ge (configuration [Ar] 3s2 3p2) to 
semicore and hence be included in the valence states. I am hoping this 
will have some convergence problems I had in the scf cycles.

However, I am not able to select an energy cutoff during lstart that 
allows me to do this:

  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
-12.5
   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
-12.0
STOP LSTART ENDS
0.266u 0.004s 11:47.23 0.0%    0+0k 0+1384io 0pf+0w

where the energy of the states are as follows:

           E-up(Ry)      E-dn(Ry)   Occupancy   q/sphere  core-state
   1S    -803.903411   -803.898620  1.00  1.00    1.0000  T
   1S    -803.903411   -803.898620  1.00  1.00    1.0000  T
   2S    -100.595366   -100.593339  1.00  1.00    1.0000  T
   2S    -100.595366   -100.593339  1.00  1.00    1.0000  T
   2P*    -89.575687    -89.572019  1.00  1.00    1.0000  T
   2P*    -89.575687    -89.572019  1.00  1.00    1.0000  T
   2P     -87.237957    -87.234392  2.00  2.00    1.0000  T
   2P     -87.237957    -87.234392  2.00  2.00    1.0000  T
   3S     -12.310095    -12.304838  1.00  1.00    1.0000  T
   3S     -12.310095    -12.304838  1.00  1.00    1.0000  T
   3P*     -8.676791     -8.670481  1.00  1.00    0.9999  F
   3P*     -8.676791     -8.670481  1.00  1.00    0.9999  F
   3P      -8.361586     -8.355077  2.00  2.00    0.9999  F
   3P      -8.361586     -8.355077  2.00  2.00    0.9999  F
   3D*     -2.187019     -2.174623  2.00  2.00    0.9944  F
   3D*     -2.187019     -2.174623  2.00  2.00    0.9944  F
   3D      -2.144032     -2.131383  3.00  3.00    0.9941  F
   3D      -2.144032     -2.131383  3.00  3.00    0.9941  F
   4S      -0.907245     -0.814618  1.00  1.00    0.6202  F
   4S      -0.907245     -0.814618  1.00  1.00    0.6202  F
   4P*     -0.336907     -0.246809  1.00  0.00    0.3584  F
   4P*     -0.336907     -0.246809  1.00  0.00    0.3584  F
   4P      -0.322467     -0.233001  1.00  0.00    0.3414  F
   4P      -0.322467     -0.233001  1.00  0.00    0.3414  F

I found an earlier question on the mailing list where the user wished to 
use a very low cutoff (~ -30 Ry) here: 
http://www.mail-archive.com/wien%40zeus.theochem.tuwien.ac.at/msg05396.html. 
What I am trying to do is far less extreme, I am only promoting one 
shell to the valence. Why can I not do this in lstart, and is there 
another way around this problem?

Thank you in advance for your help.

Best wishes,
Priyanka Seth


More information about the Wien mailing list