Wall clock time and date at job start Sat Apr 11 2020 06:10:11 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 O 1.42898 * 1 3 3 C 1.42905 * 113.99709 * 2 1 4 4 C 1.52995 * 109.47035 * 180.02562 * 3 2 1 5 5 C 1.52994 * 109.47374 * 180.02562 * 4 3 2 6 6 C 1.53003 * 109.47355 * 179.97438 * 5 4 3 7 7 C 1.50700 * 109.46866 * 180.02562 * 6 5 4 8 8 O 1.21282 * 120.00414 * 359.97438 * 7 6 5 9 9 N 1.34785 * 119.99662 * 180.02562 * 7 6 5 10 10 C 1.47155 * 121.05608 * 0.02562 * 9 7 6 11 11 C 1.53291 * 108.49670 * 128.78730 * 10 9 7 12 12 H 1.09007 * 109.50920 * 293.06245 * 11 10 9 13 13 C 1.52997 * 109.51784 * 172.95997 * 11 10 9 14 14 O 1.43087 * 109.21895 * 52.98258 * 11 10 9 15 15 C 1.42062 * 114.40171 * 298.85720 * 14 11 10 16 16 C 1.54525 * 110.32743 * 182.55269 * 15 14 11 17 17 C 1.55022 * 101.56847 * 83.95127 * 16 15 14 18 18 N 1.47493 * 104.89543 * 36.35515 * 17 16 15 19 19 C 1.36937 * 125.65696 * 155.94471 * 18 17 16 20 20 H 1.08001 * 119.99255 * 180.02562 * 19 18 17 21 21 C 1.38811 * 120.00434 * 0.03058 * 19 18 17 22 22 N 1.31613 * 119.99867 * 359.97438 * 21 19 18 23 23 Si 1.86800 * 119.99865 * 179.97438 * 21 19 18 24 24 H 1.48507 * 109.46802 * 89.99833 * 23 21 19 25 25 H 1.48496 * 109.47224 * 209.99725 * 23 21 19 26 26 H 1.48496 * 109.47423 * 330.00214 * 23 21 19 27 27 C 1.47186 * 108.67198 * 335.98457 * 18 17 16 28 28 C 1.47286 * 121.05623 * 180.02562 * 9 7 6 29 29 H 1.08997 * 109.47337 * 59.99842 * 1 2 3 30 30 H 1.08996 * 109.47262 * 179.97438 * 1 2 3 31 31 H 1.09008 * 109.47516 * 299.99673 * 1 2 3 32 32 H 1.09002 * 109.46897 * 60.00233 * 3 2 1 33 33 H 1.09003 * 109.47376 * 300.00385 * 3 2 1 34 34 H 1.08999 * 109.47047 * 59.99765 * 4 3 2 35 35 H 1.09003 * 109.47185 * 300.00297 * 4 3 2 36 36 H 1.08999 * 109.47322 * 59.99815 * 5 4 3 37 37 H 1.09003 * 109.47111 * 300.00225 * 5 4 3 38 38 H 1.09004 * 109.46587 * 60.00448 * 6 5 4 39 39 H 1.08996 * 109.47093 * 300.00438 * 6 5 4 40 40 H 1.09002 * 109.63437 * 248.48524 * 10 9 7 41 41 H 1.09003 * 109.63554 * 8.88510 * 10 9 7 42 42 H 1.09004 * 109.46807 * 60.00174 * 13 11 10 43 43 H 1.08998 * 109.47189 * 179.97438 * 13 11 10 44 44 H 1.08993 * 109.47217 * 300.00192 * 13 11 10 45 45 H 1.09002 * 111.01565 * 325.86791 * 16 15 14 46 46 H 1.09000 * 111.01102 * 202.02995 * 16 15 14 47 47 H 1.08993 * 110.35347 * 155.21204 * 17 16 15 48 48 H 1.09002 * 110.35424 * 277.50776 * 17 16 15 49 49 H 0.97003 * 119.99670 * 180.02562 * 22 21 19 50 50 H 1.09003 * 109.90608 * 120.87716 * 27 18 17 51 51 H 1.09002 * 109.90468 * 241.95008 * 27 18 17 52 52 H 1.09004 * 109.64465 * 109.65341 * 28 9 7 53 53 H 1.09003 * 109.80904 * 349.13577 * 28 9 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 8 1.4290 0.0000 0.0000 3 6 2.0102 1.3055 0.0000 4 6 3.5353 1.1848 -0.0006 5 6 4.1576 2.5824 0.0000 6 6 5.6829 2.4617 0.0000 7 6 6.2957 3.8384 0.0013 8 8 5.5829 4.8196 0.0016 9 7 7.6362 3.9794 0.0010 10 6 8.5230 2.8050 -0.0012 11 6 9.5423 2.9674 -1.1344 12 1 9.0299 2.9158 -2.0952 13 6 10.5833 1.8494 -1.0498 14 8 10.1909 4.2364 -1.0063 15 6 9.3178 5.3534 -1.0954 16 6 10.1210 6.6708 -1.0094 17 6 10.5749 6.8547 -2.4803 18 7 9.4424 6.3680 -3.2902 19 6 9.1673 6.7359 -4.5803 20 1 8.3095 6.3211 -5.0887 21 6 9.9910 7.6417 -5.2344 22 7 11.0361 8.1475 -4.6146 23 14 9.6160 8.1429 -6.9944 24 1 8.7270 9.3325 -6.9916 25 1 10.8826 8.4732 -7.6956 26 1 8.9401 7.0216 -7.6952 27 6 8.6566 5.4222 -2.4813 28 6 8.2599 5.3136 0.0028 29 1 -0.3634 0.5138 -0.8899 30 1 -0.3633 -1.0276 -0.0005 31 1 -0.3634 0.5138 0.8901 32 1 1.6885 1.8464 0.8900 33 1 1.6885 1.8465 -0.8900 34 1 3.8566 0.6440 -0.8908 35 1 3.8574 0.6437 0.8891 36 1 3.8363 3.1228 0.8904 37 1 3.8356 3.1240 -0.8895 38 1 6.0045 1.9212 -0.8903 39 1 6.0045 1.9203 0.8897 40 1 9.0445 2.7382 0.9536 41 1 7.9342 1.9017 -0.1609 42 1 10.0872 0.8833 -1.1435 43 1 11.3074 1.9653 -1.8562 44 1 11.0964 1.9026 -0.0897 45 1 10.9792 6.5628 -0.3461 46 1 9.4847 7.4963 -0.6906 47 1 10.7659 7.9072 -2.6893 48 1 11.4664 6.2599 -2.6791 49 1 11.6114 8.7807 -5.0716 50 1 7.6310 5.7787 -2.3853 51 1 8.6649 4.4377 -2.9492 52 1 8.7283 5.4979 0.9696 53 1 7.5039 6.0751 -0.1891 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 ZINC001131642637.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:10:11 Heat of formation + Delta-G solvation = -101.130730 kcal Electronic energy + Delta-G solvation = -31142.283866 eV Core-core repulsion = 27014.936250 eV Total energy + Delta-G solvation = -4127.347616 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.93 seconds Orbital eigenvalues (eV) -40.62827 -38.90968 -38.25743 -36.35412 -35.03730 -34.60288 -33.16389 -31.33241 -30.49747 -29.31069 -28.27721 -26.89078 -26.50793 -25.66898 -24.34722 -23.17644 -21.85961 -21.46404 -20.30560 -19.57854 -18.71566 -17.20850 -16.99928 -16.72788 -16.67429 -16.42330 -15.99135 -15.76729 -15.23070 -14.84381 -14.70048 -14.54202 -14.14715 -13.85002 -13.61993 -13.39573 -13.19516 -12.90776 -12.67049 -12.53135 -12.29832 -12.24848 -12.00691 -11.93941 -11.79172 -11.65858 -11.56201 -11.50764 -11.22356 -11.07171 -10.92453 -10.90585 -10.84260 -10.48442 -10.33879 -10.09752 -10.01203 -9.95801 -9.67555 -9.55748 -9.53479 -8.62587 -8.35286 -6.80962 -5.63002 -3.04362 2.61888 2.85877 3.40964 3.41701 3.46344 3.48424 3.83869 4.38180 4.49753 4.57941 4.61283 4.67354 4.85469 4.91101 4.92853 5.00431 5.01962 5.05101 5.09644 5.14576 5.20749 5.27965 5.31339 5.40160 5.41393 5.46024 5.48585 5.53641 5.57657 5.61585 5.62788 5.70463 5.88039 5.93735 5.95601 6.03499 6.07337 6.14076 6.27247 6.32654 6.48071 6.54829 6.63516 6.70851 6.80248 7.35911 7.58049 7.64202 7.65073 7.87364 8.15577 8.37506 9.24138 9.84299 10.49573 11.40996 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.011493 B = 0.002607 C = 0.002313 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2435.753620 B =10736.846116 C =12103.725591 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.032 3.968 2 O -0.395 6.395 3 C 0.060 3.940 4 C -0.118 4.118 5 C -0.102 4.102 6 C -0.138 4.138 7 C 0.515 3.485 8 O -0.539 6.539 9 N -0.614 5.614 10 C 0.079 3.921 11 C 0.067 3.933 12 H 0.068 0.932 13 C -0.141 4.141 14 O -0.353 6.353 15 C 0.095 3.905 16 C -0.130 4.130 17 C 0.157 3.843 18 N -0.600 5.600 19 C -0.221 4.221 20 H 0.075 0.925 21 C -0.011 4.011 22 N -0.921 5.921 23 Si 0.915 3.085 24 H -0.288 1.288 25 H -0.290 1.290 26 H -0.282 1.282 27 C 0.122 3.878 28 C 0.098 3.902 29 H 0.041 0.959 30 H 0.085 0.915 31 H 0.040 0.960 32 H 0.050 0.950 33 H 0.051 0.949 34 H 0.072 0.928 35 H 0.071 0.929 36 H 0.072 0.928 37 H 0.074 0.926 38 H 0.095 0.905 39 H 0.093 0.907 40 H 0.076 0.924 41 H 0.091 0.909 42 H 0.069 0.931 43 H 0.069 0.931 44 H 0.063 0.937 45 H 0.071 0.929 46 H 0.076 0.924 47 H 0.108 0.892 48 H 0.017 0.983 49 H 0.252 0.748 50 H 0.029 0.971 51 H 0.023 0.977 52 H 0.071 0.929 53 H 0.102 0.898 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -11.096 -16.386 13.694 24.065 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.065 4.065 2 O -0.314 6.314 3 C -0.017 4.017 4 C -0.156 4.156 5 C -0.140 4.140 6 C -0.177 4.177 7 C 0.303 3.697 8 O -0.415 6.415 9 N -0.348 5.348 10 C -0.044 4.044 11 C 0.009 3.991 12 H 0.086 0.914 13 C -0.198 4.198 14 O -0.272 6.272 15 C 0.054 3.946 16 C -0.169 4.169 17 C 0.033 3.967 18 N -0.326 5.326 19 C -0.349 4.349 20 H 0.093 0.907 21 C -0.207 4.207 22 N -0.570 5.570 23 Si 0.745 3.255 24 H -0.216 1.216 25 H -0.216 1.216 26 H -0.208 1.208 27 C -0.004 4.004 28 C -0.024 4.024 29 H 0.060 0.940 30 H 0.103 0.897 31 H 0.059 0.941 32 H 0.068 0.932 33 H 0.069 0.931 34 H 0.091 0.909 35 H 0.089 0.911 36 H 0.090 0.910 37 H 0.092 0.908 38 H 0.113 0.887 39 H 0.111 0.889 40 H 0.094 0.906 41 H 0.109 0.891 42 H 0.088 0.912 43 H 0.088 0.912 44 H 0.082 0.918 45 H 0.090 0.910 46 H 0.094 0.906 47 H 0.125 0.875 48 H 0.035 0.965 49 H 0.063 0.937 50 H 0.047 0.953 51 H 0.041 0.959 52 H 0.089 0.911 53 H 0.120 0.880 Dipole moment (debyes) X Y Z Total from point charges -11.097 -15.871 13.546 23.633 hybrid contribution -0.576 0.303 -0.966 1.165 sum -11.673 -15.568 12.580 23.171 Atomic orbital electron populations 1.22735 0.81027 1.02694 1.00053 1.87857 1.19742 1.28338 1.95466 1.22569 0.94197 0.84438 1.00482 1.21526 0.93053 0.97581 1.03412 1.21246 0.92783 0.95958 1.04031 1.21667 0.96865 0.94849 1.04366 1.20672 0.81474 0.91755 0.75774 1.90705 1.60236 1.40882 1.49720 1.48234 1.05480 1.07509 1.73576 1.22517 0.93011 0.88476 1.00400 1.22452 0.93717 0.85673 0.97291 0.91400 1.21691 0.98882 0.95830 1.03417 1.87866 1.33751 1.12398 1.93222 1.22317 0.90169 0.87944 0.94150 1.23011 0.99152 0.92672 1.02082 1.21946 0.88817 1.02615 0.83355 1.44498 1.28257 1.48494 1.11398 1.19055 1.08528 1.20151 0.87204 0.90671 1.24777 0.97355 0.99825 0.98696 1.69791 1.17111 1.34312 1.35788 0.86067 0.77020 0.78417 0.83965 1.21567 1.21625 1.20804 1.21711 0.94674 0.93995 0.90009 1.22229 0.96372 0.84733 0.99085 0.94021 0.89673 0.94091 0.93192 0.93074 0.90916 0.91054 0.90959 0.90763 0.88709 0.88868 0.90573 0.89081 0.91195 0.91160 0.91782 0.91034 0.90593 0.87485 0.96456 0.93722 0.95321 0.95904 0.91053 0.87959 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.03 0.21 11.01 101.05 1.11 1.32 16 2 O -0.39 -3.80 10.70 -35.23 -0.38 -4.18 16 3 C 0.06 0.50 5.94 37.16 0.22 0.72 16 4 C -0.12 -0.94 5.58 -26.74 -0.15 -1.09 16 5 C -0.10 -0.96 4.50 -26.73 -0.12 -1.08 16 6 C -0.14 -1.28 4.07 -27.88 -0.11 -1.39 16 7 C 0.51 7.13 7.71 -10.98 -0.08 7.05 16 8 O -0.54 -9.53 15.47 5.56 0.09 -9.44 16 9 N -0.61 -8.40 2.95 -172.62 -0.51 -8.91 16 10 C 0.08 0.85 5.82 -3.53 -0.02 0.83 16 11 C 0.07 0.90 3.17 -26.58 -0.08 0.82 16 12 H 0.07 1.02 6.25 -51.93 -0.32 0.70 16 13 C -0.14 -1.46 9.57 37.16 0.36 -1.10 16 14 O -0.35 -6.19 9.90 -35.23 -0.35 -6.54 16 15 C 0.09 1.75 0.81 -89.69 -0.07 1.67 16 16 C -0.13 -2.54 6.21 -24.86 -0.15 -2.70 16 17 C 0.16 3.93 5.89 -2.43 -0.01 3.92 16 18 N -0.60 -15.95 3.40 -170.20 -0.58 -16.53 16 19 C -0.22 -6.35 10.28 -15.06 -0.15 -6.51 16 20 H 0.08 2.11 7.83 -52.49 -0.41 1.70 16 21 C -0.01 -0.32 5.49 -17.43 -0.10 -0.42 16 22 N -0.92 -27.23 10.09 55.49 0.56 -26.67 16 23 Si 0.91 22.51 29.55 -169.99 -5.02 17.48 16 24 H -0.29 -7.06 7.11 56.52 0.40 -6.66 16 25 H -0.29 -7.15 7.11 56.52 0.40 -6.75 16 26 H -0.28 -6.80 7.11 56.52 0.40 -6.39 16 27 C 0.12 2.72 6.34 -3.30 -0.02 2.70 16 28 C 0.10 1.56 5.20 -3.85 -0.02 1.54 16 29 H 0.04 0.25 8.14 -51.93 -0.42 -0.18 16 30 H 0.08 0.45 8.14 -51.93 -0.42 0.02 16 31 H 0.04 0.23 8.14 -51.92 -0.42 -0.19 16 32 H 0.05 0.39 8.14 -51.93 -0.42 -0.04 16 33 H 0.05 0.42 8.14 -51.93 -0.42 -0.01 16 34 H 0.07 0.58 8.14 -51.93 -0.42 0.15 16 35 H 0.07 0.53 8.14 -51.93 -0.42 0.11 16 36 H 0.07 0.72 8.10 -51.93 -0.42 0.30 16 37 H 0.07 0.80 8.10 -51.93 -0.42 0.38 16 38 H 0.09 0.76 7.67 -51.93 -0.40 0.37 16 39 H 0.09 0.66 8.00 -51.93 -0.42 0.24 16 40 H 0.08 0.73 8.14 -51.93 -0.42 0.31 16 41 H 0.09 0.70 6.05 -51.93 -0.31 0.38 16 42 H 0.07 0.53 8.14 -51.93 -0.42 0.11 16 43 H 0.07 0.83 8.14 -51.93 -0.42 0.40 16 44 H 0.06 0.61 8.14 -51.93 -0.42 0.19 16 45 H 0.07 1.26 8.14 -51.93 -0.42 0.84 16 46 H 0.08 1.41 8.01 -51.93 -0.42 1.00 16 47 H 0.11 2.96 5.88 -51.93 -0.31 2.66 16 48 H 0.02 0.47 8.14 -51.93 -0.42 0.05 16 49 H 0.25 7.20 8.31 -40.82 -0.34 6.86 16 50 H 0.03 0.67 7.61 -51.93 -0.40 0.27 16 51 H 0.02 0.52 6.05 -51.93 -0.31 0.21 16 52 H 0.07 1.02 8.14 -51.93 -0.42 0.60 16 53 H 0.10 1.80 6.84 -51.93 -0.36 1.45 16 LS Contribution 409.66 15.07 6.17 6.17 Total: -1.00 -34.29 409.66 -8.98 -43.27 By element: Atomic # 1 Polarization: 8.60 SS G_CDS: -9.54 Total: -0.94 kcal Atomic # 6 Polarization: 5.70 SS G_CDS: 0.58 Total: 6.28 kcal Atomic # 7 Polarization: -51.58 SS G_CDS: -0.53 Total: -52.11 kcal Atomic # 8 Polarization: -19.52 SS G_CDS: -0.64 Total: -20.16 kcal Atomic # 14 Polarization: 22.51 SS G_CDS: -5.02 Total: 17.48 kcal Total LS contribution 6.17 Total: 6.17 kcal Total: -34.29 -8.98 -43.27 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. ZINC001131642637.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 -57.859 kcal (2) G-P(sol) polarization free energy of solvation -34.293 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -92.152 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.978 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -43.272 kcal (6) G-S(sol) free energy of system = (1) + (5) -101.131 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.93 seconds