DFT+U periodicity (very confused results)

Problems running VASP: crashes, internal errors, "wrong" results.


Moderators: Global Moderator, Moderator

Post Reply
Message
Author
jbliu
Newbie
Newbie
Posts: 13
Joined: Tue May 06, 2014 12:00 am
License Nr.: 5-1125

DFT+U periodicity (very confused results)

#1 Post by jbliu » Thu Mar 31, 2016 8:51 am

Dear Adimin,

The system of my present VASP job is Li-intercalated β-MnO2 grain boundary. A collinear antiferromagnetic ordering is used in all our calculations. Valence electron configurations of the potentials were taken as 1s2 2s1 for Li, 3d6 4s1 for Mn, and 2s2 2p4 for O. But I find some very confused results.

I give “Structure A” a static calculation, the resulted energy is “-1617.4093eV”. Then I put “Structure A” just translate 0.5 along direction “a” to get “Structure B”. But the resulted energy is “-1614.0948” after a static calculation! To my knowledge, “Structure A” and “Structure B” are just the same due to Periodic boundary conditions, isn’t it?

Here are my relaxation inputs:

INCAR:

SYSTEM = MnO2
LWAVE=.FALSE.
LCHARG=.FALSE.
LASPH=.TRUE.
LREAL= Auto

ISTART=0
ICHARG=2
ALGO=Normal
PREC=Normal

#LPLANE = .TRUE.
NPAR = 6
#LSCALU = .FALSE.
#NSIM = 4

ISPIN=2
ENCUT=500
ISMEAR=0
SIGMA=0.2

LDAU=.TRUE.
LDAUTYPE=2
LDAUL= -1 2 -1
LDAUU= 0.0 6.1 0.0
LDAUJ= 0.0 1.0 0.0

LORBIT=10
NELM=300
MAGMOM=160*0 40*-1 40*1 1*0

KPOINTS:

AUTO GRID
0
M
2 1 1
0 0 0

The POSCAR file for "Structure A":


