Wall clock time and date at job start Tue Mar 31 2020 21:13:20 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 N 2 2 C 1.31619 * 1 3 3 Si 1.86807 * 119.99850 * 2 1 4 4 H 1.48496 * 109.47140 * 90.00064 * 3 2 1 5 5 H 1.48495 * 109.47019 * 209.99946 * 3 2 1 6 6 H 1.48494 * 109.47049 * 330.00193 * 3 2 1 7 7 C 1.38806 * 120.00121 * 180.02562 * 2 1 3 8 8 H 1.08001 * 120.00071 * 150.71855 * 7 2 1 9 9 N 1.36938 * 120.00062 * 330.71806 * 7 2 1 10 10 C 1.47420 * 125.64663 * 345.34881 * 9 7 2 11 11 C 1.54910 * 104.83278 * 204.08974 * 10 9 7 12 12 C 1.55154 * 101.58159 * 322.99043 * 11 10 9 13 13 H 1.09002 * 110.71992 * 277.13677 * 12 11 10 14 14 C 1.53000 * 110.72565 * 154.04028 * 12 11 10 15 15 C 1.53035 * 109.46456 * 180.02562 * 14 12 11 16 16 C 1.53193 * 109.31534 * 178.60321 * 15 14 12 17 17 N 1.46923 * 108.76590 * 54.65646 * 16 15 14 18 18 C 1.34766 * 120.63409 * 126.36412 * 17 16 15 19 19 O 1.21282 * 120.00778 * 175.70843 * 18 17 16 20 20 C 1.50702 * 119.99889 * 355.70764 * 18 17 16 21 21 C 1.53001 * 115.55137 * 294.41229 * 20 18 17 22 22 C 1.52997 * 109.47323 * 300.31247 * 21 20 18 23 23 F 1.39900 * 109.47338 * 185.45550 * 22 21 20 24 24 C 1.53002 * 117.50287 * 148.74229 * 20 18 17 25 25 C 1.53007 * 117.50028 * 80.08219 * 20 18 17 26 26 C 1.46930 * 118.73769 * 306.38309 * 17 16 15 27 27 C 1.53043 * 109.45447 * 59.97222 * 14 12 11 28 28 C 1.47019 * 125.65056 * 165.03406 * 9 7 2 29 29 H 0.97007 * 120.00141 * 175.84791 * 1 2 3 30 30 H 1.09002 * 110.36645 * 322.92355 * 10 9 7 31 31 H 1.08996 * 110.37211 * 85.14365 * 10 9 7 32 32 H 1.09004 * 111.00614 * 81.06216 * 11 10 9 33 33 H 1.09004 * 111.00978 * 204.87896 * 11 10 9 34 34 H 1.08997 * 109.45843 * 300.00595 * 14 12 11 35 35 H 1.09003 * 109.49820 * 298.56871 * 15 14 12 36 36 H 1.09000 * 109.49730 * 58.66598 * 15 14 12 37 37 H 1.08999 * 109.59500 * 294.87012 * 16 15 14 38 38 H 1.08999 * 109.70148 * 174.50989 * 16 15 14 39 39 H 1.09002 * 109.47001 * 60.30911 * 21 20 18 40 40 H 1.09001 * 109.47202 * 180.30676 * 21 20 18 41 41 H 1.09007 * 109.46801 * 305.45299 * 22 21 20 42 42 H 1.08992 * 109.47311 * 65.44913 * 22 21 20 43 43 H 1.09004 * 117.49191 * 359.97438 * 24 20 18 44 44 H 1.09000 * 117.49979 * 145.02171 * 24 20 18 45 45 H 1.08993 * 117.49858 * 214.98507 * 25 20 18 46 46 H 1.08997 * 117.49430 * 359.97438 * 25 20 18 47 47 H 1.08998 * 109.58999 * 293.82296 * 26 17 16 48 48 H 1.08998 * 109.70211 * 173.46279 * 26 17 16 49 49 H 1.09005 * 109.49870 * 301.33370 * 27 14 12 50 50 H 1.08990 * 109.49804 * 61.40698 * 27 14 12 51 51 H 1.09004 * 109.88431 * 298.40148 * 28 9 7 52 52 H 1.09003 * 109.87822 * 59.56340 * 28 9 7 Note: An asterisk (*) indicates this parameter was or will be optimized. Cartesian coordinates (angstroms) NO. Atom X Y Z 1 7 0.0000 0.0000 0.0000 2 6 1.3162 0.0000 0.0000 3 14 2.2502 1.6178 0.0000 4 1 2.4977 2.0465 1.4000 5 1 3.5477 1.4403 -0.7000 6 1 1.4476 2.6527 -0.7000 7 6 2.0102 -1.2021 -0.0005 8 1 2.9868 -1.2616 -0.4580 9 7 1.4568 -2.3125 0.5790 10 6 0.2983 -2.3188 1.4907 11 6 0.4531 -3.6183 2.3196 12 6 1.0142 -4.6003 1.2575 13 1 0.2074 -5.0228 0.6585 14 6 1.8296 -5.7104 1.9236 15 6 2.3688 -6.6603 0.8518 16 6 3.2134 -7.7503 1.5193 17 7 2.4046 -8.4068 2.5554 18 6 2.2527 -9.7459 2.5621 19 8 1.6481 -10.2832 3.4657 20 6 2.8343 -10.5758 1.4466 21 6 2.2189 -10.3477 0.0645 22 6 0.7235 -10.6674 0.1130 23 9 0.1857 -10.5794 -1.1755 24 6 3.2977 -11.9933 1.7883 25 6 4.3289 -10.8968 1.5102 26 6 1.7642 -7.5757 3.5839 27 6 0.9360 -6.4875 2.8932 28 6 1.9155 -3.6965 0.3900 29 1 -0.4851 -0.8379 -0.0608 30 1 0.3250 -1.4467 2.1440 31 1 -0.6316 -2.3380 0.9224 32 1 1.1619 -3.4790 3.1359 33 1 -0.5115 -3.9587 2.6963 34 1 2.6628 -5.2701 2.4713 35 1 1.5354 -7.1200 0.3204 36 1 2.9861 -6.1020 0.1481 37 1 4.0961 -7.3012 1.9747 38 1 3.5200 -8.4850 0.7748 39 1 2.3585 -9.3070 -0.2281 40 1 2.7062 -10.9979 -0.6620 41 1 0.5801 -11.6770 0.4983 42 1 0.2204 -9.9545 0.7662 43 1 3.1824 -12.3254 2.8201 44 1 3.1645 -12.7633 1.0283 45 1 4.8737 -10.9457 0.5674 46 1 4.8916 -10.5073 2.3585 47 1 2.5287 -7.1124 4.2076 48 1 1.1128 -8.1932 4.2023 49 1 0.5334 -5.8060 3.6427 50 1 0.1161 -6.9482 2.3424 51 1 1.8250 -3.9770 -0.6594 52 1 2.9523 -3.7907 0.7130 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 ZINC001733265948.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:13:20 Heat of formation + Delta-G solvation = -54.080370 kcal Electronic energy + Delta-G solvation = -30871.024960 eV Core-core repulsion = 26785.401868 eV Total energy + Delta-G solvation = -4085.623091 eV No. of doubly occupied orbitals = 65 Molecular weight (most abundant/longest-lived isotopes) = 338.211 amu Computer time = 1.48 seconds Orbital eigenvalues (eV) -47.87355 -40.96728 -38.92287 -37.39628 -36.73039 -33.93539 -31.97079 -31.75545 -31.24609 -29.17866 -27.99872 -25.72724 -24.53782 -24.02600 -23.33455 -22.92676 -22.36425 -20.65204 -20.05100 -19.58671 -18.85900 -17.49262 -17.27708 -16.73945 -16.46699 -16.06701 -15.45257 -15.21521 -14.97190 -14.85419 -14.61449 -14.36305 -13.98790 -13.92668 -13.69670 -13.21296 -13.07600 -13.03442 -13.00544 -12.71694 -12.45660 -12.17763 -12.02359 -11.86672 -11.65180 -11.45478 -11.28684 -11.19121 -11.09437 -10.91332 -10.74568 -10.61843 -10.49718 -10.43177 -10.15932 -10.06715 -10.00015 -9.93986 -9.83983 -9.51981 -8.86911 -8.35240 -6.90294 -6.10393 -3.09331 2.23414 2.63871 3.16478 3.51603 3.54122 3.70526 3.89662 4.06915 4.25899 4.27855 4.52567 4.58209 4.60874 4.66228 4.71409 4.73879 4.78821 4.90528 4.99364 5.08505 5.14796 5.15995 5.26059 5.33059 5.35888 5.36908 5.43917 5.44894 5.50710 5.57345 5.57991 5.64009 5.66390 5.75923 5.76655 5.83660 5.91452 6.01424 6.07498 6.16056 6.29461 6.33772 6.42228 6.48664 6.65874 6.88332 7.02342 7.08832 7.12490 7.82644 7.87478 8.51603 9.05713 9.85258 10.41831 11.24852 Molecular weight = 338.21amu Principal moments of inertia in cm(-1) A = 0.019143 B = 0.002473 C = 0.002418 Principal moments of inertia in units of 10**(-40)*gram-cm**2 A = 1462.305273 B =11320.662216 C =11577.804115 Net atomic charges, atomic populations, and dipole contributions calculated with CM2 Atom NO. Type Charge No. of electrons 1 N -0.918 5.918 2 C -0.027 4.027 3 Si 1.067 2.933 4 H -0.294 1.294 5 H -0.293 1.293 6 H -0.290 1.290 7 C -0.338 4.338 8 H 0.065 0.935 9 N -0.605 5.605 10 C 0.144 3.856 11 C -0.129 4.129 12 C -0.094 4.094 13 H 0.069 0.931 14 C -0.072 4.072 15 C -0.132 4.132 16 C 0.102 3.898 17 N -0.602 5.602 18 C 0.550 3.450 19 O -0.540 6.540 20 C -0.161 4.161 21 C -0.114 4.114 22 C 0.064 3.936 23 F -0.216 7.216 24 C -0.143 4.143 25 C -0.154 4.154 26 C 0.111 3.889 27 C -0.128 4.128 28 C 0.143 3.857 29 H 0.269 0.731 30 H 0.049 0.951 31 H 0.023 0.977 32 H 0.063 0.937 33 H 0.060 0.940 34 H 0.074 0.926 35 H 0.066 0.934 36 H 0.082 0.918 37 H 0.071 0.929 38 H 0.090 0.910 39 H 0.089 0.911 40 H 0.087 0.913 41 H 0.078 0.922 42 H 0.072 0.928 43 H 0.105 0.895 44 H 0.105 0.895 45 H 0.107 0.893 46 H 0.101 0.899 47 H 0.068 0.932 48 H 0.090 0.910 49 H 0.080 0.920 50 H 0.065 0.935 51 H 0.018 0.982 52 H 0.023 0.977 For the charges calculated by CM2: Dipole moment (debyes) X Y Z Total from point charges 6.952 -24.002 4.225 25.344 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 N -0.559 5.559 2 C -0.229 4.229 3 Si 0.897 3.103 4 H -0.221 1.221 5 H -0.219 1.219 6 H -0.217 1.217 7 C -0.466 4.466 8 H 0.083 0.917 9 N -0.332 5.332 10 C 0.018 3.982 11 C -0.168 4.168 12 C -0.113 4.113 13 H 0.087 0.913 14 C -0.091 4.091 15 C -0.170 4.170 16 C -0.020 4.020 17 N -0.334 5.334 18 C 0.337 3.663 19 O -0.417 6.417 20 C -0.164 4.164 21 C -0.151 4.151 22 C 0.009 3.991 23 F -0.197 7.197 24 C -0.180 4.180 25 C -0.191 4.191 26 C -0.011 4.011 27 C -0.166 4.166 28 C 0.018 3.982 29 H 0.081 0.919 30 H 0.067 0.933 31 H 0.040 0.960 32 H 0.082 0.918 33 H 0.079 0.921 34 H 0.092 0.908 35 H 0.085 0.915 36 H 0.100 0.900 37 H 0.089 0.911 38 H 0.108 0.892 39 H 0.108 0.892 40 H 0.106 0.894 41 H 0.096 0.904 42 H 0.091 0.909 43 H 0.124 0.876 44 H 0.123 0.877 45 H 0.125 0.875 46 H 0.119 0.881 47 H 0.086 0.914 48 H 0.108 0.892 49 H 0.098 0.902 50 H 0.083 0.917 51 H 0.036 0.964 52 H 0.041 0.959 Dipole moment (debyes) X Y Z Total from point charges 6.503 -22.528 4.673 23.909 hybrid contribution 0.167 -1.487 -0.275 1.521 sum 6.670 -24.015 4.397 25.309 Atomic orbital electron populations 1.70761 1.06670 1.30202 1.48311 1.25716 0.96416 1.01019 0.99717 0.83472 0.74492 0.75827 0.76471 1.22132 1.21922 1.21666 1.19760 1.00565 0.88319 1.37919 0.91724 1.44593 1.35288 1.01198 1.52086 1.21512 0.90765 0.94646 0.91281 1.22845 0.99642 0.96832 0.97443 1.22171 0.97801 0.94287 0.96995 0.91270 1.21107 0.97895 0.94626 0.95496 1.22022 1.00153 0.96409 0.98450 1.21955 0.93699 0.95953 0.90409 1.48071 1.50345 1.05568 1.29442 1.19472 0.79404 0.81833 0.85544 1.90695 1.42427 1.71526 1.37102 1.21282 0.99299 1.00378 0.95435 1.21814 0.95185 1.04386 0.93752 1.21912 0.93335 1.03605 0.80292 1.92783 1.87556 1.97940 1.41439 1.23072 1.00893 0.90151 1.03875 1.23135 0.89422 1.02367 1.04203 1.21680 0.96932 0.92497 0.90001 1.21986 0.98273 0.96988 0.99347 1.21310 0.95504 0.85592 0.95819 0.91939 0.93293 0.95958 0.91824 0.92138 0.90806 0.91514 0.89981 0.91061 0.89213 0.89244 0.89422 0.90393 0.90920 0.87621 0.87685 0.87460 0.88104 0.91368 0.89168 0.90174 0.91664 0.96394 0.95873 Number of geometries 1 Number of calculations of the screened coulomb radii 1 The total number of SCF iterations 18. 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 N -0.92 -28.13 12.93 55.49 0.72 -27.42 16 2 C -0.03 -0.78 5.05 -17.43 -0.09 -0.87 16 3 Si 1.07 25.67 30.01 -169.99 -5.10 20.57 16 4 H -0.29 -6.90 7.11 56.52 0.40 -6.50 16 5 H -0.29 -6.65 7.11 56.52 0.40 -6.25 16 6 H -0.29 -7.03 7.11 56.52 0.40 -6.63 16 7 C -0.34 -9.66 10.02 -15.06 -0.15 -9.81 16 8 H 0.06 1.84 7.98 -52.49 -0.42 1.43 16 9 N -0.60 -15.42 3.31 -170.07 -0.56 -15.98 16 10 C 0.14 3.35 5.27 -2.53 -0.01 3.34 16 11 C -0.13 -2.24 5.37 -24.58 -0.13 -2.37 16 12 C -0.09 -1.49 2.47 -88.78 -0.22 -1.71 16 13 H 0.07 1.12 8.14 -51.93 -0.42 0.69 16 14 C -0.07 -0.87 1.57 -90.58 -0.14 -1.01 16 15 C -0.13 -1.36 5.17 -26.61 -0.14 -1.50 16 16 C 0.10 0.87 5.75 -3.71 -0.02 0.85 16 17 N -0.60 -5.98 2.97 -172.76 -0.51 -6.49 16 18 C 0.55 5.95 6.22 -10.99 -0.07 5.88 16 19 O -0.54 -7.67 16.13 5.56 0.09 -7.58 16 20 C -0.16 -1.25 0.82 -155.65 -0.13 -1.38 16 21 C -0.11 -0.83 5.20 -26.73 -0.14 -0.96 16 22 C 0.06 0.56 7.27 37.16 0.27 0.83 16 23 F -0.22 -2.41 17.36 2.25 0.04 -2.37 16 24 C -0.14 -0.95 9.77 -26.69 -0.26 -1.22 16 25 C -0.15 -0.91 9.47 -26.69 -0.25 -1.16 16 26 C 0.11 1.14 6.42 -3.70 -0.02 1.12 16 27 C -0.13 -1.38 5.16 -26.60 -0.14 -1.52 16 28 C 0.14 2.91 5.99 -3.06 -0.02 2.89 16 29 H 0.27 8.26 5.98 -40.82 -0.24 8.01 16 30 H 0.05 1.26 7.86 -51.93 -0.41 0.85 16 31 H 0.02 0.56 6.83 -51.93 -0.35 0.21 16 32 H 0.06 1.07 8.05 -51.93 -0.42 0.65 16 33 H 0.06 0.92 8.08 -51.93 -0.42 0.50 16 34 H 0.07 0.93 8.08 -51.93 -0.42 0.51 16 35 H 0.07 0.71 8.14 -51.93 -0.42 0.29 16 36 H 0.08 0.88 7.45 -51.93 -0.39 0.49 16 37 H 0.07 0.57 8.14 -51.93 -0.42 0.15 16 38 H 0.09 0.63 4.30 -51.93 -0.22 0.41 16 39 H 0.09 0.66 5.95 -51.93 -0.31 0.35 16 40 H 0.09 0.52 8.11 -51.93 -0.42 0.10 16 41 H 0.08 0.61 8.14 -51.93 -0.42 0.19 16 42 H 0.07 0.71 7.74 -51.93 -0.40 0.31 16 43 H 0.11 0.83 8.00 -51.93 -0.42 0.42 16 44 H 0.10 0.53 8.14 -51.93 -0.42 0.11 16 45 H 0.11 0.46 8.14 -51.93 -0.42 0.04 16 46 H 0.10 0.62 8.14 -51.93 -0.42 0.20 16 47 H 0.07 0.65 8.14 -51.93 -0.42 0.23 16 48 H 0.09 1.01 7.04 -51.93 -0.37 0.65 16 49 H 0.08 0.86 7.19 -51.93 -0.37 0.48 16 50 H 0.06 0.74 8.14 -51.93 -0.42 0.32 16 51 H 0.02 0.38 8.14 -51.93 -0.42 -0.04 16 52 H 0.02 0.46 7.77 -51.93 -0.40 0.06 16 LS Contribution 398.84 15.07 6.01 6.01 Total: -1.00 -33.67 398.84 -9.99 -43.65 By element: Atomic # 1 Polarization: 7.22 SS G_CDS: -9.01 Total: -1.79 kcal Atomic # 6 Polarization: -6.94 SS G_CDS: -1.66 Total: -8.60 kcal Atomic # 7 Polarization: -49.54 SS G_CDS: -0.36 Total: -49.89 kcal Atomic # 8 Polarization: -7.67 SS G_CDS: 0.09 Total: -7.58 kcal Atomic # 9 Polarization: -2.41 SS G_CDS: 0.04 Total: -2.37 kcal Atomic # 14 Polarization: 25.67 SS G_CDS: -5.10 Total: 20.57 kcal Total LS contribution 6.01 Total: 6.01 kcal Total: -33.67 -9.99 -43.65 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. ZINC001733265948.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 -10.426 kcal (2) G-P(sol) polarization free energy of solvation -33.666 kcal (3) G-ENP(sol) elect.-nuc.-pol. free energy of system -44.092 kcal (4) G-CDS(sol) cavity-dispersion-solvent structure free energy -9.988 kcal (5) G-P-CDS(sol) = G-P(sol) + G-CDS(sol) = (2) + (4) -43.654 kcal (6) G-S(sol) free energy of system = (1) + (5) -54.080 kcal 1SCF WAS SPECIFIED, SO GEOMETRY OPTIMIZATION WAS NOT USED Total computer time = 1.48 seconds