21

Two-Dimensional Grid Search

A two-dimensional grid search calculates the properties of geometries corresponding to a set of points in a grid of user-defined size in two dimensions. In many ways, it is analogous to the reaction coordinate method described above, although two variables are changed and not as much information is provided. The final result is an 11×11 matrix of energies (121 in all) corresponding to sites on the grid. A grid search is called for if the keywords STEP1=n.n and STEP2=m.m are on the keyword line and two coordinates are marked with "-1" as the optimization flags. They are related to STEP1 and STEP2 in order of their definition. The reaction described in this file is part of the surface describing the inversion of 1,5-hexadiene:

Figure 21.1. Inversion of 1,5-hexadiene

Inversion of 1,5-hexadiene

Input File (rxnpath/rxngrid1.dat) :

  am1 rhf singlet qcscf t=auto step1=0.15 step2=0.15 grad  1
1,5-HEXADIENE INVERSION
Rxn grid calculation
 C              0.000000  0    0.000000  0    0.000000  0    0    0    0
 C              1.300000  1    0.000000  0    0.000000  0    1    0    0
 C              1.540000  1  130.000000  1    0.000000  0    2    1    0
 C              1.540000 -1  120.000000  1    0.000000  1    3    2    1   2
 C              1.540000  1  120.000000  1    0.000000  1    4    3    2
 C              2.190000 -1  110.000000  1    0.000000  1    1    2    3   3
 H              1.000000  1  120.000000  1   90.000000  1    1    2    3
 H              1.000000  1  120.000000  1  180.000000  1    1    2    7
 H              1.000000  1  120.000000  1 -180.000000  1    2    3    4
 H              1.000000  1  109.500000  1  120.000000  1    3    4    5
 H              1.000000  1  109.500000  1 -120.000000  1    3    4    5
 H              1.000000  1  109.500000  1  120.000000  1    4    3    2
 H              1.000000  1  109.500000  1 -120.000000  1    4    3    2
 H              1.000000  1  120.000000  1  120.000000  1    5    4    3
 H              1.000000  1  120.000000  1   90.000000  1    6    5    4
 H              1.000000  1  120.000000  1  180.000000  1    6    5   15
 0              0.000000  0    0.000000  0    0.000000  0    0    0    0
1

A reaction grid calculation will be performed if exactly two of optimization flags are marked with a "-1". The keywords STEP1 and STEP2 are also required for a reaction grid. By default, TRUSTE will be use to perform the constrained optimization at each of the grid points but this default can be overridden by BFGS, DFP, or EF.

2

This is the first reaction coordinate and will be incremented and decremented by the value of STEP1.

3

This is the second reaction coordinate and will be incremented and decremented by the value of STEP2. These two items will be used to construct a grid of points and optimizations will be carried out at each point on the surface. The point given in the geometry definition is the center of the grid and five steps are taken in each direction from this center point. The following list shows the values that the first and second reaction coordinates will take on. The center point is indicated by ×:

                                       Step 1

	      0.79  0.94  1.09  1.24  1.39  1.54  1.69  1.84  1.99  2.14  2.29
	1.44   o     o     o     o     o     o     o     o     o     o     o
	1.59   o     o     o     o     o     o     o     o     o     o     o
	1.74   o     o     o     o     o     o     o     o     o     o     o
S	1.89   o     o     o     o     o     o     o     o     o     o     o
t	2.04   o     o     o     o     o     o     o     o     o     o     o
e	2.19   o     o     o     o     o     ×     o     o     o     o     o
p	2.34   o     o     o     o     o     o     o     o     o     o     o
	2.49   o     o     o     o     o     o     o     o     o     o     o
2	2.64   o     o     o     o     o     o     o     o     o     o     o
	2.79   o     o     o     o     o     o     o     o     o     o     o
	2.94   o     o     o     o     o     o     o     o     o     o     o

Archive File (rxnpath/rxngrid1.arc) :

