Wall clock time and date at job start Wed Apr 1 2020 06:01:43
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.52998 * 1
3 3 C 1.53003 * 109.47140 * 2 1
4 4 C 1.52995 * 109.47333 * 240.00034 * 2 1 3
5 5 C 1.50703 * 109.47167 * 119.99421 * 2 1 3
6 6 N 1.31873 * 118.94631 * 60.00742 * 5 2 1
7 7 C 1.32151 * 120.96130 * 180.02562 * 6 5 2
8 8 C 1.39945 * 118.92909 * 359.97438 * 7 6 5
9 9 C 1.47829 * 120.94278 * 180.02562 * 8 7 6
10 10 O 1.21553 * 119.99917 * 0.02562 * 9 8 7
11 11 N 1.34777 * 120.00159 * 180.02562 * 9 8 7
12 12 C 1.46498 * 119.99828 * 180.27793 * 11 9 8
13 13 H 1.09004 * 109.47372 * 34.71852 * 12 11 9
14 14 C 1.52997 * 109.47036 * 274.72131 * 12 11 9
15 15 C 1.53003 * 109.47241 * 180.02562 * 14 12 11
16 16 C 1.52999 * 109.47087 * 59.99958 * 15 14 12
17 17 C 1.52997 * 109.47114 * 300.00120 * 16 15 14
18 18 C 1.52999 * 109.46963 * 154.72290 * 12 11 9
19 19 H 1.08995 * 109.47266 * 59.99845 * 18 12 11
20 20 N 1.46503 * 109.47247 * 299.99162 * 18 12 11
21 21 C 1.36937 * 120.00056 * 154.99604 * 20 18 12
22 22 H 1.08005 * 119.99967 * 0.02562 * 21 20 18
23 23 C 1.38806 * 120.00094 * 180.02562 * 21 20 18
24 24 N 1.31618 * 120.00321 * 0.02562 * 23 21 20
25 25 Si 1.86801 * 119.99974 * 180.02562 * 23 21 20
26 26 H 1.48507 * 109.47089 * 359.97438 * 25 23 21
27 27 H 1.48494 * 109.47344 * 120.00318 * 25 23 21
28 28 H 1.48506 * 109.46807 * 240.00352 * 25 23 21
29 29 C 1.39953 * 118.11807 * 0.27864 * 8 7 6
30 30 N 1.31869 * 118.95637 * 240.02417 * 5 2 1
31 31 H 1.09000 * 109.47289 * 300.00202 * 1 2 3
32 32 H 1.08998 * 109.46994 * 60.00193 * 1 2 3
33 33 H 1.09005 * 109.47142 * 180.02562 * 1 2 3
34 34 H 1.08994 * 109.47584 * 179.97438 * 3 2 1
35 35 H 1.09007 * 109.46763 * 300.00355 * 3 2 1
36 36 H 1.09000 * 109.47195 * 59.99740 * 3 2 1
37 37 H 1.09005 * 109.47090 * 299.99579 * 4 2 1
38 38 H 1.08998 * 109.47266 * 59.99939 * 4 2 1
39 39 H 1.09000 * 109.47570 * 180.02562 * 4 2 1
40 40 H 1.07999 * 120.53451 * 180.02562 * 7 6 5
41 41 H 0.97001 * 119.99727 * 359.97438 * 11 9 8
42 42 H 1.09009 * 109.47238 * 300.00198 * 14 12 11
43 43 H 1.08997 * 109.47293 * 60.00168 * 14 12 11
44 44 H 1.09000 * 109.47191 * 180.02562 * 15 14 12
45 45 H 1.08995 * 109.46928 * 299.99922 * 15 14 12
46 46 H 1.08995 * 109.47266 * 60.00171 * 16 15 14
47 47 H 1.09003 * 109.47147 * 180.02562 * 16 15 14
48 48 H 1.08997 * 109.47293 * 180.02562 * 17 16 15
49 49 H 1.09004 * 109.46812 * 300.00187 * 17 16 15
50 50 H 0.96998 * 119.99938 * 335.01407 * 20 18 12
51 51 H 0.96999 * 120.00722 * 179.97438 * 24 23 21
52 52 H 1.07995 * 120.53723 * 179.70435 * 29 8 7
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 6 2.0400 1.4425 0.0000
4 6 2.0400 -0.7212 -1.2492
5 6 2.0323 -0.7103 1.2306
6 7 1.7012 -1.9729 1.4186
7 6 2.1191 -2.6339 2.4839
8 6 2.9247 -1.9728 3.4179
9 6 3.4175 -2.6695 4.6250
10 8 3.1243 -3.8334 4.8171
11 7 4.1918 -2.0143 5.5124
12 6 4.6851 -2.7067 6.7055
13 1 4.8892 -3.7494 6.4622
14 6 3.6274 -2.6381 7.8088
15 6 4.1428 -3.3606 9.0551
16 6 5.4290 -2.6881 9.5392
17 6 6.4867 -2.7566 8.4359
18 6 5.9715 -2.0346 7.1898
19 1 5.7673 -0.9920 7.4331
20 7 6.9844 -2.1004 6.1334
21 6 7.9926 -1.1747 6.0929
22 1 8.0414 -0.3963 6.8401
23 6 8.9520 -1.2367 5.0917
24 7 8.8924 -2.1849 4.1809
25 14 10.3278 0.0257 5.0370
26 1 10.1700 0.9791 6.1647
27 1 11.6377 -0.6643 5.1512
28 1 10.2690 0.7665 3.7513
29 6 3.2462 -0.6316 3.1796
30 7 2.7887 -0.0495 2.0851
31 1 -0.3634 0.5139 0.8900
32 1 -0.3633 0.5138 -0.8900
33 1 -0.3634 -1.0277 0.0005
34 1 3.1299 1.4426 0.0005
35 1 1.6766 1.9563 -0.8900
36 1 1.6766 1.9564 0.8900
37 1 1.6767 -1.7489 -1.2492
38 1 1.6767 -0.2074 -2.1392
39 1 3.1300 -0.7216 -1.2490
40 1 1.8441 -3.6679 2.6308
41 1 4.4260 -1.0856 5.3589
42 1 3.4232 -1.5953 8.0522
43 1 2.7109 -3.1168 7.4639
44 1 3.3890 -3.3122 9.8410
45 1 4.3472 -4.4033 8.8121
46 1 5.2245 -1.6454 9.7822
47 1 5.7964 -3.2031 10.4268
48 1 7.4032 -2.2779 8.7808
49 1 6.6911 -3.7994 8.1929
50 1 6.9407 -2.7995 5.4624
51 1 9.5626 -2.2280 3.4810
52 1 3.8612 -0.0825 3.8771
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
ZINC001085440777.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 06:01:43
Heat of formation + Delta-G solvation = -2.480055 kcal
Electronic energy + Delta-G solvation = -31085.281351 eV
Core-core repulsion = 27084.045690 eV
Total energy + Delta-G solvation = -4001.235661 eV
No. of doubly occupied orbitals = 66
Molecular weight
(most abundant/longest-lived isotopes) = 346.211 amu
Computer time = 26.12 seconds
Orbital eigenvalues (eV)
-41.43366 -39.69838 -37.92616 -36.75633 -34.29625 -34.06444 -33.60827 -31.74534
-31.18664 -29.26516 -27.75711 -27.49929 -27.46793 -25.26143 -24.39306 -23.77027
-22.37117 -21.14722 -20.71048 -19.97399 -19.65226 -17.90207 -17.70718 -17.03331
-16.48213 -15.78619 -15.43172 -15.12570 -14.63682 -14.59303 -14.42331 -14.08944
-14.00126 -13.84541 -13.68271 -13.60142 -13.40866 -13.01615 -12.75415 -12.44889
-12.36849 -12.21530 -12.15440 -11.92562 -11.80285 -11.24993 -11.18883 -11.12668
-11.08523 -10.97839 -10.92493 -10.90282 -10.80505 -10.62778 -10.37524 -10.18138
-9.87123 -9.72308 -9.66356 -9.57407 -9.32806 -8.98378 -8.71471 -7.04966
-5.83572 -3.13861 0.40072 0.98001 2.93646 3.27071 3.37859 3.43649
3.46186 3.60473 3.76336 3.96403 4.41361 4.50617 4.69231 4.73305
4.77407 4.80508 4.90868 5.00880 5.09853 5.21288 5.25587 5.43309
5.46099 5.46178 5.47639 5.62482 5.65382 5.66305 5.70369 5.74781
5.75728 5.79946 5.83362 5.89337 5.92269 5.92980 6.06886 6.09890
6.22655 6.25706 6.26244 6.39760 6.41578 6.43251 6.52757 6.75785
6.81388 6.82634 7.00948 7.62934 7.70458 7.93406 8.32992 8.35530
9.45403 10.00905 10.35956 11.38213
Molecular weight = 346.21amu
Principal moments of inertia in cm(-1)
A = 0.008807 B = 0.003375 C = 0.002622
Principal moments of inertia in units of 10**(-40)*gram-cm**2
A = 3178.650218 B = 8293.804083 C =10676.395291
Net atomic charges, atomic populations, and dipole contributions
calculated with CM2
Atom NO. Type Charge No. of electrons
1 C -0.135 4.135
2 C -0.002 4.002
3 C -0.130 4.130
4 C -0.137 4.137
5 C 0.341 3.659
6 N -0.518 5.518
7 C 0.206 3.794
8 C -0.239 4.239
9 C 0.570 3.430
10 O -0.540 6.540
11 N -0.698 5.698
12 C 0.141 3.859
13 H 0.087 0.913
14 C -0.124 4.124
15 C -0.114 4.114
16 C -0.120 4.120
17 C -0.133 4.133
18 C 0.173 3.827
19 H 0.051 0.949
20 N -0.722 5.722
21 C -0.242 4.242
22 H 0.071 0.929
23 C -0.009 4.009
24 N -0.928 5.928
25 Si 0.904 3.096
26 H -0.280 1.280
27 H -0.292 1.292
28 H -0.292 1.292
29 C 0.189 3.811
30 N -0.513 5.513
31 H 0.053 0.947
32 H 0.066 0.934
33 H 0.067 0.933
34 H 0.064 0.936
35 H 0.058 0.942
36 H 0.062 0.938
37 H 0.067 0.933
38 H 0.066 0.934
39 H 0.056 0.944
40 H 0.180 0.820
41 H 0.402 0.598
42 H 0.062 0.938
43 H 0.063 0.937
44 H 0.058 0.942
45 H 0.062 0.938
46 H 0.058 0.942
47 H 0.055 0.945
48 H 0.066 0.934
49 H 0.057 0.943
50 H 0.385 0.615
51 H 0.255 0.745
52 H 0.174 0.826
For the charges calculated by CM2:
Dipole moment (debyes) X Y Z Total
from point charges -17.014 3.655 -1.950 17.511
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.193 4.193
2 C -0.004 4.004
3 C -0.187 4.187
4 C -0.195 4.195
5 C 0.075 3.925
6 N -0.230 5.230
7 C 0.039 3.961
8 C -0.247 4.247
9 C 0.355 3.645
10 O -0.419 6.419
11 N -0.347 5.347
12 C 0.036 3.964
13 H 0.105 0.895
14 C -0.163 4.163
15 C -0.152 4.152
16 C -0.158 4.158
17 C -0.174 4.174
18 C 0.062 3.938
19 H 0.069 0.931
20 N -0.363 5.363
21 C -0.373 4.373
22 H 0.089 0.911
23 C -0.206 4.206
24 N -0.576 5.576
25 Si 0.735 3.265
26 H -0.205 1.205
27 H -0.219 1.219
28 H -0.219 1.219
29 C 0.022 3.978
30 N -0.226 5.226
31 H 0.073 0.927
32 H 0.085 0.915
33 H 0.086 0.914
34 H 0.083 0.917
35 H 0.077 0.923
36 H 0.081 0.919
37 H 0.086 0.914
38 H 0.085 0.915
39 H 0.075 0.925
40 H 0.197 0.803
41 H 0.237 0.763
42 H 0.080 0.920
43 H 0.081 0.919
44 H 0.077 0.923
45 H 0.081 0.919
46 H 0.076 0.924
47 H 0.074 0.926
48 H 0.085 0.915
49 H 0.076 0.924
50 H 0.219 0.781
51 H 0.065 0.935
52 H 0.191 0.809
Dipole moment (debyes) X Y Z Total
from point charges -17.400 3.455 -2.613 17.931
hybrid contribution 1.542 0.355 0.528 1.668
sum -15.858 3.811 -2.085 16.442
Atomic orbital electron populations
1.21840 0.92434 1.02966 1.02066 1.19445 0.96287 0.95257 0.89399
1.21735 1.02096 0.93158 1.01753 1.21854 1.01318 1.00978 0.95353
1.23333 0.88052 0.87104 0.94057 1.68276 1.24116 1.09337 1.21308
1.23458 0.86856 0.99992 0.85764 1.21847 1.05214 0.97474 1.00157
1.17695 0.78803 0.86103 0.81921 1.90759 1.58551 1.16023 1.76554
1.45808 1.50106 1.15828 1.22975 1.21377 0.93693 0.98035 0.83283
0.89483 1.21875 0.98455 1.00816 0.95156 1.21470 0.97256 0.99851
0.96599 1.21585 0.95538 1.00289 0.98387 1.22210 0.98118 1.01466
0.95595 1.20464 0.87146 0.96878 0.89308 0.93088 1.42434 1.28319
1.35253 1.30307 1.19235 1.02810 1.06090 1.09200 0.91101 1.24733
0.99280 0.98663 0.97885 1.69675 1.36819 1.34241 1.16857 0.86180
0.82362 0.80727 0.77261 1.20536 1.21886 1.21905 1.23562 0.90975
0.90736 0.92493 1.68299 1.11531 1.39282 1.03475 0.92749 0.91506
0.91374 0.91704 0.92328 0.91906 0.91355 0.91495 0.92506 0.80283
0.76287 0.91962 0.91856 0.92319 0.91904 0.92355 0.92627 0.91517
0.92407 0.78149 0.93527 0.80944
Number of geometries 1
Number of calculations of the screened coulomb radii 1
The total number of SCF iterations 14.
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.18 8.37 37.16 0.31 -0.87 16
2 C 0.00 -0.02 0.92 -155.66 -0.14 -0.16 16
3 C -0.13 -1.26 7.66 37.16 0.28 -0.97 16
4 C -0.14 -1.32 8.37 37.15 0.31 -1.01 16
5 C 0.34 4.52 6.25 -157.06 -0.98 3.53 16
6 N -0.52 -7.03 9.77 -10.39 -0.10 -7.13 16
7 C 0.21 2.85 10.88 -16.78 -0.18 2.67 16
8 C -0.24 -3.46 5.91 -104.69 -0.62 -4.08 16
9 C 0.57 9.22 7.75 -12.34 -0.10 9.13 16
10 O -0.54 -10.08 16.10 5.32 0.09 -10.00 16
11 N -0.70 -10.78 4.87 -54.76 -0.27 -11.05 16
12 C 0.14 2.26 1.98 -67.93 -0.13 2.13 16
13 H 0.09 1.61 7.57 -51.93 -0.39 1.21 16
14 C -0.12 -1.66 5.50 -26.73 -0.15 -1.80 16
15 C -0.11 -1.40 5.92 -26.73 -0.16 -1.56 16
16 C -0.12 -1.56 5.92 -26.73 -0.16 -1.72 16
17 C -0.13 -2.20 5.50 -26.76 -0.15 -2.35 16
18 C 0.17 3.09 2.29 -67.95 -0.16 2.94 16
19 H 0.05 0.91 8.08 -51.93 -0.42 0.49 16
20 N -0.72 -16.65 4.95 -53.45 -0.26 -16.91 16
21 C -0.24 -6.26 9.93 -15.06 -0.15 -6.41 16
22 H 0.07 1.78 7.83 -52.48 -0.41 1.37 16
23 C -0.01 -0.25 5.50 -17.43 -0.10 -0.34 16
24 N -0.93 -27.21 13.06 55.49 0.72 -26.48 16
25 Si 0.90 21.66 29.55 -169.99 -5.02 16.64 16
26 H -0.28 -6.48 7.11 56.52 0.40 -6.07 16
27 H -0.29 -7.13 7.11 56.52 0.40 -6.73 16
28 H -0.29 -7.09 7.11 56.52 0.40 -6.69 16
29 C 0.19 2.59 10.81 -16.78 -0.18 2.41 16
30 N -0.51 -7.15 8.98 -10.39 -0.09 -7.24 16
31 H 0.05 0.49 8.14 -51.93 -0.42 0.07 16
32 H 0.07 0.46 8.14 -51.93 -0.42 0.03 16
33 H 0.07 0.63 8.01 -51.93 -0.42 0.21 16
34 H 0.06 0.71 8.04 -51.93 -0.42 0.29 16
35 H 0.06 0.45 8.14 -51.93 -0.42 0.03 16
36 H 0.06 0.63 8.04 -51.93 -0.42 0.21 16
37 H 0.07 0.69 8.01 -51.93 -0.42 0.27 16
38 H 0.07 0.50 8.14 -51.93 -0.42 0.08 16
39 H 0.06 0.60 8.14 -51.93 -0.42 0.18 16
40 H 0.18 2.24 7.74 -52.49 -0.41 1.83 16
41 H 0.40 5.78 6.74 -40.82 -0.28 5.50 16
42 H 0.06 0.77 8.14 -51.92 -0.42 0.35 16
43 H 0.06 0.88 8.14 -51.93 -0.42 0.46 16
44 H 0.06 0.62 8.14 -51.93 -0.42 0.20 16
45 H 0.06 0.80 8.14 -51.93 -0.42 0.38 16
46 H 0.06 0.72 8.14 -51.93 -0.42 0.30 16
47 H 0.05 0.66 8.14 -51.93 -0.42 0.24 16
48 H 0.07 1.20 8.14 -51.93 -0.42 0.78 16
49 H 0.06 1.01 8.14 -51.93 -0.42 0.59 16
50 H 0.38 9.87 7.80 -40.82 -0.32 9.55 16
51 H 0.25 7.31 8.31 -40.82 -0.34 6.97 16
52 H 0.17 2.04 6.54 -52.49 -0.34 1.70 16
LS Contribution 416.60 15.07 6.28 6.28
Total: -1.00 -30.60 416.60 -9.97 -40.56
By element:
Atomic # 1 Polarization: 22.67 SS G_CDS: -8.86 Total: 13.81 kcal
Atomic # 6 Polarization: 3.97 SS G_CDS: -2.44 Total: 1.53 kcal
Atomic # 7 Polarization: -68.82 SS G_CDS: 0.00 Total: -68.82 kcal
Atomic # 8 Polarization: -10.08 SS G_CDS: 0.09 Total: -10.00 kcal
Atomic # 14 Polarization: 21.66 SS G_CDS: -5.02 Total: 16.64 kcal
Total LS contribution 6.28 Total: 6.28 kcal
Total: -30.60 -9.97 -40.56 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.
ZINC001085440777.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 38.083 kcal
(2) G-P(sol) polarization free energy of solvation -30.595 kcal
(3) G-ENP(sol) elect.-nuc.-pol. free energy of system 7.488 kcal
(4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.968 kcal
(5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -40.563 kcal
(6) G-S(sol) free energy of system = (1) + (5) -2.480 kcal
1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED
Total computer time = 26.12 seconds