New structure
1.00000000000000
5.8893198967000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 14.0908565521000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 28.1817131042000000
O Mn Li
160 80 1
Direct
0.2499999850000023 0.2604309137599969 0.9993153614700034
0.7500000000000000 0.2604309137599969 0.9993153614700034
0.0000000000000000 0.7604308247600002 0.0006845765299985
0.5000000000000000 0.7604308247600002 0.0006845765299985
0.2499999850000023 0.9535387172199989 0.0317670078599974
0.7500000000000000 0.9535387172199989 0.0317670078599974
0.0000000000000000 0.5744181912600013 0.0326922833200030
0.5000000000000000 0.5744181912600013 0.0326922833200030
0.0000000000000000 0.1181449833100032 0.0408190930699988
0.5000000000000000 0.1181449833100032 0.0408190930699988
0.2499999850000023 0.4300077190800025 0.0583602348499994
0.7500000000000000 0.4300077190800025 0.0583602348499994
0.0000000000000000 0.2750738002499986 0.0769547092800025
0.5000000000000000 0.2750738002499986 0.0769547092800025
0.0000000000000000 0.8281392854100034 0.0824713827299988
0.5000000000000000 0.8281392854100034 0.0824713827299988
0.2499999850000023 0.6702044481399980 0.0958871600100011
0.7500000000000000 0.6702044481399980 0.0958871600100011
0.2499999850000023 0.1405752276699985 0.1174443782200001
0.7500000000000000 0.1405752276699985 0.1174443782200001
0.0000000000000000 0.9835755594300011 0.1176660121600008
0.5000000000000000 0.9835755594300011 0.1176660121600008
0.0000000000000000 0.5253804691500008 0.1296273580500014
0.5000000000000000 0.5253804691500008 0.1296273580500014
0.2499999850000023 0.3674343673899969 0.1408288305599967
0.7500000000000000 0.3674343673899969 0.1408288305599967
0.2499999850000023 0.8420608161950947 0.1606329583360910
0.7500000000000000 0.8419576349941530 0.1605784805773056
0.0002971055496690 0.6755195246260213 0.1731161805975153
0.4997028944503308 0.6755195246260213 0.1731161805975153
0.0006505469152588 0.2228384026803982 0.1827273190951359
0.4993494530847414 0.2228384026803982 0.1827273190951359
0.2499999850000023 0.0567571962727391 0.1897838676862794
0.7500000000000000 0.0575895704817379 0.1900331236710447
0.2499999850000023 0.5333032025208106 0.2109229664704642
0.7500000000000000 0.5345080075642623 0.2112907492968576
-0.0006375979311096 0.3755771804143354 0.2215603738685823
0.5006375979311104 0.3755771804143354 0.2215603738685823
-0.0006449500638147 0.9210834796620457 0.2307882159450288
0.5006449500638148 0.9210834796620457 0.2307882159450288
0.2499999850000023 0.7637251282672834 0.2385963064728941
0.7500000000000000 0.7630267008509802 0.2384581761160290
0.2499999850000023 0.2404468496006257 0.2623529424465632
0.7500000000000000 0.2362062497972112 0.2612751170217582
0.0024439850343002 0.0767166995824597 0.2668567154721664
0.4975560149656996 0.0767166995824597 0.2668567154721664
0.0008614027313266 0.6212415261444986 0.2788826594566275
0.4991385972686734 0.6212415261444986 0.2788826594566275
0.2499999850000023 0.4597567637593273 0.2886836727744723
0.7500000000000000 0.4620190046964892 0.2890293108191440
0.2499999850000023 0.9420308867353459 0.3098522501200190
0.7500000000000000 0.9420132621836951 0.3098854168004734
-0.0008512404902403 0.7757784128390164 0.3169074877027884
0.5008512404902401 0.7757784128390164 0.3169074877027884
-0.0120007479172201 0.3173569750939334 0.3325693485316588
0.5120007479172195 0.3173569750939334 0.3325693485316588
0.2499999850000023 0.1541879345581658 0.3379368137490409
0.7500000000000000 0.1437268142562705 0.3354093128287287
0.2499999850000023 0.6326832332071288 0.3595766382397638
0.7500000000000000 0.6315279577694289 0.3594197522689878
0.0015371559408095 0.4751923229431999 0.3703885155218432
0.4984628440591906 0.4751923229431999 0.3703885155218432
0.0045195278578886 0.0135521988048136 0.3843985854734830
0.4954804721421112 0.0135521988048136 0.3843985854734830
0.2499999850000023 0.8556894410653192 0.3817916758165630
0.7500000000000000 0.8595182822188527 0.3827875333635560
0.2499999850000023 0.3309474363045138 0.4092602568742005
0.7500000000000000 0.3427162216387513 0.4113952013895401
-0.0107311371496977 0.1770590615156627 0.4146196468539540
0.5107311371496979 0.1770590615156627 0.4146196468539540
-0.0007235647463764 0.7253265118554928 0.4231794258398183
0.5007235647463767 0.7253265118554928 0.4231794258398183
0.2499999850000023 0.5698794736869609 0.4423095165252708
0.7500000000000000 0.5688584430483821 0.4419577574804776
0.0005596655150598 0.8821136025423965 0.4594080168732146
0.4994403344849405 0.8821136025423965 0.4594080168732146
0.0029925330184703 0.4252692603299986 0.4700633316802076
0.4970074669815297 0.4252692603299986 0.4700633316802076
0.2499999850000023 0.0485745071718454 0.4684753548873323
0.7500000000000000 0.0470573140457040 0.4691513219628882
0.0001207273837843 0.2402746369224640 0.4991804763789167
0.4998792726162157 0.2402746369224640 0.4991804763789167
0.2499999850000023 0.7389487129825475 0.5008825702703713
0.7500000000000000 0.7388168342658420 0.5010282686150459
-0.0009409898642612 0.5486889040749999 0.5327471629372299
0.5009409898642612 0.5486889040749999 0.5327471629372299
0.2499999549999998 0.9257086960048311 0.5328647004274553
0.7499999400000021 0.9256884096239244 0.5328491795872844
0.2499999549999998 0.3830304341540672 0.5425025370865157
0.7499999400000021 0.3828153074656781 0.5421271456384813
0.0002450623039668 0.0700565214205394 0.5583883085417956
0.4997549376960332 0.0700565214205394 0.5583883085417956
0.2499999549999998 0.2244516310361494 0.5759660206034765
0.7499999400000021 0.2244407893671252 0.5759663325402654
0.2499999549999998 0.6714880917998274 0.5831662427237271
0.7499999400000021 0.6718877370769896 0.5835194765452297
0.0001168681418963 0.8308067569499871 0.5960971566501523
0.4998831318581032 0.8308067569499871 0.5960971566501523
-0.0001422219101575 0.3569873034009181 0.6182750599317623
0.5001422219101574 0.3569873034009181 0.6182750599317623
0.2499999549999998 0.5158306852305059 0.6180236575079782
0.7499999400000021 0.5159832717580731 0.6181816550317707
0.2499999549999998 0.9753430820891580 0.6297425345350216
0.7499999400000021 0.9752966899954144 0.6297433605799252
-0.0000111761964053 0.1333632691188763 0.6408525622000690
0.5000111761964056 0.1333632691188763 0.6408525622000690
0.0000894396663909 0.6579173337573946 0.6609025733280827
0.4999105603336093 0.6579173337573946 0.6609025733280827
0.2499999549999998 0.8250026629960682 0.6734695584508358
0.7499999400000021 0.8249755558566991 0.6734698762337797
0.2499999549999998 0.2768652009806421 0.6832008056017919
0.7499999400000021 0.2768658582137478 0.6832012083290943
0.0001120208297647 0.4419460398165975 0.6902560700044937
0.4998879791702355 0.4419460398165975 0.6902560700044937
0.0000038733377261 0.9670185125237015 0.7106870094138881
0.4999961266622742 0.9670185125237015 0.7106870094138881
0.2499999549999998 0.1242758159289821 0.7221543301972995
0.7499999400000021 0.1242593512482375 0.7221545005651921
0.2499999549999998 0.5771713730952310 0.7302691924461743
0.7499999400000021 0.5770740951510157 0.7302317785191008
-0.0000624819751946 0.7355601111227499 0.7387120463868266
0.5000624819751945 0.7355601111227499 0.7387120463868266
-0.0000009940682968 0.2651459138216287 0.7617478053216844
0.5000009940682969 0.2651459138216287 0.7617478053216844
0.2499999549999998 0.4229194841275341 0.7694491431107383
0.7499999400000021 0.4229172004771223 0.7694498664645296
0.2499999549999998 0.8767669274389217 0.7780156426712092
0.7499999400000021 0.8767860982033757 0.7780245411562264
-0.0000079398225850 0.0340563233738224 0.7893162434504482
0.5000079398225852 0.0340563233738224 0.7893162434504482
0.0000074429937248 0.5593209527705383 0.8095098625188276
0.4999925570062759 0.5593209527705383 0.8095098625188276
0.2499999549999998 0.7238728942349867 0.8172397893727178
0.7499999400000021 0.7238711063626166 0.8172419039135788
0.2499999549999998 0.1758593527734923 0.8268588152112514
0.7499999400000021 0.1758400622275137 0.8268375827300717
0.0000018587375673 0.3431600858225707 0.8396203642696094
0.4999981412624322 0.3431600858225707 0.8396203642696094
0.0000000000000000 0.8674342783900002 0.8591710054399968
0.5000000000000000 0.8674342783900002 0.8591710054399968
0.2499999549999998 0.0253805131499973 0.8703725379499971
0.7499999400000021 0.0253805131499973 0.8703725379499971
0.2499999549999998 0.4835755594300011 0.8823338318400005
0.7499999400000021 0.4835755594300011 0.8823338318400005
0.0000000000000000 0.6405751236699970 0.8825555247800025
0.5000000000000000 0.6405751236699970 0.8825555247800025
0.0000000000000000 0.1702044921400017 0.9041126909899972
0.5000000000000000 0.1702044921400017 0.9041126909899972
0.2499999549999998 0.3281393154099987 0.9175284682699996
0.7499999400000021 0.3281393154099987 0.9175284682699996
0.2499999549999998 0.7750737102499983 0.9230451497199965
0.7499999400000021 0.7750737102499983 0.9230451497199965
0.0000000000000000 0.9300076300799986 0.9416396241499996
0.5000000000000000 0.9300076300799986 0.9416396241499996
0.2499999549999998 0.6181448713100011 0.9591807759300011
0.7499999400000021 0.6181448713100011 0.9591807759300011
0.2499999549999998 0.0744182212599966 0.9673075916799974
0.7499999400000021 0.0744182212599966 0.9673075916799974
0.0000000000000000 0.4535387172199989 0.9682328621399989
0.5000000000000000 0.4535387172199989 0.9682328621399989
0.2499999549999998 0.5138390078300006 0.0015253939999980
0.7499999400000021 0.5138390078300006 0.0015253939999980
0.2499999850000023 0.1993901048399991 0.0566408064400008
0.7500000000000000 0.1993901048399991 0.0566408064400008
0.2499999850000023 0.9060885862599974 0.0988101646000032
0.7500000000000000 0.9060885862599974 0.0988101646000032
0.2499999850000023 0.5999232770899994 0.1520907401000002
0.7500000000000000 0.5999232770899994 0.1520907401000002
0.2499999850000023 0.3008428212059368 0.2004746674020850
0.7500000000000000 0.3004686864583284 0.2004145351132569
0.2499999850000023 -0.0022830969378047 0.2491242479237537
0.7500000000000000 -0.0016729204585980 0.2494171498791961
0.2499999850000023 0.6980877640575661 0.2987960628345260
0.7500000000000000 0.6978062403791558 0.2985975161943797
0.2499999850000023 0.3965984679588764 0.3493347498196309
0.7500000000000000 0.4068340477335063 0.3480838273628763
0.2499999850000023 0.0957195721128111 0.4007939683158623
0.7500000000000000 0.0857217699034751 0.4037296176301125
0.2499999850000023 0.8003996363836081 0.4438746991202605
0.7500000000000000 0.8005344734116436 0.4441029618302292
0.2499999549999998 0.4881036069631003 0.5003589753267451
0.7499999400000021 0.4876745192723155 0.4999236507442200
0.2499999549999998 0.7974985192511311 0.5572144187986685
0.7499999400000021 0.7978723646315707 0.5572557440732899
0.2499999549999998 0.0983776311690586 0.6006267809908650
0.7499999400000021 0.0983386391650913 0.6005923660910129
0.2499999549999998 0.3997851465675423 0.6539410578201373
0.7499999400000021 0.3997465986772131 0.6539951995682143
0.2499999549999998 0.6984042008077207 0.6996648995318631
0.7499999400000021 0.6983876302623856 0.6996499176211995
0.2499999549999998 0.0002279144413511 0.7500903477642973
0.7499999400000021 0.0002415967676132 0.7500885161408217
0.2499999549999998 0.3022636030825932 0.8009843602691160
0.7499999400000021 0.3022617839027095 0.8009846982522415
0.2499999549999998 0.5998386531500017 0.8466550186300026
0.7499999400000021 0.5998386531500017 0.8466550186300026
0.2499999549999998 0.9014906408599970 0.8997363362200019
0.7499999400000021 0.9014906408599970 0.8997363362200019
0.2499999549999998 0.2017832238700024 0.9432545935800007
0.7499999400000021 0.2017832238700024 0.9432545935800007
0.0000000000000000 0.0138390678299984 0.9984746060000020
0.5000000000000000 0.0138390678299984 0.9984746060000020
0.0000000000000000 0.7017831788700022 0.0567452724200024
0.5000000000000000 0.7017831788700022 0.0567452724200024
0.0000000000000000 0.4014907298600008 0.1002635147800035
0.5000000000000000 0.4014907298600008 0.1002635147800035
0.0000000000000000 0.0998387641500003 0.1533448773700030
0.5000000000000000 0.0998387641500003 0.1533448773700030
-0.0000669106130973 0.8010524584715649 0.1992452563700227
0.5000669106130974 0.8010524584715649 0.1992452563700227
0.0013891444162654 0.4984111888015441 0.2500021820896446
0.4986108555837344 0.4984111888015441 0.2500021820896446
0.0010167997681741 0.1942262581342066 0.2984093250020118
0.4989832002318258 0.1942262581342066 0.2984093250020118
0.0005237018682229 0.8981682083708146 0.3464282580393911
0.4994762981317774 0.8981682083708146 0.3464282580393911
0.0005503219952602 0.5998279662360126 0.4005252798476874
0.4994496780047398 0.5998279662360126 0.4005252798476874
0.0012909187222209 0.2999347938791176 0.4462512577860125
0.4987090812777793 0.2999347938791176 0.4462512577860125
0.0014831492999859 0.9857843675029612 0.5022166010267808
0.4985168507000141 0.9857843675029612 0.5022166010267808
0.0000598087550793 0.3018009119267794 0.5579178611649239
0.4999401912449204 0.3018009119267794 0.5579178611649239
-0.0001990915651206 0.5942209115296437 0.5999597389852391
0.5001990915651204 0.5942209115296437 0.5999597389852391
0.0000076878842180 0.9005072562125521 0.6522621768430167
0.4999923121157820 0.9005072562125521 0.6522621768430167
0.0000132958075497 0.1995361821516842 0.7016829122572804
0.4999867041924503 0.1995361821516842 0.7016829122572804
0.0000184649159010 0.5003509789083100 0.7502682255206382
0.4999815350840990 0.5003509789083100 0.7502682255206382
-0.0000113782044715 0.8011423397644464 0.7989555324884546
0.5000113782044709 0.8011423397644464 0.7989555324884546
0.0000000000000000 0.0999233000899977 0.8479091558999983
0.5000000000000000 0.0999233000899977 0.8479091558999983
0.0000000000000000 0.4060886152600034 0.9011896864000022
0.5000000000000000 0.4060886152600034 0.9011896864000022
0.0000000000000000 0.6993900008399976 0.9433590595600023
0.5000000000000000 0.6993900008399976 0.9433590595600023
0.7500000000000000 0.2446480898961925 0.3722180202783857


