Overlay 1 IOPS:
Last Update 6/25/2001
5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 49 52 53 55 56 57 70 80 90
Overlay 1
Overlay 1 contains programs that read in geometry and optimization input, as well as control optimization and numerical frequency calculations.
IOp(5)
L103: Mode of optimization. 0 Find local minimum. 1 Find a saddle point. N Find a stationary point on the energy surface with N negative eigenvalues of the 2nd derivative matrix. L107: Mode of search: 0 Locate the maximum in the LST path. 1 Scan the LST path.
IOp(6)
L102, L103, L105, L107, L109, L112, L113, L114: Maximum number of steps (or
number of steps for an LST scan).
0 NStep = max(20,nvar+10) (L103, L112).
= min(20,nvar+10) (L102, L105, L109).
= min(40,nvar+20) (L113, L114).
N NStep = N
IOp(7)
L103, L105, L109, L112, L113, L114: Convergence on the first derivative and estimated displacement for the optimization RMS first derivative .lt. Confv, RMS est. Displacement .lt. Convx=4*convf. -1 ConvF = 1/600 Hartree/Bohr or Radian. 0 Convf = 0.0003 Hartree/Bohr or Radian. N Convf = 10 (-n). L116, L117: Convergence on electric field/charges. -1 Default value for optimizations: 10-7. 0 Default value for single-points: 10-5 in L116, 10-7 in L117. n 10-N.
IOp(8)
L103, L109, L112: Maximum step size allowed during opt.
0 DXMAXT = 0.1 Bohr or Radian (L103, estm or UnitFC).
= 0.3 Bohr or Radian (L103, read or CalcFC).
= 0.2 Bohr or Radian (L105).
= 0.3 Bohr or Radian (L113, L114).
N DXMAXT = 0.01 * N
L117: General control.
0 The type of basin to use to partition the density isosurface. The default
is 4.
1 GradVne.
2 GradRho.
3 Do not use basins, use only the center of nuclear charge.
4 Use interlocking spheres.
N0 Order of Adams-Bashforth-Moulton (ABM) predictor-corrector method to use
in solving different equations for the GradRho or Vne trajectories. Default is
4, maximum is 9.
N00 Number of small steps per ABM step to be used in starting ABM and when "slow
down" is needed in ABM. Default is 5.
N000 Which approximation to make. Default is III for Tomasi (interlocking
spheres) and IV for general surface.
1000 Approx. I—Do not do self-polarization or "compensation".
2000 Approx. II—Do self-polarization, but no "compensation".
3000 Approx. III—Do self-polarization and "compensation".
4000 Approx. IV—Do III and allow surface to "relax" in
solution if no spheres.
N0000 Whether to evaluate densities using orbitals or density matrix.
Default is to use density.
10000 Use MOs.
20000 Use density.
IOp(9)
L103: Use of trust radius. 0 Whether to update trust radius (DXMaxT, default Yes). Default is Yes for minima, No for TS. 1 No. 2 Yes. 00 Whether to scale or search the sphere when reducing the step size to the trust radius (default search for minima, scale for transition states). 10 Scale. 20 Search. L107: Whether to maintain symmetry along the search path. 0 Yes. 1 No. L117: Whether to delete points which are too close together: 0 No. 1 Yes, using a default criteria (0.05 Angstroms). -N Yes, using a (10-N Angstroms) criteria. How close to get to the isosurface in search: 0 Approx 1.0D-6 (N=20). N 2.0**-N.
IOp(10)
L103, L105, L109, L112, L113, L114: Input of initial Hessian: All values must be in atomic units (Hartree, Bohr, and Radians). 0 Use defaults. (Not valid for L109). 1 Read ((FC(I,J),J=1,I),I=1,NVAR) (8F10.6) (L103 only). 2 Read I,J,FC(I,J) (5I3,F20.0) (L103 only). End with a blank card. 3 Read from checkpoint file in internal coordinates. 4 Second derivative matrix calculated analytically. (Not valid for L109). 5 Read Cartesian forces and force constants from the checkpoint file are converted to internal coordinates. 6 Read Cartesian forces followed by Cartesian force constants (both in format 6F12.8) from input stream. Followed by a blank line. 7 Use semiempirical force constants. 8 Use unit matrix (default for L105; only recognized by 103). 9 Estimate force constants using valence force field.
IOp(11)
L103: Test of curvature. Bomb the job if the second derivative matrix has the wrong number of negative eigenvalues. 0 Default (test for Z-matrix or Cartesian TS but not for LST/QST or for minimum). 1 Don't test. 2 Test. L117: Scaling factor for determining overlaps of VDW atoms. -1 Turn off scaling. 0 Default is 1.010. N 1.000 + n*(0.001). Step size for ABM method in trudge for isodensity method. 0 0.05 (n=2). N 0.1/n.
IOp(12)
L103: Optimization control parameters. 0 Use default values. 1 Read in new values for all parameters (see InitBS).
IOp(13)
L103, L113, L114, L115: Type of Hessian update: 0 Default. (New D2Corr for L103 and L115, Powell for L113 and L114). 1 Powell. (Not in L103). 2 BFGS. (Not in L103). 3 BFGS, safeguarding positive definiteness. (Not in L103 or L115). 4 D2Corr. (New, only in L103 and L115). 5 D2Corr. (Old, only in L103 and L115). 6 D2Corr. (BFGS). 7 D2Corr. (Bofill Powell+MS for transition states). 8 D2Corr (No update, use initial Hessian).
IOp(14)
L103: Maximum number of bad steps to allow before doing a linear minimization (i.e., no quadratic step) will be attempted: 0 Default. (0 for TS, 1 for minima). N Allow N—linear only starts with the N+1st.
IOp(15)
L103, L109: Abort if derivatives too large. -1 or 0 No force test at all. N FMAXT = 0.1 * N.
IOp(16)
L103, L113, L114: Maximum allowable magnitude of the eigenvalues of the second derivative matrix. If the limit is exceeded, the size of the eigenvalue is reduced to the maximum, and processing continues. 0 EIGMAX = 25.0 Hartree/Bohr2 or Radian2. N EIGMAX = 0.1 * N.
IOp(17)
L103, L113, L114: Minimum allowable magnitude of the eigenvalues of the second derivative matrix. Similar to IOp(16). 0 EIGMIN = 0.0001. N EIGMIN = 1. / N.
IOp(18)
L103: Star only and coordinate option. 0 Proceed normally. 1 Second derivatives will be computed as directed on the variable definition cards. No optimization will occur. 10 Do optimization in Cartesian coordinates. 20 Do full optimization in redundant internal coordinates. 30 Do full optimization in pruned distance matrix coords. 40 Do optimization in Z-matrix coordinates. 50 Do full optimization in redundant internal coords with large molecular tools. 100 Read the AddRedundant input section for each structure. 1000 Do not define H-bonds 2000 Define H-bonds with no related coordinates (default) 3000 Define H-bonds and related coordinates 10000 Reduce the number of redundant internals 20000 Define all redundant internals 100000 Old definition of redundant internals. 1000000 Skip MM atoms in internal coordinate definitions. 2000000 Include MM atoms in internal coordinate definitions.
IOp(19)
L103: Search selection. 0 Default (same as 6). 2 Linear and steepest descent. 3 Steepest descent and linear only when essential. 4 Quadratic if curvature is correct; RFO if not. Linear as usual. 5 Quadratic if curvature is correct; RFO if not. No linear search. 6 RFO and linear. 7 RFO without linear. 8 Newton-Raphson and linear. 9 Newton-Raphson only. 10 GDIIS and linear. 11 GDIIS only. L113, L114: Search selection: 0 P-RFO or RFO step only (default). 1 P-RFO or RFO step for "wrong" Hessian otherwise Newton-Raphson.
IOp(20)
L101, L106, L108, L109, L110: Input units. 0 Angstroms degrees. 1 Bohrs degrees. 2 Angstroms radians. 3 Bohrs radians.
IOp(21)
L103, L113, L114: Expert switch. 0 Normal mode. 1 Expert mode. Certain cutoffs used to control the optimization will be relaxed. These include FMAXT, DXMAXT, EIGMAX and EIGMIN.
IOp(22)
L107: Whether to reorder coordinates for maximum coincidence. 0 Yes. 1 Assume reactant order equals product order. 2 Read in a re-ordering vector from the input. L115: Kind of search: 0 Both directions and generate search vector. 1 Forward direction and generate vector. 2 Backward direction and generate vector. 3 Both directions and generate vector. 4 Forward direction and read vector 8F10.6. 5 Forward direction and read vector 8F10.6. 6 Backward direction and read vector 8F10.6. 7 Both directions and read vector 8F10.6.
IOp(23)
L112: Derivative availability. 0 Energy only. 1 Energy + Forces. 2 Energy + Forces + Force constants.
IOp(24)
Whether to round tetrahedral angles. 0 Default. (Yes). 1 Yes, round angles within 0.001 degree. 2 No.
IOp(25)
Whether SCRF is used with numerical polarizability: 0 No. 1 Yes, the field in /Gen/ must be cleared each time.
IOp(26)
Accuracy of function being optimized:
- NNMM Energy 10-(NN), gradient 10-(MM).
- 1 Read in values.
0 Default (same as 1). 1 Normal accuracy for HF (energy and gradient both 1.D-7). 2 Standard grid accuracy for DFT (energy 1.D-7, gradient 1.D-5). 3 Fine grid accuracy for DFT (energy 1.D-7, gradient 1.D-6).
IOp(27)
= IJKL (i.e. 1000*I+100*J+10*K+L). Transition state searching using QST and redundant internal coordinates. L= 0,1 Input one structure, either initial guess of the minimizing structure or transition structure without QST. L= 2 Input 2 structures, the first one is the reactant, the second one is the product. The union of the two redundant coordinates is taken as the redundant coordinates for the TS. The values of the TS coordinates. are estimated by interpolating the structure of R and P. R and P are used to guide the QST optimization of the TS. L= 3 Input 3 structures. The first one is reactant the second one is the product. The third one is the initial guess of the transition structure. R and P are used to guide the QST optimization of the TS. K= 1-9 Interpolation of initial guess of TS between R and P (TS=0.1*J*R + 0.1*(10-J)*P, default J=5). J= 1 LST Constraint in internals. J= 2 QST Constraint in internals. J= 3 LST Constraint in distance matrix space. J= 4 QST Constraint in distance matrix space. I= 0-9 Control parameters for climbing phase of QST (e.g. QSTrad = 0.01*I, default QSTrad = 0.05).
IOp(28)
L101: Specification of nuclear centers. 0 By Z-matrix. L103: 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(29)
L101: Specification of nuclear centers. 0 By Z-matrix. 1 By direct coordinate input (must set IOp(29) in L202). 2 Get Z-matrix and variables from the checkpoint file. 3 Get Cartesian coordinates only from the checkpoint file. 4 By model builder, model A. 5 By model builder, model B. 6 Get Z-matrix from the checkpoint file, but read new values for some variables from the input stream. 7 Get all input (title, charge and multiplicity, structure) from the checkpoint file. 10 Print details of the model building process. 000 Default (same as 100). 100 Do not abort job if model builder generates a Z-matrix with too many variables. 200 Abort job if model builder generates a Z-matrix with too many variables. 1000 Read optimization flags in format 50L1 after the Z-matrix. 2000 Set all optimization flags to optimize. 3000 Purge flags except the frozen variables. 4000 Rebuild the coordinate system. 5000 (2+3) Purge all flags but keep the coordinate definition. 00000 Default, same as 10000. 10000 Mark Z-matrix constants as frozen variables rather than wired-in constants. 20000 Do not retain symbolic constants. 100000 Generate a symbolic z-matrix using all Cartesians if none is present on the checkpoint file (a hack to make IRCs work with Cartesian input).
IOp(30)
L103: This option is set for the last call to 103 in frequency calculations in order to preserve the values of the variables for archiving. It also suppresses error termination on large gradients. 0 Yes. 1 No.
IOp(32)
Title card punch control. 0 Do not punch. 1 Punch.
IOp(33)
L101, L102, L103, L106, L109, L110, L113, L114: Debug print. 0 Off. 1 On.
IOp(34)
L101, L102, L103: Debug + dump print. 0 Off. 1 On.
IOp(35)
Restart (L102-L112). 0 Normal optimization. 1 First point of a restart. Get geometry, wavefunction, etc. from the checkpoint file.
IOp(36)
Checkpoint. 0 Normal checkpoint of optimization. 1 Suppress checkpointing.
IOp(37)
D2E cleanup (obsolete). 0 No cleanup. 1 This is the last point at which analytic second derivatives will be done. Delete the D2E file and the buckets and truncate the read/write files.
IOp(38)
Entry control option (currently only by L106, L107, L108, L109, L110, L111, L112, and L115 but not L102, L103, and L105). 0 Continuation of run. 1 Initial entry. N>1 In L103: Initial entry of guided optimization using N levels. N0 In L106: differentiate nth derivatives once. In L110 and L111: differentiate energy N times. 000 In L106: differentiate wrt nuclear coordinates. 100 In L106: differentiate wrt electric field. 200 In L106: differentiate wrt field and nuclear.
IOp(39)
Step size control for numerical differentiation. (L106, L109, L110, L111). Path step size in L115. 0 Use internal default (0.005 Angstroms in L106 and L109, 0.01 Angstrom in L110, 0.001 au in L111). N Use step-size of 0.0001*N (Angstroms in L106, L109, L110, electric field au in L111). -1 Read stepsize in format D20.13.
IOp(40)
L113, L114: Hessian recalculation. -1 Pick up analytic second derivatives every time. 0 Just update. The default, except for CalcAll. N Recalculation the Hessian every N steps. L116: Whether to read initial E-field: 0 Start with 0.0. 1 Read from checkpoint file. 2 Read from input stream.
IOp(41)
Step number of optimization from which to take geometry. -1 for the initial geometry.
IOp(42)
L103, L115: Number of points along the reaction path in each direction. Default is 6. L117: Cutoff to be used in evaluating densities. 0 1.0D-10. N 1.0D-N.
IOp(43)
L116: Extent of reaction field. 0 Dipole. 1 Quadrupole. 2 Octapole. 3 Hexadecapole. L117: How to define radii. 0 Default is 11. 1 Use internally stored radii, centers will be on atoms. 2 Read-in centers and radii on cards. 10 Force Merz-Kollman radii (default). 20 Force CHELP (Francl) recommended radii. 30 Force CHELPG (Breneman) recommended radii. 100 Read in replacement radii for selected atom types as pairs (IAn,Rad) or (Symbol,Rad), terminated by a blank line. 200 Read in replacement radii for selected atoms as pairs (I,Rad), terminated by a blank line. Initial radius of spheres to be placed around attractors to "capture" the gradient trajectories. The final radius is then automatically optimized separately for each atom. 0 0.1. NM N.M = NM/10.
IOp(44)
IRC Type: 0 Default (same as 3). 1 Cartesian. 2 Internal. 3 Mass-weighted. L117: Maximum distance between a nucleus and its portion of the isosurface—used in Trudge only to eliminate, from the outset, points which clearly lie in another basin. This parameter should be chosen with the parameter Cont. in mind. 0 10.0 au. NM N.M au = NM/10.
IOp(45)
Read isotopes in L115. 0 Do not read isotopes. 1 Read isotopes.
IOp(46)
Order of multipoles in numerical SCRF: 0 Dipole. 1 Quadrupole. 2 Octapole. 3 Hexadecapole.
IOp(47)
Number of redundant internal coordinates to allow for. 0 Default: 50000. N N.
IOp(49)
Options to IRC path relaxation (IJKL). L 1/0 optimize reactant structure. K 1/0 optimize product structure. J 1/0 optimize TS structure (for QST input). I 1/0 bimolecular reaction.
IOp(52)
L101 and L120: Type of ONIOM calculation: 0/1 One layer, normal calculation. 2 Two layers. 3 Three layers. 00 No electrostatics included in the model systems. 10 Include electrostatics in model systems.
IOp(53)
L120: Action of each invocation of L120: 0 Do nothing 1 Create geometry of high level system. 2 Create geometry of medium level system. 3 Create geometry of low level system. 4 Integrate energy 5 Integrate energy and gradient 6 Integrate energy, gradient, and hessian N0 Save necessary information (some rwf's, energy, gradients, Hessian) of point N of the ONIOM grid. M00 Restore point M. 4--7--9* | | | 2--5--8 | | | 1--3--6
IOp(55)
L103: Options for GDIIS: ICos*1000+IChkC*100+IMix*10+Method form. L115: IRC optimization. 0 Default, use gradients to find the next geometry. 1 Use displacements to find the next geometry.
IOp(56)
Set of atom type names to parse: 0 Accept any. 1 Dreiding/UFF. 2 Amber.
IOp(57)
Whether to produce connectivity. 0 Default (2). 1 No 2 Yes, read from input stream 3 Yes, generate connectivity. 4 Yes, read from checkpoint file. 5 Yes, read from rwf file. 10 Read modifications.
IOp(70)
L118: Type of sampling (Nact) 0 Defalt (same as 3) 1 Orthant sampling 2 Microcanonical normal mode sampling 3 Fixed normal mode energy 4 Local mode sampling (now only Nact = 0 or 3 OK )
IOp(80)
.eq.1 to supress the 5th order correction after the surface hop has been made in Trajectory Surface Hopping calculations. Also needs IOp(10/80=1)
IOp(90)
Read in the velocity in Cartesian coordinates.