[gpaw-users] problems with Pt atom
Henrik Oberg
henrikob at kth.se
Thu Mar 7 15:32:59 CET 2013
Hi,
thanks a lot for the help. I was unaware of this so thanks for letting
me know.
Best regards
Henrik
On 03 7, 13, at 2:20 PM, Marcin Dulak wrote:
> Hi,
>
> GGA is unable to describe many isolated metal atoms correctly, see:
> http://dx.doi.org/10.1016/S0009-2614(96)01449-2
> https://listserv.fysik.dtu.dk/pipermail/gpaw-users/2011-December/001171.html
> I was never able to converge spin-polarized Pt to the right
> electronic configuration,
> but one should be able to get it spherically-symmetric.
> With width=0.02 eV and a symmetric, cubic box i'm getting for
> example (new 0.9 setup + cg eigensolver):
>
> Band Eigenvalues Occupancy
> 0 -54.69686 2.00000
> 1 -54.69686 2.00000
> 2 -54.69686 2.00000
> 3 -5.14838 1.90501
> 4 -5.14838 1.90501
> 5 -5.14801 1.90330
> 6 -5.14801 1.90330
> 7 -5.14801 1.90330
> 8 -5.06536 0.48007
> 9 -0.57571 0.00000
> 10 -0.57571 0.00000
> 11 -0.57571 0.00000
>
> Best regards,
>
> Marcin
>
> On 03/05/13 14:59, Henrik Oberg wrote:
>> Hi,
>>
>> I need to calculate the total energy of a Pt atom but run into
>> problems. The total energy is fluctuating wildly. Does anyone have
>> any
>> experience with metal atoms?
>> I use something like the following input:
>>
>> calc = GPAW(h=0.2, xc='PBE',
>> txt=name + '.txt',
>> stencils=(3, 3),
>> hund=True,
>> poissonsolver=PoissonSolver(eps=1e-12),
>> mixer=mixersum,
>> )
>>
>>
>> and the energy goes like:
>>
>> ....
>> iter: 139 17:11:07 +1.1 -1.2 -2.233481 0 23
>> +2.0000
>> iter: 140 17:11:28 +1.0 -0.9 -0.800716 0 22
>> +2.0000
>> iter: 141 17:11:51 -0.0 -1.0 -0.590792 0 26
>> +2.0000
>> iter: 142 17:12:11 -0.2 -1.1 -0.848061 0 22
>> +2.0000
>> iter: 143 17:12:35 +1.9 -1.1 -10.409903 0 26
>> +2.0000
>> iter: 144 17:12:57 -0.0 -0.8 0.267276 0 25
>> +2.0000
>> iter: 145 17:13:14 -0.2 -1.3 -0.758502 0 17
>> +2.0000
>> iter: 146 17:13:32 -0.9 -1.2 -0.888202 0 18
>> +2.0000
>> iter: 147 17:13:51 -1.1 -1.1 -0.998467 0 21
>> +2.0000
>> iter: 148 17:14:11 -1.2 -1.1 -0.881081 0 21
>> +2.0000
>> iter: 149 17:14:29 -1.5 -1.1 -0.829883 0 18
>> +2.0000
>>
>>
>>
>> I have tried different eigensolvers, different mixervalues etc
>> according to the covergence issues site...
>>
>> Any suggestions?
>>
>> Thanks in advance!
>>
>> Henrik Öberg
>> _______________________________________________
>> 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
>
> ***********************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://listserv.fysik.dtu.dk/pipermail/gpaw-users/attachments/20130307/0bfdd8db/attachment.html
More information about the gpaw-users
mailing list