# Determination of Molecular Diameter byPVT

Authors: Vahid Moeini; Mehri Deilam
Journal: ISRN Physical Chemistry
 (2012)

---

## Abstract

We derive an equation for calculation of molecular diameter of dense fluids, with using simultaneous Lennard-Jones (12-6) potential function and the internal pressure results. Considering the internal pressure by modeling the average configurational potential energy and then taking its derivative with respect to volume to a minimum point of potential energy has been shown that molecular diameter is function of the resultant of the forces of attraction and the forces of repulsion between the molecules in a fluid. The regularity is tested with experimental data for 10 fluids including Ar, N2, CO, CO2, CH4, C2H6, C3H8, C4H10, C6H6, and C6H5CH3. These problems have led us to try to establish a function for the accurate calculation of the molecular diameter based on the internal pressure theory for different fluids. The relationship appears to hold both compressed liquids and dense supercritical fluids.


---

## Body

## 1. Introduction

It is thus seen that all gases obey the same equation of state within the accuracy of the relation when there are no arbitrary constants specific to the individual gas. This is a statement of the principle of corresponding states. The ability of the law to predict experimental behavior is nicely shown in papers, where the reduced quantity is plotted against the reduced density at various reduced temperatures [1, 2].The principle of corresponding states, introduced by J. D. van der Waals in 1880, is valid for two or more fluids with potentials of the same form but different scale parameters (sayε, the depth of the minimum of potential, and rmin, the position of this minimum). This principle says that the thermodynamic properties of these fluids are equal when they are written in units of the ε and rmin of each fluid. The potential functions of fluids satisfying the principle of corresponding states are said to be conformal [3–5].It will have become clear that the difficulties in a statistical-mechanical treatment of pure fluids are greatest at low temperatures and high densities, that is, in the liquid region. Several models of the liquid state have been developed which provide approximate methods for calculating properties and which have some utility as bases for empirical correlations. Lattice-cell models represent one class of such theories. As the name implies, these theories are rather similar to those that might apply to the solid state, where the molecules spend most of the time at or near one location in the volume of the assemblage [6]. Despite the fact that the description of dense fluids at the molecular level is very complicated, they obey some relatively simple rules (regularities). Among these, we may refer to (1), the Taite equation that was presented by Tait about 100 years ago, (2) the linearity of the bulk modulus versus pressure for each isotherm.A general regularity was recently reported for dense fluids, according to which((∂E/∂V)T/ρRT)V2 is linear with respect to ρ2 for each isotherm, where the quantity (∂E/∂V)T is known as the internal pressure. This equation of state works very well for all types of dense fluids, for densities greater than the Boyle density. The regularity was originally suggested on the basis of a simple lattice-type model applied to a Lennard-Jones (12-6) fluid [7, 8].The purpose of this paper is to examine whether by the regularity, we will be able to predict of molecular diameter. The reason for expecting it to hold for fluids is the success of minimum point potential energy for fluid, in which fluids follow the same equation. In particular, we address the following points.(1)
Initially, Lennard-Jones (12-6) potential function and the internal pressure regularity have been used to derive a equation for calculation of molecular diameter.(2)
Experimental tests of the equation for calculation of molecular diameter.(3)
If the equation holds, determine the temperature dependences of the molecular diameter. We find that the simple model that was used to predict the temperature dependences of the interceptA and slope B for substances can also be used to predict their temperature dependences for molecular diameter, with good agreement with experiment. The results of these calculations are also compared with other calculations which obtained using the statistical mechanical theory for calculation of molecular diameter.


## 2. Theory

