Wall clock time and date at job start Sat Apr 11 2020 06:08:08 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.31000 * 1 3 3 C 1.50698 * 119.99682 * 2 1 4 4 C 1.50696 * 120.00359 * 180.02562 * 2 1 3 5 5 C 1.52995 * 109.47395 * 95.00051 * 4 2 1 6 6 N 1.46905 * 109.47174 * 180.02562 * 5 4 2 7 7 C 1.46906 * 110.99713 * 179.97438 * 6 5 4 8 8 H 1.09003 * 109.45785 * 35.02932 * 7 6 5 9 9 C 1.53038 * 109.46152 * 275.06677 * 7 6 5 10 10 C 1.53624 * 109.24306 * 178.49453 * 9 7 6 11 11 N 1.46923 * 108.54329 * 54.85499 * 10 9 7 12 12 C 1.34778 * 120.63340 * 126.36848 * 11 10 9 13 13 O 1.21278 * 120.00140 * 179.71626 * 12 11 10 14 14 C 1.50703 * 119.99723 * 359.72324 * 12 11 10 15 15 S 1.80993 * 109.47026 * 179.97438 * 14 12 11 16 16 C 1.76201 * 99.99993 * 179.97438 * 15 14 12 17 17 H 1.08002 * 120.00153 * 0.02562 * 16 15 14 18 18 C 1.38798 * 120.00523 * 179.97438 * 16 15 14 19 19 N 1.31634 * 119.99397 * 0.02562 * 18 16 15 20 20 Si 1.86798 * 120.00627 * 179.97438 * 18 16 15 21 21 H 1.48490 * 109.47339 * 180.02562 * 20 18 16 22 22 H 1.48493 * 109.47239 * 300.00544 * 20 18 16 23 23 H 1.48504 * 109.46636 * 60.00103 * 20 18 16 24 24 C 1.46932 * 118.73776 * 306.09417 * 11 10 9 25 25 C 1.53047 * 109.45599 * 154.99849 * 7 6 5 26 26 H 1.09003 * 109.52474 * 301.35549 * 25 7 6 27 27 C 1.53002 * 109.49440 * 61.43223 * 25 7 6 28 28 C 1.53003 * 109.46845 * 184.83406 * 27 25 7 29 29 H 1.07993 * 120.00567 * 0.02562 * 1 2 3 30 30 H 1.08005 * 119.99438 * 179.97438 * 1 2 3 31 31 H 1.09000 * 109.47230 * 240.00400 * 3 2 1 32 32 H 1.09003 * 109.47132 * 359.97438 * 3 2 1 33 33 H 1.08997 * 109.47322 * 120.00123 * 3 2 1 34 34 H 1.09002 * 109.47343 * 334.99489 * 4 2 1 35 35 H 1.09008 * 109.46958 * 214.99792 * 4 2 1 36 36 H 1.09003 * 109.47559 * 59.99936 * 5 4 2 37 37 H 1.09002 * 109.47186 * 299.99494 * 5 4 2 38 38 H 1.00893 * 110.99641 * 303.94886 * 6 5 4 39 39 H 1.09004 * 109.49479 * 298.57206 * 9 7 6 40 40 H 1.09002 * 109.49752 * 58.66540 * 9 7 6 41 41 H 1.09005 * 109.73479 * 295.11935 * 10 9 7 42 42 H 1.08990 * 109.66657 * 174.54899 * 10 9 7 43 43 H 1.09001 * 109.47144 * 59.99673 * 14 12 11 44 44 H 1.09006 * 109.47011 * 299.99863 * 14 12 11 45 45 H 0.97005 * 119.99723 * 179.97438 * 19 18 16 46 46 H 1.09008 * 109.58694 * 173.64463 * 24 11 10 47 47 H 1.08998 * 109.59002 * 293.93328 * 24 11 10 48 48 H 1.08990 * 109.47405 * 304.83514 * 27 25 7 49 49 H 1.09000 * 109.46813 * 64.84361 * 27 25 7 50 50 H 1.09002 * 109.47091 * 60.00231 * 28 27 25 51 51 H 1.09005 * 109.47265 * 180.02562 * 28 27 25 52 52 H 1.09002 * 109.47463 * 300.00477 * 28 27 25 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.0634 1.3051 0.0000 4 6 2.0636 -1.3050 -0.0006 5 6 2.4275 -1.6832 -1.4377 6 7 3.1616 -2.9557 -1.4383 7 6 3.5278 -3.3515 -2.8048 8 1 3.7575 -2.4619 -3.3912 9 6 2.3586 -4.0993 -3.4497 10 6 2.7328 -4.4777 -4.8908 11 7 3.9914 -5.2352 -4.8631 12 6 4.0778 -6.4450 -5.4508 13 8 5.1239 -7.0582 -5.4256 14 6 2.8745 -7.0358 -6.1395 15 16 3.3060 -8.6545 -6.8249 16 6 1.7629 -9.1013 -7.5488 17 1 0.9194 -8.4305 -7.4792 18 6 1.6318 -10.3149 -8.2096 19 7 2.6601 -11.1322 -8.2949 20 14 -0.0042 -10.7889 -8.9765 21 1 0.1174 -12.1231 -9.6168 22 1 -1.0489 -10.8396 -7.9223 23 1 -0.3818 -9.7809 -9.9995 24 6 5.1558 -4.6461 -4.1879 25 6 4.7560 -4.2635 -2.7594 26 1 5.5818 -3.7389 -2.2787 27 6 4.4236 -5.5286 -1.9658 28 6 5.6838 -6.3836 -1.8178 29 1 -0.5401 0.9352 -0.0004 30 1 -0.5399 -0.9354 -0.0004 31 1 2.6900 1.3629 -0.8900 32 1 1.3550 2.1336 -0.0005 33 1 2.6901 1.3629 0.8899 34 1 1.4387 -2.0856 0.4335 35 1 2.9743 -1.1989 0.5889 36 1 3.0524 -0.9027 -1.8717 37 1 1.5168 -1.7894 -2.0272 38 1 3.9767 -2.9031 -0.8460 39 1 2.1444 -5.0021 -2.8777 40 1 1.4777 -3.4573 -3.4604 41 1 2.8631 -3.5742 -5.4866 42 1 1.9451 -5.0936 -5.3246 43 1 2.5554 -6.3740 -6.9447 44 1 2.0637 -7.1499 -5.4200 45 1 2.5684 -11.9805 -8.7563 46 1 5.9673 -5.3733 -4.1578 47 1 5.4812 -3.7576 -4.7289 48 1 4.0545 -5.2516 -0.9784 49 1 3.6584 -6.0979 -2.4934 50 1 6.4492 -5.8144 -1.2903 51 1 5.4470 -7.2852 -1.2527 52 1 6.0531 -6.6606 -2.8053 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 ZINC001723364418.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:08:08 Heat of formation + Delta-G solvation = -34.993784 kcal Electronic energy + Delta-G solvation = -27801.629802 eV Core-core repulsion = 24119.902970 eV Total energy + Delta-G solvation = -3681.726832 eV No. of doubly occupied orbitals = 63 Molecular weight (most abundant/longest-lived isotopes) = 340.193 amu Computer time = 1.70 seconds Orbital eigenvalues (eV) -40.39641 -38.29832 -37.09174 -34.98234 -34.02994 -32.41117 -31.42773 -30.19473 -29.20471 -27.94354 -27.50903 -25.49726 -24.52595 -23.90350 -22.83836 -21.93504 -20.15590 -18.83850 -18.25141 -17.64555 -16.79204 -16.15892 -15.83477 -15.56027 -15.48571 -15.05452 -14.86835 -14.78838 -14.49275 -14.08896 -13.92693 -13.79760 -13.54399 -13.41388 -13.06312 -12.90841 -12.65734 -12.52249 -12.29342 -12.16414 -12.09240 -11.86746 -11.76800 -11.55875 -11.44737 -11.23876 -11.17274 -11.02723 -11.01825 -10.75304 -10.67147 -10.54373 -10.25085 -9.79178 -9.42847 -9.29637 -9.04268 -8.71953 -8.45830 -8.20451 -6.47768 -6.19769 -3.90086 1.83324 2.97460 3.09583 3.24549 3.27332 3.30382 3.62720 3.76919 4.22037 4.42655 4.48758 4.53417 4.57989 4.79634 4.86792 4.90004 4.95903 5.02833 5.09890 5.12753 5.15180 5.19229 5.24865 5.28138 5.33279 5.42663 5.50185 5.52871 5.52906 5.54897 5.61783 5.62260 5.67856 5.70393 5.74653 5.76392 5.90717 5.92876 5.97833 6.12682 6.14685 6.17716 6.26046 6.38615 6.48049 6.48212 6.69959 6.75159 7.07787 7.60098 7.61998 8.36654 8.61240 9.25353 10.87532 Molecular weight = 340.19amu Principal moments of inertia in cm(-1) A = 0.015004 B = 0.002269 C = 0.002061 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1865.686965 B =12337.159820 C =13585.139109 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.193 4.193 2 C -0.127 4.127 3 C -0.120 4.120 4 C -0.084 4.084 5 C 0.066 3.934 6 N -0.684 5.684 7 C 0.095 3.905 8 H 0.044 0.956 9 C -0.124 4.124 10 C 0.108 3.892 11 N -0.616 5.616 12 C 0.536 3.464 13 O -0.517 6.517 14 C -0.151 4.151 15 S -0.060 6.060 16 C -0.535 4.535 17 H 0.071 0.929 18 C 0.074 3.926 19 N -0.868 5.868 20 Si 0.918 3.082 21 H -0.283 1.283 22 H -0.274 1.274 23 H -0.274 1.274 24 C 0.117 3.883 25 C -0.102 4.102 26 H 0.075 0.925 27 C -0.110 4.110 28 C -0.157 4.157 29 H 0.096 0.904 30 H 0.098 0.902 31 H 0.064 0.936 32 H 0.060 0.940 33 H 0.064 0.936 34 H 0.078 0.922 35 H 0.066 0.934 36 H 0.028 0.972 37 H 0.073 0.927 38 H 0.340 0.660 39 H 0.080 0.920 40 H 0.070 0.930 41 H 0.064 0.936 42 H 0.091 0.909 43 H 0.088 0.912 44 H 0.088 0.912 45 H 0.269 0.731 46 H 0.099 0.901 47 H 0.059 0.941 48 H 0.068 0.932 49 H 0.067 0.933 50 H 0.047 0.953 51 H 0.052 0.948 52 H 0.064 0.936 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -0.117 24.014 16.739 29.273 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.230 4.230 2 C -0.128 4.128 3 C -0.177 4.177 4 C -0.121 4.121 5 C -0.063 4.063 6 N -0.315 5.315 7 C -0.014 4.014 8 H 0.061 0.939 9 C -0.162 4.162 10 C -0.016 4.016 11 N -0.352 5.352 12 C 0.324 3.676 13 O -0.391 6.391 14 C -0.301 4.301 15 S 0.166 5.834 16 C -0.683 4.683 17 H 0.090 0.910 18 C -0.136 4.136 19 N -0.503 5.503 20 Si 0.744 3.256 21 H -0.208 1.208 22 H -0.199 1.199 23 H -0.200 1.200 24 C -0.007 4.007 25 C -0.121 4.121 26 H 0.094 0.906 27 C -0.147 4.147 28 C -0.215 4.215 29 H 0.114 0.886 30 H 0.117 0.883 31 H 0.083 0.917 32 H 0.079 0.921 33 H 0.083 0.917 34 H 0.097 0.903 35 H 0.085 0.915 36 H 0.046 0.954 37 H 0.091 0.909 38 H 0.159 0.841 39 H 0.099 0.901 40 H 0.088 0.912 41 H 0.082 0.918 42 H 0.110 0.890 43 H 0.106 0.894 44 H 0.107 0.893 45 H 0.078 0.922 46 H 0.117 0.883 47 H 0.077 0.923 48 H 0.087 0.913 49 H 0.085 0.915 50 H 0.066 0.934 51 H 0.071 0.929 52 H 0.083 0.917 Dipole moment (debyes) X Y Z Total from point charges 0.949 22.531 16.711 28.068 hybrid contribution -1.695 1.197 -0.709 2.193 sum -0.746 23.728 16.002 28.629 Atomic orbital electron populations 1.23768 0.95222 1.01954 1.02073 1.21707 0.95468 0.95611 0.99992 1.20881 1.00261 0.94458 1.02055 1.20625 1.00287 0.94366 0.96861 1.21946 0.98592 0.89479 0.96296 1.60483 1.28993 1.31334 1.10673 1.21484 0.94690 0.96352 0.88879 0.93850 1.22036 0.97503 1.01447 0.95249 1.21816 0.86603 0.96478 0.96733 1.48041 1.11853 1.18775 1.56550 1.20289 0.89872 0.79762 0.77650 1.90679 1.30593 1.59476 1.58358 1.23553 0.98671 1.05363 1.02531 1.85384 1.17961 1.10691 1.69320 1.22819 0.97456 1.06908 1.41147 0.91040 1.24987 0.97658 0.95307 0.95627 1.70013 1.29862 1.09878 1.40526 0.86235 0.82311 0.77765 0.79260 1.20830 1.19941 1.19958 1.21573 0.89614 0.98679 0.90814 1.21729 0.96743 0.97059 0.96596 0.90629 1.21461 0.98305 0.94977 0.99959 1.21866 0.97592 0.99583 1.02425 0.88589 0.88320 0.91704 0.92073 0.91671 0.90350 0.91512 0.95379 0.90861 0.84121 0.90109 0.91186 0.91784 0.89039 0.89418 0.89329 0.92161 0.88286 0.92253 0.91341 0.91459 0.93431 0.92874 0.91741 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.19 -1.26 13.90 29.01 0.40 -0.86 16 2 C -0.13 -0.82 5.31 -101.07 -0.54 -1.36 16 3 C -0.12 -0.64 10.12 36.00 0.36 -0.28 16 4 C -0.08 -0.54 5.64 -27.89 -0.16 -0.70 16 5 C 0.07 0.46 5.46 -3.82 -0.02 0.44 16 6 N -0.68 -5.22 5.54 -106.67 -0.59 -5.81 16 7 C 0.10 0.82 2.15 -67.67 -0.15 0.67 16 8 H 0.04 0.37 7.80 -51.93 -0.40 -0.03 16 9 C -0.12 -1.20 5.26 -26.39 -0.14 -1.33 16 10 C 0.11 1.19 6.37 -3.49 -0.02 1.17 16 11 N -0.62 -8.60 2.33 -172.50 -0.40 -9.01 16 12 C 0.54 10.12 7.72 -10.98 -0.08 10.04 16 13 O -0.52 -11.63 15.27 5.56 0.08 -11.55 16 14 C -0.15 -3.07 4.50 36.01 0.16 -2.91 16 15 S -0.06 -1.64 20.57 -107.50 -2.21 -3.85 16 16 C -0.54 -15.27 9.50 30.94 0.29 -14.98 16 17 H 0.07 1.93 7.80 -52.49 -0.41 1.52 16 18 C 0.07 2.10 5.49 -17.43 -0.10 2.01 16 19 N -0.87 -26.24 13.22 55.49 0.73 -25.50 16 20 Si 0.92 21.48 29.55 -169.99 -5.02 16.45 16 21 H -0.28 -6.51 7.11 56.52 0.40 -6.11 16 22 H -0.27 -6.26 7.11 56.52 0.40 -5.85 16 23 H -0.27 -6.25 7.11 56.52 0.40 -5.85 16 24 C 0.12 1.39 5.69 -3.71 -0.02 1.37 16 25 C -0.10 -0.96 2.05 -90.50 -0.19 -1.15 16 26 H 0.08 0.61 8.10 -51.93 -0.42 0.19 16 27 C -0.11 -1.21 4.74 -26.73 -0.13 -1.33 16 28 C -0.16 -1.85 9.68 37.16 0.36 -1.49 16 29 H 0.10 0.55 7.83 -52.49 -0.41 0.14 16 30 H 0.10 0.64 7.99 -52.48 -0.42 0.22 16 31 H 0.06 0.33 8.14 -51.93 -0.42 -0.09 16 32 H 0.06 0.30 7.78 -51.93 -0.40 -0.11 16 33 H 0.06 0.31 8.14 -51.93 -0.42 -0.11 16 34 H 0.08 0.57 7.99 -51.93 -0.41 0.15 16 35 H 0.07 0.40 8.14 -51.92 -0.42 -0.03 16 36 H 0.03 0.20 8.06 -51.93 -0.42 -0.22 16 37 H 0.07 0.57 7.41 -51.93 -0.38 0.19 16 38 H 0.34 2.33 8.26 -39.30 -0.32 2.01 16 39 H 0.08 0.92 6.45 -51.93 -0.33 0.58 16 40 H 0.07 0.61 7.42 -51.93 -0.39 0.23 16 41 H 0.06 0.62 8.14 -51.93 -0.42 0.20 16 42 H 0.09 1.05 5.93 -51.93 -0.31 0.75 16 43 H 0.09 1.67 7.93 -51.91 -0.41 1.25 16 44 H 0.09 1.69 7.70 -51.91 -0.40 1.29 16 45 H 0.27 7.57 8.31 -40.82 -0.34 7.23 16 46 H 0.10 1.39 5.20 -51.92 -0.27 1.12 16 47 H 0.06 0.63 8.14 -51.93 -0.42 0.21 16 48 H 0.07 0.70 7.45 -51.93 -0.39 0.31 16 49 H 0.07 0.91 6.12 -51.93 -0.32 0.60 16 50 H 0.05 0.47 8.14 -51.93 -0.42 0.04 16 51 H 0.05 0.64 8.14 -51.93 -0.42 0.21 16 52 H 0.06 0.91 6.46 -51.93 -0.34 0.58 16 LS Contribution 416.40 15.07 6.28 6.28 Total: -1.00 -32.74 416.40 -10.34 -43.08 By element: Atomic # 1 Polarization: 9.86 SS G_CDS: -9.26 Total: 0.60 kcal Atomic # 6 Polarization: -10.73 SS G_CDS: 0.05 Total: -10.68 kcal Atomic # 7 Polarization: -40.06 SS G_CDS: -0.26 Total: -40.32 kcal Atomic # 8 Polarization: -11.63 SS G_CDS: 0.08 Total: -11.55 kcal Atomic # 14 Polarization: 21.48 SS G_CDS: -5.02 Total: 16.45 kcal Atomic # 16 Polarization: -1.64 SS G_CDS: -2.21 Total: -3.85 kcal Total LS contribution 6.28 Total: 6.28 kcal Total: -32.74 -10.34 -43.08 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. ZINC001723364418.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 8.084 kcal (2) G-P(sol) polarization free energy of solvation -32.735 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -24.652 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.342 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -43.078 kcal (6) G-S(sol) free energy of system = (1) + (5) -34.994 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.70 seconds