Wall clock time and date at job start Sat Apr 11 2020 06:10: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 O 1.45193 * 1 3 3 C 1.34225 * 116.99940 * 2 1 4 4 O 1.20826 * 120.00300 * 0.02562 * 3 2 1 5 5 C 1.50706 * 120.00086 * 179.97438 * 3 2 1 6 6 C 1.52998 * 109.47064 * 59.99515 * 5 3 2 7 7 C 1.53006 * 109.46976 * 299.99861 * 5 3 2 8 8 C 1.52998 * 109.47049 * 180.02562 * 5 3 2 9 9 H 1.08999 * 109.49555 * 299.93294 * 8 5 3 10 10 C 1.53032 * 109.50248 * 59.99932 * 8 5 3 11 11 C 1.53039 * 109.54002 * 181.34837 * 10 8 5 12 12 C 1.53187 * 109.31301 * 298.63530 * 11 10 8 13 13 N 1.46928 * 108.77730 * 54.63266 * 12 11 10 14 14 C 1.34776 * 120.62874 * 126.36878 * 13 12 11 15 15 O 1.21588 * 119.99890 * 185.22993 * 14 13 12 16 16 N 1.34768 * 120.00360 * 5.23095 * 14 13 12 17 17 C 1.46504 * 120.00353 * 175.29509 * 16 14 13 18 18 C 1.52993 * 109.47321 * 179.72443 * 17 16 14 19 19 C 1.53000 * 109.47432 * 180.02562 * 18 17 16 20 20 C 1.50697 * 109.47698 * 179.97438 * 19 18 17 21 21 H 1.08002 * 119.99560 * 359.97438 * 20 19 18 22 22 C 1.37877 * 120.00371 * 179.97438 * 20 19 18 23 23 N 1.31508 * 120.00244 * 179.97438 * 22 20 19 24 24 Si 1.86805 * 119.99686 * 0.02562 * 22 20 19 25 25 H 1.48492 * 109.47509 * 89.99537 * 24 22 20 26 26 H 1.48507 * 109.46724 * 209.99898 * 24 22 20 27 27 H 1.48497 * 109.47013 * 329.99169 * 24 22 20 28 28 C 1.46926 * 118.73872 * 306.64722 * 13 12 11 29 29 H 1.09004 * 109.47098 * 60.00016 * 1 2 3 30 30 H 1.09006 * 109.47346 * 180.02562 * 1 2 3 31 31 H 1.09003 * 109.47624 * 300.00198 * 1 2 3 32 32 H 1.09000 * 109.46861 * 179.97438 * 6 5 3 33 33 H 1.08992 * 109.47470 * 300.00851 * 6 5 3 34 34 H 1.09003 * 109.46501 * 60.00367 * 6 5 3 35 35 H 1.08996 * 109.47114 * 300.00191 * 7 5 3 36 36 H 1.09002 * 109.47131 * 59.99708 * 7 5 3 37 37 H 1.08997 * 109.47057 * 179.97438 * 7 5 3 38 38 H 1.09004 * 109.46152 * 61.32724 * 10 8 5 39 39 H 1.09003 * 109.45836 * 301.37137 * 10 8 5 40 40 H 1.09004 * 109.49826 * 178.68011 * 11 10 8 41 41 H 1.08999 * 109.52273 * 58.60011 * 11 10 8 42 42 H 1.09006 * 109.58876 * 294.84336 * 12 11 10 43 43 H 1.08998 * 109.59053 * 174.42238 * 12 11 10 44 44 H 0.97009 * 119.99885 * 355.01222 * 16 14 13 45 45 H 1.09003 * 109.46867 * 299.72054 * 17 16 14 46 46 H 1.09000 * 109.46688 * 59.72276 * 17 16 14 47 47 H 1.09008 * 109.47523 * 299.99675 * 18 17 16 48 48 H 1.08997 * 109.47786 * 59.99611 * 18 17 16 49 49 H 1.09008 * 109.46908 * 300.00025 * 19 18 17 50 50 H 1.08998 * 109.47496 * 59.99296 * 19 18 17 51 51 H 0.96997 * 120.00174 * 180.02562 * 23 22 20 52 52 H 1.09000 * 109.58564 * 173.14188 * 28 13 12 53 53 H 1.09004 * 109.58305 * 293.56850 * 28 13 12 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 8 1.4519 0.0000 0.0000 3 6 2.0613 1.1960 0.0000 4 8 1.4033 2.2093 0.0005 5 6 3.5663 1.2749 0.0006 6 6 4.1138 0.5810 -1.2482 7 6 4.1128 0.5816 1.2503 8 6 4.0000 2.7421 0.0011 9 1 3.6118 3.2365 -0.8894 10 6 3.4527 3.4362 1.2503 11 6 3.8577 4.9119 1.2348 12 6 5.3858 5.0181 1.2531 13 7 5.9279 4.2318 0.1365 14 6 6.7739 4.7919 -0.7507 15 8 7.1408 4.1605 -1.7228 16 7 7.2131 6.0517 -0.5604 17 6 8.2113 6.6243 -1.4670 18 6 8.5291 8.0575 -1.0362 19 6 9.5711 8.6557 -1.9834 20 6 9.8837 10.0676 -1.5594 21 1 9.3882 10.4957 -0.7005 22 6 10.8021 10.8130 -2.2679 23 7 11.0744 12.0453 -1.8982 24 14 11.6584 10.0727 -3.7540 25 1 10.8610 10.3522 -4.9751 26 1 13.0093 10.6730 -3.8967 27 1 11.7860 8.6046 -3.5713 28 6 5.5254 2.8252 0.0021 29 1 -0.3633 0.5138 -0.8900 30 1 -0.3634 -1.0277 0.0005 31 1 -0.3634 0.5139 0.8900 32 1 5.2023 0.6385 -1.2480 33 1 3.8047 -0.4642 -1.2482 34 1 3.7244 1.0749 -2.1385 35 1 3.7227 1.0758 2.1399 36 1 3.8037 -0.4637 1.2504 37 1 5.2013 0.6382 1.2505 38 1 3.8608 2.9570 2.1403 39 1 2.3655 3.3581 1.2622 40 1 3.4478 5.4110 2.1130 41 1 3.4709 5.3864 0.3330 42 1 5.7682 4.6261 2.1956 43 1 5.6810 6.0615 1.1424 44 1 6.8652 6.5808 0.1744 45 1 9.1207 6.0243 -1.4318 46 1 7.8183 6.6291 -2.4837 47 1 7.6198 8.6576 -1.0714 48 1 8.9220 8.0527 -0.0195 49 1 10.4806 8.0557 -1.9486 50 1 9.1779 8.6606 -2.9999 51 1 11.7203 12.5697 -2.3969 52 1 5.9126 2.4222 -0.9338 53 1 5.9193 2.2507 0.8405 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 ZINC000091743137.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:10:33 Heat of formation + Delta-G solvation = -127.680654 kcal Electronic energy + Delta-G solvation = -30585.871353 eV Core-core repulsion = 26457.372446 eV Total energy + Delta-G solvation = -4128.498907 eV No. of doubly occupied orbitals = 66 Molecular weight (most abundant/longest-lived isotopes) = 340.220 amu Computer time = 0.64 seconds Orbital eigenvalues (eV) -41.53890 -40.31750 -38.50943 -36.78986 -35.78949 -35.49187 -33.31428 -32.54855 -31.41421 -30.18343 -28.15980 -27.98832 -27.25892 -25.44145 -24.77652 -23.09481 -22.58046 -21.38726 -20.55295 -19.78558 -19.18797 -18.15449 -17.68000 -17.00442 -16.78097 -16.73672 -16.10952 -15.68680 -14.98959 -14.92010 -14.72686 -14.52699 -14.42780 -14.14731 -14.08802 -13.71503 -13.47685 -13.35657 -13.18538 -12.94111 -12.82913 -12.66912 -12.40808 -12.28486 -12.06051 -11.94092 -11.88655 -11.73276 -11.47266 -11.26224 -11.06509 -11.02666 -10.87786 -10.85689 -10.62028 -10.35560 -10.19198 -9.97349 -9.80116 -9.65476 -9.21202 -8.97601 -8.54322 -8.49849 -6.05032 -4.03611 1.79543 2.22120 2.67755 3.28588 3.32792 3.40144 3.44834 3.88838 4.07091 4.23081 4.23784 4.45869 4.52053 4.63116 4.66117 4.72414 4.76096 4.77983 4.83773 5.01426 5.06254 5.15752 5.20669 5.28688 5.29237 5.32138 5.33875 5.45797 5.48567 5.54751 5.57571 5.62804 5.67298 5.71821 5.78163 5.82316 5.85233 6.06212 6.07890 6.15965 6.34357 6.44032 6.49183 6.72242 6.84853 6.97327 7.25311 7.34778 7.38654 7.77109 8.03224 8.51554 8.69086 8.98787 9.57342 11.08183 Molecular weight = 340.22amu Principal moments of inertia in cm(-1) A = 0.027184 B = 0.001862 C = 0.001817 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1029.777394 B =15031.627593 C =15407.324412 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 C 0.034 3.966 2 O -0.357 6.357 3 C 0.470 3.530 4 O -0.498 6.498 5 C -0.045 4.045 6 C -0.147 4.147 7 C -0.148 4.148 8 C -0.083 4.083 9 H 0.084 0.916 10 C -0.128 4.128 11 C -0.132 4.132 12 C 0.110 3.890 13 N -0.631 5.631 14 C 0.706 3.294 15 O -0.593 6.593 16 N -0.742 5.742 17 C 0.132 3.868 18 C -0.113 4.113 19 C -0.011 4.011 20 C -0.522 4.522 21 H 0.072 0.928 22 C 0.066 3.934 23 N -0.903 5.903 24 Si 0.880 3.120 25 H -0.277 1.277 26 H -0.288 1.288 27 H -0.264 1.264 28 C 0.134 3.866 29 H 0.063 0.937 30 H 0.101 0.899 31 H 0.062 0.938 32 H 0.079 0.921 33 H 0.071 0.929 34 H 0.063 0.937 35 H 0.065 0.935 36 H 0.070 0.930 37 H 0.078 0.922 38 H 0.063 0.937 39 H 0.087 0.913 40 H 0.076 0.924 41 H 0.068 0.932 42 H 0.069 0.931 43 H 0.079 0.921 44 H 0.394 0.606 45 H 0.061 0.939 46 H 0.062 0.938 47 H 0.059 0.941 48 H 0.060 0.940 49 H 0.019 0.981 50 H 0.015 0.985 51 H 0.263 0.737 52 H 0.098 0.902 53 H 0.068 0.932 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges -20.964 -27.532 9.300 35.833 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.060 4.060 2 O -0.273 6.273 3 C 0.313 3.687 4 O -0.385 6.385 5 C -0.046 4.046 6 C -0.204 4.204 7 C -0.205 4.205 8 C -0.103 4.103 9 H 0.102 0.898 10 C -0.165 4.165 11 C -0.170 4.170 12 C -0.014 4.014 13 N -0.369 5.369 14 C 0.407 3.593 15 O -0.471 6.471 16 N -0.398 5.398 17 C 0.009 3.991 18 C -0.151 4.151 19 C -0.049 4.049 20 C -0.560 4.560 21 H 0.090 0.910 22 C -0.135 4.135 23 N -0.542 5.542 24 Si 0.707 3.293 25 H -0.202 1.202 26 H -0.214 1.214 27 H -0.189 1.189 28 C 0.010 3.990 29 H 0.081 0.919 30 H 0.119 0.881 31 H 0.081 0.919 32 H 0.098 0.902 33 H 0.090 0.910 34 H 0.082 0.918 35 H 0.084 0.916 36 H 0.089 0.911 37 H 0.097 0.903 38 H 0.082 0.918 39 H 0.106 0.894 40 H 0.094 0.906 41 H 0.087 0.913 42 H 0.087 0.913 43 H 0.098 0.902 44 H 0.230 0.770 45 H 0.080 0.920 46 H 0.081 0.919 47 H 0.078 0.922 48 H 0.078 0.922 49 H 0.037 0.963 50 H 0.033 0.967 51 H 0.073 0.927 52 H 0.116 0.884 53 H 0.086 0.914 Dipole moment (debyes) X Y Z Total from point charges -20.759 -27.521 8.529 35.512 hybrid contribution -0.040 0.228 -0.406 0.467 sum -20.799 -27.293 8.123 35.263 Atomic orbital electron populations 1.23281 0.76384 1.03724 1.02602 1.86503 1.23312 1.33505 1.83978 1.23742 0.91777 0.81252 0.71911 1.90624 1.67105 1.36084 1.44710 1.20506 0.86805 0.97875 0.99447 1.21962 1.02569 1.00780 0.95135 1.22024 1.02364 1.00833 0.95299 1.21364 0.96111 0.92098 1.00711 0.89786 1.21880 1.02887 0.95037 0.96672 1.21882 0.95054 0.97138 1.02951 1.21670 0.94089 0.95923 0.89711 1.47197 1.48855 1.11986 1.28858 1.15878 0.79145 0.81958 0.82313 1.90893 1.62887 1.63253 1.30091 1.45511 1.46192 1.13306 1.34785 1.21518 0.89969 0.94941 0.92707 1.21704 1.00621 0.91904 1.00865 1.19500 0.93921 0.96797 0.94663 1.21099 1.23951 0.96132 1.14779 0.91029 1.25025 0.96378 0.94670 0.97406 1.69909 1.34643 1.09039 1.40627 0.86776 0.79022 0.80131 0.83384 1.20250 1.21388 1.18863 1.21470 0.92658 0.81253 1.03586 0.91863 0.88068 0.91896 0.90192 0.91024 0.91805 0.91601 0.91086 0.90320 0.91843 0.89422 0.90576 0.91300 0.91277 0.90235 0.76952 0.92019 0.91935 0.92201 0.92160 0.96273 0.96707 0.92721 0.88403 0.91379 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 C 0.03 0.18 10.55 101.05 1.07 1.25 16 2 O -0.36 -2.73 10.37 -39.24 -0.41 -3.14 16 3 C 0.47 4.16 6.26 36.01 0.23 4.38 16 4 O -0.50 -5.61 13.41 -18.75 -0.25 -5.86 16 5 C -0.04 -0.32 0.52 -155.66 -0.08 -0.40 16 6 C -0.15 -0.98 8.06 37.16 0.30 -0.68 16 7 C -0.15 -0.79 7.71 37.16 0.29 -0.51 16 8 C -0.08 -0.73 1.32 -90.73 -0.12 -0.85 16 9 H 0.08 0.93 8.10 -51.93 -0.42 0.51 16 10 C -0.13 -1.04 4.47 -26.69 -0.12 -1.15 16 11 C -0.13 -1.07 6.08 -26.61 -0.16 -1.23 16 12 C 0.11 0.94 6.39 -3.71 -0.02 0.92 16 13 N -0.63 -7.24 2.98 -178.76 -0.53 -7.77 16 14 C 0.71 10.85 8.30 -86.93 -0.72 10.13 16 15 O -0.59 -11.31 15.77 5.29 0.08 -11.23 16 16 N -0.74 -10.82 5.46 -65.23 -0.36 -11.18 16 17 C 0.13 2.26 5.68 -4.04 -0.02 2.24 16 18 C -0.11 -2.14 5.27 -26.73 -0.14 -2.28 16 19 C -0.01 -0.26 4.04 -27.88 -0.11 -0.37 16 20 C -0.52 -14.77 9.66 -34.44 -0.33 -15.11 16 21 H 0.07 2.16 7.99 -52.49 -0.42 1.74 16 22 C 0.07 1.94 5.47 -17.82 -0.10 1.84 16 23 N -0.90 -27.85 13.96 55.43 0.77 -27.07 16 24 Si 0.88 21.88 27.47 -169.99 -4.67 17.21 16 25 H -0.28 -6.78 7.11 56.52 0.40 -6.38 16 26 H -0.29 -7.09 7.11 56.52 0.40 -6.69 16 27 H -0.26 -6.33 6.54 56.52 0.37 -5.96 16 28 C 0.13 1.29 4.97 -3.93 -0.02 1.27 16 29 H 0.06 0.35 8.13 -51.93 -0.42 -0.07 16 30 H 0.10 0.27 8.14 -51.93 -0.42 -0.15 16 31 H 0.06 0.33 8.13 -51.93 -0.42 -0.09 16 32 H 0.08 0.50 6.57 -51.93 -0.34 0.16 16 33 H 0.07 0.40 8.11 -51.93 -0.42 -0.02 16 34 H 0.06 0.48 8.14 -51.93 -0.42 0.06 16 35 H 0.07 0.35 6.46 -51.93 -0.34 0.01 16 36 H 0.07 0.34 8.11 -51.93 -0.42 -0.08 16 37 H 0.08 0.36 6.32 -51.93 -0.33 0.04 16 38 H 0.06 0.43 6.50 -51.93 -0.34 0.10 16 39 H 0.09 0.84 5.49 -51.93 -0.29 0.55 16 40 H 0.08 0.50 8.14 -51.93 -0.42 0.08 16 41 H 0.07 0.66 8.14 -51.93 -0.42 0.24 16 42 H 0.07 0.50 8.14 -51.93 -0.42 0.08 16 43 H 0.08 0.65 5.65 -51.93 -0.29 0.36 16 44 H 0.39 4.75 5.96 -40.82 -0.24 4.51 16 45 H 0.06 1.10 8.14 -51.93 -0.42 0.68 16 46 H 0.06 1.13 8.14 -51.93 -0.42 0.71 16 47 H 0.06 1.11 8.14 -51.92 -0.42 0.69 16 48 H 0.06 1.11 8.14 -51.93 -0.42 0.69 16 49 H 0.02 0.45 7.20 -51.92 -0.37 0.08 16 50 H 0.01 0.35 7.57 -51.93 -0.39 -0.05 16 51 H 0.26 7.60 8.31 -40.82 -0.34 7.26 16 52 H 0.10 1.10 5.61 -51.93 -0.29 0.81 16 53 H 0.07 0.53 6.16 -51.93 -0.32 0.21 16 LS Contribution 404.58 15.07 6.10 6.10 Total: -1.00 -35.06 404.58 -8.38 -43.44 By element: Atomic # 1 Polarization: 9.10 SS G_CDS: -9.04 Total: 0.05 kcal Atomic # 6 Polarization: -0.48 SS G_CDS: -0.08 Total: -0.55 kcal Atomic # 7 Polarization: -45.91 SS G_CDS: -0.12 Total: -46.02 kcal Atomic # 8 Polarization: -19.65 SS G_CDS: -0.58 Total: -20.23 kcal Atomic # 14 Polarization: 21.88 SS G_CDS: -4.67 Total: 17.21 kcal Total LS contribution 6.10 Total: 6.10 kcal Total: -35.06 -8.38 -43.44 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. ZINC000091743137.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 -84.242 kcal (2) G-P(sol) polarization free energy of solvation -35.059 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -119.301 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -8.380 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -43.439 kcal (6) G-S(sol) free energy of system = (1) + (5) -127.681 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 0.64 seconds