Wall clock time and date at job start Tue Mar 31 2020 19:23:33 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.53004 * 1 3 3 C 1.53002 * 109.46958 * 2 1 4 4 H 1.08994 * 109.48217 * 305.02862 * 3 2 1 5 5 C 1.53087 * 109.47836 * 185.00497 * 3 2 1 6 6 C 1.53181 * 109.16293 * 297.61460 * 5 3 2 7 7 H 1.08998 * 109.54436 * 62.97609 * 6 5 3 8 8 C 1.53000 * 109.54551 * 183.15268 * 6 5 3 9 9 C 1.50694 * 109.47392 * 185.00200 * 8 6 5 10 10 O 1.21288 * 120.00133 * 359.97438 * 9 8 6 11 11 N 1.34777 * 120.00350 * 180.02562 * 9 8 6 12 12 C 1.46496 * 120.00156 * 179.97438 * 11 9 8 13 13 C 1.50706 * 109.47027 * 180.02562 * 12 11 9 14 14 C 1.29552 * 117.66397 * 94.97376 * 13 12 11 15 15 C 1.50317 * 123.57685 * 180.78209 * 14 13 12 16 16 N 1.47008 * 109.21854 * 349.11343 * 15 14 13 17 17 C 1.36940 * 121.37480 * 222.89362 * 16 15 14 18 18 H 1.07999 * 119.99771 * 155.98353 * 17 16 15 19 19 C 1.38815 * 120.00024 * 335.98405 * 17 16 15 20 20 N 1.31615 * 119.99868 * 344.52144 * 19 17 16 21 21 Si 1.86797 * 119.99723 * 164.52270 * 19 17 16 22 22 H 1.48504 * 109.47278 * 90.00154 * 21 19 17 23 23 H 1.48507 * 109.46923 * 209.99813 * 21 19 17 24 24 H 1.48493 * 109.47392 * 329.99492 * 21 19 17 25 25 C 1.47335 * 117.25046 * 42.91594 * 16 15 14 26 26 C 1.49964 * 117.69380 * 275.00085 * 13 12 11 27 27 C 1.53181 * 109.03649 * 303.02090 * 6 5 3 28 28 C 1.53096 * 109.16091 * 56.98137 * 27 6 5 29 29 O 1.42903 * 109.41425 * 302.38776 * 28 27 6 30 30 H 1.08999 * 109.46972 * 180.02562 * 1 2 3 31 31 H 1.09006 * 109.46808 * 299.99932 * 1 2 3 32 32 H 1.09003 * 109.47306 * 59.99246 * 1 2 3 33 33 H 1.09000 * 109.46734 * 240.00751 * 2 1 3 34 34 H 1.08995 * 109.47259 * 120.00120 * 2 1 3 35 35 H 1.08998 * 109.52176 * 177.71487 * 5 3 2 36 36 H 1.08994 * 109.53054 * 57.52144 * 5 3 2 37 37 H 1.09000 * 109.46976 * 64.99878 * 8 6 5 38 38 H 1.09001 * 109.47283 * 305.00311 * 8 6 5 39 39 H 0.97000 * 119.99958 * 0.02562 * 11 9 8 40 40 H 1.09002 * 109.46710 * 60.00421 * 12 11 9 41 41 H 1.08996 * 109.47790 * 300.00189 * 12 11 9 42 42 H 1.08000 * 118.24529 * 0.86199 * 14 13 12 43 43 H 1.09003 * 109.52192 * 229.17714 * 15 14 13 44 44 H 1.09000 * 109.52148 * 109.05256 * 15 14 13 45 45 H 0.97002 * 120.00284 * 180.02562 * 20 19 17 46 46 H 1.09004 * 109.77979 * 57.81198 * 25 16 15 47 47 H 1.08995 * 109.77895 * 178.56560 * 25 16 15 48 48 H 1.09010 * 109.31574 * 40.74217 * 26 13 12 49 49 H 1.08997 * 109.16440 * 281.21642 * 26 13 12 50 50 H 1.09002 * 109.52220 * 176.88388 * 27 6 5 51 51 H 1.08993 * 109.52358 * 297.06860 * 27 6 5 52 52 H 1.09002 * 109.47226 * 182.40618 * 28 27 6 53 53 H 1.08997 * 109.48698 * 62.37933 * 28 27 6 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.5300 0.0000 0.0000 3 6 2.0400 1.4425 0.0000 4 1 1.6051 1.9818 -0.8414 5 6 3.5657 1.4446 -0.1259 6 6 4.1716 0.7376 1.0904 7 1 3.8518 -0.3043 1.1033 8 6 5.6982 0.8035 1.0129 9 6 6.2936 0.0035 2.1427 10 8 5.5683 -0.5632 2.9326 11 7 7.6320 -0.0814 2.2762 12 6 8.2108 -0.8587 3.3748 13 6 9.7145 -0.7928 3.2991 14 6 10.3138 0.1073 4.0126 15 6 11.8043 0.3017 4.0249 16 7 12.4451 -0.8262 3.3332 17 6 13.5712 -1.4230 3.8342 18 1 13.8064 -2.4426 3.5669 19 6 14.4103 -0.7175 4.6857 20 7 14.3077 0.5916 4.7749 21 14 15.6851 -1.6261 5.7048 22 1 16.9536 -1.7233 4.9388 23 1 15.9263 -0.8865 6.9698 24 1 15.1890 -2.9906 6.0161 25 6 11.8261 -1.2876 2.0783 26 6 10.4108 -1.7732 2.4030 27 6 3.6925 1.4383 2.3656 28 6 2.1618 1.4364 2.3952 29 8 1.6632 2.0814 1.2215 30 1 -0.3633 -1.0277 0.0005 31 1 -0.3633 0.5139 0.8901 32 1 -0.3634 0.5140 -0.8899 33 1 1.8933 -0.5137 -0.8901 34 1 1.8934 -0.5138 0.8899 35 1 3.9264 2.4724 -0.1668 36 1 3.8570 0.9189 -1.0352 37 1 6.0216 1.8414 1.0927 38 1 6.0318 0.3912 0.0606 39 1 8.2121 0.3727 1.6452 40 1 7.8877 -1.8966 3.2951 41 1 7.8769 -0.4464 4.3270 42 1 9.7102 0.7527 4.6335 43 1 12.1575 0.3452 5.0552 44 1 12.0546 1.2317 3.5144 45 1 14.8943 1.0848 5.3696 46 1 11.7791 -0.4642 1.3656 47 1 12.4103 -2.1062 1.6580 48 1 9.8452 -1.8804 1.4774 49 1 10.4673 -2.7391 2.9049 50 1 4.0743 0.9080 3.2380 51 1 4.0558 2.4659 2.3750 52 1 1.8149 1.9714 3.2792 53 1 1.8000 0.4087 2.4268 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 ZINC001571697462.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:23:33 Heat of formation + Delta-G solvation = -68.155139 kcal Electronic energy + Delta-G solvation = -29299.831807 eV Core-core repulsion = 25364.370171 eV Total energy + Delta-G solvation = -3935.461636 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 336.220 amu Computer time = 1.11 seconds Orbital eigenvalues (eV) -40.56358 -39.42256 -37.80899 -35.21451 -34.62010 -33.71822 -32.93530 -31.02100 -30.46785 -29.84394 -29.18548 -26.80046 -24.68532 -23.93526 -23.62327 -22.79135 -22.26626 -21.11589 -20.36999 -19.38838 -17.76420 -17.26543 -16.85768 -16.16440 -16.00140 -15.70615 -15.67483 -15.11990 -14.98382 -14.79028 -14.23715 -13.93960 -13.88548 -13.74073 -13.56125 -13.23818 -13.17703 -12.95985 -12.82631 -12.35800 -12.31477 -12.18688 -12.01387 -11.80159 -11.60568 -11.51392 -11.39046 -11.27141 -11.05898 -10.95064 -10.82879 -10.75792 -10.43693 -10.35627 -10.31005 -9.87468 -9.85777 -9.66120 -9.53092 -8.97120 -8.71985 -8.09666 -6.79140 -5.84133 -3.22904 2.55772 2.63759 3.35141 3.42306 3.44646 3.50324 3.62194 4.34502 4.41949 4.45838 4.49002 4.62785 4.65696 4.82306 4.98275 5.04079 5.08808 5.14039 5.18065 5.22146 5.27003 5.31472 5.35972 5.37695 5.45768 5.51913 5.55294 5.60426 5.64495 5.68197 5.77069 5.79290 5.89239 5.91477 5.95388 6.01965 6.07852 6.11890 6.15778 6.28528 6.42912 6.51512 6.68651 6.81680 7.01130 7.11818 7.20876 7.27939 7.57017 7.68620 7.77878 8.24292 8.33605 9.15075 9.77825 10.35609 11.26401 Molecular weight = 336.22amu Principal moments of inertia in cm(-1) A = 0.035597 B = 0.001770 C = 0.001753 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 786.391033 B =15813.143569 C =15970.741699 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.148 4.148 2 C -0.161 4.161 3 C 0.092 3.908 4 H 0.099 0.901 5 C -0.147 4.147 6 C -0.072 4.072 7 H 0.081 0.919 8 C -0.138 4.138 9 C 0.509 3.491 10 O -0.545 6.545 11 N -0.722 5.722 12 C 0.157 3.843 13 C -0.167 4.167 14 C -0.124 4.124 15 C 0.193 3.807 16 N -0.595 5.595 17 C -0.247 4.247 18 H 0.069 0.931 19 C -0.001 4.001 20 N -0.910 5.910 21 Si 0.905 3.095 22 H -0.287 1.287 23 H -0.289 1.289 24 H -0.280 1.280 25 C 0.144 3.856 26 C -0.090 4.090 27 C -0.155 4.155 28 C 0.058 3.942 29 O -0.377 6.377 30 H 0.052 0.948 31 H 0.059 0.941 32 H 0.053 0.947 33 H 0.071 0.929 34 H 0.073 0.927 35 H 0.072 0.928 36 H 0.073 0.927 37 H 0.093 0.907 38 H 0.093 0.907 39 H 0.401 0.599 40 H 0.064 0.936 41 H 0.064 0.936 42 H 0.094 0.906 43 H 0.093 0.907 44 H 0.041 0.959 45 H 0.254 0.746 46 H 0.022 0.978 47 H 0.057 0.943 48 H 0.052 0.948 49 H 0.064 0.936 50 H 0.093 0.907 51 H 0.064 0.936 52 H 0.090 0.910 53 H 0.053 0.947 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -23.788 2.130 -10.592 26.127 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.205 4.205 2 C -0.198 4.198 3 C 0.034 3.966 4 H 0.117 0.883 5 C -0.184 4.184 6 C -0.091 4.091 7 H 0.099 0.901 8 C -0.178 4.178 9 C 0.296 3.704 10 O -0.424 6.424 11 N -0.373 5.373 12 C 0.035 3.965 13 C -0.168 4.168 14 C -0.145 4.145 15 C 0.069 3.931 16 N -0.318 5.318 17 C -0.376 4.376 18 H 0.086 0.914 19 C -0.199 4.199 20 N -0.555 5.555 21 Si 0.734 3.266 22 H -0.214 1.214 23 H -0.215 1.215 24 H -0.205 1.205 25 C 0.017 3.983 26 C -0.129 4.129 27 C -0.193 4.193 28 C -0.018 4.018 29 O -0.296 6.296 30 H 0.071 0.929 31 H 0.078 0.922 32 H 0.072 0.928 33 H 0.089 0.911 34 H 0.091 0.909 35 H 0.091 0.909 36 H 0.092 0.908 37 H 0.111 0.889 38 H 0.111 0.889 39 H 0.235 0.765 40 H 0.083 0.917 41 H 0.082 0.918 42 H 0.112 0.888 43 H 0.111 0.889 44 H 0.059 0.941 45 H 0.064 0.936 46 H 0.040 0.960 47 H 0.076 0.924 48 H 0.071 0.929 49 H 0.082 0.918 50 H 0.111 0.889 51 H 0.083 0.917 52 H 0.108 0.892 53 H 0.072 0.928 Dipole moment (debyes) X Y Z Total from point charges -24.293 1.972 -9.890 26.303 hybrid contribution 1.254 -0.394 0.308 1.350 sum -23.039 1.577 -9.582 25.002 Atomic orbital electron populations 1.21701 0.95240 1.01404 1.02204 1.22217 0.96950 0.97501 1.03134 1.22182 0.93536 0.93195 0.87685 0.88271 1.22281 0.96025 1.01730 0.98393 1.20943 0.92671 1.00677 0.94766 0.90050 1.21754 0.96628 1.01359 0.98032 1.20680 0.82258 0.81156 0.86353 1.90652 1.59990 1.46999 1.44728 1.46256 1.06503 1.52837 1.31679 1.20436 0.91207 0.95604 0.89281 1.20529 0.92496 1.00619 1.03142 1.23308 1.00096 0.95517 0.95576 1.19760 0.83981 0.89974 0.99433 1.44118 1.28717 1.37738 1.21228 1.19171 0.98857 0.94856 1.24731 0.91371 1.24989 0.97895 0.95162 1.01845 1.70002 1.40620 1.04914 1.39954 0.86257 0.81369 0.79060 0.79871 1.21405 1.21543 1.20529 1.20951 0.91231 0.97075 0.89022 1.20755 0.97715 0.97570 0.96831 1.22352 0.96279 1.01163 0.99458 1.22583 0.93441 0.97910 0.87897 1.87965 1.74742 1.53206 1.13699 0.92914 0.92224 0.92819 0.91064 0.90884 0.90892 0.90848 0.88867 0.88872 0.76457 0.91721 0.91789 0.88801 0.88881 0.94105 0.93576 0.95962 0.92440 0.92903 0.91768 0.88893 0.91703 0.89158 0.92837 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.15 -1.01 9.81 37.16 0.36 -0.64 16 2 C -0.16 -1.15 4.74 -26.73 -0.13 -1.27 16 3 C 0.09 0.70 3.29 -26.68 -0.09 0.62 16 4 H 0.10 0.67 8.14 -51.93 -0.42 0.24 16 5 C -0.15 -1.05 4.93 -26.59 -0.13 -1.18 16 6 C -0.07 -0.63 1.61 -90.43 -0.15 -0.77 16 7 H 0.08 0.81 6.55 -51.93 -0.34 0.47 16 8 C -0.14 -1.17 5.18 -27.89 -0.14 -1.31 16 9 C 0.51 6.16 7.57 -10.99 -0.08 6.08 16 10 O -0.55 -8.39 14.00 5.55 0.08 -8.31 16 11 N -0.72 -8.73 5.55 -61.47 -0.34 -9.07 16 12 C 0.16 2.28 5.66 -5.19 -0.03 2.25 16 13 C -0.17 -2.77 5.23 -101.02 -0.53 -3.30 16 14 C -0.12 -2.30 9.47 -35.81 -0.34 -2.64 16 15 C 0.19 4.47 5.22 -5.10 -0.03 4.44 16 16 N -0.59 -14.53 3.01 -173.67 -0.52 -15.06 16 17 C -0.25 -6.66 9.71 -15.06 -0.15 -6.80 16 18 H 0.07 1.80 7.88 -52.49 -0.41 1.39 16 19 C 0.00 -0.03 5.01 -17.43 -0.09 -0.12 16 20 N -0.91 -25.81 11.08 55.49 0.61 -25.20 16 21 Si 0.90 21.46 29.58 -169.99 -5.03 16.43 16 22 H -0.29 -6.80 7.11 56.52 0.40 -6.40 16 23 H -0.29 -6.91 7.11 56.52 0.40 -6.51 16 24 H -0.28 -6.55 7.11 56.52 0.40 -6.14 16 25 C 0.14 2.88 6.58 -3.53 -0.02 2.85 16 26 C -0.09 -1.51 5.95 -28.18 -0.17 -1.68 16 27 C -0.16 -1.51 5.13 -26.59 -0.14 -1.64 16 28 C 0.06 0.53 6.52 37.21 0.24 0.77 16 29 O -0.38 -3.68 10.13 -35.23 -0.36 -4.04 16 30 H 0.05 0.33 8.14 -51.93 -0.42 -0.10 16 31 H 0.06 0.45 7.83 -51.93 -0.41 0.05 16 32 H 0.05 0.33 8.14 -51.93 -0.42 -0.10 16 33 H 0.07 0.44 8.14 -51.93 -0.42 0.02 16 34 H 0.07 0.60 4.98 -51.93 -0.26 0.34 16 35 H 0.07 0.52 8.14 -51.93 -0.42 0.10 16 36 H 0.07 0.43 8.14 -51.93 -0.42 0.01 16 37 H 0.09 0.69 8.14 -51.93 -0.42 0.27 16 38 H 0.09 0.63 8.14 -51.93 -0.42 0.21 16 39 H 0.40 4.21 8.46 -40.82 -0.35 3.87 16 40 H 0.06 0.97 8.11 -51.93 -0.42 0.55 16 41 H 0.06 0.97 7.78 -51.93 -0.40 0.57 16 42 H 0.09 1.57 7.74 -52.49 -0.41 1.16 16 43 H 0.09 2.41 6.41 -51.93 -0.33 2.08 16 44 H 0.04 1.01 8.14 -51.93 -0.42 0.58 16 45 H 0.25 6.99 8.31 -40.82 -0.34 6.65 16 46 H 0.02 0.44 8.14 -51.93 -0.42 0.01 16 47 H 0.06 1.12 7.96 -51.93 -0.41 0.71 16 48 H 0.05 0.76 8.14 -51.92 -0.42 0.33 16 49 H 0.06 1.13 8.14 -51.93 -0.42 0.71 16 50 H 0.09 1.10 6.07 -51.93 -0.32 0.78 16 51 H 0.06 0.62 8.14 -51.93 -0.42 0.20 16 52 H 0.09 0.75 8.14 -51.93 -0.42 0.33 16 53 H 0.05 0.48 6.23 -51.93 -0.32 0.16 16 LS Contribution 404.56 15.07 6.10 6.10 Total: -1.00 -30.46 404.56 -10.49 -40.95 By element: Atomic # 1 Polarization: 11.97 SS G_CDS: -9.43 Total: 2.54 kcal Atomic # 6 Polarization: -2.75 SS G_CDS: -1.60 Total: -4.34 kcal Atomic # 7 Polarization: -49.08 SS G_CDS: -0.25 Total: -49.33 kcal Atomic # 8 Polarization: -12.07 SS G_CDS: -0.28 Total: -12.35 kcal Atomic # 14 Polarization: 21.46 SS G_CDS: -5.03 Total: 16.43 kcal Total LS contribution 6.10 Total: 6.10 kcal Total: -30.46 -10.49 -40.95 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. ZINC001571697462.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 -27.205 kcal (2) G-P(sol) polarization free energy of solvation -30.461 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -57.665 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.490 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -40.950 kcal (6) G-S(sol) free energy of system = (1) + (5) -68.155 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.11 seconds