OJPC  Vol.9 No.4 , November 2019
Taking into Account Density Fluctuations in a Solvent in a Model of Dissolution
Show more
Abstract: Earlier it was shown by different authors that there are cavities (vacancies, holes) in any liquid. The cavities should play a prominent role in dissolution processes. Nevertheless this fact was ignored in previous model of dissolution. The sizes of the cavities in different solvents containing benzene molecules were determined using solvent induced spectral shift method. The measurements of S1←S0 benzene transition spectral shifts permit to conclude that 1) macroscopic excess volumes play an almost negligible role in processes of benzene dissolution in very different solvents and 2) the minimal size of the cavity in water able to accommodate benzene molecule coincides with the solute size. Generalization of this conclusion to other nonpolar aromatics leads to evaluation contraction of the solutes under aqueous solvent influence permits to predict the solubility values of other aromatics in water and to evaluate effect of enhancement hydrate cell around these molecules on solubility.

1. Introduction

A question on quantitative prediction of solubility stands in front of scientists almost since ancient times. Nevertheless the first attempt to use quantitative parameters for qualitative prediction of solubility was done by Hildebrand in the middle of the twentieth century [1]. He introduced parameter named density of cohesion energy, δ = E / v , where Е equals to heat of evaporation and v is molar volume. Good mutual solubility of two substances is predicted when their solubility parameters coincide. Hansen [2] improved the Hildebrand approach dividing parameter δ for components in accordance with interaction types. Both two approaches are able to predict if the solubility is good or poor, but they cannot give a quantitative answer to the question “how much is solubility?”.

An attempt to answer this question was undertaken by Ben-Naim in the middle of the second half of the twentieth century [3]. According to Ben-Naim the dissolution process is realized according to the next manner. Firstly a cavity is created in a fixed position in the solvent. Then the solute which is in the fixed point in vacuum is transferred into the cavity. Finally the solute becomes free from the fixed place in the solvent. Sum of free energies of both two processes,

μ = G c + G i = k B T ln c , (1)

is called pseudo chemical potential. It is not connected with any standard state. The term G c is free energy of creating the cavity and G i is free energy of interaction between the solute and the solvent. The term c is concentration expressed in mole shares, kB is the Boltzmann constant, and T is temperature.

This equation looks correct. Nevertheless successes of its direct application are very modest (look for example, ref. [4] and references therein). The cause is simple enough. Still between the first and the second world wars Frenkel and slightly later Schottky [5] proved that cavities (holes, vacancies) must exist even in the most ideal crystals. Naturally, they must exist in liquids. The fluctuations, namely fluctuation cavities, should participate in the process of dissolution. Luck [6] gathered a lot of indirect evidence that some kind of cavities really should exist in any liquid. A problem was how to measure these cavities, especially those of them which participate in the dissolution. The problem can be solved at least in part using the solvent induced spectral shift method. The method is in essence one of reverse spectroscopic problems when the shift of electronic spectrum of dissolved molecules serves a basis for decision of a question: how solvent molecules are distributed around the solute.

2. The Solvent-Induced Spectral Shift

2.1. Model

It is convenient to adopt the simplest model of the solution at least for a beginning. The solvent is considered as continual dielectric with dielectric constant εv and refractive index n. The solute is represened by sphere of radius r which is determined according to the Dejardin et al. procedure [7]: dependence should be built of molar volume of the solute substance in liquid state on its fluidity at different temperatures V 0 , the volume at the fluidity equal to zero, is connected with r by equation

k p a c r = ( V 0 ) 1 / 3 (2)

Here kpac is packing factor. It equals to 1.88 for molecules whose shape does not sufficiently differ from spherical (Figure 1).

2.2. Interactions in the Solution and the Spectral Shift

