Example: (non-)Linear Transit: H2O¶
In ADF2008.01 a transit calculation option has been added in the new optimization branch. This is capable of performing both linear transits, and non-linear transits, and is the default when the LINEARTRANSIT or TRANSIT sub-block is included in the ‘Geometry’ block.
The new transit code works differently to the old: the transit is represented as a sequence of constrained optimizations. A ‘Constraints’ block is used to delineate the constraints applied at each stage of the transit.
Non-linear transits are possible, and can even be combined with linear transits in other coordinates. To perform a non-linear transit in a particular coordinate, explicit values must be given.
$ADFBIN/adf << eor Title WATER Transit (non-linear), with the new optimizer branch Atoms O 0.000000 0.000000 0.000000 H 0.000000 -0.689440 -0.578509 H 0.000000 0.689440 -0.578509 End Symmetry NOSYM Constraints dist 1 2 0.8 1.0 1.25 1.5 angle 2 1 3 start=100.0 end=120.0 End Basis Type SZ Core Large End Geometry Transit 4 Optim Deloc Converge 0.0001 End End Input
In the example above, 4 values are given for the distance between atoms 1 and 2. This distance constraint will be applied simultaneously with the linear transit constraints for the angle, with other degrees of freedom optimized at each stage of the transit.
Finally, it should be pointed out that ‘partial constraints’ are used by default in the transit calculations. These constraints are not required to be fully met at each intermediate geometry, but are fully met at the converged geometries. You can use fully converged constraints by supplying the FULLCONVERGE option to the ‘Constraints’ subblock of the ‘Geometry’ block (not to be confused with the ‘Constraints’ block at root level).