Journal of Chemistry

Volume 2013, Article ID 964091, 7 pages

http://dx.doi.org/10.1155/2013/964091

## New Semiempirical Equation Describing Evaporation and Condensation in Nonassociated Liquids

^{1}Institute of Chemistry, Academy of Sciences, Academic Street 3, 2028 Kishinev, Moldova^{2}Gebze Institute of Technology, Istanbul Cd. 101, Cayirova, 41400 Kocaeli, Turkey

Received 10 August 2013; Revised 28 November 2013; Accepted 1 December 2013

Academic Editor: Pedro M. Mancini

Copyright © 2013 Mihail Yu. Gorbachev and Anatholy S. Dimoglo. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The consideration of the evaporating and condensing molecules’ interaction with the surface layer of nonassociated liquids made it possible to find an equation for relations between the saturated vapor pressure , from one side, and surface tension, critical temperature, and molar volume of the liquids, from the other side. This equation takes into account the influence of intramolecular conformational transitions of evaporating molecules on the quantity of their energy barrier. There are two types of the condensation process for the nonassociated liquids: soft and hard molecular condensation. For some vapor molecules the surface layer of the liquids behaves as an impenetrable elastic film. In the case of evaporating molecules, their one-particle potential barrier caused by the surface molecules vibration is essentially higher for the conformationally flexible molecules than that for rigid ones.

#### 1. Introduction

For nonassociated liquids in the state of phase equilibrium with their saturated vapor, both evaporation and condensation are responsible for the level of the vapor pressure. In its turn, the saturated vapor pressure of liquids () is related also to the intensity of molecular motion (i.e., molecular kinetic temperature) therein. And really, any increase in absolute temperature of a liquid leads to the corresponding increase in the motion intensity, when the number of molecules evaporating per unit time from the unit surface (and, therefore, value of ) grows [1].

Usually the dependence of on can be described by means of the following equation (Antoine’s equation [2]):

In (1) , , and are some constants which differ for various liquids. As a rule, for any liquid these constants are being defined on the basis of its experimental values. Other numerous correlations between and are similar to (1), and they are formed by means of various additional addenda depending on to its right-hand side (see [3]). However, any above-mentioned correlation applied to the calculation of for a new liquid requires the definition of its own new constants depending on the nature of this new liquid.

Furthermore, it is well known that molecules of any liquid can be divided into two groups, according to the value of their one-particle kinetic energy [1]. The first group includes the so-called “hot” molecules. Their kinetic energy is greater than or equals to some critical for jumping value . The molecules of the second group stay in their potential wells caused by the intermolecular attractive forces between these molecules. The given molecules take part in their vibration motion. The relative share of the molecules of the first group may be found by means of the Boltzmann statistics [1]: where is the Boltzmann constant and is the above-mentioned critical value of one-particle energy.

If we consider the investigated molecules as spherical particles, whose effective radius is proportional to ( is the molar volume of a liquid), then the following equation holds for the quantity [4, 5]: where denotes the surface tension of a liquid, and is a coefficient of proportionality, which is constant for many liquids (i.e., does not depend on ) [4, 5]. Then the relative share of the above-mentioned jumping molecules (the first group) can be defined by means of the following equation:

In (4) is a constant () equal to (J^{−1} K mol^{2/3}) for various nonassociated liquids.

Some part of these jumping molecules can escape the investigated liquid phase and immediately form vapor phase [1]. It has also been shown [4] that the thermal kinetic energy of a nonassociated liquid (and, therefore, its value of ) increases with the increase of value. Bearing in mind all the above reasons, the aim of the present work is to find an equation that relates the quantity with the share of jumping molecules.

#### 2. Methods of Calculations

##### 2.1. Method of Calculations Evaporation in Nonassociated Liquid

In order to find the above equation, let us consider one mole of a liquid with conformationally rigid molecules. Hereinafter we suppose that the given liquid is in thermal equilibrium with its saturated vapor (its temperature is within the interval between its normal freezing and normal boiling points). We shall consider all the molecules of the investigated liquid as spherical particles (hard spheres) due to their fast mutual reorientations. The above value of is supposed to be high enough for the consideration to be valid.

