[gpaw-users] RuO2 NiO

[Marcin Dulak]

Hi,

in some cases setting basis=dzp,
or using smaller number of density mixing steps, like  Mixer(0.05, 2), 
or even Mixer(0.05, 1),
or finally eigensolver='cg' helps.
Please attach the complete example (fine to read from a traj file, but 
preferably a clean python script).

Marcin

On 05/03/12 09:10, Niels Theis Bendtsen wrote:
> Dear GPaw users
>
> I am looking into a rutile RuO2 structure with doped Ni atoms on the surface. However when I add adsorbates (O and OH to a Ru cus site on the surface the calculations fail to converge.
>
> Have you any suggestions for me?
>
> For the clean structure which worked fine I used:
>
> calc = GPAW(h=0.20,
> 	    kpts=(2, 3, 1),
> 	    txt=name + '.txt' ,
> 	    xc='RPBE' ,
> 	    nbands=-50,
> 	    mixer=MixerSum(0.05, 5, weight=100.0),
> 	    poissonsolver=PoissonSolver(eps=1e-12),
> 	    occupations=FermiDirac(0.04),
>              spinpol=True,
> 	    maxiter = 500)
>
> and I set the 3 Ni atoms initial magnetic moment to 1.
>
> Adjusting the mixer down does not work. Increasing the number of bands slightly does not work either.
> Again I would like to state that this configuration worked for a clean surface.
> I hope you can help me.
>
> Niels Theis Bendtsen
>
>
>
>
> _______________________________________________
> gpaw-users mailing list
> gpaw-users at listserv.fysik.dtu.dk
> https://listserv.fysik.dtu.dk/mailman/listinfo/gpaw-users
>


-- 
***********************************

Marcin Dulak
Technical University of Denmark
Department of Physics
Building 307, Room 229
DK-2800 Kongens Lyngby
Denmark
Tel.: (+45) 4525 3157
Fax.: (+45) 4593 2399
email: Marcin.Dulak at fysik.dtu.dk

***********************************