huge dE value in ionic steep

Question on input files/tags, interpreting output, etc.

Please check whether the answer to your question is given in the VASP online manual or has been discussed in this forum previously!

Moderators: Global Moderator, Moderator

Locked
Message
Author
scpeter@jncasr
Newbie
Newbie
Posts: 9
Joined: Tue Apr 12, 2022 9:13 am

huge dE value in ionic steep

#1 Post by scpeter@jncasr » Wed Oct 05, 2022 4:02 am

Dear all
I am optimizing a slab model of perovskite and getting a weird E0 value though the E value at the electronic convergence step looks good. I dont know why the dE for each ionic step is coming huge and at the same time it says ionic convergence is reached. plz .have a look at the output and INCAR file.
thanku you
regards
Riyaz
Global Parameters
OSZICAR:
17 F= -.56032353E+03 E0= -.56031244E+03 d E =-.913020E-04 mag= 49.0000
N E dE d eps ncg rms rms(c)
DAV: 1 -0.550766619959E+03 0.19048E-02 -0.78478E-02 8904 0.100E+00 0.150E-01
RMM: 2 -0.550768291353E+03 -0.16714E-02 -0.17108E-02 12249 0.397E-01 0.746E-02
RMM: 3 -0.550768611867E+03 -0.32051E-03 -0.25762E-03 11439 0.968E-02 0.729E-02
RMM: 4 -0.550768628373E+03 -0.16505E-04 -0.33588E-04 10923 0.346E-02 0.588E-02
RMM: 5 -0.550768619514E+03 0.88590E-05 -0.90836E-05 8248 0.220E-02 0.374E-02
RMM: 6 -0.550768618828E+03 0.68612E-06 -0.38328E-05 6360 0.174E-02 0.328E-02
RMM: 7 -0.550768610193E+03 0.86347E-05 -0.29128E-05 5753 0.179E-02 0.182E-02
RMM: 8 -0.550768615162E+03 -0.49693E-05 -0.38522E-05 6045 0.192E-02 0.764E-03
RMM: 9 -0.550768618745E+03 -0.35829E-05 -0.11797E-05 5118 0.191E-02 0.502E-03
RMM: 10 -0.550768621619E+03 -0.28736E-05 -0.65058E-06 5048 0.149E-02
18 F= 0.41399501E+05 E0= 0.41399512E+05 d E =0.419598E+05 mag= 49.0000
N E dE d eps ncg rms rms(c)
DAV: 1 -0.550768512279E+03 0.10647E-03 -0.78460E-03 9768 0.277E-01 0.400E-02
RMM: 2 -0.550768638781E+03 -0.12650E-03 -0.12791E-03 11765 0.111E-01 0.189E-02
RMM: 3 -0.550768661156E+03 -0.22375E-04 -0.18910E-04 8958 0.323E-02 0.121E-02
RMM: 4 -0.550768664835E+03 -0.36784E-05 -0.39290E-05 6011 0.215E-02 0.119E-02
RMM: 5 -0.550768665618E+03 -0.78306E-06 -0.10012E-05 5048 0.162E-02 0.922E-03
RMM: 6 -0.550768666252E+03 -0.63475E-06 -0.62146E-06 5006 0.142E-02 0.793E-03
RMM: 7 -0.550768667270E+03 -0.10180E-05 -0.39008E-06 5005 0.148E-02 0.473E-03
RMM: 8 -0.550768668679E+03 -0.14087E-05 -0.34638E-06 4999 0.132E-02 0.306E-03
RMM: 9 -0.550768669892E+03 -0.12125E-05 -0.27761E-06 4974 0.133E-02 0.206E-03
RMM: 10 -0.550768670752E+03 -0.86061E-06 -0.13884E-06 4883 0.105E-02
19 F= 0.24167030E+03 E0= 0.24168139E+03 d E =-.411578E+05 mag= 49.0000
N E dE d eps ncg rms rms(c)
DAV: 1 -0.550767747344E+03 0.92255E-03 -0.52294E-02 8744 0.712E-01 0.101E-01
RMM: 2 -0.550768581639E+03 -0.83429E-03 -0.84832E-03 12227 0.282E-01 0.481E-02
RMM: 3 -0.550768720945E+03 -0.13931E-03 -0.11929E-03 11421 0.679E-02 0.318E-02
RMM: 4 -0.550768731184E+03 -0.10239E-04 -0.17204E-04 9395 0.250E-02 0.272E-02
RMM: 5 -0.550768731358E+03 -0.17414E-06 -0.44716E-05 6398 0.189E-02 0.238E-02
RMM: 6 -0.550768729104E+03 0.22537E-05 -0.19599E-05 5634 0.151E-02 0.216E-02
RMM: 7 -0.550768724752E+03 0.43520E-05 -0.15661E-05 5301 0.164E-02 0.150E-02
RMM: 8 -0.550768722606E+03 0.21467E-05 -0.12226E-05 5165 0.177E-02 0.609E-03
RMM: 9 -0.550768723775E+03 -0.11694E-05 -0.77207E-06 5024 0.153E-02 0.357E-03
RMM: 10 -0.550768724974E+03 -0.11986E-05 -0.44874E-06 4975 0.134E-02
20 F= 0.41141483E+04 E0= 0.41141594E+04 d E =-.372854E+05 mag= 49.0000
N E dE d eps ncg rms rms(c)
DAV: 1 -0.550768683622E+03 0.40153E-04 -0.22529E-03 10448 0.161E-01 0.248E-02
RMM: 2 -0.550768727084E+03 -0.43462E-04 -0.42143E-04 10668 0.636E-02 0.148E-02
RMM: 3 -0.550768736413E+03 -0.93290E-05 -0.63917E-05 6744 0.222E-02 0.189E-02
RMM: 4 -0.550768737235E+03 -0.82194E-06 -0.16896E-05 5290 0.203E-02 0.876E-03
RMM: 5 -0.550768738690E+03 -0.14549E-05 -0.65520E-06 4941 0.137E-02 0.569E-03
RMM: 6 -0.550768739716E+03 -0.10257E-05 -0.34901E-06 4824 0.120E-02 0.482E-03
RMM: 7 -0.550768741599E+03 -0.18830E-05 -0.18583E-06 4826 0.107E-02 0.358E-03
RMM: 8 -0.550768743277E+03 -0.16781E-05 -0.19131E-06 4900 0.118E-02 0.232E-03
RMM: 9 -0.550768744695E+03 -0.14178E-05 -0.15814E-06 4811 0.975E-03 0.181E-03
RMM: 10 -0.550768745510E+03 -0.81571E-06 -0.12021E-06 4680 0.904E-03
21 F= -.55336126E+03 E0= -.55335018E+03 d E =-.466751E+04 mag= 49.0000
N E dE d eps ncg rms rms(c)
DAV: 1 -0.550768348171E+03 0.39652E-03 -0.17658E-02 9152 0.447E-01 0.674E-02
RMM: 2 -0.550768684053E+03 -0.33588E-03 -0.33644E-03 12063 0.176E-01 0.405E-02
RMM: 3 -0.550768748485E+03 -0.64432E-04 -0.48064E-04 10612 0.464E-02 0.532E-02
RMM: 4 -0.550768744403E+03 0.40818E-05 -0.75723E-05 7243 0.219E-02 0.247E-02
RMM: 5 -0.550768743510E+03 0.89304E-06 -0.13979E-05 5228 0.179E-02 0.163E-02
RMM: 6 -0.550768743411E+03 0.99215E-07 -0.66295E-06 4993 0.142E-02 0.139E-02
RMM: 7 -0.550768742763E+03 0.64787E-06 -0.71515E-06 5085 0.163E-02 0.749E-03
RMM: 8 -0.550768743241E+03 -0.47790E-06 -0.59612E-06 5049 0.153E-02 0.380E-03
RMM: 9 -0.550768743673E+03 -0.43202E-06 -0.37312E-06 4992 0.146E-02 0.256E-03
RMM: 10 -0.550768743863E+03 -0.19011E-06 -0.23165E-06 5004 0.124E-02
22 F= -.56032375E+03 E0= -.56031268E+03 d E =-.467447E+04 mag= 49.0000



