Acta Chim. Slov. 2004, 51, 611-618. 611 Scientific Paper MASS-MODIFIED WIENER INDICES AND THE BOILING POINTS FOR LOWER CHLOROALKANES Teik-Cheng Lim Faculty of Engineering, Nanoscience and Nanotechnology Initiative, 9 Engineering Drive 1, National University of Singapore, S 117576, Republic of Singapore Received 12-10-2004 Abstract A criterion for obtaining the chemical graph of a chloroalkane has been proposed by replacing the edge vertex containing CH3 group with covalently bonded chlorine atom. By introducing ? as the group mass ratio of Cl to CH3 as well as a molecular mass factor M = M (?, n,nc), a mass-modified Wiener index has been proposed as Wm = MPW where W is the Wiener index obtainable from the chemical graph, whilst n and nc correspond to the number of chlorine and carbon atoms respectively in the chloroalkane molecule. The mass factor is defined in such a manner that Wm=W for alkanes in order to maintain consistency. Numerical trials on the boiling point versus Wiener index plane reveal that when P = 2 the chloroalkane curve shifts such that it coincides with the alkane curve. In view of the highly established empirical equation for alkanes it is suggested herein that, instead of obtaining separate empirical curves for haloalkanes, there is a good potential for boiling points of haloalkanes be obtainable based on the philosophy of mass factor in order to take advantage of the highly established empirical models of alkanes. Key words: boiling points, chloroalkanes, mass modification, Wiener index Introduction In spite of being in existence for the past 6 decades,1"5 interest in the Wiener and Wiener-like indices continue to flourish to this day.6"20 Although being initially applied for correlating with the boiling point of alkanes,1 the Wiener index has sparked interest in the general application of Quantitative Structure-Activity Relationship (QSAR) and Quantity Structure-Property Relationship (QSPR) for other types of molecules. For example: (a) Castro et al.21 applied the Wiener index, the Schultz index and the Randič connectivity index to find the correlation between the boiling point of saturated alcohols and these three indices;17'21 (b) Roy et al.22 used a type of Quantity Structure-Toxicity Relationship (QSTR) to model the toxicity of nitrobenzenes; T.-C. Lim: Mass-Modified Wiener Indices and the Boiling Points for Lower Chloroalkanes 612 Acta Chim. Slov. 2004, 51, 611-618. (c) Farkas et al.23 applied Quantitative Structure-Retention Relationship (QSRR) to predict gas chromatographic retention indices for saturated O-, N- and S-heterocyclic compounds; and (d) A Wiener index for cyclic structures was proposed by Randič.10 Recently Charton et al.24 appreciated the fact that topological parameters for haloalkanes and haloalkenes are function of some combination of polarizability, dipole moment, the number of halogen atoms of each type, and the ratio of the number of branches to the number of atoms along the main polymer chain. In a series of paper, Yang et al.25"28 demonstrated the usefulness of the vertex electronegativity, the normalized bond length, and the difference of interaction of electrons in the two vertices. The results exhibit good prediction of the retention index of gas chromatography of hydrocarbons, the standard formation enthalpy of methyl halides, halogen-silicon and inorganic compounds containing transition metals. In this paper, we introduce a mass-modified Wiener index for chloroalkanes whereby the atomic mass of chlorine atoms are taken into consideration for correlating with the boiling points of lower chloroalkanes. Analysis Although the Wiener index was first defined by Wiener to obtain the sum of distances between carbon atoms in saturated hydrocarbons,1 it was Hosoya who introduced the formula for Wiener index 1 " W = ^HSv (1) using the distance matrix T.-C. Lim: Mass-Modified Wiener Indices and the Boiling Points for Lower Chloroalkanes Acta Chim. Slov. 2004, 51, 611-618. 613 Customarily, the vertices refer to carbon atoms but not hydrogen atoms. Although halogen compounds are those whereby at least one hydrogen atom is replaced by a halogen atom, these halogen atoms are considered herein as vertices, unlike hydrogen atoms. The justification for such decision is that the halogen atoms are heavier than hydrogen and even the carbon atom. It follows that the halogen atoms need to be differentiated from the carbon atoms. Since it is customary that the carbon atom is considered as a vertex without any weight attached to it, we herein let the carbon atom be of unit weight such that the chlorine atoms are proportionately weighted. As such, we introduce the mass modified Wiener index Wm =MPW (3) with the molecular mass ratio defined as an + nc M =---------c- (4) n + nc where «=number of chlorine atoms, «c=number of carbon atoms and a is the mass ratio of the chlorine atom to the CH3 group, _ m a~mc+3mH ^ in which m =chlorine atomic mass, mc =carbon atomic mass and mH =hydrogen atomic mass. We appreciate the fact that hydrogen mass appears in Equation (5) because this equation considers the replacement of CH3 group with chlorine atom, whilst the hydrogen atom is not considered in Equation (4) due to the neglect of hydrogen atoms as vertices. Table 1 shows the molecular structure and boiling points of chloromethanes and chloroethanes with their corresponding alkanes on the basis of similarity in chemical graphs.30"33 Results and discussion Substituting m = 35.453, mc =12.011 and mH =1.0079 into Equation (5) gives the Cl-to-CH3 group mass ratio as a = 2.358. Based on this ratio with the number of chlorine and carbon atoms as shown on the left column of Table 1, the molecular mass ratio M is obtained from Equation (4). Needless to say, the molecular mass ratio for alkanes is M = \. Plots of the boiling point in absolute temperature, Tb, versus the T.-C. Lim: Mass-Modified Wiener Indices and the Boiling Points for Lower Chloroalkanes 614 Acta Chim. Slov. 2004, 51, 611-618. fourth root of the Wiener index for both chloroalkanes and alkanes are furnished in Figure 1. Table 1. Chloroalkanes and alkanes arranged according to molecular graph with their corresponding boiling points.30-33 Chloroalkanes Alkanes Structure Name Tb (K) Structure Name Tb (K) 0-9 Chloromethane 248.95 Dichloromethane 312.95 Trichloromethane 334.85 Tetrachloromethane 349.65 Chloroethane 285.45 1,1-dichloroethane 330.45 1,2-dichloroethane 356.65 Ethane 184.55 Propane 231.05 2-methylpropane 261.45 2,2-dimethylpropane 282.65 Propane 2-methylpropane Butane Pentachloroethane 435.15 Hexachloroethane 231.05 261.45 272.65 301.05 1,1,2-trichloroethane 386.95 >—•—• 2-methylbutane tr hl'r'2 th 41935 J^~^l 2,3-dimethylbutane 331.15 2,2,3-trimethylbutane 354.05 2,2,3,3-tetramethylbutane 373.96 We observe that data from these two classes of molecules give two separate and distinct data points, whereby for equal Wiener index, chloroalkanes exhibit higher boiling points. It is interesting to note that when we select P = 2 in Equation (3), the resulting mass-modified Wiener indices for chloroalkanes increases in such a manner that the corresponding data points fall on the same linear regression as those of alkanes, as shown in Figure 2. The best fit line gives T 92.846^My[W + 111.31 (6) T.-C. Lim: Mass-Modified Wiener Indices and the Boiling Points for Lower Chloroalkanes Acta Chim. Slov. 2004, 51, 611-618. 615 Figure 1. Boiling point of chloroalkanes (circles) and alkanes (triangles) plotted against fourth root of Wiener index. Figure 2. Boiling point of chloroalkanes (circles) and alkanes (triangles) plotted against product of square root of molecular mass ratio and fourth root of Wiener index. With reference to Figure 3, an improved curve-fit can be obtained when the boiling points are plotted against W°2283 (in which the power 0.2283 has been suggested by Mihalič and Trinajstič),32 whereby the power of P = 2 in Equation (3) is retained. Not surprisingly, a lower index would result in a higher slope, as can be seen in the following curve-fit result, T.-C. Lim: Mass-Modified Wiener Indices and the Boiling Points for Lower Chloroalkanes 616 Acta Chim. Slov. 2004, 51, 611-618. M2W\ +91.524 (7) That Equation (7) appears to exhibit close resemblance to the empirical equation by Mihalič and Trinajstič,32 Tb =108.6fF02283 +95.15 (8) in absolute temperature further suggests that empirical equations of alkanes be used for chloroalkanes by incorporating M. Both exponents of the Wiener index, i.e. 0.25 and 0.2283, were selected due to their capability in producing linear [plots of boiling points versus W1'4 or W°2283. The former exponent is convenient whilst the latter was obtained via regression analysis from numerous data of alkanes. Figure 3. Boiling point of chloroalkanes (circles) and alkanes (triangles) plotted against Wm°'22*3 whereby Wm =M2W. Conclusions and Suggestions A molecular mass modification to the Wiener index has been introduced whereby a factor M is the ratio of the chloroalkane mass in reference to a graph-similar alkane mass. Whether raised to a power of 0.25 or 0.2283, multiplication of the Wiener index with the square of M enables the collapse of chloroalkane data points onto the alkane regression fit. This implies that instead of obtaining a separate empirical equation for haloalkanes, there exists the possibility of extending the empirical equation of alkanes - T.-C. Lim: Mass-Modified Wiener Indices and the Boiling Points for Lower Chloroalkanes Acta Chim. Slov. 2004, 51, 611-618. 617 which are highly established - for application in haloalkanes by simply factoring the molecular mass ratio M raised to a certain power. In the present paper, a power of 2 has been shown to be suitable for the special čase of lower chloroalkanes. Further refmements can be made by considering higher chloroalkanes and their corresponding graph-similar alkanes. Having demonstrated the suitability of the molecular mass ratio in aligning the Tb versus M2W curve for chloroalkanes, it is suggested that similar work be performed for other haloalkanes in obtaining the parameters P for aligning their Tb versus M2W curves to that of alkanes. References 1. H. Wiener, J. Am. Chem. Soc. 1947, 69, 17-20. 2. H. Wiener, J. Am. Chem. Soc. 1947, 69, 2636-2638. 3. H. Wiener, J. Chem. Phys. 1947, 15,166-166. 4. H. Wiener, J. Phys. Chem. 1948, 52, 425-430. 5. H. Wiener, J. Phys. Chem. 1948, 52, 1082-1089. 6. B. Zhou, I. Gutman, Chem. Phys. Lett. 2004, 394, 93-95. 7. I. Gutman, Croat. Chem. Acta 2004, 77, 61-64. 8. I. Gutman, D. Vukicevic, J. Zerovnik, Croat. Chem. Acta 2004, 77, 103-109. 9. M. Randič, J. Chem. Inf. Comput. Sci. 2004, 44, 373-377. 10. M. Randič, Acta Chim. Slov. 2002, 49, 483-496. 11. S. Nikolič, N. Trinajstič, M. Randič, Chem. Phys. Lett. 2001, 333, 319-321. 12. D. Plavšič, M. Šoškic, I. Landeka, N. Trinajstič, J. Chem. Inf. Comput. Sci. 1996, 36, 1123-1126. 13. S. Nikolič, N. Trinajstič, Z. Mihalič, Croat. Chem. Acta 1995, 68, 105-129. 14. D. Plavšič, S. Nikolič, N. Trinajstič, D. J. Klein, Croat. Chem. Acta 1993, 66, 345-353. 15. A. Perdih, B. Ferdih, Acta Chim. Slov. 2002, 49, 67-110. 16. A. Perdih, B. Perdih, Acta Chim. Slov. 2003, 50, 513-538. 17. M. Kompany-Zareh,v4cta Chim. Slov. 2003, 50, 259-273. 18. P.V. Khadikar, S. Karmarkar, R. G. Varma,v4cto Chim. Slov. 2002, 49, 755-771. 19. M. V. Diudea, I. Gutman, Croat. Chem. Acta 1998, 71, 21-51. 20. A. A. Dobiynin, R. Entringer, I. Gutman, Acta Appl. Math. 2001, 66, 211-249. 21. E. A. Castro, M. Tueros, http://preprint.chemweb.com/cps/physchem/0110012. 22. K. Roy, G. Ghosh, QSAR Comb. Sci. 2004, 23, 99-108. 23. O. Farkas, K. Heberger, I. G. Zenkevich, Chemom. Intell. Lab. Syst. 2004, 72, 173-184. 24. M. Charton, B. I. Charton, J. Comput.-Aided Mol. Des. 2003, 17, 211-221. 25. F. Yang, Z. D. Wang, J. Sun, M. D. Luo, S. S. Qu, Acta Chim. Sinica 2002, 60, 764-768. 26. F. Yang, Z. D. Wang, Y. P. Huang, X. R. Ding, P. J. Zhou, J. Chem. Inf. Comput. Sci. 2003, 43, 753-756. 27. F. Yang, Z. D. Wang, Y. P. Huang, P. J. Zhou, J. Chem. Inf. Comput. Sci. 2003, 43, 1337-1341. 28. H. L. Yang, Z. D. Wang, Y. P. Huang, F. Yang, J. Wuhan. Univ. Technol. Mater. Sci. Ed. 2004, 19, 41-43. 29. H. Hosoya, Buli. Chem. Soc. Jpn. 1971, 44, 2332-2339. 30. E. Althaus, M. Jakubith, Chemistry and Chemical Engineering; VCH: Weinheim, 1993. 31. J. A. Dean, Handbook of Organic Chemistry; McGraw-Hill: New York, 1987. 32. Z. Mihalič, N. Trinajstič, J. Mol. Struct. THEOCHEM 1991, 232, 65-78. 33. D. W. Scott, D. R. Douslin, M. E. Gross, G. D. Oliver, H. M. Huff, J. Am. Chem. Soc. 1952, 74, 883-887. T.-C. Lim: Mass-Modifed JViener Indices and the Boiling Points for Lower Chloroalkanes 618 Acta Chim. Slov. 2004, 51, 611-618. Povzetek Za pripravo kemijskega grafa kloroalkana se kovalentno vezan klorov atom predstavi kot skupino CH3. Uvede se a kot masno razmerje med Cl in CH3, molekulski masni faktor M M = M(a, n ,nc) ter masno spremenjeni Wienerjev indeks Wm= hfw, kjer je W Wienerjev indeks, ki ga dobimo iz kemijskega grafa, medtem ko sta n in nc število klorovih odnosno ogljikovih atomov v molekuli kloroalkana. Masni faktor je definiran tako, da je za alkane Wm= W. Če izberemo P=2, se na diagramu vrelišč proti Wienerjevemu indeksu krivulja kloroalkanov premakne tako, da sovpada s krivuljo alkanov. Ker je empirična enačba za vrelišča alkanov dobro uveljavljena, so dobri tudi obeti za ugotavljanje vrelišč haloalkanov s pomočjo masnega faktorja in empiričnih modelov za vrelišča alkanov. T.-C. Lim: Mass-Modified Wiener Indices and the Boiling Points for Lower Chloroalkanes