Regularities, which were originally devised for normal dense fluids, are based on the cell theory and consider only nearest adjacent interaction. Regularities are applied well to Lennard-Jones fluids which their interaction potential can be modeled rather accurately by the (12-6) powers of inverse intermolecular distance. Lennard-Jones potential function suitably describes the interactions between the molecules of a fluid under the condition that it behaves as a normal fluid. These regularities were attempted to calculate the thermodynamic quantities by modeling the average configurationally potential energy and then taking its derivative with respect to volume [9, 10].In the present work, using the experimentalPVT, the molecular diameter will be calculated in the following stages.Stage 1.
It assumes that any kinetic energy contribution to the internal energyE will vanish on taking the derivative, since the temperature is held constant. We approximate the average potential energy by summing the contribution from nearest neighbors only, assuming single inverse powers for the effective repulsion and attraction, that is,
(1)U=N2z(ρ)(Cir̅i-Cjr̅j),
where U is the total potential energy among N molecules, z(ρ) is the average number of nearest neighbors, r̅ is average distance between nearest neighbors, and Ci and Cj are constants. It is wellknown that z(ρ) to be proportional to ρ, as is the case for liquid argon, rubidium, cesium, and taken v≈r̅3, so that U can be written as [11]
(2)UN=Kiv(i/3)+1-Kjv(j/3)+1,
where Ki and Kj are constants. Carrying out the differentiation, we obtain the internal pressure for dense fluids. 							(3)pint⁡=(∂E∂V)T≈[∂(U/N)∂v]T=A1ρ(i/3)+2-B1ρ(j/3)+2,
where A1 and B1 are constants and in order to model the experimental results, it should take j=3 and i=9. Combining the foregoing results, it finds the internal pressure for dense fluids based on Lennard-Jones (12-6) potential function. It shows that the isotherm ((∂E/∂V)T/ρRT)V2 is linear function of ρ2 and so predicted the temperature dependence of its parameters [8]. The final result is therefore of the form:
(4)((∂E/∂V)TρRT)V2=A+Bρ2,
where the intercept A and the slope B significantly are related to attraction and repulsion, respectively, and both depend on temperature. In practice, it has been shown that A=A1/RT and B=B2-B1/RT. Thus, the intercept A is proportional to 1/T and the slope B is functional in 1/T. We now consider the density problem where the system is prepared in a minimum state. Carrying out the differentiation (∂E/∂V)T=0, we obtain the density for each fluid at the position of this minimum.
(5)A+B(ρmin⁡ρC)2=0.
The minimum value of this density can also be written in the familiar form:
(6)ρmin⁡=(-AB)1/2ρC.Stage 2.
The Lennard-Jones potential function has been most widely used for thermodynamic and transport property prediction of substances in both pure and mixture states and in all phases of matter. However, the Lennard-Jones potential function can be written as pair-intermolecular potential to be used in thermophysical property computations, that is,(7)U(r)=Aε[(σr)m-(σr)n],
where A is a constant, ε is the potential well depths, and σ is the diameter of molecule. The potential function obeys boundary conditions such that [12]
(8)A=(mm-n)(mn)n/(m-n).
By differentiation of the expression for the Lennard-Jones potential show that the distance rmin at the minimum is
(9)rmin⁡=(mn)1/(m-n)σ,
in order to model the experimental results, it should take m=12 and n=6. Where the position of this minimum of the system is given by
(10)rmin⁡=(1ρmin⁡)1/3,σ=121/6rmin⁡.
At present work, with using simultaneous of (6) of the Stage 1 and (10) of the Stage 2, the molecular diameter of each fluid can be obtained which provides a simple method to determine the molecular diameter. A typical curve of molecular diameter as a function of temperature for Ar is shown in Figure 1. This result is quite remarkable. It means that the molecular diameter of each fluid is given by the intercept A and the slope B significantly that are related to attraction and repulsion, respectively, and both depend on temperature. It is clear that the result is also valid for each fluid in any temperature.Figure 1
Typical curve of molecular diameter versusT/TC for Ar at selected.


## 3. Internal Pressure

Whatever the model chosen for the liquid state, the cohesive forces are of primary importance. Ignoring, for the time being, the origin of these forces, we can estimate their magnitudes from thermodynamic internal pressure. In the case of an ideal gas, the internal pressure is zero because intermolecular forces are absent. In the case of imperfect gas, the internal pressure becomes appreciable, and in the case of a liquid, it may become much greater than the external pressure [8].The internal pressure is the resultant of the forces of attraction and the forces of repulsion between the molecules in a liquid. It therefore depends markedly on the volumeV and thus on the external pressure p. This effect is considered in the following for diethyl ether at 298 K. For moderate increases in p, the internal pressure decreases only slightly, but as p exceeds 500 MPa, the internal pressure begins to decrease rapidly and goes to large negative values as the liquid is further compressed. This behavior reflects on a larger scale the law of force between individual molecules; at high compressions; the repulsive forces become predominant.Joel Hildebrand was the first to point out the significance of the internal pressure of liquids in determining solubility relationships. A solution of two liquids that differ considerably in the internal pressure usually exhibits considerable positive deviation from ideality, that is, a tendency toward lowered mutual solubility [12–17].


## 4. Experimental Test