input: ISTART = 1
ICHARG=2
ISPIN = 2
LMAXMIX = 4
LASPH = .TRUE.
MAGMOM = 16*4.5 33*0 11*0 12*0
LDAU = .TRUE.
LDAUL = 2 -1 -1 -1
LDAUU = 3.9 0.0 0.0 0.0
LDAUJ = 0.0 0.0 0.0 0.0
LCHARG = .TRUE.
ADDGRID= .TRUE.
#LSCALAPACK = .FALSE.
ALGO = Fast
IDIPOL = 3
#LDIPOL = .TRUE.
LREAL = Auto
Electronic Relaxation
LORBIT = 11
ISMEAR = 0
SIGMA = 0.05
ENCUT = 450
NELM = 150
NELMIN = 10
EDIFF = 1E-05
Ionic Relaxation
NSW = 100
IBRION = 2
ISIF = 2
IVDW = 11
NCORE = 4

martin.schlipf
Global Moderator
Global Moderator
Posts: 458
Joined: Fri Nov 08, 2019 7:18 am

Re: huge dE value in ionic steep

#2 Post by martin.schlipf » Wed Oct 05, 2022 7:34 am

Can you provide a complete set of input (INCAR, POSCAR, KPOINTS, POTCAR) and output files (OUTCAR, standard output), please?

