[Wannier] Spin polarized calculations with f-electrons
Mario Burbano
burbanom at tcd.ie
Thu Mar 29 11:09:13 CEST 2012
Dear all,
Over the last few months I have been working on doped (Y, Gd, Sm, La, etc)
and reduced ceria systems. My goal is to perform DFT calculations on these
compounds and then extract the ionic forces and dipoles, which I will then
use to fit an interionic potential model. To this end, I wanted to use the
WCs to extract the dipoles. I have carried out this work using mainly
VASP/Wannier90. I also used CPMD (which can output the WCs) in some cases
to check that the results are consistent. My motivation for using
VASP/Wannier90 is the implementation of the Hubbard (U) correction and the
availability of pseudopotentials for lanthanoid elements.
The results I have obtained so far for yttrium-doped CeO2 and
scandium-doped CeO2 have been consistent with each other using either DFT
program. By this I mean that the forces and the dipole moments obtained
from the DFT codes show excellent agreement between them. However, it is
when I carry out simulations with elements that have f-electrons (Ce3+,
Gd3+, Sm3+) that I start having some serious problems. In the case of
Gd-doped ceria, when I visualize the wannier90_centres.xyz file from the
Wannier90 analysis of the data, I can see that the f-electrons that should
have been assigned to Gd, have instead been localized near the nucleus of
the Ce atoms, in addition to the expected 8 valence electrons. I have
encountered similar problems with systems that include reduced cerium
(Ce3+). In this case the single f-electrons from each Ce3+ are assigned to
a single Ce atom (which according to the VASP output should not be
reduced). These wannier centres are conspicuous among the others given that
their spreads are typically higher (>4 Bohr^2) than those of all the other
WCs in the system (2.5 – 3 Bohr^2). In addition, the up and down spin
components from the wannier90 calculation are the same, which, I think, is
unexpected. Finally, a comparison with the results from CPMD is not
available here because CPMD does not have the DFT+U method implemented.
I am quite confident that the VASP calculations I performed are fine.
Indeed, I have plotted the electronic density of states for said
calculations from VASP and I found them to be in accordance with what would
be expected, i.e. in the case of Gd-doped ceria, the Gd f-states are deep
in the valence band. The calculation also showed charges consistent with Ce
being 4+ and Gd being 3+ and the magnetization of the Gd ions is 7, which
is consistent with the literature.
My question is if there is a specific set of keywords that I need to use
for these f-electron systems. To this end, I attach my wannier90.win,
wannier90.wout, INCAR file from VASP. A plot of the EDOS for Gd-doped ceria
is also available upon request. I would greatly appreciate some help with
this issue.
Mario Burbano
INCAR
SYSTEM = CeO2
Start parameter for this Run:
NWRITE = 2
LPETIM=F write-flag & timer
ISTART = 1 job : 0-new 1-cont 2-samecut
ICHARG = 1
Electronic Relaxation 1
PREC = normal
ENMAX = 500.00 eV Use no more than 400 for hyb
NELM = 150 Or use 150 depending on system
INIWAV = 1 random initial wavefunctions
EDIFF = 1E-05
Ionic Relaxation
EDIFFG = -0.01 stopping-criterion for IOM
NSW = 0 number of steps for IOM
NBLOCK = 1; KBLOCK = 40 inner block and outer block
IBRION = 2 ionic relax: 0-MD 1-quasi-New 2-CG
!POTIM = 0.2
!TEBEG = 1200
ISIF = 2 stress and relaxation
IWAVPR = 1 prediction: 0-non 1-charg 2-wave 3-comb
ISYM = 0 0-nonsym 1-usesym
LCORR = F Harris-correction to forces
ISMEAR = 0
LMAXMIX = 6
!NPAR = 4
ALGO = Normal
LREAL = A
ROPT = 1E-04 1E-04 1E-04
ISPIN = 2
!MAGMOM = 63*0.0 2*7.0 30*0.0
LORBIT = 12
!NBANDS = 576
!LWANNIER90 = .TRUE.
LWANNIER90_RUN = .TRUE.
LWRITE_MMN_AMN = .TRUE.
LASPH = .TRUE.
LORBIT = 12
LDAU = .TRUE. DFT+U Parameters
LDAUL = -1 3 3
LDAUU = 0 5 7
LDAUJ = 0 0 0
LMAXMIX = 6
LASPH = .TRUE.
wannier90.win
num_wann = 394
exclude_bands : 395-516
translate_home_cell = T
write_xyz = T
num_iter = 10000
!restart = default
length_unit = bohr
conv_tol = 1E-05
conv_window = 2
use_bloch_phases = T
begin unit_cell_cart
11.0999000 0.0000000 0.0000000
0.0000000 11.0999000 0.0000000
0.0000000 0.0000000 11.0999000
end unit_cell_cart
begin atoms_cart
O 1.0910258 1.5555279 1.1133891
O 1.6168928 1.4308410 6.8344120
O 1.2726312 7.1209920 1.5239751
O 6.7235477 1.3062118 1.5243693
O 6.9849921 1.2755700 7.0070873
O 7.0117723 6.9623947 1.3763876
O 6.9734855 6.8747087 7.0389660
O 1.2816080 4.2024615 1.5141458
O 1.3287112 4.3624672 6.9362300
O 1.2972118 9.4526459 1.3483589
O 1.3706034 9.6482123 6.9524353
O 7.1090675 4.1821279 1.1530493
O 7.1452833 4.1399094 6.6534944
O 6.8159130 9.5397849 1.4425250
O 7.0278822 9.4519766 6.9553476
O 3.8301972 1.6956449 1.5085080
O 4.1800092 1.2534059 6.6640303
O 4.2504328 6.8128734 1.5556019
O 3.9247379 6.7899324 6.9399874
O 9.5556851 1.1386685 1.2527318
O 9.6036212 1.1291998 6.7260156
O 9.5866608 7.0162440 1.2348022
O 9.7441521 6.6994064 6.7313684
O 4.0548661 4.3197004 1.1723938
O 4.1273180 4.3148634 6.8252434
O 4.1772558 9.8811025 1.4115724
O 4.2169280 9.6342031 7.0420338
O 9.7559647 4.2100936 1.2570907
O 9.8481326 4.0300227 6.9914025
O 9.5395371 9.6699812 1.3797357
O 9.7535514 9.5559178 7.0318532
O 1.6712692 1.4300929 4.0238828
O 1.2543316 1.3619023 9.5704645
O 1.4847537 6.7609402 4.3163762
O 1.2548616 6.9356260 9.7869100
O 7.2461518 1.3818800 4.2487517
O 6.7620253 1.5077195 9.8902004
O 7.1182846 7.1713511 4.3022937
O 7.0775316 6.9802018 9.7163932
O 1.3355604 3.9342011 4.1256461
O 1.3149627 4.1328975 9.5825493
O 1.5207152 9.3671966 4.1077831
O 1.2353957 9.4353610 9.7722335
O 6.9859142 4.3444563 3.9066656
O 7.0419841 4.2985905 9.5477929
O 6.9607615 9.9409150 4.4807489
O 6.8498114 9.9006835 9.7935517
O 4.2730595 1.7229243 4.1831688
O 3.9713144 1.4241965 9.7405612
O 4.1850513 6.8978010 4.0957753
O 4.0475432 6.8265599 9.7106643
O 10.0315959 1.3967283 3.8460075
O 9.5551204 1.2485580 9.3996756
O 9.7999703 6.8490059 3.9202658
O 9.6494379 6.9913007 9.6843607
O 4.2028441 4.1931154 4.0845843
O 4.1862072 4.3245559 9.5876039
O 4.1470296 9.7356322 4.3733518
O 4.1740254 9.6720670 9.7716207
O 9.7654108 4.1122673 4.1680258
O 9.7818919 3.9942592 9.7749854
O 9.5273401 9.6357585 4.1431350
O 9.6338072 9.5554361 9.5372859
Gd 5.7302426 2.7687112 2.8334760
Gd 10.9789137 10.9654572 0.0105097
Ce 0.1778309 0.3185873 5.4495023
Ce 11.0260132 5.4730070 11.0700061
Ce 10.9607076 5.5356176 5.4386890
Ce 5.4440578 0.0834043 0.1366895
Ce 5.6382141 0.1742708 5.4959768
Ce 5.4885401 5.5045347 0.0919183
Ce 5.7278683 5.5908857 5.5776579
Ce 2.7531946 2.6635696 11.0489633
Ce 2.7925504 2.9913267 5.4315247
Ce 2.7046285 8.3944282 0.0673516
Ce 2.8253220 8.2792851 5.6985709
Ce 8.2519614 2.7117986 11.0442467
Ce 8.4126651 2.7737429 5.5618419
Ce 8.1952205 8.1584337 0.0519898
Ce 8.3428229 8.4431157 5.6094422
Ce 3.0058818 11.0854546 2.9571521
Ce 2.7320840 0.0080979 8.2708000
Ce 2.9056957 5.4776607 2.7579451
Ce 2.6843252 5.6239897 8.3753906
Ce 8.4102373 11.0299493 2.7425168
Ce 8.2087693 10.9612438 8.3471079
Ce 8.3530073 5.5853908 2.6058968
Ce 8.3022411 5.5315589 8.3088875
Ce 11.0936636 2.8089634 2.6070143
Ce 11.0195643 2.5974887 8.4131393
Ce 0.0012041 8.3390840 2.8207236
Ce 0.0256974 8.3442472 8.3057507
Ce 5.3643603 2.7681623 8.3381502
Ce 5.5863572 8.2839246 2.8753859
Ce 5.5999122 8.4307771 8.3534403
end atoms_cart
mp_grid = 1 1 1
begin kpoints
0.000000000000 0.000000000000 0.000000000000
end kpoints
Truncated wannier90.wout
WF centre and spread 315 ( 2.955411, 7.977270, -8.320092 )
3.10397663
WF centre and spread 316 ( 1.983168, -7.084571, 3.330386 )
2.82780519
WF centre and spread 317 ( 8.389708, -8.287388, 8.121573 )
3.04872989
WF centre and spread 318 ( -7.973836, 2.911651, 2.511068 )
3.06077796
WF centre and spread 319 ( -3.401918, 2.461755, 1.944609 )
2.99807308
WF centre and spread 320 ( -2.747481, -7.649929, 7.787093 )
3.10885119
WF centre and spread 321 ( -2.943842, 7.137560, -2.093789 )
2.88281541
WF centre and spread 322 ( -1.987234, -8.089698, -7.978904 )
3.10697570
WF centre and spread 323 ( 2.773243, -8.264631, -2.946507 )
2.99032824
WF centre and spread 324 ( -7.568257, 3.066973, 7.682292 )
2.93448218
WF centre and spread 325 ( 2.813919, 8.344642, 3.277612 )
2.91319501
WF centre and spread 326 ( 10.053932, -5.344354, -4.650336 )
2.47476242
WF centre and spread 327 ( 8.056433, -7.708514, -3.028838 )
3.09028990
WF centre and spread 328 ( -2.112588, 7.406764, 7.518556 )
2.92660667
WF centre and spread 329 ( -7.769204, -3.388896, 2.335361 )
2.87424976
WF centre and spread 330 ( -3.338727, 1.952770, -3.644332 )
2.75161953
WF centre and spread 331 ( -2.122343, 8.371526, 1.936664 )
2.79772379
WF centre and spread 332 ( -2.123271, 8.223344, 8.318030 )
2.95807968
WF centre and spread 333 ( 7.649308, -7.742433, 2.587871 )
3.09400540
WF centre and spread 334 ( -8.553601, -3.383708, 3.220682 )
2.89142249
WF centre and spread 335 ( 2.170253, 7.862902, 8.177933 )
3.01993865
WF centre and spread 336 ( -8.170698, 2.009430, -8.168123 )
2.87459996
WF centre and spread 337 ( -3.378961, -3.125420, 2.321571 )
2.97235965
WF centre and spread 338 ( 7.314661, 7.718848, 1.837712 )
2.82518456
WF centre and spread 339 ( 7.601711, 3.611174, 8.276649 )
2.86500869
WF centre and spread 340 ( 8.260864, -2.292446, -2.068143 )
2.98322396
WF centre and spread 341 ( 1.735624, 3.413602, 2.563879 )
2.83866685
WF centre and spread 342 ( 1.996838, -7.902926, 2.514744 )
3.08720409
WF centre and spread 343 ( -2.013425, 7.984178, -8.201266 )
2.99172689
WF centre and spread 344 ( 7.425085, 8.514192, -7.639445 )
2.93024986
WF centre and spread 345 ( 2.480124, -3.691701, 7.358987 )
2.80989381
WF centre and spread 346 ( 8.017905, -8.470672, -2.185893 )
3.00688507
WF centre and spread 347 ( 8.303022, 8.539391, -8.446310 )
2.95251901
WF centre and spread 348 ( -8.030027, -8.106592, -2.203690 )
3.06248896
WF centre and spread 349 ( 4.841020, 5.442349, 0.532952 )
2.44329218
WF centre and spread 350 ( -3.010810, -2.561720, -7.310560 )
2.92096297
WF centre and spread 351 ( -8.281606, -1.688036, 8.869493 )
2.77204423
WF centre and spread 352 ( -3.208049, 1.779769, -7.877091 )
2.84982433
WF centre and spread 353 ( -1.530521, 2.207644, 7.554694 )
2.96087159
WF centre and spread 354 ( -8.181692, 7.923215, 6.975434 )
2.79038804
WF centre and spread 355 ( 7.982252, 8.693198, 1.757743 )
2.77611846
WF centre and spread 356 ( -7.291098, -8.247983, 2.951081 )
2.82219938
WF centre and spread 357 ( -7.448289, -8.447889, -7.259313 )
2.86047293
WF centre and spread 358 ( 3.242955, -8.553512, 7.805895 )
3.03592398
WF centre and spread 359 ( -2.346036, -8.256709, -2.318301 )
3.06070397
WF centre and spread 360 ( -2.542093, -3.172905, 7.484967 )
2.90458533
WF centre and spread 361 ( 3.155946, -7.818175, 8.713861 )
3.07195605
WF centre and spread 362 ( 2.895759, -2.745856, 2.889512 )
3.18059131
WF centre and spread 363 ( 2.363604, 2.550049, 2.458465 )
3.30535165
WF centre and spread 364 ( -3.396423, -3.160862, 8.253963 )
2.80642521
WF centre and spread 365 ( 7.623483, 3.579213, 3.094698 )
2.82291720
WF centre and spread 366 ( 2.111921, 8.714749, -8.232192 )
3.05464852
WF centre and spread 367 ( -2.083752, 7.098527, -2.888565 )
2.87670157
WF centre and spread 368 ( -2.512735, -8.164707, 1.935589 )
2.93885801
WF centre and spread 369 ( -2.080538, 7.957223, -2.080406 )
2.85354471
WF centre and spread 370 ( 2.422977, -7.808186, 7.833933 )
3.10047156
WF centre and spread 371 ( 8.240706, 7.585315, 7.308912 )
2.86138806
WF centre and spread 372 ( -2.618142, -2.336775, 2.192220 )
2.83746790
WF centre and spread 373 ( 2.017842, 8.321124, 2.425184 )
2.96342095
WF centre and spread 374 ( -3.338213, -7.343168, 1.909590 )
2.77140536
WF centre and spread 375 ( 2.584014, 3.353781, 1.732777 )
2.85955278
WF centre and spread 376 ( 2.896469, 8.242927, -3.268844 )
2.94442950
WF centre and spread 377 ( -2.535725, -3.013976, 3.065033 )
3.01193621
WF centre and spread 378 ( -2.424804, -2.518216, -2.609676 )
2.96225890
WF centre and spread 379 ( -3.149585, -3.322517, -2.601553 )
3.05900124
WF centre and spread 380 ( 2.822494, -7.100193, 2.435167 )
2.98763778
WF centre and spread 381 ( 3.505806, -6.903623, -8.067028 )
11.46708747
WF centre and spread 382 ( 10.176170, 5.241469, -5.221244 )
2.42578292
WF centre and spread 383 ( 5.235784, -0.065968, -5.354240 )
3.69196155
WF centre and spread 384 ( 10.112686, 5.240451, -5.221938 )
2.92049603
WF centre and spread 385 ( 10.141411, 5.227326, -5.245713 )
2.85520635
WF centre and spread 386 ( 9.722165, 4.820787, -5.230830 )
6.16646019
WF centre and spread 387 ( 1.361169, 0.518060, 9.417393 )
9.99612809
WF centre and spread 388 ( 10.110982, 5.170869, -5.222982 )
2.96430550
WF centre and spread 389 ( 10.112519, 5.273404, -5.282557 )
4.20368371
WF centre and spread 390 ( 4.723585, -0.037301, -6.030209 )
9.74295322
WF centre and spread 391 ( 0.392734, 0.529797, 10.290223 )
2.79897493
WF centre and spread 392 ( 0.557381, 0.533053,-10.451580 )
5.65979199
WF centre and spread 393 ( 0.487947, 0.548060, 10.254358 )
4.40782260
WF centre and spread 394 ( 3.062461, 0.220112, -8.035170 )
14.14469343
Sum of centres and spreads ( 2.577830, 49.036143, -1.714244 )
1137.06428007
Spreads (Bohr^2) Omega I = 949.037551765
================ Omega D = 0.000000000
Omega OD = 188.026728302
Final Spread (Bohr^2) Omega Total = 1137.064280067
------------------------------------------------------------------------------
Wannier centres written to file wannier90.up_centres.xyz
Writing checkpoint file wannier90.up.chk... done
Time for wannierise 1946.970 (sec)
Total Execution Time 1947.812 (sec)
*===========================================================================*
| TIMING
INFORMATION |
*===========================================================================*
| Tag Ncalls Time
(s)|
|---------------------------------------------------------------------------|
|kmesh: get : 2
0.013|
|wann: main : 1
1946.963|
*---------------------------------------------------------------------------*
All done: wannier90 exiting
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/wannier/attachments/20120329/4cc39337/attachment.html>
More information about the Wannier
mailing list