Wall clock time and date at job start Sat Apr 11 2020 06:11:36 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.53002 * 1 3 3 O 1.42895 * 109.47512 * 2 1 4 4 C 1.42898 * 114.00166 * 179.97438 * 3 2 1 5 5 H 1.08995 * 110.79891 * 26.70671 * 4 3 2 6 6 C 1.54943 * 110.73930 * 263.41882 * 4 3 2 7 7 N 1.48391 * 103.37303 * 204.38735 * 6 4 3 8 8 C 1.34778 * 127.17020 * 219.23990 * 7 6 4 9 9 O 1.21282 * 120.00036 * 179.79493 * 8 7 6 10 10 C 1.50701 * 119.99878 * 359.79084 * 8 7 6 11 11 C 1.53000 * 109.47013 * 174.24459 * 10 8 7 12 12 O 1.42907 * 112.84903 * 299.51466 * 11 10 8 13 13 C 1.42902 * 113.99437 * 296.51217 * 12 11 10 14 14 C 1.53783 * 113.61670 * 168.31500 * 11 10 8 15 15 C 1.53954 * 87.02087 * 220.04620 * 14 11 10 16 16 C 1.53785 * 113.61497 * 70.70707 * 11 10 8 17 17 C 1.47487 * 105.79382 * 39.06747 * 7 6 4 18 18 C 1.54005 * 107.24668 * 334.78088 * 17 7 6 19 19 H 1.09007 * 110.32383 * 120.10078 * 18 17 7 20 20 C 1.52996 * 110.32527 * 242.31668 * 18 17 7 21 21 N 1.46500 * 109.47264 * 293.54519 * 20 18 17 22 22 C 1.36939 * 120.00051 * 180.02562 * 21 20 18 23 23 H 1.07999 * 119.99544 * 0.02562 * 22 21 20 24 24 C 1.38812 * 119.99704 * 180.02562 * 22 21 20 25 25 N 1.31617 * 120.00457 * 359.97438 * 24 22 21 26 26 Si 1.86804 * 119.99613 * 179.97438 * 24 22 21 27 27 H 1.48503 * 109.47087 * 0.02562 * 26 24 22 28 28 H 1.48505 * 109.46753 * 120.00052 * 26 24 22 29 29 H 1.48496 * 109.47441 * 240.00249 * 26 24 22 30 30 H 1.09003 * 109.47072 * 300.00391 * 1 2 3 31 31 H 1.09001 * 109.46728 * 59.99381 * 1 2 3 32 32 H 1.08997 * 109.47551 * 179.97438 * 1 2 3 33 33 H 1.09003 * 109.47072 * 239.99609 * 2 1 3 34 34 H 1.08997 * 109.46851 * 120.00145 * 2 1 3 35 35 H 1.08994 * 110.66459 * 322.72720 * 6 4 3 36 36 H 1.09005 * 110.65978 * 85.84383 * 6 4 3 37 37 H 1.08994 * 109.47585 * 294.24431 * 10 8 7 38 38 H 1.09004 * 109.46899 * 54.24756 * 10 8 7 39 39 H 1.08995 * 109.47306 * 63.27103 * 13 12 11 40 40 H 1.09006 * 109.46810 * 183.27195 * 13 12 11 41 41 H 1.08997 * 109.47352 * 303.27028 * 13 12 11 42 42 H 1.09002 * 113.61396 * 334.57046 * 14 11 10 43 43 H 1.09004 * 113.61103 * 105.46630 * 14 11 10 44 44 H 1.08993 * 113.54141 * 270.81359 * 15 14 11 45 45 H 1.08998 * 113.73878 * 139.93598 * 15 14 11 46 46 H 1.08993 * 113.61279 * 254.53021 * 16 11 10 47 47 H 1.09003 * 113.61262 * 25.43341 * 16 11 10 48 48 H 1.08996 * 109.90925 * 215.32452 * 17 7 6 49 49 H 1.09001 * 109.91141 * 94.23993 * 17 7 6 50 50 H 1.08998 * 109.47204 * 53.55137 * 20 18 17 51 51 H 1.09008 * 109.46622 * 173.55050 * 20 18 17 52 52 H 0.96998 * 119.99993 * 0.02562 * 21 20 18 53 53 H 0.96996 * 120.00084 * 179.97438 * 25 24 22 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 8 2.0064 1.3472 0.0000 4 6 3.4310 1.4600 0.0006 5 1 3.8888 0.5831 -0.4570 6 6 3.9642 1.6690 1.4403 7 7 5.2659 2.3494 1.2287 8 6 6.4119 2.1441 1.9077 9 8 7.4075 2.7748 1.6216 10 6 6.4560 1.1310 3.0224 11 6 7.8936 0.9901 3.5268 12 8 8.7970 0.5475 2.5118 13 6 8.5130 -0.7606 2.0114 14 6 8.0105 0.2018 4.8421 15 6 9.0748 1.2289 5.2695 16 6 8.3960 2.2170 4.3062 17 6 5.0486 3.2952 0.1180 18 6 3.8783 2.7514 -0.7224 19 1 3.0621 3.4736 -0.7470 20 6 4.3481 2.4300 -2.1425 21 7 4.7202 3.6710 -2.8264 22 6 5.1734 3.6317 -4.1180 23 1 5.2561 2.6857 -4.6326 24 6 5.5265 4.8075 -4.7658 25 7 5.4252 5.9603 -4.1388 26 14 6.1454 4.7537 -6.5275 27 1 6.1654 3.3468 -7.0024 28 1 5.2438 5.5552 -7.3935 29 1 7.5183 5.3163 -6.5873 30 1 -0.3633 0.5139 0.8900 31 1 -0.3633 0.5139 -0.8900 32 1 -0.3634 -1.0276 -0.0005 33 1 1.8934 -0.5139 -0.8900 34 1 1.8933 -0.5138 0.8900 35 1 3.2889 2.3055 2.0119 36 1 4.1071 0.7111 1.9403 37 1 5.8150 1.4616 3.8396 38 1 6.1049 0.1678 2.6520 39 1 7.5360 -0.7606 1.5282 40 1 9.2763 -1.0462 1.2875 41 1 8.5103 -1.4727 2.8365 42 1 7.1065 0.2308 5.4504 43 1 8.4029 -0.8070 4.7130 44 1 10.0874 0.9557 4.9727 45 1 9.0053 1.5228 6.3168 46 1 9.0990 2.8326 3.7453 47 1 7.5975 2.7980 4.7678 48 1 5.9473 3.3603 -0.4952 49 1 4.7965 4.2795 0.5126 50 1 5.2118 1.7666 -2.0981 51 1 3.5421 1.9412 -2.6899 52 1 4.6456 4.5206 -2.3643 53 1 5.6723 6.7819 -4.5913 RHF calculation, no. of doubly occupied orbitals= 66 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 ZINC001264526883.mol2 53 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:11:36 Heat of formation + Delta-G solvation = -83.196733 kcal Electronic energy + Delta-G solvation = -31264.465716 eV Core-core repulsion = 27137.895777 eV Total energy + Delta-G solvation = -4126.569940 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.77 seconds Orbital eigenvalues (eV) -40.66633 -39.65117 -38.14668 -36.50863 -35.69933 -34.60641 -32.32309 -31.35643 -30.72486 -29.30181 -28.20429 -27.23080 -25.78322 -25.43510 -23.86779 -23.01855 -22.84935 -21.80000 -20.09141 -19.62057 -19.19047 -17.38635 -17.20295 -16.69474 -16.56266 -16.11941 -15.72599 -15.39063 -15.23389 -14.92992 -14.64601 -14.38469 -14.19908 -14.12793 -13.71915 -13.53098 -13.35698 -13.19757 -12.91978 -12.70697 -12.55245 -12.45780 -12.27014 -11.94262 -11.75244 -11.66241 -11.51485 -11.50787 -11.19720 -11.11250 -11.01162 -11.00121 -10.73937 -10.71405 -10.39594 -10.29379 -10.04184 -9.82855 -9.72935 -9.50312 -9.43303 -8.85622 -8.51946 -6.89706 -5.72532 -3.01709 2.76783 2.89718 3.47005 3.51340 3.52577 3.54632 3.65053 4.16489 4.36733 4.45382 4.49263 4.55945 4.63999 4.68424 4.82494 4.93107 5.01521 5.02022 5.08969 5.11000 5.13594 5.21293 5.22063 5.28411 5.31674 5.37965 5.44419 5.45392 5.50812 5.59314 5.59653 5.66346 5.68025 5.73417 5.78554 5.95147 5.98978 6.00970 6.04594 6.09375 6.12013 6.25701 6.36858 6.42208 6.70808 6.78240 7.15456 7.20657 7.64239 7.73647 8.33030 8.45557 9.57097 10.13054 10.48048 11.50288 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.008285 B = 0.003311 C = 0.002575 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 3378.868087 B = 8455.122375 C =10871.242109 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.151 4.151 2 C 0.055 3.945 3 O -0.369 6.369 4 C 0.084 3.916 5 H 0.078 0.922 6 C 0.082 3.918 7 N -0.623 5.623 8 C 0.530 3.470 9 O -0.512 6.512 10 C -0.154 4.154 11 C 0.110 3.890 12 O -0.343 6.343 13 C 0.026 3.974 14 C -0.165 4.165 15 C -0.136 4.136 16 C -0.116 4.116 17 C 0.114 3.886 18 C -0.104 4.104 19 H 0.098 0.902 20 C 0.164 3.836 21 N -0.731 5.731 22 C -0.240 4.240 23 H 0.070 0.930 24 C -0.013 4.013 25 N -0.931 5.931 26 Si 0.907 3.093 27 H -0.281 1.281 28 H -0.292 1.292 29 H -0.291 1.291 30 H 0.062 0.938 31 H 0.065 0.935 32 H 0.068 0.932 33 H 0.050 0.950 34 H 0.045 0.955 35 H 0.082 0.918 36 H 0.075 0.925 37 H 0.102 0.898 38 H 0.089 0.911 39 H 0.038 0.962 40 H 0.087 0.913 41 H 0.041 0.959 42 H 0.089 0.911 43 H 0.072 0.928 44 H 0.083 0.917 45 H 0.068 0.932 46 H 0.085 0.915 47 H 0.089 0.911 48 H 0.094 0.906 49 H 0.067 0.933 50 H 0.019 0.981 51 H 0.023 0.977 52 H 0.385 0.615 53 H 0.255 0.745 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -3.545 -16.543 21.226 27.144 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.209 4.209 2 C -0.022 4.022 3 O -0.288 6.288 4 C 0.024 3.976 5 H 0.096 0.904 6 C -0.042 4.042 7 N -0.360 5.360 8 C 0.318 3.682 9 O -0.386 6.386 10 C -0.194 4.194 11 C 0.069 3.931 12 O -0.261 6.261 13 C -0.070 4.070 14 C -0.202 4.202 15 C -0.173 4.173 16 C -0.153 4.153 17 C -0.008 4.008 18 C -0.124 4.124 19 H 0.116 0.884 20 C 0.037 3.963 21 N -0.375 5.375 22 C -0.371 4.371 23 H 0.088 0.912 24 C -0.209 4.209 25 N -0.580 5.580 26 Si 0.738 3.262 27 H -0.206 1.206 28 H -0.219 1.219 29 H -0.218 1.218 30 H 0.081 0.919 31 H 0.084 0.916 32 H 0.087 0.913 33 H 0.069 0.931 34 H 0.063 0.937 35 H 0.100 0.900 36 H 0.094 0.906 37 H 0.120 0.880 38 H 0.108 0.892 39 H 0.057 0.943 40 H 0.105 0.895 41 H 0.059 0.941 42 H 0.107 0.893 43 H 0.091 0.909 44 H 0.102 0.898 45 H 0.087 0.913 46 H 0.104 0.896 47 H 0.107 0.893 48 H 0.112 0.888 49 H 0.085 0.915 50 H 0.036 0.964 51 H 0.041 0.959 52 H 0.218 0.782 53 H 0.065 0.935 Dipole moment (debyes) X Y Z Total from point charges -2.811 -15.998 20.885 26.457 hybrid contribution 0.006 -0.234 -0.866 0.897 sum -2.804 -16.231 20.018 25.924 Atomic orbital electron populations 1.21739 0.93758 1.02059 1.03298 1.22531 0.95645 0.84029 0.99969 1.87973 1.19944 1.26148 1.94701 1.23146 0.81197 0.99129 0.94102 0.90439 1.22970 0.85414 0.97682 0.98100 1.48818 1.09528 1.41681 1.35956 1.20003 0.81456 0.83915 0.82785 1.90776 1.32766 1.48491 1.66578 1.22200 0.99282 0.99532 0.98398 1.22222 0.87374 0.96880 0.86620 1.88074 1.51548 1.31984 1.54520 1.22752 1.01362 0.84591 0.98318 1.23602 1.00317 0.99123 0.97190 1.22813 0.98664 0.96622 0.99244 1.22903 1.00938 0.95635 0.95868 1.22135 0.99660 0.91771 0.87202 1.22406 0.97606 0.96314 0.96025 0.88379 1.20671 0.96346 0.90761 0.88560 1.42333 1.78899 1.03197 1.13069 1.19233 1.35603 0.93062 0.89182 0.91203 1.24740 1.02307 0.95357 0.98468 1.69665 1.54473 0.97119 1.36693 0.86121 0.78198 0.77176 0.84711 1.20598 1.21879 1.21823 0.91934 0.91590 0.91296 0.93127 0.93666 0.89980 0.90632 0.87972 0.89222 0.94306 0.89461 0.94078 0.89281 0.90908 0.89799 0.91317 0.89638 0.89289 0.88802 0.91494 0.96361 0.95917 0.78151 0.93534 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.15 -1.25 10.28 37.16 0.38 -0.87 16 2 C 0.05 0.52 6.50 37.16 0.24 0.77 16 3 O -0.37 -4.70 10.08 -35.23 -0.36 -5.06 16 4 C 0.08 1.08 3.11 -24.86 -0.08 1.00 16 5 H 0.08 0.97 7.55 -51.93 -0.39 0.58 16 6 C 0.08 0.91 6.68 -1.97 -0.01 0.89 16 7 N -0.62 -8.95 3.30 -169.54 -0.56 -9.51 16 8 C 0.53 7.98 7.67 -10.98 -0.08 7.89 16 9 O -0.51 -10.00 15.27 -7.92 -0.12 -10.12 16 10 C -0.15 -1.59 3.92 -27.88 -0.11 -1.70 16 11 C 0.11 1.21 0.38 -89.80 -0.03 1.17 16 12 O -0.34 -4.66 9.82 -48.71 -0.48 -5.14 16 13 C 0.03 0.27 10.55 101.05 1.07 1.33 16 14 C -0.16 -1.31 7.22 -25.83 -0.19 -1.49 16 15 C -0.14 -1.18 8.02 -25.82 -0.21 -1.38 16 16 C -0.12 -1.28 6.58 -25.91 -0.17 -1.45 16 17 C 0.11 1.96 6.30 -2.98 -0.02 1.94 16 18 C -0.10 -1.72 2.83 -89.25 -0.25 -1.97 16 19 H 0.10 1.69 8.08 -51.92 -0.42 1.27 16 20 C 0.16 3.37 5.34 -4.04 -0.02 3.35 16 21 N -0.73 -18.86 5.33 -59.91 -0.32 -19.18 16 22 C -0.24 -6.87 9.99 -15.06 -0.15 -7.02 16 23 H 0.07 1.95 7.83 -52.49 -0.41 1.54 16 24 C -0.01 -0.37 5.50 -17.43 -0.10 -0.47 16 25 N -0.93 -28.76 13.06 55.49 0.72 -28.03 16 26 Si 0.91 22.96 29.55 -169.99 -5.02 17.94 16 27 H -0.28 -6.87 7.11 56.52 0.40 -6.47 16 28 H -0.29 -7.37 7.11 56.52 0.40 -6.97 16 29 H -0.29 -7.43 7.11 56.52 0.40 -7.03 16 30 H 0.06 0.50 8.14 -51.93 -0.42 0.08 16 31 H 0.06 0.61 8.14 -51.93 -0.42 0.19 16 32 H 0.07 0.43 8.14 -51.93 -0.42 0.01 16 33 H 0.05 0.48 8.12 -51.93 -0.42 0.06 16 34 H 0.05 0.35 8.14 -51.93 -0.42 -0.08 16 35 H 0.08 0.86 8.14 -51.93 -0.42 0.44 16 36 H 0.08 0.62 7.09 -51.93 -0.37 0.26 16 37 H 0.10 0.81 8.14 -51.93 -0.42 0.39 16 38 H 0.09 0.77 6.50 -51.93 -0.34 0.44 16 39 H 0.04 0.39 6.94 -51.93 -0.36 0.03 16 40 H 0.09 0.83 8.14 -51.93 -0.42 0.41 16 41 H 0.04 0.32 6.73 -51.93 -0.35 -0.03 16 42 H 0.09 0.55 8.10 -51.93 -0.42 0.13 16 43 H 0.07 0.54 6.84 -51.93 -0.36 0.19 16 44 H 0.08 0.77 8.14 -51.93 -0.42 0.35 16 45 H 0.07 0.49 8.14 -51.93 -0.42 0.06 16 46 H 0.08 1.16 8.14 -51.93 -0.42 0.73 16 47 H 0.09 0.86 8.10 -51.93 -0.42 0.44 16 48 H 0.09 1.97 7.48 -51.93 -0.39 1.58 16 49 H 0.07 1.14 8.14 -51.93 -0.42 0.71 16 50 H 0.02 0.39 8.14 -51.93 -0.42 -0.03 16 51 H 0.02 0.47 8.14 -51.92 -0.42 0.05 16 52 H 0.38 10.40 7.51 -40.82 -0.31 10.10 16 53 H 0.25 7.60 8.31 -40.82 -0.34 7.26 16 LS Contribution 419.61 15.07 6.32 6.32 Total: -1.00 -34.99 419.61 -9.12 -44.11 By element: Atomic # 1 Polarization: 16.25 SS G_CDS: -9.58 Total: 6.67 kcal Atomic # 6 Polarization: 1.73 SS G_CDS: 0.27 Total: 2.00 kcal Atomic # 7 Polarization: -56.57 SS G_CDS: -0.15 Total: -56.73 kcal Atomic # 8 Polarization: -19.36 SS G_CDS: -0.95 Total: -20.32 kcal Atomic # 14 Polarization: 22.96 SS G_CDS: -5.02 Total: 17.94 kcal Total LS contribution 6.32 Total: 6.32 kcal Total: -34.99 -9.12 -44.11 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. ZINC001264526883.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 -39.084 kcal (2) G-P(sol) polarization free energy of solvation -34.989 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -74.073 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.123 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -44.113 kcal (6) G-S(sol) free energy of system = (1) + (5) -83.197 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.77 seconds