There are two sorts of interactions in a dilute aqueous solution of nonpolar substance: solute-solvent and solvent-solvent interactions. Only the first one affects directly the spectral shift whereas the second of them affects indirectly participating in organization distribution solvent molecules around the solute. An electronic spectrum can be used for study the solution structure rather than vibrational one because in contrast to vibrational spectrum it belongs to whole molecule rather than to some of its fragments. Let the simplest case of the solute will be considered, when the solute is nonpolar molecule, and let electronic absorption of the solute is far from the solvent edge of the solvent absorption and let solute electronic states are mutually independent. Then the shift of a purely electronic or electronic-vibrational (vibronic) band in the transfer of a molecule from the gas phase to the solution may be considered as sum of the different contributions:

Δ ν = Δ ν d i s p + Δ ν e l s t + Δ ν c h e m (3)

Here Δ ν is the shift of the spectral band expressed in wave numbers, Δ ν d i s p is contribution from dispersion interactions of the solute with solvent molecules, Δ ν e l s t is that part of the shift which is determined by interaction with sources of constant electric fields (ions, dipoles,…) in the solvent, and Δ ν c h e m is the term which corresponds to those chemical interactions including hydrogen bonding that do not change individuality of the solute molecule.

Figure 1. Model, A solute (sphere of radius r) in a spherical cavity of radius R created in the solvent.

Minimum of free energy of dispersion interaction take place when the solute touches the cavity border. The share of the shift stipulated by this type of interactions is

Δ ν d i s p = C φ ( R , r ) f ( n ) , (4)

where C is a positive coefficient depending on properties of the transition in consideration, φ ( R , r ) = R 3 / [ r 3 ( 2 R r ) 3 ] is geometrical factor where R is radius of the cavity containing the solute molecule whose radius is r, f ( n ) = ( n 2 1 ) / ( n 2 + 2 ) , n is refraction index [8].

It was shown in ref. [9] that there are neither electrical nor chemical interactions between the solute and the solvent in the aqueous solution of benzene. This fact will be used below at construction the simplest version of solubility model in which density fluctuations in the solute are taken into account.

This consideration is related to mutually independent electronic transitions. When electronic states are connected by a vibration, then the low which describes the shift suffers changes [10] [11]. So the shift of S1 − S0 benzene transition is approximately expressed as

Δ ν ( 10 ) = C 1 | Δ ν 10 | 1.910 , (5)

where C1 is the positive coefficient, Δ ν 10 is the shift in the absence of electronic-vibrational coupling of the S1 electronic state with other ones, and index 1.910 is a correction for this coupling [12].

Solvents whose aromatic molecules contain oxygen make exciplexes with high-energy states of aromatic solutes [12]. This fact permits to solve reverse spectroscopy problem using only the most low-energy transitions at handling with such solvents.

2.3. Experimental Data

Experimental details including purification of substances, recording and measuring spectral shifts were done in refs. [9] [10] [11]. Data concerning solubilities are cited below.

3. Solubility

3.1. Microscopic Balance of Volumes

The average size of cavities in the solvent able to participate in the dissolution process can be found from balance of volumes at dissolution:

V 1 u v = V a v = V 1 u + V b v + V E (6)

Here V 1 u v is average volume of that cavity in the solvent, (superscript v) which contains one solute molecule (subscript 1u), V a v is the same volume in the solvent obtained after removal a of its molecules, V 1 u is average volume per one solute molecule in the solute substance, V b v is average volume of those fluctuating cavities in the solvent which participate in the dissolution process, and V E is excess volume of mixing the solution components.

Solutions of benzene in different solvents can be considered as an instructive example. One obtains from Equation (5)

| Δ ν d i s p ( 10 ) | 1 / 910 = k f ( n ) φ ( R , r ) (7)

where k = 2 824.9695422 cm 1 and r = 2.72 × 10 10 m . The results of calculation microscopic balance of volumes at benzene dissolution in different solvents are listed in Table 1.

It is readily seen from the Table that benzene forms simple intermixing only with carbon tetrachloride, R b v = 0 , not with all other considered solvents where R b v > 0 . The positive values of R b v witness that just cavities participate in dissolution process, at least, in considered cases.

