[Wien] spin-polarized calculations

foyevtsova foyevtsova at phas.ubc.ca
Tue Dec 11 23:14:00 CET 2018


Thanks, it's working now! case.scfm and other files do look like 
spin-polarized.

But is it normal that there is a small jump in energy:

:ENE  : ********** TOTAL ENERGY IN Ry =      -386334.44392673
:ENE  : ********** TOTAL ENERGY IN Ry =      -386334.44396356
:ENE  : ********** TOTAL ENERGY IN Ry =      -386334.44396524
:ENE  : ********** TOTAL ENERGY IN Ry =      -386334.44847247

and also the :DIS value suddenly becomes quite larger?

:DIS  :  CHARGE DISTANCE       ( 0.0075857 for atom   12 spin 1)      
0.0006768
:DIS  :  CHARGE DISTANCE       ( 0.1267154 for atom   12 spin 1)      
0.0420162

Kateryna



On 2018-12-11 12:08, Peter Blaha wrote:
>> is there a way to use files from a converged non-magnetic calculation 
>> in a new spin-polarized calculation (but with no spin-polarization, of 
>> course)?
>> 
>> I tried the following:
>> 
>> 1) run_lapw
>> 2) cp case.vector case.vectorup(dn)
>>     cp case.energy case.energyup(dn)
>>     cp case.vsp case.vspup(dn)
>>     cp case.vns case.vnsup(dn)
>> 3) x lapw2 -up
>>     x lapw2 -dn
>> 4) x lcore
> 
> You need:  x lcore -up / -dn
> and   rm case.clmsum_old
> 
>> 5) x mixer with
>>     PRATT  0.0   YES  (BROYD/PRATT, BG charge (-1 for core hole), 
>> norm)
>>     1.00            mixing FACTOR for BROYD/PRATT scheme
> 
> check case.scfm. Does it look like a spin-polarized scf file (up and
> dn charges, ....)
> check case.clmup/dn/sum  are they identical /twice as large ???
> 
> And case.clmsum should basically be the same as the save clmsum file
> from run_lapw.
> 
>> 6) clean
>> 7) runsp_lapw
>> 
>> lapw0 finished successfully, but lapw1 gave a
>>   'SELECT' - no energy limits found for atom   6  L= 0
>>   'SELECT' - E-bottom    2.76451   E-top -200.00000
>> error.
>> 
>> Would you have any advice?
>> 
>> Thank you,
>> Kateryna
>> _______________________________________________
>> 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


More information about the Wien mailing list