Thank you very much!

Jian Xu

admin
Administrator
Administrator
Posts: 2921
Joined: Tue Aug 03, 2004 8:18 am
License Nr.: 458

Re: DFT+U periodicity (very confused results)

#2 Post by admin » Thu Mar 31, 2016 11:35 am

Shifted and not shifted cells represent the same structure.
Because of large cells the numerical noise can cause
that calculated energies are not exactly the same, but
estimated deviations are of the order of 1meV.
Check your cells, there must be some (non-negligible) difference.

jbliu
Newbie
Newbie
Posts: 13
Joined: Tue May 06, 2014 12:00 am
License Nr.: 5-1125

Re: DFT+U periodicity (very confused results)

#3 Post by jbliu » Thu Mar 31, 2016 4:48 pm

admin wrote:Shifted and not shifted cells represent the same structure.
Because of large cells the numerical noise can cause
that calculated energies are not exactly the same, but
estimated deviations are of the order of 1meV.
Check your cells, there must be some (non-negligible) difference.
Thank you for your reply. Would you mind check those two strucures for static calculation, respectively? I would appreciate it very much!

The POSCAR file for "Structure B" is presented here. The INCAR, KPOINTS, and POTCAR (PAW-PBE)files are just the same with the static calculation for "Structure A".Valence electron configurations of the potentials were taken as 1s2 2s1 for Li, 3d6 4s1 for Mn, and 2s2 2p4 for O.

