Wall clock time and date at job start Tue Mar 31 2020 21:15:16 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.53000 * 1 3 3 C 1.53005 * 109.47071 * 2 1 4 4 C 1.52999 * 109.47305 * 119.99572 * 2 1 3 5 5 H 1.08994 * 109.47878 * 179.41571 * 4 2 1 6 6 N 1.46900 * 109.46782 * 59.41952 * 4 2 1 7 7 C 1.46902 * 111.00179 * 85.36621 * 6 4 2 8 8 C 1.53000 * 109.47158 * 175.55423 * 7 6 4 9 9 H 1.09007 * 109.46927 * 304.99349 * 8 7 6 10 10 O 1.42898 * 109.47434 * 64.99673 * 8 7 6 11 11 C 1.53006 * 109.46815 * 185.00084 * 8 7 6 12 12 N 1.46905 * 109.46702 * 179.97438 * 11 8 7 13 13 C 1.46899 * 111.00082 * 293.95136 * 12 11 8 14 14 C 1.46904 * 110.99527 * 169.99609 * 12 11 8 15 15 C 1.50692 * 109.46960 * 189.99624 * 14 12 11 16 16 H 1.08007 * 120.00475 * 0.02562 * 15 14 12 17 17 C 1.37883 * 120.00124 * 180.02562 * 15 14 12 18 18 N 1.31510 * 119.99597 * 0.02562 * 17 15 14 19 19 Si 1.86802 * 120.00038 * 180.02562 * 17 15 14 20 20 H 1.48498 * 109.47027 * 359.97438 * 19 17 15 21 21 H 1.48502 * 109.46995 * 120.00220 * 19 17 15 22 22 H 1.48499 * 109.47017 * 239.99914 * 19 17 15 23 23 C 1.50705 * 109.47198 * 299.42478 * 4 2 1 24 24 C 1.38238 * 119.97558 * 79.72502 * 23 4 2 25 25 C 1.38204 * 120.02380 * 179.73043 * 24 23 4 26 26 C 1.38458 * 119.97360 * 0.51808 * 25 24 23 27 27 F 1.35100 * 120.02882 * 179.73257 * 26 25 24 28 28 C 1.38463 * 119.94854 * 359.74784 * 26 25 24 29 29 C 1.38201 * 119.97061 * 0.02562 * 28 26 25 30 30 H 1.08995 * 109.47094 * 300.00248 * 1 2 3 31 31 H 1.09001 * 109.46644 * 60.00131 * 1 2 3 32 32 H 1.09000 * 109.47062 * 179.97438 * 1 2 3 33 33 H 1.08991 * 109.47773 * 239.99989 * 2 1 3 34 34 H 1.09002 * 109.47277 * 174.20790 * 3 2 1 35 35 H 1.09001 * 109.47060 * 294.21245 * 3 2 1 36 36 H 1.09007 * 109.46597 * 54.20926 * 3 2 1 37 37 H 1.00901 * 111.00070 * 209.32034 * 6 4 2 38 38 H 1.08998 * 109.47228 * 55.55157 * 7 6 4 39 39 H 1.09004 * 109.46825 * 295.54914 * 7 6 4 40 40 H 0.96697 * 113.99815 * 60.00448 * 10 8 7 41 41 H 1.08997 * 109.46752 * 299.99282 * 11 8 7 42 42 H 1.09001 * 109.46995 * 59.99785 * 11 8 7 43 43 H 1.08995 * 109.46923 * 60.00319 * 13 12 11 44 44 H 1.08999 * 109.46893 * 180.02562 * 13 12 11 45 45 H 1.09006 * 109.46881 * 300.00486 * 13 12 11 46 46 H 1.09000 * 109.47371 * 310.00126 * 14 12 11 47 47 H 1.09005 * 109.46525 * 70.00222 * 14 12 11 48 48 H 0.97006 * 119.99807 * 180.02562 * 18 17 15 49 49 H 1.08005 * 119.99126 * 359.97119 * 24 23 4 50 50 H 1.07999 * 120.01081 * 180.29389 * 25 24 23 51 51 H 1.07998 * 120.01427 * 179.97438 * 28 26 25 52 52 H 1.08005 * 119.98539 * 179.97438 * 29 28 26 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 6 2.0400 -0.7211 1.2493 5 1 3.1299 -0.7303 1.2441 6 7 1.5389 -2.1020 1.2564 7 6 2.4163 -2.9845 0.4757 8 6 1.9258 -4.4287 0.5973 9 1 1.8598 -4.7025 1.6503 10 8 0.6365 -4.5464 -0.0074 11 6 2.9089 -5.3625 -0.1117 12 7 2.4383 -6.7491 0.0055 13 6 2.4734 -7.2006 1.4030 14 6 3.2202 -7.6451 -0.8570 15 6 2.5623 -9.0001 -0.9013 16 1 1.6539 -9.1779 -0.3448 17 6 3.1200 -10.0130 -1.6524 18 7 4.2263 -9.7967 -2.3297 19 14 2.3041 -11.6925 -1.7079 20 1 1.0813 -11.6768 -0.8655 21 1 1.9363 -12.0181 -3.1093 22 1 3.2465 -12.7179 -1.1926 23 6 1.5494 -0.0025 2.4797 24 6 0.2569 -0.2026 2.9272 25 6 -0.1923 0.4516 4.0586 26 6 0.6493 1.3161 4.7378 27 9 0.2096 1.9604 5.8409 28 6 1.9420 1.5200 4.2855 29 6 2.3904 0.8610 3.1565 30 1 -0.3633 0.5138 0.8899 31 1 -0.3633 0.5138 -0.8900 32 1 -0.3633 -1.0277 -0.0005 33 1 1.8934 -0.5138 -0.8899 34 1 3.1251 1.4443 0.1037 35 1 1.7637 1.9256 -0.9373 36 1 1.5944 1.9854 0.8337 37 1 1.4306 -2.4426 2.2000 38 1 3.4351 -2.9135 0.8566 39 1 2.3979 -2.6815 -0.5712 40 1 0.6203 -4.3177 -0.9468 41 1 2.9750 -5.0887 -1.1647 42 1 3.8924 -5.2726 0.3496 43 1 1.8304 -6.5610 2.0076 44 1 2.1204 -8.2302 1.4611 45 1 3.4958 -7.1456 1.7770 46 1 3.2668 -7.2306 -1.8641 47 1 4.2296 -7.7441 -0.4577 48 1 4.6185 -10.5092 -2.8584 49 1 -0.4010 -0.8736 2.3947 50 1 -1.2007 0.2914 4.4108 51 1 2.5993 2.1944 4.8142 52 1 3.3986 1.0203 2.8032 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) F: (AM1): M.J.S. DEWAR AND E. G. ZOEBISCH, THEOCHEM, 180, 1 (1988). 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 ZINC000576089175.mol2 52 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 21:15:16 Heat of formation + Delta-G solvation = -72.891802 kcal Electronic energy + Delta-G solvation = -29520.664798 eV Core-core repulsion = 25434.225982 eV Total energy + Delta-G solvation = -4086.438816 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 338.211 amu Computer time = 1.58 seconds Orbital eigenvalues (eV) -49.19595 -40.27669 -38.24347 -36.92620 -35.43978 -33.62509 -32.50767 -31.57327 -31.01854 -29.36914 -27.33354 -27.15473 -26.35850 -24.85926 -23.61779 -22.93362 -22.37226 -21.26800 -20.20721 -19.63291 -17.93805 -16.88032 -16.72626 -16.43945 -15.98550 -15.84626 -15.55956 -15.44104 -15.06122 -14.71403 -14.58424 -14.31956 -14.07622 -13.96953 -13.71846 -13.54607 -13.40409 -13.28593 -12.72736 -12.53573 -12.49409 -12.23516 -12.04351 -11.97114 -11.87813 -11.36660 -11.31413 -11.20776 -11.08654 -11.03751 -10.88938 -10.75719 -10.68329 -10.64563 -10.28888 -9.96974 -9.56564 -9.41505 -9.10453 -8.87161 -8.57571 -8.44130 -7.20033 -6.02992 -4.03727 0.76944 0.88893 3.21254 3.31389 3.38747 3.44461 3.79627 3.93165 4.25391 4.39863 4.43555 4.55276 4.64668 4.75733 4.83407 4.87419 4.95407 5.01661 5.07576 5.14819 5.23317 5.30680 5.34986 5.47518 5.50323 5.55439 5.58316 5.63387 5.64957 5.70222 5.74277 5.77988 5.79998 5.89267 5.96043 6.05121 6.12745 6.15781 6.19284 6.28700 6.31715 6.39946 6.54079 6.76255 6.82350 7.17037 7.26746 7.32230 7.80727 8.04324 8.14275 8.20828 8.55446 9.04727 9.60543 11.08014 Molecular weight = 338.21amu Principal moments of inertia in cm(-1) A = 0.020535 B = 0.001955 C = 0.001881 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1363.198248 B =14321.277864 C =14881.832391 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.145 4.145 2 C -0.079 4.079 3 C -0.149 4.149 4 C 0.131 3.869 5 H 0.052 0.948 6 N -0.669 5.669 7 C 0.027 3.973 8 C 0.117 3.883 9 H 0.091 0.909 10 O -0.535 6.535 11 C 0.024 3.976 12 N -0.528 5.528 13 C 0.037 3.963 14 C 0.192 3.808 15 C -0.504 4.504 16 H 0.065 0.935 17 C 0.053 3.947 18 N -0.900 5.900 19 Si 0.898 3.102 20 H -0.271 1.271 21 H -0.284 1.284 22 H -0.285 1.285 23 C -0.094 4.094 24 C -0.069 4.069 25 C -0.157 4.157 26 C 0.094 3.906 27 F -0.143 7.143 28 C -0.157 4.157 29 C -0.084 4.084 30 H 0.053 0.947 31 H 0.044 0.956 32 H 0.074 0.926 33 H 0.072 0.928 34 H 0.051 0.949 35 H 0.055 0.945 36 H 0.055 0.945 37 H 0.343 0.657 38 H 0.036 0.964 39 H 0.077 0.923 40 H 0.369 0.631 41 H 0.070 0.930 42 H 0.026 0.974 43 H 0.046 0.954 44 H 0.074 0.926 45 H 0.007 0.993 46 H 0.028 0.972 47 H -0.012 1.012 48 H 0.260 0.740 49 H 0.146 0.854 50 H 0.135 0.865 51 H 0.133 0.867 52 H 0.129 0.871 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -4.207 26.078 9.911 28.214 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.202 4.202 2 C -0.098 4.098 3 C -0.207 4.207 4 C 0.021 3.979 5 H 0.070 0.930 6 N -0.301 5.301 7 C -0.103 4.103 8 C 0.057 3.943 9 H 0.108 0.892 10 O -0.339 6.339 11 C -0.106 4.106 12 N -0.252 5.252 13 C -0.112 4.112 14 C 0.067 3.933 15 C -0.542 4.542 16 H 0.083 0.917 17 C -0.148 4.148 18 N -0.540 5.540 19 Si 0.726 3.274 20 H -0.196 1.196 21 H -0.210 1.210 22 H -0.212 1.212 23 C -0.094 4.094 24 C -0.087 4.087 25 C -0.175 4.175 26 C 0.075 3.925 27 F -0.121 7.121 28 C -0.176 4.176 29 C -0.103 4.103 30 H 0.072 0.928 31 H 0.063 0.937 32 H 0.093 0.907 33 H 0.091 0.909 34 H 0.070 0.930 35 H 0.074 0.926 36 H 0.074 0.926 37 H 0.162 0.838 38 H 0.054 0.946 39 H 0.095 0.905 40 H 0.215 0.785 41 H 0.088 0.912 42 H 0.044 0.956 43 H 0.065 0.935 44 H 0.093 0.907 45 H 0.025 0.975 46 H 0.046 0.954 47 H 0.006 0.994 48 H 0.070 0.930 49 H 0.164 0.836 50 H 0.153 0.847 51 H 0.151 0.849 52 H 0.147 0.853 Dipole moment (debyes) X Y Z Total from point charges -5.189 25.893 9.758 28.153 hybrid contribution 1.645 -0.621 0.419 1.807 sum -3.545 25.272 10.177 27.474 Atomic orbital electron populations 1.21892 0.94437 1.03076 1.00808 1.21023 0.97136 0.94965 0.96641 1.21893 1.00633 0.96317 1.01875 1.20541 0.96718 0.88530 0.92077 0.93042 1.60480 1.50307 1.00368 1.18910 1.22805 0.94100 0.94073 0.99327 1.21856 0.84245 0.90336 0.97872 0.89166 1.86582 1.28504 1.91266 1.27522 1.22725 0.95431 0.91987 1.00428 1.62017 1.55635 1.00038 1.07542 1.22245 0.97722 1.01663 0.89538 1.19747 0.90785 0.91190 0.91544 1.21210 1.09769 0.94815 1.28454 0.91722 1.25007 0.95066 0.99142 0.95582 1.69970 1.19815 1.33143 1.31097 0.86434 0.79416 0.83718 0.77824 1.19598 1.21006 1.21174 1.19964 0.94401 0.98394 0.96654 1.21099 0.95430 0.97932 0.94282 1.21227 1.00628 0.98554 0.97128 1.17778 0.93721 0.95515 0.85479 1.91498 1.89867 1.79906 1.50870 1.21088 0.94785 1.02942 0.98754 1.20902 0.99771 0.95565 0.94027 0.92776 0.93688 0.90665 0.90921 0.93023 0.92598 0.92563 0.83758 0.94575 0.90475 0.78507 0.91179 0.95632 0.93491 0.90686 0.97472 0.95380 0.99384 0.93029 0.83598 0.84721 0.84922 0.85275 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.19 7.96 37.16 0.30 -0.90 16 2 C -0.08 -0.62 2.84 -90.62 -0.26 -0.87 16 3 C -0.15 -1.04 8.59 37.16 0.32 -0.72 16 4 C 0.13 1.03 2.13 -68.86 -0.15 0.89 16 5 H 0.05 0.41 7.58 -51.93 -0.39 0.02 16 6 N -0.67 -6.06 4.89 -107.95 -0.53 -6.58 16 7 C 0.03 0.27 5.17 -3.82 -0.02 0.25 16 8 C 0.12 1.57 2.98 -26.73 -0.08 1.49 16 9 H 0.09 1.26 6.53 -51.93 -0.34 0.92 16 10 O -0.54 -7.97 13.15 -35.23 -0.46 -8.43 16 11 C 0.02 0.40 4.62 -3.82 -0.02 0.39 16 12 N -0.53 -11.28 4.62 -244.36 -1.13 -12.41 16 13 C 0.04 0.73 9.14 60.06 0.55 1.27 16 14 C 0.19 5.01 4.72 -4.98 -0.02 4.98 16 15 C -0.50 -14.24 8.46 -34.44 -0.29 -14.53 16 16 H 0.06 1.83 7.12 -52.48 -0.37 1.45 16 17 C 0.05 1.51 5.46 -17.82 -0.10 1.41 16 18 N -0.90 -26.78 13.29 55.43 0.74 -26.05 16 19 Si 0.90 21.76 29.54 -169.99 -5.02 16.74 16 20 H -0.27 -6.47 7.11 56.52 0.40 -6.07 16 21 H -0.28 -6.82 7.11 56.52 0.40 -6.41 16 22 H -0.29 -6.83 7.11 56.52 0.40 -6.43 16 23 C -0.09 -0.72 3.57 -104.62 -0.37 -1.10 16 24 C -0.07 -0.55 6.80 -39.59 -0.27 -0.82 16 25 C -0.16 -1.28 10.03 -39.51 -0.40 -1.67 16 26 C 0.09 0.84 7.30 -39.42 -0.29 0.55 16 27 F -0.14 -1.41 17.27 2.25 0.04 -1.37 16 28 C -0.16 -1.20 10.03 -39.51 -0.40 -1.60 16 29 C -0.08 -0.61 8.86 -39.59 -0.35 -0.96 16 30 H 0.05 0.44 5.57 -51.93 -0.29 0.15 16 31 H 0.04 0.36 8.14 -51.93 -0.42 -0.07 16 32 H 0.07 0.71 7.59 -51.93 -0.39 0.32 16 33 H 0.07 0.62 7.58 -51.93 -0.39 0.23 16 34 H 0.05 0.33 8.14 -51.93 -0.42 -0.09 16 35 H 0.05 0.38 8.14 -51.93 -0.42 -0.05 16 36 H 0.06 0.41 7.02 -51.93 -0.36 0.05 16 37 H 0.34 3.19 7.91 -39.29 -0.31 2.88 16 38 H 0.04 0.36 7.95 -51.93 -0.41 -0.05 16 39 H 0.08 0.77 7.58 -51.93 -0.39 0.38 16 40 H 0.37 4.64 9.05 45.56 0.41 5.05 16 41 H 0.07 1.16 7.93 -51.93 -0.41 0.75 16 42 H 0.03 0.42 8.12 -51.93 -0.42 0.00 16 43 H 0.05 0.80 6.47 -51.93 -0.34 0.46 16 44 H 0.07 1.60 6.76 -51.93 -0.35 1.25 16 45 H 0.01 0.13 8.14 -51.93 -0.42 -0.29 16 46 H 0.03 0.77 7.93 -51.93 -0.41 0.36 16 47 H -0.01 -0.32 8.12 -51.93 -0.42 -0.75 16 48 H 0.26 7.29 8.31 -40.82 -0.34 6.95 16 49 H 0.15 1.16 6.85 -52.48 -0.36 0.80 16 50 H 0.13 0.92 8.06 -52.49 -0.42 0.50 16 51 H 0.13 0.81 8.06 -52.49 -0.42 0.39 16 52 H 0.13 0.73 8.04 -52.48 -0.42 0.31 16 LS Contribution 411.44 15.07 6.20 6.20 Total: -1.00 -30.77 411.44 -10.07 -40.84 By element: Atomic # 1 Polarization: 11.07 SS G_CDS: -8.06 Total: 3.01 kcal Atomic # 6 Polarization: -10.09 SS G_CDS: -1.84 Total: -11.93 kcal Atomic # 7 Polarization: -44.12 SS G_CDS: -0.92 Total: -45.04 kcal Atomic # 8 Polarization: -7.97 SS G_CDS: -0.46 Total: -8.43 kcal Atomic # 9 Polarization: -1.41 SS G_CDS: 0.04 Total: -1.37 kcal Atomic # 14 Polarization: 21.76 SS G_CDS: -5.02 Total: 16.74 kcal Total LS contribution 6.20 Total: 6.20 kcal Total: -30.77 -10.07 -40.84 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. ZINC000576089175.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 -32.054 kcal (2) G-P(sol) polarization free energy of solvation -30.768 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -62.822 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -10.069 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -40.838 kcal (6) G-S(sol) free energy of system = (1) + (5) -72.892 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.58 seconds