scpeter@jncasr
Newbie
Newbie
Posts: 9
Joined: Tue Apr 12, 2022 9:13 am

Re: huge dE value in ionic steep

#3 Post by scpeter@jncasr » Wed Oct 05, 2022 12:04 pm

thanks for the reply, yes
plz. have a look, I have shared only the order for POTCAR, not the whole file due to the word limit.

INCAR:
Global Parameters
ISTART = 1
ICHARG=2
ISPIN = 2
LMAXMIX = 4
LASPH = .TRUE.
MAGMOM = 16*4.5 33*0 11*0 12*0
LDAU = .TRUE.
LDAUL = 2 -1 -1 -1
LDAUU = 3.9 0.0 0.0 0.0
LDAUJ = 0.0 0.0 0.0 0.0
LCHARG = .TRUE.
ADDGRID= .TRUE.
#LSCALAPACK = .FALSE.
ALGO = Fast
IDIPOL = 3
#LDIPOL = .TRUE.
LREAL = Auto
Electronic Relaxation
LORBIT = 11
ISMEAR = 0
SIGMA = 0.05
ENCUT = 450
NELM = 150
NELMIN = 10
EDIFF = 1E-05
Ionic Relaxation
NSW = 100
IBRION = 2
ISIF = 2
IVDW = 11
NCORE = 4


