Wall clock time and date at job start Tue Mar 31 2020 19:49:48 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.31002 * 1 3 3 C 1.50699 * 119.99865 * 2 1 4 4 N 1.46494 * 109.47528 * 124.99662 * 3 2 1 5 5 C 1.46497 * 120.00130 * 270.00191 * 4 3 2 6 6 C 1.50703 * 109.46987 * 269.99459 * 5 4 3 7 7 H 1.08003 * 120.00123 * 0.02562 * 6 5 4 8 8 C 1.37876 * 120.00195 * 179.97438 * 6 5 4 9 9 N 1.31510 * 119.99841 * 0.02562 * 8 6 5 10 10 Si 1.86804 * 120.00030 * 180.02562 * 8 6 5 11 11 H 1.48504 * 109.47058 * 0.02562 * 10 8 6 12 12 H 1.48493 * 109.47269 * 120.00203 * 10 8 6 13 13 H 1.48506 * 109.47023 * 240.00321 * 10 8 6 14 14 C 1.34784 * 120.00111 * 89.99123 * 4 3 2 15 15 O 1.21280 * 120.00110 * 180.02562 * 14 4 3 16 16 C 1.50699 * 119.99691 * 0.02562 * 14 4 3 17 17 C 1.53004 * 109.46636 * 179.97438 * 16 14 4 18 18 H 1.08997 * 109.47451 * 55.00332 * 17 16 14 19 19 C 1.53005 * 109.46705 * 174.99847 * 17 16 14 20 20 N 1.46494 * 109.46982 * 294.99501 * 17 16 14 21 21 C 1.34784 * 119.99635 * 155.00372 * 20 17 16 22 22 O 1.21279 * 119.99925 * 359.97438 * 21 20 17 23 23 C 1.50700 * 119.99961 * 180.02562 * 21 20 17 24 24 C 1.53003 * 109.46986 * 300.00090 * 23 21 20 25 25 C 1.53005 * 109.47219 * 180.02562 * 24 23 21 26 26 C 1.52994 * 109.46906 * 59.99762 * 25 24 23 27 27 C 1.52998 * 109.47366 * 299.99803 * 26 25 24 28 28 C 1.52992 * 109.47257 * 180.02562 * 23 21 20 29 29 H 1.07997 * 120.00294 * 179.97438 * 1 2 3 30 30 H 1.08005 * 119.99925 * 0.02562 * 1 2 3 31 31 H 1.07991 * 119.99877 * 179.97438 * 2 1 3 32 32 H 1.09000 * 109.47187 * 4.99453 * 3 2 1 33 33 H 1.09001 * 109.46976 * 244.99878 * 3 2 1 34 34 H 1.09003 * 109.47336 * 29.99358 * 5 4 3 35 35 H 1.09003 * 109.47462 * 149.99970 * 5 4 3 36 36 H 0.97006 * 119.99718 * 179.97438 * 9 8 6 37 37 H 1.09002 * 109.47418 * 59.99392 * 16 14 4 38 38 H 1.09001 * 109.46692 * 299.99243 * 16 14 4 39 39 H 1.08997 * 109.47041 * 180.02562 * 19 17 16 40 40 H 1.09005 * 109.46855 * 299.99723 * 19 17 16 41 41 H 1.08996 * 109.47503 * 59.99963 * 19 17 16 42 42 H 0.96995 * 120.00292 * 335.27805 * 20 17 16 43 43 H 1.08997 * 109.47392 * 59.99562 * 23 21 20 44 44 H 1.09005 * 109.46932 * 299.99942 * 24 23 21 45 45 H 1.08998 * 109.47453 * 59.99968 * 24 23 21 46 46 H 1.09000 * 109.46608 * 299.99989 * 25 24 23 47 47 H 1.08999 * 109.46783 * 179.97438 * 25 24 23 48 48 H 1.09005 * 109.46734 * 180.02562 * 26 25 24 49 49 H 1.08992 * 109.47456 * 60.00368 * 26 25 24 50 50 H 1.09000 * 109.47603 * 300.00663 * 27 26 25 51 51 H 1.08998 * 109.47444 * 180.02562 * 27 26 25 52 52 H 1.09004 * 109.47330 * 299.99977 * 28 23 21 53 53 H 1.09005 * 109.47124 * 59.99773 * 28 23 21 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.3100 0.0000 0.0000 3 6 2.0635 1.3051 0.0000 4 7 2.9937 1.3320 1.1314 5 6 2.5587 1.8651 2.4248 6 6 2.8441 3.3438 2.4803 7 1 3.3032 3.8373 1.6363 8 6 2.5183 4.0670 3.6080 9 7 1.9588 3.4663 4.6354 10 14 2.8714 5.9001 3.6764 11 1 3.4994 6.3312 2.4016 12 1 1.6008 6.6418 3.8774 13 1 3.7945 6.1846 4.8045 14 6 4.2497 0.8662 0.9825 15 8 5.0196 0.8881 1.9193 16 6 4.6972 0.3182 -0.3481 17 6 6.1537 -0.1395 -0.2472 18 1 6.2550 -0.8511 0.5722 19 6 6.5697 -0.8091 -1.5585 20 7 7.0135 1.0197 0.0034 21 6 8.2036 0.8563 0.6147 22 8 8.5624 -0.2503 0.9576 23 6 9.0880 2.0489 0.8731 24 6 9.4439 2.7174 -0.4564 25 6 10.3423 3.9279 -0.1941 26 6 11.6235 3.4700 0.5056 27 6 11.2676 2.8015 1.8350 28 6 10.3695 1.5915 1.5725 29 1 -0.5400 -0.9353 -0.0004 30 1 -0.5400 0.9354 -0.0004 31 1 1.8500 -0.9352 0.0004 32 1 1.3585 2.1316 0.0895 33 1 2.6213 1.4026 -0.9314 34 1 1.4886 1.6968 2.5458 35 1 3.0988 1.3606 3.2260 36 1 1.7299 3.9751 5.4290 37 1 4.6125 1.0950 -1.1081 38 1 4.0680 -0.5286 -0.6223 39 1 7.6072 -1.1356 -1.4866 40 1 5.9299 -1.6717 -1.7450 41 1 6.4684 -0.0977 -2.3780 42 1 6.7247 1.9053 -0.2669 43 1 8.5613 2.7611 1.5082 44 1 9.9707 2.0053 -1.0916 45 1 8.5311 3.0433 -0.9551 46 1 9.8159 4.6401 0.4413 47 1 10.5962 4.4038 -1.1413 48 1 12.2633 4.3324 0.6928 49 1 12.1499 2.7578 -0.1298 50 1 10.7412 3.5137 2.4704 51 1 12.1804 2.4757 2.3337 52 1 10.1164 1.1149 2.5196 53 1 10.8961 0.8796 0.9368 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 ZINC000292287158.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:49:48 Heat of formation + Delta-G solvation = -85.255017 kcal Electronic energy + Delta-G solvation = -30426.277173 eV Core-core repulsion = 26490.074031 eV Total energy + Delta-G solvation = -3936.203142 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 0.81 seconds Orbital eigenvalues (eV) -40.18865 -39.14624 -37.62176 -35.90738 -35.22738 -33.21568 -32.16837 -31.98190 -30.80270 -30.34372 -28.62196 -26.92800 -25.38923 -24.02146 -22.91436 -22.61221 -21.16674 -20.29066 -19.98519 -19.61824 -18.43381 -17.04509 -16.60811 -16.32125 -16.16126 -15.55314 -15.21948 -14.80360 -14.69916 -14.43844 -14.27218 -14.16446 -14.04822 -13.47736 -13.38556 -12.98852 -12.91838 -12.64449 -12.41288 -12.32569 -12.29450 -12.10040 -11.82653 -11.71898 -11.64150 -11.31335 -11.27172 -11.21862 -11.04392 -10.90836 -10.87295 -10.71014 -10.61889 -10.47984 -10.21821 -10.18647 -9.65143 -9.38710 -9.21929 -8.96503 -8.73167 -8.46232 -7.69803 -6.06637 -4.11578 2.41565 2.96082 3.26393 3.32839 3.35384 3.37811 3.47446 3.94854 4.41359 4.58093 4.60051 4.65240 4.84309 4.94434 5.11202 5.12554 5.17432 5.20764 5.27279 5.27851 5.33656 5.35861 5.38883 5.45034 5.46884 5.52426 5.55113 5.65226 5.66066 5.73094 5.76448 5.90707 5.97838 5.98061 6.00584 6.06315 6.08488 6.14801 6.16547 6.27113 6.29084 6.38392 6.75310 6.76383 6.88016 7.07035 7.43957 7.53212 7.77381 7.89654 8.21148 8.39572 8.53774 8.88000 9.08097 9.53586 11.01486 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.008696 B = 0.003819 C = 0.002891 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 3218.973177 B = 7330.562633 C = 9683.524268 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.182 4.182 2 C -0.169 4.169 3 C 0.137 3.863 4 N -0.595 5.595 5 C 0.248 3.752 6 C -0.525 4.525 7 H 0.057 0.943 8 C 0.061 3.939 9 N -0.893 5.893 10 Si 0.896 3.104 11 H -0.274 1.274 12 H -0.284 1.284 13 H -0.283 1.283 14 C 0.501 3.499 15 O -0.529 6.529 16 C -0.153 4.153 17 C 0.158 3.842 18 H 0.094 0.906 19 C -0.160 4.160 20 N -0.711 5.711 21 C 0.528 3.472 22 O -0.547 6.547 23 C -0.108 4.108 24 C -0.108 4.108 25 C -0.119 4.119 26 C -0.120 4.120 27 C -0.121 4.121 28 C -0.097 4.097 29 H 0.094 0.906 30 H 0.099 0.901 31 H 0.104 0.896 32 H 0.086 0.914 33 H 0.073 0.927 34 H 0.031 0.969 35 H 0.042 0.958 36 H 0.262 0.738 37 H 0.091 0.909 38 H 0.100 0.900 39 H 0.061 0.939 40 H 0.064 0.936 41 H 0.055 0.945 42 H 0.397 0.603 43 H 0.099 0.901 44 H 0.065 0.935 45 H 0.066 0.934 46 H 0.066 0.934 47 H 0.063 0.937 48 H 0.060 0.940 49 H 0.061 0.939 50 H 0.063 0.937 51 H 0.063 0.937 52 H 0.069 0.931 53 H 0.067 0.933 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 9.931 -3.106 -15.265 18.474 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.220 4.220 2 C -0.188 4.188 3 C 0.013 3.987 4 N -0.326 5.326 5 C 0.129 3.871 6 C -0.564 4.564 7 H 0.075 0.925 8 C -0.141 4.141 9 N -0.532 5.532 10 Si 0.724 3.276 11 H -0.198 1.198 12 H -0.210 1.210 13 H -0.209 1.209 14 C 0.287 3.713 15 O -0.405 6.405 16 C -0.193 4.193 17 C 0.053 3.947 18 H 0.112 0.888 19 C -0.219 4.219 20 N -0.361 5.361 21 C 0.314 3.686 22 O -0.424 6.424 23 C -0.130 4.130 24 C -0.145 4.145 25 C -0.157 4.157 26 C -0.157 4.157 27 C -0.159 4.159 28 C -0.135 4.135 29 H 0.112 0.888 30 H 0.117 0.883 31 H 0.122 0.878 32 H 0.104 0.896 33 H 0.091 0.909 34 H 0.049 0.951 35 H 0.060 0.940 36 H 0.072 0.928 37 H 0.109 0.891 38 H 0.118 0.882 39 H 0.080 0.920 40 H 0.083 0.917 41 H 0.074 0.926 42 H 0.231 0.769 43 H 0.117 0.883 44 H 0.084 0.916 45 H 0.084 0.916 46 H 0.084 0.916 47 H 0.082 0.918 48 H 0.079 0.921 49 H 0.079 0.921 50 H 0.082 0.918 51 H 0.082 0.918 52 H 0.088 0.912 53 H 0.085 0.915 Dipole moment (debyes) X Y Z Total from point charges 10.518 -3.850 -14.307 18.169 hybrid contribution -0.380 1.471 0.022 1.520 sum 10.138 -2.378 -14.284 17.677 Atomic orbital electron populations 1.23821 0.95879 1.01755 1.00565 1.23433 0.96409 0.99062 0.99919 1.20567 0.93148 0.95272 0.89761 1.48242 1.10520 1.60379 1.13499 1.19250 0.95306 0.93720 0.78869 1.21442 1.40453 0.93304 1.01163 0.92509 1.24932 0.95099 0.98660 0.95383 1.69946 1.42056 1.36040 1.05175 0.86478 0.77938 0.85107 0.78079 1.19849 1.20981 1.20892 1.20283 0.82420 0.78973 0.89650 1.90718 1.52334 1.55727 1.41732 1.22125 0.96821 1.02725 0.97678 1.20626 0.87351 0.88683 0.98030 0.88834 1.22177 1.02302 1.00496 0.96891 1.45838 1.18799 1.10481 1.60946 1.20453 0.81581 0.90412 0.76146 1.90646 1.73718 1.26741 1.51346 1.21354 0.96623 0.96231 0.98765 1.21504 0.99969 0.97445 0.95631 1.21536 0.96182 0.97369 1.00566 1.21550 0.97252 1.00097 0.96848 1.21566 0.99530 0.97557 0.97221 1.21299 0.93917 0.97190 1.01069 0.88801 0.88309 0.87804 0.89573 0.90850 0.95144 0.93986 0.92778 0.89108 0.88225 0.91993 0.91722 0.92578 0.76943 0.88293 0.91602 0.91573 0.91570 0.91840 0.92145 0.92070 0.91785 0.91843 0.91177 0.91475 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 12. 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.18 -2.27 14.32 29.01 0.42 -1.86 16 2 C -0.17 -2.28 9.74 -36.03 -0.35 -2.63 16 3 C 0.14 2.08 5.65 -5.20 -0.03 2.05 16 4 N -0.59 -11.94 2.44 -176.44 -0.43 -12.37 16 5 C 0.25 6.36 5.16 -5.19 -0.03 6.33 16 6 C -0.52 -14.49 9.39 -34.44 -0.32 -14.81 16 7 H 0.06 1.55 7.81 -52.48 -0.41 1.14 16 8 C 0.06 1.73 5.46 -17.82 -0.10 1.63 16 9 N -0.89 -26.41 13.29 55.43 0.74 -25.67 16 10 Si 0.90 21.48 29.54 -169.99 -5.02 16.46 16 11 H -0.27 -6.40 7.11 56.52 0.40 -6.00 16 12 H -0.28 -6.66 7.11 56.52 0.40 -6.25 16 13 H -0.28 -6.82 7.11 56.52 0.40 -6.42 16 14 C 0.50 9.72 7.53 -10.98 -0.08 9.63 16 15 O -0.53 -12.77 14.93 5.56 0.08 -12.69 16 16 C -0.15 -2.02 3.79 -27.88 -0.11 -2.13 16 17 C 0.16 2.21 2.14 -67.93 -0.15 2.06 16 18 H 0.09 1.63 7.39 -51.93 -0.38 1.25 16 19 C -0.16 -1.73 9.48 37.16 0.35 -1.38 16 20 N -0.71 -10.17 4.11 -53.80 -0.22 -10.39 16 21 C 0.53 8.34 7.12 -10.98 -0.08 8.26 16 22 O -0.55 -10.47 15.64 5.56 0.09 -10.39 16 23 C -0.11 -1.39 2.53 -91.77 -0.23 -1.62 16 24 C -0.11 -1.07 5.10 -26.73 -0.14 -1.21 16 25 C -0.12 -1.06 5.92 -26.73 -0.16 -1.22 16 26 C -0.12 -1.05 5.92 -26.74 -0.16 -1.21 16 27 C -0.12 -1.22 5.92 -26.75 -0.16 -1.38 16 28 C -0.10 -1.20 4.49 -26.76 -0.12 -1.32 16 29 H 0.09 1.01 8.06 -52.49 -0.42 0.59 16 30 H 0.10 1.22 7.84 -52.48 -0.41 0.81 16 31 H 0.10 1.31 7.83 -52.49 -0.41 0.90 16 32 H 0.09 1.39 7.72 -51.93 -0.40 0.99 16 33 H 0.07 0.89 6.44 -51.93 -0.33 0.55 16 34 H 0.03 0.81 8.07 -51.93 -0.42 0.39 16 35 H 0.04 1.24 7.42 -51.93 -0.39 0.86 16 36 H 0.26 7.33 8.31 -40.82 -0.34 6.99 16 37 H 0.09 0.99 7.40 -51.93 -0.38 0.61 16 38 H 0.10 1.10 7.78 -51.93 -0.40 0.70 16 39 H 0.06 0.75 8.14 -51.93 -0.42 0.33 16 40 H 0.06 0.60 8.14 -51.93 -0.42 0.17 16 41 H 0.06 0.51 8.14 -51.93 -0.42 0.09 16 42 H 0.40 5.08 8.03 -40.82 -0.33 4.75 16 43 H 0.10 1.34 8.14 -51.93 -0.42 0.92 16 44 H 0.07 0.65 8.14 -51.93 -0.42 0.23 16 45 H 0.07 0.61 7.57 -51.93 -0.39 0.21 16 46 H 0.07 0.61 8.14 -51.93 -0.42 0.19 16 47 H 0.06 0.49 8.14 -51.93 -0.42 0.07 16 48 H 0.06 0.47 8.14 -51.93 -0.42 0.05 16 49 H 0.06 0.52 8.14 -51.93 -0.42 0.10 16 50 H 0.06 0.66 8.14 -51.93 -0.42 0.24 16 51 H 0.06 0.58 8.14 -51.93 -0.42 0.16 16 52 H 0.07 0.99 8.10 -51.93 -0.42 0.57 16 53 H 0.07 0.85 8.10 -51.93 -0.42 0.43 16 LS Contribution 424.37 15.07 6.40 6.40 Total: -1.00 -34.35 424.37 -9.52 -43.87 By element: Atomic # 1 Polarization: 15.30 SS G_CDS: -9.71 Total: 5.59 kcal Atomic # 6 Polarization: 0.63 SS G_CDS: -1.44 Total: -0.81 kcal Atomic # 7 Polarization: -48.52 SS G_CDS: 0.08 Total: -48.43 kcal Atomic # 8 Polarization: -23.25 SS G_CDS: 0.17 Total: -23.08 kcal Atomic # 14 Polarization: 21.48 SS G_CDS: -5.02 Total: 16.46 kcal Total LS contribution 6.40 Total: 6.40 kcal Total: -34.35 -9.52 -43.87 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. ZINC000292287158.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 -41.387 kcal (2) G-P(sol) polarization free energy of solvation -34.348 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -75.734 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.521 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -43.868 kcal (6) G-S(sol) free energy of system = (1) + (5) -85.255 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.81 seconds