Input Block 'String'


Click here for full thread
Gets Around
You can also go back and forth between neb and string calculations, e.g.

This run never really converges well because the two H's that are doing all the real work pop off to form H2 and do a little dance and then come back to the CH2 group they left behind. Optimizing the details of the little dance is about as likely has finding blinker fluid at the auto parts store

Note that this sends H1 right through the Carbon - atom avoidance code must work


echo
start ch4_neb

permanent_dir ./permh
scratch_dir ./permh

geometry nocenter
 symmetry c1
C 0.000000 0.000000 0.000000
H 0.000000 0.000000 1.089000
H 1.026719 0.000000 -0.363000
H -0.513360 -0.889165 -0.363000
H -0.513360 0.889165 -0.363000
end

geometry endgeom nocenter
 symmetry c1
C 0.000000 0.000000 0.000000
H 0.000000 0.000000 -1.089000
H 1.026719 0.000000 0.363000
H -0.513360 -0.889165 0.363000
H -0.513360 0.889165 0.363000
end

basis
 * library STO-3G
end

scf
  thresh 1.0d-10
maxiter 5001
end

neb
  impose
nhist 10
nbeads 10
maxiter 100
stepsize 0.10d0
kbeads 0.1
end
task scf neb ignore

neb
  kbeads 1.0
nbeads 20
end
task scf neb ignore

string
 nbeads 20
nhist 10
maxiter 100
stepsize 1.0
freeze1 .true.
freezeN .true.
interpol 1
end

task scf string ignore