Dimer method using IOPT 1 breaks the adsorbate molecule

Vasp transition state theory tools

Moderator: moderators

Post Reply
sxd375
Posts: 21
Joined: Tue Feb 13, 2018 3:30 pm
Location: State College

Dimer method using IOPT 1 breaks the adsorbate molecule

Post by sxd375 » Sun May 06, 2018 5:01 pm

Dr. Graeme

I have been running dimer for a transition state with LBFGS optimiser, with default values other than INVCURV = 0.005 (I read in the forum, it makes the optimiser rather conservative in taking steps). However, the adsorbate molecule breaks off and drifts apart in some 200th iteration leading to very high forces. I donot understand why. I have attached the dimer calculation here.

Please let me know what I am doing wrong in here.

P.S. also as I pointed out about LBFGS optimiser during dimer, the DIMCAR doesnot update/rotate for the second time even after 300 steps and as ever I get following DIMCAR(for other cases even though similar happened, but dimer converged still.)

Step Force Torque Energy Curvature Angle
1 0.18801 1.99353 -802.98099 -4.85432 1.72500
1 0.18801 1.15474 -802.98099 -4.95597 3.34210
1 0.18801 0.95404 -802.98099 -4.98347 1.63867

Thanks
Attachments
dim.rar
(17.17 MiB) Downloaded 216 times

Post Reply