Wall clock time and date at job start Wed Apr 1 2020 04:05:48
AMSOL-version 7.1
by G. D. Hawkins, D. J. Giesen, G. C. Lynch, C. C. Chambers,
I. Rossi, J. W. Storer, J. Li, J. D. Thompson, P. Winget,
B. J. Lynch, D. Rinaldi, D. A. Liotard, C. J. Cramer, and D. G. Truhlar
Copyright 2004 by Regents of the University of Minnesota.
All rights reserved.
Notice: recipients of this code are asked to comply
with the user agreement in Section 1 of the documentation file.
*******************************************************************************
* 1SCF - SCF CALCULATION WITHOUT GEOMETRY OPTIMIZATION
* - USE EF ROUTINE FOR MINIMUM SEARCH (DEFAULT)
* GEO-OK - OVERRIDE INTERATOMIC DISTANCE CHECK
* TLIMIT= - A TIME OF 15. SECONDS REQUESTED
* CHARGE - CHARGE ON SYSTEM= -1
* AM1 - THE AM1 HAMILTONIAN TO BE USED
* DEV - DEVELOPER OPTIONS ARE ALLOWED
* SM5.42R - SM5.42R CALCULATIONS WILL BE PERFORMED
* - CHARGE MODEL 2 WILL BE USED (CM2A)
* - THE SOLVENT IS USER-SPECIFIED
* DIELEC - THE SOLVENT DIELECTRIC CONSTANT IS 2.06
* IOFR - THE SOLVENT INDEX OF REFRACTION IS 1.4345
* ALPHA - THE SOLVENT ALPHA IS 0.00
* BETA - THE SOLVENT BETA IS 0.00
* GAMMA - THE SOLVENT MACROSCOPIC SURFACE TENSION IS
* 38.93 CAL MOL^-1 ANGSTROM^-2
* FACARB - THE FRACTION OF AROMATIC CARBONS IS 0.00
* FEHALO - THE FRACTION OF ELECTRONEGATIVE HALOGENS IS 0.00
*******************************************************************************
Atom NO. Chemical Bond length Bond angle Dihedral angle
number (k) symbol (angstroms) (degrees) (degrees)
(I) NA:I NB:NA:I NC:NB:NA:I NA NB NC
1 1 C
2 2 C 1.53000 * 1
3 3 H 1.09013 * 109.46847 * 2 1
4 4 C 1.50696 * 109.46954 * 240.00296 * 2 1 3
5 5 O 1.21287 * 119.99669 * 239.99240 * 4 2 1
6 6 N 1.34777 * 120.00417 * 60.00337 * 4 2 1
7 7 C 1.47199 * 126.70005 * 359.67511 * 6 4 2
8 8 C 1.54034 * 106.26369 * 148.64194 * 7 6 4
9 9 C 1.53001 * 112.77318 * 141.67694 * 8 7 6
10 10 C 1.56817 * 109.45009 * 268.28686 * 8 7 6
11 11 N 1.48111 * 105.10238 * 88.54223 * 10 8 7
12 12 C 1.36938 * 127.84525 * 219.46101 * 11 10 8
13 13 H 1.08000 * 119.99711 * 15.65982 * 12 11 10
14 14 C 1.38809 * 120.00225 * 195.65939 * 12 11 10
15 15 N 1.31619 * 120.00179 * 9.11562 * 14 12 11
16 16 Si 1.86799 * 120.00103 * 189.11791 * 14 12 11
17 17 H 1.48500 * 109.47018 * 0.02562 * 16 14 12
18 18 H 1.48503 * 109.46853 * 119.99941 * 16 14 12
19 19 H 1.48501 * 109.47116 * 239.99790 * 16 14 12
20 20 C 1.48116 * 104.31203 * 39.46244 * 11 10 8
21 21 C 1.54316 * 104.71058 * 320.22419 * 20 11 10
22 22 C 1.53000 * 112.07168 * 146.94444 * 21 20 11
23 23 C 1.47194 * 126.70470 * 180.02562 * 6 4 2
24 24 C 1.50696 * 109.47098 * 119.99870 * 2 1 3
25 25 C 1.35316 * 126.12997 * 85.02798 * 24 2 1
26 26 C 1.39860 * 107.92434 * 179.75840 * 25 24 2
27 27 N 1.30840 * 108.20313 * 0.24295 * 26 25 24
28 28 N 1.34955 * 126.12819 * 265.20027 * 24 2 1
29 29 C 1.46498 * 126.03479 * 0.02562 * 28 24 2
30 30 H 1.08999 * 109.47059 * 180.02562 * 1 2 3
31 31 H 1.08999 * 109.47203 * 300.00111 * 1 2 3
32 32 H 1.09001 * 109.47490 * 60.00129 * 1 2 3
33 33 H 1.09000 * 110.07280 * 267.95443 * 7 6 4
34 34 H 1.09004 * 110.06807 * 29.48705 * 7 6 4
35 35 H 1.08998 * 109.46931 * 301.38018 * 9 8 7
36 36 H 1.09003 * 109.46949 * 61.38255 * 9 8 7
37 37 H 1.09003 * 109.47082 * 181.38065 * 9 8 7
38 38 H 1.08996 * 110.47408 * 329.46723 * 10 8 7
39 39 H 1.09004 * 110.02653 * 207.38647 * 10 8 7
40 40 H 0.96997 * 120.00729 * 179.97438 * 15 14 12
41 41 H 1.08995 * 110.38825 * 201.41880 * 20 11 10
42 42 H 1.08995 * 110.38791 * 79.03024 * 20 11 10
43 43 H 1.09005 * 109.47251 * 182.00953 * 22 21 20
44 44 H 1.09001 * 109.47563 * 302.00833 * 22 21 20
45 45 H 1.09000 * 109.47234 * 62.01129 * 22 21 20
46 46 H 1.09005 * 110.07063 * 330.52259 * 23 6 4
47 47 H 1.09000 * 110.07504 * 91.99541 * 23 6 4
48 48 H 1.08002 * 126.03223 * 0.02562 * 25 24 2
49 49 H 1.08002 * 125.90433 * 180.21931 * 26 25 24
50 50 H 1.09005 * 109.46658 * 270.25042 * 29 28 24
51 51 H 1.08991 * 109.47849 * 30.25455 * 29 28 24
52 52 H 1.09002 * 109.47114 * 150.25771 * 29 28 24
Note: An asterisk (*) indicates this parameter was or will be optimized.
Cartesian coordinates (angstroms)
NO. Atom X Y Z
1 6 0.0000 0.0000 0.0000
2 6 1.5300 0.0000 0.0000
3 1 1.8933 1.0278 0.0000
4 6 2.0323 -0.7103 -1.2305
5 8 2.7297 -1.6964 -1.1192
6 7 1.7068 -0.2498 -2.4546
7 6 0.8700 0.9221 -2.7598
8 6 1.4055 1.5018 -4.0826
9 6 0.2918 2.0105 -5.0001
10 6 2.4889 2.5927 -3.7741
11 7 3.7690 1.8523 -3.6900
12 6 4.7997 2.0534 -2.8112
13 1 4.6419 2.6418 -1.9194
14 6 6.0473 1.5010 -3.0665
15 7 6.2923 0.9387 -4.2310
16 14 7.3785 1.5554 -1.7572
17 1 6.8541 2.2362 -0.5461
18 1 7.7866 0.1708 -1.4084
19 1 8.5528 2.3018 -2.2760
20 6 3.6811 0.8451 -4.7724
21 6 2.2193 0.3523 -4.7314
22 6 1.6963 -0.0125 -6.1222
23 6 2.1126 -0.8247 -3.7474
24 6 2.0323 -0.7104 1.2304
25 6 2.2089 -2.0452 1.3649
26 6 2.6928 -2.2876 2.6545
27 7 2.7978 -1.1401 3.2741
28 7 2.3835 -0.1295 2.3969
29 6 2.3435 1.3049 2.6918
30 1 -0.3633 -1.0277 0.0005
31 1 -0.3633 0.5138 0.8900
32 1 -0.3634 0.5138 -0.8900
33 1 -0.1706 0.6189 -2.8756
34 1 0.9582 1.6621 -1.9642
35 1 -0.3949 1.1944 -5.2249
36 1 -0.2501 2.8147 -4.5023
37 1 0.7270 2.3848 -5.9267
38 1 2.2733 3.0870 -2.8269
39 1 2.5247 3.3225 -4.5829
40 1 7.1641 0.5530 -4.4096
41 1 4.3658 0.0193 -4.5792
42 1 3.9003 1.3037 -5.7366
43 1 0.6531 -0.3199 -6.0490
44 1 1.7744 0.8542 -6.7787
45 1 2.2888 -0.8316 -6.5299
46 1 3.0790 -1.3195 -3.6499
47 1 1.3620 -1.5356 -4.0929
48 1 2.0150 -2.7911 0.6082
49 1 2.9372 -3.2551 3.0678
50 1 3.2974 1.7580 2.4220
51 1 1.5439 1.7719 2.1169
52 1 2.1595 1.4520 3.7560
RHF calculation, no. of doubly occupied orbitals= 66
REFERENCES FOR PARAMETERS IN GAS-PHASE HAMILTONIAN:
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)
N: (AM1): M.J.S. DEWAR ET AL, J. AM. CHEM. SOC. 107 3902-3909 (1985)
O: (AM1): M.J.S. DEWAR ET AL, J. AM. CHEM. SOC. 107 3902-3909 (1985)
Si: (AM1): M.J.S.DEWAR, C. JIE, ORGANOMETALLICS, 6, 1486-1490 (1987).
REFERENCE FOR CHARGE MODEL 2:
J. Li, J. Xing, C. J. Cramer, and D. G. Truhlar, J. Chem. Phys. 111 (1999) 885.
REFERENCE FOR THE PARAMETERIZATION OF THE SM5.42R SOLVATION MODEL:
J. Li, T. Zhu, G. D. Hawkins, P. Winget, D. A. Liotard, C. J. Cramer, and
D. G. Truhlar, Theor. Chem. Acc. 103 (1999) 9-63
P. Winget, J. D. Thompson C. J. Cramer,and D. G. Truhlar, J. Phys. Chem. B. submitted.
CHARGE=-1 AM1 1SCF TLIMIT=15 GEO-OK SM5.42R
& SOLVNT=GENORG IOFR=1.4345 ALPHA=0.00 BETA=0.00 GAMMA=38.93
& DIELEC=2.06 FACARB=0.00 FEHALO=0.00 DEV
ZINC000982282090.mol2 52
1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED
SCF FIELD WAS ACHIEVED
Wall clock time and date at job start Wed Apr 1 2020 04:05:48
Heat of formation + Delta-G solvation = 82.934624 kcal
Electronic energy + Delta-G solvation = -33360.418902 eV
Core-core repulsion = 29362.887099 eV
Total energy + Delta-G solvation = -3997.531803 eV
No. of doubly occupied orbitals = 66
Molecular weight
(most abundant/longest-lived isotopes) = 346.211 amu
Computer time = 1.28 seconds
Orbital eigenvalues (eV)
-41.24720 -40.29608 -38.32690 -34.75472 -34.15931 -33.00542 -32.77225 -31.59561
-30.74710 -28.96454 -26.85152 -26.59299 -26.02495 -25.80858 -24.86481 -24.20147
-22.15250 -21.07266 -20.21261 -20.01511 -18.74770 -17.81850 -16.98632 -16.77603
-15.96984 -15.61722 -15.35531 -15.12771 -14.65067 -14.46809 -14.06557 -13.89043
-13.61027 -13.47981 -13.39046 -13.27465 -13.11842 -12.88462 -12.69448 -12.44046
-12.33288 -12.29109 -11.97842 -11.78781 -11.68952 -11.55896 -11.43941 -11.39013
-11.23255 -10.85389 -10.82074 -10.61023 -10.44250 -10.12786 -10.04413 -9.79122
-9.69824 -9.55318 -9.34747 -9.07454 -8.69873 -8.21869 -8.16692 -6.83372
-5.79585 -3.13037 2.01403 2.61351 2.99967 3.12592 3.38392 3.42698
3.44601 3.47702 4.42849 4.46118 4.59681 4.69344 4.77656 4.97712
5.03974 5.06812 5.17473 5.22272 5.30556 5.39134 5.43038 5.44698
5.51344 5.52220 5.55722 5.64968 5.66036 5.69746 5.86045 5.88008
5.94689 6.01258 6.03289 6.07202 6.11656 6.15623 6.25324 6.30119
6.44114 6.48239 6.49333 6.60146 6.65137 6.79273 6.90562 6.92607
7.16792 7.28943 7.52943 7.71183 7.89455 7.98289 8.21207 8.58459
9.10522 9.74687 10.50481 11.34317
Molecular weight = 346.21amu
Principal moments of inertia in cm(-1)
A = 0.009102 B = 0.005272 C = 0.003959
Principal moments of inertia in units of 10**(-40)*gram-cm**2
A = 3075.394283 B = 5309.914375 C = 7069.966264
Net atomic charges, atomic populations, and dipole contributions
calculated with CM2
Atom NO. Type Charge No. of electrons
1 C -0.142 4.142
2 C -0.042 4.042
3 H 0.105 0.895
4 C 0.547 3.453
5 O -0.526 6.526
6 N -0.595 5.595
7 C 0.121 3.879
8 C -0.092 4.092
9 C -0.122 4.122
10 C 0.142 3.858
11 N -0.588 5.588
12 C -0.222 4.222
13 H 0.070 0.930
14 C 0.001 3.999
15 N -0.912 5.912
16 Si 0.904 3.096
17 H -0.281 1.281
18 H -0.289 1.289
19 H -0.290 1.290
20 C 0.172 3.828
21 C -0.097 4.097
22 C -0.123 4.123
23 C 0.130 3.870
24 C 0.094 3.906
25 C -0.230 4.230
26 C 0.038 3.962
27 N -0.302 5.302
28 N -0.373 5.373
29 C 0.099 3.901
30 H 0.066 0.934
31 H 0.063 0.937
32 H 0.075 0.925
33 H 0.060 0.940
34 H 0.070 0.930
35 H 0.055 0.945
36 H 0.051 0.949
37 H 0.061 0.939
38 H 0.042 0.958
39 H 0.023 0.977
40 H 0.253 0.747
41 H 0.087 0.913
42 H 0.021 0.979
43 H 0.050 0.950
44 H 0.059 0.941
45 H 0.054 0.946
46 H 0.080 0.920
47 H 0.052 0.948
48 H 0.153 0.847
49 H 0.180 0.820
50 H 0.075 0.925
51 H 0.084 0.916
52 H 0.087 0.913
For the charges calculated by CM2:
Dipole moment (debyes) X Y Z Total
from point charges -16.619 0.093 3.782 17.044
Note: The Mulliken population analysis charges presented below were not
used in the solvation calculation but are provided for completeness.
The chosen solvation model uses CM2 partial charges for calculating
solvation energies.
Net atomic charges, atomic populations, and dipole contributions
using Mulliken population analysis
Atom NO. Type Charge No. of electrons
1 C -0.200 4.200
2 C -0.066 4.066
3 H 0.123 0.877
4 C 0.334 3.666
5 O -0.402 6.402
6 N -0.327 5.327
7 C -0.003 4.003
8 C -0.095 4.095
9 C -0.180 4.180
10 C 0.015 3.985
11 N -0.310 5.310
12 C -0.350 4.350
13 H 0.088 0.912
14 C -0.197 4.197
15 N -0.558 5.558
16 Si 0.735 3.265
17 H -0.206 1.206
18 H -0.215 1.215
19 H -0.217 1.217
20 C 0.046 3.954
21 C -0.100 4.100
22 C -0.181 4.181
23 C 0.007 3.993
24 C -0.033 4.033
25 C -0.258 4.258
26 C -0.147 4.147
27 N -0.124 5.124
28 N -0.157 5.157
29 C -0.043 4.043
30 H 0.085 0.915
31 H 0.082 0.918
32 H 0.094 0.906
33 H 0.079 0.921
34 H 0.089 0.911
35 H 0.074 0.926
36 H 0.070 0.930
37 H 0.080 0.920
38 H 0.061 0.939
39 H 0.041 0.959
40 H 0.063 0.937
41 H 0.105 0.895
42 H 0.039 0.961
43 H 0.069 0.931
44 H 0.078 0.922
45 H 0.073 0.927
46 H 0.099 0.901
47 H 0.071 0.929
48 H 0.171 0.829
49 H 0.197 0.803
50 H 0.093 0.907
51 H 0.102 0.898
52 H 0.106 0.894
Dipole moment (debyes) X Y Z Total
from point charges -15.593 0.214 5.148 16.422
hybrid contribution 0.761 -0.184 -0.134 0.794
sum -14.832 0.031 5.014 15.656
Atomic orbital electron populations
1.21887 0.91610 1.03483 1.03030 1.19678 0.97436 1.00096 0.89351
0.87717 1.19379 0.80056 0.84449 0.82710 1.90715 1.37368 1.26684
1.85463 1.48012 1.50526 1.29200 1.05003 1.22135 0.93149 0.87703
0.97272 1.22878 0.95632 0.97378 0.93610 1.21594 0.97021 1.00281
0.99109 1.22094 0.86681 0.91966 0.97770 1.44720 1.14707 1.34860
1.36718 1.18837 0.87388 1.29777 0.98990 0.91175 1.24792 0.96893
1.01422 0.96574 1.69871 1.22142 1.43037 1.20732 0.86177 0.81972
0.77517 0.80854 1.20625 1.21529 1.21687 1.21613 0.89272 0.93881
0.90615 1.22567 1.00057 0.93658 0.93686 1.21633 0.99682 1.01341
0.95453 1.21977 1.00470 0.95123 0.81700 1.20272 1.01526 0.93526
0.87939 1.21375 1.11514 0.95036 0.97923 1.23812 0.99947 0.96899
0.94062 1.77276 1.22525 0.95818 1.16754 1.48353 1.51281 1.06722
1.09342 1.21827 1.03364 0.75741 1.03354 0.91485 0.91752 0.90634
0.92148 0.91129 0.92636 0.92968 0.91958 0.93948 0.95910 0.93662
0.89469 0.96142 0.93131 0.92224 0.92699 0.90108 0.92943 0.82947
0.80274 0.90654 0.89751 0.89397
Number of geometries 1
Number of calculations of the screened coulomb radii 1
The total number of SCF iterations 16.
Note: The number of geometries may not
correspond to the number of cycles due to rejected geometry changes.
In the following table subtotal= G_P + SS G_CDS.
Atom Chem. CM2 G_P Area Sigma k SS G_CDS Subtotal M
number symbol chg. (kcal) (Ang**2) cal/(Ang**2) (kcal) (kcal) value
1 C -0.14 -1.50 8.94 37.16 0.33 -1.17 16
2 C -0.04 -0.62 2.06 -92.93 -0.19 -0.81 16
3 H 0.10 1.48 7.23 -51.92 -0.38 1.11 16
4 C 0.55 9.96 7.17 -10.99 -0.08 9.88 16
5 O -0.53 -11.56 13.97 5.55 0.08 -11.48 16
6 N -0.60 -10.01 3.31 -171.45 -0.57 -10.58 16
7 C 0.12 1.63 4.31 -3.12 -0.01 1.62 16
8 C -0.09 -1.36 1.15 -150.81 -0.17 -1.53 16
9 C -0.12 -1.37 7.65 37.16 0.28 -1.08 16
10 C 0.14 2.66 6.37 -1.09 -0.01 2.65 16
11 N -0.59 -14.03 3.38 -167.41 -0.57 -14.60 16
12 C -0.22 -6.02 10.35 -15.06 -0.16 -6.18 16
13 H 0.07 1.87 7.85 -52.49 -0.41 1.46 16
14 C 0.00 0.02 5.40 -17.43 -0.09 -0.08 16
15 N -0.91 -26.73 10.87 55.49 0.60 -26.12 16
16 Si 0.90 22.91 29.56 -169.99 -5.03 17.89 16
17 H -0.28 -6.93 7.11 56.52 0.40 -6.53 16
18 H -0.29 -7.57 7.11 56.52 0.40 -7.17 16
19 H -0.29 -7.28 7.11 56.52 0.40 -6.88 16
20 C 0.17 3.89 4.91 -2.46 -0.01 3.87 16
21 C -0.10 -1.63 1.06 -152.32 -0.16 -1.79 16
22 C -0.12 -1.67 7.50 37.16 0.28 -1.39 16
23 C 0.13 2.28 6.08 -3.26 -0.02 2.26 16
24 C 0.09 1.53 5.73 -82.26 -0.47 1.06 16
25 C -0.23 -3.85 9.14 -40.02 -0.37 -4.21 16
26 C 0.04 0.59 12.02 -17.37 -0.21 0.38 16
27 N -0.30 -5.05 12.58 30.62 0.39 -4.67 16
28 N -0.37 -5.80 4.01 -63.89 -0.26 -6.06 16
29 C 0.10 1.11 11.06 59.85 0.66 1.77 16
30 H 0.07 0.73 8.14 -51.93 -0.42 0.31 16
31 H 0.06 0.57 8.14 -51.93 -0.42 0.15 16
32 H 0.07 0.70 5.15 -51.93 -0.27 0.43 16
33 H 0.06 0.66 7.33 -51.93 -0.38 0.28 16
34 H 0.07 0.91 6.22 -51.93 -0.32 0.58 16
35 H 0.05 0.55 7.09 -51.93 -0.37 0.18 16
36 H 0.05 0.51 8.14 -51.93 -0.42 0.08 16
37 H 0.06 0.72 7.17 -51.93 -0.37 0.35 16
38 H 0.04 0.76 7.53 -51.93 -0.39 0.37 16
39 H 0.02 0.44 8.00 -51.93 -0.42 0.02 16
40 H 0.25 7.27 8.31 -40.82 -0.34 6.93 16
41 H 0.09 2.31 5.79 -51.93 -0.30 2.01 16
42 H 0.02 0.49 7.97 -51.93 -0.41 0.08 16
43 H 0.05 0.57 7.10 -51.93 -0.37 0.20 16
44 H 0.06 0.79 7.16 -51.93 -0.37 0.41 16
45 H 0.05 0.74 8.14 -51.93 -0.42 0.32 16
46 H 0.08 1.72 7.14 -51.93 -0.37 1.35 16
47 H 0.05 0.80 7.89 -51.93 -0.41 0.39 16
48 H 0.15 2.61 6.80 -52.49 -0.36 2.25 16
49 H 0.18 2.14 8.06 -52.49 -0.42 1.72 16
50 H 0.07 0.87 8.14 -51.93 -0.42 0.45 16
51 H 0.08 0.73 7.56 -51.93 -0.39 0.33 16
52 H 0.09 0.83 8.14 -51.93 -0.42 0.40 16
LS Contribution 396.08 15.07 5.97 5.97
Total: -1.00 -34.66 396.08 -8.16 -42.82
By element:
Atomic # 1 Polarization: 9.97 SS G_CDS: -8.38 Total: 1.59 kcal
Atomic # 6 Polarization: 5.64 SS G_CDS: -0.40 Total: 5.25 kcal
Atomic # 7 Polarization: -61.63 SS G_CDS: -0.40 Total: -62.03 kcal
Atomic # 8 Polarization: -11.56 SS G_CDS: 0.08 Total: -11.48 kcal
Atomic # 14 Polarization: 22.91 SS G_CDS: -5.03 Total: 17.89 kcal
Total LS contribution 5.97 Total: 5.97 kcal
Total: -34.66 -8.16 -42.82 kcal
The number of atoms in the molecule is 52
The average number of expansion shells was 16.00
The maximum number of expansion shells was 16
The minimum number of expansion shells was 16
**** NOTA BENE ****
This is the net solvation energy for this exact molecular structure
(nuclear and electronic)! The standard-state solvation energy should be
obtained as the difference between the heat of formation plus delta-G solvation
for the relaxed solvated system and that for the relaxed gas-phase system.
ZINC000982282090.mol2 52
1SCF run
This is a breakdown of the solvation energy calculated
without geometric relaxation in solution:
(1) E-EN(sol) electronic-nuclear energy of solute 125.758 kcal
(2) G-P(sol) polarization free energy of solvation -34.663 kcal
(3) G-ENP(sol) elect.-nuc.-pol. free energy of system 91.095 kcal
(4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.161 kcal
(5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -42.824 kcal
(6) G-S(sol) free energy of system = (1) + (5) 82.935 kcal
1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED
Total computer time = 1.28 seconds