POSCAR:
Wien97 struct file generated by XCrySDe
1.00000000000000
7.7521495819000004 0.0000000000000000 0.0000000000000000
0.0000000000000000 7.7521495819000004 0.0000000000000000
0.0000000000000000 0.0000000000000000 27.1471996307000012
Mn O N La
16 33 11 12
Selective dynamics
Direct
0.4984756248774939 0.0009386686122902 0.7276062549904201 T T T
0.9987641468479783 0.0009500515402375 0.7258555810754509 T T T
0.4993208936418975 0.5009649003253309 0.7277180818819950 T T T
0.9988751173661683 0.5007962654878245 0.7216414644997889 T T T
0.9991607044276339 0.0013256797260156 0.5775899244571542 T T T
0.4994418125271170 0.0010032504177420 0.5788129005151070 T T T
0.9958111627767176 0.5025513451265416 0.5773152497852010 T T T
0.4974136100972720 0.5011735154689472 0.5784770535955174 T T T
0.0023881780602090 0.4996052572023554 0.4326068851197603 T T T
0.5024419472044456 0.5017595107305274 0.4324439210055306 T T T
0.0007767565054100 0.0002242689009406 0.4311437845324114 T T T
0.4992437574869733 0.0014426381149231 0.4309118853661640 T T T
0.9980714540505260 0.5029699411820612 0.2882987603271742 T T T
0.4982381365625942 0.5031132288510801 0.2838979541948280 T T T
0.4979311586302221 0.0031907165895431 0.2819585714357674 T T T
0.9992035745895149 0.0036555384484857 0.2821844441270070 T T T
0.2492689609104135 0.5106161585657356 0.7178339440762694 T T T
0.7494348943744314 0.4926428486961293 0.7150270956719886 T T T
0.2485640831851211 0.9930118366054566 0.7164870310665180 T T T
0.7487449273413472 0.0090704784304589 0.7170331334286438 T T T
0.0077906590620987 0.2506403657521084 0.7189306517091941 T T T
0.9888112035075933 0.7505358350497242 0.7177642428068719 T T T
0.9969645216190829 0.0022253300237282 0.6468477004917140 T T T
0.5012089557991932 -0.0003755331762534 0.6490946705989419 T T T
0.4939237835145441 0.5013747479695850 0.6483169551443719 T T T
0.9422835445127948 0.2497698025247356 0.5770851499971493 T T T
0.0535086674937056 0.7538549188906059 0.5774960358572593 T T T
0.5529328630272687 0.2507725074579472 0.5775437004192037 T T T
0.4438833860650941 0.7513658424910281 0.5774790124978798 T T T
0.2489819373174706 0.4483647929935294 0.5760043041217736 T T T
0.7464521611160345 0.5566348791947426 0.5779698842017986 T T T
0.4999998810000008 0.4999998810000008 0.5058692100000002 F F F
0.4999998810000008 0.0000000000000000 0.5052256580000005 F F F
0.0000000000000000 0.0000000000000000 0.5048177839999965 F F F
0.0513237796373138 0.2493071423528148 0.4323088090183805 T T T
0.9527733810611453 0.7519239226998498 0.4337376403226259 T T T
0.2502521617809297 0.9512579770231189 0.4339021454677154 T T T
0.7502824465782857 0.0507443176342831 0.4305951001937353 T T T
0.2532293871541245 0.5523404107309651 0.4308601799231061 T T T
0.7502169081555002 0.4483687585565788 0.4344765916449893 T T T
0.5098697713516229 0.5013164016791513 0.3629937999797546 T T T
0.0099508550903918 -0.0017787367818005 0.3613649956260242 T T T
0.4871832541072508 0.0032532451212651 0.3605949789756238 T T T
0.9910021355925250 0.2537629697888356 0.2946258359777232 T T T
0.0089641317913192 0.7528542446035000 0.2924555695475985 T T T
0.5050891801805415 0.2534883296765553 0.2930022332101508 T T T
0.4893305204935494 0.7527144687284975 0.2930730538106243 T T T
0.2478882905778342 0.4941203680556961 0.2942741345545674 T T T
0.7477040702697276 0.5127692704849154 0.2887018037301980 T T T
0.4923826309345909 0.2515688699351660 0.7326870468245467 T T T
0.5058010412444087 0.7508159101303999 0.7327044559530067 T T T
0.0075035777360152 0.4963468531079207 0.6479776291948620 T T T
0.2485999827258811 0.0493488097958656 0.5772034160529498 T T T
0.7494264217428814 0.9532580810866649 0.5767343418966365 T T T
0.0000000000000000 0.4999998810000008 0.5055875179999987 F F F
0.4550627708353019 0.2518457890958791 0.4329634074908301 T T T
0.5471047722556647 0.7508249731605641 0.4325723249859217 T T T
0.9840121400801353 0.5085682316963817 0.3620403306074791 T T T
0.2481660734386711 0.0101535506869808 0.2753179791831822 T T T
0.7490469661538360 0.9961728932221988 0.2795727374067022 T T T
0.2449106589099752 0.2534983459751270 0.6557009933620980 T T T
0.7567566649011345 0.2486606891843558 0.6541553028740792 T T T
0.2417801190571889 0.7514367472948839 0.6544186430486983 T T T
0.7515729728777427 0.7521149327323433 0.6553753064137651 T T T
0.2499277000000006 0.2479888800000012 0.5039203759999964 F F F
0.7500722999999994 0.2479888800000012 0.5039203759999964 F F F
0.2499277000000006 0.7520111199999988 0.5039203759999964 F F F
0.7500722999999994 0.7520111199999988 0.5039203759999964 F F F
0.2510166510889366 0.2543688768539046 0.3550454817319292 T T T
0.7469272405051859 0.2630324854421234 0.3552159172707253 T T T
0.2495763462139748 0.7471274855114487 0.3555463851005752 T T T
0.7497242735329948 0.7495776729810868 0.3532216392661911 T T T

