19.01.2015 Views

MOLPRO

MOLPRO

MOLPRO

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

43 VIBRATIONAL FREQUENCIES (FREQUENCIES) 334<br />

43 VIBRATIONAL FREQUENCIES (FREQUENCIES)<br />

FREQUENCIES,options,<br />

Calculate harmonic vibrational frequencies and normal modes. The hessian is calculated analytically<br />

or numerically by finite differences in 3N cartesian coordinates (Z-Matrix coordinates<br />

will be destroyed on entry). If analytic gradients are available these are differentiated once to<br />

build the hessian, otherwise the energy is differentiated twice. If for the wavefunction method<br />

dipole moments are available, the dipole derivatives and the IR intensities are also calculated.<br />

Note that numerical hessians cannot be computed when dummy atoms holding basis functions<br />

are present. To get reasonable results it is necessary to do a geometry optimization before using<br />

the frequency calculation.<br />

43.1 Options<br />

The following options are available:<br />

ANALYTICAL<br />

CENTRAL<br />

NUMERICAL<br />

FORWARD<br />

SYMM=AUTO|NO<br />

AUTO<br />

NOAUTO|NOSYM<br />

Use analytical second derivatives of the energy. At present, analytical<br />

second derivatives are only possible for closed shell Hartree-Fock<br />

(HF) and MCSCF wavefunctions without symmetry. It is not yet possible<br />

to calculate IR-intensities analytically. Note that, due to technical<br />

reasons, the analytical MCSCF second derivatives have to be computed<br />

in the MCSCF-program using e.g. multi; cpmcscf,hess<br />

(see MULTI) before they can be used in FREQUENCIES. If analytical<br />

MCSCF second derivatives have been computed using multi;<br />

cpmcscf,hess, FREQUENCIES will use them by default.<br />

Use central differences/high quality force constants (default).<br />

Differentiate the energy twice, using central differences.<br />

Use forward differences/low quality force constants.<br />

During the numerical calculation of the hessian, the symmetry of the<br />

molecule may be lowered. Giving SYMM=AUTO the program uses<br />

the maximum possible symmetry of the molecular wavefunction in<br />

each energy/gradient calculation, and this option therefore minimizes<br />

the computational effort. With SYMM=NO no symmetry is used during<br />

the frequency calculation (default). For single reference calculations<br />

like HF, MP2, CCSD, RCCSD the AUTO option can be safely<br />

used and is recommended. However, the AUTO option cannot be used<br />

for multireference methods (MCSCF/MRCI/ACPF/AQCC/RS2). If<br />

given, the option is disabled in these cases. For these methods frequency<br />

calculations are only possible without symmetry. Symmetry<br />

is turned off atomatically if the state symmetry is 1. Note that this may<br />

fail if there are lower states in other symmetries. Use of RESTRICT,<br />

SELECT, REF, PROJECT, LOCAL, state-averaged MCSCF will lead<br />

on errors unless the calculation is performed in C 1 symmetry In such<br />

cases the whole calculation must be done without symmetry.<br />

Same as SYMM=AUTO, see above.<br />

Same as SYMM=NO, see above.<br />

HESSREC|SAVE=record Save hessian to record. By default the hessian is saved on record<br />

5300.2.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!