Condensed Matter Physics, 2014, Vol. 17, No 4, 43001: 1–27 DOI: 10.5488/CMP.17.43001 http://www.icmp.lviv.ua/journal The phase transition of the ﬁrst order in the critical region of the gas–liquid system arXiv:1501.02325v1 [cond-mat.stat-mech] 10 Jan 2015 I.R. Yukhnovskii Institute for Condensed Matter Physics of the National Academy of Sciences of Ukraine, 1 Svientsitskii St., 79011 Lviv, Ukraine Received July 15, 2014, in ﬁnal form October 1, 2014 This is a summarising investigation of the events of the phase transition of the ﬁrst order that occur in the critical region below the liquid–gas critical point. The grand partition function has been completely integrated in the phase-space of the collective variables. The basic density measure is the quartic one. It has the form of the exponent function with the ﬁrst, second, third and fourth degree of the collective variables. The problem has been reduced to the Ising model in an external ﬁeld, the role of which is played by the generalised chemical potential µ∗ . The line µ∗ (η) = 0, where η is the density, is the line of the phase transition. We consider the isothermal compression of the gas till the point where the phase transition on the line µ∗ (η) = 0 is reached. When the path of the pressing reaches the line µ∗ = 0 in the gas medium, a droplet of liquid springs up. The work for its formation is obtained, the surface-tension energy is calculated. On the line µ∗ = 0 we have a two-phase system: the gas and the liquid (the droplet) one. The equality of the gas and of the liquid chemical potentials is proved. The process of pressing is going on. But the pressure inside the system has stopped, two ﬁxed densities have arisen: one for the gas-phase η g = η c (1 − d /2) and the other for the liquid-phase η l = η c (1 + d /2) (symmetrically to the rectlinear diameter), where η c = 0.13044 is the critical density. Starting from that moment the external pressure works as a latent work of pressure. Its value is obtained. As a result, the gas-phase disappears and the whole system turns into liquid. The jump of the density is equal to η c d , p where d = D/2G ∼ τν/2 . D and G are coeﬃcients of the Hamiltonian in the last cell connected with the renormalisation-group symmetry. The equation of state is written. Key words: critical point, phase transition of the ﬁrst order, grand partition function, quartic density measure, collective variables PACS: 05.70.Jk, 64.70.F-, 64.60.F- 1. Introduction The present research follows our previous consideration concerning the liquid–gas critical point. In this work we continue the investigations of the system of interacting particles in the vicinity of the critical point [1–4]. The great partition function was calculated. The potential of the interaction between particles consists of the sum of the van der Waals interaction and of the short-range interaction between the hardspheres. We used the collective variables (CV) method. The Jacobian of the transition from the Cartesiancoordinates space to the space of the collective variables was averaged with the statistical weight of a distribution function of the hard-spheres. The problem was reduced to the Ising model in an external ﬁeld, presented by the generalised chemical potential µ∗ . As a basic density measure, the quartic distribution function is used. It has the form of the exponent function with the ﬁrst, second, third and fourth degree of the collective variables. Detailed calculations were performed in [2] for the case of τ Ê 0, where τ = (T − Tc )/Tc . Therein, the coordinates of the critical point Tc and η c were obtained, where Tc is the critical temperature, η c is the critical density η c = (N /V )c (πσ3 /6), σ is the diameter of a hard sphere. © I.R. Yukhnovskii, 2014 43001-1 I.R. Yukhnovskii Table 1. The critical indexes for the quartic and sextic density measures in the CV method [5]. Critical parameters ν α β γ Approximation ρ 4 Approximation ρ 6 0.605 0.637 0.185 0.088 0.303 0.319 1.210 1.275 The critical indexes correspond to that of the Ising model. Their values, the sizes of the particles, the critical temperatures and densities are given in tables 1 and 2 [1]. The ﬁrst results of the investigations in the area of temperatures below the critical point, τ É 0, are given in [1–4]. It has been shown that the events connected with the phase transition of the ﬁrst order from the gas phase to the liquid phase (and inversely) are described by the integral over the variable ρ 0 for k = 0. This integral embraces the results of the integration over all collective variables with the exception of the variable ρ 0 . The collective variable ρ 0 is connected with the distribution of the number of particles. The integration over the variable ρ 0 was performed in the steepest descent method. It is shown that the discriminant Q of the equation concerning the maximum of integrand can take positive or negative values, or be equal to zero. When the discriminant Q is positive we have pure monophasic states of the gas or of the liquid. In the area Q < 0, the maximization problem possesses three real roots, ρ max = b cos 0 ϕ + 2πn , 3 n = 0, 1, 2, 0 É ϕ É π. However, only two of them, actually for n = 0 and n = 1, correspond to a maximum of the integrand. The areas of the principal and of the secondary maxima were determined. It was shown that there is a break between the end of the principal maximum for the gas phase and the beginning of the principal maximum for the liquid phase. Those were the main results of the papers [1, 2]. In this work we would like to pass the whole route of the phase transition of the ﬁrst order. The work consists of the introduction and seven sections. In Introduction we deﬁne the problem and give the main results of the mentioned works. In section 2 the description of the monophasic states is given, which corresponds to the positive values of the discriminant Q . Section 3 is devoted to the coordinate form of the reference system. We work within the framework of the quartic density measure. The reference system is included in our problem through the cumulants M1 , M2 , M3 and M4 . In accordance with the accuracy of the consideration of the whole problem, it is needed to take a suitable expression for the pressure P 0 , the chemical potential µ0 and the free energy F 0 of the reference system. Section 4 is an essential one in this work. The domain of the phase transition is considered. Here, the discriminant is negative Q < 0. In this region we have two maxima — the principal and the secondary one. The principal maxima correspond to the pure gas or to the pure liquid phases. The secondary maxima are connected with the stable density ﬂuctuations, which are either of the liquid type in the gas phase, or of the gas type in the liquid phase. Table 2. The critical temperatures Tc and the effective hard sphere diameter σ for some systems. System CO-CO Ar-Ar Kr-Kr Xe-Xe O2 -O2 N2 -N2 43001-2 P c , atm 34.532 47.964 54.182 57.537 49.77 33.55 Tc , °C (exp.) Tc , °C σ/σ0 ε/k B , K [1] σ0 , Å (exp.) σ, Å [1] 37.92 51.75 56.76 65.88 54.76 44.08 –140.23 –122.65 –63.1 16.62 –118.84 –147.05 –138.46 –123.27 –67.84 16.84 –110.8 –150.02 3.76 3.405 3.6 4.1 3.58 3.698 3.37 3.14 3.367 3.71 3.18 3.365 0.898 0.922 0.935 0.905 0.89 0.91 100.2 119.8 171 221 117.5 95.05 P c , atm Phase transition of the ﬁrst order Under consideration is the equation for the envelope curve of two maxima, the principal and the secondary ones: PV − ln Ξρ 0 = 0, Θ ∂ ln Ξρ 0 ∂βµ = N. Having opened the second equation we get the probabilities of the system to be in the states corresponding to the principal or to the secondary maxima. We have got analytical expressions for the probabilities to be in the gas or in the liquid phases. The central considerations concern the events along the line µ∗ = 0. Here, when the gas system is pressed, a droplet of liquid appears inside it, that is, a new phase arises — the phase of liquid. In the appearance of the droplet lies the essence of the phase transition of the ﬁrst order (as a result of isothermal pressure on the gas system). The work for the formation of a droplet — the surface-tension energy — was calculated. A fundamental quantity ∆ is introduced. It describes an inverse to the usual order parameter, characterizing the weight of the pure monophase state in the system. On the rectangle vertex of the Q = 0, the weight is equal to unity for the gas phase and zero for the liquid phase or, on the contrary, dependent on which of them the maternal is. The quantity ∆ is equal to zero on the line of the phase transition. This means that on that line the gas phase (maternal) and the liquid phase (in a droplet-form) have equal statistical weights. In its initial determination ∆ was a function of the reference system cumulants M2 , M3 , M4 and, therefore, was a function of the compressibility coeﬃcients and its derivatives. The equality ∆ = 0 means that on the line µ∗ = 0, the compressibility coeﬃcient equals zero. We have some metastable state — the state of phase transition of the ﬁrst order. In section 6, the equality of chemical potential of the gas phase and of the liquid phase on the line µ∗ = 0 is shown. In section 7, the gas-condensation process is analyzed. An analytical expression for the surface-tension energy is obtained. The latent work A lat in the isothermal pressure process (an analogue to the latent heat in an isobaric process), which acts along the line µ∗ = 0 during the gas-liquid phase transition, was obtained. And, in the end, in section 8, the equation of state was obtained. The equation for the isotherm P = P (η)τ=const has three parts: 1) a monophasic gas branch: P = P (η, τ) for densities η É η g = η c (1 − d /2), d = coeﬃcients in Hamiltonian at the second and at the fourth degree of ρ 0 ); p D/2G (D and G are 2) a horizontal line µ∗ (η) = 0 for the two-phase system, here P = 0, and η g É η É η l , η l = η c (1 + d /2). The transition of gas into liquid is takes place at the expense of the latent work A lat. acting from without; 3) the liquid-branch for η > η l = η c (1 + d /2). Thus, during the gas-condensation process we get a jump of the density of the amount of η l − η g = η c d ∼ τν/2 . As we see, the magnitudes of the densities of the beginning and of the end of the transition process from the gas state to the liquid state are located symmetrically to the rectilinear diameter. When the system approaches to the critical point, the speciﬁc surface-tension energy disappears as τ5/2ν , because D ∼ τ2ν , G ∼ τν , and the jump of the density η is proportional to τν/2 . In the preamble we stressed that within the framework of the Gibbsian statistics one can describe the phase transition of the ﬁrst order [6]. As for other works, we beneﬁted from the monographs of R. Balescu [7], L.D. Landau and E.M. Lifshits [8], J.-P. Hansen and I.R. McDonald [9] as well as from [10–12]. In this research we made a number of assumptions, especially concerning the reference system, the dependence of the coeﬃcients D and G on the density and others. In fact, by means of the method of collective variables it is possible to describe the critical point theory of the liquid-gas system with a more suﬃcient accuracy. Now we get to the theme. 43001-3 I.R. Yukhnovskii 1.1. How it should look like from the very start 3 Let a big volume V0 in the form of a sphere with the radius R 0 , V0 = 4π 3 R 0 and with the number of particles N , N = n 0 V0 be picked out in a gas medium with the density n 0 , pressure P at the temperature T . We consider the phase transition of the ﬁrst order as a result of the quasistatic pressure on the gas sphere. The process takes place at temperatures close to the critical point Tc . The dimensionless temperature τ is introduced τ = (T − Tc )/Tc , and |τ| É 0, 01, τ < 0. The pressure inside the sphere of radius R 0 is equal to P +δP , where δP is the additional pressure, connected with the spherical form of the volume V . The gas energy inside the sphere: E 0 = (P + δP )V0 − αS 0 , (1.1) where S 0 is the area of the sphere, S 0 = 4πR 02 , α is the surface-tension coeﬃcient. The minimum of the difference E 0 − PV0 means that δ(E 0 − PV0 )/δR = 0 or δP 0 = 2α . R0 (1.2) Such a picture may be useful to remember when discussing the isothermal phase transition in the gasliquid system illustrated in this work. We consider the grand partition function Ξ= ∞ zN X ZN , N =0 N ! (1.3) where N is the number of particles, z N is the activity N "µ z = mk B T 2π ¶3/2 1 ħ3 #N exp(βµN ), (1.4) m is the mass of a particle, k B is Bolzman’s constant, T is temperature, ħ is Plank’s constant, β = (k B T )−1 , µ is the chemical potential, Z N is the conﬁguration integral of N particles [1, 2]. The collective-variables method is used {ρ k } [3–5] to calculate Ξ. The Jacobian of transition from the Cartesian coordinates to the collective variables was calculated with the statistical weight of the system of hard cores. The latter was taken as a reference system. As a result, for Ξ we get the expression [1, 2]: Ξ = Ξ0 Ξ1 Ξp 0 , (1.5) where Ξ0 is the partition function of the reference system (RS), ln Ξ0 = βµ0 〈N 〉 − βF 0 = p 0 V0 , Θ (1.6) µ0 is the chemical potential of RS, F 0 is free energy of RS. The quantity Ξ1 is the partition function [1], which includes all interaction effects, with the exception of the effects connected with the formation of a new phase, i.e., the liquid phase within the maternal sphere of gas in the volume V0 : £ ¤ Ξ1 = Ξg exp −β (F CR + F IGR ) . (1.7) Here, Ξg is the partition function in which the effects appearing after integration over the collective variables ρ k for k > B are taken into account, where B is the point of the ﬁrst zero of the Fourier transform of the attraction potential; the integration is performed with the Gaussian density measure; F CR + F IGR are the free energies arising after integration over ρ k with 0 < k É B : F CR in the critical regime, where the renormalization group (RG) symmetry between the block-Hamiltonians exists and 2πn τ /a nτ É k É B , n τ is the number and a nτ is the periods of the block-lattice, on which, given τ, the renorm-group symmetry between the coeﬃcients of the block-Hamiltonians terminates; the integration is performed in the quartic density measure [3–5]; F IGR appears as a result of the integration over ρ k performed here in an inverseGaussian regime [3] for 0 < k É 2πn τ /a nτ . We suppose that all expressions in Ξ1 are known, and that all of them are ﬁxed quantities within the temperature region τ < 0.01. In this work, we consider only one 43001-4 Phase transition of the ﬁrst order part, Ξρ 0 , of the Ξ. It is connected with the density ﬂuctuations of the number of particles inside of the sphere of the radius R 0 . The phenomena of the phase transition of the ﬁrst order at τ < 0 are described by the expression Ξρ 0 . The chemical potential µ of the whole system is only within Ξρ 0 . As a result of calculations carried out in [1, 2], the initial form of the partition function Ξρ 0 for integration over ρ 0 is as follows: h iZ Ξρ 0 = exp µ∗ (1 − ∆)N exp[N E (ρ 0 )]dρ 0 , (1.8) where µ∗ is the generalized chemical potential of the system1 , µ∗ = β(µ − µ0 ) + ξ − |α(0)|(1 − ∆), (1.9) µ0 is the chemical potential of the RS, M3 , | M4 | ξ= (1.10) µ ¶ 1 ∆ = − ξ M2 + M3 ξ 2 , 3 (1.11) M2 , M3 , M4 are the cumulants of the Jacobian of transition [3, 4]; at small densities, ∆ is a linear function of density η, see ﬁgure 2. ¶ ˜ N Φ(k) α(0) = , V Θ µ (1.12) ˜ where Φ(k) is the Fourier transform of the attraction potential; E (ρ 0 ) = µ∗ ρ 0 + Dρ 20 −Gρ 40 , (1.13) G = G 0 |τ|ν (1.14) D = D 0 |τ|2ν , are coeﬃcients of the limiting block-structure Hamiltonian Hnτ at the variable ρ 0 , D 0 ' 1.19; G 0 = 1.67; ν = ln s ∗ / ln E 1 = 0.605 is the correlation length critical exponent, E 1 is the greater of two eigenvalues of the matrix for linearized recurrent relations of Wilson’s type [3–5, 13–16] for coeﬃcients of blockstructure Hamiltonians, s ∗ is the optimal parameter in the scale transformations s ∗ = 3.58 . . ., (s ∗ ) < E 1 < (s ∗ )2 . We work in a critical region of temperatures and densities2 , τ < τ∗ , η < η∗ . For the limit interval of the temperature τ∗ we have the expression: [4] " C˜1 τ∗ ˜ r ∗ + βΦ(0)η #ν 6 πσ3 " = 1− a2 6 ˜ βΦ(0)η πσ3 #1/2 . (1.15) Here, r ∗ is one of the ﬁxed-point coordinates, a 2 is a coeﬃcient in the Jacobian of transition. On both sides of (1.15) there are inverse radia of correlation in the quartic and in Gaussian density measure, correspondingly [4]. For the limiting point of the interval of densities we have: s ∗ |η − η c | = η c D 0 ∗ ν/2 (τ ) , G0 (1.16) 3 N πσ ∗ where η = V 6 , σ is a diameter of a particle (hard core), η is the limiting value of density for the critical regime. In general, coeﬃcients D 0 and G 0 in (1.13) depend on density. In this work we take them as η = η c , where η c = 0.13044 is the critical density determined from (1.11) by the condition ∆ = 0. The Hamiltonian E (ρ 0 ) in (1.8) is obtained as a result of integration in (1.5) over all variables ρ k , with the exception of ρ 0 3 . 1 While investigating the characteristics of the integrand, the origin of the coordinates is transferred to the point ∆ = 0, µ∗ = 0, and thus the change of ρ 0 may accept both positive and negative values (see ﬁgure 1). 2 “critical” means that within this region there exists the RG symmetry. 3 The variable ρ is introduced in Ξ in (1.5) as a result of integration over d Γ and summation over N of the expression δ(ρ − 0 0 N p N / 〈N 〉), where N is the independent variable of the number of particles. In order to integrate over ρ 0 in (1.8) we surely have to take the grand partition function. 43001-5 I.R. Yukhnovskii The chemical potential µ appears only in Ξρ 0 . Thus, the variable ρ 0 is connected with the density and with the distribution of the number of particles. By deﬁnition ln Ξ = PV Θ ∂ ln Ξ ∂ ln Ξρ 0 = = 〈N 〉. ∂(βµ) ∂µ∗ and (1.17) Based on these two expressions we get the equation of state: P = p(τ, η), (1.18) connecting such variables as pressure, temperature, density, p , τ, η. The curve of equations (1.17) is an envelope curve for pV /θ = ln Ξ, the latter depends on p , τ, η, µ. We work within the framework of the equilibrium thermodynamics. The number N is large. We use the presence of the thermodynamic limit: N → ∞, N = const. V V → ∞; (1.19) In this connection, when we integrate over ρ 0 in (1.8) we use the steepest descent method. The integration in (1.8) was carried out and discussed in [1, 2]. We give here only the ﬁnal conclusions. R In the integral exp[N E (ρ 0 )]dρ 0 in (1.8) we looked for the maximum of the subintegral function. The derivative ∂E (ρ 0 ) = µ∗ + 2Dρ 0 − 4Gρ 30 = 0 ∂ρ 0 (1.20) gives us the cubic equations: ρ 30 + V ρ 0 + W = 0, on condition of E¨ (ρ 0 ) < 0: V =− where 2D − 12Gρ 20 < 0 1D , 2G W =− µ∗ , 4G 2 3(ρ max 0 ) + V > 0. or (1.21) (1.22) The equation (1.21) has three roots, depending on the sign of the discriminant Q : Q= W2 V3 + . 4 27 (1.23) Here, the ﬁrst addend is positive, and the second is always negative. Therefore, the discriminant Q of equation (1.21) may be positive, equal to zero and negative. The case of Q = 0 separates the roots. At p p 3 3 ∗ /G Q = 0 (W /2)2 = − (V /3)3 , we have three real roots, among them only the root ρ max = 2 −W /2 = µ 0 meets the requirement (1.22)4 . From equation (1.23) at Q = 0 we get Ãs µ∗ = ±G 2D 3G !3 s ρ max =± 0 and thus 2D . 3G (1.24) The roots (1.24) are the vertex of the rectangle on the plane µ∗ , ρ 0 . Let us introduce the notation s 3 a = Gb , b= 2D . 3G (1.25) Out of the rectangle we have the region Q > 0, inside Q < 0. In the case of Q > 0, equation (1.20), or (1.21), has only one root [1] ρ max = ρ1 = q 0 4 Two other roots, 43001-6 ρ 2 and ρ 3 = − 21 µ µ∗ G ¶1/3 p 3 µ∗ /G do not obey (1.22). , µ∗ = ±G µ ρ1 q ¶3 , (1.26) Phase transition of the ﬁrst order m * , 1 0 - 5 Q > 0 2 Q < 0 -b b /2 -d d b r m a x 0 -0 ,0 5 -b /2 -0 ,1 0 0 ,0 0 0 ,0 5 0 ,1 0 -2 Figure 1. The curve of the generalized chemical potential µ∗ as functions of ρ max 0 . On this curve, the integrand exp[N E (ρ)] in the region of (1.8) has maxima [solutions (1.26) for Q > 0 and (1.27) for Q < 0]. ¡p p p Here, a = G 2D/3G ∼ τ5/2ν ; b = 2D/3G ∼ τν/2 ; d = D/2G ∼ τν/2 . The thick lines correspond to principal maxima of the function E (µ∗ , ρ 0 ), thin lines correspond to secondary maxima, and the very thin line is a line of minimum. The origin of the system’s coordinates is located in the point µ∗ = 0, ∆ = 0. The abscissa axis is along the line µ∗ = 0. The ordinate axis is along the rectilinear diameter ∆ = 0. ¢3 where ( q =4 −1/3 µ |γ| 1+ 4 ) ¶1/3 +... , q = 4−1/3 at γ = 0, q =1 at γ = 1, µ ¶ µ ¶ V 3. W 2 γ= − < 1, 3 2 which corresponds to (1.24) at Q = 0. In the case of Q < 0, we have three real roots: µ∗ = a cos γ, n = 0, 1, 2. Only the roots and ϕ 3 ρ 1,2,3 = b cos γ + 2kπ . 3 (1.27) ϕ + 2π (1.28) 3 meet the requiremets (1.22). The principal maximum of exp[N E (ρ 0 )] for the root ρ 02 takes place for π π 2 É ϕ É π (the thick line in ﬁgure 1), and for the root ρ 01 , when 0 É ϕ É 2 . The secondary maxima of π π exp[N E (ρ 0 )] exist for ρ 02 in the region 0 É ϕ < 2 and for ρ 01 in the region 2 É ϕ É π. The curve for the principal maxima of the generalized chemical potential µ∗ has a jump in ρ 0 from −d to +d on the axis µ∗ = 0. A smooth transition of the solution ρ 1 into ρ 02 at the point (−b, −a) and ρ 1 into ρ 01 at the point (b, a) ρ 01 = b cos and ρ 02 = b cos takes place. Herein below we shall consider an isothermal process of a quasistatic pressure on the gas system at τ = const, τ < 0, |τ| < τ∗ . 2. Monophasic states of liquid and gas The solution ρ 0max = ρ 1 in (1.26), describes the values of densities of absolute maxima for the function exp[N E (ρ)] in integral (1.8) for Ξρ 0 . Considering N as a very big number, according to the steepestdescent method, we write: ln Ξρ 0 = N (1 − ∆)µ∗ + N E (ρ 1 ). (2.1) 43001-7 I.R. Yukhnovskii Eliminating µ∗ by means of (1.17) ∂ ln Ξρ 0 ∂µ∗ =N (2.2) one gets (1 − ∆)N + N ρ 1 = N , and ρ 1 = ∆. (2.3) Then, in correspondence with (1.20) µ∗ = µ∗ (∆) = −2D∆ + 4G∆3 ∗ 2 4 E (ρ max 0 ) = µ ∆ + D∆ −G∆ . and (2.4) Relationships (1.20)–(2.2) together with (1.5)–(1.8) give us the equation of state PV = P 0 V + P 1 V + P ρ 0 V, (2.5) P 0 V = Θ ln Ξ0 (2.6) P 1 V = Θ ln Ξ1 (2.7) where — contribution from the reference system, is the contribution from the region k Ê B and from the integrals over ρ k for 0 < k É B . The latter includes integrals in the critical regime and in the inverse Gaussian regime, and the free energy of the RS up to the cumulant M4 . Finally, P ρ 0 V = Θ ln Ξρ 0 , (2.8) ln Ξρ 0 = N E (∆), (2.9) where in accordance with (2.1)–(2.4) E (∆) = µ∗ (∆) + D∆2 −G∆4 , (2.10) 3 ∗ µ (∆) = −2D∆ + 4G∆ , ¡ ¢ P ρ 0 V = ΘN −2D∆ + 4G∆3 + D∆2 −G∆4 . (2.11) (2.12) The sum of (2.6), (2.7) and (2.12) in (2.5) presents the gas-state isotherm when ∆ É −b , and the liquid-state isotherm when ∆ Ê b . Our main task is to consider events in the region |∆| É b . In this region, the discriminant (1.23) of equation (1.21) is negative. Note that we restrict ourselves to τ < 0, |τ| < τ∗ ' 0, 02 and s |η − η c | < η c D 0 ∗ ν/2 (τ ) . G0 (2.13) By deﬁnition, the coordinates of the critical point are τ = 0, ∆(η c ) = 0, η c = 0.13044. The expression for ∆ is given in (1.11). The quantity ∆ is also determined in (2.3). Here, for ρ 1 we have its value (1.26), the quantities D and G being connected with the attractive forces. Thus, owing to expression (2.2), both the connection of the effects of the long-range and of the short-range (of the RS) interactions takes place. In the critical region, we can also use for ∆ the following expansion: µ d∆ ∆(η) = ∆(η c ) − dη ¡ ¶ ηc (η − η c ). (2.14) ¢ Here, the following expressions hold d∆/dη η=η = 1/η c . According to the deﬁnition of the critical point c ∆|η=ηc = 0, and for (2.14) we get (see ﬁgure 2): ∆= 43001-8 1 (η − η c ) ηc (2.15) Phase transition of the ﬁrst order 0 ,1 5 D D 0 ,1 0 D lin 0 ,0 5 h 0 ,0 0 0 ,1 1 5 0 ,1 2 0 0 ,1 2 5 0 ,1 3 0 0 ,1 3 5 0 ,1 4 0 0 ,1 4 5 0 ,1 5 0 -0 ,0 5 -0 ,1 0 -0 ,1 5 Figure 2. The region of the linear approximation for ∆ as a function of the density and vice versa. Our discussion concerns the critical region of densities [see equation (1.16)] |η − η c | É 0.02 ∆ = −(M2 ξ + 13 M3 ξ2 ), η = (N /V )(πσ3 /6), η c = 0.13044. and η = ηc ∆ + ηc . (2.16) It follows that for gaseous state ρ 1 É −b and for liquid state ρ 1 Ê b in (2.15), one has some restrictions on values of ∆ and η. Besides, we should take into account the restrictions on density in (2.13), as well as those connected with cumulants M4 [2, 5]. For µ∗ we have expressions (2.4) and (1.9). The latter includes the chemical potential of the reference system, and (2.5) has a contribution into the pressure of the whole system. 3. On the reference system Let us turn back to formula (1.5) for the grand partition function Ξ. The logarithm of Ξ gives the whole equation of state accounting for external shell, thus accordingt to (1.1) ln Ξ = P sph Vsph − αS sph , where according to (1.2), P sph = p + δP 0 . In this section, we omit the terms connected with the existence of the external shell assuming ln Ξ = PV /θ . In accordance with (1.5) ln Ξ = ln Ξ0 + ln Ξ1 + ln Ξρ 0 . (3.1) Our aim is to separate in the right-hand side the terms belonging only to the RS. At ﬁrst it is ln Ξ0 = P 0 V /θ , where P 0 is the pressure of the RS. Then, from (2.9), (2.10) and (1.9) for ln Ξρ 0 we may write: h i ln Ξρ 0 = N β(µ − µ0 ) + ξ − |α(0)|(1 − ∆) + N (D∆2 −G∆4 ). (3.2) Thus, in the right-hand side we have two addends connected with RS, namely, N (ξ−βµ0 ). They do not disappear, when excluding the attractive interaction, taking α = 0. Of course, in the right-hand side of (3.1) we should exclude some additional expressions, completing together the free energy F 0 of the RS. We are talking about the factor5 lnQ −1 (−ξ, a 2 , a 4 ) that forms a part in the expression for Ξ1 . Collecting together all terms remaining in the expression (3.1) at α = 0, we can write (3.1) in the form: PV P 0 V − = N {µ − µ0 } − (F − F 0 ), Θ Θ 5 It is shown in [2], Q(−ξ, a 2 , a 4 ) = (3.3) ³ ´ iN B a2 2 a4 4 , where a 2 , a 4 are the coeﬃcients in the transition −∞ exp −ξρ l − 2 ρ l − 4! ρ l dρ l hR ∞ Jacobian into the phase space of collective variables. 43001-9 I.R. Yukhnovskii where P 0 V /Θ = ln Ξ0 , and F − F0 = h ¡ ¢i N −|α(0)|(1 − ∆) + N D∆2 −G∆4 + · ¸ NB + ln Ξg − β(F CR + F IGR ) − N −ξ + lnQ(ξ, a 2 , a 4 ) . N (3.4) As a result, we may write F 0 in the form F0 = N ξ − NB lnQ(ξ, a 2 , a 4 ). Θ (3.5) Herein below in our computations, in the critical region we shall consider the quantities entering (3.4) as the known quantities. We have presented here the expressions for F 0 because we say that the pressure P 0 , and the chemical potential µ0 ought to be determined from the expression for F 0 . The latter is obtained here using the quartic density measure while calculating the Jacobian. Thus, for the chemical potential of the RS we shall take the expression µ0 = ∂F 0 ' ξ. ∂N (3.6) The extraction of the RS while calculating the equation of state according to (3.3), will not be considered here, and will be the subject of another research. 4. The ﬁeld of the phase transition We consider the partition function (1.5) in the case of negative values of the discriminant (1.23) of equation (1.20). The part of the expression (1.5) will be analysed in detail, which is connected with the chemical potential via Ξρ 0 according to (1.17). For Ξρ 0 , we have expressions (1.8)–(1.13). We shall remember here the main form h i Z∞ h i ∗ Ξρ 0 = exp N µ (1 − ∆) exp N E (ρ) dρ. (4.1) −∞ In the region Q < 0, equation (1.20) for extremum of exp N E (ρ) has the roots of (1.27). We shall investigate the equation (1.17) for the envelope curve in the case of Q < 0. Let us consider more in detail the second of the two equations ∂ ln Ξρ 0 ∂µ∗ − N = 0. (4.2) Taking into account (4.1), one has R∞ N ρ exp(N E (ρ))dρ −∞ N (1 − ∆) + ∞ R = N, exp(N E (ρ))dρ −∞ or R∞ ρ exp[N E (ρ)]dρ −∞ R∞ = ∆. (4.3) exp[N E (ρ)]dρ −∞ On the right we have ∆ = ∆(η), given in (1.11), where the cumulants of the RS are present. On the lefthand side of equation (4.3) there are the magnitudes connected with the long-range effects of the Van der Waals forces. This important equation connects together the long-range effects accumulated in the left-hand side, and the short-range effects accumulated in ∆. As we remember, long-range forces are described in the collective variables phase space and the short-range ones in the Cartesian-phase space. 43001-10 Phase transition of the ﬁrst order R∞ We devide the integral −∞ ρ exp[N E (ρ)]dρ into two parts Z∞ Z0 Z∞ ρ exp[N E (ρ)]dρ + ρ exp[N E (ρ)]dρ. −∞ 0 ρ exp[N E (ρ)]dρ = −∞ In the ﬁrst one, along with ﬁgure 1, the maximum of the subintegral function exp[N E (ρ)] takes place for ϕ+2π ρ = ρ 02 = b cos 3 , 0 É ϕ É π. This is the domain of the gas state of the matter. In the second term, the ϕ maximum for exp[N E (ρ)] takes place for ρ = ρ 01 = b cos 3 — this is the liquid domain. In the whole integral (4.3), the principal maximum for the gas-phase trips together with the secondary maximum for liquid, and vice-versa the principal maximum for the liquid phase trips together with the secondary maximum for the gas. Let us go to a more detailed consideration of expressions (4.3). So, we carry out the integrals in (4.3), and denote by I= Z∞ ρ exp[N E (ρ)]dρ, K = and by −∞ Z∞ exp[N E (ρ)]dρ. (4.4) −∞ Taking into account that the maxima of the subintegral function go along the lines ρ = ρ 02 and ρ = ρ 01 and that the curve ρ 02 , where −b É ρ 02 É −b/2, is connected with the gas-phase, and the curve ρ 01 , where b/2 É ρ 01 É b , is connected with the liquid phase, we shall write down integrals (4.4) as the sums: I = Ig + Il K = Kg + Kl, and (4.5) where Z0 Ig = Z∞ I l = ρ exp[N E (ρ)]dρ, ρ exp[N E (ρ)]dρ, −∞ Kg = (4.6) 0 Z0 Kl = exp[N E (ρ)]dρ, −∞ Z∞ exp[N E (ρ)]dρ. (4.7) 0 All of them are of the same type, and we shall calculate I g precisely. We expand the function E (ρ) at the ϕ+2π 3 . maximum point ρ = ρ 02 = b cos Then, I g = exp[N E (ρ 02 )](I11 ρ 02 + I12 ), (4.8) where I11 Z0 = −∞ = I12 · ¸ N exp − |E¨ (ρ 02 )|(ρ − ρ 02 )2 dρ = 2 pπ · N 2 Z0 = −∞ 2 + 0 ¸ N exp − |E¨ (ρ 02 )|x 2 dx , 2 · · ¸ N 2 ¨ (ρ − ρ 02 ) exp − |E (ρ 02 )|(ρ − ρ 02 ) dρ 02 = 2 1 = |E¨ (ρ 02 )| |ρ Z 02 | ¸−1/2 |N E¨ (ρ 02 )| 1 exp − N |E¨ (ρ 02 )|ρ 202 . 2 For the second integral in I11 , we introduce the probability integral and get6 |ρ Z 02 | 0 6 While calculating p · ¸ N π 1 exp − E¨ (ρ 02 )x 2 dx = (erfc z) q , 2 2 N ¨ |E (ρ 02 )| 2 h i I11 , I12 and K we used exp N E (ρ 02 ) − 12 |E¨ (ρ 02 )ρ 202 ' exp[N E (0)], E (0) = 0. 43001-11 I.R. Yukhnovskii r z= · ¸ ∞ X 1 1 −z 2 m 1.3 . . . (2m − 1) 1+ (−1) . erfc z = p e (2z 2 )m πz m=1 N |E¨ (ρ 02 )||ρ 02 | , 2 Substituting into (4.8) and gathering together all the terms after shortening in the limit of big N , we ﬁnd I g = ρ 02 eN E (ρ 02 ) q p π . (4.9) N ¨ 2 |E (ρ 02 )| The integral I l when calculated undergoes the same transformations and we obtain I l = ρ 01 eN E (ρ 01 ) q p π . (4.10) N ¨ 2 |E (ρ 01 )| The integrals K g and K l are actually calculated while calculating I11 . For K g and K l , we have: K g = exp[N E (ρ 02 )] q K l = exp[N E (ρ 01 )] q p π , N ¨ 2 |E (ρ 02 )| p π . (4.11) N ¨ 2 |E (ρ 01 )| After all these calculations, we return to expressions (4.2), (4.3), that together with (4.1) describe the envelope-curve equation, or, in other words, to the equation of state in the region of Q < 0. Expression (4.3) now has the form: ∆= Ig + Il Kg + Kl , ∆ = ρ 02 w g + ρ 01 w l , or explicitly (4.12) where wg = wl = c g exp[N E (ρ 02 )] c g exp[N E (ρ 02 )] + c l exp[N E (ρ 01 )] = Kg Kg + Kl , c l exp[N E (ρ 01 )] Kl = , c g exp[N E (ρ 02 )] + c l exp[N E (ρ 01 )] K g + K l c g = (|E¨ (ρ 02 )|)−1/2 , c l = |E¨ (ρ 01 )|−1/2 . Here, the functions w g and w l play the role of the probabilities of the gas or of the liquid state, accordingly. wl = 1 1 2 0 at at at ϕ = 0, ϕ = π2 , ϕ = π; wg = 0 1 2 1 at at at ϕ = 0, ϕ = π2 , ϕ = π. As we see, w g + w l = 1. (4.13) Herein, −b, ϕ = π, ϕ + 2π −d , ϕ = π2 , ρ 02 = b cos = b 3 − 2 , ϕ = 0. b , ϕ = π, ϕ 2 ρ 01 = b cos = d , ϕ = π2 , 3 b, ϕ = 0. (4.14) (4.15) For ∆, we have: ∆ = ∆g + ∆l , 43001-12 (4.16) Phase transition of the ﬁrst order Figure 3. Probabilities w l and w g for single-phase systems, of liquid (w l ) and of gas (w g ), respectively, as functions of the angle ϕ. Two bottom axes are to represent values ρ 01 (ϕ) and ρ 02 (ϕ). −b, ϕ = π, −d /2, ϕ = π2 , 0, ϕ = 0. ϕ = π, 0, d /2, ϕ = π2 , ∆l = w l ρ 01 = b, ϕ = 0. ∆g = w g ρ 02 = q (4.17) (4.18) q 3D 1D ν/2 ν/2 Here, b = ;d= . 2 G ∼τ 2 G ∼τ In formulas (4.13), (4.17) and (4.18) we give important points for functions w g , w l and ∆g , ∆l . In fact, we should be driven by expressions (4.12) or even by more precise solutions of equation (4.3). In ﬁgure 3, the graphics of w g and w l are presented as functions of the angle ϕ. Attention should be paid to the fact that the function ∆ as a sum ∆ = w g ρ 02 + w l ρ 01 = ∆g + ∆l turns into zero at ∆g = −d /2 and ∆l = d /2. In ﬁgure 4, the line (−d /2 ÷ d /2) is interpreted as the line of the gasliquid, ﬁrst-order phase transition. The points ∆g |ϕ=π/2 = −d /2 and ∆l |ϕ=π/2 = d /2 do not lie on the curve b 0 ,1 0 0 D r 0 2 w G 0 ,0 7 5 d r 0 1 w L 0 ,0 5 0 j G w 0 2 N = 1 0 7 d /2 0 ,0 0 0 0 1 D , r p p /2 w L , r 0 ,0 2 5 - 3 t = - 1 0 -0 ,0 2 5 -0 ,0 5 0 -d /2 -d -0 ,0 7 5 -0 ,1 0 0 -b ϕ+2π ϕ Figure 4. Plot for ∆ = ∆g + ∆l and ∆g = w g ρ 02 ; ∆l = w l ρ 01 ; ρ 02 = b cos 3 ; ρ 01 = b cos 3 ; 0 É ϕ É π. The vertical fragment d /2 ÷ −d /2 joins the point of the largest densities of gas with that of the smallest density of liquid in the two-phase system. 43001-13 I.R. Yukhnovskii 0 ,0 0 0 5 -ffL -b -0 ,2 r fL -d -fG 0 ,0 0 0 0 -0 ,1 0 ,0 0 ,1 d ffG b m a x 0 ,2 m * E e t = - 1 0 - 2 -0 ,0 0 0 5 ϕ+2πn Figure 5. The quantities E , E and µ∗ as functions of ρ max = b cos , n = 1, 2, 3 within the region 3 Q < 0 [see equations (1.23), (1.27) and (1.28)]. Thick lines represent the principal maxima for E (ρ). Thin lines represent the secondary maxima. For gaseous branch n = 1, the points −b, −d , −b/2 correspond to values ϕ = 0, π/2, π; for liquid branch n = 0, the points b, d and b/2 correspond to the values ϕ = 0, π/2, π, respectively. ∆(ϕ)|ϕ=π/2 . The sum of their values corresponds to the point ∆(ϕ) = 0 for ϕ = π2 . The expression ∆g = −d /2 consists of two factors: ρ 02 = −d is the end of the principal maxima of the function exp[N E (ρ 0 )] located upon the gaseous branch ρ 02 (ϕ), and w g = 1/2 that characterises the probability of arising of the value ρ 02 . For the point ∆l = d /2 to w l = 1/2 and the value ρ 01 = d , which starts the curve £ , it corresponds ¤ of the principal maxima of exp N E (ρ 0 ) located upon the liquid branch ρ 01 (ϕ). Note, that the density of either gaseous or liquid phase is independent of either ρ 01 or ρ 02 , but depends on ∆, according to (4.16)– (4.18). A jump of the density at the phase transition is connected with the points ρ 01 = d and ρ 02 = −d , given in ﬁgure 1 for µ∗ = µ∗ (ρ 0 ), but its magnitude depends on the values of products w g ρ 02 and w l ρ 01 in these points, in other words, on the probabilities w g (−d ) = 12 and w l (−d ) = 21 as well. Thus, we have η g |ϕ=π/2 = η c (−d /2) + η c and η l |ϕ=π/2 = η c (d /2) + η c (4.19) and the jump of density equals η c d . Also note that E¨ (ρ 0 ), which enter the coeﬃcients c g and c l in (4.12), are even functions of ρ 0 , E¨ (ρ 0 ) = 2D −12Gρ 20 . Thus, E¨ (b) = E¨ (−b), E¨ (d ) = E¨ (−d ). In general, for the points ρ 02 and ρ 01 , which are symmetric relative to the rectilinear diameter ∆ = 0, there is E¨ (ρ 02 )/E¨ (ρ 01 ) = 1. In ﬁgure 5, there are presented the curves E = ln Ξρ 0 , E (ρ 0 ) and µ∗ (ρ 0 ) as functions of ρ max = ϕ+2πn b cos 3 for n = 0, 1, ρ max correspond to maximum values of the function E (ρ) in equation (4.1) for ln Ξρ 0 . The very thin line represents the minima of the function E (ρ), which is the region of thermodynamic instability. The quantity exp[N E (ρ)] characterises the probability measure for a state to exist. The points f f l and f f g , which are the ends of f l ÷ f f l and f g ÷ f f g determine the limits for stable density ﬂuctuations of liquid type arising in the gaseous phase and of gaseous type in a liquid. The points −b/2 and b/2 determine the limits of secondary maxima for gaseous and liquid phases, respectively. The points −d /2 and d /2 are the points of ﬁxed densities for gas in the gaseous phase and for liquid in the liquid stage in the two-phase system, which arise at the ﬁrst-order phase transition. In ﬁgure 5, the closeness of the curves µ∗ (ρ max ) and E (ρ max ) is obvious, which makes it possible, in what follows, to write E (∆) ≈ µ∗ (∆). Let us return to the initial expression for Ξρ 0 , given in (4.1) for Q < 0. Similarly to the case of Q > 0, 43001-14 Phase transition of the ﬁrst order in equations (2.8)–(2.12) we write ln Ξρ 0 in the form: ln Ξρ 0 = N E (∆), (4.20) keeping the notation of equations (2.8)–(2.12). However, now expression (4.12) is taken for ∆. The quantity E (∆), where ∆ is a solution to equation (4.2) speciﬁed in (4.12), is a curve that envelops the principal and the secondary maxima of the function E (ρ). It follows from ﬁgure 5 that the curves E (ρ max ) and µ∗ (ρ max ) stretch along each other. This is not surprising because E (ρ max ) consists of the sum µ∗ (ρ max ) and expressions Dρ 2max − Gρ 4max . In the critical region µ∗ ∼ τ5/2ν , and Dρ 2max −Gρ 4max ∼ τ3ν . Therefore, in what follows, we shall write E (∆) = µ∗ (∆) + D∆2 −G∆4 , µ∗ (∆) = −2D∆ + 4G∆3 , E (∆) ' µ∗ (∆) (4.21) and for the equation of state, taking into account (4.20) and (4.21): Θ ln Ξρ 0 , V N P ρ 0 = Θ µ∗ (∆), V P ρ0 = −b É ∆ É b. (4.22) The value of ∆ is given in (4.12). For µ∗ we may take two equivalent expressions. In the ﬁrst one, according to (1.27) µ∗ = a cos ϕ, 0 É ϕ É π, for the second one according to (1.20) µ∗ = −2Dρ max + 4Gρ 3max , where, based on (1.27), ρ max may take on two values: ρ max = ρ 02 = b cos ϕ + 2π , 3 0 É ϕ É π, and this means the branch of the gaseous type, and ρ max = ρ 01 = b cos ϕ , 3 0 É ϕ É π, describes the branch of the liquid type. 5. Properties of functions E (∆) and µ∗ (∆) in the region of Q < 0 In the region of Q < 0 for 0 É ϕ É π, there are two solutions to equations (1.20), (1.21): the solution of the gaseous type ρ 02 and that of the liquid type ρ 01 . Moreover, when ρ 02 (ϕ) for π/2 É ϕ É π describes the principal maximum, in the same interval of values ϕ the solution ρ 01 (ϕ) describes the secondary maximum. Inverse situation is seen in the region where 0 É ϕ É π/2: the principal maximum is due to ρ 01 , and the secondary one is due to ρ 02 . In general, this is governed by the probabilities w g and w l given in (4.12). In order to take into account both maxima, in the expression for µ∗ (ρ max ) we should write: µ∗ = µ∗ (∆(ϕ)) = −2D∆ + 4G∆3 , (5.1) where ∆ = ∆g + ∆l , ∆g = w g ρ 02 , ∆l = w l ρ 01 . 43001-15 I.R. Yukhnovskii 1 ,5 µ, 10 −5 1 ,0 0 ,5 ϕ 0 ,0 π -0 ,5 µ = a c o s (ϕ) 3 7 3 8 −2D ∆+4G ∆ f o r N = 1 0 -1 ,0 −2D ∆+4G ∆ f o r N = 1 0 -1 ,5 Figure 6. The comparison of the curves µ∗ = a cos ϕ and µ∗ = −2D∆(ϕ) + 4G∆(ϕ)3 as functions of the angle ϕ, 0 É ϕ É π; N = 107 and N = 108 . As N grows, the attitudes of maximum and minimum remain unchanged, the widths go to zero. See also ﬁgures 3 and 4 for w g and w l as well as for ∆g , ∆l and ∆ = ∆g + ∆l . Figure 6 presents the curves for the generalized chemical potential µ∗ (ϕ) as functions of ∆(ϕ): one curve is for µ∗ = a cos ϕ and the other is for µ∗ = −2D∆(ϕ) + 4G∆(ϕ)3 . We are going to work with the latter one. They almost fully coincide except for the narrow interval of values ϕ near ϕ = π/2. It follows from ﬁgures 1 and 6, in the region of Q < 0, that in the process of isothermal compression of the gas along ρ 02 (ϕ) for π/2 É ϕ É π, the system via ∆l enters the region of liquid-like values of µ∗ (∆g + ∆l ). In other words, at compression of the gaseous phase in the region Q < 0, the density ﬂuctuations of the liquid type emerge. Ultimately, at ϕ = π/2, a droplet of liquid appears in the gaseous phase. What to do next depends on the type of the shell the droplet possesses when the two-phase system arises. The simplest case is when the shell of the droplet is presented as a geometric sphere conﬁning the droplet. The shell does not contain particles. The question arises whether this is physically justiﬁed. In fact, the shell of the droplet has a volume, and is of the form of a spherical layer with a non-neglected thickness. All of them, i.e., two phases and the surface layer, should be present in the right-hand side of (5.1). The left-hand side of (5.1), which initially is of the form PV /Θ = ln Ξρ 0 , is written as a sum of terms dependent on the right-hand side of (5.1), in the expression ln Ξρ 0 . Returning to (5.1), for ∆(ϕ), in accordance with (4.12), we write ∆ = ∆g + ∆l , ∆g = w g ρ 02 , ∆l = w l ρ 01 . Then, E (∆) = E (∆g + ∆l ) = µ∗ (∆g + ∆l ) + D(∆g + ∆l )2 −G(∆g + ∆l )4 . (5.2) Let us write this as a sum E (∆g + ∆l ) = E g (∆g , ∆l ) + E l (∆g , ∆l ), µ∗ (∆g + ∆l ) = µ∗g (∆g , ∆l ) + µ∗l (∆g , ∆l ). (5.3) By doing so, we make no changes compared to either (5.1) or (5.2). The curves for ∆g and ∆l are presented in ﬁgure 3. Take the point ∆g = −d /2 on the curve ∆g , and the point ∆l = d /2 on the curve ∆l . Both correspond to equal probabilities w g = w l = 1/2 for w g and w l in ﬁgure 3. Obviously, ∆g (−d /2) + ∆l (d /2) = ∆(π/2) = 0. (5.4) In ﬁgures 1, 5, 6 the points ∆(π/2) = 0 and µ∗ (∆(ϕ))∆=0 = 0 correspond to the transition from the gaseous part µ∗ (∆(ϕ)) to the liquid part µ∗ (∆(ϕ)). At this point, ∆ = 0, 43001-16 µ∗ (∆) = 0, E (∆) = 0. (5.5) Phase transition of the ﬁrst order Our aim is to show that the liquid-gas phase transition occurs along the line ∆g = −d /2÷∆l = d /2, and that condition (5.5) is obeyed starting at ∆g = −d /2 and ﬁnishing at ∆l = d /2, and that this process can be described within the Gibbs statistics. Since the ﬁrst-order phase transition at a constant temperature τ occurs due to the external latent work of the pressure, we should demonstrate that at ϕ = π/2, and w l = w g = 1/2, the following is true ¯ ¯ µ∗g (ϕ)|ϕ=π/2 = µ∗g (∆g , ∆l )¯ ∆g =−d /2 = 0; ∆l =d /2 ¯ ¯ µ∗l (ϕ)|ϕ=π/2 = µ∗l (∆g , ∆l )¯ ∆g =−d /2 = 0; ∆l =d /2 ¯ ¯ E g (ϕ)|ϕ=π/2 = E g (∆g , ∆l )¯ ∆g =−d /2 = 0; ∆l =d /2 ¯ ¯ E l (ϕ)|ϕ=π/2 = E l (∆g , ∆l )¯ ∆g =−d /2 = 0. (5.6) ∆l =d /2 For this reason, we return to (5.2) and (5.3). Let us show that equations (5.5), (5.6) do take place along the line µ∗ (∆g , ∆l ) = 0 for ∆ in the interval ∆g = −d /2 ÷ ∆l = d /2. We write the expressions µ∗ (∆g + ∆l ) n and E (∆g + ∆l ), selecting the mixed products ∆m g ∆l and separating them for µg and µl and for E g and E l µ∗g (∆g , ∆l ) = −2D∆g + 4G∆3g + 12G∆g ∆2l , µ∗l (∆g , ∆l ) = −2D∆l + 4G∆3l + 12G∆l ∆2g . (5.7) Obviously, the sum of these expressions is equal to the initial form of µ∗ (∆g + ∆l ) in (5.2). Next, E g (∆g , ∆l ) = µ∗g (∆g , ∆l ) + D∆2g −G∆4g + D∆g ∆l h ³ ´ i −G 2 ∆3g ∆l + ∆g ∆3l + 3∆2g ∆2l , E l (∆g , ∆l ) = µ∗l (∆g , ∆l ) + D∆2l −G∆4l + D∆g ∆l h ³ ´ i −G 2 ∆3g ∆l + ∆g ∆3l + 3∆2g ∆2l . (5.8) The sum of these expressions equals E g (∆g , ∆l ) in (5.2). By substitution ∆g = w g ρ 02 and ∆l = w l ρ 01 we make sure that on both ends of the line of phase transition ∆g = −d /2 ÷ ∆l = d /2, ϕ = π/2, (5.9) (see ﬁgures 3 and 4) equalities (5.6) hold: µ∗g (−d /2, d /2) = 0, µ∗l (−d /2, d /2) = 0, E g (−d /2, d /2) = 0, E l (−d /2, d /2) = 0. (5.10) We conclude that the line −d /2 ÷ d /2 has the properties of the line of the ﬁrst-order phase transition in the gas-liquid system, at arbitrarily large N . ∗ The expressions for generalized chemical potentials µ∗ g and µl as well as for E g and E l consist of both “pure” phase terms such as −D∆g + 4G∆3g , and “mixed” products such as ∆g ∆2l , ∆2g ∆l , ∆g ∆l etc. On the line −d /2 ÷ d /2, the “mixed” products are of the same order of magnitude as the “pure” ones. Note that ∗ 5ν/2 in (5.8), quantities µ∗ , and the expression D∆2 −G∆4 is proportional to τ3ν . g , µl are proportional to τ In ﬁgure 7, there are presented the curves for the mixed products ∆g ∆2l , ∆2g ∆l . Their non-zero values are located near the limiting value ϕ = π/2 and characterise the width of transition layer between gas and n liquid. The size of the region of values ∆m g ∆l depends on N . At N increasing, it narrows and the region ∗ of the surface shell joins the surface µ = 0, for the angle ϕ = π/2. Hence, the magnitude of the maxima keep their values unchanged. As a result, the effect of the mixed products on the requirements (5.5), (5.6) remains unchanged. 43001-17 I.R. Yukhnovskii -5 4 2 ∆l ∆g + ∆l ∆g , 1 0 2 2 ϕ 0 0 ,0 0 ,5 1 ,0 -2 fo r N = 1 0 fo r N = 1 0 1 ,5 2 ,0 2 ,5 3 ,0 7 8 -4 Figure 7. The region of the surface layer described by the products ∆g ∆2l , ∆2g ∆l in the expressions for n E (∆g , ∆l ) and µ∗ (∆g , ∆l ) in (5.7), (5.8). The width of this region of values ∆m g ∆l depends on N . At N ∗ increasing, this region narrows and joins the surface ϕ = π/2, µ = 0. In addition, the magnitudes of maxima remain constant, and the effect of the mixed products on the requirement (5.5), (5.6) does not change. It is seen in equation (5.7) that the mixed product 12G∆g ∆2l forces the gaseous branch µ∗ g (∆g , ∆l ) into the point µ∗ = 0, ∆g = −d /2, and the product 12G∆g ∆2l nulliﬁes the liquid branch in the point µ∗ = 0, ∆ = d /2, and we get the curve presented in ﬁgure 8. It remains to clarify the problem of the horizontal branch in the equation of state, the fact of originating of the liquid phase in the gaseous one at isothermal compression of the gas, as well as the form of the left-hand side in equation (4.22) for the equation of state. µ, 1 0 -5 n This way, the existence of mixed products ∆m g ∆l in equations (5.7), (5.8) leads to the appearance ∗ ∗ ∗ of a shelf on the curve µ (∆) = µg (∆g , ∆l ) + µg (∆g , ∆l ) in the equation of state. Conditions (5.6) have a fundamental signiﬁcance because condensation into the liquid starts at density of gas η g = η c (1 − d /2), and terminates at density of liquid η l = η c (1 + d /2). In other words, the gas isotherm ends at the point 1 ,5 1 ,0 0 ,5 0 ,0 -0 ,0 7 5 -0 ,0 5 0 -0 ,0 2 5 0 ,0 0 0 0 ,0 2 5 0 ,0 5 0 0 ,0 7 5 -0 ,5 -1 ,0 ∆ ( ϕ) w l ρ0 1 ( ϕ) w g ρ0 2 ( ϕ) -1 ,5 Figure 8. The curve of µ∗ (∆) as a relult of combination of µ∗ = a cos ϕ and ∆ = w g ρ 02 (ϕ)+ w l ρ 01 (ϕ). The inﬂuence of the mixed products 12G∆g ∆2l on nullifying the branch µ∗ g (∆g , ∆l ) in the point ∆g = −d /2 and of 12G∆g ∆2l on nullifying the branch µ∗ (∆g , ∆l ) in the point ∆l = d /2 is shown. l 43001-18 Phase transition of the ﬁrst order ρ 02 = d , ∆g = −d /2 at density η g = η c (1 − d /2), and the liquid-phase isotherm ends at the point ρ 01 = d , ∆l = d /2 at density η l = η c (1 + d /2). On the line η c (1 − d /2) ÷ η c (1 + d /2), the phase transition occurs ∗ from the gas into liquid. This phase transition from gas into liquid along µ∗ g (−d /2; d /2) ÷ µl (−d /2; d /2) occurs due to the external latent work of the pressure. Therefore, according to (5.7), there are contributions into the energy density of the surface layer due ∗ to µ∗ g (∆g , ∆l ) and µl (∆g , ∆l ) h i µ∗g (∆g , ∆l ) = 12G∆g ∆2l , surf ¤ £ ∗ µl (∆g , ∆l ) surf = 12G∆g ∆2l , (5.11) 3 µ∗l (−d /2, d /2) − µ∗g (−d /2, d /2) = − Dd . 2 (5.12) Since ∆g = w g ρ 02 , ∆l = w l ρ 01 , and the probabilities w g and w l depend on the number of particles N and the angle ϕ, as is seen in ﬁgure 8. At N increasing, the products of w g and w l degenerate into δ−functions placed along ϕ = π/2. This means that in the limit of large N , the two-phase liquid-gas system degenerates into the system gas — a droplet of liquid divided by the surface without particles. The energy of the phase layer transforms into the energy associated with the creation of a geometric surface of division. The two-phase system emerges: gas and a droplet of liquid of radius R dr . If for the gaseous phase P g is the pressure, Vg is the volume, P g Vg¢is the energy, then for the droplet of liquid the ¡ 3 R , the pressure inside it P = P + 2α/R volume is 4π l g dr , where α is the surface tension coeﬃcient, the 3 ¡dr ¢¡ ¢ 3 2 energy P l Vl = P g + 2α/R dr 4πR dr /3 and plus the surface tension energy −αS dr = −α4πR dr . Hence, the formation of the droplet with a geometric interface is concerned with an additional energy of the surface 3 3 α 4πR dr 2α 4πR dr − α4πR 2 = − , R dr 3 R dr 3 (5.13) where the surface-energy density is −α/R dr . The magnitude of the surface energy is connected with quantity (5.11). Also, the change of the volume for the gaseous phase takes place, from Vl to Vg = V − (4πR 3 /3)(n l /n g ). 6. The equality of the chemical potentials of gas and liquid on the condensation-boiling line µ∗ (∆g , ∆l ) For µ∗ , by deﬁnition, we have equation (1.9): µ∗ = β(µ − µ0 ) + ξ − |α(0)|(1 − ∆). (6.1) In the approximation accepted in this research, the chemical potential of the reference system is βµ0 = ξ. (6.2) Then, for the chemical potential µ we get the following expression: µ = −Θ |α(0)|(1 − ∆), ∆ = w g ρ 02 + w l ρ 01 . (6.3) On the line of phase transition the conditions (5.5), (5.6), µ∗ (∆g , ∆l ) = 0, hold together with Ng + Nl = N , dNg = −dNl , (6.4) where Ng and Nl are the number of particles in the gas phase and in the liquid phase, respectively, µg , µl are the chemical potentials of the liquid and of the gas, respectively. According to (5.6) βµg = |α(0)|g (1 − ∆g ), βµl = |α(0)|l (1 − ∆l ). (6.5) 43001-19 I.R. Yukhnovskii According to deﬁnition (1.12) N 1 ˜ |Φ(0)|, V Θ 6 1 ˜ |α(0)|g = η g |Φ(0)|, πσ3 Θ 6 1 ˜ |α(0)|l = η l |Φ(0)|, πσ3 Θ |α(0)| = then η g — gas density, η l — liquid density in a droplet. (6.6) We consider the beginning of condensation (boiling) ∆ = ∆g + ∆l , 1 ∆g = ρ 02 , 2 1 ∆l = ρ 01 , 2 ρ 02 = −d , ρ 01 = d , 1 ∆g = − d , 2 1 ∆l = d . 2 (6.7) The densities η g and η l are expressed via the value η c in the critical point and via ∆ based on the relationship (2.15): η = ηc ∆ + ηc , ³ 1 ´ ηg = ηc 1 − d , 2 ³ 1 ´ ηl = ηc 1 + d . 2 (6.8) Then, for (5.4) we obtain ³ 6 d2 ´ η 1 − , c πσ3 4 ³ 6 d2 ´ ˜ η . µl = |Φ(0)| 1 − c πσ3 4 ˜ µg = |Φ(0)| (6.9) Thus, µg = µl . (6.10) It was proved that on the line µ = 0 there exists a two-phase system of a gas phase and of a liquid phase: the maternal sphere of the gas and the droplet of liquid inside it and the chemical potentials of both phases are the same. Let us imagine that conditions (6.7) and (6.8) necessary for equalities (6.9) and (6.10), may be also valid for other quantities ∆g and ∆l , not only for that of (6.7), which are located symmetrically to the point of the intersection of the rectlinear diameter with the line µ∗ = 0 that is ∆g = −∆l . Such points are located at the ends of the intervals which pass through the origin µ∗ = 0 and touch the lines ρ 02 and ρ 01 that are the curves of the maxima for gaseous and liquid phases, respectively, see ﬁgure 1. For each pair of points, denoted by ρ˜ 02 and ρ˜ 01 , the equalities E (ρ˜ 02 ) = E (ρ˜ 01 ) take place, since ρ˜ 02 = −ρ˜ 01 , and E (ρ 0 ) is an even function. Thus, for each pair of end-point values there is w(ρ˜ 02 ) = w(ρ˜ 01 ). Hence, for ˜ = w(ρ˜ 02 )ρ˜ 02 + w(ρ˜ 01 )ρ˜ 01 = 0, µ∗ = 0, each pair, the following is true: w(ρ˜ 02 )ρ˜ 02 = −w(ρ˜ 01 )ρ˜ 01 and ∆ ˜ ˜ E (0) = 0, and also µg (∆g ) = µl (∆l ). However, conditions (5.6) are not obeyed if ∆g = −d /2, ∆l = d /2 are ˜ g = −∆ ˜ l. substituted by any of ∆ Therefore, neither approaching nor moving away from the densities of gas and liquid compared to (6.7) occur. The conclusion is that densities for gaseous and liquid phases do not change at the ﬁrstorder phase transtition. The density of gas remains ﬁxed ∗ ³ ´ η g = η c ∆g + η c = η c 1 − d /2 . (6.11) ³ ´ η l = η c ∆l + η c = η c 1 + d /2 . (6.12) The density of liquid remains ﬁxed The chemical potentials of both phases are the same µg = µl . This remains true as long as the phase transition takes place. As a result, under the action of latent work of the pressure, all the gas condenses into liquid. A jump-like change of the gas density takes place ηl − ηg = ηc d . 43001-20 (6.13) Phase transition of the ﬁrst order 7. The process of condensation of gas into liquid (the boiling of liquid) This section begins with (4.22) PV = N E (∆), Θ ˜ (∆), P = ΘηE ˜ =Θ Θ 6 . πσ3 (7.1) Let us follow the process of isothermal quasistatic compression in the gaseous phase (see section 1). The process takes place along the line ∆ = w g ρ 02 + w l ρ 01 . Starting at ∆ = ∆g = w g ρ˜02 : ρ 02 = −b, w g = 1, ∆ = −b, w l = 0, (7.2) where (see ﬁgures 1–4) there is a homogeneous gaseous system (see ﬁgures 1 and 5). Next, after the quantities ρ 02 and ∆ increase, we get ρ 02 > −b , ∆ > −b , and the probabilities are as follows: w g < 1, w l > 0, w g + w l = 1. In a gaseous system, there emerge density ﬂuctuations of a liquid type. During the quasistatic process for −b É ρ 02 É f f l ﬂuctuations are unstable, because, according to ﬁgure 5, under the line f f l − f l E (∆) < 0 and exp[N E (∆)] → 0. Being further compressed, while crossing the line f f l − f l (see ﬁgure 5) in a gaseous system, there emerge stable density ﬂuctuations of a liquid type because E (∆) > 0 therein. The argument ∆ is determined by equation (4.12). Here, both terms in ∆ become of non-zero value. However, the term w g ρ 02 in ∆ is a principal one. As earlier, we have only gas, and there is a single-phase system (see ﬁgure 4). When the external pressure increases, the pressure in the gaseous system also increases. The density of a gaseous system grows, in accordance with (2.15). The gaseous system evolves along the same isotherm µ∗ g (ρ 02 ). Now we get to the limiting densities η g = η c (1 + w g ρ 02 ), η l = η c (1 + w l ρ 01 ), where ρ 02 = −d , ρ 01 = d , E (d ) = E (−d ), w g = w l = 1/2. (7.3) Based on (4.12), (5.6), (5.8)–(5.10), it follows ∆ = ∆g + ∆l = 0, µ∗ = µ∗g + µ∗l = 0, E = E g + E l = 0. (7.4) This state describes the end of the curve for the principal maxima for a gaseous phase and the beginning of the principal maxima for a liquid phase. Here, we arrive at (5.10) and (7.1). When work is done on the gaseous system, within the gas system there arise processes that evolve along the curve ∆l = w l ρ 01 , although under the conditions of a secondary maximum, as is seen in ﬁgures 1 and 5. The system does an internal work to create stable density ﬂuctuations of liquid type, which in the case of ρ 02 = −d , ∆g = −d /2, ρ 01 = d , ∆l = d /2 causes the creation of a droplet. In order to determine the properties of a droplet, it should have exact boundaries. Equations (5.7) ∗ and (5.8) for µ∗ , µ∗ g , µl , E , E g , E l are “bulk” formulas. In these formulas, the interface is a phase layer with a number of particles. Our further calculations are performed for the case of very large N . The model of a two-phase region is considered as gas and a droplet of liquid arising in gas, the interface n ∗ is a geometric surface without particles. The mixed products ∆m g ∆l , entering µ (∆g , ∆l ) and E (∆g , ∆l ), presented in (5.7) and (5.8) are replaced with the energy of additional pressure in the droplet of liquid 3 2 (4πR dr /3)(2α/R dr ) and the surface energy −α4πR dr . The two-phase system exists only along the line µ∗ = 0. On this line, the chemical potentials of gaseous and liquid phases are the same. Equality (5.11) is invariant of the process of condensation (i.e., boiling). This is true for the case when the thermodynamic potential Φρ 0 from the beginning to the end of the process of condensation (i.e., boiling) remains constant Φρ 0 = Ng µg + Nl µl = N µg = N µl , (7.5) since Ng +Nl = N remains constant and is equal to the initial number of particles, and dNg = −dNl , which, in fact, holds. The amount of gas decreases, while the amount of liquid increases till all the gas transforms into liquid. Liquid phase originates in the form of a droplet conﬁned by the interface S of some surface energy. The sharp bend of the isotherm in ﬁgure 8 is related to the instant of the creation of a droplet, and the sum of the moduli of the mixed products (5.11) is equal to the energy of the droplet. 43001-21 I.R. Yukhnovskii 7.1. The energy of the surface tension According to (4.12), (6.5)–(6.7), along the line µ∗ = 0 we have, ∆ = ∆g + ∆l , ∆g = w g (−d ) = −d /2, ∆l = w l d = d /2, w g = w l = 1/2. For the basic relationship (7.1), the left-hand side takes the form P g Vg + P l Vl − αS l , (7.6) where P g , P l are the pressures in the gas and in the liquid phases, Vg and Vl are the volumes of the gas and of the liquid phases, S l is the surface of the droplet of liquid, α is the surface-tension coeﬃcient. In the right-hand side of (7.1) we obtain: Ng E g (∆g ) + Nl E l (∆l ) ' Ng µ∗g (∆g ) + Nl µ∗l (∆l ). (7.7) ∗ According to the conclusions drawn in (6.5) and (6.6), the quantities µ∗ g , µl take on ﬁxed values. In accordance with (6.5)–(6.7): E g (∆g ) = µ∗ (∆g ) = −2D∆g + 4G∆3g = 3Dd /4, E l (∆l ) = µ∗ (∆l ) = −2D∆l + 4G∆3l = −3Dd /4, (7.8) and the process of the phase transition is reduced to redistribution of Ng and Nl . The width of the interface between gas and liqud droplet is considered to be zero. Liquid phase is a droplet of radius R that originates in the gaseous phase. The droplet has strict boundaries. Combining (7.5)–(7.8), for (7.1) on the line µ∗ = 0 we obtain: P g Vg + P l Vl − αS l = Θ{Ng E g + Nl E l }. (7.9) The emergence of a liquid droplet, restricted by a shell, in a gaseous phase is the essence of the ﬁrst order phase transition. The droplet is of a spherical form. The surface tension energy is minimal. For (7.9), we get: 4π 3 η l R , 3 dr η g µ ¶ 2α 4π 3 R , P l Vl = P g + R dr 3 dr P g Vg = P g V − P g (7.10) where V is the initial volume, R dr is the radius of the ﬁrst droplet, η l , η g are the densities of gas and liquid, respectively. For the energy of the droplet surface, we have: 2 − αS = −α4πR dr =− 3 3α 4πR dr . R dr 3 (7.11) In the right-hand side of (7.9), it is written Θ[Ng E g + Nl E l ] = Θ[N E g + Nl (E l − E g )]. (7.12) Collecting (7.8)–(7.10) in equation (7.7), we obtain: P g Vg + P l Vl − αS = P g V − ¶ 3 3 µ 4πR dr α 4πR dr ηl − Pg −1 , R dr 3 3 ηg ¶# 3 µ 3 4πR dr α 4πR dr ηl V− −1 − = N ΘE g + Nl (E l − E g ). 3 ηg R dr 3 " Pg (7.13) 4π 3 3 The change of gas volume by 4π 3 R dr η l /η g − 1 , that takes place in the product P g V − 3 R 0 η l /η g − 1 in (7.13) corresponds to those changes in the quantity E g (∆g ) which were observed at the transition from E g (∆g ) at ∆g = −b to E g (∆g ) at ∆g = w g (−d ) = −d /2. We compare ¡ ¢ µ · ¶¸ 4πR 3 η l Pg V − − 1 = N ΘE g . 3 ηg 43001-22 £ ¡ ¢¤ (7.14) Phase transition of the ﬁrst order Hence, for the speciﬁc7 surface tension energy of the droplet of radius R , we have − α = n l (E l − E g )Θ. R (7.15) 3 Let us replace E l and E g with (7.8), and n l with n l = Nl /(4πR dr /3) = η l 6/πσ3 . In the right-hand side of (7.15), we obtain: Θn l (E l − E g ) = −η l 9Θ Dd , πσ3 (7.16) which conforms to (5.10). Therefore, based on (7.15), there is a contribution to speciﬁc pressure P total from the surface tension energy of the liquid droplet of radius R : P total = − 3 6 3 α ˜ l Dd , = − Θ 3 η l Dd = − Θη R 2 πσ 2 (7.17) ˜ = Θ 6/πσ3 . where Θ The quantity P total is proportional to τ5/2ν . It governs the rate of decreasing the surface tension energy at approximatiion to the critical point τ = 0. 7.2. Latent work of condensation (latent heat of boiling) Condensation begins at the moment of the creation of the droplet with the surface energy (7.17). As was mentioned above, condensation of a gas starts from the state N E (∆) = N µ∗g (∆g ), 1 ∆g = − d , 2 ρ 02 = −d , 1 ηg = ηc − d ηc 2 and terminates in a liquid system in the state N E (∆) = N µ∗ (∆l ), ρ 01 = d , ∆l = 21 d and η l = η c + 12 d η c . At l the moment when condensation starts, the pressure is equal to P g = n g ΘE g (∆g ) = n g Θµ∗ (∆g ). (7.18) Replacing ∆g with ∆g = w g (−d ) = − 12 d and 1 3 µ∗ (∆g ) = d D − d 3G = Dd . 2 4 (7.19) 3 P g = n g Θ Dd . 4 (7.20) Then, The initial volume of gas is a sphere of radius R 0 , of volume V0 = 4πR 03 /3. The ﬁnal volume of gas after transferring into a liquid during condensation at a ﬁxed pressure is the volume of liquid Vl in the form of a sphere of radius R l : Vl = 4π 3 R . 3 l All N particles of a gaseous sphere V0 passed into a sphere of liquid Vl . Thus, 1 1 V0 = Vl , vg vl Vl = vl V0 , vg the total change of the volume is µ (Vl − V0 ) = ¶ vl − 1 V0 . vg (7.21) The work performed is the latent work (i.e., heat) of the process of condensation: µ −P g (Vl − V0 ) = −P g ¶ vl 4π 3 −1 R . vg 3 l 7 By a volume unity. 43001-23 I.R. Yukhnovskii Having substituted (7.20), the speciﬁc latent work of condensation is8 µ ¶ vl 3 3 L = − 1− n g Θ Dd = −(v g − v l ) ΘDd n g2 . vg 4 4 (7.22) Condensation at isothermal compressibility is a dynamical process. Equation (7.9) concerns the states of a system in which two phases have already been created: the initial gaseous phase and a new liquid phase in the form of a droplet. The conditions ∆ = 0, µ∗ = 0, E (∆) = 0 require an external pressure to be applied in order to do the latent work of condensation. As a result, the gaseous phase disappears (or remains as a nascent bubble of gas), and we write: P l Vl = L(Vg − Vl ). (7.23) [For a nascent bubble of gas, if any, P g Vg nascent −αS gnascent is energy, and we can write: P l Vl +P g Vg nascent − αS gnascent = L(Vg − Vg nascent − Vl )]. For the two-phase system along the line µ∗ = 0, E (∆) = 0 for external pressure we write down: P µ∗ =0 = L ³ vg vl ´ −1 . (7.24) Using L , we get µ ¶ vl 3 P µ∗ =0 = 1 − n g ΘDd . vg 4 (7.25) On the curve of the equation of state, this expression represents a horizontal fragment of the transition from a gaseous phase into a liquid phase under the isothermal compression of a gaseous phase. 8. The equation of state Now we shall summarize the results of our calculations in the critical region of τ < 0, |τ| < τ∗ and |η − η c | < |η∗ − η c |, where τ∗ and |η∗ − η c | are the boundaries of the critical region presented in expressions (1.15), (1.16) and (2.13). The condensation of gas at isothermal compression is a process analogous to the boiling at an isothermal decrease of pressure of a liquid. This is a quasiequilibrial directed process. Equation (7.9) concerns the states of a system in which two phases have already been formed: the initial gaseous phase and a new liquid one, in the form of a droplet. (The droplet has got a shell that separates liquid from gas, and its creation is the essence of the ﬁrst-order phase transition.) Conditions ∆ = 0, µ∗ = 0, E (∆) = 0 require an external pressure to be applied in order to do the latent work of condensation (7.22). As a result, in (7.9) the gaseous phase disappears (or remains as a nascent bubble of gas). The conditions existing on the line of condensation µ∗ = 0, E (∆) = 0, ∆g = −d /2 ÷ ∆l = d /2, η g = η c (1 − d /2) ÷ η c (1 + d /2) = η l are determined by (6.9), i.e., by the equality of the chemical potentials of different phases. The thermodynamic potential Φρ 0 of a system on the line of condensation remains constant Φρ 0 = Ng µg + Nl µl = N µg = N µl , (8.1) according to (6.9) ˜ Φρ 0 = N |Φ(0)| 6 η c (1 − d 2 /4). πσ3 (8.2) On the line of the ﬁrst-order phase transition, equations (7.9) hold. To obtain these equations, the mixed n products ∆m g ∆l were nulliﬁed in the exact relationships (5.7)–(5.8), in particular those in (5.11). Instead of them, the surface tension energy of a droplet was introduced as well as an additional internal pressure of a droplet. This allowed us to describe the magnitude of energy associated with the creation of the liquid droplet. The droplet has got a strictly deﬁned interface. We have described the quantity of the latent work of the pressure. 8 Dimension of 43001-24 L coincides with the dimension of pressure. Phase transition of the ﬁrst order N=10 7 τ = - 10 -5 1,5x10 µ* -3 -5 1,0x10 -6 5,0x10 ∆G -d/2 -0,08 -0,06 -0,04 0,0 -0,02 0,00 d/2 0,02 0,04 -6 -5,0x10 0,06 0,08 ∆L -5 -1,0x10 -5 -1,5x10 ∗ Figure 9. The graph of E (∆g ) ' µ∗ g (∆g ) and E (∆l ) ' µl (∆l ); circles denote densities of the gaseous phase ∗ = η (1 + d /2), which remains constant during the liquid-gas η∗ = η (1 − d /2) and of the liquid phase η c c g l ϕ+4π phase transition. Both solutions are in the region of instability ρ 03 = b cos 3 (see ﬁgures 1 and 5). The jump of the density at phase transition is η l − η g = η c d . The slope of the gaseous branch ∆g while approaching the point ∆g = −d /2 is equal to 0. The slope of the liquid branch while turning away from the point ∆l = d /2 is equal to 0. The task of our further research will be to clarify the problem of the thermodynamic limit and to compare the theory with experimental results. We will need a total expression for ln Ξ, not only ln Ξρ 0 . These urgent problems are easier to consider farther from the critical point, since there are no requirements following from the renormalization-group symmetry near T = Tc . Thus, the general expression for the equation of state takes the form P ρ0 = n ˜ η g µ∗g (∆g , ∆l )∆ É−d /2, µ∗ <−h/2 Θ g + [“shelf”, the region of constant pressure]−h/2Éµ∗ Éh/2, −d /2É∆Éd /2 ª + η l µ∗l (∆g , ∆l )∆Êd , µ∗ >h/2 , (8.3) where µ∗l (∆g , ∆l ) = −2D∆l + 4G∆3l + 12G∆l ∆2g , h → 0, µ∗g (∆g , ∆l ) = −2D∆g + 4G∆3g + 12G∆g ∆2l , h → 0. On the¡ “shelf”, ¢the external pressure is equal to the speciﬁc work of condensation according to (7.9), ˜ 3 Dd , and the thermodynamic potential remains constant, Φρ = N µg (−d /2, d /2) = P µ∗ = 1 − v l /v g η g Θ 0 4 N µl (−d /2, d /2). The phase transition is a dynamical process. It takes place without any breaks as a quasistatic isothermal compression of a gaseous phase. There are ﬁve stages of the process. The ﬁrst one is a gaseous phase p with density η = η c (−b) + η c . Here, w g (ϕ) = w g (π) = 1, w l = 0, ∆ = ∆g = −b , b = 2D/(3G), ∆l = 0. When moving along ∆ = ∆g = w g (ϕ)ρ 02 (ϕ), there are density ﬂuctuations of a liquid type in the gaseous state along the secondary maximum ∆l = w l ρ 01 . The region of a sharp bend of the isotherm µ∗ g (∆g , ∆l ) is reached due to the terms in (5.11) (see ﬁgure 9). In a gaseous phase, a liquid droplet appears with an additional, in comparison with the gaseous system, energy density (7.15)–(7.18). A liquid droplet that emerges means the end of the ﬁrst stage and the beginning of the second stage, being the main stage of the ﬁrst-order phase transition. The thickness of interface, accordint to ﬁgure 7 is equal to half-width of the maximum for the exn pressions 12∆g ∆2l , 12∆l ∆2g , and 2∆g ∆l . In general, the interface includes all the mixed products 2∆m g ∆l , ∗ via µ (∆g , ∆l ) and E (∆g , ∆l ). All terms of the mixed products distribute accordingly in equations (5.7) 43001-25 I.R. Yukhnovskii and (5.8) between E g and E l in such a way that in the ﬁxed limiting points of a phase transition ∆g = −d /2 and ∆l = d /2, the fundamental equalities (5.6) hold. The third stage is considered to be the result of the action of the external latent work of pressure on the two-phase system, presented in (7.23)–(7.24). Pressure and temperature do not change, accoridng to (6.5)–(6.6). Gas transfers into liquid according to (7.5). The external latent work of pressure has been done. The density jump ∆η = η l − η g = η c d takes place. Everything ends with the fourth stage, which is a liquid with a nascent bubble of gas. Finally, the ﬁfth stage is a liquid of density growing from the value η l = η c d /2 + η c to η l = η c b + η c . Thus, the events in the region Q < 0 ﬁnish. The region of liquid densities begins, where the discriminant Q > 0. This is a homogeneous liquid phase. We started with an initial sphere of gas in a boundless system. By isothermal compression, we reached a point of the ﬁrst-order phase transition. Gas transferred into liquid, condensed into a droplet. The quantities that describe the process were determined. We worked within the framework of the grand canonical distribution for the Gibbs statistics. The successful solution of the problem is ensured by the choice of a related phase space of collective variables, among which there is a variable ρ 0 . This variable is responsible for a description of the phase transition and the density jump. Calculations are carried out using the quartic basic measure density. The zero-approximation performed in this work provides good results on both the critical indexes and the description of the mechanism of density jump at the condensation of gas into liquid. We have gone through all the events of the ﬁrst-order phase transition below the critical point. The curve for the equation of state has a form presented in ﬁgure 9. A peculiar thing is the density jump p between the points η g = η c (1 − d /2) and η l = η c (1 + d /2), where d = D/2G , where D and G are the coeﬃcients at the second and fourth power of ρ 0 , respectively, in the initial Hamiltonian E (ρ 0 ). Figure 1 resembles the Van der Waals isotherm. The theory put forward by a great physicist [17], was the basis for numerous researches devoted to the gas-liquid phase transition, including the present research. At the same time, this work is not void of numerous evident drawbacks. Among them we note: • we worked in a very narrow regions of temperatures and densities, near the critical point; • when calculating the Jacobian of transition to the phase space of collective variables, the statistical weight was introduced, the distribution over short-range interaction between hard spheres. The cumulants of the Jacobian, M1 , M2 , M3 and M4 , were taken at k = 0. This was connected with the fact that they have wide “shelves” at k = 0. The cumulants were determined via compressibility of the reference system. We neglected the terms proportional to k 2 , which were proportional to the binary correlation function for the system of hard spheres; • we neglected the effect of integration over ρ k for k , 0 on the form of the isotherm of phase transition. • we restricted the expression for the chemical potential of the reference system to βµ0 = ξ. We postulated that the dependence E (ρ max ) = E (∆), being accurate in the region of Q > 0, remains also accurate for Q < 0, if the sum ∆ = ∆g + ∆l is assumed as ∆. All the above demands a more detailed consideration. The mutually-coordinated short-range and long-range interaction being taken into account will be one of the main tasks in further research. Acknowledgements The author expresses his gratitude to I.M. Mryglod and O.L. Ivankiv for help, to M.P. Kozlovskii and all the researchers of the Department of the Statistical Theory of Condensed Matter for useful discussions; to O.V. Patsahan for the data on the property of the parameter ∆ in the critical region; to A.D. Trokhimchuk for discussions concerning the short-range order in condensed systems; to Yu.V. Holovach for proofreading the paper and consistent remarks. Our special thanks go to R. Romanik for valuable considerations and the graphs, and to V.O. Kolomiets for collaboration. 43001-26 Phase transition of the ﬁrst order References 1. Yukhnovskii I., Kolomiets V., Idzyk I., Condens. Matter Phys., 2013, 16, 23604; doi:10.5488/CMP.16.23604. 2. Yukhnovskii I., Idzyk I., Kolomiets V., J. Stat. Phys., 1995, 80, 405; doi:10.1007/BF02178366. 3. Yukhnovskii I.R., Phase Transitions of the Second Order: Collective Variables Method, World scientiﬁc, Singapore, 1987. 4. Yukhnovs’kii I., Riv. Nuovo Cimento, 1989, 12, No. 1, 1; doi:10.1007/BF02740597. 5. Yukhnovskii I.R., Kozlovskii M.P., Pylyuk I.V., Microscopic Theory of Phase Transitions in Three-Dimentional Systems, Eurosvit, Lviv, 2001 (in Ukrainian). 6. Martynov G.A., Phys-Usp., 1999, 42, No. 6, 517; doi:10.1070/PU1999v042n06ABEH000543. 7. Balescu R., Equilibrium and Non-Equilibrium Statistical Mechanics, John Wiley & Sons, New York ,1975. 8. Landau L.D., Lifshitz E.M., Statistical Physics, A Course of Theoretical physics Vol. 5, 3-rd Edn., Nauka, Moscow, 1976 (in Russian). 9. Hansen J.P., McDonald I.R., Theory of Simple Liquids, 3-rd Edn., Academic Press, London, 2006. 10. Derzhko O., Mygal’ V., Some Problems of Non-Homogeneous Classical-Fluid Theories: Lecture Notes, Lviv University Press, Lviv, 1999 (in Ukrainian). 11. De Pablo J.J., Yan Q., Escobedo F.A., Annu. Rev. Phys. Chem., 1999, 50, No. 1, 377; doi:10.1146/annurev.physchem.50.1.377. 12. Anisimov M.A., Condens. Matter Phys., 2013, 16, 23603; doi:10.5488/CMP.16.23603. 13. Widom B., J. Chem. Phys., 1965, 43, No. 11, 3898; doi:10.1063/1.1696618. 14. Kadanoff L.P., Physics, 1966, 2, 263. 15. Wilson K.G., Phys. Rev. B, 1971, 4, No. 9, 3174; doi:10.1103/PhysRevB.4.3174. 16. Wilson K.G., Phys. Rev. B, 1971, 4, No. 9, 3184; doi:10.1103/PhysRevB.4.3184. 17. Van der Waals J.D., In: Nobel Lectures in Physics 1901–1921, World Scientiﬁc, Singapore, 1998, 254-265. Фазовий перехiд 1-го роду в областi критичної точки газ–рiдина I.Р. Юхновський Iнститут фiзики конденсованих систем НАН України, вул. I. Свєнцiцького, 1, 79011 Львiв, Україна Розглядається поведiнка системи взаємодiючих частинок в областi температур нижче критичної точки T É Tc . Завершується розрахунок великої статистичної суми, початий у попереднiх роботах у методi колективних змiнних. За базову густину мiри береться четвiрний (а не Гаусовий) розподiл. Описанi подiї пов’язанi з фазовим переходом 1-го роду, що вiдбуваються в результатi iзотермiчного квазiстатичного стиснення газу. Видiлена лiнiя µ∗ (η) = 0, на якiй у газовiй фазi пiд дiєю тиску виникає крапля рiдини. Має мiсце рiвнiсть хiмiчних потенцiалiв газової i рiдкої (у краплi) фаз; знайдено величину поверхневої енергiї краплi, розраховано скриту роботу конденсацiї, визначено скачок густини, написано рiвняння стану. Робота становить певну кiнцеву стадiю дослiджень в областi температур i густин, що включає у собi критичну точку рiдина–газ. Ключовi слова: критична точка, фазовий перехiд першого роду, велика статистична сума, колективнi змiннi 43001-27

1/--pages