Next interesting consequence is seen from the Table. The excess volume plays negligible role In all considered cases. This fact supports conclusion made above that dissolution realizes through fluctuation cavities even in cases of infinite solubility, as for example in alkanes.

3.2. Minimal Radius of the Suitable Cavity

The minimal radius of the cavity yet participating in the process of benzene dissolution can be evaluated from the function of cavity size distribution in the solvent [5]:

Table 1. Solvents, their refractive indices, n, shifts of benzene S 1 S 0 transition spectrum, Δ ν ( 10 ) , average radii of cavities containing a benzene molecule, R, excess volumes of mixing, V E , averaged over different sources, and R b v is radius of average fluctuation cavity participating in dissolution and R b v is the same without taking in consideration excess volume.

d p = d R R e x p [ G b v c ( R ) / k B T ] (8)

where G b v c is average free energy of the fluctuation cavity surface which participate in dissolution (the same indexes are used in Equation (6)). It may be expressed through microscopic surface tension [13] [14]:

G c = κ γ σ . (9)

Here σ is the cavity surface area, γ is the macroscopic surface tension, and κ is the coefficient correcting the macroscopic surface tension to the microscopic one. It is expressed as [13]:

κ 1 + ( σ 1 / σ ) ( κ 1 1 ) , (10)

where σ 1 is the area of the surface of the cavity created in the liquid as a result of removal of one of its molecules, and

κ 1 ( γ σ 1 ) 1 k B T ln ( k B T / P s V 1 ) . (11)

Here P s is pressure of saturated vapor and V 1 is the volume per one molecule in the liquid. For associated liquids

κ 1 ( γ σ 1 ) 1 k B T ln ( k B T / ξ P s V 1 ) , (12)

where ξ is average degree of association of vapor molecules [15].

Now Vbv can be expressed as

( R b v ) 2 = R m d R R 3 exp [ G c ( R ) / k B T ] R m d R R exp [ G c ( R ) / k B T ] = ( R m ) 2 + k B T 4 π κ γ (13)

Here R m is the radius of that minimum cavity which is still good for acceptance the solute, γ = 71.95 mN / m [16], κ 1 [15]. Hence R m = 2.81 × 10 10 m . This value exceeds r = 2.72 × 10 10 m [7] adopted here for benzene molecular radius only about 3%. This is too low difference for our crude model of solution. Therefore we may think that the minimal size of the cavity in water able to take the benzene molecule coincides with benzene molecule size. This conclusion is extended further to other big nonpolar solutes.

3.3. Approaches to Solubility

3.3.1. Fluctuation Approach

Let us consider low solubility of substance consisting of big nonpolar molecules, so low solubility that solute molecules do not touch each other. Let N is full amount of solvent molecules and n is amount of cavities able to accept a solute molecule. Then solubility is determined by

μ * / k T = ln ( n / N ) = ln c (14)

If the similar coincidence takes place also for other nonpolar big rigid nonpolar molecules, one has

ln ( n e / n ) = ln r e d R R exp ( 4 π γ R 2 / k B T ) r d R R exp ( 4 π γ R 2 / k B T ) = 4 π γ ( r e 2 r 2 ) / k B T . (15)

One obtains after substituting the right side of Equation (15) into Equation (14):

μ * = μ e * + 4 π γ ( r 2 r e 2 ) (16)

Here prim numbers the approach to evaluation the solubility. The results obtained with this approach are given in the third column of Table 2. They are not very significantly deviated from empirical data. Note that this approach is not connected even with phase states of solution components.

3.3.2. Energetic Approach

Free energy of molecular transfer out of the fixed position in the substance which will be dissolved into the fixed position in vacuum and then into the fixed position in the solvent, μ * , is considered in the second approach called energetic one. In the idealized case when solute properties do not change at these transitions,