New structure
1.00000000000000
5.8893198967000000 0.0000000000000000 0.0000000000000000
0.0000000000000000 14.0908565521000000 0.0000000000000000
0.0000000000000000 0.0000000000000000 28.1817131042000000
O Mn Li
160 80 1
Direct
-0.2500000150000000 0.2604309137600000 0.9993153614700000
0.2500000000000000 0.2604309137600000 0.9993153614700000
-0.5000000000000000 0.7604308247600000 0.0006845765300000
0.0000000000000000 0.7604308247600000 0.0006845765300000
-0.2500000150000000 0.9535387172200000 0.0317670078600000
0.2500000000000000 0.9535387172200000 0.0317670078600000
-0.5000000000000000 0.5744181912600000 0.0326922833200000
0.0000000000000000 0.5744181912600000 0.0326922833200000
-0.5000000000000000 0.1181449833100000 0.0408190930700000
0.0000000000000000 0.1181449833100000 0.0408190930700000
-0.2500000150000000 0.4300077190800000 0.0583602348500000
0.2500000000000000 0.4300077190800000 0.0583602348500000
-0.5000000000000000 0.2750738002500000 0.0769547092800000
0.0000000000000000 0.2750738002500000 0.0769547092800000
-0.5000000000000000 0.8281392854100000 0.0824713827300000
0.0000000000000000 0.8281392854100000 0.0824713827300000
-0.2500000150000000 0.6702044481400000 0.0958871600100000
0.2500000000000000 0.6702044481400000 0.0958871600100000
-0.2500000150000000 0.1405752276700000 0.1174443782200000
0.2500000000000000 0.1405752276700000 0.1174443782200000
-0.5000000000000000 0.9835755594300000 0.1176660121600000
0.0000000000000000 0.9835755594300000 0.1176660121600000
-0.5000000000000000 0.5253804691500000 0.1296273580500000
0.0000000000000000 0.5253804691500000 0.1296273580500000
-0.2500000150000000 0.3674343673900000 0.1408288305600000
0.2500000000000000 0.3674343673900000 0.1408288305600000
-0.2500000150000000 0.8420608162000000 0.1606329583400000
0.2500000000000000 0.8419576349900000 0.1605784805800000
-0.4997028944500000 0.6755195246300000 0.1731161806000000
-0.0002971055500000 0.6755195246300000 0.1731161806000000
-0.4993494530800000 0.2228384026800000 0.1827273191000000
-0.0006505469200000 0.2228384026800000 0.1827273191000000
-0.2500000150000000 0.0567571962700000 0.1897838676900000
0.2500000000000000 0.0575895704800000 0.1900331236700000
-0.2500000150000000 0.5333032025200000 0.2109229664700000
0.2500000000000000 0.5345080075600000 0.2112907493000000
-0.5006375979300000 0.3755771804100000 0.2215603738700000
0.0006375979300000 0.3755771804100000 0.2215603738700000
-0.5006449500600000 0.9210834796600000 0.2307882159500000
0.0006449500600000 0.9210834796600000 0.2307882159500000
-0.2500000150000000 0.7637251282700000 0.2385963064700000
0.2500000000000000 0.7630267008500000 0.2384581761200000
-0.2500000150000000 0.2404468496000000 0.2623529424500000
0.2500000000000000 0.2362062498000000 0.2612751170200000
-0.4975560149700000 0.0767166995800000 0.2668567154700000
-0.0024439850300000 0.0767166995800000 0.2668567154700000
-0.4991385972700000 0.6212415261400000 0.2788826594600000
-0.0008614027300000 0.6212415261400000 0.2788826594600000
-0.2500000150000000 0.4597567637600000 0.2886836727700000
0.2500000000000000 0.4620190047000000 0.2890293108200000
-0.2500000150000000 0.9420308867400000 0.3098522501200000
0.2500000000000000 0.9420132621800000 0.3098854168000000
-0.5008512404900000 0.7757784128400000 0.3169074877000000
0.0008512404900000 0.7757784128400000 0.3169074877000000
-0.5120007479200000 0.3173569750900000 0.3325693485300000
0.0120007479200000 0.3173569750900000 0.3325693485300000
-0.2500000150000000 0.1541879345600000 0.3379368137500000
0.2500000000000000 0.1437268142600000 0.3354093128300000
-0.2500000150000000 0.6326832332100000 0.3595766382400000
0.2500000000000000 0.6315279577700000 0.3594197522700000
-0.4984628440600000 0.4751923229400000 0.3703885155200000
-0.0015371559400000 0.4751923229400000 0.3703885155200000
-0.4954804721400000 0.0135521988000000 0.3843985854700000
-0.0045195278600000 0.0135521988000000 0.3843985854700000
-0.2500000150000000 0.8556894410700000 0.3817916758200000
0.2500000000000000 0.8595182822200000 0.3827875333600000
-0.2500000150000000 0.3309474363000000 0.4092602568700000
0.2500000000000000 0.3427162216400000 0.4113952013900000
-0.5107311371500000 0.1770590615200000 0.4146196468500000
0.0107311371500000 0.1770590615200000 0.4146196468500000
-0.5007235647500000 0.7253265118600000 0.4231794258400000
0.0007235647500000 0.7253265118600000 0.4231794258400000
-0.2500000150000000 0.5698794736900000 0.4423095165300000
0.2500000000000000 0.5688584430500000 0.4419577574800000
-0.4994403344800000 0.8821136025400000 0.4594080168700000
-0.0005596655200000 0.8821136025400000 0.4594080168700000
-0.4970074669800000 0.4252692603300000 0.4700633316800000
-0.0029925330200000 0.4252692603300000 0.4700633316800000
-0.2500000150000000 0.0485745071700000 0.4684753548900000
0.2500000000000000 0.0470573140500000 0.4691513219600000
-0.4998792726200000 0.2402746369200000 0.4991804763800000
-0.0001207273800000 0.2402746369200000 0.4991804763800000
-0.2500000150000000 0.7389487129800000 0.5008825702700000
0.2500000000000000 0.7388168342700000 0.5010282686200000
-0.5009409898600000 0.5486889040700000 0.5327471629400000
0.0009409898600000 0.5486889040700000 0.5327471629400000
-0.2500000450000000 0.9257086960000000 0.5328647004300000
0.2499999400000000 0.9256884096200000 0.5328491795900000
-0.2500000450000000 0.3830304341500000 0.5425025370900000
0.2499999400000000 0.3828153074700000 0.5421271456400000
-0.4997549377000000 0.0700565214200000 0.5583883085400000
-0.0002450623000000 0.0700565214200000 0.5583883085400000
-0.2500000450000000 0.2244516310400000 0.5759660206000000
0.2499999400000000 0.2244407893700000 0.5759663325400000
-0.2500000450000000 0.6714880918000000 0.5831662427200000
0.2499999400000000 0.6718877370800000 0.5835194765500000
-0.4998831318600000 0.8308067569500000 0.5960971566500000
-0.0001168681400000 0.8308067569500000 0.5960971566500000
-0.5001422219100000 0.3569873034000000 0.6182750599300000
0.0001422219100000 0.3569873034000000 0.6182750599300000
-0.2500000450000000 0.5158306852300000 0.6180236575100000
0.2499999400000000 0.5159832717600000 0.6181816550300000
-0.2500000450000000 0.9753430820900000 0.6297425345400000
0.2499999400000000 0.9752966900000000 0.6297433605800000
-0.5000111762000000 0.1333632691200000 0.6408525622000000
0.0000111762000000 0.1333632691200000 0.6408525622000000
-0.4999105603300000 0.6579173337600000 0.6609025733300000
-0.0000894396700000 0.6579173337600000 0.6609025733300000
-0.2500000450000000 0.8250026630000000 0.6734695584500000
0.2499999400000000 0.8249755558600000 0.6734698762300000
-0.2500000450000000 0.2768652009800000 0.6832008056000000
0.2499999400000000 0.2768658582100000 0.6832012083300000
-0.4998879791700000 0.4419460398200000 0.6902560700000000
-0.0001120208300000 0.4419460398200000 0.6902560700000000
-0.4999961266600000 0.9670185125200000 0.7106870094100000
-0.0000038733400000 0.9670185125200000 0.7106870094100000
-0.2500000450000000 0.1242758159300000 0.7221543302000000
0.2499999400000000 0.1242593512500000 0.7221545005700000
-0.2500000450000000 0.5771713731000000 0.7302691924500000
0.2499999400000000 0.5770740951500000 0.7302317785200000
-0.5000624819800000 0.7355601111200000 0.7387120463900000
0.0000624819800000 0.7355601111200000 0.7387120463900000
-0.5000009940700000 0.2651459138200000 0.7617478053200000
0.0000009940700000 0.2651459138200000 0.7617478053200000
-0.2500000450000000 0.4229194841300000 0.7694491431100000
0.2499999400000000 0.4229172004800000 0.7694498664600000
-0.2500000450000000 0.8767669274400000 0.7780156426700000
0.2499999400000000 0.8767860982000000 0.7780245411600000
-0.5000079398200000 0.0340563233700000 0.7893162434500000
0.0000079398200000 0.0340563233700000 0.7893162434500000
-0.4999925570100000 0.5593209527700000 0.8095098625200000
-0.0000074429900000 0.5593209527700000 0.8095098625200000
-0.2500000450000000 0.7238728942300000 0.8172397893700000
0.2499999400000000 0.7238711063600000 0.8172419039100000
-0.2500000450000000 0.1758593527700000 0.8268588152100000
0.2499999400000000 0.1758400622300000 0.8268375827300000
-0.4999981412600000 0.3431600858200000 0.8396203642700000
-0.0000018587400000 0.3431600858200000 0.8396203642700000
-0.5000000000000000 0.8674342783900000 0.8591710054400000
0.0000000000000000 0.8674342783900000 0.8591710054400000
-0.2500000450000000 0.0253805131500000 0.8703725379500000
0.2499999400000000 0.0253805131500000 0.8703725379500000
-0.2500000450000000 0.4835755594300000 0.8823338318400000
0.2499999400000000 0.4835755594300000 0.8823338318400000
-0.5000000000000000 0.6405751236700000 0.8825555247800000
0.0000000000000000 0.6405751236700000 0.8825555247800000
-0.5000000000000000 0.1702044921400000 0.9041126909900000
0.0000000000000000 0.1702044921400000 0.9041126909900000
-0.2500000450000000 0.3281393154100000 0.9175284682700000
0.2499999400000000 0.3281393154100000 0.9175284682700000
-0.2500000450000000 0.7750737102500000 0.9230451497200000
0.2499999400000000 0.7750737102500000 0.9230451497200000
-0.5000000000000000 0.9300076300800000 0.9416396241500000
0.0000000000000000 0.9300076300800000 0.9416396241500000
-0.2500000450000000 0.6181448713100000 0.9591807759300000
0.2499999400000000 0.6181448713100000 0.9591807759300000
-0.2500000450000000 0.0744182212600000 0.9673075916800000
0.2499999400000000 0.0744182212600000 0.9673075916800000
-0.5000000000000000 0.4535387172200000 0.9682328621400000
0.0000000000000000 0.4535387172200000 0.9682328621400000
-0.2500000450000000 0.5138390078300000 0.0015253940000000
0.2499999400000000 0.5138390078300000 0.0015253940000000
-0.2500000150000000 0.1993901048400000 0.0566408064400000
0.2500000000000000 0.1993901048400000 0.0566408064400000
-0.2500000150000000 0.9060885862600000 0.0988101646000000
0.2500000000000000 0.9060885862600000 0.0988101646000000
-0.2500000150000000 0.5999232770900000 0.1520907401000000
0.2500000000000000 0.5999232770900000 0.1520907401000000
-0.2500000150000000 0.3008428212100000 0.2004746674000000
0.2500000000000000 0.3004686864600000 0.2004145351100000
-0.2500000150000000 -0.0022830969400000 0.2491242479200000
0.2500000000000000 -0.0016729204600000 0.2494171498800000
-0.2500000150000000 0.6980877640600000 0.2987960628300000
0.2500000000000000 0.6978062403800000 0.2985975161900000
-0.2500000150000000 0.3965984679600000 0.3493347498200000
0.2500000000000000 0.4068340477300000 0.3480838273600000
-0.2500000150000000 0.0957195721100000 0.4007939683200000
0.2500000000000000 0.0857217699000000 0.4037296176300000
-0.2500000150000000 0.8003996363800000 0.4438746991200000
0.2500000000000000 0.8005344734100000 0.4441029618300000
-0.2500000450000000 0.4881036069600000 0.5003589753300000
0.2499999400000000 0.4876745192700000 0.4999236507400000
-0.2500000450000000 0.7974985192500000 0.5572144188000000
0.2499999400000000 0.7978723646300000 0.5572557440700000
-0.2500000450000000 0.0983776311700000 0.6006267809900000
0.2499999400000000 0.0983386391700000 0.6005923660900000
-0.2500000450000000 0.3997851465700000 0.6539410578200000
0.2499999400000000 0.3997465986800000 0.6539951995700000
-0.2500000450000000 0.6984042008100000 0.6996648995300000
0.2499999400000000 0.6983876302600000 0.6996499176200000
-0.2500000450000000 0.0002279144400000 0.7500903477600000
0.2499999400000000 0.0002415967700000 0.7500885161400000
-0.2500000450000000 0.3022636030800000 0.8009843602700000
0.2499999400000000 0.3022617839000000 0.8009846982500000
-0.2500000450000000 0.5998386531500000 0.8466550186300000
0.2499999400000000 0.5998386531500000 0.8466550186300000
-0.2500000450000000 0.9014906408600000 0.8997363362200000
0.2499999400000000 0.9014906408600000 0.8997363362200000
-0.2500000450000000 0.2017832238700000 0.9432545935800000
0.2499999400000000 0.2017832238700000 0.9432545935800000
-0.5000000000000000 0.0138390678300000 0.9984746060000000
0.0000000000000000 0.0138390678300000 0.9984746060000000
-0.5000000000000000 0.7017831788700000 0.0567452724200000
0.0000000000000000 0.7017831788700000 0.0567452724200000
-0.5000000000000000 0.4014907298600000 0.1002635147800000
0.0000000000000000 0.4014907298600000 0.1002635147800000
-0.5000000000000000 0.0998387641500000 0.1533448773700000
0.0000000000000000 0.0998387641500000 0.1533448773700000
-0.5000669106100000 0.8010524584700000 0.1992452563700000
0.0000669106100000 0.8010524584700000 0.1992452563700000
-0.4986108555800000 0.4984111888000000 0.2500021820900000
-0.0013891444200000 0.4984111888000000 0.2500021820900000
-0.4989832002300000 0.1942262581300000 0.2984093250000000
-0.0010167997700000 0.1942262581300000 0.2984093250000000
-0.4994762981300000 0.8981682083700000 0.3464282580400000
-0.0005237018700000 0.8981682083700000 0.3464282580400000
-0.4994496780000000 0.5998279662400000 0.4005252798500000
-0.0005503220000000 0.5998279662400000 0.4005252798500000
-0.4987090812800000 0.2999347938800000 0.4462512577900000
-0.0012909187200000 0.2999347938800000 0.4462512577900000
-0.4985168507000000 0.9857843675000000 0.5022166010300000
-0.0014831493000000 0.9857843675000000 0.5022166010300000
-0.4999401912400000 0.3018009119300000 0.5579178611600000
-0.0000598087600000 0.3018009119300000 0.5579178611600000
-0.5001990915700000 0.5942209115300000 0.5999597389900000
0.0001990915700000 0.5942209115300000 0.5999597389900000
-0.4999923121200000 0.9005072562100000 0.6522621768400000
-0.0000076878800000 0.9005072562100000 0.6522621768400000
-0.4999867041900000 0.1995361821500000 0.7016829122600000
-0.0000132958100000 0.1995361821500000 0.7016829122600000
-0.4999815350800000 0.5003509789100000 0.7502682255200000
-0.0000184649200000 0.5003509789100000 0.7502682255200000
-0.5000113782000000 0.8011423397600000 0.7989555324900000
0.0000113782000000 0.8011423397600000 0.7989555324900000
-0.5000000000000000 0.0999233000900000 0.8479091559000000
0.0000000000000000 0.0999233000900000 0.8479091559000000
-0.5000000000000000 0.4060886152600000 0.9011896864000000
0.0000000000000000 0.4060886152600000 0.9011896864000000
-0.5000000000000000 0.6993900008400000 0.9433590595600000
0.0000000000000000 0.6993900008400000 0.9433590595600000
0.2500000000000000 0.2446480899000000 0.3722180202800000