In examining the accuracy of this work, we wish to address the following specific questions. (i) For what types of fluids is the equation valid? (ii) Over what range of density does the equations hold for these isotherms?To see whether the regularity is limited to the certain type of fluid or is generally true, we have obtained the molecular diameter for a number of different fluids at different temperatures [18, 19]. The various fluids tested are listed in Table 1, together with the intercept, slope, pressure range of the experimental data, and the molecular diameter for a typical isotherm. We conclude that the equation is a property of many types of dense fluids. It turns out that all the fluid isotherms hold for ρ≥ρB mol L−1 that this is nearly equal to the Boyle density. The Boyle temperature TB and Boyle volume VB are defined in terms of the second virial coefficient B2 as follows: B2(TB)=0 and VB=TBB2′(TB), where B2′=dB2/dT.Table 1
Comparison between the calculatedσP<?brm?>int<?erm?> and experimental values of molecular diameters ση for different fluids.
Fluidρ/mol L−1T/KABΔp/MPaσPint,Åση,ÅaAr13.40742000.9343−0.095350–6003.82493.434bN211.1839702.8413−0.23190.06–303.91513.722cCO10.74002000.8460−0.088445–1004.13513.706dCO210.63854001.0995−0.101940–1004.06613.881eCH410.07202500.9904−0.104650–2004.23213.774eC2H66.76004001.05723−0.101030–704.75354.48eC3H84.93004001.2486−0.091020–1005.04815.061eC4H103.92305001.1896−0.082430–705.40165.339fC6H63.86206401.1561−0.074730–905.36735.270gC6H5CH33.16603005.8742−0.46480.1013–1005.9318932.5Referencesa[20], b[21], c[22], d[23], e[24], f[25], g[26].To investigate the liquid region, we turn first with fluid in its liquid range. (By liquid we meanT≺TC) Except near the critical point, all the isotherms are found to be quite linear over almost the entire density range from the vaporization line to the freezing line. The general conclusion, which is borne out for the other substances examined, is that the linear relation holds in the liquid region for ρ≥ρB. In summary, we can tell, the upper density limit for the linearity in liquid region is the freezing line.Thus, the general rule in both the liquid and supercritical regions is valid that the density holds forρ≥ρB.


## 5. Discussion

An important parameter related to the nature of fluids is the molecular diameterσ, given in terms of experimental quantities as second virial coefficient σB and viscosity coefficient ση. Although authors had used several microscopic arguments in a heuristic way in obtaining universal parameters, its actual basis amounts only to the use of virial expansions for equation of state and truncation of these expansions is the only approximation. As justification for the accuracy of the truncation, other than good agreement with experimental data. However, there is no direct indication of any upper density limit on the molecular parameters [15, 16].The molecular parameters are particularly useful in determination of molecular structure. Qualitatively, nature of molecule has been indicated by molecular diameter at given temperatures. In this paper, the internal pressure has been used to calculate molecular diameter. By use of the temperature dependence of internal pressure parameters, the molecular diameter has also been calculated at each temperature. In this paper, we have also shown that internal pressure is closely related to the molecular diameter.The cohesive forces, which are the resultant of forces of attraction and forces of repulsion between liquid molecules, holding a liquid together create a pressure within the liquid which has been termed the internal pressure.Internal pressure, fundamental liquid property, is closely related to the different properties of liquids such as ultrasonic velocity, free volume, viscosity, surface tension, solubility parameter, and latent heat of vaporization in the liquid phase.The precise meaning of the internal pressure is contained in a generalized manner in the following well-known thermodynamic equation. United forces of external and internal pressure equalize the thermal pressure which tries to expand the matter. If the internal pressure of a fluid is available, then the molecular diameter can be calculated easily. The main objective of the present research is to compute the molecular diameter.Therefore, the internal pressure was used to calculate the molecular diameter of quite different fluids for which the reliable experimental values are available. The calculated and experimental values of molecular diameter are found to present a close agreement with each other.By careful examination of the present results, it seems to us that this work contains point of particular interest; the physical interpretation in terms of molecular forces is as follows.If the shape of intermolecular potential function,U(r), for molecules in fluids is considered, it is seen that U(r) can be divided into two parts, attractive UA and repulsive UR parts, namely,(11)U(r)=UA+UR.
Differentiation at constant temperature with respect to volume gives(12)(∂U(r)/∂V)T=(∂UA/∂V)T+(∂UR/∂V)T,
namely,(13)Pi=Pi,A+Pi,R,
where attractive internal pressure, Pi,A is seen to be positive that it is related to the A term and repulsive internal pressure, Pi,R is negative, therefore, it is proportional to B term the internal pressure. By application of the regularity, we suggest that the molecular diameter is due to a balance between intermolecular attractive and repulsive forces [8].A simple model had been used to predict the temperature dependences of the molecular diameter and shown to agree with experiment. At the simplest level, we now have a new and useful way of molecular diameter on the isotherms of compressed dense fluids. Moreover, the temperature dependences of molecular diameter that characterize each linear isotherm are well described by the interceptA that is proportional to 1/T and the slope B is functional in 1/T.Finally, we expect the regularity to hold for dense mixtures as well as for single substances. The reason for this expectation is the wellknown success of so-called one-fluid approximations for mixtures. Such approximations assume that mixtures obey the same EOS as single substances, but with parameters that depend on the composition of the mixture. Since the present linearity has been shown to be consistent with an EOS, adjusting the parameters of the EOS for compositions will not alter this result. Testing of this expectation and investigation of the composition dependences of the parametersA and B remain for future work [7].



---

*Source: 521827-2012-03-15.xml*