Here it is to be noted that the share of the evaporating molecules of this liquid is less than that of the molecules jumping within its volume. It has been shown [6] that the surface layer of any nonassociated liquid is formed of its vibrating molecules.

Here, vibrating molecules are molecules of liquids with weight centres staying for a while in their potential wells formed by neighbouring molecules. Vibrating molecules may become jumping ones and *vice versa*, but they have also a period of “settled life,” which allows for considering them as vibrating molecules.

Thus, molecules abandoning the liquid phase must have high enough values of their one-particle kinetic energy in order to get over the potential barrier caused by the vibrating molecules belonging to the surface layer. As far as the energy well of a vibrating molecule is proportional to (see (3)), hereinafter we shall consider that the quantity may be described by means of the following equation:

Equation (5) is similar to (3); however, there is a principal difference between them. The coefficient in (5) may depend on and on the nature of the investigated liquid phase as well. In other words, the quantity is proportional to the energy of the one-particle potential well and also may depend on the structure of the surface layer, namely, its thickness which is being defined by the chemical nature of its vibrating molecules and the value of of the given liquid phase. If one-particle kinetic energy of a jumping molecule is higher than the above mentioned its critical value , the molecule will be capable of escaping liquid phase. The relative share of such molecules may be found by means of the Boltzmann statistics [1]: where , in its turn, depends on and the nature of the liquid phase.

For different non-associated liquids it was shown [6] that molecules in their near-surface layers form a peculiar quasi crystal phase. This phase is being formed for the account of jumping molecules quitting the near-surface layers when passing to the gas phase. By this, the temperature of these layers becomes lower than that in the whole volume of the liquid.

The molecules of near-surface layers participate mainly in vibrational but not in jumping motion. Only they are responsible, from one part, for the surface tension phenomena in such liquids, and from the other part, namely, they form a potential barrier for evaporating or condensing molecules. The order in the given quasi crystal phase grows while nearing the surface layer of the liquid. From the liquid part, this area is messy enough, and its border is diffusive. In this work we suppose that the volume concentration of molecules in this area is practically equal to their concentration in the whole volume.

Now, let the total number of all jumping molecules be equal to ( is Avogadro’s number) and let their volume concentration in the given liquid phase be . Then the number of collisions (per unit time) among the above-mentioned jumping molecules with energies greater than or equal to and the vibrating ones situated on some unit area of the surface layer may be defined from the equation below [1]: where is the molecular mass of one molecule and the quantity is being defined by means of (4). At the same time, will be also equal to the number of evaporating jumping molecules, which are capable of quitting the investigated liquid through its unit surface per unit time (because their energy is high enough for getting over the surface layer).

##### 2.2. Describing Condensation in Nonassociated Liquids

Now let us consider the inverse process of molecular condensation as passing from the saturated vapor phase into the considered liquid one. Hereinafter we shall suppose that the process can be divided into two simultaneous processes. The first one is the condensation of those vapor molecules, whose one-particle kinetic energy is less than or equal to the above-mentioned critical energy (see (3)). In this case, attractive forces between vapour molecules and vibrating surface molecules are capable of fixing the former ones on the surface of the liquid. The number of the vapor molecules that are being fixed on the area unit of this surface per unit time may be easily found by means of the following equation [1]: where is the value of the total volume’s concentration for the molecules of the liquid under study in its saturated vapor phase. As far as for the considered fixed vapor molecules the critical value of one-particle kinetic energy is being defined by (3); our (8) contains as its parameter (see (4)). However, it should be noted that in (8) denotes the relative share of those vapor molecules which can be fixed on the liquid surface due to the intermolecular attractive forces. The condensation of such “slow” molecules may be called as “soft” condensation.

The presence of quasi-crystal area consisting of a few surface layers of the liquid (see [6]) is the reason of the fact that molecules moving to the surface layer (from either gas phase or liquid phase) must necessarily be influenced by a potential that has been formed by the molecules of the quasi-crystal area. This potential (from both sides of the very surface layer) will look like a Morse curve. However, the Morse curve, be it nearby the layer or in departure of it, can be represented by two parabolas, one for the area of minimum and the other for the more distant area. Thus, the energy of the quitting the surface level molecules will be described by two squared members being kinetic and potential energies of their interaction with quasi crystal molecules. According to the introduced approximation, potential energy of a molecule will quadratically depend on its distance till the very surface layer.