alex
Hero Member
Hero Member
Posts: 584
Joined: Tue Nov 16, 2004 2:21 pm
License Nr.: 5-67
Location: Germany

Re: DFT+U periodicity (very confused results)

#4 Post by alex » Fri Apr 01, 2016 4:09 pm

Hi jbliu,

my first guess would be that one of your two calculations did not converge or converged to a different electronic structure. Check for that, please.
The structures look identical, but the symmetry the program finds might differ for both cases. Check for that, too. This might result in issues stated above.

Question: It looks like there is some kind of step in the middle of the cell. Is this intended?

Cheers,

alex

admin
Administrator
Administrator
Posts: 2921
Joined: Tue Aug 03, 2004 8:18 am
License Nr.: 458

Re: DFT+U periodicity (very confused results)

#5 Post by admin » Fri Apr 01, 2016 4:21 pm

The structures look OK.
The energies I have calculated (GGA_PBE) are:
energy (converged in steps) magn.moment
strA -1622.019 33 -0.35
strB (produced by myself) -1622.019 32 -0.34
strB (yours) -1622.019 32 -0.34

You should use smaller SIGMA.

jbliu
Newbie
Newbie
Posts: 13
Joined: Tue May 06, 2014 12:00 am
License Nr.: 5-1125

Re: DFT+U periodicity (very confused results)

