MOLCAS manual:

Next: 6.31 mknemo Up: 6. Programs Previous: 6.29 mckinley


6.30 mclr

The mclr program in MOLCAS performs response calculations on single and multiconfigurational SCF wave functions with the technique described in [63]. The right hand side (RHS) and thus the perturbation has to be defined through a preceding MCKINLEY calculation. Second order derivatives are obtained from a MCKINLEY and a consecutive mclr calculation, with a geometrical displacement as the external perturbation. If the response of a geometrical perturbation is calculated, harmonic frequencies corresponding to the most abundant masses are printed. MCLR also calculates isotope shifted frequencies. Per default, vibrational frequencies are calculated for all possible single isotopic substitutions. MCLR always calculates the response of an electric field and prints the polarizability.

Note that the user should not normally, for frequency calculations, request the execution of this module since this will be automatic.

The mclr code also calculates the Lagrangian multipliers required for a SA-MCSCF single state gradient [64], where the RHS is generated by the program itself. Through an mclr and a consecutive alaska calculation, analytical gradients of a SA-MCSCF state may be obtained. Note that alaska will automatically run the mclr module! Thus, with structure geometry optimizations of excited MCSCF states can be performed (See input example below).

The MCLR program is based on the split GUGA formalism. However, it uses determinant based algorithms to solve the configuration interaction problem [65], in analogy to how it is done in the RASSCF. For spin symmetric wave function (MS=0) the time reversal symmetry is used, and the innermost loops are performed in combinations instead of determinant.

The upper limit to the size of the CI wave function that can be handled with the present program is about the same as for the RASSCF. The present version of the code is just able to handle CASSCF wave function, RASSCF wave function will soon be included.

The orbital handling is based on a one index transformation technique. The integrals is the transformed to occupied orbitals in two indexes, this can be done directly or from disk based integrals generated by seward.

6.30.1 Dependencies

To start the MCLR module the one-electron integrals generated by SEWARD have to be available. Moreover, MCLR requires the wave function from a SCF or RASSCF calculation and apart from in an evaluation of SA-MCSCF gradients, it also requires the differentiated integrals from MCKINLEY.

6.30.2 Files Input files

MCLR will use the following input files: ONEINT, ORDINT,RUNFILE, ABDATA, RYSRW, JOBIPH (for more information see [*]). Output files

MCKINTCommunication file between mclr and mckinley and rassi
UNSYMASCII file where all essential information, like geometry, Hessian normal modes and dipole derivatives are stored.
MLDNFQMolden input file for harmonic frequency analysis. Scratch files

TEMP0xx=1,8 used for for integral transformation and storing half transformed integrals.
REORDUsed for storing data used in the transformation of CI vectors from determinant base to CSF base.
TEMPCIVExchange file for temporary storing the CI vectors during the PCG.
RESPBinary bile where the solution of the response equations are stored.
JOPRUsed for half transformed integrals in direct mode.
KOPRUsed for half transformed integrals in direct mode.

6.30.3 Input

This section describes the input to the MCLR program in the MOLCAS program system. The input for each module is preceded by its name like:

  &MCLR Optional keywords

A list of these keywords is given below:
SALAMakes MCLR compute the Lagrangian multipliers for a state average MCSCF wave function. These multipliers are required by ALASKA to obtain analytical gradients for an excited state, when the excited state is determined by a SA optimization. SALA has to be followed by an integer on the next line, specifying the excited state for which the gradient is required. This integer specifies which root in the SA root the multipliers are calculated for.
EXPDimensionHere follows the dimension of the explicit Hamiltonian used as preconditioner in the Preconditioned conjugate gradient algorithm. Default 100.
ITERationsSpecify the maximum number of iterations in the PCG. Default 50.
LOWMemoryLowers the amount of memory used, by paging out the CI vectors on disk. This will lower the performance, but the program will need less memory.
PRINtRaise the print level, default 0.
RASSiThis keyword is used for transforming the CI vectors to split GUGA representation, and transforming the orbital rotations to AO basis, to make the response accessible for state interaction calculations.
SEWArdSpecify one particle operators, used as right hand side, form the ONEINT file constructed by SEWARD The keyword is followed by one row for each perturbation: LABEL symmetry Component
EndSewardMarks the end of perturbation specifications read from SEWARD ONEINT file.
THREsholdSpecify the convergence threshold for the PCG. Default is 1.0e-04.
DISOTOPECalculates frequencies modified for double isotopic substitution.
THERmochemistryRequest an user specified thermochemical analysis. The keyword must be followed by a line containing the Rotational Symmetry Number, a line containing the Pressure (in atm), and lines containing the Temperatures (in K) for which the thermochemistry will be calculated. The section is ended by the keyword "End of PT".

TIMECalculates the time dependent response of an electric periodic perturbation. The frequency of the perturbation should be specified on the following line. Used to calculated time dependent polarizabilities and required in a RASSI calculation of two photon transition moments.
MASSUsed to generate single and double (in conjunction with DISO) isotope shifted frequencies, with the isotope masses specified by the user. This implementation can be useful for example in calculating intermolecular frequencies which are contaminated by the BSSE. By setting the corresponding masses to the very large numbers, ghost orbitals can be used in the frequency calculation. MASS needs the atomic label and the new mass in units of u (real), for each element of the molecule. Input example

A default input for a harmonic frequency calculation.


An input for a harmonic frequency calculation with modified isotopic masses for hydrogen and oxygen.

H  =  2.0079
O  =  150000.998

Thermochemistry for an asymmetric top (Rotational Symmetry Number = 1), at 1.0 atm and 273.15, 298.15, 398.15 and 498.15 K.

  273.15  ;  298.15  ;  398.15  ;  498.15
End  of  PT

The time dependent response is calculated for a perturbation of frequency 0.2 au.

TIME  =  0.2
The time dependent response is calculated for a perturbation of frequency 0.2 au.

SALA  =  2
Computes the Lagrangian multipliers for state number 2 in the SA root. Note, that 2 refers to the SA root. Thus, if the ground state is not included in the SA, the numbering of roots in the ci root and SA root differ. With the following RASSCF input

  2  3
  2  3
  1  1
RlxRoot  =  2

SALA 2 yields the gradient for ci root number 3. Geometry optimization of an excited SA-CASSCF state can be done normally using AUTO, and requires the use of the RLXR keyword in the RASSCF input to specify the selected root to be optimized. An explicit input to MCLR is not required but can be specified if default options are not appropriate.

next up previous contents
Next: 6.31 mknemo Up: 6. Programs Previous: 6.29 mckinley