μ * = G 1 u i v + G 1 u c v G 1 u i G 1 u c (17)

Here G 1 u i v is free energy of interaction (superscript i) of one solute molecule (subscript 1u) with the solvent (superscript v), G 1 u c v is free energy of creation the cavity in the solvent (superscript c) where the solute molecule can be placed, G 1 u i is free energy of interaction between the solute and its environment in the solute substance, and G 1 u c is free energy of creation the cavity instead removed solute molecule. The detailed Equation (17) looks as

Table 2. Geometrical characteristics, quasi-chemical potentials of dissolution and increase of aqueous molecules in the first hydrate cell of some substituted benzenes in water at 293 K.

μ * = G a v i + Δ G i + G a v c G b v i G b v c G 1 u i G 1 u c (18)

where G a v i is free energy of interaction between content of sphere which radius is R 1 u v before removal а solvent molecules out of it and the rest solvent, G a v c is free energy of this cavity surface, terms with subscripts bv and 1u describe similar characteristics of original vacancies in the solvent and components of that free energy which must be spend for removal one molecule out of the substance to be dissolved, respectively, and Δ G i is a correction which must be introduced into the process description after replacement claster of a solvent molecules for one solute molecule. According to [13] [14],

G i = 2 G c (19)

and Equation (18) becomes simplified to

μ * = G a v c + G b v c + Δ G i + G 1 u c (20)

When size of the cavity which contains the solute molecule is held rather by induced electrostatic forces than by collisions at thermal movement, then equilibrium takes place at

G a v c + G b v c + Δ G i = 0 (21)


μ * = G 1 u c (22)

In essence, μ * is pseudo chemical potential of transfere a molecule out of the condensed substance liable to dissolution into the solvent. It is described with the next equation:

μ * / k T = ln C . (23)

Thus it is the quasi chemical potential. The double prim numbers the approach to evaluation μ * . The data on calculated solubilities are given in the fourth column of Table 2. They also are close enough to empirical results.

When the solute molecule is transferred out of a solid phase into the solvent then Equation (23) should be specified. Zhang and Gobas [17] supposed that a surface molecule of a solid substance dissolving in a liquid is bound with other ones in that manner as molecules of super cooled liquid. Then

μ * = k T ln ( R T / p s V 1 u ) Δ G f (24)

Here Δ G f is a change of free energy at conversation the solute substance into the state of super cooled liquid which equals to

Δ G f = Δ H f ( 1 T / T m ) + T T m Δ C p d T T T T m Δ C p d T T (25)

Here Δ H f is enthalpy of solute substance fusion at melting point T m and Δ C p is change of specific heat at changing the phase state. For example, Eq. (25) may be used for evaluation that characteristic of naphthalene which is solid at the room temperature. One can evaluate function Δ C p ( T ) applying square-law extrapolation specific heat of the liquid phase taking necessary values from ref. [18]: Δ C p = 4.184 ( 1.00286 × 10 4 Т 2 2.2563 Т + 66.1362 ) J / mol K .

3.3.3. United Approach

Corrections to solute size changes should be introduced in both two above described approaches. Correcting term δ μ to the quasi-chemical potential μ * is caused by the molecular size decreasing because of pressing by reaction field forces [9]. It can be found from the next expression:

μ * = μ e * + δ μ = μ e * + 4 π [ ( r + δ r ) 2 r e 2 ] γ μ e * + 4 π ( r 2 + 2 r δ r r e 2 ) γ (26)


δ μ 8 π γ r δ r . (27)

Such correction is the negative value because its sign coincides with the sign of δ r .

Correction to μ in energetic approach looks as

μ = μ * + δ μ * (28)

where δ μ * is correction which takes into account reversible positive work making by forces of hydrophobic (electric) repulsion which compress the solute molecule.

δ μ * 4 π γ ( R 1 2 R 2 ) (29)

