Improved first order mean spherical approximation for simple fluids

A perturbation approach based on the first-order mean spherical approximation (FMSA) is proposed. It consists in adopting a hard-sphere plus short-range attractive Yukawa fluid as the novel reference system, over which the perturbative solution of the Ornstein-Zernike equation is performed. A choice...

Повний опис

Збережено в:

Бібліографічні деталі
Дата:	2011
Автори:	Hlushak, S., Trokhymchuk, A., Nezbeda, I.
Формат:	Стаття
Мова:	English
Опубліковано:	Інститут фізики конденсованих систем НАН України 2011
Назва видання:	Condensed Matter Physics
Онлайн доступ:	http://dspace.nbuv.gov.ua/handle/123456789/120011
Теги:	Додати тег Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:	Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:	Improved first order mean spherical approximation for simple fluids / S. Hlushak, A. Trokhymchuk, I. Nezbeda // Condensed Matter Physics. — 2011. — Т. 14, № 3. — С. 33004: 1-8. — Бібліогр.: 19 назв. — англ.

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine

id	irk-123456789-120011
record_format	dspace
spelling	irk-123456789-1200112017-06-11T03:04:53Z Improved first order mean spherical approximation for simple fluids Hlushak, S. Trokhymchuk, A. Nezbeda, I. A perturbation approach based on the first-order mean spherical approximation (FMSA) is proposed. It consists in adopting a hard-sphere plus short-range attractive Yukawa fluid as the novel reference system, over which the perturbative solution of the Ornstein-Zernike equation is performed. A choice of the optimal range of the reference attraction is discussed. The results are compared against conventional FMSA/HS theory and Monte-Carlo simulation data for compressibility factor and vapor-liquid phase diagrams of the medium-ranged Yukawa fluid. Proposed theory keeps the same level of simplicity and transparency, as the conventional FMSA/HS approach does, but shows to be more accurate. Пропонується новий п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льш точною. 2011 Article Improved first order mean spherical approximation for simple fluids / S. Hlushak, A. Trokhymchuk, I. Nezbeda // Condensed Matter Physics. — 2011. — Т. 14, № 3. — С. 33004: 1-8. — Бібліогр.: 19 назв. — англ. 1607-324X PACS: 01.65.Q DOI:10.5488/CMP.14.33004 arXiv:1202.4260 http://dspace.nbuv.gov.ua/handle/123456789/120011 en Condensed Matter Physics Інститут фізики конденсованих систем НАН України
institution	Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection	DSpace DC
language	English
description	A perturbation approach based on the first-order mean spherical approximation (FMSA) is proposed. It consists in adopting a hard-sphere plus short-range attractive Yukawa fluid as the novel reference system, over which the perturbative solution of the Ornstein-Zernike equation is performed. A choice of the optimal range of the reference attraction is discussed. The results are compared against conventional FMSA/HS theory and Monte-Carlo simulation data for compressibility factor and vapor-liquid phase diagrams of the medium-ranged Yukawa fluid. Proposed theory keeps the same level of simplicity and transparency, as the conventional FMSA/HS approach does, but shows to be more accurate.
format	Article
author	Hlushak, S. Trokhymchuk, A. Nezbeda, I.
spellingShingle	Hlushak, S. Trokhymchuk, A. Nezbeda, I. Improved first order mean spherical approximation for simple fluids Condensed Matter Physics
author_facet	Hlushak, S. Trokhymchuk, A. Nezbeda, I.
author_sort	Hlushak, S.
title	Improved first order mean spherical approximation for simple fluids
title_short	Improved first order mean spherical approximation for simple fluids
title_full	Improved first order mean spherical approximation for simple fluids
title_fullStr	Improved first order mean spherical approximation for simple fluids
title_full_unstemmed	Improved first order mean spherical approximation for simple fluids
title_sort	improved first order mean spherical approximation for simple fluids
publisher	Інститут фізики конденсованих систем НАН України
publishDate	2011
url	http://dspace.nbuv.gov.ua/handle/123456789/120011
citation_txt	Improved first order mean spherical approximation for simple fluids / S. Hlushak, A. Trokhymchuk, I. Nezbeda // Condensed Matter Physics. — 2011. — Т. 14, № 3. — С. 33004: 1-8. — Бібліогр.: 19 назв. — англ.
series	Condensed Matter Physics
work_keys_str_mv	AT hlushaks improvedfirstordermeansphericalapproximationforsimplefluids AT trokhymchuka improvedfirstordermeansphericalapproximationforsimplefluids AT nezbedai improvedfirstordermeansphericalapproximationforsimplefluids
first_indexed	2025-07-08T17:04:23Z
last_indexed	2025-07-08T17:04:23Z
_version_	1837099138242576384
fulltext	Condensed Matter Physics, 2011, Vol. 14, No 3, 33004: 1–8 DOI: 10.5488/CMP.14.33004 http://www.icmp.lviv.ua/journal Improved first order mean-spherical approximation for simple fluids S. Hlushak1, A. Trokhymchuk1,2, I. Nezbeda3 1 Institute for Condensed Matter Physics of the National Academy of Sciences of Ukraine, 1 Svientsitskii Str., 79011 Lviv, Ukraine 2 Department of Chemistry and Biochemistry, Brigham Young University, Provo, UT 84602, USA 3 Faculty of Science, J.E. Purkinje University, 400 96 Usti n. Lab., Czech Republic Received July 1, 2011, in final form July 25, 2011 A perturbation approach based on the first-order mean-spherical approximation (FMSA) is proposed. It con- sists in adopting a hard-sphere plus short-range attractive Yukawa fluid as the novel reference system, over which the perturbative solution of the Ornstein-Zernike equation is performed. A choice of the optimal range of the reference attraction is discussed. The results are compared against conventional FMSA/HS theory and Monte-Carlo simulation data for compressibility factor and vapor-liquid phase diagrams of the medium-ranged Yukawa fluid. The proposed theory keeps the same level of simplicity and transparency as the conventional FMSA/HS approach does, but turns out to be more accurate. Key words: hard-core Yukawa fluid, mean-spherical approximation, perturbation theory, Monte Carlo simulations PACS: 01.65.Q 1. Introduction The most common and universal method in theories of fluids is the perturbation approach. Nonetheless, without analytic results of the integral equation theory for simple systems the per- turbation expansion could not be accomplished. To be specific, the analytic results of the mean- spherical approximation (MSA) for the fluid of hard spheres (for which the MSA is identical to the Percus-Yevick theory) [1, 2] and the Yukawa fluid [3, 4] are used to describe the properties of the reference system. In addition to these conventional MSA results there is also a perturbed version of the MSA called the first-order mean-spherical approximation (FMSA) [5]. When compared to computer simulation data for the Yukawa fluid model, the FMSA theory has been shown to be reasonably good but still less accurate than the common MSA [6]. However, an advantage of the FMSA theory, hereinafter referred to as FMSA/HS since it employs the hard sphere (HS) fluid for the reference, is its simplicity and transparency. These are the factors causing its increasing attractiveness and potential applicability in the liquid state theory and also, in general, in the soft condensed matter theory. The most important ideas of the FMSA/HS theory are outlined in the following section. The goal of this paper is to present a modification of the FMSA/HS theory. The modification pursues the idea that there is an alternative to the HS reference, namely a mean field theory based on a Yukawa fluid reference [7]. In section 3 the possible choices in this direction are considered. The results and discussions are the subject of section 4. The study is summarized in section 5. 2. Conventional FMSA/HS theory The FMSA/HS is a theory based on the solution of the Ornstein-Zernike (OZ) equation, h̃ (k) = c̃ (k) + ρh̃ (k) c̃ (k) (1) c© S. Hlushak, A. Trokhymchuk, I. Nezbeda, 2011 33004-1 http://dx.doi.org/10.5488/CMP.14.33004 http://www.icmp.lviv.ua/journal S. Hlushak, A. Trokhymchuk, I. Nezbeda by employing the perturbation expansions limited to the first order correction term, h̃ (k) = h̃o (k) + ∆h̃ (k) , c̃ (k) = c̃o (k) + ∆c̃ (k) , Q̂ (k) = Q̂o (k) + ∆Q̂ (k) , (2) for the total, h, and direct, c, correlation functions combined with Baxter factorization function, Q, respectively. Here and in what follows all symbols with the tilde will denote the three-dimensional Fourier transforms, while all symbols with the hats will denote the one-dimensional Fourier trans- forms or the Laplace transforms. The terms with subscript “o” denote the contribution of a reference system, and ∆h, ∆c, ∆Q are the respective corrections. The FMSA/HS theory has been developed by Tang and Lu [5, 8, 9] for the hard-core based fluid models u (r) = { ∞, r < σ, φ(r), r > σ, (3) with σ being the diameter of a hard core and φ(r) being the potential function for out-of-core pair interaction between molecules. In their studies Tang and Lu were considering the fluid of hard spheres (HS) as the reference (or unperturbed) system for which all three functions that enter expansions (2) are well known. In particular, for the total correlation function of the HS reference system, ho ≡ hHS , Tang and Lu used the expression ĝo (s) ≡ ĝHS (s) = L (sσ) e−sσ (1− η) 2 Q̂HS (sσ) s2 , (4) that is the Laplace transform of the radial distribution function gHS(r) = hHS(r)+1 resulting from the solution of the OZ equation for HS fluid in the Percus-Yevick approximation [1]. The function Q̂HS (t) refers to the Laplace transform of the Baxter factorization function of the HS fluid, Q̂o (t) ≡ Q̂HS (t) = S (t) + 12ηL (t) e−t (1− η) 2 t3 , (5) where S (t) = (1− η)2 t3 + 6η (1− η) t2 + 18η2t− 12η (1 + 2η) , (6) L (t) = ( 1 + η 2 ) t+ 1 + 2η (7) and η = πρσ3/6 is the packing fraction. The main result obtained by Tang and Lu concerns an expression for correction ∆h to the total correlation function (or radial distribution function) resulting from the out-of-core attraction φ(r). It reads [5] ∆ĥ (k) = P (ikσ) Q̂2 HS (ikσ) , (8) where P (ikσ) = ∆U (k) 2Q̂2 HS (−ikσ) − e−ikσ 2iπ ∞ ∫ −∞ ∆U (y) eiyσ (y − k) Q̂2 HS (−iyσ) dy , (9) with function ∆U (k) defined as ∆U (k) = ∞ ∫ σ r∆c (r) e−ikrdr. (10) In accordance with this approach the necessary closure to be used in equation (10) reads as ∆cFMSA/HS (r) = −βφ (r) , for r > σ, (11) 33004-2 Improved MSA theory for simple fluids with β = 1/kBT and T being the temperature. In general, without any relation to the above version of the FMSA theory, a key point of its application to particular fluid model [determined by the out-of-core interaction potential φ(r)] concerns an evaluation of the integral in equation (9). So far, within the framework of the FMSA/HS theory this has been done for the Yukawa, Lennard-Jones, Kihara and sticky fluids by Tang and Lu [8] and for the square-well fluids by Tang and Lu [9] and by Hlushak et al. [10]. In the particular case of the Yukawa (Y) fluid model, φY (r) = −ǫσ e−z(r−σ) r , (12) for which the function ∆U (k) ∼ e−ikσ, an integration contour in the right-hand side of equation (9) can be closed in the upper complex half-plane and evaluation of the integral is rather simple and does not require the calculation of the residues at zeroes of the function Q̂HS (−ikσ). Then, the Laplace transform of the correction term (8) for the total correlation function of the Yukawa fluid within the FMSA/HS theory reads ∆ĥFMSA/HS (s) = βǫσe−sσ (s+ z) Q̂2 HS (sσ) Q̂ 2 HS (zσ) . (13) Once the radial distribution function ĝ = ĝo+∆ĥ is known, the thermodynamics of the system can be calculated. 3. FMSA theory based on Yukawa reference system The first order mean-spherical approximation theory that we are dealing with is a kind of perturbation theory approach. Usually, within the perturbation theory, in order to improve the performance of the first order approximation one should calculate the second order correction term. As an alternative to this common way here we propose to qualitatively modify the reference system over which the perturbation is calculated, and as a result to continue working within the same first order approximation. This is of particular importance since our aim is to keep the simplicity and transparency of the improved FMSA theory at the level established by the conventional FMSA/HS approach. 3.1. Division of the potential Similar to the classical perturbation theory approach we proceed by dividing the initial inter- action potential u(r) into two parts u (r) = uo (r) + ∆u (r) , (14) that are the reference and residual contributions, respectively. However, for the reference system potential uo (r) we will require that besides the hard-core repulsion it should (i) include a piece of the attractive tail of the same strength ǫ as the total potential, and (ii) extend over the range that is somewhat shorter than the range of the total interaction. The function that allows us to satisfy and control these requirements in an easy and a natural way is the Yukawa (Y0) potential. Thus, we may define the desired reference fluid as follows: uo (r) ≡ uY0(r) = { ∞, r < σ, −εσe−zo(r−σ)/r, r > σ . (15) The improved FMSA closure, referred to as FMSA/Y0, reads as follows: ∆cFMSA/Y0 (r) = −βφ (r)− βǫσe−zo(r−σ)/r, for r > σ , (16) where zo is the so far undefined parameter. 33004-3 S. Hlushak, A. Trokhymchuk, I. Nezbeda 3.2. Yukawa reference system We note here that some authors have already suggested to utilize the non-HS reference system in liquid state theory (e.g., see [11, 12]). Moreover, the short-range attractive Yukawa fluid has been already considered as an alternative to the HS reference system by Melnyk et al. [7, 13–16] in their studies within the framework of an augmented van der Waals theory for simple fluids. It is very important that, like in the case of the HS reference potential, all properties for the Yukawa reference potential (15), including the Baxter function Qo ≡ QY0 , are available in the literature. They can be obtained either within the MSA theory for Yukawa fluid [e.g., see Blum and Hoye [4], Kalyuzhnyi et al. [17, 18]] or within the conventional FMSA/HS theory due to Tang and Lu as it is described in previous section 2. The MSA theory is more accurate than FMSA/HS theory, but the latter is simpler. For the purpose of the present study we decided to sacrifice the accuracy in order to maintain simplicity and transparency. Thus, in what follows the FMSA/HS theory will be employed to describe properties of the Y0 reference system. And, as we will see hereinafter, despite a less accurate description of the reference system, the FMSA/Y0 theory still shows a notable improvement against the FMSA/HS theory. Following the conventional FMSA/HS approach, the Laplace transform of the radial distribution function of the Y0 reference system reads ĝo (s) ≡ ĝY0 (s; zo) = ĝHS (s) + βǫσe−sσ (s+ zo) Q̂2 HS (sσ) Q̂ 2 HS (zoσ) , (17) where the first term corresponds to the contribution of a hard-sphere repulsion, while the second one is the contribution due to the short-range attraction attributed to the Y0 reference system in accordance with equation (15). Similarly, the FMSA/HS result for the Baxter factorization function of the Y0 reference is Q̂o (t) ≡ Q̂Y0 (t; zoσ) = Q̂HS (t) + 12βεη zoσQ̂HS (−t) e−t − (t+ zoσ) Q̂HS (t) tzoσ (t+ zoσ) Q̂2 HS (zoσ) . (18) 3.3. Full system within the FMSA/Y0 theory After all properties of interest for the Y0 reference fluid are specified, we turn to the entire system determined by the interaction potential u(r) or more precisely by the out-of-core potential function φ(r). Although this function can be substituted by any potential function used in literature to represent the simple fluids [e.g., Lennard-Jones, Sutherland, Yukawa or Kihara potentials, etc.] for the purpose of present study we proceed with the Yukawa (Y) fluid model already defined according to equation (12). By introducing Yukawa potential φY(r) into the FMSA/Y0 closure (16) we note that for both functions that enter ∆cFMSA/Y0(r), the function ∆U (k) is proportional to e−ikσ and the integral in equation (9) can be easily evaluated in the way it was already discussed in Introduction section. Then, the correction term ∆ĥFMSA/Y0 that is necessary to evaluate the radial distribution function of the full system, ĝ (s; zo, z) = ĝo (s; zo) + ∆ĥ (s; zo, z) , (19) is given by ∆ĥFMSA/Y0 (s; zo, z) = − βǫσe−sσ (s+ zo) Q̂2 Y0 (sσ; zoσ) Q̂ 2 Y0 (zoσ; zoσ) + βǫσe−sσ (s+ z) Q̂2 Y0 (sσ; zoσ) Q̂ 2 Y0 (zσ; zoσ) . (20) The calculations of thermodynamics within the FMSA/Y0 approach can be made through the energy route, i.e., in the way that is quite similar to how it was done by Tang et al. within the 33004-4 Improved MSA theory for simple fluids conventional FMSA/HS theory [6]. It consists in evaluating the internal energy using its definition through the radial distribution function, U NkT = 2πρβ ∞ ∫ 0 drr2g (r; zo, z)u (r) = 12ηβεσezσ ĝ (z; zo, z) . (21) This result is used to derive the Helmholtz free energy in the form A−Aid NkT = a0 + a1 + a2 , (22) where a0 = 4η − 3η2 (1− η) 2 , (23) a1 = −12η βεL (zσ) (1− η) 2 Q̂HS (zσ) (zσ) 2 , (24) a2 = −6ηβ2ε2 [ 1 (zσ + zoσ) Q̂2 HS (zσ) Q̂ 2 HS (zoσ) − 1 (zσ + zoσ) Q̂2 Y0 (zσ; zoσ) Q̂ 2 Y0 (zoσ; zoσ) + 1 2zσQ̂2 Y0 (zσ; zoσ) Q̂ 2 Y0 (zσ; zoσ) ] . (25) In the limit zo = z the above formulas reduce to the conventional FMSA/HS theory by Tang and Lu [6]. The chemical potential and pressure are obtained employing standard thermodynamic relations. 4. Results and discussions 4.1. How much of the attraction should be treated as the refer ence? 0.0 0.1 0.2 0.3 0.4 0.5 �� 3 0 2 4 6 8 10 z o ,m in z=0.5 z=1 z=1.8 z=3 z=4 Figure 1. Density dependence of the Y0 refer- ence decay parameter zo,min that minimizes the free energy for Yukawa fluid models with differ- ent value of the decay exponent z as it is spec- ified for each curve. One can see that at low density zo,min → 2z. The above outlined FMSA/Y0 theory con- tains the parameter zo that determines the range of reference attraction and needs to be specified. The most straightforward approach to proceed deals with the minimization of the free energy (22) of the Yukawa system over the range of possible values of zo . Noting that zo enters only a2 contribution, we deduce that the outcome of the minimization, i.e., zo,min , de- pends only on the density of the fluid but not on the temperature. The corresponding density dependence of zo,min for several Yukawa fluids is presented in figure 1. It is evident that in the limit of a vanishing density and for all systems considered, the parameter zo,min tends to the value of 2z. This rather remarkable result for choosing the Y0 reference system in such a strict form is quite probably limited to the FMSA theory only. Nevertheless, it is quite consistent with our earlier findings for the Y0 reference decay parameter within the framework of the Yukawa based van der Waals theory for the simple fluids [7, 13–16]. First of all, there is a requirement for 33004-5 S. Hlushak, A. Trokhymchuk, I. Nezbeda zo,min to be larger than z providing in this way the range of the reference attraction which will be shorter than that in the parent fluid. Secondly, it is obvious that the reference fluid should be in a one-phase region for the set of density and temperature parameters used in the studies of the parent fluid. In the case of a hard sphere reference, this requirement is satisfied since HS fluid does not exhibit the liquid/vapor transition. By adding an attraction one provides the possibility for the phase coexistence. However, the critical temperature is always getting lower if the range of attraction is shorter. In [14] we presented a collection of the Monte Carlo generated liquid/vapor envelopes for Yukawa fluid with different values of the decay parameter (see figure 1 in [14]) from which it follows that the values 3 < zoσ < 6 can be used as a reference decay for the Lennard- Jones-like Yukawa fluid (zσ = 1.8) since the critical point temperature in such a reference fluid will always be lower than the triple point temperature in parent fluid. On the other hand, analyzing the dependence of the critical point coordinates of the Lennard- Jones-like Yukawa fluid (zσ = 1.8) on the reference system decay parameter zo , in [13] we have shown that the best agreement with computer simulation data can be reached if 3 < zoσ < 4. Thus, the result zo = 2z agrees with these findings and in all our subsequent calculations for the Lennard-Jones-like Yukawa fluid we impose zoσ = 3.6, making it independent of the fluid density. 4.2. Comparison between FMSA/HS and FMSA/Y0 solutions To illustrate the improvements that the FMSA/Y0 theory yields over the conventional FMSA/HS theory, we compare the predictions made up by these two theoretical approaches for compressibility factor and vapor/liquid phase diagram of the most popular Yukawa fluid model defined by the decay parameter zσ = 1.8. Figure 2. Isotherms T ∗ = 1.5 and T ∗ = 1 (from the top to the bottom) for compressibil- ity factor, Z = PV/NkT , of the Yukawa fluid with zσ = 1.8 as obtained from the conven- tional FMSA/HS theory (black dashed lines) and from the improved FMSA/Y0 theory with zoσ = 3.6 (red solid lines). The symbols denote computer simulation data by Schukla [19]. Figure 3. Vapor-liquid phase diagram of the Yukawa fluid with zσ = 1.8. The notations are the same as in figure 2. The compressibility factors, Z = PV/NkT , of the Yukawa fluid with zσ = 1.8 and for two temperatures, T ∗ = 1.5 and 1, that result from both FMSA approaches are presented in figure 2. To estimate the effectiveness of the theory, the computer simulation data due to Schukla [19] are shown as well. As it was already pointed out by Tang [6], the conventional FMSA/HS theory does a good job in predicting the compressibility factor of the Yukawa fluid. Nevertheless, the proposed 33004-6 Improved MSA theory for simple fluids FMSA/Y0 theory shows to be even more accurate, especially at lower temperature and in the region of intermediate densities. While the enhancement may seem to be minor in the case of compressibility factor, it is more pronounced for the vapor/luqiud phase diagram shown in figure 3. The results of the FMSA/Y0 theory lead to the shrinkage of coexisting densities envelope in the region near the critical point, approaching the computer simulation data. The estimated improvement in comparison with the FMSA/HS predictions is around 30%. 5. Conclusions Being compared with computer simulation data for the Lennard-Jones-like Yukawa fluid (zσ = 1.8), the first order mean-spherical approximation (FMSA/HS) due to Tang et Lu has shown to be only slightly less accurate in the calculations of the thermodynamics and liquid/vapor phase coexistence than the full MSA theory (details of this comparison and a corresponding discussion can be found in the [6]). At the same time, the FMSA/HS is a much simpler theory. Abbreviation HS underlines here that the theory is based on the hard-sphere reference system. In the present study we reported a further improvement of this approach, which we called a conventional FMSA/HS theory by introducing the short-ranged Yukawa fluid (with decay parameter zo = 2z) as a new reference system. Consequently, we are referring to this theory as the FMSA/Y0 theory. Numerical calculations that we performed for the same model discussed by Tang [6] showed that the proposed modifications of the reference system make the first order mean-spherical approximation even more accurate. It is important that in order to treat the novel reference we are employing the conventional FMSA/HS approach. Due to this, the level of simplicity and transparency of the improved FMSA/Y0 theory is kept at the same level as that of the conventional FMSA/HS theory. In particular, it is easy to see that the resulting new expressions for main ingredients of the FMSA ideology – the Laplace transforms of the radial distribution function, equation (17), and the Baxter factorization function, equation (18), of the new reference system, as well as the equation (20) for the correction term – are only slightly longer than their original counterparts [see equations (4), (5) and (13)] and both being composed of the same variables, functions and model parameters. Only the parameter of the FMSA/Y0 theory appears to be the decay parameter zo, for which in the case of Yukawa fluid we obtained zo = 2z, and which determines an amount of the out-of-core attraction that should be attributed to the Y0 reference system. Acknowledgement This work was supported by the Grant Agency of the Academy of Sciences of the Czech Republic (Grant No. IAA400720710) and the Czech-Ukrainian Bilateral Cooperative Program. References 1. Wertheim M.S., Phys. Rev. Lett., 1963, 10, 321; doi:10.1103/PhysRevLett.10.321. 2. Baxter R.J., Aust. J. Phys., 1968, 21, 563; doi:10.1071/PH680563. 3. Waisman E., Mol. Phys., 1973, 25, 45; doi:10.1080/00268977300100061. 4. Hoye J. S., Blum L., J. Stat. Phys., 1977, 16, 399; doi:10.1007/BF01013184. 5. Tang Y., Lu B. C.-Y., J. Chem. Phys., 1993, 99, 9828; doi:10.1063/1.465465. 6. Tang Y., J. Chem. Phys., 2003, 118, 4140; doi:10.1063/1.1541615. 7. Nezbeda I., Melnyk R., Trokhymchuk A., J. Supercritical Fluids., 2010, 55, 448; doi:10.1016/j.supflu.2010.10.041. 8. Tang Y., Lu B. C.-Y., J. Chem. Phys., 1994, 100, 3079; doi:10.1063/1.466449. 9. Tang Y., Lu B. C.-Y., J. Chem. Phys., 1994, 100, 6665; doi:10.1063/1.466449. 10. Hlushak S., Sokolowsli S., Trokhymchuk A., J. Chem. Phys., 2009, 130, 234511; doi:10.1063/1.3154583. 11. Sowers G.M., Sandler S.I., Fluid Phase Equlibria., 1991, 63, 1; doi:10.1016/0378-3812(91)80017-P. 12. Kahl G., J. Non-Cristalline Solids, 1990, 117, 104; doi:10.1016/0022-3093(90)90889-T. 33004-7 http://dx.doi.org/10.1103/PhysRevLett.10.321 http://dx.doi.org/10.1071/PH680563 http://dx.doi.org/10.1080/00268977300100061 http://dx.doi.org/10.1007/BF01013184 http://dx.doi.org/10.1063/1.465465 http://dx.doi.org/10.1063/1.1541615 http://dx.doi.org/10.1016/j.supflu.2010.10.041 http://dx.doi.org/10.1063/1.466449 http://dx.doi.org/10.1063/1.466449 http://dx.doi.org/10.1063/1.3154583 http://dx.doi.org/10.1016/0378-3812(91)80017-P http://dx.doi.org/10.1016/0022-3093(90)90889-T S. Hlushak, A. Trokhymchuk, I. Nezbeda 13. Melnyk R., Moucka F., Nezbeda I., Trokhymchuk A., J. Chem. Phys., 2007, 127, 094510; doi:10.1063/1.2766937. 14. Melnyk R., Nezbeda I., Henderson D., Trokhymchuk A., Fluid Phase Equlibria., 2009, 279, 1; doi:10.1016/j.fluid.2008.12.004. 15. Melnyk R., Orea P., Nezbeda I., Trokhymchuk A., J. Chem. Phys., 2010, 132, 134504; doi:10.1063/1.3371710. 16. Melnyk R., Nezbeda I., Trokhymchuk A., Mol. Phys., 2011, 109, 113; doi:10.1080/00268976.2010.542034. 17. Kalyuzhnyi Yu. V., Blum L., Rescic J., Stell G., J. Chem. Phys., 2000, 113, 1135; doi:10.1063/1.481892. 18. Hlushak S., Kalyuzhnyi Yu. V. Preprint arXiv:cond-mat.soft/0805.0688v1, 2008. 19. Schukla K. P., J. Chem. Phys., 2000, 112, 10358; doi:10.1063/1.481673. Удосконалене середньосферичне наближення першого порядку для простих рiдин С. Глушак1, А. Трохимчук1,2, I. Незбеда3 1 Iнститут фiзики конденсованих систем НАН України, вул. I. Свєнцiцького, 1, 79011 Львiв, Україна 2 Факультет хiмiї та бiохiмiї, Унiверситет Бригам Янг, Прово, США 3 Факультет природничих наук, Ун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 як i звичайне ССНПП/ТС, але є бiльш точною. Ключовi слова: рiдина Юкави, середньосферичне наближення, теорiя збурень, моделювання Монте-Карло 33004-8 http://dx.doi.org/10.1063/1.2766937 http://dx.doi.org/10.1016/j.fluid.2008.12.004 http://dx.doi.org/10.1063/1.3371710 http://dx.doi.org/10.1080/00268976.2010.542034 http://dx.doi.org/10.1063/1.481892 http://dx.doi.org/10.1063/1.481673 Introduction Conventional FMSA/HS theory FMSA theory based on Yukawa reference system Division of the potential Yukawa reference system Full system within the FMSA/Y0 theory Results and discussions How much of the attraction should be treated as the reference? Comparison between FMSA/HS and FMSA/Y0 solutions Conclusions

Improved first order mean spherical approximation for simple fluids

Репозитарії

Схожі ресурси