In this way, square roots will appear nowhere (including the equation for the particles balance as well). It becomes clear that the minimum for the Morse potential related to the gas phase is vacant, and a vapor molecule with low energy can be settled there, finishing in this way the very surface layer. This minimum lies in a space inside the gas phase but not inside the layer mentioned. The surface layer lies in the area of potential energy of the intermolecular interaction maximum. However, a softly condensed molecule, when getting in the minimum, imbeds anyway into the surface layer, sooner or later, without getting additional energy.

As a result, quasi crystal phase presence and low kinetic energy of its vibrating molecules cause soft condensation instead of evaporation. The energy of the softly condensed molecules redistributes at once by vibrating molecules of the quasi crystal area, and it is much more harder for them to quit the surface of the liquid; it is easier and more favorably to embed into it.

The second above-mentioned condensation process is the condensation of those vapor molecules, whose one-particle kinetic energy is greater than or equal to the value of (see (5)). This energy is high enough for getting over the potential barrier formed by the molecules vibrating on the surface of the liquid. The number of these “quick” vapor molecules penetrating into the volume of the liquid phase under study (per unit time and through a unit area of its surface) is defined by the equation below [1]: where may depend on and on the nature of the investigated liquid phase (the coefficient is explained before, relative to (5)). The condensation of such “quick” molecules may be called as “hard” condensation.

Thus, just the vibrating molecules that form the surface layer of any liquid can cause the two above-considered condensation processes: condensation of molecules, whose energy is less than or equal to ; condensation of those molecules, whose energy is greater than or equal to . For molecules with their one-particle kinetic energy ranging from up to , the surface layer of the liquid phase is an elastic interfacial molecular film, which reflects the given molecules back into the saturated vapor phase.

For the saturated vapor phase, which is in equilibrium state with the corresponding liquid phase, the number must be equal to the sum of numbers and . In this way, we have from (7)–(9) (after cancelling their common multiplier ) the following equality:

Next, if it has been supposed that the saturated vapor phase obeys the perfect gas law, then the concentration can be expressed by means of the following equation:
where is the saturated vapor pressure [1]. From (10) and (11) one can easily obtain the following expression used to calculate values of :
where is the gas constant ( Jmol^{−1} K^{−1}), and is being defined by means of (4). In order to use (12) in practice, we should define the dependence of the quantity on and the nature of the investigated liquid phase.

#### 3. Results and Discussion

##### 3.1. Definition of the Temperature Dependence of

To define the dependence of on , let us consider the compounds presented in Table 1. (Here it should be reminded that we consider liquids, whose molecules are conformationally rigid.)

Experimental values of , , and for all these compounds that participate in the calculation of are taken from [3, 7–11]. The comparison of numerical values of calculated with their use in (12) shows that for all compounds from Table 1 is directly proportional to , where denotes the reduced temperature of the compounds. Hereinafter denotes critical temperatures of investigated liquids. For the liquids presented in Table 1 their values of were taken from [7]. The data in the last column of Table 1 show that the ratio remains constant and does not depend on and the nature of the liquid under study. Thus, the expression for the quantity entering (12) looks like

Then (12) will have the form

Here it is worth to be noted that the value of in (14) has been defined by means of (4) with the symbol (but not ) in the name of the exponent.

##### 3.2. Towards the Evaluation of the Constant

In order to verify (13) and (14) and to define the quantity more exactly, the compounds presented in Table 2 were investigated. Experimental values of , , , and were taken from [3, 7–13].

The data in Table 2 show that does remain constant for these compounds, although their chemical structures differ significantly. The average numerical value of the constant defined from the data sets in Tables 1 and 2 is equal to (J^{−1} K mol^{−2/3}).

Now let us consider the case when our molecules have energetically degenerated (or quasi-degenerated) conformational states in the liquid phase. It has been shown [14] that the intramolecular conformational transitions occurring among these states will lead to the corresponding increase in the total number of bimolecular collisions between neighboring molecules. This increase can be described (see [14]) by means of the function , whose analytical form is