This file is abbreviated to conserve space. The archive file contains the geometries corresponding to each grid point (in the order they are computed) and so contains 121 geometries. Only the first and last geometry are shown here separated by "=========" to indicate the missing geometry. A complete copy of the file is found in the test suite results directory.

 Timestamp: 2004-02-12-14-45-20-0000008014-hpux

                     SUMMARY OF   AM1   CALCULATION
                                                       Feb-12-2004
                          AMPAC Version 8.13
                             Presented by:

                        Semichem, Inc.
                        PO Box 1649
                        Shawnee KS 66222
                        (913)268-3271
                        (913)268-3445 (fax)

 FORMULA: C6H10
 1,5-HEXADIENE INVERSION
 Rxn grid calculation

     GEOMETRY OPTIMISED : ENERGY MINIMISED
     SCF FIELD WAS ACHIEVED

          FINAL HEAT OF FORMATION =      35.206020 kcal
                                  =     147.337193 kJ
          ELECTRONIC ENERGY       =   -3921.652079 eV
          CORE-CORE REPULSION     =    3017.264491 eV
          TOTAL ENERGY            =    -904.387588 eV
          GRADIENT NORM           =       0.157153
          RMS GRADIENT NORM       =       0.024848
          UNSTABLE MODE(S)        =       0 ( ESTIMATE  )
          DIPOLE                  =       0.117254 debyes
          NO. OF FILLED LEVELS    =      17 (OCC = 2)
          KOOPMAN IONISATION POTENTIAL =      8.82 eV
          MOLECULAR POINT GROUP   = C2    0.100000
          MOLECULAR WEIGHT        =      82.145
          COMPUTATION TIME        =       0.02     seconds

          FINAL GEOMETRY OBTAINED                                 CHARGE  1
 AM1 RHF SINGLET QCSCF T=AUTO STEP1=0.15 STEP2=0.15 GRAD
 1,5-HEXADIENE INVERSION
 Rxn grid calculation
  C     0.000000  0    0.000000  0    0.000000  0    0   0   0   -0.1854
  C     1.344434  1    0.000000  0    0.000000  0    1   0   0   -0.1969
  C     1.478627  1  121.661128  1    0.000000  0    2   1   0   -0.1251
  C     1.540000 -1  109.338471  1  -68.221670  1    3   2   1   -0.1251
  C     1.478644  1  109.342139  1   48.195073  1    4   3   2   -0.1970
  C     2.190000 -1   95.333280  1   65.292648  1    1   2   3   -0.1854
  H     1.099488  1  121.392710  1  172.273579  1    1   2   3    0.1054
  H     1.100141  1  123.026296  1  166.927675  1    1   2   7    0.1136
  H     1.100202  1  116.414007  1 -257.849915  1    2   3   4    0.1190
  H     1.119799  1  109.032404  1  169.890197  1    3   4   5    0.0840
  H     1.121482  1  108.702828  1  -72.950161  1    3   4   5    0.0854
  H     1.119844  1  109.026704  1  169.881233  1    4   3   2    0.0840
  H     1.121538  1  108.704126  1  -72.945162  1    4   3   2    0.0854
  H     1.100293  1  116.407406  1  102.189776  1    5   4   3    0.1190
  H     1.099455  1  121.393780  1  172.194620  1    6   5   4    0.1054
  H     1.100119  1  123.047266  1  167.035349  1    6   5  15    0.1136
  0     0.000000  0    0.000000  0    0.000000  0    0   0   0

