Wall clock time and date at job start Tue Mar 31 2020 21:38:53 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.52993 * 1 3 3 H 1.08997 * 109.55477 * 2 1 4 4 C 1.53296 * 109.51193 * 239.87513 * 2 1 3 5 5 N 1.47159 * 108.40379 * 187.70036 * 4 2 1 6 6 C 1.34783 * 120.98626 * 230.61232 * 5 4 2 7 7 O 1.21276 * 119.99574 * 179.97438 * 6 5 4 8 8 C 1.50690 * 120.00164 * 359.97438 * 6 5 4 9 9 H 1.09008 * 115.55423 * 325.58655 * 8 6 5 10 10 C 1.53001 * 117.49711 * 111.25329 * 8 6 5 11 11 C 1.53002 * 117.50215 * 179.91337 * 8 6 5 12 12 C 1.53005 * 117.49378 * 0.02562 * 11 8 6 13 13 C 1.52999 * 117.50298 * 145.02013 * 11 8 6 14 14 C 1.47224 * 117.97278 * 50.58611 * 5 4 2 15 15 C 1.53258 * 108.43057 * 309.53966 * 14 5 4 16 16 C 1.53004 * 109.36248 * 172.08309 * 15 14 5 17 17 C 1.53184 * 109.34109 * 181.17051 * 16 15 14 18 18 N 1.46923 * 108.77234 * 305.41379 * 17 16 15 19 19 C 1.36934 * 120.64213 * 233.69277 * 18 17 16 20 20 H 1.08003 * 119.99962 * 27.67107 * 19 18 17 21 21 C 1.38812 * 120.00102 * 207.67174 * 19 18 17 22 22 N 1.31614 * 120.00506 * 17.36646 * 21 19 18 23 23 Si 1.86805 * 119.99922 * 197.36772 * 21 19 18 24 24 H 1.48507 * 109.47007 * 89.99838 * 23 21 19 25 25 H 1.48501 * 109.46937 * 209.99813 * 23 21 19 26 26 H 1.48495 * 109.47310 * 329.99984 * 23 21 19 27 27 C 1.46933 * 118.71924 * 53.67213 * 18 17 16 28 28 C 1.52979 * 109.40034 * 292.17295 * 15 14 5 29 29 O 1.42891 * 109.32243 * 52.15338 * 15 14 5 30 30 H 1.09004 * 109.47390 * 180.12431 * 1 2 3 31 31 H 1.09004 * 109.47602 * 300.12105 * 1 2 3 32 32 H 1.08998 * 109.47435 * 60.12403 * 1 2 3 33 33 H 1.09004 * 109.68801 * 307.44022 * 4 2 1 34 34 H 1.09006 * 109.68563 * 67.96214 * 4 2 1 35 35 H 1.09004 * 117.49994 * 214.97742 * 10 8 6 36 36 H 1.08997 * 117.50372 * 359.97438 * 10 8 6 37 37 H 1.08999 * 109.47504 * 59.99546 * 12 11 8 38 38 H 1.08997 * 109.47187 * 179.97438 * 12 11 8 39 39 H 1.09003 * 109.46795 * 299.99602 * 12 11 8 40 40 H 1.09001 * 109.46634 * 60.00243 * 13 11 8 41 41 H 1.08993 * 109.47234 * 180.02562 * 13 11 8 42 42 H 1.08998 * 109.47108 * 300.00887 * 13 11 8 43 43 H 1.09003 * 109.66418 * 189.80458 * 14 5 4 44 44 H 1.09004 * 109.70384 * 69.29898 * 14 5 4 45 45 H 1.08996 * 109.49044 * 301.13594 * 16 15 14 46 46 H 1.09007 * 109.49009 * 61.21147 * 16 15 14 47 47 H 1.09006 * 109.58825 * 185.62504 * 17 16 15 48 48 H 1.09007 * 109.58516 * 65.20040 * 17 16 15 49 49 H 0.96997 * 120.00937 * 179.97438 * 22 21 19 50 50 H 1.08997 * 109.58230 * 186.55251 * 27 18 17 51 51 H 1.08998 * 109.58268 * 66.13247 * 27 18 17 52 52 H 1.09010 * 109.49167 * 298.82704 * 28 15 14 53 53 H 1.08996 * 109.49348 * 58.90356 * 28 15 14 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.5299 0.0000 0.0000 3 1 1.8947 1.0271 0.0000 4 6 2.0419 -0.7252 -1.2498 5 7 3.5014 -0.8748 -1.1349 6 6 4.3185 -0.5006 -2.1394 7 8 5.5188 -0.6316 -2.0253 8 6 3.7407 0.0804 -3.4040 9 1 2.7779 -0.3277 -3.7119 10 6 4.0026 1.5614 -3.6849 11 6 4.7204 0.4876 -4.5064 12 6 6.2119 0.2731 -4.2405 13 6 4.2737 0.2913 -5.9565 14 6 4.0412 -1.4450 0.1105 15 6 3.4305 -0.6877 1.2947 16 6 3.8205 -1.3861 2.5990 17 6 3.1858 -0.6458 3.7803 18 7 3.5729 0.7700 3.7142 19 6 4.1561 1.3851 4.7896 20 1 4.7337 0.8062 5.4951 21 6 4.0049 2.7527 4.9732 22 7 3.0493 3.4000 4.3407 23 14 5.1580 3.6746 6.1179 24 1 6.3341 4.1586 5.3512 25 1 4.4437 4.8329 6.7123 26 1 5.6128 2.7651 7.2000 27 6 3.3227 1.5131 2.4715 28 6 3.9555 0.7492 1.3041 29 8 2.0075 -0.6763 1.1655 30 1 -0.3634 -1.0277 0.0022 31 1 -0.3634 0.5157 0.8889 32 1 -0.3634 0.5119 -0.8911 33 1 1.5765 -1.7083 -1.3204 34 1 1.8016 -0.1406 -2.1379 35 1 3.2123 2.1276 -4.1780 36 1 4.5924 2.1156 -2.9548 37 1 6.4009 -0.7857 -4.0638 38 1 6.7865 0.6051 -5.1052 39 1 6.5107 0.8469 -3.3632 40 1 3.3728 0.8762 -6.1414 41 1 5.0659 0.6213 -6.6285 42 1 4.0644 -0.7638 -6.1325 43 1 5.1256 -1.3351 0.1244 44 1 3.7786 -2.5009 0.1754 45 1 3.4637 -2.4158 2.5809 46 1 4.9053 -1.3785 2.7054 47 1 3.5406 -1.0783 4.7159 48 1 2.1005 -0.7309 3.7244 49 1 2.9433 4.3556 4.4692 50 1 3.7681 2.5053 2.5432 51 1 2.2486 1.6032 2.3098 52 1 5.0393 0.7397 1.4206 53 1 3.6945 1.2380 0.3655 RHF calculation, no. of doubly occupied orbitals= 65 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 ZINC001098586715.mol2 53 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Tue Mar 31 2020 21:38:53 Heat of formation + Delta-G solvation = -28.439664 kcal Electronic energy + Delta-G solvation = -31952.147067 eV Core-core repulsion = 28018.407624 eV Total energy + Delta-G solvation = -3933.739443 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 1.27 seconds Orbital eigenvalues (eV) -40.83809 -39.45113 -37.19276 -36.34476 -34.53749 -33.53907 -31.47909 -30.24193 -29.91132 -29.45672 -27.59358 -27.34541 -26.27279 -24.64304 -23.20848 -22.54903 -21.65952 -20.75981 -19.68648 -18.97365 -18.52611 -17.39678 -16.93597 -15.99314 -15.56799 -15.46568 -14.97720 -14.62544 -14.60055 -14.26844 -14.11245 -13.85543 -13.55350 -13.46222 -13.22190 -13.18915 -12.96544 -12.73143 -12.54040 -12.37561 -12.12738 -11.99449 -11.94353 -11.68597 -11.48530 -11.42479 -11.40444 -11.35699 -11.00118 -10.73734 -10.53371 -10.38076 -10.32080 -10.23369 -10.05345 -9.94883 -9.79817 -9.58374 -9.38117 -9.27088 -8.42352 -8.31025 -6.61307 -5.74549 -3.15853 2.47381 2.85741 3.40452 3.47306 3.48284 3.58977 3.92067 4.29230 4.45042 4.53512 4.58002 4.74611 4.81841 4.89869 4.97686 5.14265 5.15028 5.29154 5.35212 5.38633 5.42025 5.44302 5.50192 5.51993 5.57139 5.64559 5.67398 5.71562 5.74820 5.83259 5.84646 5.90326 5.92694 6.03779 6.05823 6.10833 6.22475 6.31522 6.37399 6.43371 6.45892 6.49829 6.59653 6.69912 6.85540 7.09326 7.38297 7.42430 7.61786 7.72539 7.83369 8.19876 8.49337 9.26125 9.90156 10.47783 11.33916 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.013882 B = 0.003515 C = 0.003235 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2016.532629 B = 7964.597555 C = 8653.635952 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.141 4.141 2 C 0.071 3.929 3 H 0.070 0.930 4 C 0.080 3.920 5 N -0.607 5.607 6 C 0.546 3.454 7 O -0.535 6.535 8 C -0.189 4.189 9 H 0.121 0.879 10 C -0.145 4.145 11 C -0.095 4.095 12 C -0.103 4.103 13 C -0.115 4.115 14 C 0.090 3.910 15 C 0.119 3.881 16 C -0.117 4.117 17 C 0.152 3.848 18 N -0.591 5.591 19 C -0.250 4.250 20 H 0.066 0.934 21 C 0.001 3.999 22 N -0.903 5.903 23 Si 0.901 3.099 24 H -0.288 1.288 25 H -0.292 1.292 26 H -0.279 1.279 27 C 0.181 3.819 28 C -0.159 4.159 29 O -0.356 6.356 30 H 0.061 0.939 31 H 0.071 0.929 32 H 0.067 0.933 33 H 0.073 0.927 34 H 0.090 0.910 35 H 0.104 0.896 36 H 0.101 0.899 37 H 0.066 0.934 38 H 0.047 0.953 39 H 0.065 0.935 40 H 0.058 0.942 41 H 0.065 0.935 42 H 0.060 0.940 43 H 0.097 0.903 44 H 0.068 0.932 45 H 0.059 0.941 46 H 0.073 0.927 47 H 0.052 0.948 48 H 0.027 0.973 49 H 0.252 0.748 50 H 0.084 0.916 51 H 0.021 0.979 52 H 0.062 0.938 53 H 0.048 0.952 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -3.018 -7.848 -18.972 20.752 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.199 4.199 2 C 0.012 3.988 3 H 0.087 0.913 4 C -0.043 4.043 5 N -0.340 5.340 6 C 0.334 3.666 7 O -0.412 6.412 8 C -0.210 4.210 9 H 0.139 0.861 10 C -0.183 4.183 11 C -0.095 4.095 12 C -0.160 4.160 13 C -0.172 4.172 14 C -0.033 4.033 15 C 0.079 3.921 16 C -0.157 4.157 17 C 0.025 3.975 18 N -0.314 5.314 19 C -0.380 4.380 20 H 0.083 0.917 21 C -0.197 4.197 22 N -0.547 5.547 23 Si 0.731 3.269 24 H -0.215 1.215 25 H -0.218 1.218 26 H -0.204 1.204 27 C 0.055 3.945 28 C -0.198 4.198 29 O -0.276 6.276 30 H 0.080 0.920 31 H 0.090 0.910 32 H 0.086 0.914 33 H 0.091 0.909 34 H 0.108 0.892 35 H 0.122 0.878 36 H 0.119 0.881 37 H 0.085 0.915 38 H 0.066 0.934 39 H 0.084 0.916 40 H 0.077 0.923 41 H 0.084 0.916 42 H 0.079 0.921 43 H 0.115 0.885 44 H 0.086 0.914 45 H 0.077 0.923 46 H 0.091 0.909 47 H 0.071 0.929 48 H 0.045 0.955 49 H 0.062 0.938 50 H 0.102 0.898 51 H 0.039 0.961 52 H 0.081 0.919 53 H 0.067 0.933 Dipole moment (debyes) X Y Z Total from point charges -2.532 -7.839 -18.733 20.464 hybrid contribution 0.472 0.607 0.380 0.858 sum -2.060 -7.232 -18.352 19.833 Atomic orbital electron populations 1.21692 0.92662 1.02680 1.02820 1.22376 0.95123 0.96368 0.84925 0.91258 1.22539 0.80192 1.01618 0.99998 1.48079 1.06792 1.62836 1.16281 1.19642 0.86182 0.77302 0.83450 1.90687 1.14548 1.55133 1.80843 1.22361 1.01885 1.00515 0.96259 0.86147 1.23246 1.03058 0.89246 1.02728 1.21602 0.94632 1.01542 0.91763 1.21354 0.90832 1.03069 1.00740 1.21389 1.01202 1.01852 0.92800 1.22459 1.00800 0.97007 0.83010 1.21709 0.80619 0.95968 0.93774 1.21997 1.01055 0.98225 0.94438 1.20787 0.96474 0.85944 0.94284 1.44034 1.70646 1.05748 1.10955 1.19211 1.23977 0.91324 1.03446 0.91675 1.25018 0.97467 0.98077 0.99157 1.70035 1.30415 1.08340 1.45948 0.86311 0.80597 0.79534 0.80475 1.21534 1.21818 1.20443 1.20635 0.96552 0.95200 0.82101 1.22694 1.01286 0.95257 1.00609 1.88001 1.20917 1.76363 1.42271 0.91998 0.90993 0.91384 0.90918 0.89201 0.87816 0.88072 0.91493 0.93367 0.91554 0.92301 0.91639 0.92104 0.88468 0.91410 0.92270 0.90891 0.92914 0.95480 0.93764 0.89790 0.96080 0.91903 0.93319 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 17. 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.60 9.57 37.15 0.36 -1.24 16 2 C 0.07 0.98 3.02 -26.58 -0.08 0.90 16 3 H 0.07 1.09 6.29 -51.93 -0.33 0.76 16 4 C 0.08 0.92 5.76 -3.53 -0.02 0.89 16 5 N -0.61 -8.44 2.53 -172.23 -0.44 -8.87 16 6 C 0.55 7.80 7.48 -10.99 -0.08 7.72 16 7 O -0.53 -9.18 12.58 5.56 0.07 -9.11 16 8 C -0.19 -2.04 4.67 -91.77 -0.43 -2.47 16 9 H 0.12 1.05 6.91 -51.92 -0.36 0.69 16 10 C -0.15 -1.45 9.53 -26.69 -0.25 -1.71 16 11 C -0.09 -0.96 2.70 -154.47 -0.42 -1.37 16 12 C -0.10 -1.20 8.09 37.16 0.30 -0.90 16 13 C -0.12 -0.86 9.85 37.16 0.37 -0.50 16 14 C 0.09 1.38 5.30 -3.51 -0.02 1.36 16 15 C 0.12 2.07 0.68 -90.49 -0.06 2.01 16 16 C -0.12 -2.09 5.10 -26.63 -0.14 -2.23 16 17 C 0.15 3.26 6.28 -3.72 -0.02 3.24 16 18 N -0.59 -15.05 2.99 -172.32 -0.52 -15.56 16 19 C -0.25 -6.95 9.59 -15.06 -0.14 -7.09 16 20 H 0.07 1.78 7.89 -52.48 -0.41 1.37 16 21 C 0.00 0.03 4.91 -17.43 -0.09 -0.05 16 22 N -0.90 -26.45 11.19 55.49 0.62 -25.83 16 23 Si 0.90 22.00 29.59 -169.99 -5.03 16.98 16 24 H -0.29 -7.05 7.11 56.52 0.40 -6.65 16 25 H -0.29 -7.18 7.11 56.52 0.40 -6.78 16 26 H -0.28 -6.67 7.11 56.52 0.40 -6.27 16 27 C 0.18 4.40 5.08 -3.71 -0.02 4.38 16 28 C -0.16 -3.17 4.33 -26.64 -0.12 -3.29 16 29 O -0.36 -6.14 8.64 -35.23 -0.30 -6.44 16 30 H 0.06 0.65 8.14 -51.93 -0.42 0.22 16 31 H 0.07 0.94 8.14 -51.93 -0.42 0.52 16 32 H 0.07 0.59 8.14 -51.93 -0.42 0.17 16 33 H 0.07 0.76 8.14 -51.93 -0.42 0.34 16 34 H 0.09 0.82 5.83 -51.93 -0.30 0.51 16 35 H 0.10 0.81 8.05 -51.93 -0.42 0.39 16 36 H 0.10 1.17 8.05 -51.93 -0.42 0.75 16 37 H 0.07 0.87 7.26 -51.93 -0.38 0.49 16 38 H 0.05 0.47 7.92 -51.93 -0.41 0.06 16 39 H 0.07 0.90 7.02 -51.93 -0.36 0.54 16 40 H 0.06 0.37 7.90 -51.93 -0.41 -0.04 16 41 H 0.06 0.45 7.92 -51.93 -0.41 0.03 16 42 H 0.06 0.44 8.14 -51.93 -0.42 0.01 16 43 H 0.10 1.62 7.08 -51.93 -0.37 1.25 16 44 H 0.07 0.96 8.14 -51.93 -0.42 0.54 16 45 H 0.06 0.92 8.14 -51.93 -0.42 0.49 16 46 H 0.07 1.30 8.14 -51.93 -0.42 0.87 16 47 H 0.05 1.10 7.93 -51.93 -0.41 0.69 16 48 H 0.03 0.61 7.82 -51.93 -0.41 0.20 16 49 H 0.25 7.15 8.31 -40.82 -0.34 6.82 16 50 H 0.08 2.31 6.03 -51.93 -0.31 1.99 16 51 H 0.02 0.54 7.83 -51.93 -0.41 0.14 16 52 H 0.06 1.33 8.14 -51.92 -0.42 0.91 16 53 H 0.05 0.93 6.16 -51.93 -0.32 0.61 16 LS Contribution 396.24 15.07 5.97 5.97 Total: -1.00 -31.73 396.24 -9.86 -41.59 By element: Atomic # 1 Polarization: 11.02 SS G_CDS: -9.38 Total: 1.64 kcal Atomic # 6 Polarization: 0.51 SS G_CDS: -0.86 Total: -0.36 kcal Atomic # 7 Polarization: -49.94 SS G_CDS: -0.33 Total: -50.27 kcal Atomic # 8 Polarization: -15.32 SS G_CDS: -0.23 Total: -15.55 kcal Atomic # 14 Polarization: 22.00 SS G_CDS: -5.03 Total: 16.98 kcal Total LS contribution 5.97 Total: 5.97 kcal Total: -31.73 -9.86 -41.59 kcal The number of atoms in the molecule is 53 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. ZINC001098586715.mol2 53 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 13.153 kcal (2) G-P(sol) polarization free energy of solvation -31.728 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -18.575 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.864 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -41.592 kcal (6) G-S(sol) free energy of system = (1) + (5) -28.440 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.27 seconds