Wall clock time and date at job start Tue Mar 31 2020 19:47:17 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.42895 * 1 3 3 C 1.42903 * 113.99424 * 2 1 4 4 H 1.08993 * 109.47634 * 338.26250 * 3 2 1 5 5 C 1.50702 * 109.46994 * 218.25854 * 3 2 1 6 6 O 1.21282 * 120.00049 * 355.00444 * 5 3 2 7 7 N 1.34774 * 120.00392 * 175.00368 * 5 3 2 8 8 C 1.47176 * 125.65634 * 0.02562 * 7 5 3 9 9 C 1.53698 * 107.10907 * 181.47240 * 8 7 5 10 10 C 1.54517 * 103.64807 * 21.54082 * 9 8 7 11 11 C 1.47482 * 125.66471 * 179.72042 * 7 5 3 12 12 H 1.08999 * 110.35365 * 274.87475 * 11 7 5 13 13 C 1.52992 * 110.35497 * 37.16068 * 11 7 5 14 14 C 1.52275 * 110.63345 * 263.43223 * 9 8 7 15 15 C 1.53320 * 109.41580 * 296.23300 * 14 9 8 16 16 N 1.46918 * 108.71080 * 306.49289 * 15 14 9 17 17 C 1.36936 * 120.56403 * 234.25600 * 16 15 14 18 18 H 1.07999 * 120.00236 * 30.01712 * 17 16 15 19 19 C 1.38814 * 120.00074 * 210.01469 * 17 16 15 20 20 N 1.31615 * 119.99848 * 13.10781 * 19 17 16 21 21 Si 1.86803 * 119.99990 * 193.10250 * 19 17 16 22 22 H 1.48501 * 109.47018 * 90.00359 * 21 19 17 23 23 H 1.48503 * 109.47106 * 210.00219 * 21 19 17 24 24 H 1.48496 * 109.47178 * 330.00266 * 21 19 17 25 25 C 1.46935 * 118.87569 * 54.23200 * 16 15 14 26 26 C 1.52273 * 110.62432 * 139.65864 * 9 8 7 27 27 C 1.52998 * 109.46831 * 98.26472 * 3 2 1 28 28 C 1.52999 * 117.50141 * 351.10790 * 27 3 2 29 29 C 1.53001 * 117.50071 * 59.73293 * 27 3 2 30 30 H 1.09001 * 109.47335 * 305.44480 * 1 2 3 31 31 H 1.09002 * 109.47662 * 65.44993 * 1 2 3 32 32 H 1.09006 * 109.46982 * 185.44751 * 1 2 3 33 33 H 1.09001 * 109.90674 * 300.83664 * 8 7 5 34 34 H 1.08999 * 109.90859 * 61.91360 * 8 7 5 35 35 H 1.08998 * 111.01084 * 83.24916 * 10 9 8 36 36 H 1.08998 * 111.01452 * 207.11884 * 10 9 8 37 37 H 1.09006 * 109.47208 * 58.55030 * 13 11 7 38 38 H 1.09004 * 109.47165 * 178.54575 * 13 11 7 39 39 H 1.09004 * 109.47372 * 298.54831 * 13 11 7 40 40 H 1.09004 * 109.47939 * 56.21573 * 14 9 8 41 41 H 1.08998 * 109.48454 * 176.24937 * 14 9 8 42 42 H 1.08997 * 109.59570 * 186.72313 * 15 14 9 43 43 H 1.09002 * 109.72323 * 66.33710 * 15 14 9 44 44 H 0.97005 * 120.00321 * 179.97438 * 20 19 17 45 45 H 1.08998 * 109.60198 * 185.99776 * 25 16 15 46 46 H 1.08996 * 109.72775 * 65.59543 * 25 16 15 47 47 H 1.09002 * 109.47802 * 183.73877 * 26 9 8 48 48 H 1.08991 * 109.48011 * 303.75822 * 26 9 8 49 49 H 1.09004 * 115.55831 * 205.42246 * 27 3 2 50 50 H 1.08998 * 117.49779 * 214.98251 * 28 27 3 51 51 H 1.09002 * 117.49556 * 0.02562 * 28 27 3 52 52 H 1.08993 * 117.50127 * 359.97438 * 29 27 3 53 53 H 1.09000 * 117.49415 * 144.99734 * 29 27 3 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.4289 0.0000 0.0000 3 6 2.0101 1.3055 0.0000 4 1 1.2858 2.0257 -0.3806 5 6 3.2336 1.3108 -0.8798 6 8 3.6267 0.2772 -1.3780 7 7 3.8910 2.4639 -1.1131 8 6 3.5255 3.7839 -0.5745 9 6 4.5347 4.7972 -1.1376 10 6 5.0906 4.0965 -2.3976 11 6 5.0997 2.6115 -1.9453 12 1 5.9926 2.3970 -1.3582 13 6 5.0253 1.6851 -3.1606 14 6 5.6558 5.0498 -0.1386 15 6 5.0761 5.6948 1.1257 16 7 4.3068 6.8847 0.7374 17 6 4.5765 8.1030 1.3016 18 1 5.5787 8.3361 1.6298 19 6 3.5619 9.0384 1.4515 20 7 2.3057 8.6770 1.2975 21 14 3.9794 10.8099 1.8721 22 1 4.1696 11.5844 0.6194 23 1 2.8704 11.4063 2.6593 24 1 5.2300 10.8472 2.6719 25 6 3.2433 6.7455 -0.2669 26 6 3.8350 6.0939 -1.5222 27 6 2.4068 1.6871 1.4275 28 6 1.9451 0.7667 2.5591 29 6 3.4175 0.7890 2.1438 30 1 -0.3634 0.5960 0.8372 31 1 -0.3634 0.4270 -0.9347 32 1 -0.3633 -1.0231 0.0976 33 1 2.5171 4.0484 -0.8929 34 1 3.5789 3.7684 0.5140 35 1 4.4290 4.2446 -3.2511 36 1 6.0999 4.4407 -2.6231 37 1 4.1355 1.9219 -3.7440 38 1 5.9125 1.8244 -3.7784 39 1 4.9750 0.6492 -2.8254 40 1 6.1311 4.1038 0.1211 41 1 6.3941 5.7180 -0.5818 42 1 5.8876 5.9857 1.7927 43 1 4.4219 4.9860 1.6335 44 1 1.5966 9.3306 1.4026 45 1 2.8465 7.7291 -0.5180 46 1 2.4440 6.1186 0.1281 47 1 4.5544 6.7732 -1.9797 48 1 3.0360 5.8791 -2.2316 49 1 2.4488 2.7566 1.6341 50 1 1.6832 1.2306 3.5100 51 1 1.3589 -0.1110 2.2865 52 1 3.7994 -0.0738 1.5983 53 1 4.1245 1.2678 2.8213 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 ZINC001086920255.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 19:47:17 Heat of formation + Delta-G solvation = -26.905565 kcal Electronic energy + Delta-G solvation = -32013.945849 eV Core-core repulsion = 28080.272930 eV Total energy + Delta-G solvation = -3933.672920 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 1.42 seconds Orbital eigenvalues (eV) -40.77708 -39.25197 -37.74582 -36.02792 -35.81816 -33.12224 -32.07104 -30.93450 -29.94459 -28.75547 -27.23584 -26.20200 -25.64635 -24.27403 -22.75002 -22.39469 -21.94749 -20.98191 -20.63449 -19.06210 -18.26293 -17.08007 -16.81151 -16.70753 -16.21046 -15.85146 -15.08719 -14.86259 -14.74340 -14.30335 -14.20809 -13.89571 -13.75427 -13.35146 -13.30476 -13.02684 -12.99486 -12.64849 -12.56976 -12.11229 -12.09336 -11.73181 -11.66812 -11.54322 -11.46429 -11.39985 -11.24311 -11.13922 -10.89372 -10.78616 -10.55830 -10.49898 -10.38113 -10.24205 -10.20879 -9.93860 -9.86873 -9.57437 -9.38679 -9.13486 -8.51708 -8.46456 -6.67893 -5.83198 -3.29064 2.67027 2.80835 3.34012 3.41415 3.42004 3.64277 3.87775 4.44941 4.48756 4.58435 4.70114 4.77979 4.89570 5.00516 5.06276 5.12268 5.22738 5.25642 5.30350 5.32585 5.34457 5.36520 5.44246 5.50822 5.58084 5.65324 5.67283 5.75208 5.80729 5.85003 5.93041 5.95791 6.03265 6.06626 6.12824 6.18650 6.21508 6.26034 6.27541 6.47495 6.55642 6.57609 6.62562 6.69400 6.79615 6.91004 6.92441 7.30623 7.59800 7.67232 7.75250 8.30144 8.42403 9.17793 9.79799 10.36652 11.24663 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.013521 B = 0.003709 C = 0.003236 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 2070.319183 B = 7548.378404 C = 8651.361108 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.027 3.973 2 O -0.357 6.357 3 C 0.106 3.894 4 H 0.089 0.911 5 C 0.536 3.464 6 O -0.528 6.528 7 N -0.608 5.608 8 C 0.106 3.894 9 C -0.054 4.054 10 C -0.124 4.124 11 C 0.138 3.862 12 H 0.074 0.926 13 C -0.164 4.164 14 C -0.120 4.120 15 C 0.145 3.855 16 N -0.590 5.590 17 C -0.252 4.252 18 H 0.064 0.936 19 C 0.005 3.995 20 N -0.905 5.905 21 Si 0.902 3.098 22 H -0.287 1.287 23 H -0.291 1.291 24 H -0.279 1.279 25 C 0.170 3.830 26 C -0.128 4.128 27 C -0.184 4.184 28 C -0.163 4.163 29 C -0.165 4.165 30 H 0.047 0.953 31 H 0.042 0.958 32 H 0.088 0.912 33 H 0.066 0.934 34 H 0.080 0.920 35 H 0.073 0.927 36 H 0.080 0.920 37 H 0.052 0.948 38 H 0.058 0.942 39 H 0.084 0.916 40 H 0.051 0.949 41 H 0.063 0.937 42 H 0.053 0.947 43 H 0.017 0.983 44 H 0.254 0.746 45 H 0.089 0.911 46 H 0.017 0.983 47 H 0.060 0.940 48 H 0.050 0.950 49 H 0.122 0.878 50 H 0.095 0.905 51 H 0.101 0.899 52 H 0.097 0.903 53 H 0.095 0.905 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -0.505 -13.611 -4.229 14.262 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.070 4.070 2 O -0.274 6.274 3 C 0.044 3.956 4 H 0.107 0.893 5 C 0.324 3.676 6 O -0.405 6.405 7 N -0.343 5.343 8 C -0.017 4.017 9 C -0.054 4.054 10 C -0.162 4.162 11 C 0.035 3.965 12 H 0.092 0.908 13 C -0.222 4.222 14 C -0.159 4.159 15 C 0.019 3.981 16 N -0.313 5.313 17 C -0.382 4.382 18 H 0.082 0.918 19 C -0.194 4.194 20 N -0.549 5.549 21 Si 0.733 3.267 22 H -0.214 1.214 23 H -0.217 1.217 24 H -0.204 1.204 25 C 0.045 3.955 26 C -0.167 4.167 27 C -0.202 4.202 28 C -0.200 4.200 29 C -0.202 4.202 30 H 0.066 0.934 31 H 0.061 0.939 32 H 0.107 0.893 33 H 0.084 0.916 34 H 0.098 0.902 35 H 0.092 0.908 36 H 0.099 0.901 37 H 0.071 0.929 38 H 0.077 0.923 39 H 0.103 0.897 40 H 0.070 0.930 41 H 0.081 0.919 42 H 0.072 0.928 43 H 0.035 0.965 44 H 0.064 0.936 45 H 0.107 0.893 46 H 0.035 0.965 47 H 0.079 0.921 48 H 0.068 0.932 49 H 0.140 0.860 50 H 0.113 0.887 51 H 0.119 0.881 52 H 0.115 0.885 53 H 0.114 0.886 Dipole moment (debyes) X Y Z Total from point charges -0.389 -14.309 -4.434 14.985 hybrid contribution -0.155 1.388 0.241 1.417 sum -0.544 -12.921 -4.194 13.595 Atomic orbital electron populations 1.22914 0.80737 1.02978 1.00355 1.87985 1.19595 1.25213 1.94616 1.21334 0.94157 0.86509 0.93560 0.89310 1.19893 0.84347 0.82417 0.80962 1.90623 1.67372 1.31455 1.51033 1.48567 1.28679 1.07559 1.49506 1.22370 0.98808 0.81005 0.99480 1.20927 0.95328 0.94743 0.94428 1.23021 1.02036 0.93094 0.98012 1.21955 0.87352 0.95597 0.91565 0.90833 1.22176 1.01303 1.01948 0.96735 1.21929 0.96661 0.99610 0.97682 1.20949 0.94673 0.89446 0.93067 1.43995 1.39340 1.03716 1.44216 1.19195 0.93161 0.90642 1.35176 0.91815 1.25084 0.95073 0.99142 1.00087 1.69976 0.98141 1.33107 1.53697 0.86289 0.77394 0.85063 0.78001 1.21449 1.21741 1.20449 1.20826 0.89087 1.00317 0.85306 1.22137 0.99191 0.94893 1.00504 1.22746 1.07187 0.98535 0.91751 1.22920 0.97924 1.01590 0.97567 1.22887 0.92435 1.01814 1.03034 0.93402 0.93924 0.89316 0.91553 0.90191 0.90793 0.90115 0.92921 0.92336 0.89706 0.93044 0.91853 0.92806 0.96531 0.93634 0.89263 0.96482 0.92087 0.93159 0.85974 0.88695 0.88056 0.88501 0.88650 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.03 0.26 10.11 101.05 1.02 1.28 16 2 O -0.36 -4.80 8.48 -55.09 -0.47 -5.27 16 3 C 0.11 1.40 2.43 -27.88 -0.07 1.33 16 4 H 0.09 1.07 7.72 -51.93 -0.40 0.67 16 5 C 0.54 8.22 6.68 -10.99 -0.07 8.15 16 6 O -0.53 -9.56 14.16 -14.30 -0.20 -9.77 16 7 N -0.61 -8.63 3.25 -163.07 -0.53 -9.16 16 8 C 0.11 1.52 4.25 -3.31 -0.01 1.51 16 9 C -0.05 -0.79 0.53 -154.09 -0.08 -0.88 16 10 C -0.12 -1.57 5.42 -24.75 -0.13 -1.70 16 11 C 0.14 1.82 3.45 -66.21 -0.23 1.60 16 12 H 0.07 0.98 7.74 -51.93 -0.40 0.58 16 13 C -0.16 -2.10 8.71 37.15 0.32 -1.78 16 14 C -0.12 -1.89 4.78 -26.94 -0.13 -2.01 16 15 C 0.15 2.79 6.22 -3.65 -0.02 2.77 16 16 N -0.59 -13.66 2.98 -172.17 -0.51 -14.17 16 17 C -0.25 -6.57 9.58 -15.06 -0.14 -6.72 16 18 H 0.06 1.65 7.86 -52.49 -0.41 1.24 16 19 C 0.00 0.12 5.00 -17.43 -0.09 0.04 16 20 N -0.90 -25.01 11.08 55.49 0.61 -24.40 16 21 Si 0.90 20.98 29.57 -169.99 -5.03 15.96 16 22 H -0.29 -6.68 7.11 56.52 0.40 -6.28 16 23 H -0.29 -6.87 7.11 56.52 0.40 -6.46 16 24 H -0.28 -6.39 7.11 56.52 0.40 -5.98 16 25 C 0.17 3.70 5.01 -3.64 -0.02 3.68 16 26 C -0.13 -2.18 4.70 -26.94 -0.13 -2.31 16 27 C -0.18 -2.29 4.94 -90.62 -0.45 -2.74 16 28 C -0.16 -1.82 9.73 -26.73 -0.26 -2.08 16 29 C -0.16 -2.04 9.81 -26.73 -0.26 -2.30 16 30 H 0.05 0.41 8.01 -51.93 -0.42 -0.01 16 31 H 0.04 0.38 8.14 -51.93 -0.42 -0.04 16 32 H 0.09 0.70 8.14 -51.93 -0.42 0.28 16 33 H 0.07 0.96 7.59 -51.93 -0.39 0.57 16 34 H 0.08 1.22 4.09 -51.93 -0.21 1.00 16 35 H 0.07 0.92 7.93 -51.93 -0.41 0.51 16 36 H 0.08 0.94 8.01 -51.93 -0.42 0.53 16 37 H 0.05 0.69 8.14 -51.93 -0.42 0.26 16 38 H 0.06 0.63 8.14 -51.93 -0.42 0.21 16 39 H 0.08 1.27 5.89 -51.93 -0.31 0.96 16 40 H 0.05 0.74 7.47 -51.93 -0.39 0.35 16 41 H 0.06 1.03 8.14 -51.93 -0.42 0.61 16 42 H 0.05 1.04 7.93 -51.93 -0.41 0.63 16 43 H 0.02 0.32 6.54 -51.93 -0.34 -0.02 16 44 H 0.25 6.83 8.31 -40.82 -0.34 6.49 16 45 H 0.09 2.24 6.04 -51.93 -0.31 1.93 16 46 H 0.02 0.38 7.47 -51.93 -0.39 0.00 16 47 H 0.06 1.08 8.14 -51.93 -0.42 0.66 16 48 H 0.05 0.78 8.10 -51.93 -0.42 0.36 16 49 H 0.12 1.57 6.53 -51.93 -0.34 1.24 16 50 H 0.09 0.96 8.14 -51.93 -0.42 0.53 16 51 H 0.10 1.11 6.54 -51.93 -0.34 0.77 16 52 H 0.10 1.34 8.14 -51.93 -0.42 0.91 16 53 H 0.10 1.11 8.14 -51.93 -0.42 0.68 16 LS Contribution 395.25 15.07 5.96 5.96 Total: -1.00 -29.67 395.25 -10.17 -39.84 By element: Atomic # 1 Polarization: 12.42 SS G_CDS: -9.25 Total: 3.16 kcal Atomic # 6 Polarization: -1.40 SS G_CDS: -0.75 Total: -2.15 kcal Atomic # 7 Polarization: -47.30 SS G_CDS: -0.43 Total: -47.73 kcal Atomic # 8 Polarization: -14.36 SS G_CDS: -0.67 Total: -15.03 kcal Atomic # 14 Polarization: 20.98 SS G_CDS: -5.03 Total: 15.96 kcal Total LS contribution 5.96 Total: 5.96 kcal Total: -29.67 -10.17 -39.84 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. ZINC001086920255.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 12.932 kcal (2) G-P(sol) polarization free energy of solvation -29.666 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -16.734 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.171 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -39.837 kcal (6) G-S(sol) free energy of system = (1) + (5) -26.906 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.42 seconds