Toward a Quantitative Theory of Ultrasmall Liquid Droplets and VaporsLiquid Nucleation
|
|
- Phebe Parker
- 5 years ago
- Views:
Transcription
1 4988 Ind. Eng. Chem. Res. 2008, 47, Toward a Quantitative Theory of Ultrasmall Liquid Droplets and VaporsLiquid Nucleation Zhidong Li and Jianzhong Wu* Department of Chemical and EnVironmental Engineering, UniVersity of California, RiVerside, California Thermodynamic properties of small systems can be drastically different from those corresponding to their macroscopic counterparts due to the surface and fluctuation effects. While the equations of state for macroscopic systems are well advanced, quantitative predictions of the structural and thermodynamic properties of small systems from a molecular perspective remain a daunting challenge. This article illustrates applications of a nonmean-field density functional theory to two types of ultrasmall liquid droplets: one is stabilized in a container of finite size and the other is unstable as appeared during vaporsliquid nucleation. For small systems of simple fluids represented by the Lennard-Jones model, theoretical predictions are compared with results from molecular simulations for the microscopic structure, the droplet size, and the free energy of formation over a broad range of conditions. The numerical agreement of theory with simulation data is comparable to that for the corresponding macroscopic systems. While the Tolman length, a correlation of curvature on surface tension, is negligible at least for small droplets of simple fluids, the vaporsliquid interfacial tension declines with the droplet size approximately proportional to the Gaussian curvature. Surprisingly, the Laplace equation for pressure change across a curved surface remains accurate even for a liquid droplet with the radius only a few times the molecular diameter. I. Introduction Small systems are ubiquitous in nature and in modern technology. 1,2 Common examples include cellular compartments, molecular motors, and viral nucleocapsids in biology, various nanoassemblies of amphiphilic molecules such as micelles, vesicles, and microemulsions, and solidlike metal or semiconductor clusters or quantum dots. Small systems may also emerge as unstable or metastable transition states during various stages of phase transitions in macroscopic systems such as vaporsliquid or liquidssolid nucleation. Whereas it has been well documented that the thermodynamic properties of small systems can be drastically different from those corresponding to their macroscopic counterparts and they may vary at different environments (open, closed, isothermal, isobaric, etc.), it is not until recently that the properties of small systems have become amenable to direct experimental measurements. In contrast to statistical thermodynamics for macroscopic systems, few theoretical tools are available for quantitative analysis of the properties of small systems in terms of their chemical constituents, the intermolecular forces, and interactions with the surrounding environment. In this work, we intend to apply a density functional theory (DFT) to represent the properties of ultrasmall liquid droplets quantitatively. For direct comparison with simulation results, we are concerned here only with a simple liquid model where the intermolecular force can be represented by the Lennard-Jones (LJ) potential. It has been long recognized that for macroscopic systems specified by constant temperature and chemical potential, a liquid droplet is thermodynamically unstable relative to a bulk state. Because of the fractionation effect, however, the same liquid droplet can be stabilized in a closed system of finite size. Toward a better understanding of vaporsliquid transition, formation of liquid droplets in a supersaturated vapor has been scrutinized by extensive molecular simulations. 3 8 While the properties of the critical nucleus, a tiny unstable liquid droplet, * To whom all correspondence should be addressed. jwu@ engr.ucr.edu. Tel.: Fax: have been examined by molecular simulations in great details, theoretical interpretation of the simulation data often relies on the classical nucleation theory (CNT), 9 14 proposed by J. W. Gibbs over 100 years ago. Without any specific knowledge on the microscopic structure of the critical nucleus, CNT presumes that the thermodynamic properties of a tiny liquid droplet are identical to those corresponding to a macroscopic liquid system. Although CNT captures qualitatively the free-energy barrier of a vaporsliquid phase transition far from the spinodal line, its numerical performance is often unsatisfactory and deteriorates at conditions where the properties of the droplet deviate from those of the bulk phase. Applications of DFT to vaporsliquid nucleation have been well documented. In particular, DFT has been used to describe the interfacial structure, nucleus size, surface tension, free energy of nucleation, and kinetic and mechanical properties In comparison to CNT, DFT surpasses a number of unrealistic approximations including the assumption of a sharp vapors liquid interface and utilization of the bulk interfacial tension. In addition, DFT accounts for the microscopic structure of the critical nucleus explicitly. Regrettably, most previous applications of DFT rely either on the local density approximation (LDA) for the free energy due to the short-range intermolecular repulsion and/or on the van der Waals-like mean-field approximation (MFA) for longer-range intermolecular attractions. Because the MFA is not accurate for the thermodynamic properties of bulk phases, the performance of an MFA-based DFT is at most semiquantitative. In this work, we intend to improve the numerical performance by using a modified fundamental measure theory (MFMT) 22,23 for the short-range repulsion and a quadratic density expansion (QDE) 24,25 for longrange attractions. In our previous publications, 26,27 we have demonstrated that the nonmean-field DFT provides accurate predictions of the properties of bulk phases and planar interfacial tensions. Because DFT is directly applicable to both bulk and inhomogeneous systems on the basis of the same set of /ie070578i CCC: $ American Chemical Society Published on Web 07/20/2007
2 molecular parameters, we expect that it remains reliable for representing the properties of small systems. A second purpose of this work is to examine the interfacial properties of a liquid droplet in terms of its size. The size or curvature-dependence of surface tension is conventionally described by the Tolman equation, 28 which provides a firstorder correction of the surface tension of a liquid droplet with respect to that in the planar limit. The Tolman length can be defined as the linear coefficient in a Taylor expansion of the droplet surface tension as a function of the curvature; it can also be related to the separation between locations of the equimolar dividing surface and of the thermodynamically defined surface in the planar limit. Whereas the Tolman length has been used extensively in surface chemistry, its magnitude and sign remain controversial even for simple fluids. For example, DFT predicts that the Tolman length is negative for LJ fluids and the magnitude is about 15% to 20% of the LJ diameter. 16,20 However, alternative theories for LJ fluids indicate that the Tolman length can be positive at similar conditions. 29 On the basis of DFT calculations, McGraw and Laaksonen 30 asserted that the Tolman length for the LJ fluid is zero and the size dependence of surface tension is mainly attributed to the second-order correction of the curvature. A small Tolman length appears consistent with results from recent molecular simulations. 7,31 Vaporsliquid transition is a dynamic process that involves growing of a liquid droplet that cannot be directly addressed by equilibrium DFT. To represent the free-energy landscape of nucleation, Talanquer and Oxtoby extended DFT to unstable systems. 19,20 The finite size effect is accounted for by assuming that the liquid droplet in contact with a vapor has a density equal to that at the boundary of the container. The density profile and subsequently thermodynamic properties of the critical nucleus can be obtained by using a pseudo grand canonical ensemble. With this method, the properties of an arbitrary liquid droplet can be obtained by inserting different stable clusters into the same supersaturated vapor and by assuming that the internal structures of these stable clusters are fixed. While the pseudo grand canonical ensemble is efficient for calculation of the free-energy profile of vaporsliquid nucleation, the hypothetical boundary condition is sometimes difficult to justify, i.e., the stable cluster may exist in a very small container where the surrounding vapor density is not equal to that at the boundary of the container. Besides, it remains to be validated whether the internal structure of stable clusters depends on the properties of the coexisting vapor. In this work, we propose a new method to obtain the nucleation free-energy landscape directly from numerical iteration. This method has been tested by predicting the density profiles of nonequilibrium nuclei. II. Density Functional Theory We consider ultrasmall liquid droplets of a LJ fluid surrounded by a supersaturated vapor. The system is thermodynamically stable in the canonical ensemble where the total number of molecules is fixed, but it is unstable in the grand canonical ensemble that allows fluctuation of the overall number density. As in Monte Carlo simulations, 7,8 the intermolecular potential is represented by a truncated-and-shifted LJ potential u(r) ) { ulj (r) - u LJ (r c ), r e r c 0, r > r c } (1) where u LJ (r) is u LJ (r) ) 4ɛ [( σ r ) 12 - ( σ r ) 6 ] (2) Ind. Eng. Chem. Res., Vol. 47, No. 15, and σ and ɛ are the size and energy parameters, respectively. Throughout this work, the cutoff distance r c is fixed at 2.5σ. To apply DFT to the LJ fluid, we divide u LJ (r) into a shortrange repulsion and a long-range attraction as introduced originally in the BakersHenderson (BH) perturbation theory 32 u rep (r) ) { ulj (r), r e σ 0, r > σ } (3) u att (r) ) { 0, r e σ u LJ (r), r > σ } (4) For the truncated-and-shifted LJ potential, the intermolecular energy also includes a residual interaction potential u res (r) ) { -ulj (r c ), r e r c -u LJ (r), r > r c } (5) The repulsive part of the LJ potential can be approximated by a hard-sphere (HS) potential with an effective diameter T / d ) T / T /2σ (6) where T* ) k B T/ɛ with k B being the Boltzmann constant and T the absolute temperature. Alternatively, the LJ potential can be divided into a repulsive part and an attractive part following the WeekssChandlersAnderson (WCA) theory. 34 The nonmean-field DFT (NMF-DFT) for the LJ fluids has been reported in our previous publications. 27,35 37 Briefly, the intrinsic Helmholtz functional F is decomposed into four components F ) F id + F ex hs + F ex ex att + F res (7) where F id corresponds to that for a system of ideal-gas molecules, F ex hs represents contribution to the Helmholtz energy due to the HS repulsion, F ex att is that due to van der Waals attraction, and F ex res arises from the residual part of the intermolecular potential. The ideal-gas term is known exactly F id ) k B T dr F(r)[ln F(r) - 1] (8) The excess Helmholtz energy due to the HS repulsion F ex hs can be described accurately by a modified fundamental measure theory (MFMT) 22,23 F ex hs ) k B T Φ hs [n R (r)] dr (9) in which the reduced excess energy density Φ hs depends on six weighted densities n R (r) ) F(r )ω (R) ( r - r ) dr derived on the basis of hard-sphere geometry 38 Φ hs )-n 0 ln(1 - n 3 ) + n 1 n 2 - n V1 n V n π[ 2] n 3 ln(1 - n 3 ) + n 2 3 (n 3 2-3n 2 n V2 n V2 ) (10) 3 (1 - n 3 ) n 3 For the LJ potential, the effective HS diameter d instead of the LJ diameter σ should be used to calculate the weighted densities. 35,36 The excess Helmholtz energy due to the van der Waals attraction F ex att is represented by a quadratic density expansion 15,27,39 F ex att ) dr F(r)f att [n 3 (r)] + k B T 4 dr dr C att ( r - r, F)[F(r) -F(r )] 2 (11)
3 4990 Ind. Eng. Chem. Res., Vol. 47, No. 15, 2008 with Fj )[F(r) +F(r )]/2. The attractive Helmholtz energy per molecule in the bulk f att b depends on the packing density η ) πf b d 3 /6, which is replaced by n 3 (r) for inhomogeneous systems. att Analytical expressions for f b and for the direct correlation function (DCF) C att (r) are derived from the first-order meanspherical approximation (MSA). 26,35 37 Finally, the Helmholtz energy due to the weak but long-ranged residual potential is accounted for simply by using the mean-field approximation F ex res ) 1 dr dr u res ( r - r )F(r)F(r ) (12) 2 A typical mean-field theory for LJ fluids applies eq 12 to the entire attractive part of the intermolecular potential defined by the WCA perturbation theory. 34 In this case, the attractive part of the Helmholtz energy is given by F ex att ) 1 dr dr u att ( r - r )F(r)F(r ) (13) 2 with u rep (r) ) { ulj (r) + ɛ, r e σ 0, r > σ }, u att (r) ) { -ɛ, r e σ u LJ (r), r > σ } (14) To facilitate a comparison between mean-field and nonmeanfield approximations, we also perform calculations by using the DFT with mean-field approximation (MFA-DFT) for the attractive potential. For consistency, F ex hs is represented by MFMT for both cases. In MFA-DFT, the effective HS diameter is calculated from d ) dr{1 - exp[-u rep (r)/k B T]}. At equilibrium, the density profile satisfies the Eulers Lagrange equation F(r) ) exp{βµ - δf(r)} δβfex (15) where β ) 1/k B T. In eq 15, the chemical potential µ is specified by the coexisting bulk phase for an open system or for a closed system with N particles, it is replaced by the normalization condition F(r) ) N exp{- δf(r)} δβfex dr exp{- δf(r)} δβfex (16) Equations 15 and 16 are solved iteratively by using the Picard method with a step length 0.02σ. The initial guess for the density profile at the vaporsliquid interface is provided by a step function, F(r) )F L for r e r 0 and F(r) )F V for r > r 0, where r 0 is the radius of the liquid droplet. III. Results and Discussion A. Bulk Properties. While the limitations of a mean-field theory for representing the thermodynamic properties of a bulk fluid are well recognized, much less known is its reliability for inhomogeneous systems. In fact, a theory for inhomogeneous systems including DFT must reduce to an equation of state for bulk systems when the local density profile F(r) is replaced by the average bulk density F b. Therefore, we can test the numerical performance of different versions of DFT by comparison with the simulation results in the bulk limit (even though such comparison does not guarantee that a theory good for bulk fluids is necessarily accurate at inhomogeneous conditions). Figure 1 shows the phase diagrams for truncated-and-shifted LJ fluids predicted by our nonmean-field theory and by the mean-field theory according to BH or WCA separation of the attractive Figure 1. Comparison of phase diagrams and saturation pressures for the Lennard-Jones fluid predicted from the mean-field and nonmean-field theories and from Monte Carlo simulation. 42,44 In part a, the dashed lines represent the spinodal curve predicted from the nonmean-field theory. and repulsive LJ potentials. A comparison with the simulation results indicates that, except near the critical point, the nonmeanfield theory faithfully reproduces the coexisting vapor and liquid densities as well as the saturation pressure while predictions from the mean-field theory are at most semiquantitative. Evidently, the WCA separation of the LJ potential overestimates the intermolecular attraction, which is responsible for the small values of the coexisting vapor/liquid densities and the saturation pressure at a given temperature. On the other hand, the BH separation of the LJ potential underestimates the attraction, which leads to less accurate mean-field properties. The dashed lines in Figure 1a show the spinodal curve predicted by the nonmean-field theory. These lines set the boundaries of the metastable region and are important in study of vaporsliquid nucleation. B. Stable Liquid Droplets. We now consider the properties of a liquid droplet coexisting with its own vapor in a closed spherical container. The system was studied by Monte Carlo simulations. 8 Because of the finite size, the liquid droplet is thermodynamically stable and can be represented by a canonical ensemble (i.e., NVT ensemble) along with the coexisting vapor phase. As specified in Monte Carlo simulations, 8 we assume that the external potential due to the cage surface is represented by V ext (R B - r) ) { ulj (R B - r) + ɛ, R B - r e σ 0, R B - r > σ } (17)
4 Ind. Eng. Chem. Res., Vol. 47, No. 15, Figure 2. Density profiles of a stable liquid droplet coexisting with a vapor in a spherical cavity with N Lennard-Jones molecules. (a) T* ) 0.71 and (b) T* ) Solid and dashed lines are predictions of NVTsNMFsDFT and of NVTsMFAsDFT, respectively, and the symbols are simulation results from ref 8. where R B is the radius of the container, which is about 10 times the molecular diameter. Because the container is in direct contact with the vapor phase, the external potential has negligible effect on the properties of the liquid droplet. Figures 2 and 3 present the liquidsvapor interfacial density profiles predicted by mean-field and nonmean-field DFTs in comparison with the simulation results. 8 Partitioning of the LJ molecules among the coexisting liquid and vapor phases ensures that at a given temperature the liquid droplet grows with the total number of particles within the container or shrinks with temperature if the total number of particles is fixed. Because most of the volume is occupied by the vapor phase, theoretical predictions are extremely sensitive to its density, which must be calculated concomitant with the properties of the small liquid droplet. As expected, the performance of NMF-DFT is superior to the MFA-DFT, and in most cases, the nonmean-field theory is able to reproduce the simulation data quantitatively. In particular, the NMF-DFT works the best at low temperature and large droplet size. The discrepancies between NMF-DFT and simulation at high temperature are possibly due to the fact that a small deviation in the vapor density leads to a magnified error for the droplet density. In general, MFA-DFT overestimates the droplet size and liquid density for small droplets or at high temperature, and the trend is opposite for large droplets or at low temperature. Besides, the MFA-DFT yields oscillation of the liquid density at low temperature. Because WCA Figure 3. Same as Figure 2 but at different temperatures for (a) N ) 256 and (b) N ) overestimates the intermolecular attraction, this oscillation is possibly due to the fact that the droplets show some ordered structure. C. Unstable Liquid Droplets. In contrast to a finite system, a liquid droplet is thermodynamically unstable if it is surrounded by a bulk supersaturated vapor phase at constant pressure or chemical potential. In this case, coexistence of the liquid droplet with its surrounding vapor reflects a free-energy maximum. As a result, the liquid droplet may either grow or shrink depending on whether its size is larger or smaller than that of the critical nucleus. To calculate the properties of a unstable droplet, we follow a method originally proposed by Zeng and Oxtoby. 21 At a given temperature, the initial guesses of the densities for the liquid and vapor phases, F L, F V, are obtained from the densities of the equilibrium liquid F 0 L and that of the supersaturated vapor F n V. The saddle point is identified by trying different values of the droplet radius r 0, on the basis of the fact that a tiny change of r 0 will lead to either shrinkage or growth of the droplet size. Figure 4 shows variation of the free-energy barrier during iteration near the saddle point at temperature T* ) and supersaturation S ) Here, the supersaturation ratio is defined as S ) p n V /p 0 with p n V and p 0 being the pressure of the supersaturated vapor and that at equilibrium, respectively. The free-energy barrier is defined as Ω ) Ω[F(r)] + p n V V with V being the system volume. Figure 4 indicates that Ω bifurcates around r 0 /σ ) The droplet eventually shrinks when r 0 /σ ) 4.8 and grows when r 0 /σ ) Both curves exhibit a flat plateau lasting for several hundreds of iteration steps within
5 4992 Ind. Eng. Chem. Res., Vol. 47, No. 15, 2008 and the corresponding vapor spinodal density is Thus, the densities of the supersaturated vapors under study (F b / ) for S ) ) are much lower than that at the spinodal point. In the DFT calculations, the critical nucleation barrier is represented by Ω* and the critical excess number of particles n* is given by n ) 4π [F(r) -F n V ]r 2 dr (18) Figure 4. Nucleation barrier calculated from µvtsnmfsdft as a function of the iteration steps at T* ) and S ) The thin dashed line marks the critical nucleation barrier. Figure 5. (a) Free-energy barrier and (b) excess number of particles of a critical nucleus as a function of supersaturation at T* ) Solid and dashed lines represent predictions of µvtsnmfsdft and CNT, respectively, and the symbols are simulation data from ref 7. which the droplet structure remains nearly unchanged. This plateau marks the critical nucleation barrier Ω*. Figure 5 compares predictions from NMF-DFT with simulation results for vaporsliquid nucleation of the LJ fluid at T* ) in the NPT ensemble. 7 At this temperature, the coexistence vapor and liquid densities are and , respectively, As indicated by Lee et al., 17 Ω in the µvt ensemble is equal to G in the NPT ensemble. In the µvt ensemble, we consider the interactions between the molecules in the system and those in the reservoir explicitly. However, such interaction is absent in the NPT ensemble, which could be significant as the nucleus size increases. Figure 5 indicates that, far from the spindal point, agreement of the theory with the simulation is good for both Ω* and n*. The nucleation barrier falls monotonically with the degree of supersaturation S, which eventually leads to a spontaneous nucleation. Far away from the spinodal point, the droplet size falls with the degree of supersaturation, which is consistent with predictions of the Kelvin equation F 0 L k B T ln(s) ) 2γ/R where γ is the surface tension and R is the droplet radius. Approaching the spinodal point, however, the droplet size may exhibit a shallow minimum. 40,41 In this case, the Kelvin equation becomes inappropriate because at higher density the vapor phase cannot be treated as an ideal gas and the structure of the nucleus becomes very different from that of a macroscopic droplet. We did not calculate the droplet size at high supersaturation in this work because the low free-energy barrier makes the critical nucleus hardly discernible. For comparison, we have also calculated the properties of critical nuclei by using the classical nucleation theory (CNT). According to CNT, the critical nucleation barrier Ω* and droplet size n* are given by, respectively, Ω / ) 16πγ 3 (19) 3 p 2 n / ) 32πγ 3 3 p (F L 3 0 -F V n ) (20) where γ is the surface tension of a planar vaporsliquid interface and p ) p L n - p V n with p L n being the pressure of a hypothetical liquid droplet that has the same chemical potential as that of the vapor. By applying the NMF-DFT to a planar interface, we obtained that γ / ) γ σ 2 /ɛ ) 0.527, which is close to the simulation result In consistency with the previous work, 7,41 CNT predicts the size of critical nuclei accurately. For the free-energy barrier, the performance of CNT is comparable to that of NMF-DFT at low supersaturation. But, contrary to NMF-DFT, the accuracy of CNT deteriorates as the degree of supersaturation increases. Figure 6 shows the density profiles of the critical nuclei at three values of supersaturation (S ) 1.53, 1.79, and 2.04) and T* ) The lines are predicted by the NMF-DFT, and the symbols are the simulation results. 7 Unlike that for closed systems, the performance of the DFT is near-perfect. For sufficiently large droplet (S ) 1.53), the density of the liquid core is close to that of the bulk liquid at the same chemical potential of the surrounding vapor. In that case, both the size of the critical droplet and liquid density falls as the vapor pressure increases. D. Tolman Length. The Tolman length measures the extent by which the surface tension of a small liquid drop deviates from its planar value. It can be formally defined as the leading
6 Ind. Eng. Chem. Res., Vol. 47, No. 15, Figure 6. Density profiles for the critical nuclei at different degrees of supersaturation S at T* ) The thick lines and symbols correspond to predictions from µvtsnmfsdft and from simulation, 7 respectively. The thin horizontal lines represent the corresponding density in the bulk liquid that has the same chemical potential as that of the vapor. coefficient in the Taylor expansion of the surface tension in terms of the curvature: γ ) γ (1-2δ ) TL R + (21) where δ TL is the Tolman length and 1/R is the curvature of a spherical liquid droplet. At T* ) 0.7, the coexistence vapor and liquid densities are and , respectively. The vapor spinodal density is which is much higher than the density of supersaturated vapor used in the work (F / b ) for S ) ). By using the nonmeanfield DFT, we predict γ / ) 0.599, which agrees well with the simulation result To define the droplet size or curvature, we use an equimolar dividing surface to partition the liquid droplet and its surrounding vapor. The radius of the equimolar dividing surface is calculated from R 3 3 e ) V [F(r) -F V n ]r 2 3 n / dr ) F(0) -F n 4π[F(0) -F V n ] (22) In terms of R e and γ e, the free energy of nucleus formation is given by Ω / )- 4 3 πr e3 p + 4πR e 2 γ e (23) On the basis of the droplet radius R e and the free energy of formation Ω* calculated from DFT, we obtain the surface tension γ e. Subsequently, the Tolman length can be estimated by fitting γ e in terms of the Helfrich formula 43 γ e ) γ - 2δ TL γ + k s (24) R 2 e R e where k s is the rigidity coefficient. Equation 24 can be obtained from eq 21 by replacing R and γ by R e and γ e, respectively, and by neglecting the high-order terms. Figure 7a shows γ e calculated from eq 23 as a function of 1/R e at T* ) 0.7 and By fitting the data with eq 24, we estimate that δ TL /σ ) and βk s ) for T* ) 0.7, and δ TL /σ ) and βk s ) for T* ) Figure 7. (a) Surface tension γ e and (b) pressure drop p versus 1/R e for the critical nuclei at T* ) 0.7 and Symbols are from µvtsnmfsdft. The thick lines are obtained by the best fitting (a) of eq 24 and (b) of the Laplace equation p ) 2γ /R e, respectively. The thin horizontal line in part a marks the surface tension at the planar limit. Because the Tolman length δ TL is close to zero, the first-order correction to the surface tension appears insignificant at least for the LJ fluid. Our DFT calculations support an earlier conjecture by McGraw and Laaksonen assuming that the Tolman length is zero. 30 Figure 7a indicates that the surface tension rises with the droplet size and it reaches the planar limit asymptotically. Minimization of Ω* with respect to R e gives the generalized Laplace equation 30,43 p ) 2γ e R e + γ e R e ) 2γ R e - 2δ TL γ R e 2 (25) Because the Tolman length is close to zero, eq 25 suggest that the original Laplace equation is accurate even for a very small liquid droplet. Figure 7b shows an excellent linear relation as predicted by the original Laplace equation p ) 2γ /R e. The good fix also suggests that δ TL can be ignored for the LJ fluids. E. Free-Energy Profile of VaporsLiquid Nucleation. Different from the method of the pseudo grand canonical ensemble, we find that the nucleation free-energy landscape can be directly obtained from numerical iteration in the grand canonical ensemble. To illustrate, Figure 8 shows variation of the free energy for different values of r 0 and plots the freeenergy barrier Ω as a function of the excess number of particles n during the iteration at T* ) and S ) Except the initial iteration steps, the shrinkage (r 0 /σ e 4.8) or
7 4994 Ind. Eng. Chem. Res., Vol. 47, No. 15, 2008 Figure 8. Free-energy landscape obtained by using different initial guesses of r 0 at T* ) and S ) 1.53 calculated from µvtsnmfsdft. Figure 10. Density profiles for several precritical nuclei at T* ) and S ) Lines represent results from µvtsnmfsdft, and symbols are simulation data from ref 7. free-energy landscape during liquid droplet formation. Figure 9 shows that the energy landscape from the pseudo µvt ensemble is systematically lower than that from the exact result. The difference is probably due to the assumption that the internal structures of the stable clusters are independent of the supersaturated vapor where they are inserted. Figure 10 presents the density profiles of several precritical nuclei in comparison with the simulation results 7 at T* ) and S ) Except for the two smallest nuclei, the structures of all the other nuclei are faithfully reproduced by the NMF- DFT. The good agreement also supports the rationality of the new method for determining the free-energy landscape. Figure 9. Free-energy landscape β Ω at different degrees of supersaturation S at T* ) Solid and dashed lines are predictions of µvtsnmfsdft and pseudo µvtsnmfsdft, respectively. the growth (r 0 /σ g 4.82) of the nuclei follow a unique route independent of r 0. This unique pathway corresponds to the freeenergy change at different nucleation stages, i.e., the free-energy landscape. Figure 9 shows the free-energy landscape of vaporsliquid nucleation at T* ) and different degrees of supersaturation. Several essential features are captured by these results. First, the energy barrier is zero when small nuclei disappear in the vapor phase and approaches a negative infinity when the large nuclei become the liquid phase. Second, near the spinodal point, the free-energy landscape narrows and the critical barrier becomes smaller. We speculate that the energy barrier disappears at the spinodal point. For comparison, we have also carried out NMF-DFT calculations in a pseudo µvt ensemble as proposed by Talanquer and Oxtoby. 19,20 For these calculations, we encapsulate a fixed number of particles in an imaginary spherical container and the density of surrounding vapor is taken into account by assuming it to be identical to that of the container boundary. A set of stable clusters are obtained by minimization of the Helmholtz energy with respect to the container volume. By assuming that the internal structures of the stable clusters remain unchanged when they are inserted into a given supersaturated vapor, we are able to calculate the IV. Conclusions We have investigated the properties of ultrasmall liquid droplets and vaporsliquid nucleation for the Lennard-Jones fluids by using a nonlocal density functional theory. Unlike a typical mean-field approach, here the excess Helmholtz freeenergy functional is represented by a modified fundamental measure theory for the short-range repulsion and by a quadratic density expansion for the long-range attraction. By extensive comparison with simulation results, we find that the DFT can be used for predicting the microscopic structure of liquid droplets, the size, and the free-energy barriers of critical nuclei. We have also calculated the Tolman length of the Lennard- Jones fluids and studied the influence of the droplet curvature on the surface tension. Our calculation suggests that the Tolman length is at least close to zero but the surface tension is still curvature-dependent. We find that the McGrawsLaaksonen assumption is quite reasonable. As the droplet becomes smaller, the surface tension falls gradually from the planar limit. We are able to obtain the free-energy landscape of vapors liquid nucleation directly from numerical iteration. This new method allows us to capture structural and thermodynamic properties of ultrasmall liquid droplets at either equilibrium or nonequilibrium conditions. Acknowledgment We thank Drs. Dong Fu and Vicente A. Talanquer for fruitful discussions. This research is sponsored by the U.S. Department of Energy (DE-FG02-06ER46296) and uses the computational
8 resources from the National Energy Research Scientific Computing Center (NERSC), which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC03-76SF Literature Cited (1) Bustamante, C.; Liphardt, J.; Ritort, F. The nonequilibrium thermodynamics of small systems. Phys. Today 2005, 58, 43. (2) Hill, T. L. Thermodynamics of Small Systems; Dover: New York, (3) Laasonen, K.; Wonczak, S.; Strey, R.; Laaksonen, A. Molecular dynamics simulations of gas-liquid nucleation of Lennard-Jones fluid. J. Chem. Phys. 2000, 113, (4) Neimark, A. V.; Vishnyakov, A. Monte Carlo simulation study of droplet nucleation. J. Chem. Phys. 2005, 122, (5) Neimark, A. V.; Vishnyakov, A. Vapor-to-droplet transition in a Lennard-Jones fluid: Simulation study of nucleation barriers using the ghost field method. J. Phys. Chem. B 2005, 109, (6) Shen, V. K.; Debenedetti, P. G. A computational study of homogeneous liquid-vapor nucleation in the Lennard-Jones fluid. J. Chem. Phys. 1999, 111, (7) ten Wolde, P. R.; Frenkel, D. Computer simulation study of gasliquid nucleation in a Lennard-Jones system. J. Chem. Phys. 1998, 109, (8) Thompson, S. M.; Gubbins, K. E.; Walton, J.; Chantry, R. A. R.; Rowlinson, J. S. A Molecular-Dynamics Study of Liquid-Drops. J. Chem. Phys. 1984, 81, 530. (9) Frenkel, J. Kinetic Theory of Liquids; Dover: New York, (10) Becker, R.; Doring, W. Kinetic treatment of germ formation in supersaturated vapour. Ann. Phys. Berlin 1935, 24, 719. (11) Volmer, M.; Weber, A. Z. Phys.Chem. 1925, 199, 277. (12) Frenkel, J. A general theory of heterophase fluctuations and pretransition phenomena. J. Chem. Phys. 1939, 7, 538. (13) Frenkel, J. Statistical theory of condensation phenomena. J. Chem. Phys. 1939, 7, 200. (14) Drossinos, Y.; Kevrekidis, P. G. Classical nucleation theory revisited. Phys. ReV. E2003, 67, (15) Fu, D. Nucleation for Lennard-Jones fluid by density functional theory. Chin. Phys. Lett. 2005, 22, (16) Koga, K.; Zeng, X. C.; Shchekin, A. K. Validity of Tolman s equation: How large should a droplet be. J. Chem. Phys. 1998, 109, (17) Lee, D. J.; Dagama, M. M. T.; Gubbins, K. E. A Microscopic Theory for Spherical Interfaces - Liquid-Drops in the Canonical Ensemble. J. Chem. Phys. 1986, 85, 490. (18) Oxtoby, D. W.; Evans, R. Nonclassical Nucleation Theory for the Gas-Liquid Transition. J. Chem. Phys. 1988, 89, (19) Talanquer, V.; Oxtoby, D. W. Dynamical Density-Functional Theory of Gas-Liquid Nucleation. J. Chem. Phys. 1994, 100, (20) Talanquer, V.; Oxtoby, D. W. Density-Functional Analysis of Phenomenological Theories of Gas-Liquid Nucleation. J. Phys. Chem. 1995, 99, (21) Zeng, X. C.; Oxtoby, D. W. Gas-Liquid Nucleation in Lennard- Jones Fluids. J. Chem. Phys. 1991, 94, (22) Yu, Y. X.; Wu, J. Z. Structures of hard-sphere fluids from a modified fundamental-measure theory. J. Chem. Phys. 2002, 117, (23) Roth, R.; Evans, R.; Lang, A.; Kahl, G. Fundamental measure theory for hard-sphere mixtures revisited: the White Bear version. J. Phys.: Condens. Matter 2002, 14, Ind. Eng. Chem. Res., Vol. 47, No. 15, (24) Ebner, C.; Saam, W. F.; Stroud, D. Density-Functional Theory of Simple Classical Fluids 0.1. Surfaces. Phys. ReV. A1976, 14, (25) Ebner, C.; Saam, W. F. New Phase-Transition Phenomena in Thin Argon Films. Phys. ReV. Lett. 1977, 38, (26) Tang, Y. P.; Tong, Z. F.; Lu, B. C. Y. Analytical equation of state based on the Ornstein-Zernike equation. Fluid Phase Equilib. 1997, 134, 21. (27) Fu, D.; Wu, J. Z. A self-consistent approach for modelling the interfacial properties and phase diagrams of Yukawa, Lennard-Jones and square-well fluids. Mol. Phys. 2004, 102, (28) Tolman, R. C. The Effect of Droplet Size on Surface Tension. J. Chem. Phys. 1949, 17, 333. (29) Kashchiev, D. Determining the curvature dependence of surface tension. J. Chem. Phys. 2003, 118, (30) McGraw, R.; Laaksonen, A. Interfacial curvature free energy, the Kelvin relation, and vapor-liquid nucleation rate. J. Chem. Phys. 1997, 106, (31) Lei, Y. A.; Bykov, T.; Yoo, S.; Zeng, X. C. The Tolman length: Is it positive or negative. J. Am. Chem. Soc. 2005, 127, (32) Cotterman, R. L.; Schwarz, B. J.; Prausnitz, J. M. Molecular Thermodynamics for Fluids at Low and High-Densities 0.1. Pure Fluids Containing Small or Large Molecules. AIChE J. 1986, 32, (33) Tang, Y. P. Role of the Barker-Henderson diameter in thermodynamics. J. Chem. Phys. 2002, 116, (34) Weeks, J. D.; Chandler, D.; Andersen, H. C. Role of Repulsive Forces in Determining Equilibrium Structure of Simple Liquids. J. Chem. Phys. 1971, 54, (35) Li, Z. D.; Cao, D. P.; Wu, J. Z. Layering, condensation, and evaporation of short chains in narrow slit pores. J. Chem. Phys. 2005, 122, (36) Li, Z. D.; Cao, D. P.; Wu, J. Z. Density-functional theory and Monte Carlo simulation for the surface structure and correlation functions of freely jointed Lennard-Jones polymeric fluids. J. Chem. Phys. 2005, 122, (37) Tang, Y. P.; Wu, J. Z. Modeling inhomogeneous van der Waals fluids using an analytical direct correlation function. Phys. ReV. E2004, 70, (38) Rosenfeld, Y. Free-Energy Model for the Inhomogeneous Hard- Sphere Fluid Mixture and Density-Functional Theory of Freezing. Phys. ReV. Lett. 1989, 63, 980. (39) Fu, D.; Wu, J. Z. Vapor-liquid equilibria and interfacial tensions of associating fluids within a density functional theory. Ind. Eng. Chem. Res. 2005, 44, (40) Cahn, J. W.; Hilliard, J. E. Free Energy of a Nonuniform System 0.3. Nucleation in a 2-Component Incompressible Fluid. J. Chem. Phys. 1959, 31, 688. (41) Wood, S. M.; Wang, Z. G. Nucleation in binary polymer blends: A self-consistent field study. J. Chem. Phys. 2002, 116, (42) Trokhymchuk, A.; Alejandre, J. Computer simulations of liquid/ vapor interface in Lennard-Jones fluids: Some questions and answers. J. Chem. Phys. 1999, 111, (43) Blokhuis, E. M.; Bedeaux, D. Pressure Tensor of a Spherical Interface. J. Chem. Phys. 1992, 97, (44) Smit, B. Phase-Diagrams of Lennard-Jones Fluids. J. Chem. Phys. 1992, 96, ReceiVed for review April 25, 2007 ReVised manuscript received June 7, 2007 Accepted June 8, 2007 IE070578I
of Nebraska - Lincoln
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Xiao Cheng Zeng Publications Published Research - Department of Chemistry 10-1-2006 Homogeneous nucleation at high supersaturation
More informationThermodynamic expansion of nucleation free-energy barrier and size of critical nucleus near the vapor-liquid coexistence
JOURNAL OF CHEMICAL PHYSICS VOLUME 110, NUMBER 7 15 FEBRUARY 1999 hermodynamic expansion of nucleation free-energy barrier and size of critical nucleus near the vapor-liquid coexistence Kenichiro Koga
More informationA patching model for surface tension of spherical droplet and Tolman length. II
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Xiao Cheng Zeng Publications Published Research - Department of Chemistry -5-999 A patching model for surface tension of
More informationCrystal nucleation for a model of globular proteins
JOURNAL OF CHEMICAL PHYSICS VOLUME 120, NUMBER 17 1 MAY 2004 Crystal nucleation for a model of globular proteins Andrey Shiryayev and James D. Gunton Department of Physics, Lehigh University, Bethlehem,
More informationSolvation of a Spherical Cavity in Simple Liquids: Interpolating between the Limits
Article Subscriber access provided by - Access paid by the UC Riverside Libraries Solvation of a Spherical Cavity in Simple Liquids: Interpolating between the Limits Jianzhong Wu J. Phys. Chem. B, 2009,
More informationDirect determination of the Tolman length from the bulk pressures of liquid drops via molecular dynamics simulations
Direct determination of the Tolman length from the bulk pressures of liquid drops via molecular dynamics simulations Alan E. van Giessen, and Edgar M. Blokhuis Citation: The Journal of Chemical Physics
More informationHeterogeneous nucleation on mesoscopic wettable particles: A hybrid thermodynamic/densityfunctional
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Xiao Cheng Zeng Publications Published esearch - Department of Chemistry 7--00 Heterogeneous nucleation on mesoscopic wettable
More informationBinary Hard-Sphere Mixtures Within Spherical Pores
Journal of the Korean Physical Society, Vol. 35, No. 4, October 1999, pp. 350 354 Binary Hard-Sphere Mixtures Within Spherical Pores Soon-Chul Kim Department of Physics, Andong National University, Andong
More informationA density-functional theory for bulk and inhomogeneous Lennard-Jones fluids from the energy route
JOURNAL OF CHEMICAL PHYSICS VOLUME 119, NUMBER 14 8 OCTOBER 2003 A density-functional theory for bulk and inhomogeneous Lennard-Jones fluids from the energy route Yiping Tang a) Honeywell Hi-Spec Solutions,
More informationA patching model for surface tension and the Tolman length
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Xiao Cheng Zeng Publications Published Research - Department of Chemistry 8-22-1999 A patching model for surface tension
More informationVapor-to-Droplet Transition in a Lennard-Jones Fluid: Simulation Study of Nucleation Barriers Using the Ghost Field Method
5962 J. Phys. Chem. B 2005, 109, 5962-5976 Vapor-to-Droplet Transition in a Lennard-Jones Fluid: Simulation Study of Nucleation Barriers Using the Ghost Field Method Alexander V. Neimark* and Aleksey Vishnyakov
More informationA computational study of homogeneous liquid vapor nucleation in the Lennard-Jones fluid
JOURNAL OF CHEMICAL PHYSICS VOLUME 111, NUMBER 8 22 AUGUST 1999 A computational study of homogeneous liquid vapor nucleation in the Lennard-Jones fluid Vincent K. Shen and Pablo G. Debenedetti a) Department
More informationhydrate systems Gránásy Research Institute for Solid State Physics & Optics H-1525 Budapest, POB 49, Hungary László
Phase field theory of crystal nucleation: Application to the hard-sphere and CO 2 Bjørn Kvamme Department of Physics, University of Bergen Allégaten 55, N-5007 N Bergen, Norway László Gránásy Research
More informationMonte Carlo Calculations of Effective Surface Tension for Small Clusters
Monte Carlo Calculations of Effective Surface Tension for Small Clusters Barbara N. Hale Physics Department and Cloud and Aerosol Science Laboratory, University of Missouri- Rolla, Rolla, MO 65401, USA
More informationNucleation rate (m -3 s -1 ) Radius of water nano droplet (Å) 1e+00 1e-64 1e-128 1e-192 1e-256
Supplementary Figures Nucleation rate (m -3 s -1 ) 1e+00 1e-64 1e-128 1e-192 1e-256 Calculated R in bulk water Calculated R in droplet Modified CNT 20 30 40 50 60 70 Radius of water nano droplet (Å) Supplementary
More informationA kinetic theory of homogeneous bubble nucleation
JOURNAL OF CHEMICAL PHYSICS VOLUME 118, NUMBER 2 8 JANUARY 2003 A kinetic theory of homogeneous bubble nucleation Vincent K. Shen and Pablo G. Debenedetti a) Department of Chemical Engineering, Princeton
More informationYang-Xin Yu* and Feng-Qi You. Yiping Tang. Guang-Hua Gao and Yi-Gui Li
334 J. Phys. Chem. B 2006, 110, 334-341 Structure and Adsorption of A Hard-Core Multi-Yukawa Fluid Confined in A Slitlike Pore: Grand Canonical Monte Carlo Simulation and Density Functional Study Yang-Xin
More informationCluster formation and bulk phase behavior of colloidal dispersions
Cluster formation and bulk phase behavior of colloidal dispersions Tao Jiang and Jianzhong Wu* Department of Chemical and Environmental Engineering, University of California, Riverside, California 92521-0444,
More informationChemical Potential. Combining the First and Second Laws for a closed system, Considering (extensive properties)
Chemical Potential Combining the First and Second Laws for a closed system, Considering (extensive properties) du = TdS pdv Hence For an open system, that is, one that can gain or lose mass, U will also
More informationMelting line of the Lennard-Jones system, infinite size, and full potential
THE JOURNAL OF CHEMICAL PHYSICS 127, 104504 2007 Melting line of the Lennard-Jones system, infinite size, and full potential Ethan A. Mastny a and Juan J. de Pablo b Chemical and Biological Engineering
More informationSurface property corrected modification of the classical nucleation theory
CCP5 Annual Meeting Surface property corrected modification of the classical nucleation theory Sheffield Hallam University, September 15, 2010 Martin Horsch, Hans Hasse, and Jadran Vrabec The critical
More informationSurface Tension of the Vapor Liquid Interface with Finite Curvature
Colloid Journal, Vol. 65, No. 4, 00, pp. 44045. Translated from Kolloidnyi Zhurnal, Vol. 65, No. 4, 00, pp. 480494. Original Russian Text Copyright 00 by Zhukhovitskii. Surface Tension of the VaporLiquid
More informationWeighted density functional theory of the solvophobic effect
PHYSICAL REVIEW E, VOLUME 64, 021512 Weighted density functional theory of the solvophobic effect Sean X. Sun Department of Chemistry, University of California, Berkeley, California 94720 Received 24 July
More informationDensity functional theory of inhomogeneous liquids. III. Liquid-vapor nucleation
THE JOURNAL OF CHEMICAL PHYSICS 129, 24451 28 Density functional theory of inhomogeneous liquids. III. Liquid-vapor nucleation James F. Lutsko a Center for Nonlinear Phenomena and Complex Systems CP 231,
More informationSurface analysis algorithms in the mardyn program and the ls1 project
Surface analysis algorithms in the mardyn program and the ls1 project Stuttgart, 15 th December 1 M. T. Horsch Surface tension The virial route Bakker-Buff equation: γ R 2 out in dz z Normal pressure decays
More informationDensity-functional theory for the structures and thermodynamic properties of highly asymmetric electrolyte and neutral component mixtures
PHYSICAL REVIEW E 70, 031109 (2004) Density-functional theory for the structures and thermodynamic properties of highly asymmetric electrolyte and neutral component mixtures Zhidong Li and Jianzhong Wu*
More informationDependence of the surface tension on curvature
Dependence of the surface tension on curvature rigorously determined from the density profiles of nanodroplets Athens, st September M. T. Horsch,,, S. Eckelsbach, H. Hasse, G. Jackson, E. A. Müller, G.
More informationVIII. Phase Transformations. Lecture 38: Nucleation and Spinodal Decomposition
VIII. Phase Transformations Lecture 38: Nucleation and Spinodal Decomposition MIT Student In this lecture we will study the onset of phase transformation for phases that differ only in their equilibrium
More informationA New Uniform Phase Bridge Functional: Test and Its Application to Non-uniform Phase Fluid
Commun. Theor. Phys. (Beijing, China) 39 (2003) pp. 231 237 c International Academic Publishers Vol. 39, No. 2, February 15, 2003 A New Uniform Phase Bridge Functional: Test and Its Application to Non-uniform
More informationGibbs ensemble simulation of phase equilibrium in the hard core two-yukawa fluid model for the Lennard-Jones fluid
MOLECULAR PHYSICS, 1989, VOL. 68, No. 3, 629-635 Gibbs ensemble simulation of phase equilibrium in the hard core two-yukawa fluid model for the Lennard-Jones fluid by E. N. RUDISILL and P. T. CUMMINGS
More informationVapor-Liquid Equilibria and Interfacial Tensions of Associating Fluids within a Density Functional Theory
1120 Ind. Eng. Chem. Res. 2005, 44, 1120-1128 Vapor-Liquid Equilibria and Interfacial Tensions of Associating Fluids within a Density Functional Theory Dong Fu and Jianzhong Wu*, Department of Environmental
More informationComparative study on methodology in molecular dynamics simulation of nucleation
THE JOURNAL OF CHEMICAL PHYSICS 126, 224517 2007 Comparative study on methodology in molecular dynamics simulation of nucleation Jan Julin, Ismo Napari, and Hanna Vehkamäki Department of Physical Sciences,
More informationIntroduction Statistical Thermodynamics. Monday, January 6, 14
Introduction Statistical Thermodynamics 1 Molecular Simulations Molecular dynamics: solve equations of motion Monte Carlo: importance sampling r 1 r 2 r n MD MC r 1 r 2 2 r n 2 3 3 4 4 Questions How can
More informationThermodynamics of Three-phase Equilibrium in Lennard Jones System with a Simplified Equation of State
23 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 28 (2014) Thermodynamics of Three-phase Equilibrium in Lennard Jones System with a Simplified Equation of State Yosuke
More informationHard-sphere fluids in contact with curved substrates
Hard-sphere fluids in contact with curved substrates P. Bryk, 1,2,3 R. Roth, 2,3 K. R. Mecke, 2,3 and S. Dietrich 2,3 1 Department for the Modeling of Physico-Chemical Processes, Maria Curie-Skłodowska
More informationDensity Functional Theory Study on the Structure and Capillary Phase Transition of a Polymer Melt in a Slitlike Pore: Effect of Attraction
14418 J. Phys. Chem. B 2006, 110, 14418-14425 Density Functional Theory Study on the Structure and Capillary Phase Transition of a Polymer elt in a Slitlike Pore: Effect of Attraction Yang-Xin Yu,* Guang-Hua
More informationDENSITY FUNCTIONAL THEORY FOR STUDIES OF MULTIPLE STATES OF INHOMOGENEOUS FLUIDS AT SOLID SURFACES AND IN PORES.
J. Smith, D. Stroud, MRS Symposium Proceedings Series, v.49, p.7-33, 998. DENSITY FUNCTIONAL THEORY FOR STUDIES OF MULTIPLE STATES OF INHOMOGENEOUS FLUIDS AT SOLID SURFACES AND IN PORES. A.. NEIMARK, and
More informationPhase transitions of quadrupolar fluids
Phase transitions of quadrupolar fluids Seamus F. O Shea Department of Chemistry, University of Lethbridge, Lethbridge, Alberta, Canada, T1K 3M4 Girija S. Dubey Brookhaven National Laboratory, Upton, New
More informationMOLECULAR DYNAMICS SIMULATION OF HETEROGENEOUS NUCLEATION OF LIQUID DROPLET ON SOLID SURFACE
MOLECULAR DYNAMICS SIMULATION OF HETEROGENEOUS NUCLEATION OF LIQUID DROPLET ON SOLID SURFACE Tatsuto Kimura* and Shigeo Maruyama** *Department of Mechanical Engineering, The University of Tokyo 7-- Hongo,
More informationAn Extended van der Waals Equation of State Based on Molecular Dynamics Simulation
J. Comput. Chem. Jpn., Vol. 8, o. 3, pp. 97 14 (9) c 9 Society of Computer Chemistry, Japan An Extended van der Waals Equation of State Based on Molecular Dynamics Simulation Yosuke KATAOKA* and Yuri YAMADA
More informationDensity functional theory for molecular orientation of hard rod fluids in hard slits
Vol 16 No 8, August 2007 c 2007 Chin. Phys. Soc. 1009-1963/2007/16(08)/2296-08 Chinese Physics and IOP Publishing Ltd Density functional theory for molecular orientation of hard rod fluids in hard slits
More informationDensity Functional Theory of the Interface between Solid and Superfluid Helium 4
Density Functional Theory of the Interface between Solid and Superfluid Helium 4 Frédéric Caupin and Tomoki Minoguchi Laboratoire de Physique Statistique de l Ecole Normale Supérieure associé aux Universités
More informationUsing Fundamental Measure Theory to Treat the Correlation Function of the Inhomogeneous Hard-Sphere Fluid
Using Fundamental Measure Theory to Treat the Correlation Function of the Inhomogeneous Hard-Sphere Fluid Jeff B. Schulte, Patrick A. Kreitzberg, Chris V. Haglund, and David Roundy Department of Physics,
More informationHill s nano-thermodynamics is equivalent with Gibbs thermodynamics for curved surfaces
Hill s nano-thermodynamics is equivalent with Gibbs thermodynamics for curved surfaces arxiv:1806.03444v1 [physics.chem-ph] 9 Jun 2018 Dick Bedeaux and Signe Kjelstrup Porelab, Department of Chemistry,
More informationSelf-consistent field theory for the nucleation of micelles
JOURNAL OF CHEMICAL PHYSICS VOLUME 110, NUMBER 11 15 MARCH 1999 Self-consistent field theory for the nucleation of micelles N. A. M. Besseling a) and M. A. Cohen Stuart Laboratory for Physical Chemistry
More informationContents. 1 Introduction and guide for this text 1. 2 Equilibrium and entropy 6. 3 Energy and how the microscopic world works 21
Preface Reference tables Table A Counting and combinatorics formulae Table B Useful integrals, expansions, and approximations Table C Extensive thermodynamic potentials Table D Intensive per-particle thermodynamic
More informationTitle Super- and subcritical hydration of Thermodynamics of hydration Author(s) Matubayasi, N; Nakahara, M Citation JOURNAL OF CHEMICAL PHYSICS (2000), 8109 Issue Date 2000-05-08 URL http://hdl.handle.net/2433/50350
More informationGas-liquid phase separation in oppositely charged colloids: stability and interfacial tension
7 Gas-liquid phase separation in oppositely charged colloids: stability and interfacial tension We study the phase behaviour and the interfacial tension of the screened Coulomb (Yukawa) restricted primitive
More informationDensity functional theory for inhomogeneous mixtures of polymeric fluids
JOURNAL OF CHEICAL PHYSICS VOLUE 117, NUBER 5 1 AUGUST 2002 Density functional theory for inhomogeneous mixtures of polymeric fluids Yang-Xin Yu a) and Jianzhong Wu b) Department of Chemical and Environmental
More informationImproved association in a classical density functional theory for water. Abstract
Improved association in a classical density functional theory for water Eric J. Krebs, Jeff B. Schulte, and David Roundy Department of Physics, Oregon State University, Corvallis, OR 97331, USA Abstract
More informationJacco Snoeijer PHYSICS OF FLUIDS
Jacco Snoeijer PHYSICS OF FLUIDS dynamics dynamics freezing dynamics freezing microscopics of capillarity Menu 1. surface tension: thermodynamics & microscopics 2. wetting (statics): thermodynamics & microscopics
More informationComparison of different mixing rules for prediction of density and residual internal energy of binary and ternary Lennard Jones mixtures
Fluid Phase Equilibria 178 (2001) 87 95 Comparison of different mixing rules for prediction of density and residual internal energy of binary and ternary Lennard Jones mixtures Jian Chen a,, Jian-Guo Mi
More informationSurface thermodynamics of planar, cylindrical, and spherical vapour-liquid interfaces of water
Surface thermodynamics of planar, cylindrical, and spherical vapour-liquid interfaces of water Gabriel V. Lau, Ian J. Ford, Patricia A. Hunt, Erich A. Müller, and George Jackson Citation: The Journal of
More informationCARBON 2004 Providence, Rhode Island. Adsorption of Flexible n-butane and n-hexane on Graphitized Thermal Carbon Black and in Slit Pores
CARBON Providence, Rhode Island Adsorption of Flexible n-butane and n-hexane on Graphitized Thermal Carbon Black and in Slit Pores D. D. Do* and H. D. Do, University of Queensland, St. Lucia, Qld 7, Australia
More informationLOCAL PRESSURE OF CONFINED FLUIDS INSIDE NANOSLIT PORES (A DENSITY FUNCTIONAL THEORY PREDICTION)
International Journal of Theoretical Physics, Group Theory, and Nonlinear Optics Volume 17, Numer, pp.101-119 (01) ISSN 154 4674 LOCAL PRESSURE OF CONFINED FLUIDS INSIDE NANOSLIT PORES Fatemeh Heidari
More informationMolecular dynamics simulation of atomic clusters in equilibrium with a vapour
Molecular Simulation, Vol. 33, No. 3, 15 March 2007, 245 251 Molecular dynamics simulation of atomic clusters in equilibrium with a vapour M. SALONEN*, I. NAPARI and H. VEHKAMÄKI Department of Physical
More informationPORE SIZE DISTRIBUTION OF CARBON WITH DIFFERENT PROBE MOLECULES
PORE SIZE DISTRIBUTION OF CARBON WITH DIFFERENT PROBE MOLECULES Atichat Wongkoblap*, Worapot Intomya, Warangkhana Somrup, Sorod Charoensuk, Supunnee Junpirom and Chaiyot Tangsathitkulchai School of Chemical
More informationMolecular dynamics simulation of the liquid vapor interface: The Lennard-Jones fluid
Molecular dynamics simulation of the liquid vapor interface: The Lennard-Jones fluid Matthias Mecke and Jochen Winkelmann a) Institut für Physikalische Chemie, Universität Halle-Wittenberg, Geusaer Str.,
More informationMolecular Driving Forces
Molecular Driving Forces Statistical Thermodynamics in Chemistry and Biology SUBGfittingen 7 At 216 513 073 / / Ken A. Dill Sarina Bromberg With the assistance of Dirk Stigter on the Electrostatics chapters
More informationsubmitted to Physical Review E Numerical Tests of Nucleation Theories for the Ising Models Seunghwa Ryu 1 and Wei Cai 2
submitted to Physical Review E Numerical Tests of Nucleation Theories for the Ising Models Seunghwa Ryu 1 and Wei Cai 2 1 Department of Physics, Stanford University, Stanford, California 94305 2 Department
More informationStatistical Mechanics of a Thin Film on a Solid Substrate
Statistical Mechanics of a Thin Film on a Solid Substrate arxiv:1406.0698v2 [cond-mat.stat-mech] 6 Jan 2017 Diploma Thesis, revised version Andreas Nold submitted at the Technische Universität Darmstadt
More informationNucleation: Theory and Applications to Protein Solutions and Colloidal Suspensions
Nucleation: Theory and Applications to Protein Solutions and Colloidal Suspensions Richard P. Sear Department of Physics, University of Surrey Guildford, Surrey GU2 7XH, United Kingdom r.sear@surrey.ac.uk
More informationJOURNAL REVIEW. Density Functional Theory for Chemical Engineering: From Capillarity to Soft Materials
THERMODYNAMICS JOURNAL REVIEW Density Functional Theory for Chemical Engineering: From Capillarity to Soft Materials Jianzhong Wu Dept. of Chemical and Environmental Engineering, University of California,
More informationPhase Equilibria of binary mixtures by Molecular Simulation and PR-EOS: Methane + Xenon and Xenon + Ethane
International Journal of ChemTech Research CODEN( USA): IJCRGG ISSN : 0974-4290 Vol.5, No.6, pp 2975-2979, Oct-Dec 2013 Phase Equilibria of binary mixtures by Molecular Simulation and PR-EOS: Methane +
More informationNUCLEATION THEORY USING EQUATIONS OF STATE ABDALLA A. OBEIDAT ATHESIS. Presented to the Faculty of the Graduate School of the
NUCLEATION THEORY USING EQUATIONS OF STATE by ABDALLA A. OBEIDAT ATHESIS Presented to the Faculty of the Graduate School of the UNIVERSITY OF MISSOURI-ROLLA in Partial Fulfillment of the Requirements for
More informationAdsorption properties of a colloid-polymer mixture confined in a slit pore
PHYSICAL REVIEW E, VOLUME 64, 041507 Adsorption properties of a colloid-polymer mixture confined in a slit pore Soon-Chul Kim* Department of Physics, Andong National University, Andong, 760-749 Korea Peter
More informationEQUATION OF STATE DEVELOPMENT
EQUATION OF STATE DEVELOPMENT I. Nieuwoudt* & M du Rand Institute for Thermal Separation Technology, Department of Chemical Engineering, University of Stellenbosch, Private bag X1, Matieland, 760, South
More informationCH.7 Fugacities in Liquid Mixtures: Models and Theories of Solutions
CH.7 Fugacities in Liquid Mixtures: Models and Theories of Solutions The aim of solution theory is to express the properties of liquid mixture in terms of intermolecular forces and liquid structure. The
More informationChapter 7. Pickering Stabilisation ABSTRACT
Chapter 7 Pickering Stabilisation ABSTRACT In this chapter we investigate the interfacial properties of Pickering emulsions. Based upon findings that indicate these emulsions to be thermodynamically stable,
More informationLennard-Jones as a model for argon and test of extended renormalization group calculations
JOURNAL OF CHEMICAL PHYSICS VOLUME 111, NUMBER 2 22 NOVEMBER 1999 Lennard-Jones as a model for argon and test of extended renormalization group calculations John A. White Department of Physics, American
More informationStatistical mechanics of nucleation: Incorporating translational and rotational free energy into thermodynamics of a microdroplet
PHYSICAL REVIEW E 7, 0607 2006 Statistical mechanics of nucleation: Incorporating translational and rotational free energy into thermodynamics of a microdroplet Isamu Kusaka The Koffolt Laboratories, The
More informationfiziks Institute for NET/JRF, GATE, IIT-JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES
Content-Thermodynamics & Statistical Mechanics 1. Kinetic theory of gases..(1-13) 1.1 Basic assumption of kinetic theory 1.1.1 Pressure exerted by a gas 1.2 Gas Law for Ideal gases: 1.2.1 Boyle s Law 1.2.2
More information510 Subject Index. Hamiltonian 33, 86, 88, 89 Hamilton operator 34, 164, 166
Subject Index Ab-initio calculation 24, 122, 161. 165 Acentric factor 279, 338 Activity absolute 258, 295 coefficient 7 definition 7 Atom 23 Atomic units 93 Avogadro number 5, 92 Axilrod-Teller-forces
More informationEstimation of Saturated Vapor Pressure from Nucleation Data
CHINESE JOURNAL OF CHEMICAL PHYSICS VOLUME 19, NUMBER 5 OCTOBER 27, 2006 Estimation of Saturated Vapor Pressure from Nucleation Data Mei-li Wang a, Ji-ming Yang a, Xi-sheng Luo b a. Department of Modern
More informationSupplemental Material for Temperature-sensitive colloidal phase behavior induced by critical Casimir forces
Supplemental Material for Temperature-sensitive colloidal phase behavior induced by critical Casimir forces Minh Triet Dang, 1 Ana Vila Verde, 2 Van Duc Nguyen, 1 Peter G. Bolhuis, 3 and Peter Schall 1
More informationEntropic wetting and the fluid fluid interface of a model colloid polymer mixture
INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS: CONDENSED MATTER J. Phys.: Condens. Matter 4 (22) L L8 PII: S953-8984(2)28952-3 LETTER TO THE EDITOR Entropic wetting and the fluid fluid interface of
More informationINHOMOGENEOUS DENSITY THEORIES COUPLED INTO A MOLECULAR EQUATION OF STATE FOR THE DESCRIPTION OF THE FLUID-FLUID INTERFACE
12th Joint European Thermodynamics Conference Brescia, July 1-5, 2013 INHOMOGENEOUS DENSITY THEORIES COUPLED INTO A MOLECULAR EQUATION OF STATE FOR THE DESCRIPTION OF THE FLUID-FLUID INTERFACE Felix Llovell*,
More information2. Thermodynamics. Introduction. Understanding Molecular Simulation
2. Thermodynamics Introduction Molecular Simulations Molecular dynamics: solve equations of motion r 1 r 2 r n Monte Carlo: importance sampling r 1 r 2 r n How do we know our simulation is correct? Molecular
More informationSupplemental Material for Curvature-dependence of the liquid-vapor surface tension beyond the Tolman approximation
Supplemental Material for Curvature-dependence of the liquid-vapor surface tension beyond the Tolman approximation Nicolas Bruot and Frédéric Caupin Institut Lumière Matière, UMR536 Université Claude Bernard
More informationWetting Transitions at Fluid Interfaces and Related Topics
Wetting Transitions at Fluid Interfaces and Related Topics Kenichiro Koga Department of Chemistry, Faculty of Science, Okayama University Tsushima-Naka 3-1-1, Okayama 7-853, Japan Received April 3, 21
More informationInterfaces and interfacial energy
Interfaces and interfacial energy 1/14 kinds: l/g }{{ l/l } mobile s/g s/l s/s Example. Estimate the percetage of water molecules on the surface of a fog droplet of diameter (i) 0.1 mm (naked eye visibility
More informationOn the local and nonlocal components of solvation thermodynamics and their relation to solvation shell models
JOURNAL OF CHEMICAL PHYSICS VOLUME 109, NUMBER 12 22 SEPTEMBER 1998 On the local and nonlocal components of solvation thermodynamics and their relation to solvation shell models Nobuyuki Matubayasi Institute
More informationNUCLEATION IN FINITE SYSTEMS: THEORY AND COMPUTER SIMULATION*t. M. RAO and B. J. BERNE. Chemistry Dept, Columbia University, New York, U.S.A.
NUCLEATION IN FINITE SYSTEMS: THEORY AND COMPUTER SIMULATION*t M. RAO and B. J. BERNE Chemistry Dept, Columbia University, New York, U.S.A. Abstract. The free energy of formation of a droplet in a finite
More informationDirect calculation of interfacial tensions from computer simulation: Results for freely jointed tangent hard sphere chains
Direct calculation of interfacial tensions from computer simulation: Results for freely jointed tangent hard sphere chains Luis G. MacDowell Departamento de Química Física, Facultad de Ciencias Químicas,
More informationPhase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany
Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals Phase Equilibria Phase diagrams and classical thermodynamics
More informationMolecular Modeling of Matter
Molecular Modeling of Matter Keith E. Gubbins Lecture 1: Introduction to Statistical Mechanics and Molecular Simulation Common Assumptions Can treat kinetic energy of molecular motion and potential energy
More informationMonte Carlo simulations of confined hard sphere fluids
Physics Procedia 3 (2010) 1493 1497 www.elsevier.com/locate/procedia Monte Carlo simulations of confined hard sphere fluids G.E. Schröder-Turk, S. Kuczera, D. Hörndlein, F. Krach, T. Hiester, K. Mecke
More informationUB association bias algorithm applied to the simulation of hydrogen fluoride
Fluid Phase Equilibria 194 197 (2002) 249 256 UB association bias algorithm applied to the simulation of hydrogen fluoride Scott Wierzchowski, David A. Kofke Department of Chemical Engineering, University
More informationTheoretical Studies of the Correlations in Moderately Asymmetric Binary Hard-Sphere Solid Mixtures
Ames Laboratory Publications Ames Laboratory 5-008 Theoretical Studies of the Correlations in Moderately Asymmetric Binary Hard-Sphere Solid Mixtures Vadim B. Warshavsky Iowa State University Xueyu Song
More informationPrinciples of Equilibrium Statistical Mechanics
Debashish Chowdhury, Dietrich Stauffer Principles of Equilibrium Statistical Mechanics WILEY-VCH Weinheim New York Chichester Brisbane Singapore Toronto Table of Contents Part I: THERMOSTATICS 1 1 BASIC
More informationDensity dependence of dielectric permittivity of water and estimation of the electric field for the breakdown inception
Journal of Physics: Conference Series PAPER OPEN ACCESS Density dependence of dielectric permittivity of water and estimation of the electric field for the breakdown inception To cite this article: D I
More informationPerturbation theory calculations of model pair potential systems
Graduate Theses and Dissertations Graduate College 2016 Perturbation theory calculations of model pair potential systems Jianwu Gong Iowa State University Follow this and additional works at: http://lib.dr.iastate.edu/etd
More informationCH 240 Chemical Engineering Thermodynamics Spring 2007
CH 240 Chemical Engineering Thermodynamics Spring 2007 Instructor: Nitash P. Balsara, nbalsara@berkeley.edu Graduate Assistant: Paul Albertus, albertus@berkeley.edu Course Description Covers classical
More informationStep 1. Step 2. g l = g v. dg = 0 We have shown that over a plane surface of water. g v g l = ρ v R v T ln e/e sat. this can be rewritten
The basic question is what makes the existence of a droplet thermodynamically preferable to the existence only of water vapor. We have already derived an expression for the saturation vapor pressure over
More informationPressure Dependent Study of the Solid-Solid Phase Change in 38-Atom Lennard-Jones Cluster
University of Rhode Island DigitalCommons@URI Chemistry Faculty Publications Chemistry 2005 Pressure Dependent Study of the Solid-Solid Phase Change in 38-Atom Lennard-Jones Cluster Dubravko Sabo University
More informationEnergetics of small n-pentanol clusters from droplet nucleation rate data
JOURNAL OF CHEMICAL PHYSICS VOLUME 112, NUMBER 12 22 MARCH 2000 Energetics of small n-pentanol clusters from droplet nucleation rate data Michael Knott, Hanna Vehkamäki, a) and Ian J. Ford Department of
More informationPre-yield non-affine fluctuations and a hidden critical point in strained crystals
Supplementary Information for: Pre-yield non-affine fluctuations and a hidden critical point in strained crystals Tamoghna Das, a,b Saswati Ganguly, b Surajit Sengupta c and Madan Rao d a Collective Interactions
More informationModeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase. Andy Long Kridsanaphong Limtragool
Modeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase Andy Long Kridsanaphong Limtragool Motivation We want to study the spontaneous formation of micelles and vesicles Applications
More informationStatistical geometry of cavities in a metastable confined fluid
PHYSICAL REVIEW E VOLUME 62, NUMBER 1 JULY 2000 Statistical geometry of cavities in a metastable confined fluid Aleksey Vishnyakov, 1 Pablo G. Debenedetti, 2 and Alexander V. Neimark 1, * 1 TRI/Princeton,
More informationGases and the Virial Expansion
Gases and the irial Expansion February 7, 3 First task is to examine what ensemble theory tells us about simple systems via the thermodynamic connection Calculate thermodynamic quantities: average energy,
More information