========

 Timestamp: 2004-02-12-14-45-20-0000008014-hpux

                     SUMMARY OF   AM1   CALCULATION
                                                       Feb-12-2004
                          AMPAC Version 8.13
                             Presented by:

                        Semichem, Inc.
                        PO Box 1649
                        Shawnee KS 66222
                        (913)268-3271
                        (913)268-3445 (fax)

 FORMULA: C6H10
 1,5-HEXADIENE INVERSION
 Rxn grid calculation

     GEOMETRY OPTIMISED : ENERGY MINIMISED
     SCF FIELD WAS ACHIEVED

          FINAL HEAT OF FORMATION =     154.169060 kcal
                                  =     645.197515 kJ
          ELECTRONIC ENERGY       =   -3845.725386 eV
          CORE-CORE REPULSION     =    2946.496423 eV
          TOTAL ENERGY            =    -899.228963 eV
          GRADIENT NORM           =       0.098920
          RMS GRADIENT NORM       =       0.015641
          UNSTABLE MODE(S)        =       0 ( ESTIMATE  )
          DIPOLE                  =       1.239565 debyes
          NO. OF FILLED LEVELS    =      17 (OCC = 2)
          KOOPMAN IONISATION POTENTIAL =      7.72 eV
          MOLECULAR POINT GROUP   = C2    0.100000
          MOLECULAR WEIGHT        =      82.145
          COMPUTATION TIME        =       9.34     seconds

          FINAL GEOMETRY OBTAINED                                 CHARGE  2
 AM1 RHF SINGLET QCSCF T=AUTO STEP1=0.15 STEP2=0.15 GRAD
 1,5-HEXADIENE INVERSION
 Rxn grid calculation
  C     0.000000  0    0.000000  0    0.000000  0    0   0   0   -0.1350
  C     1.525687  1    0.000000  0    0.000000  0    1   0   0   -0.2051
  C     1.450381  1  113.086785  1    0.000000  0    2   1   0   -0.0847
  C     2.290000 -1  106.042680  1  -40.526018  1    3   2   1   -0.0851
  C     1.450485  1  106.022529  1   33.582836  1    4   3   2   -0.2050
  C     1.440000 -1  114.358636  1   66.035586  1    1   2   3   -0.1350
  H     1.128058  1  107.146059  1  188.031302  1    1   2   3    0.0688
  H     1.126727  1  108.026342  1  114.107632  1    1   2   7    0.0856
  H     1.127706  1  111.933656  1 -276.885012  1    2   3   4    0.0907
  H     1.093766  1  111.697633  1  164.356002  1    3   4   5    0.0657
  H     2.084619  1  134.840888  1   19.349081  1    3   4   5    0.1144
  H     1.093772  1  111.532113  1  164.210983  1    4   3   2    0.0656
  H     2.084393  1  134.831308  1   19.362668  1    4   3   2    0.1144
  H     1.127689  1  111.964925  1   83.091637  1    5   4   3    0.0902
  H     1.128056  1  107.149290  1  188.076712  1    6   5   4    0.0688
  H     1.126735  1  108.023720  1  114.110585  1    6   5  15    0.0857
  0     0.000000  0    0.000000  0    0.000000  0    0   0   0
1

This is the first geometry in the archive file and corresponds to the central geometry of the grid.

2

This is the last geometry in the archive file.