Here R 1 is radius of the cavity containing the solute in the case if it is not subjected to deformation, and R is the same after deformation. R can be expressed in quasi-spherical approximation as R = [ R 1 3 r 3 + ( r + δ r ) 3 ] 1 / 3 . We get after expansion R in the Tailor series and taking into account that R 3 r / 2 , and neglecting the infinitesimal terms of decomposition, that

δ μ * 8 π γ r δ r (30)

One obtains comparing Equations. (27) and (30) that

δ μ * = δ μ * (31)

It follows from Equation (31) that difference between values μ * and μ * is caused only by contraction of substituent size under reaction (reactive field) of the aqueous solvent on the solute. Hence

μ * = ( μ * + μ * ) / 2 (32)

Only deformation of the solute is taken into consideration at calculating μ * . The value of the solute size contraction is

δ r Δ μ / ( 16 π γ r ) (33)

where Δ μ = μ * μ * .

The values of μ * with values of δ r are given in the fifth column of Table 2. Firstly, it is readily seen from the Table that the predicted values of μ * are closer to the measured ones then predicted by any of above approaches. So the main points of presented consideration look right. Secondly, the substance with rigid molecules, namely naphthalene, does not show any size contraction under aqueous solvent influence. Contractions show molecules containing alkyl substituents. The more branched is a substituent, the more is contraction. This fact is evidently caused by facility of ordinary bond deformation.

3.3.4. Taking into Account Solvent Shell Strengthening

Interaction between water molecules in hydrate shell of benzene molecule is more strong then in pure water. Really, ions K+ and Cl- destroy water structure, i.e., weaken interaction between the molecule and other water molecules [19]. Nevertheless, addition salt KCl into aqueous solution of benzene does not destroys its hydrate cell, in contrast to addition such salts as RbCl and CsCl [9] which are more actively then KCl [19].

One can evaluate contribution of this enhancement into the μ* value comparing calculated values with measured ones. The corresponding correction equals to

Δ μ * / k T = 0.214 Δ N 1 + 0.055 (34)

Here N 1 is amount of water molecules in the first hydrate shell of benzene and Δ N 1 is the change of this amount after transition to another solute. Approximately

N 1 = 4 π [ R + ( V 1 u ) 1 / 3 / 2 ] 2 / ( V 1 u ) 2 / 3 , (35)

where R 1.5 r . The correlation factor of dependence (34) equals to ρ = 0.943 , of root mean square deviation of coefficient at Δ N 1 σ A = 0.080 and of free term σ В = 0.028 . One can readily see from Equation (34) that interaction between water molecules in the first hydrate shell is enhanced owing to interaction with the solute and lowering the solvent free energy per one its molecule in the solvent shell equals approximately 0.2kT. The low value of the free term in the right side of Equation (34) witness that adopted approximation is correct (see Table 3).

4. Conclusions

The above consideration clearly shows that fluctuations of density such as vacancies (holes, cavities) in diverse solvents should be taken into account at evaluations solubility of different solutes. This fact leads to paradoxy at the first glance conclusion that the excess volume plays a very modest role in microscopic balance of volumes at dissolution process, at least, in considered cases.

The calculated values of solubility based on values of cavities obtained from spectral shifts data are close enough to empirical ones, and evaluated size

Table 3. Geometrical characteristics and quasi-chemical potentials of dissolution of some alkyl substituted benzenes in water at 293 K and solute contractions.

changes, and change of strength of hydrogen bonds in the solvent near the solutes look right. So, the main points of the above consideration are correct.

In essence, the called conclusions are obtained owing to taking into account, evidently or not evidently, effect of electric field described in ref. [9]. Hence such effect should be taken into consideration in advanced models of solubility, including computer simulation.

Cite this paper: Ar’ev, I. (2019) Taking into Account Density Fluctuations in a Solvent in a Model of Dissolution. Open Journal of Physical Chemistry, 9, 204-215. doi: 10.4236/ojpc.2019.94012.