#6 Post by jbliu » Sat Apr 02, 2016 1:50 pm

admin wrote:The structures look OK.
The energies I have calculated (GGA_PBE) are:
energy (converged in steps) magn.moment
strA -1622.019 33 -0.35
strB (produced by myself) -1622.019 32 -0.34
strB (yours) -1622.019 32 -0.34

You should use smaller SIGMA.
Thank you very much! However, I feel disappointed because I haven't got the same results as yours. I have tried to use smaller SIGMA in INCAR, but the energy and magnetic moment of StructureA and StructureB are still different.

would you mind send your calculation files which contains all the Input and Output files to me? My address is: jbliu@tsinghua.edu.cn

Thank you very much for your help!

The energies I have calculated (GGA_PBE) are:
energy (converged in steps) magn.moment
strA (SIGMA=0.05) -1618.4407 185 DAVsteps -0.6728
strB (SIGMA=0.05) -1617.2765 154 DAVsteps 1.6357

strA (SIGMA=0.02) -1618.3588 319 DAVsteps -0.0589
strB (SIGMA=0.02) -1618.3065 389 DAVsteps -0.1379

jbliu
Newbie
Newbie
Posts: 13
Joined: Tue May 06, 2014 12:00 am
License Nr.: 5-1125

Re: DFT+U periodicity (very confused results)

