Overlay 7 IOPS
Last Update 6/25/2001
5 6 7 8 9 10 11 12 13 14 15 16 18 25 28 29 30 31 32 40 41 42 43 44 45 52 60 63 70
Overlay 7
Overlay 7 consists of programs used in the calculation of first and second derivatives of energy with respect to nuclear coordinates.
IOp(5)
Operation of Link 707. 0 Compute both Born-Oppenheimer corrections and one-electron integral second derivative terms. 1 Compute only Born-Oppenheimer corrections.
IOp(6)
The operation of link 705 (NYI). 0 Default (12). 1 Do not form the ecp contribution to the SCF forces. 2 Form the ecp contribution to the SCF forces. 10 Do not form 1e derivative matrices. 20 Increment the 1e derivative matrices with ecp terms.
IOp(7)
Use of internal coordinates. 0 Yes. 1 No. 2 Yes, but neglect first derivatives in conversion of second derivatives to internal coordinates.
IOp(8)
Harmonic frequency calculation. 0 No. 1 Yes, with most common isotopes. 2 Yes, with read-in isotopes. 10 Print higher precision normal modes. Nxx Default scale factor is #N (1=HF, 1/1.12, (2=CBS4=0.91671, 3=CBSQ=0.91844). I. If M=1, only harmonic thermochemistry. I. If M=2, do hindered rotor analysis. II. If M=3, read hindered rotor parameters from input.
IOp(9)
Whether to rotate derivatives back to the Z-matrix orientation: 0/1 Yes/No.
IOp(10)
First/second derivative control. 0 Do only first derivatives. 1 Do only second derivatives. 2 Do both.
IOp(11)
Control of integral derivative algorithm: 0 Default. Use IsAlg to decide. 2 Scalar Rys SPDF. 3 Berny SP, Scalar Rys DF. 4 Old vector Rys SPDF. 5 Berny SP, old vector Rys DF. 6 FoFDir: Rys spdf. 7 Berny SP, FoFDir Rys df. 8 FoFDir: HGP sp, Rys df. 9 Berny SP, FoFDir Rys df (same as 7). 10 FoFDir: HGP spd, Rys f. 11 Berny SP, FoFDir HGP d Rys f. 12 FoFDir: HGP spdf. 13 Berny SP, FoFDir HGP df. 14 FoFDir: PRISM spdf. 15 FoFDir: Berny SP, PRISM df.
IOp(12)
Selection of density matrix. 0 Usual SCF density. N Use generalized density number N for both the one-electron integral derivatives and the corresponding 2PDM terms.
IOp(13)
Contraction with two-particle density matrices: 0 Default (same as 1). 1 Use HF 2PDM. 2 Use external 2PDM. 3 Use both HF and external 2PDM. 4 Generate 2PDM from CIS square 1PDM (for debugging). 5 Generate 2PDM from CIS square 1PDM and use HF/z 2PDM as well. 6 Contract with external 2PDM derivatives. The types of derivatives are given by IOp(15). 7 Form derivative 2PDM from CIS and HF derivative density matrices. The types of derivatives are given by IOp(15). 10 Leave the external 2PDM on the disk instead of deleting it. 0-5 imply use of the generalized density in L701, while 6-7 imply use of the generalized density derivatives in L701.
IOp(14)
The state for CIS gradients. Defaults to 1.
IOp(15)
The nature of the perturbation(s). 0 Default (1st order nuclear and electric field). IJK Nuclear Kth order, electric field Jth order, magnetic field Ith order. 1000 Generate simulated density derivatives. Only 1, 10, and 11 are valid in Overlay 7.
IOp(16)
Number of translations and rotations to remove during redundant coordinate transformations: -2 0. -1 Normal (6 or 5 for linear molecules). 0 Default, same as -1. N N.
IOp(18)
Derivative accuracy option: 0 Compute to 10(-8) accuracy. 1 Do as accurately as possible in L702. 2 Use the original 'Berny' values in L702. 10 Do as accurately as possible in L703. 20 Use sleazier cutoffs in L703. 100 Do as accurately as possible in l708. 200 Use sleazier cutoffs in l708.
IOp(25)
Type of derivatives available. 0 First. 1 Second. 2 Third. 10 Read derivatives from checkpoint file (in Z-matrix orientation).
IOp(28)
Skip option to defer integral evaluation to L703. 0 Compute as normal. 2 Do all gradient integrals in L703.
IOp(29)
Mode of use of l716. 0 Normal. 1 Normal + generate estimated initial force constants. 6 Nuclear repulsion only (useful for testing). 00 Default method for initial force constants.
IOp(30)
Use of symmetry in Overlay 7. 0 Use (subject to availability). 1 Do not use.
IOp(31)
Handling of forces contributions. 0 Just use the forces in IRWFX. 1 Compute HF forces from D2E file and increment both FX and FXYZ (non-O11 PSCF grad and HF freq). 00 Use FX in conversion of force constants to internal coordinates. (HF freq, PSCF Freq=Numer). 10 Use FXYZ in conversion of forces constants to internal coordinates (PSCF opt with HF 2nd derivative).
IOp(32)
Punch option. 0 None. 1 Punch energy in format D24.16, forces and lower triangular force constants in format 6F12.8. 2 Punch nuclear coordinate derivatives. Forces are punched in 3D20.12 format, one card per atom. Force constants and third derivatives are punched in 4E20.12 format in compressed form.
IOp(40)
Neglect of integrals (only option 1 works in Overlay 7): 0 Keep all integrals. 1 Neglect four center integrals. 2 Neglect three center two-electron integrals as well. 3 Neglect 2e integrals with diatomic differential overlap. 10 Neglect three center one-electron integrals. 20 Neglect 1e integrals with diatomic differential overlap. 30 Do only overlap and not other 1e integrals.
IOp(41)
NDDO flag. 0 Evaluate integrals correctly. 1 Apply NDDO approximation.
IOp(42)
1PDM: 0 Use SCF total density. N Use generalized density N.
IOp(43)
Second order simultaneous optimization flag. 0 .false. 1 .true (other second derivative options must also be set appropriately). 2 .true (debugging option: compute fifth order WG and GG terms in L715).
IOp(44)
Handling of an applied electric field. -1 Do not add electric field terms to forces. 0 Update forces for a uniform electric field. 1 Update forces for the self-consistent reaction field (SCRF) method. 2 Update forces for a uniform electric field, with forces done the usual way for CIS or MP2 second derivatives.
IOP(45)
Controlling the projection of the reaction path. 0 Do not project. The point is a stationary point. 1 Project the reaction path and compute 3N-7 frequencies.
IOp(52)
Whether ECP integrals should be done in L701 as usual. 0 Yes. 1 No.
IOp(60-62)
Override standard values of IRadAn, IRanWt, and IRanGd.
IOp(63)
Whether to do FMM. 0 Use global default. 1 Turn off FMM here regardless.
Iop(70)
Whether to allow cavity to move in PCM derivatives. 0 Default (No). 1 Yes. 2 No.