[1]   Hildebrand, J.H. (1916) Solubility. Journal of the American Chemical Society, 38, 1452-1473.

[2]   Hansen, C.M. (2000) Hansen Solubility Parameters. A User’s Handbook. CRC Press, Boca Raton, London, New York, Washington DC.

[3]   Ben-Naim, A. (2006) Molecular Theory of Solutions. Oxford University Press, Oxford.

[4]   Graziano, G.J. (1998) On the Size Dependence of Hydrophobic Hydration. Journal of the Chemical Society, Faraday Transactions, 94, 3345-3352.

[5]   Frenkel, Y.I. (1975) Kinetic Theory of Liquids. Nauka, Leningrad.

[6]   Luck, W.A.P. (1981) Einfache Modelle für kalorische Eigenschaften der Flüssigkeiten. Berichte der Bunsengesellschaft für physikalische Chemie, 83, 859.

[7]   Dejardin, L.J., Marrony, R., Delseny, C., Brunet, S. and Berge, R. (1981) ne nouvelle définition des volumes libres dans les liquides non associés: Application à la détermination des diamètres moléculaires. Rheologica Acta, 20, 497-500.

[8]   Ar’ev, I.A. (1987) Investigation the benzene hydration with spectral shift method. Theoretical and Experimental Chemistry, 23, 329.

[9]   Ar’ev, I.A. and Chernova, L.G. (2011) Interaction of Benzene with Water and Aqueous Solutions of Alkali Metal Chlorides. Russian Journal of Physical Chemistry A, 85, 1592.

[10]   Ar’ev, I.A., Dyadyusha, G.G. and Makhlinets, N.V. (1983) Effect of Spin-Orbital Interaction in Mono- and p-Dihalogenobezene Molecule on Their Electronic Spectrum Shifts under Solvent Action. Optics and Spectroscopy, 55, 285.

[11]   Ar’ev, I.A., Dyadyusha, G.G., Klimusheva, G.V. and Soroka, G.M. (1983) Effect of Spin-Orbital Interaction in Mono- and p-Dihalogenobenzenes in the S1 State under the Action of Environment. Optics and Spectroscopy, 55, 653.

[12]   Ar’ev, I.A., Lebovka, N.I. and Solovieva, E.A. (2013) Effects of Partial Charge-Transfer Solute—Solvent Interactions in Absorption Spectra of Aromatic Hydrocarbons in Aqueous and Alcoholic Solutions. Molecular Physics, 111, 3077-3080.

[13]   Sinanoglu, O. (1981) Microscopic Surface Tension down to Molecular Dimensions and Microthermodynamic Surface Areas of Molecules or Clusters. The Journal of Chemical Physics, 75, 463.

[14]   Sinanoglu, O. (1981) What Size Cluster Is Like a Surface? Chemical Physics Letters, 81, 188-190.

[15]   Ar’ev, I.A. (2014) Correcting the Microscopic Coefficient of Surface Tension of Associated Liquids. Russian Journal of Physical Chemistry A, 88, 173-174.

[16]   Vargaftik, N.B., Volyak, L.D. and Volkov, B.N. (1975) Surface Tension of Water at Temperatures from 0 up to 370 Degrees of Celsium. In: Surface Phenomena in Liquids, Leningrad University, Leningrad, 180-192. (In Russian)

[17]   Zhang, X. and Gobas, F.A.P.C. (1995) A Thermodynamic Analysis of the Relationships between Molecular Size, Hydrophobicity, Aqueous Solubility and Octanol=Water Partitioning of Organic Chemicals. Chemosphere, 31, 3501.

[18]   Nikolskii, B.P. (1966) Handbook for Chemist. Vol. 1, 2, Khimia, Moscow, Leningrad. (In Russian)

[19]   Samoylov, O.Ya. (1957) The Structure of Aqueous Solutions of Electrolytes and Hydration of Ions. Acad. Sci. of USSR, Moscow. (In Russian)