Wall clock time and date at job start Tue Mar 31 2020 06:11:52 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 * SM5.42R - SM5.42R CALCULATIONS WILL BE PERFORMED * - CHARGE MODEL 2 WILL BE USED (CM2A) * - THE SOLVENT IS WATER ******************************************************************************* 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 C 1.53003 * 109.46597 * 2 1 4 4 C 1.53006 * 109.47054 * 179.97438 * 3 2 1 5 5 N 1.46496 * 109.47238 * 179.97438 * 4 3 2 6 6 N 1.28387 * 119.10522 * 270.25893 * 5 4 3 7 7 C 1.31400 * 122.26740 * 180.27644 * 6 5 4 8 8 C 1.47887 * 119.89378 * 179.74512 * 7 6 5 9 9 O 1.21551 * 120.00163 * 158.33982 * 8 7 6 10 10 N 1.34775 * 119.99940 * 338.34113 * 8 7 6 11 11 C 1.46506 * 119.99992 * 178.22569 * 10 8 7 12 12 C 1.50700 * 109.47084 * 269.99882 * 11 10 8 13 13 H 1.08003 * 120.00105 * 359.97438 * 12 11 10 14 14 C 1.37878 * 119.99879 * 179.97438 * 12 11 10 15 15 N 1.31511 * 119.99894 * 179.97438 * 14 12 11 16 16 Si 1.86801 * 120.00118 * 0.02562 * 14 12 11 17 17 H 1.48498 * 109.46998 * 89.99754 * 16 14 12 18 18 H 1.48499 * 109.47140 * 209.99906 * 16 14 12 19 19 H 1.48508 * 109.46794 * 329.99651 * 16 14 12 20 20 C 1.46501 * 120.00179 * 358.22449 * 10 8 7 21 21 C 1.52997 * 109.47288 * 225.50014 * 20 10 8 22 22 C 1.53002 * 109.46952 * 65.05270 * 21 20 10 23 23 C 1.52999 * 109.46859 * 287.69130 * 22 21 20 24 24 C 1.53001 * 109.47013 * 185.04594 * 21 20 10 25 25 C 1.53005 * 109.46965 * 65.00601 * 24 21 20 26 26 C 1.41573 * 120.21521 * 359.47090 * 7 6 5 27 27 C 1.36328 * 118.23589 * 0.25127 * 26 7 6 28 28 C 1.34615 * 119.09690 * 90.00524 * 5 4 3 29 29 O 1.21972 * 120.22706 * 0.02562 * 28 5 4 30 30 H 1.09005 * 109.47016 * 179.97438 * 1 2 3 31 31 H 1.09004 * 109.47273 * 299.99490 * 1 2 3 32 32 H 1.08993 * 109.47135 * 60.00149 * 1 2 3 33 33 H 1.08995 * 109.47525 * 240.00397 * 2 1 3 34 34 H 1.08993 * 109.47135 * 119.99851 * 2 1 3 35 35 H 1.09004 * 109.46893 * 60.00277 * 3 2 1 36 36 H 1.08999 * 109.47661 * 299.99727 * 3 2 1 37 37 H 1.09004 * 109.46711 * 59.99952 * 4 3 2 38 38 H 1.08999 * 109.46857 * 300.00456 * 4 3 2 39 39 H 1.08998 * 109.46743 * 30.00026 * 11 10 8 40 40 H 1.09003 * 109.46919 * 149.99830 * 11 10 8 41 41 H 0.97000 * 119.99884 * 179.97438 * 15 14 12 42 42 H 1.09000 * 109.47043 * 345.50411 * 20 10 8 43 43 H 1.09001 * 109.46852 * 105.50048 * 20 10 8 44 44 H 1.08994 * 109.47148 * 305.04940 * 21 20 10 45 45 H 1.08997 * 109.47093 * 47.69066 * 22 21 20 46 46 H 1.09002 * 109.47037 * 167.69315 * 22 21 20 47 47 H 1.09003 * 109.47373 * 71.14024 * 23 22 21 48 48 H 1.09004 * 109.47235 * 191.13913 * 23 22 21 49 49 H 1.09001 * 109.47136 * 311.13972 * 23 22 21 50 50 H 1.09000 * 109.47205 * 184.99726 * 24 21 20 51 51 H 1.09003 * 109.47524 * 305.00631 * 24 21 20 52 52 H 1.09002 * 109.47122 * 179.97438 * 25 24 21 53 53 H 1.09005 * 109.47100 * 299.99830 * 25 24 21 54 54 H 1.09003 * 109.47351 * 59.99304 * 25 24 21 55 55 H 1.07995 * 120.88552 * 180.28027 * 26 7 6 56 56 H 1.08005 * 121.03655 * 180.02562 * 27 26 7 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.0399 1.4426 0.0000 4 6 3.5700 1.4427 0.0006 5 7 4.0582 2.8239 0.0000 6 7 4.2620 3.4137 -1.1220 7 6 4.7041 4.6490 -1.1950 8 6 4.9258 5.2755 -2.5161 9 8 4.9464 6.4865 -2.6187 10 7 5.1050 4.5038 -3.6064 11 6 5.2858 5.1241 -4.9213 12 6 3.9427 5.3048 -5.5805 13 1 3.0443 4.9888 -5.0711 14 6 3.8609 5.8736 -6.8338 15 7 2.6888 6.0308 -7.4092 16 14 5.4148 6.4195 -7.7151 17 1 5.9414 5.2963 -8.5314 18 1 5.1037 7.5710 -8.5996 19 1 6.4346 6.8252 -6.7146 20 6 5.1189 3.0446 -3.4768 21 6 6.2909 2.4726 -4.2769 22 6 7.6087 2.9453 -3.6597 23 6 7.8095 4.4315 -3.9626 24 6 6.2323 0.9441 -4.2452 25 6 4.9741 0.4667 -4.9732 26 6 4.9632 5.3840 -0.0132 27 6 4.7468 4.7737 1.1865 28 6 4.2760 3.4417 1.1760 29 8 4.0687 2.8555 2.2253 30 1 -0.3633 -1.0277 -0.0005 31 1 -0.3634 0.5138 0.8900 32 1 -0.3633 0.5138 -0.8899 33 1 1.8934 -0.5137 -0.8900 34 1 1.8933 -0.5138 0.8899 35 1 1.6765 1.9563 0.8900 36 1 1.6766 1.9564 -0.8900 37 1 3.9337 0.9287 -0.8892 38 1 3.9329 0.9289 0.8908 39 1 5.7656 6.0955 -4.8018 40 1 5.9120 4.4834 -5.5422 41 1 2.6312 6.4313 -8.2908 42 1 5.2295 2.7742 -2.4267 43 1 4.1837 2.6368 -3.8604 44 1 6.2288 2.8173 -5.3090 45 1 7.5789 2.7953 -2.5805 46 1 8.4342 2.3734 -4.0835 47 1 7.1019 5.0202 -3.3788 48 1 8.8270 4.7216 -3.7002 49 1 7.6428 4.6116 -5.0246 50 1 7.1144 0.5364 -4.7390 51 1 6.2048 0.6028 -3.2103 52 1 4.9326 -0.6223 -4.9510 53 1 4.0919 0.8745 -4.4795 54 1 5.0015 0.8080 -6.0081 55 1 5.3224 6.4014 -0.0606 56 1 4.9292 5.2923 2.1161 RHF calculation, no. of doubly occupied orbitals= 70 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=WATER ZINC000075583916.mol2 56 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 06:11:52 Heat of formation + Delta-G solvation = -22.428916 kcal Electronic energy + Delta-G solvation = -35659.849202 eV Core-core repulsion = 31377.180734 eV Total energy + Delta-G solvation = -4282.668468 eV No. of doubly occupied orbitals = 70 Molecular weight (most abundant/longest-lived isotopes) = 363.231 amu Computer time = 1.23 seconds Orbital eigenvalues (eV) -43.37830 -41.17411 -38.93378 -38.62799 -37.19639 -35.71365 -35.07959 -34.54551 -33.07098 -32.60685 -30.88788 -30.04796 -28.66516 -27.01984 -25.99167 -25.58676 -24.74847 -23.78501 -22.75562 -22.18057 -21.64554 -20.51876 -19.58134 -19.06208 -18.80240 -17.94902 -17.72085 -17.19254 -17.03529 -16.82176 -16.39523 -15.99902 -15.78339 -15.27638 -15.02162 -14.84583 -14.79494 -14.66733 -14.53905 -14.32474 -14.15372 -14.09543 -13.88144 -13.62047 -13.55636 -13.45395 -13.36625 -13.01321 -12.95136 -12.56424 -12.52416 -12.38714 -12.34693 -12.23872 -12.10132 -12.01127 -11.79114 -11.66901 -11.54778 -11.49266 -11.40824 -11.31267 -11.20748 -11.11758 -10.76963 -10.60408 -9.70789 -9.29092 -8.30251 -6.34164 -0.69389 0.01137 0.97162 1.33070 1.37467 1.47386 1.71209 2.08547 2.33159 2.36020 2.80225 3.04268 3.46963 3.54051 3.68896 3.75596 3.82459 3.87609 3.91176 4.02906 4.15489 4.31901 4.39183 4.43157 4.44342 4.48298 4.49057 4.60349 4.64953 4.65865 4.68866 4.72491 4.80070 4.84771 4.87696 4.89942 4.90773 4.96619 5.02229 5.04023 5.14247 5.18179 5.19616 5.21024 5.26624 5.28825 5.29276 5.40892 5.47476 5.50517 5.54294 5.60213 5.63735 6.24949 6.33865 6.55540 6.67371 6.80004 7.13898 7.31367 8.86780 Molecular weight = 363.23amu Principal moments of inertia in cm(-1) A = 0.008966 B = 0.003585 C = 0.003363 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 3122.285364 B = 7808.270361 C = 8323.607654 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C -0.155 4.155 2 C -0.113 4.113 3 C -0.131 4.131 4 C 0.103 3.897 5 N -0.335 5.335 6 N -0.168 5.168 7 C 0.014 3.986 8 C 0.569 3.431 9 O -0.569 6.569 10 N -0.564 5.564 11 C 0.230 3.770 12 C -0.533 4.533 13 H 0.027 0.973 14 C 0.017 3.983 15 N -0.945 5.945 16 Si 0.934 3.066 17 H -0.253 1.253 18 H -0.284 1.284 19 H -0.248 1.248 20 C 0.107 3.893 21 C -0.113 4.113 22 C -0.103 4.103 23 C -0.168 4.168 24 C -0.102 4.102 25 C -0.157 4.157 26 C -0.028 4.028 27 C -0.184 4.184 28 C 0.473 3.527 29 O -0.554 6.554 30 H 0.069 0.931 31 H 0.049 0.951 32 H 0.047 0.953 33 H 0.080 0.920 34 H 0.077 0.923 35 H 0.054 0.946 36 H 0.058 0.942 37 H 0.124 0.876 38 H 0.117 0.883 39 H 0.026 0.974 40 H 0.058 0.942 41 H 0.267 0.733 42 H 0.117 0.883 43 H 0.058 0.942 44 H 0.062 0.938 45 H 0.064 0.936 46 H 0.079 0.921 47 H 0.045 0.955 48 H 0.059 0.941 49 H 0.044 0.956 50 H 0.077 0.923 51 H 0.077 0.923 52 H 0.070 0.930 53 H 0.048 0.952 54 H 0.039 0.961 55 H 0.174 0.826 56 H 0.194 0.806 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 7.928 -14.609 14.442 22.020 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.212 4.212 2 C -0.151 4.151 3 C -0.169 4.169 4 C -0.015 4.015 5 N -0.156 5.156 6 N 0.001 4.999 7 C -0.156 4.156 8 C 0.354 3.646 9 O -0.447 6.447 10 N -0.293 5.293 11 C 0.108 3.892 12 C -0.572 4.572 13 H 0.046 0.954 14 C -0.181 4.181 15 N -0.585 5.585 16 Si 0.757 3.243 17 H -0.177 1.177 18 H -0.209 1.209 19 H -0.172 1.172 20 C -0.016 4.016 21 C -0.132 4.132 22 C -0.140 4.140 23 C -0.225 4.225 24 C -0.140 4.140 25 C -0.214 4.214 26 C -0.051 4.051 27 C -0.214 4.214 28 C 0.273 3.727 29 O -0.435 6.435 30 H 0.088 0.912 31 H 0.068 0.932 32 H 0.066 0.934 33 H 0.099 0.901 34 H 0.096 0.904 35 H 0.073 0.927 36 H 0.077 0.923 37 H 0.142 0.858 38 H 0.135 0.865 39 H 0.044 0.956 40 H 0.075 0.925 41 H 0.076 0.924 42 H 0.134 0.866 43 H 0.076 0.924 44 H 0.080 0.920 45 H 0.082 0.918 46 H 0.098 0.902 47 H 0.064 0.936 48 H 0.078 0.922 49 H 0.063 0.937 50 H 0.096 0.904 51 H 0.096 0.904 52 H 0.089 0.911 53 H 0.067 0.933 54 H 0.058 0.942 55 H 0.191 0.809 56 H 0.212 0.788 Dipole moment (debyes) X Y Z Total from point charges 7.472 -15.006 14.319 22.046 hybrid contribution -0.151 1.542 -0.421 1.606 sum 7.321 -13.463 13.898 20.689 Atomic orbital electron populations 1.21776 0.95591 1.02410 1.01454 1.21402 0.96180 0.93908 1.03617 1.21919 0.93262 0.99780 1.01959 1.22446 0.95700 0.76370 1.07002 1.47285 1.50376 1.15680 1.02261 1.71517 1.06735 1.10101 1.11547 1.21701 1.07107 0.90887 0.95904 1.17621 0.77895 0.86030 0.83005 1.90770 1.53828 1.14594 1.85508 1.47824 1.66835 1.07799 1.06798 1.19564 0.95326 0.93942 0.80341 1.21562 0.93226 1.40523 1.01841 0.95419 1.25731 0.99409 0.95028 0.97974 1.70047 1.26130 1.45828 1.16505 0.85907 0.80100 0.79064 0.79248 1.17712 1.20937 1.17208 1.22426 0.98690 0.77779 1.02676 1.21708 0.94639 0.97039 0.99786 1.21466 0.96670 0.95304 1.00598 1.22226 1.01255 0.97730 1.01287 1.21429 0.98705 0.92195 1.01685 1.21849 0.98064 1.01460 1.00028 1.22292 0.91412 1.01775 0.89619 1.22895 1.03936 0.93606 1.00930 1.18388 0.80823 0.89540 0.83926 1.90936 1.54653 1.67783 1.30117 0.91203 0.93164 0.93430 0.90125 0.90409 0.92730 0.92343 0.85759 0.86515 0.95553 0.92465 0.92412 0.86562 0.92397 0.92004 0.91777 0.90199 0.93632 0.92237 0.93713 0.90434 0.90398 0.91065 0.93253 0.94210 0.80886 0.78844 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 15. 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 -2.29 9.91 71.98 0.71 -1.58 16 2 C -0.11 -1.56 5.93 30.60 0.18 -1.38 16 3 C -0.13 -2.44 5.54 30.60 0.17 -2.27 16 4 C 0.10 1.84 5.55 86.38 0.48 2.32 16 5 N -0.34 -7.96 3.22 -472.47 -1.52 -9.49 16 6 N -0.17 -4.56 7.03 -48.20 -0.34 -4.90 16 7 C 0.01 0.41 6.45 42.59 0.27 0.69 16 8 C 0.57 22.45 6.31 86.80 0.55 23.00 16 9 O -0.57 -26.59 16.15 -3.88 -0.06 -26.65 16 10 N -0.56 -22.09 1.33 -829.51 -1.11 -23.19 16 11 C 0.23 10.63 1.89 85.63 0.16 10.79 16 12 C -0.53 -30.53 9.94 25.60 0.25 -30.28 16 13 H 0.03 1.79 8.06 -2.91 -0.02 1.77 16 14 C 0.02 0.96 5.47 84.66 0.46 1.43 16 15 N -0.94 -58.96 13.96 21.45 0.30 -58.66 16 16 Si 0.93 43.14 27.47 68.60 1.88 45.02 16 17 H -0.25 -10.75 7.11 99.48 0.71 -10.04 16 18 H -0.28 -13.07 7.11 99.48 0.71 -12.36 16 19 H -0.25 -10.86 6.54 99.48 0.65 -10.21 16 20 C 0.11 3.26 4.49 86.38 0.39 3.65 16 21 C -0.11 -2.93 1.86 -10.79 -0.02 -2.95 16 22 C -0.10 -2.38 5.40 30.59 0.17 -2.21 16 23 C -0.17 -4.81 8.18 71.98 0.59 -4.22 16 24 C -0.10 -1.96 5.38 30.60 0.16 -1.80 16 25 C -0.16 -3.34 9.67 71.98 0.70 -2.64 16 26 C -0.03 -0.68 9.68 22.65 0.22 -0.46 16 27 C -0.18 -4.00 10.11 22.57 0.23 -3.77 16 28 C 0.47 12.07 7.47 85.01 0.63 12.71 16 29 O -0.55 -16.33 16.62 -5.22 -0.09 -16.42 16 30 H 0.07 0.85 8.14 -2.38 -0.02 0.83 16 31 H 0.05 0.77 8.14 -2.39 -0.02 0.75 16 32 H 0.05 0.79 8.14 -2.39 -0.02 0.77 16 33 H 0.08 0.94 8.14 -2.39 -0.02 0.92 16 34 H 0.08 0.91 8.14 -2.39 -0.02 0.89 16 35 H 0.05 1.19 8.14 -2.39 -0.02 1.17 16 36 H 0.06 1.29 8.14 -2.39 -0.02 1.27 16 37 H 0.12 1.69 8.14 -2.38 -0.02 1.67 16 38 H 0.12 1.90 7.36 -2.39 -0.02 1.89 16 39 H 0.03 1.32 6.30 -2.39 -0.02 1.31 16 40 H 0.06 2.46 3.65 -2.39 -0.01 2.45 16 41 H 0.27 15.54 8.31 -92.71 -0.77 14.77 16 42 H 0.12 3.07 4.44 -27.75 -0.12 2.95 16 43 H 0.06 1.95 6.34 -2.39 -0.02 1.94 16 44 H 0.06 1.94 6.02 -2.39 -0.01 1.93 16 45 H 0.06 1.35 8.05 -2.39 -0.02 1.33 16 46 H 0.08 1.52 8.05 -2.39 -0.02 1.50 16 47 H 0.04 1.52 4.51 -2.39 -0.01 1.51 16 48 H 0.06 1.47 8.14 -2.38 -0.02 1.45 16 49 H 0.04 1.41 6.13 -2.39 -0.01 1.39 16 50 H 0.08 1.28 8.09 -2.39 -0.02 1.26 16 51 H 0.08 1.21 8.14 -2.39 -0.02 1.19 16 52 H 0.07 1.19 8.14 -2.39 -0.02 1.17 16 53 H 0.05 1.15 6.46 -2.38 -0.02 1.14 16 54 H 0.04 0.98 8.14 -2.39 -0.02 0.96 16 55 H 0.17 3.78 7.86 -2.91 -0.02 3.76 16 56 H 0.19 3.03 8.06 -2.91 -0.02 3.01 16 Total: -1.00 -75.01 431.17 6.08 -68.93 By element: Atomic # 1 Polarization: 23.62 SS G_CDS: 0.70 Total: 24.32 kcal Atomic # 6 Polarization: -5.28 SS G_CDS: 6.31 Total: 1.03 kcal Atomic # 7 Polarization: -93.57 SS G_CDS: -2.67 Total: -96.24 kcal Atomic # 8 Polarization: -42.92 SS G_CDS: -0.15 Total: -43.07 kcal Atomic # 14 Polarization: 43.14 SS G_CDS: 1.88 Total: 45.02 kcal Total: -75.01 6.08 -68.93 kcal The number of atoms in the molecule is 56 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. ZINC000075583916.mol2 56 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 46.505 kcal (2) G-P(sol) polarization free energy of solvation -75.012 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -28.507 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy 6.078 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -68.934 kcal (6) G-S(sol) free energy of system = (1) + (5) -22.429 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.23 seconds