In (15) its dimensionless constant does not depend on and the nature of the liquid phase. It is equal to 0.05681 [14]. (For liquids with conformationally rigid molecules and .)

Now let us consider the process of piercing of the evaporating conformationally flexible “hot” molecules through the vibrating molecules of the surface layer. Hereinafter we shall suppose that these vibrating molecules are fixed in one of their conformational states, due to the action of the attractive force field existing between them and also due to the above-mentioned absence of the translational jumping motion for the molecules of the surface layer [6]. In other words, it is presupposed that the presence of uncompensated attractive force (the latter causes the surface tension in the investigated liquid and is directed into its volume) leads to the fixation of each molecule of the surface layer in one of its conformations. Thus, the evaporating “hot” molecules, which are capable of making transitions among their conformational states, will pierce the surface layer formed by the vibrating molecules existing in only one conformational state.

##### 3.3. Evaporating Conformationally Flexible Molecules

Further, let us consider interaction of two neighboring molecules of the liquid under study. If each of these molecules can have conformational states, then the increase in the number of their bimolecular collisions is being described by means of the function (see (15)). Now let us imagine that for one of these molecules (hereinafter, the first molecule) its conformational transitions are stopped (i.e., the molecule is fixed in one of its conformations). In this case, the increase in the number of bimolecular collisions between the considered molecules will be caused by conformational transitions of the second (conformationally flexible) molecule. As the situation under consideration is symmetrically reversible with respect to both molecules (namely, the first molecule is considered as conformationally flexible and the second one is conformationally rigid), so the above increase is to be equal to the square root of the function :

And in the case, when both neighboring molecules interact while doing their conformational transitions, the total increase in the number of their bimolecular collisions will be equal to (and, therefore, to ).

When the evaporating molecules of the liquid phase pierce through the vibrating molecules of its surface layer, we have the realization of the above-considered situation in practice. Thinking the evaporating molecules as piercing the surface layer, we deal with the translational motion of a conformationally flexible molecule through the layer of conformationally rigid ones. Hereafter we shall consider that the one-particle potential barrier (see (5)) which must be got over by this molecule is also proportional to the number of its bimolecular collisions with the corresponding molecules of the surface layer. Thus, in order to calculate the barrier adequately, we must take into account the increase in the number of the bimolecular collisions which are caused by the conformational transitions of this evaporating molecule. Bearing in mind the above-mentioned reasons related to the quantity describing this increase (see (16)), we can rewrite (5) in the following form:

Thus, for the saturated vapor pressure of those liquids, whose molecules can have their conformational transitions, the following final expression holds:

In order to verify (18) let us therein denote the product as . The values of calculated on the basis of (18) are presented in Table 3 for the liquids, whose molecules have conformational transitions.

Experimental values of , , , and required for the calculation of were taken from [3, 7, 8, 10, 15]. Here the definition of (the number of one-particle conformational states for the molecules of the investigated liquids) is to be mentioned. For any such molecule, the value of in (15) and (18) is supposed to be equal to the total number of local minima on the surface of its full internal potential energy, which arise on this surface under all possible internal rotations and intramolecular inversions that change the conformation of its skeleton [5, 16]. All its atoms are considered hereinafter as non-equivalent. Thus, in the case of ethane (see Table 3) its molecule may pass through three its possible (hindered) conformations, which arise from the internal rotations around C–C bond.

In the case of trimethylamine, the internal rotation around each of its C–N bonds can be characterized by the presence of three local minima on the surface of its internal one-particle potential energy [17]. As the one-particle number of the above C–N bonds is equal to three, so the number of the minima connected with three internal C–N rotations will be 3^{3}. For the molecule of trimethylamine we have also to take into account the spatial inversion of its nitrogen atom [17]. The given type of internal one-particle motion leads to the total number of the above mentioned minima equal to .

1,4-Dioxane may have intra-molecular conformational transitions among four its possible conformations: two conformations of the type “chair” and two conformations of the type “bath” [17]. So, for 1,4-dioxane . The usage of the above considered values of for the calculation of the ratio shows (see Table 3) that this ratio is the constant found above (cf. the data of Tables 2 and 3).

