PAPER First- and second-order wetting transitions at liquid vapor interfaces K. Koga,* a J. O. Indekeu* b and B. Widom* c Received 8th December 2009, Accepted 13th January 2010 DOI: /b925671g Wetting

•0 views•6 pages

Previous

Next

All materials on our website are shared by users. If you have any questions about copyright issues, please report us to resolve them. We are always happy to assist you.

Documenttranscript

PAPER First- and second-order wetting transitions at liquid vapor interfaces K. Koga,* a J. O. Indekeu* b and B. Widom* c Received 8th December 2009, Accepted 13th January 2010 DOI: /b925671g Wetting transitions, in which one liquid wets, or spreads at, the interface between a second liquid and their common vapor, are defined and first- and second-order transitions are distinguished. The mean-field density-functional models of fluid interfaces are recalled. A criterion is noted for determining when the wetting transitions in those models are required to be of first order or may be of second order. It is seen how two examples of such density-functional models that have been treated in the past, one leading to a first-order and the other to a secondorder wetting transition, provide examples of the application of the criterion. 1. Introduction Faraday Discussions We consider two liquid phases (called phases 2 and 3) in equilibrium with their common vapor (phase 1). The three phases may meet at a line of common contact with three non-vanishing contact angles, or alternatively one of the phases, say 2, may wet (spread at) the interface between 1 and 3. Let the three interfacial tensions be denoted s 12, s 23, and s 13, and suppose the liquid vapor tension s 13 to be the greatest of the three. Then the condition for the three phases to meet at a line of common contact ( non-wetting ) is the triangle inequality s 13 s 12 + s 23 (non-wetting). (1) The alternative ( wetting ) is that phase 2 spreads as a macroscopically thick film at the 1,3 interface, with s 13 ¼ s 12 + s 23 (wetting). (2) A transition between these two modes of three-phase contact is termed a wetting transition. Such a transition might occur as the temperature or a chemical potential varies, leading to two continuous ranges of three-phase states (triple points), with (1) holding over one range and (2) over the other. Between those two ranges of triple points there would then be a boundary at which (1) goes over continuously to (2) or the reverse. According to the phase rule, to have a continuum of three-phase states instead of just isolated triple points requires that the system be a mixture of at least two components. In a one-component system there are only isolated triple a Department of Chemistry, Faculty of Science, Okayama University, Okayama, , Japan. b Institute for Theoretical Physics, Katholieke Universiteit Leuven, BE-3001 Leuven, Belgium. c Department of Chemistry, Baker Laboratory, Cornell University, Ithaca, NY, , USA. This journal is ª The Royal Society of Chemistry 2010 Faraday Discuss., 2010, 146, points, where either one of the phases wets the interface between the other two, or not. There is no transition between those two circumstances because no field variable can then vary without one or two of the three phases disappearing. For simplicity we consider here systems with the minimum number, two, of components. By the phase rule, the system of three phases and two components then has one degree of freedom; i.e., there is then a one-dimensional manifold of triple points parameterized by a single field variable, which we henceforth call b, and which we may think of as some function of the temperature and the two chemical potentials. The wetting transition between (1) and (2) occurs when b takes some value b w. We shall take b b w to be the range of non-wet states, where (1) is then satisfied, and b b w to be the range of states where the 1,3 interface is wet by 2, so that (2) holds. In reality, many such liquid liquid vapor three-phase equilibria occur in systems of c 2 components. To model such systems and still have the simplicity of having only one variable field we would imagine holding fixed c 2 of the system s c +1 field variables. Then with the two phase-equilibrium constraints we would still be considering only a one-dimensional manifold of triple points, parameterized by a single varying field. The distinction between first- and second-order wetting transitions is illustrated in Fig. 1, where s 12 + s 23 s 13, which we call Ds, is shown schematically as a function of b. In the first-order transition Ds vanishes linearly with b b w, and has a metastable extension, shown by the dashed line in Fig. 1(a), in which the non-wet interface for b b w is at a local but not global minimum of the free energy. At the global minimum, Ds ¼ 0 over the range b b w. At a classical second-order wetting transition (in the sense of the Ehrenfest classification), Ds vanishes proportionally to (b b w ) 2 as b approaches b w from above [Fig. 1(b)], and does not have a metastable extension to states with b b w. A wetting transition may be very weakly first order, when the curve of Ds in Fig. 1(a) crosses the b-axis at b ¼ b w so shallowly that it may be difficult to locate b w or to distinguish the slope of the curve at b w from Mean-field density-functional theories recalled In each of the three interfaces, the densities r 1 and r 2 of the two components vary in the direction z perpendicular to that interface. We define an excess free energy of inhomogeneity J[r 1 (z), r 2 (z);b], which is a functional of the two spatially varying densities and depends on the single field variable b. We take it to be of the classical square-gradient form J ¼ F½r 1 ðzþ; r 2 ðzþ; bš þ 1 h r ðzþ2 þ r 0 2 ðzþ2i (3) Fig. 1 (a) First-order transition; (b) second-order transition. 218 Faraday Discuss., 2010, 146, This journal is ª The Royal Society of Chemistry 2010 where F(r 1, r 2 ;b) is such that it, along with vf/vr 1 and vf/vr 2, vanish when, for the given b, the densities r 1 and r 2 are those of any of the three coexisting bulk phases; while F 0 otherwise. The tension s of each of the three interfaces is then s ¼ min r 1 ðzþ; r 2 ðzþ ð N N J dz (4) i.e., the variational minimum of the integrated free-energy density, minimized with respect to the densities r 1 (z), r 2 (z) subject to their taking their respective bulk-phase values at z ¼N. The minimizing r 1 (z) and r 2 (z) satisfy the Euler Lagrange equations vf ¼ d2 r 1 vr 1 dz ; vf ¼ d2 r 2 (5) 2 vr 2 dz 2 These are to be solved subject to the boundary conditions at z ¼N noted above. When z is eliminated between r 1 (z) and r 2 (z) for any of the interfaces the result is a trajectory in the r 1, r 2 plane that describes how r 1 and r 2 vary with each other through that interface. There are three such trajectories, one for each of the three interfaces, as shown schematically in Fig. 2. The density variables are taken to have the values r 1 ¼ 1, r 2 ¼ 0 in phase 1; r 1 ¼ 0, r 2 ¼ b in phase 2; and r 1 ¼ 1, r 2 ¼ 0 in phase 3. Here r 1 and r 2 are not to be taken literally as densities but rather as two independent functions of the physical densities, chosen to make the notation and the representation in Fig. 2 simple. The curve in Fig. 2 extending from 1,0 to 0, b is the 1,2 interfacial trajectory; that from 0, b to 1,0 is the 2,3 interfacial trajectory; that from 1,0 to 1,0 that does not go via 0, b is that of a 1,3 interface that is not wet by phase 2; while the composite trajectory that does go from 1,0 to 1,0 via 0, b is that of a 1,3 interface that is wet by phase 2. Each trajectory corresponds to a solution r 1 (z), r 2 (z) of the Euler Lagrange eqns (5) with the appropriate boundary conditions. When these are substituted in the integrand on the right-hand side of (4) the result is the tension s of the corresponding interface. The three tensions are functions of b. The equilibrium 1,3 interface is or is not wet by phase 2 according to whether the tension s 13 that corresponds to the direct 1,3 trajectory is respectively greater or less than the sum s 12 + s 23, which is the tension of the wet 1,3 interface. The wetting transition occurs at that b( ¼ b w ) for which s 13 of the non-wet interface equals s 12 + s 23. If the transition is of second order, there is no solution of the Euler Lagrange equations corresponding to a non-wet 1,3 interface when b b w [cf. Fig. 1(b)], hence no direct 1,3 trajectory in the r 1, r 2 plane, only the indirect trajectory via 0, b. The latter is always present because it is a composite of the 1,2 and 2,3 interfacial Fig. 2 Wetting and non-wetting trajectories in the r 1, r 2 plane. The bulk-phase densities are r 1 ¼ 1, r 2 ¼ 0 (phase 1); r 1 ¼ 0, r 2 ¼ b (phase 2); r 1 ¼ 1, r 2 ¼ 0 (phase 3). This journal is ª The Royal Society of Chemistry 2010 Faraday Discuss., 2010, 146, trajectories, and those interfaces are present irrespective of whether phase 2 wets the 1,3 interface or not. 3. First- and second-order transitions Wetting transitions, which here are transitions in the structure of the 1,3 interface, are in many ways analogous to the transitions between bulk phases in three dimensions. For the latter, while the free energy is necessarily continuous at the transition, its first derivatives with respect to the thermodynamic fields, viz., the energy, entropy, and mass densities and the chemical composition, are discontinuous for first-order transitions (except for accidental azeotropies), but are still continuous at transitions of higher order. Discontinuity in the thermodynamic densities reflects discontinuity in the structures of the phases; at higher-order transitions the structures are continuous. The relevant free energy in the wetting transitions is Ds. Its derivative with respect to the field variable b is discontinuous at b ¼ b w for the first-order transition in Fig. 1(a) but continuous for the second-order transition in Fig. 1(b). In our meanfield density-functional models the structures of the two alternative 1,3 surface phases are the alternative chemical composition profiles r 1 (z), r 2 (z) of the wet and non-wet 1,3 interfaces. These correspond to the two alternative 1,3 trajectories in the r 1, r 2 plane of Fig. 2, one going directly from 1,0 to 1,0 (non-wet) and the other indirectly, via 0, b (wet). If the transition is of first order the structure is discontinuous; i.e., the two alternative trajectories are still distinct at b ¼ b w. If it is of second order the two trajectories coincide at the transition, where there is then only the one, which goes via 0, b w. As noted earlier, the function F(r 1, r 2 ; b) in (3) and its derivatives vf/vr 1 and vf/r 2 are 0 at the points r 1, r 2 in Fig. 2 representing the compositions of the bulk phases, while F 0 everywhere else. Near those points contours of constant F are ellipses. For simplicity we suppose F(r 1, r 2 ; b) to be an even function of r 1, so the trajectories are symmetric in r 1, as in the example in Fig. 2. Then in the immediate neighborhood of the phase-2 point at 0, b we will have Fðr 1 ; r 2 ; bþ r1 a b r2 þ (6) pffiffiffiffi pffiffiffiffi where a 1 F and a2 F are the semi-axes of the elliptical contours of constant F in the directions of the r 1 and r 2 axes, respectively, with a 1 and a 2, in general, functions of b. In related work 2 it is shown from the Euler Lagrange eqns (5) and their first integral, and from (6), that when the wetting transition is of higher than first order, d 2 r 2 /dr 12 evaluated at the maximum on the non-wetting 1,3 trajectory in Fig. 2 must diverge at the wetting transition, while d 2 r 2 /dr 12 on the wetting trajectory near 0, b diverges or remains finite as b/b w according to whether, at b ¼ b w, a 1 /a 2 2 or a 1 /a 2 $ 2, respectively. Then since the wetting and non-wetting trajectories coincide at a higher-order wetting transition, a necessary condition for the transition to be of higher order is that a 1 /a 2 2 at b ¼ b w, and then a sufficient condition that it be of first order is that a 1 /a 2 $ 2. We shall now recall two wetting-transition models that have been treated in the past, 1,3 in one of which the transition is of first order and in the other of which it is of second order, and we shall see how the stated criterion is satisfied. In the first, 1 F(r 1, r 2 ;b) ¼ 16r 22 (r 2 b) 2 + [(r 2 br 1 ) 2 b 2 ] 2 + [(r 2 + br 1 ) 2 b 2 ] 2 (7) The coexisting phases are again of compositions ( 1,0), (0, b), and (1,0), as in Fig. 2. In the immediate p neighborhood of 0, b this F is again of the form (6), now with a 1 =a 2 ¼ ffiffi 3 =b. There is a first-order wetting transition at b ¼ bw x0.51, 1 where a Faraday Discuss., 2010, 146, This journal is ª The Royal Society of Chemistry 2010 then a 1 /a 2 x3.4. Thus, the wetting transition in this model had to be of first order, as it is. Indeed, the non-wetting trajectory in this model just runs along the r 1 axis itself, from 1,0 to 1,0, nowhere near the wetting trajectory, which passes through 0, b. In the second of the two earlier models, 3 F(r 1, r 2 ;b) ¼ [(r 1 + 1) 2 + r 22 ][(r 1 1) 2 + r 22 ][r 12 +(r 2 b) 2 ] (8) This, too, is of the form (6) near 0, b, now with a 1 /a 2 ¼ 1. It thus satisfies the necessary condition for its wetting transition to be of higher order, and indeed what is found is a classical second-order transition as in Fig. 1(b). (In this case the diverging second derivative d 2 r 2 /dr 12 at b ¼ b w is manifested as a discontinuous first derivative: b r 2 r 1.) From (6) it is apparent that the ratio a 1 /a 2 plays the role of an asymmetry, or anisotropy parameter, a 1 /a 2 ¼ 1 corresponding to the symmetric (isotropic) model. In general, the thermodynamic singularities at wetting are expected to be nonuniversal with respect to this ratio, in the sense that the critical exponents associated with higher than first-order wetting transitions vary continuously with a 1 /a 2. Indeed, in earlier work on models of wall wetting with a two-component order parameter 4 or anisotropic vector order parameter 5 this non-universality was demonstrated. Particularly interesting cases arise when the ratio a 1 /a 2 reflects the ratio of characteristic physical length scales of the problem. These two lengths can already appear in a single-component order parameter model, in the form of decay lengths of the wall-fluid and the fluid fluid potential, respectively. 6 In this model, the non-universality of the critical exponents at wetting, found at the classical (van der Waals theory) level, was found to be robust when thermal fluctuations are taken into account using functional renormalization group calculations. 7 In the entirely different physical system of a type-i superconductor, in which a wetting or interface delocalization transition was uncovered theoretically 8 and verified experimentally, 9 these two lengths are the superconducting coherence length x and the magnetic penetration depth l. For the case in which the theoretically predicted interface delocalization transition is critical, the critical exponent of the surface excess free energy was shown 10 to vary continuously with the Ginzburg Landau parameter k, which is the ratio of l to x. Now, returning to our model defined through (3) and (6), higher than first-order wetting transitions are found in the entire interval 0 a 1 /a 2 1. 2 This is again consistent with the criterion we exposed. Moreover, for these transitions, too, the critical exponents of the surface free energy at wetting are found to be non-universal and vary continuously with the asymmetry or anisotropy parameter a 1 /a 2 (ref. 2). It is of interest that the order of the wetting transition in mean-field density-functional models with excess free-energy density J of the form (3), with two spatially varying densities r 1 (z) and r 2 (z), depends on the structure of the function F only in the immediate neighborhood of the composition r 1, r 2 of the wetting phase. Acknowledgements We acknowledge support by a Grant-in-Aid for Scientific Research and the Next Generation Super Computing Project from MEXT, Japan; by grant G , FWO-Vlaanderen; and by the U.S. National Science Foundation. References 1 I. Szleifer and B. Widom, Mol. Phys., 1992, 75, K. Koga, J. O. Indekeu and B. Widom, Phys. Rev. Lett., 2010, 104, K. Koga and B. Widom, J. Chem. Phys., 2008, 128, E. H. Hauge, Phys. Rev. B: Condens. Matter, 1986, 33, C. J. Walden, B. L. Gyorffy and A. O. Parry, Phys. Rev. B: Condens. Matter, 1990, 42, T. Aukrust and E. H. Hauge, Phys. Rev. Lett., 1985, 54, This journal is ª The Royal Society of Chemistry 2010 Faraday Discuss., 2010, 146, 7 E. H. Hauge and K. Olaussen, Phys. Rev. B: Condens. Matter, 1985, 32, J. O. Indekeu and J. M. J. van Leeuwen, Phys. Rev. Lett., 1995, 75, 1618; J. O. Indekeu and J. M. J. van Leeuwen, Physica C, 1995, 251, V. F. Kozhevnikov, M. J. Van Bael, P. K. Sahoo, K. Temst, C. Van Haesendonck, A. Vantomme and J. O. Indekeu, New J. Phys., 2007, 9, J. M. J. van Leeuwen and E. H. Hauge, J. Stat. Phys., 1997, 87, Faraday Discuss., 2010, 146, This journal is ª The Royal Society of Chemistry 2010

We Need Your Support

Thank you for visiting our website and your interest in our free products and services. We are nonprofit website to share and download documents. To the running of this website, we need your help to support us.

Thanks to everyone for your continued support.

No, Thanks

SAVE OUR EARTH

We need your sign to support Project to invent "SMART AND CONTROLLABLE REFLECTIVE BALLOONS" to cover the Sun and Save Our Earth.

More details...Sign Now!

We are very appreciated for your Prompt Action!

x