Wall clock time and date at job start Tue Mar 9 2021 15:15:50 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= 0 * 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 O 2 2 C 1.42906 * 1 3 3 Si 1.86296 * 109.47126 * 2 1 4 4 H 1.48504 * 109.47239 * 299.99608 * 3 2 1 5 5 H 1.48502 * 109.47352 * 59.99647 * 3 2 1 6 6 H 1.48505 * 109.46927 * 179.97438 * 3 2 1 7 7 C 1.50702 * 109.47242 * 239.99329 * 2 1 3 8 8 C 1.38638 * 119.93570 * 74.98152 * 7 2 1 9 9 F 1.35099 * 119.87894 * 0.04271 * 8 7 2 10 10 C 1.38671 * 120.23549 * 179.74266 * 8 7 2 11 11 C 1.37815 * 120.12108 * 0.56052 * 10 8 7 12 12 C 1.39671 * 119.87688 * 359.44200 * 11 10 8 13 13 C 1.47724 * 120.11860 * 180.30566 * 12 11 10 14 14 O 1.21572 * 119.99882 * 179.97438 * 13 12 11 15 15 N 1.34771 * 120.00124 * 359.97438 * 13 12 11 16 16 C 1.46502 * 120.00266 * 180.02562 * 15 13 12 17 17 C 1.52560 * 108.99424 * 276.16231 * 16 15 13 18 18 C 1.54011 * 109.45738 * 162.04442 * 17 16 15 19 19 Si 1.86517 * 107.16034 * 58.32144 * 18 17 16 20 20 C 1.86518 * 110.13757 * 75.81571 * 19 18 17 21 21 C 1.53975 * 107.20634 * 164.96093 * 20 19 18 22 22 C 1.52542 * 109.47508 * 301.76670 * 21 20 19 23 23 C 1.52551 * 111.91409 * 77.26631 * 22 21 20 24 24 C 1.53954 * 109.49305 * 282.74567 * 23 22 21 25 25 C 1.86536 * 108.14374 * 315.42978 * 19 18 17 26 26 C 1.52558 * 108.99406 * 153.84063 * 16 15 13 27 27 C 1.37851 * 119.94046 * 254.99987 * 7 2 1 28 28 H 0.96703 * 113.99156 * 60.00357 * 1 2 3 29 29 H 1.09008 * 109.46768 * 119.99830 * 2 1 3 30 30 H 1.08004 * 119.94312 * 180.30137 * 10 8 7 31 31 H 1.08002 * 120.06547 * 179.72341 * 11 10 8 32 32 H 0.97001 * 119.99997 * 0.02562 * 15 13 12 33 33 H 1.08996 * 108.97468 * 35.00380 * 16 15 13 34 34 H 1.08991 * 109.46923 * 282.12380 * 17 16 15 35 35 H 1.09004 * 109.46098 * 42.11406 * 17 16 15 36 36 H 1.09006 * 109.96945 * 298.89514 * 18 17 16 37 37 H 1.08997 * 109.96694 * 177.75132 * 18 17 16 38 38 H 1.08994 * 109.90669 * 284.38928 * 20 19 18 39 39 H 1.09003 * 109.91013 * 45.54089 * 20 19 18 40 40 H 1.09001 * 109.47282 * 181.76027 * 21 20 19 41 41 H 1.08998 * 109.46726 * 61.76675 * 21 20 19 42 42 H 1.09008 * 108.99111 * 316.64205 * 22 21 20 43 43 H 1.08998 * 108.92815 * 197.85269 * 22 21 20 44 44 H 1.09007 * 109.46607 * 162.74575 * 23 22 21 45 45 H 1.08993 * 109.46932 * 42.75300 * 23 22 21 46 46 H 1.09002 * 109.90492 * 298.77804 * 24 23 22 47 47 H 1.08997 * 109.91006 * 177.63103 * 24 23 22 48 48 H 1.08994 * 109.91078 * 285.11402 * 25 19 18 49 49 H 1.09005 * 109.97023 * 164.00228 * 25 19 18 50 50 H 1.09002 * 109.47190 * 317.90136 * 26 16 15 51 51 H 1.09000 * 109.47087 * 77.89783 * 26 16 15 52 52 H 1.07998 * 120.06454 * 0.02562 * 27 7 2 Note: An asterisk (*) indicates this parameter was or will be optimized. Cartesian coordinates (angstroms) NO. Atom X Y Z 1 8 0.0000 0.0000 0.0000 2 6 1.4291 0.0000 0.0000 3 14 2.0501 1.7564 0.0000 4 1 1.5551 2.4565 -1.2126 5 1 1.5550 2.4565 1.2125 6 1 3.5351 1.7564 0.0006 7 6 1.9314 -0.7106 -1.2304 8 6 1.8685 -0.0838 -2.4654 9 9 1.3576 1.1631 -2.5621 10 6 2.3339 -0.7298 -3.6008 11 6 2.8512 -2.0039 -3.5086 12 6 2.9169 -2.6399 -2.2668 13 6 3.4765 -4.0029 -2.1608 14 8 3.5331 -4.5543 -1.0788 15 7 3.9244 -4.6351 -3.2636 16 6 4.4789 -5.9871 -3.1586 17 6 5.9474 -5.8916 -2.7565 18 6 6.4192 -7.2521 -2.2101 19 14 6.1438 -8.5118 -3.5577 20 6 7.4517 -8.3381 -4.8761 21 6 7.4039 -9.6144 -5.7361 22 6 7.6416 -10.8335 -4.8506 23 6 6.3996 -11.2035 -4.0457 24 6 6.2512 -10.2369 -2.8566 25 6 4.4527 -8.2200 -4.2888 26 6 4.3015 -6.7017 -4.4947 27 6 2.4523 -1.9826 -1.1258 28 1 -0.3932 0.4417 -0.7651 29 1 1.7924 -0.5139 0.8901 30 1 2.2871 -0.2352 -4.5598 31 1 3.2093 -2.5079 -4.3941 32 1 3.8796 -4.1950 -4.1268 33 1 3.9351 -6.5344 -2.3887 34 1 6.5454 -5.6209 -3.6265 35 1 6.0639 -5.1308 -1.9846 36 1 5.8388 -7.5202 -1.3272 37 1 7.4784 -7.2046 -1.9572 38 1 7.2415 -7.4656 -5.4946 39 1 8.4338 -8.2386 -4.4139 40 1 8.1776 -9.5654 -6.5023 41 1 6.4264 -9.6957 -6.2115 42 1 8.4591 -10.6181 -4.1624 43 1 7.9213 -11.6784 -5.4798 44 1 6.4979 -12.2236 -3.6745 45 1 5.5188 -11.1323 -4.6837 46 1 7.1189 -10.3181 -2.2019 47 1 5.3427 -10.4689 -2.3009 48 1 3.6851 -8.5800 -3.6039 49 1 4.3670 -8.7348 -5.2458 50 1 3.3101 -6.4861 -4.8932 51 1 5.0589 -6.3529 -5.1967 52 1 2.5002 -2.4702 -0.1634 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) 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=0 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 m_22_15739362_12496548.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 9 2021 15:15:50 Heat of formation + Delta-G solvation = -139.103251 kcal Electronic energy + Delta-G solvation = -30815.288248 eV Core-core repulsion = 26623.639390 eV Total energy + Delta-G solvation = -4191.648857 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 365.174 amu Computer time = 0.77 seconds Orbital eigenvalues (eV) -50.47257 -41.46613 -39.88470 -37.95726 -37.57207 -37.40888 -34.28145 -32.80907 -32.15292 -31.99941 -31.83245 -27.73761 -27.07531 -26.66290 -24.54700 -23.71596 -23.52935 -23.34919 -21.57869 -21.32076 -21.12688 -19.23863 -18.90625 -18.10318 -17.84158 -17.23853 -16.98238 -16.50330 -16.30672 -16.18138 -15.93333 -15.55704 -15.47547 -15.20840 -15.15712 -15.07463 -14.93974 -14.73423 -14.43570 -14.08897 -13.85842 -13.79519 -13.49712 -13.39300 -13.35996 -13.20572 -12.94776 -12.85785 -12.64085 -12.40098 -12.37884 -12.24344 -12.10662 -12.06565 -11.98188 -11.67099 -11.52830 -11.24963 -11.04271 -10.66968 -10.42532 -10.27678 -10.20769 -10.18207 -10.09177 -9.86002 -0.55663 -0.07298 1.08998 1.46440 1.49969 1.63619 1.71428 1.73871 1.79778 1.94256 2.04188 2.54490 2.55345 2.90144 3.22034 3.45331 3.49684 3.52806 3.67680 3.81527 3.92793 3.97401 4.03115 4.04133 4.11218 4.20841 4.26757 4.28388 4.37570 4.40193 4.52411 4.52666 4.55196 4.68703 4.69110 4.71635 4.77555 4.80040 4.84853 4.92810 4.93128 4.99492 5.02838 5.12455 5.14127 5.15871 5.22101 5.26213 5.30455 5.34574 5.38839 5.39709 5.46526 5.51735 5.55824 6.40413 6.42085 7.14455 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 O -0.532 6.532 2 C -0.080 4.080 3 Si 0.936 3.064 4 H -0.229 1.229 5 H -0.251 1.251 6 H -0.251 1.251 7 C -0.157 4.157 8 C 0.116 3.884 9 F -0.123 7.123 10 C -0.153 4.153 11 C -0.056 4.056 12 C -0.144 4.144 13 C 0.560 3.440 14 O -0.528 6.528 15 N -0.723 5.723 16 C 0.135 3.865 17 C -0.117 4.117 18 C -0.405 4.405 19 Si 0.997 3.003 20 C -0.405 4.405 21 C -0.109 4.109 22 C -0.133 4.133 23 C -0.109 4.109 24 C -0.408 4.408 25 C -0.403 4.403 26 C -0.116 4.116 27 C -0.010 4.010 28 H 0.376 0.624 29 H 0.125 0.875 30 H 0.147 0.853 31 H 0.139 0.861 32 H 0.399 0.601 33 H 0.093 0.907 34 H 0.066 0.934 35 H 0.068 0.932 36 H 0.079 0.921 37 H 0.078 0.922 38 H 0.073 0.927 39 H 0.075 0.925 40 H 0.061 0.939 41 H 0.062 0.938 42 H 0.065 0.935 43 H 0.061 0.939 44 H 0.061 0.939 45 H 0.063 0.937 46 H 0.075 0.925 47 H 0.074 0.926 48 H 0.080 0.920 49 H 0.082 0.918 50 H 0.071 0.929 51 H 0.071 0.929 52 H 0.151 0.849 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 1.530 -0.393 -5.061 5.302 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 O -0.337 6.337 2 C -0.190 4.190 3 Si 0.756 3.244 4 H -0.151 1.151 5 H -0.175 1.175 6 H -0.174 1.174 7 C -0.158 4.158 8 C 0.097 3.903 9 F -0.101 7.101 10 C -0.171 4.171 11 C -0.075 4.075 12 C -0.147 4.147 13 C 0.347 3.653 14 O -0.405 6.405 15 N -0.375 5.375 16 C 0.031 3.969 17 C -0.155 4.155 18 C -0.492 4.492 19 Si 1.191 2.809 20 C -0.491 4.491 21 C -0.146 4.146 22 C -0.170 4.170 23 C -0.147 4.147 24 C -0.494 4.494 25 C -0.489 4.489 26 C -0.154 4.154 27 C -0.029 4.029 28 H 0.223 0.777 29 H 0.142 0.858 30 H 0.165 0.835 31 H 0.156 0.844 32 H 0.234 0.766 33 H 0.111 0.889 34 H 0.084 0.916 35 H 0.087 0.913 36 H 0.098 0.902 37 H 0.097 0.903 38 H 0.092 0.908 39 H 0.094 0.906 40 H 0.080 0.920 41 H 0.081 0.919 42 H 0.083 0.917 43 H 0.080 0.920 44 H 0.080 0.920 45 H 0.081 0.919 46 H 0.094 0.906 47 H 0.093 0.907 48 H 0.098 0.902 49 H 0.101 0.899 50 H 0.090 0.910 51 H 0.090 0.910 52 H 0.169 0.831 Dipole moment (debyes) X Y Z Total from point charges 1.942 -0.226 -3.729 4.211 hybrid contribution -0.273 0.215 -1.063 1.118 sum 1.669 -0.012 -4.792 5.074 Atomic orbital electron populations 1.86572 1.16566 1.81470 1.49054 1.24991 0.85487 1.11044 0.97479 0.86135 0.80747 0.75513 0.81967 1.15134 1.17528 1.17425 1.20338 1.04213 0.96614 0.94618 1.17916 0.97348 0.80630 0.94419 1.91594 1.84465 1.37550 1.96539 1.21107 1.01788 0.95956 0.98291 1.21197 0.95648 0.93781 0.96846 1.20051 1.04781 0.95463 0.94438 1.18106 0.77596 0.85202 0.84345 1.90820 1.55694 1.66912 1.27056 1.46056 1.63508 1.17736 1.10181 1.21135 0.94738 0.82009 0.98994 1.21619 0.94492 0.98187 1.01228 1.26216 1.04435 1.06740 1.11792 0.79483 0.67295 0.66781 0.67314 1.26249 1.12094 0.99868 1.10882 1.21238 1.01880 0.93687 0.97813 1.21506 0.98581 0.98099 0.98825 1.21234 0.97819 0.99258 0.96346 1.26251 1.04596 1.16212 1.02362 1.26178 1.17889 0.97247 1.07582 1.21520 1.02867 0.95298 0.95750 1.21061 0.91271 0.91347 0.99203 0.77701 0.85756 0.83498 0.84365 0.76592 0.88905 0.91569 0.91287 0.90185 0.90279 0.90827 0.90597 0.92038 0.91905 0.91671 0.92019 0.92018 0.91873 0.90568 0.90685 0.90163 0.89943 0.90994 0.91049 0.83130 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 12. 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 O -0.53 -0.54 12.88 -35.23 -0.45 -1.00 16 2 C -0.08 -0.01 2.80 36.01 0.10 0.09 16 3 Si 0.94 0.26 26.09 -169.99 -4.43 -4.18 16 4 H -0.23 -0.14 5.28 56.52 0.30 0.16 16 5 H -0.25 -0.39 7.11 56.52 0.40 0.01 16 6 H -0.25 -0.28 7.11 56.52 0.40 0.12 16 7 C -0.16 -0.06 4.78 -104.61 -0.50 -0.56 16 8 C 0.12 -0.01 6.63 -39.27 -0.26 -0.27 16 9 F -0.12 -0.05 12.57 2.25 0.03 -0.02 16 10 C -0.15 0.26 10.01 -39.58 -0.40 -0.13 16 11 C -0.06 0.12 9.54 -39.20 -0.37 -0.25 16 12 C -0.14 -0.04 5.87 -104.96 -0.62 -0.65 16 13 C 0.56 0.95 7.56 -12.39 -0.09 0.86 16 14 O -0.53 -2.86 15.99 5.30 0.08 -2.78 16 15 N -0.72 0.78 4.98 -55.23 -0.28 0.50 16 16 C 0.14 -0.10 2.02 -68.38 -0.14 -0.24 16 17 C -0.12 0.06 5.20 -26.43 -0.14 -0.08 16 18 C -0.41 0.40 4.97 37.69 0.19 0.59 16 19 Si 1.00 -2.58 9.18 -169.99 -1.56 -4.15 16 20 C -0.40 0.79 4.97 37.67 0.19 0.98 16 21 C -0.11 0.15 5.64 -26.46 -0.15 0.00 16 22 C -0.13 0.14 5.40 -27.20 -0.15 -0.01 16 23 C -0.11 0.11 5.63 -26.46 -0.15 -0.04 16 24 C -0.41 0.53 4.97 37.65 0.19 0.72 16 25 C -0.40 0.95 4.97 37.66 0.19 1.14 16 26 C -0.12 0.29 5.08 -26.45 -0.13 0.16 16 27 C -0.01 -0.01 9.21 -39.20 -0.36 -0.37 16 28 H 0.38 -0.69 8.41 45.56 0.38 -0.31 16 29 H 0.12 -0.12 7.77 -51.92 -0.40 -0.53 16 30 H 0.15 -0.55 8.06 -52.48 -0.42 -0.98 16 31 H 0.14 -0.73 6.40 -52.49 -0.34 -1.07 16 32 H 0.40 -1.59 6.62 -40.82 -0.27 -1.86 16 33 H 0.09 0.11 7.29 -51.93 -0.38 -0.27 16 34 H 0.07 -0.13 7.71 -51.93 -0.40 -0.53 16 35 H 0.07 0.11 7.95 -51.93 -0.41 -0.30 16 36 H 0.08 -0.01 7.98 -51.93 -0.41 -0.42 16 37 H 0.08 -0.09 7.98 -51.93 -0.41 -0.51 16 38 H 0.07 -0.17 7.98 -51.93 -0.41 -0.59 16 39 H 0.08 -0.12 7.98 -51.93 -0.41 -0.54 16 40 H 0.06 -0.07 8.14 -51.93 -0.42 -0.49 16 41 H 0.06 -0.12 7.72 -51.93 -0.40 -0.52 16 42 H 0.06 -0.06 7.89 -51.92 -0.41 -0.47 16 43 H 0.06 -0.07 8.14 -51.93 -0.42 -0.49 16 44 H 0.06 -0.04 8.14 -51.93 -0.42 -0.46 16 45 H 0.06 -0.10 7.72 -51.93 -0.40 -0.50 16 46 H 0.08 -0.08 7.98 -51.93 -0.41 -0.50 16 47 H 0.07 -0.09 7.98 -51.93 -0.41 -0.51 16 48 H 0.08 -0.16 7.98 -51.93 -0.41 -0.57 16 49 H 0.08 -0.22 7.98 -51.93 -0.41 -0.64 16 50 H 0.07 -0.19 8.13 -51.93 -0.42 -0.61 16 51 H 0.07 -0.26 7.71 -51.93 -0.40 -0.66 16 52 H 0.15 0.14 7.58 -52.49 -0.40 -0.26 16 LS Contribution 401.68 15.07 6.05 6.05 Total: 0.00 -6.60 401.68 -11.32 -17.92 By element: Atomic # 1 Polarization: -6.13 SS G_CDS: -8.16 Total: -14.29 kcal Atomic # 6 Polarization: 4.54 SS G_CDS: -2.61 Total: 1.93 kcal Atomic # 7 Polarization: 0.78 SS G_CDS: -0.28 Total: 0.50 kcal Atomic # 8 Polarization: -3.41 SS G_CDS: -0.37 Total: -3.78 kcal Atomic # 9 Polarization: -0.05 SS G_CDS: 0.03 Total: -0.02 kcal Atomic # 14 Polarization: -2.33 SS G_CDS: -6.00 Total: -8.32 kcal Total LS contribution 6.05 Total: 6.05 kcal Total: -6.60 -11.32 -17.92 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. m_22_15739362_12496548.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 -121.184 kcal (2) G-P(sol) polarization free energy of solvation -6.600 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -127.783 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -11.320 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -17.920 kcal (6) G-S(sol) free energy of system = (1) + (5) -139.103 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.78 seconds