Thus, bearing in mind that , we can use (18) for the calculations of the quantity for various nonassociated liquids formed by both conformationally rigid and conformationally flexible molecules.

##### 3.4. The Calculation of the Quantity Vapor Pressure of the Inorganic and Organic Liquids

The values of the saturated vapor pressure of the inorganic and organic liquids presented in Table 4 were calculated by means of (18). Experimental values of , , needed for this calculation were taken from [1, 3, 7–11, 15] and also from [18–27]. The values of were calculated on the base of data presented in [17, 28]. Note that in the case of the benzene derivatives (see Table 4) the internal rotation around any of their C–C bonds (where one of carbon atoms belongs to the benzene ring) corresponds to six local minima of their one-particle potential energy [17].

Thus, for instance, the total number of 1,3,5-trimethylbenzene will be 6^{3}. The comparison of the calculated values of and the experimental ones taken from [3, 8, 9, 15, 21, 23] shows that (18) describes the quantity adequately. The average absolute value (module) of calculation error for is 1.6% (Table 4).

However, such low-boiling quantum liquids as helium and hydrogen do not obey (18). This equation does not feet also for liquid metals, ionic liquids, and the chemical compounds, whose molecules are strongly associated in their liquid phase: water, alcohols, and so on. Thus, for instance, for quinoline being the associated liquid (and considered at its normal boiling point—510.85 K [15]), experimental values of , , and are equal to (Nm^{−1}), (m^{3} mol^{−1}), and 782.2 K, respectively [15, 29, 30]. The value of calculated by means of (18) equals 117305 Pa (instead of 101325 Pa). At the same time, naphthalene, whose molecules are nonpolar, obeys (18) at its normal boiling point (see Table 4).

#### 4. Conclusion

The following conclusions may be drawn. Equation (18) allows us to perform the calculation of the saturated vapor pressure for various nonassociated liquids. There are two types of the condensation process of the nonassociated liquids: soft and hard molecular condensation. For a part of vapor molecules their surface layer behaves as some impenetrable elastic film. For the evaporating conformationally flexible molecules their one-particle potential barrier caused by the vibrating surface molecules is essentially higher than that for the evaporating conformationally rigid ones. The increase in the value of this barrier can be properly described by means of (17).

#### Conflict of Interests

The authors declare that there are no competing financial interests related to the publication of this paper.

#### References