Ouput File (rxnpath/rxngrid1.out) :

 Timestamp: 2004-02-12-14-45-20-0000008014-hpux
 *******************************************************************************
                             AM1 CALCULATION RESULTS
 *******************************************************************************
 *                             AMPAC Version 8.13
 *                                Presented by:
 *
 *                           Semichem, Inc.
 *                           PO Box 1649
 *                           Shawnee KS 66222
 *                           (913)268-3271
 *                           (913)268-3445 (fax)
 *
 *  STEP1    - 2-D GRID CALCULATION TO BE PERFORMED
 *  T=AUTO   - AUTOMATIC DETERMINATION OF ALLOWED TIME
 *  GRADIENTS- ALL GRADIENTS TO BE PRINTED
 *  SINGLET  - IS THE REQUIRED SPIN MULTIPLICITY
 *  AM1      - THE AM1 HAMILTONIAN TO BE USED
 *  QCSCF    - USE QUADRATICALLY CONVERGENT SCF METHOD
 *******************************************************************************
 AM1 RHF SINGLET QCSCF T=AUTO STEP1=0.15 STEP2=0.15 GRAD
 1,5-HEXADIENE INVERSION
 Rxn grid calculation
    ATOM    CHEMICAL   BOND LENGTH    BOND ANGLE    TWIST ANGLE
   NUMBER   SYMBOL     (ANGSTROMS)     (DEGREES)     (DEGREES)
    (I)                   NA:I          NB:NA:I      NC:NB:NA:I   NA  NB  NC
      1     C
      2     C          1.30000 *                                   1
      3     C          1.54000 *      130.00000 *                  2   1
      4     C          1.54000        120.00000 *     0.00000 *    3   2   1
      5     C          1.54000 *      120.00000 *     0.00000 *    4   3   2
      6     C          2.19000        110.00000 *     0.00000 *    1   2   3
      7     H          1.00000 *      120.00000 *    90.00000 *    1   2   3
      8     H          1.00000 *      120.00000 *   180.00000 *    1   2   7
      9     H          1.00000 *      120.00000 *  -180.00000 *    2   3   4
     10     H          1.00000 *      109.50000 *   120.00000 *    3   4   5
     11     H          1.00000 *      109.50000 *  -120.00000 *    3   4   5
     12     H          1.00000 *      109.50000 *   120.00000 *    4   3   2
     13     H          1.00000 *      109.50000 *  -120.00000 *    4   3   2
     14     H          1.00000 *      120.00000 *   120.00000 *    5   4   3
     15     H          1.00000 *      120.00000 *    90.00000 *    6   5   4
     16     H          1.00000 *      120.00000 *   180.00000 *    6   5  15

   MOLECULAR POINT GROUP            SYMMETRY CRITERIA
            C1                          0.10000000

          SINGLET STATE CALCULATION

        **  REFERENCES TO PARAMETERS  **

  H  (AM1): M.J.S. DEWAR ET AL, J. AM. CHEM. SOC. 107 3902-3909 (1985)
  C  (AM1): M.J.S. DEWAR ET AL, J. AM. CHEM. SOC. 107 3902-3909 (1985)

          CARTESIAN COORDINATES
    NO.       ATOM         X         Y         Z
     1         6        0.0000    0.0000    0.0000
     2         6        1.3000    0.0000    0.0000
     3         6        2.2899    1.1797    0.0000
     4         6        1.7632    2.6268    0.0000
     5         6        0.2466    2.8943    0.0000
     6         6       -0.7490    2.0579    0.0000
     7         1       -0.5000    0.0000   -0.8660
     8         1       -0.5000   -0.0000    0.8660
     9         1        1.6420   -0.9397    0.0000
    10         1        2.8470    1.0272   -0.8164
    11         1        2.8470    1.0272    0.8164
    12         1        2.0919    3.1017    0.8164
    13         1        2.0919    3.1017   -0.8164
    14         1       -0.1706    3.4075    0.7500
    15         1       -1.1319    1.7363    0.8660
    16         1       -1.1319    1.7363   -0.8660


 "GRID" PROCEDURE TO BUILD A 2-D SURFACE FROM A  42-D HYPERSURFACE.
 (VERSION 2, NOVEMBER 94).
 ATTENTION: SUCH A METHOD CAN YIELD SPURIOUS DISCONTINUITIES AS OPTIMIZED
 VARIABLES MAY RESPOND DISCONTINUOUSLY TO SMOOTH CHANGES IN THE DRIVING COORDS.
 THE PRESENT ALGORITHM IS TUNED TO INCREASE CHANCES OF GETTING A RELIABLE
 RESPONSE IN -AT LEAST- THE CENTRAL REGION OF THE GRID.

 STANDARD DEVIATION ON ENERGY   (KCAL)       0.00000055521
 STANDARD DEVIATION ON GRADIENT (KCAL/A,RD,RD)  0.00010624 0.00020591 0.00019285

 CENTRAL POINT IS LABELED (0,0)    1
   LABELS   1ST COORD  2ND COORD  HEAT OF FORM.   RMS-GRAD  TIME(sec)
 (  0,  0)      1.540      2.190    35.206020      0.157    0.47

 OPTIMIZED GEOMETRY OF CENTRAL POINT FOLLOWS.    2
 IT DEFINES THE "GERM" FOR SUBSEQUENT OPTIMIZATIONS.
    ATOM    CHEMICAL   BOND LENGTH    BOND ANGLE    TWIST ANGLE
   NUMBER   SYMBOL     (ANGSTROMS)     (DEGREES)     (DEGREES)
    (I)                   NA:I          NB:NA:I      NC:NB:NA:I   NA  NB  NC
      1     C
      2     C          1.34443 *                                   1
      3     C          1.47863 *      121.66113 *                  2   1
      4     C          1.54000        109.33847 *   -68.22167 *    3   2   1
      5     C          1.47864 *      109.34214 *    48.19507 *    4   3   2
      6     C          2.19000         95.33328 *    65.29265 *    1   2   3
      7     H          1.09949 *      121.39271 *   172.27358 *    1   2   3
      8     H          1.10014 *      123.02630 *   166.92768 *    1   2   7
      9     H          1.10020 *      116.41401 *  -257.84992 *    2   3   4
     10     H          1.11980 *      109.03240 *   169.89020 *    3   4   5
     11     H          1.12148 *      108.70283 *   -72.95016 *    3   4   5
     12     H          1.11984 *      109.02670 *   169.88123 *    4   3   2
     13     H          1.12154 *      108.70413 *   -72.94516 *    4   3   2
     14     H          1.10029 *      116.40741 *   102.18978 *    5   4   3
     15     H          1.09946 *      121.39378 *   172.19462 *    6   5   4
     16     H          1.10012 *      123.04727 *   167.03535 *    6   5  15

 (  1,  0)      1.690      2.190    40.715086      0.092    0.10   3
 (  1,  1)      1.690      2.340    35.569657      0.162    0.07
 (  0,  1)      1.540      2.340    29.284640      0.074    0.07
 ( -1,  1)      1.390      2.340    39.472880      0.045    0.05
 ( -1,  0)      1.390      2.190    46.049829      0.113    0.07
 ( -1, -1)      1.390      2.040    55.017053      0.092    0.05
 (  0, -1)      1.540      2.040    43.383839      0.038    0.06
 (  1, -1)      1.690      2.040    47.827082      0.064    0.06
 (  2, -1)      1.840      2.040    56.244542      0.131    0.04
 (  2,  0)      1.840      2.190    50.784948      0.061    0.05
 (  2,  1)      1.840      2.340    46.660623      0.018    0.06
 (  2,  2)      1.840      2.490    44.370286      0.103    0.05
 (  1,  2)      1.690      2.490    32.564778      0.176    0.06
 (  0,  2)      1.540      2.490    25.690258      0.061    0.06
 ( -1,  2)      1.390      2.490    35.349386      0.098    0.05
 ( -2,  2)      1.240      2.490    80.022567      0.144    0.05
 ( -2,  1)      1.240      2.340    84.646131      0.096    0.07
 ( -2,  0)      1.240      2.190    91.814084      0.100    0.05
 ( -2, -1)      1.240      2.040   101.426172      0.013    0.06
 ( -2, -2)      1.240      1.890   111.343234      0.147    0.06
 ( -1, -2)      1.390      1.890    64.209727      0.096    0.05
 (  0, -2)      1.540      1.890    51.534259      0.068    0.06
 (  1, -2)      1.690      1.890    54.293969      0.139    0.05
 (  2, -2)      1.840      1.890    59.588052      0.031    0.07
 (  3, -2)      1.990      1.890    59.881860      0.053    0.06
 (  3, -1)      1.990      2.040    62.164062      0.113    0.06
 (  3,  0)      1.990      2.190    59.758149      0.184    0.06
 (  3,  1)      1.990      2.340    57.177253      0.064    0.06
 (  3,  2)      1.990      2.490    55.843204      0.101    0.06
 (  3,  3)      1.990      2.640    55.294123      0.026    0.08
 (  2,  3)      1.840      2.640    43.183662      0.036    0.06
 (  1,  3)      1.690      2.640    30.884486      0.027    0.07
 (  0,  3)      1.540      2.640    23.594410      0.094    0.05
 ( -1,  3)      1.390      2.640    32.871190      0.031    0.06
 ( -2,  3)      1.240      2.640    77.168857      0.026    0.06
 ( -3,  3)      1.090      2.640   182.248547      0.055    0.06
 ( -3,  2)      1.090      2.490   185.472388      0.029    0.07
 ( -3,  1)      1.090      2.340   190.565660      0.023    0.06
 ( -3,  0)      1.090      2.190   198.268740      0.149    0.06
 ( -3, -1)      1.090      2.040   208.431612      0.054    0.07
 ( -3, -2)      1.090      1.890   218.901830      0.073    0.05
 ( -3, -3)      1.090      1.740   226.074854      0.066    0.08
 ( -2, -3)      1.240      1.740   118.038558      0.016    0.07
 ( -1, -3)      1.390      1.740    70.130497      0.023    0.07
 (  0, -3)      1.540      1.740    56.087957      0.106    0.06
 (  1, -3)      1.690      1.740    56.309770      0.180    0.06
 (  2, -3)      1.840      1.740    57.143421      0.103    0.06
 (  3, -3)      1.990      1.740    52.820088      0.116    0.06
 (  4, -3)      2.140      1.740    46.049474      0.130    0.06
 (  4, -2)      2.140      1.890    55.601675      0.081    0.06
 (  4, -1)      2.140      2.040    62.735013      0.074    0.07
 (  4,  0)      2.140      2.190    65.665444      0.042    0.07
 (  4,  1)      2.140      2.340    65.798849      0.170    0.07
 (  4,  2)      2.140      2.490    65.888239      0.114    0.06
 (  4,  3)      2.140      2.640    66.219725      0.102    0.07
 (  4,  4)      2.140      2.790    66.505399      0.092    0.07
 (  3,  4)      1.990      2.790    55.012373      0.077    0.07
 (  2,  4)      1.840      2.790    42.488310      0.050    0.06
 (  1,  4)      1.690      2.790    29.872941      0.019    0.06
 (  0,  4)      1.540      2.790    22.320478      0.078    0.07
 ( -1,  4)      1.390      2.790    31.354102      0.039    0.07
 ( -2,  4)      1.240      2.790    75.402313      0.153    0.07
 ( -3,  4)      1.090      2.790   180.217575      0.055    0.06
 ( -4,  4)      0.940      2.790   380.493252      0.145    0.06
 ( -4,  3)      0.940      2.640   382.783957      0.159    0.07
 ( -4,  2)      0.940      2.490   386.341387      0.058    0.06
 ( -4,  1)      0.940      2.340   391.835668      0.039    0.06
 ( -4,  0)      0.940      2.190   399.982589      0.041    0.06
 ( -4, -1)      0.940      2.040   410.595964      0.016    0.07
 ( -4, -2)      0.940      1.890   421.512537      0.075    0.07
 ( -4, -3)      0.940      1.740   429.028108      0.141    0.07
 ( -4, -4)      0.940      1.590   431.820637      0.085    0.12
 ( -3, -4)      1.090      1.590   229.264380      0.138    0.09
 ( -2, -4)      1.240      1.590   121.315792      0.081    0.08
 ( -1, -4)      1.390      1.590    73.038927      0.145    0.06
 (  0, -4)      1.540      1.590    57.847484      0.124    0.07
 (  1, -4)      1.690      1.590    55.724895      0.039    0.07
 (  2, -4)      1.840      1.590    53.335676      0.025    0.06
 (  3, -4)      1.990      1.590    46.578555      0.088    0.05
 (  4, -4)      2.140      1.590    38.387051      0.127    0.06
 (  5, -4)      2.290      1.590    31.876958      0.174    0.06
 (  5, -3)      2.290      1.740    40.473556      0.101    0.07
 (  5, -2)      2.290      1.890    51.376128      0.046    0.06
 (  5, -1)      2.290      2.040    60.809929      0.094    0.06
 (  5,  0)      2.290      2.190    67.902045      0.052    0.07
 (  5,  1)      2.290      2.340    72.161152      0.196    0.06
 (  5,  2)      2.290      2.490    74.504151      0.096    0.06
 (  5,  3)      2.290      2.640    76.102134      0.126    0.06
 (  5,  4)      2.290      2.790    77.174452      0.032    0.07
 (  5,  5)      2.290      2.940    77.837854      0.098    0.07
 (  4,  5)      2.140      2.940    66.666752      0.066    0.08
 (  3,  5)      1.990      2.940    54.814965      0.044    0.07
 (  2,  5)      1.840      2.940    42.036659      0.089    0.06
 (  1,  5)      1.690      2.940    29.236332      0.075    0.07
 (  0,  5)      1.540      2.940    21.539412      0.147    0.06
 ( -1,  5)      1.390      2.940    30.443276      0.123    0.08
 ( -2,  5)      1.240      2.940    74.352067      0.124    0.06
 ( -3,  5)      1.090      2.940   179.001726      0.076    0.08
 ( -4,  5)      0.940      2.940   379.090230      0.031    0.08
 ( -5,  5)      0.790      2.940   720.524130      0.115    0.07
 ( -5,  4)      0.790      2.790   722.098432      0.089    0.09
 ( -5,  3)      0.790      2.640   724.584260      0.116    0.05
 ( -5,  2)      0.790      2.490   728.360401      0.069    0.05
 ( -5,  1)      0.790      2.340   734.092028      0.191    0.06
 ( -5,  0)      0.790      2.190   742.488665      0.160    0.08
 ( -5, -1)      0.790      2.040   753.379197      0.099    0.07
 ( -5, -2)      0.790      1.890   764.647140      0.160    0.07
 ( -5, -3)      0.790      1.740   772.282145      0.046    0.13
 ( -5, -4)      0.790      1.590   521.937832      0.094    0.34
 ( -5, -5)      0.790      1.440   523.243635      0.155    0.26
 ( -4, -5)      0.940      1.440   223.363107      0.068    0.07
 ( -3, -5)      1.090      1.440    65.779094      0.132    0.05
 ( -2, -5)      1.240      1.440     1.769214      0.156    0.04
 ( -1, -5)      1.390      1.440    -5.174507      0.152    0.06
 (  0, -5)      1.540      1.440    18.250951      0.190    0.05
 (  1, -5)      1.690      1.440    53.416845      0.085    0.06
 (  2, -5)      1.840      1.440    89.180280      0.149    0.07
 (  3, -5)      1.990      1.440   120.319189      0.123    0.15
 (  4, -5)      2.140      1.440   141.573539      0.087    0.11
 (  5, -5)      2.290      1.440   154.169060      0.099    0.08

 HORIZONTAL STEP=  0.150 ON COORD. 1 OF ATOM   4
   VERTICAL STEP=  0.150 ON COORD. 1 OF ATOM   6
 UPPER LEFT  CORNER COORDINATES:    0.790    2.940 AT (-5, 5)
 LOWER RIGHT CORNER COORDINATES:    2.290    1.440 AT ( 5,-5)
 THE GRID IS RELATIVE TO THE CENTRAL HEAT OF FORMATION=    35.206020

          WHOLE OF (11x11) GRID (Kcal/Mol), SUITABLE FOR PLOTTING:  4

  5 685.32 343.88 143.80  39.15  -4.76 -13.67  -5.97   6.83  19.61  31.46  42.63
  4 686.89 345.29 145.01  40.20  -3.85 -12.89  -5.33   7.28  19.81  31.30  41.97
  3 689.38 347.58 147.04  41.96  -2.33 -11.61  -4.32   7.98  20.09  31.01  40.90
  2 693.15 351.14 150.27  44.82   0.14  -9.52  -2.64   9.16  20.64  30.68  39.30
  1 698.89 356.63 155.36  49.44   4.27  -5.92   0.36  11.45  21.97  30.59  36.96
  0 707.28 364.78 163.06  56.61  10.84   0.00   5.51  15.58  24.55  30.46  32.70
 -1 718.17 375.39 173.23  66.22  19.81   8.18  12.62  21.04  26.96  27.53  25.60
 -2 729.44 386.31 183.70  76.14  29.00  16.33  19.09  24.38  24.68  20.40  16.17
 -3 737.08 393.82 190.87  82.83  34.92  20.88  21.10  21.94  17.61  10.84   5.27
 -4 486.73 396.61 194.06  86.11  37.83  22.64  20.52  18.13  11.37   3.18  -3.33
 -5 488.04 188.16  30.57 -33.44 -40.38 -16.96  18.21  53.97  85.11 106.37 118.96
      -5     -4     -3     -2     -1      0      1      2      3      4      5

 STATIONARY POINT(s) LOCATED ON THE GRID: 5

 SADDLE-POINT ( 1,-4)    HEAT OF FORMATION=     55.7249
    ATOM    CHEMICAL   BOND LENGTH    BOND ANGLE    TWIST ANGLE
   NUMBER   SYMBOL     (ANGSTROMS)     (DEGREES)     (DEGREES)
    (I)                   NA:I          NB:NA:I      NC:NB:NA:I   NA  NB  NC
      1     C
      2     C          1.44147 *                                   1
      3     C          1.41308 *      119.41827 *                  2   1
      4     C          1.69000        104.91213 *   -60.57178 *    3   2   1
      5     C          1.41307 *      104.90891 *    53.77366 *    4   3   2
      6     C          1.59000        107.40400 *    59.88983 *    1   2   3
      7     H          1.11756 *      113.19321 *   179.15800 *    1   2   3
      8     H          1.12009 *      113.36263 *   124.50449 *    1   2   7
      9     H          1.09715 *      119.78746 *  -256.90337 *    2   3   4
     10     H          1.11265 *      106.18778 *   176.59113 *    3   4   5
     11     H          1.11505 *      101.73105 *   -67.71685 *    3   4   5
     12     H          1.11266 *      106.18595 *   176.58984 *    4   3   2
     13     H          1.11503 *      101.72935 *   -67.71724 *    4   3   2
     14     H          1.09714 *      119.78945 *   103.09617 *    5   4   3
     15     H          1.11757 *      113.19069 *   179.16461 *    6   5   4
     16     H          1.12009 *      113.36251 *   124.50029 *    6   5  15


   MAXIMUM    (-5,-3)    HEAT OF FORMATION=    772.2821
    ATOM    CHEMICAL   BOND LENGTH    BOND ANGLE    TWIST ANGLE
   NUMBER   SYMBOL     (ANGSTROMS)     (DEGREES)     (DEGREES)
    (I)                   NA:I          NB:NA:I      NC:NB:NA:I   NA  NB  NC
      1     C
      2     C          1.38691 *                                   1
      3     C          1.62955 *      115.73844 *                  2   1
      4     C          0.79000        124.93643 *   -33.14728 *    3   2   1
      5     C          1.62956 *      124.93584 *    17.63290 *    4   3   2
      6     C          1.74000        102.65410 *    47.98593 *    1   2   3
      7     H          1.11176 *      116.44041 *   166.28250 *    1   2   3
      8     H          1.11378 *      118.91472 *   136.05718 *    1   2   7
      9     H          1.08657 *      116.35559 *  -229.32693 *    2   3   4
     10     H          1.22612 *      124.34594 *   144.66484 *    3   4   5
     11     H          1.23931 *      125.32768 *   -98.38562 *    3   4   5
     12     H          1.22612 *      124.34607 *   144.66985 *    4   3   2
     13     H          1.23932 *      125.32893 *   -98.37928 *    4   3   2
     14     H          1.08656 *      116.35644 *   130.66660 *    5   4   3
     15     H          1.11176 *      116.44161 *   166.28999 *    6   5   4
     16     H          1.11378 *      118.91380 *   136.05657 *    6   5  15


     FULL COMPUTATION TIME :      9.41 SECONDS
 Process Info: 9.3u 0.6s 0:11 89%
1

This is a description of the center point of the grid, including heat of formation and gnorm.

2

This is the geometry of the center point of the grid.

3

This is a list of the AM1 energies of the optimized geometries at each pair of STEP1, STEP2 coordinates. Due to the size of the output, the geometries are not printed. Coordinates and information about each geometry can be found in the visualization (*.vis) file. The AMPAC GUI can be used to examine these geometries.

4

This is an 11×11 matrix of the energies at the points in the order described above. Future versions of the AMPAC GUI will present this data in graphical form.

5

A list of (constrained) critical points and geometries that were found during the grid search. Constrained critical points may not actually correspond to a true critical point. Critical points of interest should be reoptimized with all optimization flags set to on and a FORCE calculation performed on the resulting geometry.