next up previous contents index
Next: Trajectories and Dynamics Up: Locating stationary points Previous: Transition states following a

Level-shifted mode-following

 

tex2html_wrap9437

The input needed for doing a level-shifted mode following  is very similar to the input for following a gradient extremal , and the keyword that is needed in order to invoke this kind of calculation is .MODFOL  . As for gradient extremals, we need to specify which mode we follow. However, a mode following does not use mass-weighted molecular coordinates  as default, and isotopic composition of the molecule is therefore not needed. Note, however, that mass-weighted coordinates can be requested through the keyword .MASSES   as described in the input section for the *WALK   module. A typical input following the third mode will thus look like:

**DALTON INPUT
.WALK
*WALK
.START
.MODFOL
.INDEX
 1
.MODE
 3
**WAVE FUNCTIONS
.HF
*END OF INPUT

The level-shifted mode-following uses an algorithm similar to the one used in the ordinary geometry optimization of a molecule, but whereas one in minimizations chooses a step so that the level shift parameter is less than the lowest eigenvalue of the molecular Hessian , this level shift parameter is chosen to be in-between the eigenvalues tex2html_wrap_inline9431 and tex2html_wrap_inline9433 if we are following mode number t. This approach is described was pionered by Cerjan and Miller [17], and is also described in Ref. [11]. As for the gradient extremal approach, higher-order transition states  can be requested throguh the use of the keyword .INDEX  .

Note that it may often be necessary to start the mode-following calculation by stepping out of the stationary point along the mode of interest using the keyword .EIGEN   in the .WALK   module. We refer to the reference manual for a further description of this option.

The index  of the critical point -- that is, the number of negative Hessian eigenvalues  -- sought, need to be specified (by the keyword .INDEX  ), as the calculation would otherwise continue until a critical point with index zero (corresponding to a minimum), is found.


next up previous contents index
Next: Trajectories and Dynamics Up: Locating stationary points Previous: Transition states following a

Kenneth Ruud
Sat Apr 5 10:26:29 MET DST 1997