CI-NEB: nearly there but not quite converged
Jörg Saßmannshausen
j.sassma... at ucl.ac.uk
Sun Apr 30 21:58:16 UTC 2017
Dear all,
I am battling a problem with some CI-NEB calculations I am doing.
They are nearly converged but not quite:
RMS DISPLACEMENT = 0.00010 [ 0.00010] (YES)
MAX DISPLACEMENT = 0.00051 [ 0.00020] ( NO)
RMS FORCE = 0.00019 [ 0.00050] (YES)
MAX FORCE = 0.00102 [ 0.00100] ( NO)
If I do a frequency calculation on this geometry for the transition state I usually find the imaginary frequency of the
transition state but also one or two more which corespond to say a rotation of a Me group or a rotation of a Cp ring of
my metallocene.
Now, I could say: I ignore it but I really would like to get that done 'properly'.
I am already using quite tight parameters, for example:
EPS_SCF 1.0E-8
EPS_DEFAULT 1.0E-16
I am doing a PBE-D3/DZVP level of theory calculations.
Also, if I am selecting such a step for a restart, it seems to jump off too far:
RMS DISPLACEMENT = 0.00130 [ 0.00010] ( NO)
MAX DISPLACEMENT = 0.00893 [ 0.00020] ( NO)
RMS FORCE = 0.00253 [ 0.00050] ( NO)
MAX FORCE = 0.01721 [ 0.00100] ( NO)
These are the parameters I am using:
BAND_TYPE CI-NEB
K_SPRING 0.02
ROTATE_FRAMES T
&CI_NEB
NSTEPS_IT 5
&END
and
&OPTIMIZE_BAND
OPTIMIZE_END_POINTS T
OPT_TYPE DIIS
&DIIS
MAX_STEPS 600
N_DIIS 7
NO_LS
STEPSIZE 0.5
MAX_STEPSIZE 1.0
&END
I am happy to post the input,the force_eval and subsys files if required.
Does anybody have some good advice here?
All the best from London
Jörg
--
*************************************************************
Dr. Jörg Saßmannshausen, MRSC
University College London
Department of Chemistry
20 Gordon Street
London
WC1H 0AJ
email: j.sassma... at ucl.ac.uk
web: http://sassy.formativ.net
Please avoid sending me Word or PowerPoint attachments.
See http://www.gnu.org/philosophy/no-word-attachments.html
-------------- next part --------------
Q29udGVudC1UeXBlOiBhcHBsaWNhdGlvbi9wZ3Atc2lnbmF0dXJlOyBuYW1lPXNpZ25hdHVyZS5h
c2MgDQpDb250ZW50LURlc2NyaXB0aW9uOiBUaGlzIGlzIGEgZGlnaXRhbGx5IHNpZ25lZCBtZXNz
YWdlIHBhcnQuDQoNCi0tLS0tQkVHSU4gUEdQIFNJR05BVFVSRS0tLS0tDQpWZXJzaW9uOiBHbnVQ
RyB2MS40LjEyIChHTlUvTGludXgpDQoNCmlGNEVBQkVJQUFZRkFsa0dYZmdBQ2drUVRDbi9MeTlT
L3djR3JBRC9SZXFjZzN4RDdBWVhCWWl6aDNVM1ZBYzcNCmdiWlIyTmhqUm5uV0FselpRUVVBLzB5
RmViUDRMUExaeTd1ZjBFaDZOWUsyL3E1NS9NUmJhU002a091WHpaS0cNCj1Xd1p1DQotLS0tLUVO
RCBQR1AgU0lHTkFUVVJFLS0tLS0NCg==
More information about the CP2K-user
mailing list