#7 Post by jbliu » Tue Apr 05, 2016 4:52 pm

admin wrote:The structures look OK.
The energies I have calculated (GGA_PBE) are:
energy (converged in steps) magn.moment
strA -1622.019 33 -0.35
strB (produced by myself) -1622.019 32 -0.34
strB (yours) -1622.019 32 -0.34

You should use smaller SIGMA.
Thank you very much for your help! I have got the same results as yours by using your POTCAR files. We find that there are some differences between our POTCAR files. Could you tell me which version of POTCAR do you use?

Thank you!

admin
Administrator
Administrator
Posts: 2921
Joined: Tue Aug 03, 2004 8:18 am
License Nr.: 458

Re: DFT+U periodicity (very confused results)

#8 Post by admin » Wed Apr 06, 2016 9:19 am

Standard GGA_PBE, nothing special.

alex
Hero Member
Hero Member
Posts: 584
Joined: Tue Nov 16, 2004 2:21 pm
License Nr.: 5-67
Location: Germany

Re: DFT+U periodicity (very confused results)

#9 Post by alex » Wed Apr 06, 2016 1:10 pm

Hell jbliu,

I got the same as admin:
A:
DAV: 34 -0.162201921405E+04 -0.30559E-05 -0.38580E-04 3540 0.898E-02
1 F= -.16220192E+04 E0= -.16217778E+04 d E =-.482814E+00 mag= -0.3559
B:
DAV: 32 -0.162201923641E+04 -0.75992E-04 -0.93062E-04 4224 0.945E-02
1 F= -.16220192E+04 E0= -.16217768E+04 d E =-.484872E+00 mag= -0.3542

