Wall clock time and date at job start Sat Apr 11 2020 06:10:00 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 N 2 2 C 1.31625 * 1 3 3 Si 1.86802 * 120.00039 * 2 1 4 4 H 1.48507 * 109.46723 * 239.99770 * 3 2 1 5 5 H 1.48501 * 109.47091 * 359.97438 * 3 2 1 6 6 H 1.48503 * 109.47395 * 120.00063 * 3 2 1 7 7 C 1.38794 * 119.99838 * 179.97438 * 2 1 3 8 8 H 1.08004 * 120.00112 * 180.02562 * 7 2 1 9 9 S 1.76200 * 120.00152 * 0.02562 * 7 2 1 10 10 C 1.81000 * 99.99943 * 180.02562 * 9 7 2 11 11 C 1.50702 * 109.47145 * 179.97438 * 10 9 7 12 12 O 1.21278 * 119.99963 * 0.02562 * 11 10 9 13 13 N 1.34778 * 120.00089 * 179.97438 * 11 10 9 14 14 C 1.46505 * 119.99661 * 180.02562 * 13 11 10 15 15 C 1.52999 * 109.46520 * 180.02562 * 14 13 11 16 16 H 1.09001 * 112.85560 * 43.79994 * 15 14 13 17 17 C 1.53787 * 113.60945 * 174.99821 * 15 14 13 18 18 C 1.53788 * 87.08033 * 139.97856 * 17 15 14 19 19 H 1.09002 * 113.60895 * 89.19615 * 18 17 15 20 20 N 1.46894 * 113.61389 * 220.01506 * 18 17 15 21 21 C 1.46906 * 111.00093 * 204.99715 * 20 18 17 22 22 C 1.52997 * 109.46945 * 180.02562 * 21 20 18 23 23 C 1.53006 * 109.46646 * 180.02562 * 22 21 20 24 24 C 1.54266 * 110.03373 * 292.16642 * 23 22 21 25 25 C 1.54883 * 104.19342 * 217.10711 * 24 23 22 26 26 C 1.54884 * 102.74976 * 37.94325 * 25 24 23 27 27 C 1.53873 * 110.02977 * 175.00113 * 23 22 21 28 28 C 1.53781 * 113.61298 * 272.49211 * 15 14 13 29 29 H 0.97003 * 119.99987 * 0.02562 * 1 2 3 30 30 H 1.09006 * 109.47624 * 60.00210 * 10 9 7 31 31 H 1.09003 * 109.47514 * 299.99653 * 10 9 7 32 32 H 0.96997 * 120.00047 * 0.02562 * 13 11 10 33 33 H 1.08991 * 109.47406 * 60.00065 * 14 13 11 34 34 H 1.08996 * 109.46815 * 299.99559 * 14 13 11 35 35 H 1.08997 * 113.61988 * 25.43409 * 17 15 14 36 36 H 1.08998 * 113.61417 * 254.52509 * 17 15 14 37 37 H 1.00897 * 111.00130 * 328.95359 * 20 18 17 38 38 H 1.09004 * 109.46292 * 60.00465 * 21 20 18 39 39 H 1.08990 * 109.47009 * 300.00615 * 21 20 18 40 40 H 1.08995 * 109.47505 * 60.00699 * 22 21 20 41 41 H 1.09000 * 109.47175 * 300.00215 * 22 21 20 42 42 H 1.09000 * 110.01709 * 53.61060 * 23 22 21 43 43 H 1.08998 * 110.48638 * 335.78849 * 24 23 22 44 44 H 1.09000 * 110.48060 * 98.42546 * 24 23 22 45 45 H 1.09004 * 110.75629 * 279.62669 * 25 24 23 46 46 H 1.09001 * 110.91126 * 156.33259 * 25 24 23 47 47 H 1.09000 * 110.48596 * 203.37518 * 26 25 24 48 48 H 1.08999 * 110.60178 * 80.79430 * 26 25 24 49 49 H 1.09006 * 110.03296 * 238.57431 * 27 23 22 50 50 H 1.09003 * 110.03149 * 359.94824 * 27 23 22 51 51 H 1.09000 * 113.61530 * 105.47921 * 28 15 14 52 52 H 1.08998 * 113.61439 * 334.56926 * 28 15 14 Note: An asterisk (*) indicates this parameter was or will be optimized. Cartesian coordinates (angstroms) NO. Atom X Y Z 1 7 0.0000 0.0000 0.0000 2 6 1.3163 0.0000 0.0000 3 14 2.2503 1.6177 0.0000 4 1 3.1041 1.6963 -1.2126 5 1 1.2853 2.7465 -0.0006 6 1 3.1041 1.6965 1.2125 7 6 2.0102 -1.2020 0.0005 8 1 3.0902 -1.2021 0.0010 9 16 1.1292 -2.7279 0.0005 10 6 2.5157 -3.8914 0.0021 11 6 1.9873 -5.3027 0.0030 12 8 0.7910 -5.5023 0.0021 13 7 2.8440 -6.3431 0.0037 14 6 2.3302 -7.7151 0.0040 15 6 3.5023 -8.6985 0.0054 16 1 4.2920 -8.4149 0.7012 17 6 3.0717 -10.1697 0.1287 18 6 4.1664 -10.4945 -0.9015 19 1 5.1383 -10.7036 -0.4545 20 7 3.7637 -11.5042 -1.8895 21 6 4.9354 -12.1980 -2.4407 22 6 4.4769 -13.2390 -3.4640 23 6 5.6973 -13.9620 -4.0376 24 6 6.5335 -12.9892 -4.8944 25 6 7.0383 -13.8585 -6.0727 26 6 5.8234 -14.7728 -6.3679 27 6 5.2458 -15.0960 -4.9746 28 6 4.0147 -9.0482 -1.4017 29 1 -0.4850 0.8401 -0.0004 30 1 3.1237 -3.7310 0.8926 31 1 3.1251 -3.7318 -0.8875 32 1 3.8008 -6.1835 0.0040 33 1 1.7220 -7.8755 -0.8861 34 1 1.7213 -7.8748 0.8938 35 1 2.0614 -10.3608 -0.2329 36 1 3.2580 -10.5981 1.1135 37 1 3.1097 -12.1605 -1.4899 38 1 5.4763 -12.6940 -1.6347 39 1 5.5903 -11.4751 -2.9268 40 1 3.9362 -12.7428 -4.2699 41 1 3.8217 -13.9620 -2.9782 42 1 6.3060 -14.3660 -3.2287 43 1 5.9115 -12.1729 -5.2616 44 1 7.3741 -12.6006 -4.3195 45 1 7.9033 -14.4503 -5.7733 46 1 7.2759 -13.2404 -6.9384 47 1 6.1457 -15.6858 -6.8685 48 1 5.0858 -14.2467 -6.9740 49 1 5.6360 -16.0504 -4.6207 50 1 4.1573 -15.1319 -5.0198 51 1 4.9623 -8.5713 -1.6523 52 1 3.2621 -8.9327 -2.1816 RHF calculation, no. of doubly occupied orbitals= 63 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) 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 ZINC001688900135.mol2 52 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:00 Heat of formation + Delta-G solvation = -51.181715 kcal Electronic energy + Delta-G solvation = -25728.971369 eV Core-core repulsion = 22046.542576 eV Total energy + Delta-G solvation = -3682.428793 eV No. of doubly occupied orbitals = 63 Molecular weight (most abundant/longest-lived isotopes) = 340.193 amu Computer time = 0.69 seconds Orbital eigenvalues (eV) -39.95270 -39.04561 -38.62986 -35.58082 -33.94597 -32.40490 -31.52906 -29.73160 -29.15781 -27.38432 -26.18538 -25.80775 -25.38098 -22.64246 -21.79996 -21.02065 -20.48134 -20.00244 -19.49756 -17.80267 -16.67878 -16.41145 -16.21004 -15.76010 -15.65530 -15.50322 -15.10857 -14.78280 -14.72140 -14.52600 -13.88814 -13.83042 -13.76747 -13.28797 -13.10249 -12.68729 -12.57548 -12.55828 -12.47706 -11.91818 -11.91226 -11.66895 -11.55383 -11.45985 -11.35567 -11.09453 -11.02658 -11.01993 -11.01311 -10.98938 -10.78741 -10.74723 -10.66320 -10.41885 -9.79463 -9.32765 -9.03908 -8.82245 -8.76356 -8.20618 -6.46961 -6.19016 -3.89432 3.08272 3.16646 3.25200 3.29346 3.33757 3.70018 3.78330 4.11723 4.18336 4.21984 4.23066 4.32742 4.42881 4.54083 4.59279 4.71807 4.80843 4.91664 4.99945 5.03307 5.10989 5.15824 5.16145 5.18303 5.23278 5.23733 5.26532 5.36188 5.44801 5.46264 5.50582 5.54238 5.61538 5.65690 5.66422 5.67825 5.71676 5.74789 5.84774 5.89998 5.92116 5.93947 5.99866 6.03423 6.16083 6.18071 6.69691 6.76364 7.00102 7.61587 7.78568 8.37462 8.68366 9.26062 10.88385 Molecular weight = 340.19amu Principal moments of inertia in cm(-1) A = 0.019825 B = 0.001378 C = 0.001324 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1412.032379 B =20317.840702 C =21139.703164 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 N -0.868 5.868 2 C 0.074 3.926 3 Si 0.918 3.082 4 H -0.274 1.274 5 H -0.283 1.283 6 H -0.274 1.274 7 C -0.536 4.536 8 H 0.072 0.928 9 S -0.060 6.060 10 C -0.153 4.153 11 C 0.530 3.470 12 O -0.519 6.519 13 N -0.738 5.738 14 C 0.134 3.866 15 C -0.123 4.123 16 H 0.095 0.905 17 C -0.132 4.132 18 C 0.058 3.942 19 H 0.071 0.929 20 N -0.678 5.678 21 C 0.058 3.942 22 C -0.111 4.111 23 C -0.091 4.091 24 C -0.121 4.121 25 C -0.127 4.127 26 C -0.126 4.126 27 C -0.121 4.121 28 C -0.119 4.119 29 H 0.269 0.731 30 H 0.089 0.911 31 H 0.089 0.911 32 H 0.397 0.603 33 H 0.068 0.932 34 H 0.064 0.936 35 H 0.093 0.907 36 H 0.068 0.932 37 H 0.340 0.660 38 H 0.026 0.974 39 H 0.073 0.927 40 H 0.073 0.927 41 H 0.059 0.941 42 H 0.068 0.932 43 H 0.070 0.930 44 H 0.062 0.938 45 H 0.064 0.936 46 H 0.065 0.935 47 H 0.062 0.938 48 H 0.065 0.935 49 H 0.061 0.939 50 H 0.065 0.935 51 H 0.063 0.937 52 H 0.093 0.907 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 15.244 -31.855 -6.533 35.914 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 N -0.503 5.503 2 C -0.135 4.135 3 Si 0.745 3.255 4 H -0.200 1.200 5 H -0.208 1.208 6 H -0.200 1.200 7 C -0.685 4.685 8 H 0.090 0.910 9 S 0.166 5.834 10 C -0.304 4.304 11 C 0.316 3.684 12 O -0.394 6.394 13 N -0.392 5.392 14 C 0.010 3.990 15 C -0.142 4.142 16 H 0.113 0.887 17 C -0.171 4.171 18 C -0.052 4.052 19 H 0.089 0.911 20 N -0.309 5.309 21 C -0.071 4.071 22 C -0.148 4.148 23 C -0.110 4.110 24 C -0.159 4.159 25 C -0.164 4.164 26 C -0.163 4.163 27 C -0.159 4.159 28 C -0.157 4.157 29 H 0.079 0.921 30 H 0.107 0.893 31 H 0.107 0.893 32 H 0.231 0.769 33 H 0.086 0.914 34 H 0.083 0.917 35 H 0.111 0.889 36 H 0.087 0.913 37 H 0.159 0.841 38 H 0.044 0.956 39 H 0.091 0.909 40 H 0.091 0.909 41 H 0.078 0.922 42 H 0.087 0.913 43 H 0.088 0.912 44 H 0.081 0.919 45 H 0.083 0.917 46 H 0.083 0.917 47 H 0.081 0.919 48 H 0.084 0.916 49 H 0.079 0.921 50 H 0.084 0.916 51 H 0.082 0.918 52 H 0.111 0.889 Dipole moment (debyes) X Y Z Total from point charges 13.016 -31.646 -7.049 34.936 hybrid contribution 2.773 0.755 0.811 2.987 sum 15.790 -30.891 -6.238 35.248 Atomic orbital electron populations 1.70009 1.05951 1.27531 1.46820 1.24981 0.95013 0.98983 0.94544 0.86233 0.79242 0.81284 0.78782 1.19953 1.20807 1.19956 1.22836 0.94880 0.91726 1.59014 0.91002 1.85392 1.05014 0.96270 1.96713 1.23533 1.06499 0.95520 1.04853 1.20301 0.86552 0.86315 0.75211 1.90698 1.15404 1.85648 1.47615 1.45994 1.09842 1.05562 1.77791 1.21242 0.95231 0.80291 1.02247 1.22632 0.98325 0.94528 0.98711 0.88722 1.23197 0.99664 0.95361 0.98865 1.23007 0.97426 0.92806 0.91958 0.91104 1.60609 1.15069 1.18119 1.37115 1.22053 0.92580 0.95910 0.96581 1.21703 0.97176 0.97104 0.98858 1.21760 0.96295 0.95417 0.97561 1.22474 0.99191 0.97868 0.96317 1.22501 0.97717 0.97951 0.98275 1.22353 0.97587 0.99434 0.96940 1.22212 1.00778 0.97993 0.94910 1.23161 1.01722 0.94485 0.96367 0.92147 0.89307 0.89265 0.76925 0.91411 0.91733 0.88891 0.91308 0.84106 0.95582 0.90889 0.90853 0.92218 0.91327 0.91158 0.91875 0.91725 0.91651 0.91938 0.91608 0.92054 0.91625 0.91774 0.88887 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 N -0.87 -26.39 13.22 55.49 0.73 -25.65 16 2 C 0.07 2.13 5.49 -17.43 -0.10 2.03 16 3 Si 0.92 21.54 29.55 -169.99 -5.02 16.51 16 4 H -0.27 -6.27 7.11 56.52 0.40 -5.87 16 5 H -0.28 -6.51 7.11 56.52 0.40 -6.11 16 6 H -0.27 -6.27 7.11 56.52 0.40 -5.86 16 7 C -0.54 -15.44 9.50 30.94 0.29 -15.15 16 8 H 0.07 1.95 7.80 -52.48 -0.41 1.54 16 9 S -0.06 -1.67 20.57 -107.50 -2.21 -3.88 16 10 C -0.15 -3.15 6.16 36.01 0.22 -2.93 16 11 C 0.53 10.14 7.85 -10.98 -0.09 10.05 16 12 O -0.52 -12.11 15.86 5.56 0.09 -12.02 16 13 N -0.74 -10.04 5.56 -60.29 -0.33 -10.38 16 14 C 0.13 1.58 5.51 -4.04 -0.02 1.56 16 15 C -0.12 -1.01 3.75 -89.80 -0.34 -1.35 16 16 H 0.09 0.71 8.14 -51.93 -0.42 0.28 16 17 C -0.13 -0.94 7.21 -25.91 -0.19 -1.13 16 18 C 0.06 0.41 4.14 -66.90 -0.28 0.13 16 19 H 0.07 0.49 7.87 -51.93 -0.41 0.08 16 20 N -0.68 -4.79 6.56 -105.86 -0.69 -5.48 16 21 C 0.06 0.36 5.38 -3.82 -0.02 0.34 16 22 C -0.11 -0.63 4.60 -26.73 -0.12 -0.75 16 23 C -0.09 -0.48 2.93 -89.50 -0.26 -0.74 16 24 C -0.12 -0.67 5.53 -25.08 -0.14 -0.81 16 25 C -0.13 -0.66 6.91 -24.74 -0.17 -0.83 16 26 C -0.13 -0.64 6.83 -25.08 -0.17 -0.81 16 27 C -0.12 -0.62 6.46 -25.62 -0.17 -0.79 16 28 C -0.12 -0.99 7.43 -25.92 -0.19 -1.18 16 29 H 0.27 7.60 8.31 -40.82 -0.34 7.26 16 30 H 0.09 1.68 8.14 -51.91 -0.42 1.26 16 31 H 0.09 1.69 8.14 -51.91 -0.42 1.27 16 32 H 0.40 4.23 8.47 -40.82 -0.35 3.88 16 33 H 0.07 0.89 7.84 -51.93 -0.41 0.48 16 34 H 0.06 0.83 8.14 -51.93 -0.42 0.40 16 35 H 0.09 0.70 8.10 -51.93 -0.42 0.28 16 36 H 0.07 0.44 8.14 -51.93 -0.42 0.01 16 37 H 0.34 2.11 8.55 -39.29 -0.34 1.78 16 38 H 0.03 0.16 8.13 -51.93 -0.42 -0.26 16 39 H 0.07 0.51 7.51 -51.93 -0.39 0.12 16 40 H 0.07 0.50 8.00 -51.93 -0.42 0.08 16 41 H 0.06 0.32 8.14 -51.93 -0.42 -0.10 16 42 H 0.07 0.36 8.14 -51.93 -0.42 -0.06 16 43 H 0.07 0.43 7.76 -51.93 -0.40 0.03 16 44 H 0.06 0.34 8.14 -51.93 -0.42 -0.09 16 45 H 0.06 0.32 8.14 -51.93 -0.42 -0.10 16 46 H 0.06 0.33 8.14 -51.93 -0.42 -0.10 16 47 H 0.06 0.29 8.14 -51.93 -0.42 -0.13 16 48 H 0.07 0.35 8.14 -51.93 -0.42 -0.07 16 49 H 0.06 0.29 8.14 -51.93 -0.42 -0.13 16 50 H 0.06 0.34 7.67 -51.93 -0.40 -0.06 16 51 H 0.06 0.51 8.14 -51.93 -0.42 0.09 16 52 H 0.09 0.92 7.90 -51.93 -0.41 0.51 16 LS Contribution 426.24 15.07 6.42 6.42 Total: -1.00 -33.83 426.24 -12.57 -46.40 By element: Atomic # 1 Polarization: 10.24 SS G_CDS: -9.82 Total: 0.42 kcal Atomic # 6 Polarization: -10.60 SS G_CDS: -1.73 Total: -12.34 kcal Atomic # 7 Polarization: -41.21 SS G_CDS: -0.30 Total: -41.51 kcal Atomic # 8 Polarization: -12.11 SS G_CDS: 0.09 Total: -12.02 kcal Atomic # 14 Polarization: 21.54 SS G_CDS: -5.02 Total: 16.51 kcal Atomic # 16 Polarization: -1.67 SS G_CDS: -2.21 Total: -3.88 kcal Total LS contribution 6.42 Total: 6.42 kcal Total: -33.83 -12.57 -46.40 kcal The number of atoms in the molecule is 52 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. ZINC001688900135.mol2 52 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 -4.782 kcal (2) G-P(sol) polarization free energy of solvation -33.825 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -38.607 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -12.574 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -46.400 kcal (6) G-S(sol) free energy of system = (1) + (5) -51.182 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.70 seconds