KPOINTS
K-Spacing Value to Generate K-Mesh: 0.040
0
Monkhorst-Pack
3 3 1
0.0 0.0 0.0

POTCAR
PAW_PBE Mn_pv 02Aug2007
13.0000000000000

PAW_PBE O 08Apr2002
6.00000000000000

PAW_PBE N 08Apr2002
5.00000000000000

PAW_PBE La 06Sep2000
11.0000000000000

scpeter@jncasr
Newbie
Newbie
Posts: 9
Joined: Tue Apr 12, 2022 9:13 am

Re: huge dE value in ionic steep

#4 Post by scpeter@jncasr » Wed Oct 05, 2022 12:26 pm

https://drive.google.com/file/d/1rIREQc ... sp=sharing

you can also get file from my drive via this link.
Thank you

martin.schlipf
Global Moderator
Global Moderator
Posts: 458
Joined: Fri Nov 08, 2019 7:18 am

Re: huge dE value in ionic steep

#5 Post by martin.schlipf » Thu Oct 06, 2022 7:47 am

The link is dead
Sorry, the file you have requested does not exist.

Make sure that you have the correct URL and the file exists.
Please make sure to include the output of the calculation as well.

scpeter@jncasr
Newbie
Newbie
Posts: 9
Joined: Tue Apr 12, 2022 9:13 am

Re: huge dE value in ionic steep

#6 Post by scpeter@jncasr » Thu Oct 06, 2022 12:56 pm

https://drive.google.com/file/d/1_LBzfJ ... sp=sharing

plz try this, sry for the inconvenience

martin.schlipf
Global Moderator
Global Moderator
Posts: 458
Joined: Fri Nov 08, 2019 7:18 am

Re: huge dE value in ionic steep

#7 Post by martin.schlipf » Fri Oct 07, 2022 7:23 am

Sorry I still cannot access it. Now I get an error message
You need access
Request access, or switch to an account with access.

scpeter@jncasr
Newbie
Newbie
Posts: 9
Joined: Tue Apr 12, 2022 9:13 am

Re: huge dE value in ionic steep

#8 Post by scpeter@jncasr » Fri Oct 07, 2022 1:25 pm


martin.schlipf
Global Moderator
Global Moderator
Posts: 458
Joined: Fri Nov 08, 2019 7:18 am

Re: huge dE value in ionic steep

#9 Post by martin.schlipf » Mon Oct 10, 2022 9:52 am

Thank you for the input and output files.

It seems that something is not working as intended with the dipole correction. If you look for the term dipol+quadrupol energy correction in your OUTCAR file, you will see that it wildly oscillates. This in turn causes the jumps in the total energy. It appears that the dipole correction is not accounted for when judging whether the calculation should stop. I will inquire whether that is by design or by mistake.

One thing you could try is to fix the dipole center by hand. You can do this with the DIPOL tag. Please choose a position in the vacuum region and check whether that solves your issue.

martin.schlipf
Global Moderator
Global Moderator
Posts: 458
Joined: Fri Nov 08, 2019 7:18 am

Re: huge dE value in ionic steep

#10 Post by martin.schlipf » Mon Oct 10, 2022 1:45 pm

I got some more input on this issue: With LDIPOL = F, only the total energy is corrected not the potential. You should also verify that your system is converged with respect to the vacuum size.

If fixing the position of DIPOL does not help to solve the issue, converge the calculation without dipole corrections and try to include them only in the last step.

scpeter@jncasr
Newbie
Newbie
Posts: 9
Joined: Tue Apr 12, 2022 9:13 am

Re: huge dE value in ionic steep

#11 Post by scpeter@jncasr » Mon Oct 10, 2022 3:05 pm

Thank you for the suggestion, I will surely try it. As you said with LDIPOL = F the total energy is corrected so my question is if the energy value and the structural parameters are reliable without the LDIPOL tag.

Locked