Wall clock time and date at job start Sat Mar 6 2021 20:15:25 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= 0 * AM1 - THE AM1 HAMILTONIAN TO BE USED * SM5.42R - SM5.42R CALCULATIONS WILL BE PERFORMED * - CHARGE MODEL 2 WILL BE USED (CM2A) * - THE SOLVENT IS WATER ******************************************************************************* 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.52995 * 1 3 3 C 1.52997 * 109.47209 * 2 1 4 4 C 1.53003 * 109.46902 * 239.99646 * 2 1 3 5 5 Si 1.86306 * 109.47294 * 120.00201 * 2 1 3 6 6 C 1.86305 * 109.46912 * 299.99840 * 5 2 1 7 7 C 1.86298 * 109.47317 * 179.97438 * 5 2 1 8 8 O 1.63099 * 109.46988 * 59.99474 * 5 2 1 9 9 C 1.42905 * 113.99927 * 180.02562 * 8 5 2 10 10 C 1.53041 * 109.49402 * 150.00366 * 9 8 5 11 11 C 1.53049 * 109.53178 * 181.31786 * 10 9 8 12 12 C 1.53183 * 109.31217 * 298.63202 * 11 10 9 13 13 N 1.46923 * 108.89650 * 54.56840 * 12 11 10 14 14 C 1.34776 * 120.63833 * 126.45820 * 13 12 11 15 15 O 1.21650 * 120.00764 * 175.73525 * 14 13 12 16 16 C 1.47050 * 119.99407 * 355.74182 * 14 13 12 17 17 C 1.36224 * 123.49000 * 232.47781 * 16 14 13 18 18 Cl 1.73599 * 124.60490 * 0.02812 * 17 16 14 19 19 S 1.70443 * 110.78836 * 179.97438 * 17 16 14 20 20 C 1.75493 * 92.46745 * 0.02562 * 19 17 16 21 21 C 1.50702 * 124.68584 * 179.97438 * 20 19 17 22 22 O 1.42896 * 109.47161 * 244.97146 * 21 20 19 23 23 Si 1.86296 * 109.47181 * 4.97401 * 21 20 19 24 24 H 1.48507 * 109.47246 * 59.99916 * 23 21 20 25 25 H 1.48505 * 109.47401 * 179.97438 * 23 21 20 26 26 H 1.48503 * 109.47085 * 299.99833 * 23 21 20 27 27 C 1.32690 * 110.63014 * 359.70624 * 20 19 17 28 28 C 1.46968 * 118.71228 * 306.64057 * 13 12 11 29 29 H 1.09004 * 109.46865 * 299.99828 * 1 2 3 30 30 H 1.08998 * 109.47201 * 60.00058 * 1 2 3 31 31 H 1.09002 * 109.47474 * 180.02562 * 1 2 3 32 32 H 1.09007 * 109.47104 * 180.02562 * 3 2 1 33 33 H 1.08998 * 109.46992 * 300.00071 * 3 2 1 34 34 H 1.08995 * 109.47063 * 60.00590 * 3 2 1 35 35 H 1.08990 * 109.47500 * 300.00498 * 4 2 1 36 36 H 1.08998 * 109.46810 * 60.00680 * 4 2 1 37 37 H 1.09000 * 109.46724 * 179.97438 * 4 2 1 38 38 H 1.08998 * 109.47030 * 179.97438 * 6 5 2 39 39 H 1.08993 * 109.46857 * 300.00035 * 6 5 2 40 40 H 1.09000 * 109.46832 * 60.00107 * 6 5 2 41 41 H 1.09004 * 109.46905 * 300.00010 * 7 5 2 42 42 H 1.08993 * 109.47208 * 59.99502 * 7 5 2 43 43 H 1.08998 * 109.46855 * 180.02562 * 7 5 2 44 44 H 1.09000 * 109.49529 * 29.93597 * 9 8 5 45 45 H 1.08996 * 109.46476 * 61.29918 * 10 9 8 46 46 H 1.08999 * 109.46178 * 301.33398 * 10 9 8 47 47 H 1.09007 * 109.49167 * 178.68010 * 11 10 9 48 48 H 1.09000 * 109.49465 * 58.58902 * 11 10 9 49 49 H 1.09007 * 109.58735 * 174.42455 * 12 11 10 50 50 H 1.08996 * 109.59039 * 294.70474 * 12 11 10 51 51 H 1.09002 * 109.46915 * 124.97258 * 21 20 19 52 52 H 0.96698 * 114.00198 * 60.00014 * 22 21 20 53 53 H 1.07998 * 123.45700 * 180.24894 * 27 20 19 54 54 H 1.08997 * 109.46595 * 293.56743 * 28 13 12 55 55 H 1.09006 * 109.59336 * 173.21457 * 28 13 12 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.0399 -0.7214 -1.2492 5 14 2.1510 -0.8783 1.5211 6 6 1.5300 0.0000 3.0423 7 6 4.0140 -0.8790 1.5207 8 8 1.6073 -2.4160 1.5211 9 6 2.0285 -3.1814 2.6520 10 6 0.9700 -4.2382 2.9763 11 6 1.4049 -5.0369 4.2073 12 6 2.7257 -5.7520 3.9066 13 7 3.7105 -4.7610 3.4518 14 6 4.9138 -4.6635 4.0511 15 8 5.7492 -3.8933 3.6166 16 6 5.2169 -5.4884 5.2301 17 6 5.6829 -4.9788 6.4043 18 17 5.9899 -3.2951 6.6956 19 16 5.9224 -6.2186 7.5492 20 6 5.3856 -7.4594 6.4302 21 6 5.3143 -8.9303 6.7503 22 8 6.2453 -9.6415 5.9322 23 14 5.7355 -9.1962 8.5455 24 1 4.7680 -8.4571 9.3957 25 1 5.6647 -10.6456 8.8612 26 1 7.1091 -8.6979 8.8107 27 6 5.0456 -6.9229 5.2651 28 6 3.3565 -3.8755 2.3336 29 1 -0.3633 0.5138 0.8900 30 1 -0.3633 0.5138 -0.8900 31 1 -0.3634 -1.0277 0.0005 32 1 3.1300 1.4424 -0.0005 33 1 1.6766 1.9563 -0.8900 34 1 1.6767 1.9562 0.8900 35 1 1.6766 -1.7489 -1.2493 36 1 1.6766 -0.2075 -2.1392 37 1 3.1299 -0.7209 -1.2494 38 1 1.8930 -0.5141 3.9323 39 1 1.8934 1.0275 3.0423 40 1 0.4400 -0.0001 3.0423 41 1 4.3769 -1.3931 0.6307 42 1 4.3777 0.1485 1.5205 43 1 4.3773 -1.3932 2.4105 44 1 2.1588 -2.5216 3.5098 45 1 0.8576 -4.9122 2.1272 46 1 0.0178 -3.7481 3.1792 47 1 0.6396 -5.7737 4.4517 48 1 1.5413 -4.3600 5.0506 49 1 3.0902 -6.2414 4.8098 50 1 2.5689 -6.4962 3.1258 51 1 4.3061 -9.2960 6.5557 52 1 7.1647 -9.3674 6.0531 53 1 4.6818 -7.4995 4.4276 54 1 3.2516 -4.4649 1.4227 55 1 4.1371 -3.1267 2.1990 RHF calculation, no. of doubly occupied orbitals= 71 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). S: (AM1): M.J.S. DEWAR, Y.-C. YUAN, INORGANIC CHEM., 29, 589 (1990) Cl: (AM1): M.J.S. DEWAR AND E. G. ZOEBISCH, THEOCHEM, 180, 1 (1988). 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=0 AM1 1SCF TLIMIT=15 GEO-OK SM5.42R & SOLVNT=WATER s_22_10954526_13431040.mol2 55 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Sat Mar 6 2021 20:15:25 Heat of formation + Delta-G solvation = -129.751824 kcal Electronic energy + Delta-G solvation = -35438.154171 eV Core-core repulsion = 30944.185962 eV Total energy + Delta-G solvation = -4493.968209 eV No. of doubly occupied orbitals = 71 Molecular weight (most abundant/longest-lived isotopes) = 419.131 amu Computer time = 1.36 seconds Orbital eigenvalues (eV) -41.66673 -39.57268 -38.78779 -38.09562 -37.14219 -36.87952 -36.82724 -33.47284 -31.87912 -31.14751 -29.37224 -29.16014 -27.89831 -27.50471 -27.16846 -26.95680 -25.02822 -23.88248 -22.20634 -21.34881 -20.42774 -20.12444 -19.31476 -18.15447 -17.96988 -17.71738 -17.34032 -16.69939 -16.26753 -16.09324 -15.96113 -15.91771 -15.52737 -15.43704 -14.93792 -14.80041 -14.77151 -14.71147 -14.70150 -14.32137 -14.23287 -14.12920 -13.94745 -13.84275 -13.66012 -13.43974 -13.22311 -13.06386 -12.98087 -12.95372 -12.93487 -12.71479 -12.46699 -12.42006 -12.27292 -12.25201 -12.10422 -12.04602 -11.91920 -11.79699 -11.72134 -11.66878 -11.60050 -11.15991 -11.09404 -10.56236 -10.51734 -10.10693 -9.96886 -9.79101 -9.62692 -0.61967 -0.25411 0.40062 0.78256 0.93822 1.09245 1.28001 1.41762 1.49695 1.63900 1.72254 1.91852 1.95408 2.00175 2.40845 3.21683 3.50119 3.66225 3.75406 3.79345 3.84527 3.91954 3.96972 4.09605 4.10428 4.12749 4.35119 4.42940 4.43765 4.49878 4.58108 4.61329 4.62167 4.66650 4.68272 4.70926 4.75842 4.79682 4.80562 4.82950 4.89394 4.92422 4.98536 5.01336 5.07750 5.07833 5.08491 5.18489 5.20442 5.32256 5.37749 5.38767 5.43437 5.51806 5.56201 5.62878 6.16305 6.63434 6.90285 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.137 4.137 2 C -0.374 4.374 3 C -0.143 4.143 4 C -0.138 4.138 5 Si 1.308 2.692 6 C -0.516 4.516 7 C -0.522 4.522 8 O -0.618 6.618 9 C 0.104 3.896 10 C -0.120 4.120 11 C -0.132 4.132 12 C 0.087 3.913 13 N -0.594 5.594 14 C 0.575 3.425 15 O -0.562 6.562 16 C -0.150 4.150 17 C -0.146 4.146 18 Cl -0.008 7.008 19 S 0.233 5.767 20 C -0.217 4.217 21 C -0.027 4.027 22 O -0.542 6.542 23 Si 0.897 3.103 24 H -0.233 1.233 25 H -0.235 1.235 26 H -0.238 1.238 27 C -0.068 4.068 28 C 0.071 3.929 29 H 0.060 0.940 30 H 0.052 0.948 31 H 0.060 0.940 32 H 0.059 0.941 33 H 0.055 0.945 34 H 0.062 0.938 35 H 0.060 0.940 36 H 0.052 0.948 37 H 0.055 0.945 38 H 0.062 0.938 39 H 0.077 0.923 40 H 0.077 0.923 41 H 0.070 0.930 42 H 0.075 0.925 43 H 0.057 0.943 44 H 0.059 0.941 45 H 0.084 0.916 46 H 0.081 0.919 47 H 0.112 0.888 48 H 0.063 0.937 49 H 0.105 0.895 50 H 0.098 0.902 51 H 0.155 0.845 52 H 0.391 0.609 53 H 0.182 0.818 54 H 0.091 0.909 55 H 0.089 0.911 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -4.550 -5.423 3.187 7.763 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.194 4.194 2 C -0.418 4.418 3 C -0.200 4.200 4 C -0.195 4.195 5 Si 1.438 2.562 6 C -0.618 4.618 7 C -0.623 4.623 8 O -0.576 6.576 9 C 0.046 3.954 10 C -0.157 4.157 11 C -0.170 4.170 12 C -0.034 4.034 13 N -0.327 5.327 14 C 0.362 3.638 15 O -0.441 6.441 16 C -0.170 4.170 17 C -0.292 4.292 18 Cl 0.021 6.979 19 S 0.500 5.500 20 C -0.335 4.335 21 C -0.141 4.141 22 O -0.349 6.349 23 Si 0.718 3.282 24 H -0.156 1.156 25 H -0.158 1.158 26 H -0.161 1.161 27 C -0.096 4.096 28 C -0.052 4.052 29 H 0.079 0.921 30 H 0.071 0.929 31 H 0.079 0.921 32 H 0.078 0.922 33 H 0.074 0.926 34 H 0.081 0.919 35 H 0.079 0.921 36 H 0.071 0.929 37 H 0.075 0.925 38 H 0.081 0.919 39 H 0.096 0.904 40 H 0.097 0.903 41 H 0.089 0.911 42 H 0.094 0.906 43 H 0.076 0.924 44 H 0.077 0.923 45 H 0.103 0.897 46 H 0.100 0.900 47 H 0.130 0.870 48 H 0.082 0.918 49 H 0.123 0.877 50 H 0.116 0.884 51 H 0.172 0.828 52 H 0.241 0.759 53 H 0.199 0.801 54 H 0.109 0.891 55 H 0.107 0.893 Dipole moment (debyes) X Y Z Total from point charges -4.060 -5.486 5.011 8.467 hybrid contribution 0.577 1.722 -1.718 2.500 sum -3.483 -3.764 3.293 6.095 Atomic orbital electron populations 1.21545 0.93679 1.02118 1.02024 1.25771 1.00790 1.02754 1.12476 1.21619 1.01109 0.94825 1.02423 1.21549 1.00821 1.00276 0.96827 0.76388 0.61693 0.53414 0.64729 1.27088 1.07292 1.09716 1.17696 1.27408 1.26321 1.04607 1.03916 1.85414 1.78341 1.46073 1.47791 1.21901 0.93623 0.93545 0.86333 1.21724 0.97961 0.96974 0.99080 1.21965 0.97198 1.00400 0.97439 1.22394 0.88027 0.90582 1.02359 1.47978 1.14587 1.36578 1.33598 1.17819 0.82028 0.81623 0.82346 1.90838 1.47738 1.40145 1.65385 1.20244 1.07769 0.95313 0.93666 1.25427 1.09455 0.93243 1.01106 1.98417 1.93915 1.10141 1.95418 1.83402 1.62152 0.92131 1.12360 1.25013 1.09334 1.01074 0.98056 1.24398 0.99318 0.84148 1.06242 1.86538 1.20953 1.61793 1.65621 0.86197 0.82851 0.84016 0.75116 1.15644 1.15783 1.16093 1.21110 0.99977 0.89451 0.99097 1.22498 0.97825 0.95151 0.89689 0.92119 0.92940 0.92081 0.92244 0.92599 0.91855 0.92105 0.92888 0.92550 0.91898 0.90367 0.90348 0.91118 0.90603 0.92380 0.92334 0.89695 0.90014 0.87009 0.91835 0.87686 0.88385 0.82780 0.75929 0.80065 0.89094 0.89271 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 15. 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 -0.06 7.85 71.98 0.56 0.50 16 2 C -0.37 -0.29 0.53 -10.79 -0.01 -0.29 16 3 C -0.14 -0.06 7.85 71.98 0.56 0.51 16 4 C -0.14 -0.20 7.85 71.98 0.56 0.37 16 5 Si 1.31 0.76 8.07 68.60 0.55 1.31 16 6 C -0.52 -0.53 8.56 113.37 0.97 0.44 16 7 C -0.52 -2.24 7.53 113.37 0.85 -1.38 16 8 O -0.62 -0.83 8.28 -148.98 -1.23 -2.07 16 9 C 0.10 0.05 2.30 30.67 0.07 0.12 16 10 C -0.12 0.48 5.54 30.62 0.17 0.65 16 11 C -0.13 0.91 6.08 30.67 0.19 1.09 16 12 C 0.09 -0.50 5.75 86.36 0.50 -0.01 16 13 N -0.59 -1.08 2.98 -823.11 -2.45 -3.53 16 14 C 0.58 4.30 7.04 86.55 0.61 4.91 16 15 O -0.56 -8.85 16.12 -4.19 -0.07 -8.92 16 16 C -0.15 -0.47 4.95 -19.84 -0.10 -0.57 16 17 C -0.15 -0.57 6.87 67.35 0.46 -0.11 16 18 Cl -0.01 -0.05 28.46 -2.72 -0.08 -0.13 16 19 S 0.23 -0.36 20.19 -56.49 -1.14 -1.50 16 20 C -0.22 1.04 5.65 24.52 0.14 1.18 16 21 C -0.03 0.21 2.94 71.24 0.21 0.42 16 22 O -0.54 3.05 12.88 -148.98 -1.92 1.14 16 23 Si 0.90 -5.30 26.20 68.60 1.80 -3.50 16 24 H -0.23 1.21 7.11 99.48 0.71 1.91 16 25 H -0.24 0.78 7.11 99.48 0.71 1.49 16 26 H -0.24 1.12 7.11 99.48 0.71 1.83 16 27 C -0.07 0.23 7.82 22.70 0.18 0.40 16 28 C 0.07 0.20 5.75 86.35 0.50 0.70 16 29 H 0.06 -0.05 7.79 -2.39 -0.02 -0.07 16 30 H 0.05 0.05 8.14 -2.39 -0.02 0.03 16 31 H 0.06 0.02 7.96 -2.39 -0.02 0.00 16 32 H 0.06 -0.01 7.80 -2.38 -0.02 -0.03 16 33 H 0.05 0.04 8.14 -2.39 -0.02 0.02 16 34 H 0.06 -0.07 7.79 -2.39 -0.02 -0.09 16 35 H 0.06 0.08 7.96 -2.39 -0.02 0.06 16 36 H 0.05 0.09 8.14 -2.39 -0.02 0.07 16 37 H 0.06 0.06 7.79 -2.39 -0.02 0.04 16 38 H 0.06 0.11 7.10 -2.39 -0.02 0.09 16 39 H 0.08 0.04 7.79 -2.39 -0.02 0.02 16 40 H 0.08 0.01 7.79 -2.39 -0.02 -0.01 16 41 H 0.07 0.32 7.79 -2.39 -0.02 0.30 16 42 H 0.07 0.26 7.79 -2.39 -0.02 0.24 16 43 H 0.06 0.44 6.34 -2.39 -0.02 0.42 16 44 H 0.06 0.06 6.15 -2.39 -0.01 0.04 16 45 H 0.08 -0.41 8.14 -2.39 -0.02 -0.43 16 46 H 0.08 -0.37 8.14 -2.39 -0.02 -0.38 16 47 H 0.11 -1.25 8.14 -2.38 -0.02 -1.27 16 48 H 0.06 -0.32 8.14 -2.39 -0.02 -0.33 16 49 H 0.11 -0.80 4.63 -2.38 -0.01 -0.81 16 50 H 0.10 -0.90 8.14 -2.39 -0.02 -0.92 16 51 H 0.15 -1.85 7.96 -2.39 -0.02 -1.86 16 52 H 0.39 -4.21 8.70 -74.06 -0.64 -4.85 16 53 H 0.18 -1.32 7.59 -2.91 -0.02 -1.34 16 54 H 0.09 0.10 8.14 -2.39 -0.02 0.08 16 55 H 0.09 0.63 3.84 -2.38 -0.01 0.62 16 Total: 0.00 -16.31 449.21 2.90 -13.41 By element: Atomic # 1 Polarization: -6.15 SS G_CDS: 1.01 Total: -5.14 kcal Atomic # 6 Polarization: 2.50 SS G_CDS: 6.43 Total: 8.93 kcal Atomic # 7 Polarization: -1.08 SS G_CDS: -2.45 Total: -3.53 kcal Atomic # 8 Polarization: -6.63 SS G_CDS: -3.22 Total: -9.85 kcal Atomic # 14 Polarization: -4.54 SS G_CDS: 2.35 Total: -2.19 kcal Atomic # 16 Polarization: -0.36 SS G_CDS: -1.14 Total: -1.50 kcal Atomic # 17 Polarization: -0.05 SS G_CDS: -0.08 Total: -0.13 kcal Total: -16.31 2.90 -13.41 kcal The number of atoms in the molecule is 55 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. s_22_10954526_13431040.mol2 55 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 -116.337 kcal (2) G-P(sol) polarization free energy of solvation -16.314 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -132.651 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy 2.899 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -13.415 kcal (6) G-S(sol) free energy of system = (1) + (5) -129.752 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.36 seconds