Maybe your installation is weird.

Cheers,

alex

jbliu
Newbie
Newbie
Posts: 13
Joined: Tue May 06, 2014 12:00 am
License Nr.: 5-1125

Re: DFT+U periodicity (very confused results)

#10 Post by jbliu » Thu Apr 07, 2016 9:31 am

admin wrote:Standard GGA_PBE, nothing special.
Thank you for your reply. However, We find there indeed exists some differences between your POTCAR files and our POTCAR files (version 5.2 and version 5.4) which are downloaded in the official VASP websites. Would you mind tell us where can download POTCAR versions you used? Or could you send your complete POTCAR files which contain all of available elements besides Mn, O and Li to us? My address is: jbliu@tsinghua.edu.cn

Thank you very much!

The POTCAR file (version 5.2) We used:
PAW_PBE O 08Apr2002
6.00000000000000
parameters from PSCTR are:
VRHFIN =O: s2p4
LEXCH = PE
EATOM = 432.3788 eV, 31.7789 Ry

TITEL = PAW_PBE O 08Apr2002
LULTRA = F use ultrasoft PP ?
IUNSCR = 1 unscreen: 0-lin 1-nonlin 2-no
RPACOR = 1.200 partial core radius
POMASS = 16.000; ZVAL = 6.000 mass and valenz
RCORE = 1.520 outmost cutoff radius
RWIGS = 1.550; RWIGS = 0.820 wigner-seitz radius (au A)
ENMAX = 400.000; ENMIN = 300.000 eV
ICORE = 2 local potential
LCOR = T correct aug charges
LPAW = T paw PP
EAUG = 605.392
DEXC = 0.000
RMAX = 1.553 core radius for proj-oper
RAUG = 1.300 factor for augmentation sphere
RDEP = 1.550 radius for radial grids
RDEPT = 1.329 core radius for aug-charge

Atomic configuration
4 entries
n l j E occ.
1 0 0.50 -514.6923 2.0000
2 0 0.50 -23.9615 2.0000
2 1 0.50 -9.0305 4.0000
3 2 1.50 -9.5241 0.0000
Description
l E TYP RCUT TYP RCUT
0 -23.9615318 23 1.200
0 -9.5240782 23 1.200
1 -9.0304911 23 1.520
1 8.1634956 23 1.520
2 -9.5240782 7 1.500
Error from kinetic energy argument (eV)
NDATA = 100
STEP = 20.000 1.050

The POTCAR file (version 5.2) you used:
PAW_PBE O 08Apr2002
6.00000000000000000
parameters from PSCTR are:
VRHFIN =O: s2p4
LEXCH = PE
EATOM = 432.3788 eV, 31.7789 Ry

TITEL = PAW_PBE O 08Apr2002
LULTRA = F use ultrasoft PP ?
IUNSCR = 0 unscreen: 0-lin 1-nonlin 2-no
RPACOR = .000 partial core radius
POMASS = 16.000; ZVAL = 6.000 mass and valenz
RCORE = 1.520 outmost cutoff radius
RWIGS = 1.550; RWIGS = .820 wigner-seitz radius (au A)
ENMAX = 400.000; ENMIN = 300.000 eV
ICORE = 2 local potential
LCOR = T correct aug charges
LPAW = T paw PP
EAUG = 605.392
DEXC = .000
RMAX = 2.264 core radius for proj-oper
RAUG = 1.300 factor for augmentation sphere
RDEP = 1.550 radius for radial grids
QCUT = -5.520; QGAM = 11.041 optimization parameters

Description
l E TYP RCUT TYP RCUT
0 .000 23 1.200
0 -.700 23 1.200
1 .000 23 1.520
1 .600 23 1.520
2 .000 7 1.500
Error from kinetic energy argument (eV)
NDATA = 100
STEP = 20.000 1.050


Post Reply