Wall clock time and date at job start Sat Apr 11 2020 06:12:13 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.52996 * 1 3 3 C 1.53123 * 109.57495 * 2 1 4 4 C 1.53034 * 109.27992 * 177.60435 * 3 2 1 5 5 N 1.46847 * 109.51345 * 300.80438 * 4 3 2 6 6 C 1.46895 * 110.99982 * 185.84805 * 5 4 3 7 7 C 1.50700 * 109.47354 * 169.99549 * 6 5 4 8 8 O 1.21282 * 120.00036 * 0.02562 * 7 6 5 9 9 N 1.34778 * 120.00099 * 179.97438 * 7 6 5 10 10 C 1.46505 * 119.99960 * 180.02562 * 9 7 6 11 11 C 1.52997 * 109.47074 * 180.02562 * 10 9 7 12 12 N 1.46899 * 109.47323 * 179.97438 * 11 10 9 13 13 C 1.46900 * 110.99953 * 66.03848 * 12 11 10 14 14 C 1.46893 * 111.00249 * 189.99581 * 12 11 10 15 15 C 1.50701 * 109.47113 * 169.99899 * 14 12 11 16 16 O 1.21298 * 119.99877 * 359.97438 * 15 14 12 17 17 N 1.34775 * 120.00096 * 180.02562 * 15 14 12 18 18 C 1.37099 * 119.99964 * 179.97438 * 17 15 14 19 19 H 1.07996 * 120.00314 * 359.97438 * 18 17 15 20 20 C 1.38953 * 119.99768 * 180.02562 * 18 17 15 21 21 N 1.31413 * 120.00317 * 359.97438 * 20 18 17 22 22 Si 1.86802 * 119.99787 * 180.02562 * 20 18 17 23 23 H 1.48501 * 109.46790 * 90.00029 * 22 20 18 24 24 H 1.48497 * 109.47060 * 209.99478 * 22 20 18 25 25 H 1.48494 * 109.47272 * 329.99945 * 22 20 18 26 26 C 1.46843 * 111.20124 * 61.74903 * 5 4 3 27 27 C 1.53042 * 109.52278 * 298.25596 * 26 5 4 28 28 H 1.09000 * 109.46968 * 300.29446 * 1 2 3 29 29 H 1.09000 * 109.47430 * 60.29800 * 1 2 3 30 30 H 1.09002 * 109.46935 * 180.29882 * 1 2 3 31 31 H 1.09003 * 109.52281 * 239.83230 * 2 1 3 32 32 H 1.08997 * 109.51083 * 57.65029 * 3 2 1 33 33 H 1.08998 * 109.48980 * 297.53960 * 3 2 1 34 34 H 1.09002 * 109.46271 * 180.79158 * 4 3 2 35 35 H 1.08994 * 109.46242 * 60.81457 * 4 3 2 36 36 H 1.09000 * 109.47118 * 289.99436 * 6 5 4 37 37 H 1.09003 * 109.47534 * 49.99526 * 6 5 4 38 38 H 0.97001 * 120.00093 * 0.02562 * 9 7 6 39 39 H 1.09000 * 109.46742 * 300.00416 * 10 9 7 40 40 H 1.08993 * 109.47197 * 60.00027 * 10 9 7 41 41 H 1.08998 * 109.46708 * 299.99818 * 11 10 9 42 42 H 1.09003 * 109.47648 * 59.99696 * 11 10 9 43 43 H 1.09004 * 109.46775 * 60.00345 * 13 12 11 44 44 H 1.09000 * 109.47452 * 179.97438 * 13 12 11 45 45 H 1.09001 * 109.47645 * 300.00796 * 13 12 11 46 46 H 1.08994 * 109.47219 * 289.99809 * 14 12 11 47 47 H 1.09003 * 109.47051 * 50.00116 * 14 12 11 48 48 H 0.97002 * 119.99746 * 359.97438 * 17 15 14 49 49 H 0.96991 * 120.00307 * 180.02562 * 21 20 18 50 50 H 1.08993 * 109.46313 * 58.27251 * 26 5 4 51 51 H 1.09000 * 109.46480 * 178.25361 * 26 5 4 52 52 H 1.08997 * 109.50758 * 179.11703 * 27 26 5 53 53 H 1.08998 * 109.50897 * 299.20258 * 27 26 5 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.0430 1.4427 0.0000 4 6 3.5721 1.4353 0.0604 5 7 4.0116 0.7408 1.2773 6 6 5.4685 0.8344 1.4399 7 6 5.8548 0.3320 2.8071 8 8 5.0005 -0.0655 3.5708 9 7 7.1497 0.3233 3.1809 10 6 7.5252 -0.1657 4.5099 11 6 9.0427 -0.0680 4.6792 12 7 9.4192 -0.5577 6.0120 13 6 9.1448 -1.9950 6.1416 14 6 10.8286 -0.2626 6.3025 15 6 11.1138 -0.5526 7.7536 16 8 10.2297 -0.9642 8.4749 17 7 12.3508 -0.3536 8.2501 18 6 12.6101 -0.6170 9.5703 19 1 11.8230 -0.9834 10.2126 20 6 13.8854 -0.4113 10.0823 21 7 14.8431 0.0351 9.3010 22 14 14.2388 -0.7708 11.8810 23 1 14.6775 -2.1821 12.0268 24 1 15.3125 0.1337 12.3649 25 1 13.0075 -0.5500 12.6812 26 6 3.5710 -0.6600 1.2779 27 6 2.0417 -0.7151 1.2536 28 1 -0.3633 0.5184 0.8873 29 1 -0.3634 0.5092 -0.8926 30 1 -0.3633 -1.0277 0.0054 31 1 1.8942 -0.5163 -0.8882 32 1 1.6470 1.9699 0.8679 33 1 1.7172 1.9445 -0.9111 34 1 3.9400 2.4612 0.0747 35 1 3.9676 0.9199 -0.8148 36 1 5.9588 0.2282 0.6781 37 1 5.7800 1.8735 1.3333 38 1 7.8330 0.6409 2.5700 39 1 7.0352 0.4404 5.2718 40 1 7.2136 -1.2047 4.6162 41 1 9.5329 -0.6739 3.9172 42 1 9.3541 0.9712 4.5730 43 1 9.7172 -2.5432 5.3932 44 1 9.4332 -2.3317 7.1374 45 1 8.0808 -2.1770 5.9901 46 1 11.4669 -0.8849 5.6754 47 1 11.0298 0.7886 6.0957 48 1 13.0577 -0.0241 7.6733 49 1 15.7333 0.1783 9.6583 50 1 3.9666 -1.1655 0.3970 51 1 3.9378 -1.1558 2.1766 52 1 1.7152 -1.7549 1.2381 53 1 1.6456 -0.2224 2.1415 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 ZINC001663304493.mol2 53 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED SCF FIELD WAS ACHIEVED Wall clock time and date at job start Sat Apr 11 2020 06:12:13 Heat of formation + Delta-G solvation = -67.122747 kcal Electronic energy + Delta-G solvation = -29082.913498 eV Core-core repulsion = 24989.333330 eV Total energy + Delta-G solvation = -4093.580168 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.226 amu Computer time = 0.75 seconds Orbital eigenvalues (eV) -40.31234 -38.98278 -38.55691 -36.78966 -34.93886 -34.56848 -33.80265 -31.25689 -31.16074 -30.37050 -28.97543 -27.70502 -26.25924 -25.08348 -23.69201 -23.51815 -22.17807 -21.23603 -20.50516 -20.00500 -19.16975 -17.85434 -17.26561 -16.68238 -16.15172 -16.04625 -15.69417 -15.47608 -14.98741 -14.81077 -14.72747 -14.60906 -14.47491 -14.21217 -14.03016 -13.57465 -13.46622 -13.32848 -13.21686 -12.94641 -12.72916 -12.43917 -12.17164 -12.12991 -11.91208 -11.81737 -11.52936 -11.50260 -11.39297 -11.27620 -11.06197 -11.03854 -10.94314 -10.68577 -10.62540 -10.38626 -10.21771 -9.85813 -9.47801 -9.17395 -8.64430 -8.48143 -7.98961 -7.89179 -6.41832 -3.82502 2.49062 2.76411 3.18070 3.26915 3.33127 3.59281 3.89984 4.05476 4.26570 4.38302 4.41916 4.46028 4.53687 4.80376 4.83223 4.88589 4.91535 4.99397 5.03901 5.10090 5.12381 5.17086 5.24680 5.28171 5.30946 5.37986 5.39049 5.49888 5.62436 5.63605 5.68394 5.68787 5.78177 5.84152 5.85184 5.88621 5.91018 5.92881 6.01282 6.11459 6.25635 6.29195 6.55381 6.60834 6.74881 7.19048 7.30444 7.35183 7.55501 8.06443 8.09772 8.61815 9.18334 9.54234 10.06190 10.75946 Molecular weight = 340.23amu Principal moments of inertia in cm(-1) A = 0.041062 B = 0.001408 C = 0.001385 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 681.734281 B =19880.398754 C =20217.415367 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.146 4.146 2 C -0.096 4.096 3 C -0.111 4.111 4 C 0.056 3.944 5 N -0.520 5.520 6 C 0.037 3.963 7 C 0.512 3.488 8 O -0.518 6.518 9 N -0.732 5.732 10 C 0.123 3.877 11 C 0.045 3.955 12 N -0.515 5.515 13 C 0.018 3.982 14 C 0.047 3.953 15 C 0.443 3.557 16 O -0.576 6.576 17 N -0.564 5.564 18 C -0.299 4.299 19 H 0.103 0.897 20 C 0.023 3.977 21 N -0.903 5.903 22 Si 0.930 3.070 23 H -0.278 1.278 24 H -0.284 1.284 25 H -0.269 1.269 26 C 0.050 3.950 27 C -0.109 4.109 28 H 0.056 0.944 29 H 0.052 0.948 30 H 0.053 0.947 31 H 0.064 0.936 32 H 0.073 0.927 33 H 0.066 0.934 34 H 0.075 0.925 35 H 0.032 0.968 36 H 0.062 0.938 37 H 0.104 0.896 38 H 0.398 0.602 39 H 0.079 0.921 40 H 0.069 0.931 41 H 0.035 0.965 42 H 0.079 0.921 43 H 0.015 0.985 44 H 0.098 0.902 45 H 0.063 0.937 46 H 0.048 0.952 47 H 0.089 0.911 48 H 0.392 0.608 49 H 0.271 0.729 50 H 0.025 0.975 51 H 0.095 0.905 52 H 0.065 0.935 53 H 0.074 0.926 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -19.009 2.368 -24.936 31.445 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.203 4.203 2 C -0.114 4.114 3 C -0.149 4.149 4 C -0.071 4.071 5 N -0.245 5.245 6 C -0.094 4.094 7 C 0.297 3.703 8 O -0.393 6.393 9 N -0.385 5.385 10 C 0.000 4.000 11 C -0.084 4.084 12 N -0.239 5.239 13 C -0.129 4.129 14 C -0.084 4.084 15 C 0.229 3.771 16 O -0.460 6.460 17 N -0.201 5.201 18 C -0.428 4.428 19 H 0.121 0.879 20 C -0.179 4.179 21 N -0.544 5.544 22 Si 0.757 3.243 23 H -0.204 1.204 24 H -0.209 1.209 25 H -0.194 1.194 26 C -0.077 4.077 27 C -0.147 4.147 28 H 0.075 0.925 29 H 0.071 0.929 30 H 0.073 0.927 31 H 0.082 0.918 32 H 0.092 0.908 33 H 0.085 0.915 34 H 0.094 0.906 35 H 0.050 0.950 36 H 0.080 0.920 37 H 0.122 0.878 38 H 0.231 0.769 39 H 0.098 0.902 40 H 0.087 0.913 41 H 0.054 0.946 42 H 0.097 0.903 43 H 0.033 0.967 44 H 0.116 0.884 45 H 0.082 0.918 46 H 0.066 0.934 47 H 0.107 0.893 48 H 0.225 0.775 49 H 0.081 0.919 50 H 0.043 0.957 51 H 0.113 0.887 52 H 0.083 0.917 53 H 0.092 0.908 Dipole moment (debyes) X Y Z Total from point charges -20.545 2.313 -23.069 30.978 hybrid contribution 2.130 -0.535 -1.153 2.481 sum -18.414 1.778 -24.222 30.478 Atomic orbital electron populations 1.21806 0.95067 1.01690 1.01766 1.21214 0.95718 0.96174 0.98324 1.21636 0.94238 0.96704 1.02272 1.22118 0.96330 0.98751 0.89944 1.61566 1.09745 1.15579 1.37655 1.21940 0.90128 1.02876 0.94504 1.20458 0.82175 0.78586 0.89098 1.90704 1.48997 1.45819 1.53764 1.45944 1.08029 1.66658 1.17909 1.21314 0.94842 1.00422 0.83449 1.22251 0.94513 1.00988 0.90656 1.61437 1.19999 1.17262 1.25232 1.22625 1.00483 0.86336 1.03492 1.21864 0.88156 1.01696 0.96678 1.19711 0.83728 0.85022 0.88684 1.90497 1.45101 1.50326 1.60078 1.41857 1.09027 1.58120 1.11072 1.19709 0.96201 1.41150 0.85763 0.87882 1.25252 0.95279 0.97154 1.00208 1.69607 1.02696 1.45791 1.36303 0.85860 0.78236 0.78138 0.82087 1.20446 1.20938 1.19381 1.22281 0.96133 0.87869 1.01445 1.21611 0.94545 1.00228 0.98279 0.92481 0.92859 0.92746 0.91759 0.90832 0.91540 0.90613 0.94955 0.92049 0.87783 0.76866 0.90236 0.91269 0.94650 0.90292 0.96675 0.88374 0.91849 0.93381 0.89319 0.77516 0.91898 0.95657 0.88702 0.91666 0.90763 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 13. 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.15 -0.84 9.20 37.16 0.34 -0.50 16 2 C -0.10 -0.63 2.79 -90.49 -0.25 -0.88 16 3 C -0.11 -0.71 5.38 -26.65 -0.14 -0.85 16 4 C 0.06 0.38 5.46 -3.84 -0.02 0.35 16 5 N -0.52 -4.66 4.40 -236.12 -1.04 -5.69 16 6 C 0.04 0.28 5.45 -4.98 -0.03 0.25 16 7 C 0.51 5.49 7.50 -10.98 -0.08 5.41 16 8 O -0.52 -7.81 14.32 5.56 0.08 -7.73 16 9 N -0.73 -6.64 5.56 -60.29 -0.33 -6.98 16 10 C 0.12 1.45 5.60 -4.04 -0.02 1.43 16 11 C 0.04 0.52 5.07 -3.83 -0.02 0.50 16 12 N -0.52 -8.38 4.50 -244.38 -1.10 -9.48 16 13 C 0.02 0.30 8.96 60.06 0.54 0.84 16 14 C 0.05 0.82 5.45 -4.98 -0.03 0.80 16 15 C 0.44 10.36 7.47 -10.99 -0.08 10.27 16 16 O -0.58 -15.46 13.91 5.54 0.08 -15.38 16 17 N -0.56 -14.15 5.29 -10.96 -0.06 -14.21 16 18 C -0.30 -8.29 9.68 -14.93 -0.14 -8.43 16 19 H 0.10 3.00 7.16 -52.49 -0.38 2.63 16 20 C 0.02 0.60 5.50 -17.47 -0.10 0.51 16 21 N -0.90 -24.38 13.05 55.50 0.72 -23.66 16 22 Si 0.93 21.00 29.55 -169.99 -5.02 15.98 16 23 H -0.28 -6.21 7.11 56.52 0.40 -5.81 16 24 H -0.28 -6.28 7.11 56.52 0.40 -5.88 16 25 H -0.27 -6.10 7.11 56.52 0.40 -5.70 16 26 C 0.05 0.46 5.11 -3.83 -0.02 0.44 16 27 C -0.11 -0.89 5.38 -26.64 -0.14 -1.04 16 28 H 0.06 0.35 8.14 -51.93 -0.42 -0.08 16 29 H 0.05 0.27 8.14 -51.93 -0.42 -0.16 16 30 H 0.05 0.30 8.14 -51.93 -0.42 -0.12 16 31 H 0.06 0.42 8.14 -51.93 -0.42 -0.01 16 32 H 0.07 0.55 8.14 -51.93 -0.42 0.13 16 33 H 0.07 0.34 8.14 -51.93 -0.42 -0.08 16 34 H 0.08 0.43 8.02 -51.93 -0.42 0.02 16 35 H 0.03 0.20 8.14 -51.93 -0.42 -0.22 16 36 H 0.06 0.36 8.12 -51.93 -0.42 -0.07 16 37 H 0.10 0.58 7.99 -51.93 -0.41 0.17 16 38 H 0.40 2.33 8.47 -40.82 -0.35 1.98 16 39 H 0.08 1.12 8.14 -51.93 -0.42 0.70 16 40 H 0.07 0.90 6.55 -51.93 -0.34 0.55 16 41 H 0.04 0.37 8.12 -51.93 -0.42 -0.05 16 42 H 0.08 0.83 7.93 -51.93 -0.41 0.41 16 43 H 0.01 0.22 8.14 -51.93 -0.42 -0.21 16 44 H 0.10 1.94 5.82 -51.93 -0.30 1.64 16 45 H 0.06 0.93 6.47 -51.93 -0.34 0.59 16 46 H 0.05 0.76 8.12 -51.93 -0.42 0.34 16 47 H 0.09 1.37 7.93 -51.93 -0.41 0.96 16 48 H 0.39 9.36 7.90 -40.82 -0.32 9.04 16 49 H 0.27 6.94 8.31 -40.82 -0.34 6.60 16 50 H 0.03 0.21 8.14 -51.93 -0.42 -0.21 16 51 H 0.09 1.11 5.78 -51.93 -0.30 0.81 16 52 H 0.06 0.49 8.14 -51.93 -0.42 0.06 16 53 H 0.07 0.72 8.14 -51.93 -0.42 0.29 16 LS Contribution 416.30 15.07 6.27 6.27 Total: -1.00 -33.41 416.30 -10.05 -43.46 By element: Atomic # 1 Polarization: 17.77 SS G_CDS: -9.45 Total: 8.32 kcal Atomic # 6 Polarization: 9.30 SS G_CDS: -0.20 Total: 9.10 kcal Atomic # 7 Polarization: -58.21 SS G_CDS: -1.81 Total: -60.02 kcal Atomic # 8 Polarization: -23.27 SS G_CDS: 0.16 Total: -23.12 kcal Atomic # 14 Polarization: 21.00 SS G_CDS: -5.02 Total: 15.98 kcal Total LS contribution 6.27 Total: 6.27 kcal Total: -33.41 -10.05 -43.46 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. ZINC001663304493.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 -23.665 kcal (2) G-P(sol) polarization free energy of solvation -33.407 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -57.072 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.051 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -43.458 kcal (6) G-S(sol) free energy of system = (1) + (5) -67.123 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.76 seconds