- E. A. Moelwyn-Hughes,
*Physical Chemistry*, Pergamon Press, London, UK, 1961. - C. Antoine, “Tension des vapors: nouvelle relation entre les tension et les temperatures,”
*Comptes Rendus*, vol. 107, pp. 681–684, 1888. View at Google Scholar - R. C. Reid, J. M. Prausnitz, and T. K. Sherwood,
*The Properties of Gases and Liquids*, McGraw-Hill, New York, NY, USA, 1977. - M. Y. Gorbachev, “Thermal kinetic energy of liquids with conformationally rigid molecules,”
*Physics and Chemistry of Liquids*, vol. 44, no. 2, pp. 145–152, 2006. View at Publisher · View at Google Scholar · View at Scopus - M. Y. Gorbachev, “Influence of excited molecules on the process of frictional flow of nonassociated liquids,”
*Physics and Chemistry of Liquids*, vol. 42, no. 3, pp. 269–278, 2004. View at Publisher · View at Google Scholar · View at Scopus - M. Y. Gorbachev, “Dependence of surface tension of near-boiling non-associated liquids on their molar volume and some critical constants,”
*Physics and Chemistry of Liquids*, vol. 39, no. 3, pp. 315–325, 2001. View at Google Scholar · View at Scopus - J. A. Dean, Ed.,
*Handbook of Organic Chemistry*, McGraw-Hill, New York, NY, USA, 1987. - N. B. Vargaftik,
*Tables of Thermophysical Properties of Liquids and Gases*, Halsted Press, New York, NY, USA, 1975. - A. G. Morachevsk and I. B. Sladkov,
*Phisiko-Khimicheskie Svoistva Molekuljarnykh Neorganicheskikh Soedinenii*, Khimia, Sankt Peterburg, Russia, 1996. - B. P. Nikolsky, Ed.,
*Spravochnik Khimika. Tom 1*, Khimia, Moscow, Russia, 1963. - G. G. Devjatykh and A. D. Zorin,
*Letuchie Neorganicheskie Gidridy Osoboi Chistoty*, Nauka, Moscow, Russia, 1974. - S. D. Razumovskii,
*Kislorod-Elementarnye Formy i Svoistva*, Khimia, Moscow, Russia, 1979. - A. A. Ravdel and A. M. Ponomareva, Eds.,
*Kratky Spravochnik Phisikohimicheskikh Velichin*, Khimia, Leningrad, Russia, 1983. - M. Y. Gorbachev, “Influence of transitions among degenerated states of atoms and molecules on some physico-chemical quantities of liquids: Vaporization entropies of liquid 3d-metals,”
*Physics and Chemistry of Liquids*, vol. 38, no. 1, pp. 17–23, 2000. View at Google Scholar · View at Scopus - V. A. Rabinovich and Z. Ja. Khavin,
*Kratky Khimichesky Spravochnik*, Khimia, Leningrad, Russia, 1991. - M. Y. Gorbachev, “Influence of intra-molecular conformational transitions and some physico-chemical quantities of liquids on their thermal conductivity,”
*Physics and Chemistry of Liquids*, vol. 40, no. 4, pp. 395–403, 2002. View at Publisher · View at Google Scholar · View at Scopus - V. G. Dashevsky,
*Conformacii Organicheskikh Molecul*, Khimija, Moscow, Russia, 1974. - B. S. Almeida and M. M. Telo da Gama, “Surface tension of simple mixtures: comparison between theory and experiment,”
*Journal of Physical Chemistry*, vol. 93, no. 10, pp. 4132–4138, 1989. View at Google Scholar · View at Scopus - C. A. Croxton,
*Liquid State Physics—A Statistical Mechanical Introduction*, Cambridge University Press, Cambridge, UK, 1974. - A. K. Saroha, “Safe handling of chlorine,”
*Journal of Chemical Health and Safety*, vol. 13, no. 2, pp. 5–11, 2006. View at Publisher · View at Google Scholar · View at Scopus - N. B. Vargaftik,
*Spravochnik po teplophisicheskim svoistvam gazov i zhidkostei*, Gosudarstvennoe Izdatelstvo Phiziko-Matematicheskoi Lietratury, Moscow, Russia, 1963. - G. J. Janz, G. R. Lakshminaryanan, R. P. T. Tomkins, and J. Wong, “Molten Salts, v. 2 Section 2. Surface Tension Data,” National Bureau of Standarts 28, Washington, DC, USA, 1969.
- I. L. Knunjanc, Ed.,
*Khimicheskaja Enciklopedija. Tom 1*, Sovetskaja Enciklopedija, Moscow, Russia, 1988. - I. L. Knunjanc, Ed.,
*Khimicheskaja Enciklopedija. Tom 3*, Bolishaja Rossiiskaja Enciklopedija, Moscow, Russia, 1992. - I. L. Knunjanc, Ed.,
*Kratkaja Khimicheskaja Enciklopedija. Tom 4*, Sovetskaja Enciklopedija, Moscow, Russia, 1965. - I. T. Goronovsky, P. Yu. Nazarenko, and E. F. Nekrjach,
*Kratky Spravochnik Po Khimii*, Naukova Dumka, Kiev, Ukraine, 1987. - C. L. Yaw and W. Andrew,
*Thermophysical Properties of Chemicals and Hydrocarbons*, Norwich, New York, NY, USA, 2008. - S. G. Frankiss, “Vibrational spectra and structures of S
_{2}Cl_{2}, S_{2}Br_{2}, Se_{2}Cl_{2}and Se_{2}Br_{2},”*Journal of Molecular Structure*, vol. 2, no. 4, pp. 271–279, 1968. View at Google Scholar · View at Scopus - K. S. Birdi, Ed.,
*Handbook of Surface and Colloid Chemistry*, CRC Press, New York, NY, USA, 3rd edition, 2008. - D. Ambrose, “Critical temperatures of some phenols and other organic compounds,”
*Transactions of the Faraday Society*, vol. 59, pp. 1988–1993, 1963. View at Google Scholar · View at Scopus