# 3.3.16. ALASKA and SLAPAF — A Molecular Structure Optimization¶

One of the most powerful functions of ab initio calculations is geometry predictions. The minimum energy structure of a molecule for a given method and basis set is instructive especially when experiment is unable to determine the actual geometry. Molcas performs a geometry optimization with analytical gradients at the SCF or RASSCF level of calculation, and with numerical gradients at the CASPT2 level.

In order to perform geometry optimization an input file must contain a loop, which includes several calls: calculation of integrals (SEWARD), calculation of energy (SCF, RASSCF, CASPT2), calculation of gradients (ALASKA), and calculation of the new geometry (SLAPAF).

This is an example of such input

&GATEWAY
coord= file.xyz
basis= ANO-S-MB
>> EXPORT MOLCAS_MAXITER=25
>> Do While <<
&SEWARD
&SCF
&SLAPAF
>> EndDo <<


The initial coordinates will be taken from xyz file file.xyz, and the geometry will be optimized at the SCF level in this case. After the wave function calculation, calculation of gradients is required, although code ALASKA is automatically called by Molcas. SLAPAF in this case required the calculation of an energy minimum (no input). Other options are transition states (TS), minimum energy paths (MEP-search), etc The loop will be terminated if the geometry converges, or maximum number of iterations (MOLCAS_MAXITER) will be reached (the default value is 50).

There are several EMIL commands (see Section 4.1.3.2) which can be used to control geometry optimization. For example, it is possible to execute some Molcas modules only once:

>> IF ( ITER = 1 )
* this part of the input will be executed only during the first iteration
>> ENDIF


Program SLAPAF is tailored to use analytical or numerical gradients produced by ALASKA to relax the geometry of a molecule towards an energy minimum (default, no input required then) or a transition state. The program is also used for finding inter state crossings (ISC), conical intersections (CI), to compute reaction paths, intrinsic reaction coordinate (IRC) paths, etc.

## 3.3.16.1. SLAPAF — Basic and Most Common Keywords¶

TS

Computing a transition state

FindTS

Computing a transition state with a constraint

MEP-search

Computing a steepest-descent minimum reaction path

ITER

Number of iterations

INTErnal

Definition of the internal coordinates