Wedge filling and interface delocalization in finite Ising lattices with antisymmetric surface fields

Critical effects and scaling at meniscus osculation transitions

We propose a simple scaling theory describing critical effects at rounded meniscus osculation transitions which occur when the Laplace radius of a condensed macroscopic drop of liquid coincides with the local radius of curvature R_{w} in a confining parabolic geometry. We argue that the exponent β_{osc} characterizing the scale of the interfacial height ℓ_{0}∝R_{w}^{β_{osc}} at osculation, for large R_{w}, falls into two regimes representing fluctuation-dominated and mean-field-like behavior, respectively. These two regimes are separated by an upper critical dimension, which is determined here explicitly and depends on the range of the intermolecular forces. In the fluctuation-dominated regime, representing the universality class of systems with short-range forces, the exponent is related to the value of the interfacial wandering exponent ζ by β_{osc}=3ζ/(4-ζ). In contrast, in the mean-field regime, which was not previously identified and which occurs for systems with longer-range forces (and higher dimensions), the exponent β_{osc} takes the same value as the exponent β_{s}^{co} for complete wetting, which is determined directly by the intermolecular forces. The prediction β_{osc}=3/7 in d=2 for systems with short-range forces (corresponding to ζ=1/2) is confirmed using an interfacial Hamiltonian model which determines the exact scaling form for the decay of the interfacial height probability distribution function. A numerical study in d=3, based on a microscopic model density-functional theory, determines that β_{osc}≈β_{s}^{co}≈0.326 close to the predicted value of 1/3 appropriate to the mean-field regime for dispersion forces.

Meniscus osculation and adsorption on geometrically structured walls

We study the adsorption of simple fluids at smoothly structured, completely wet walls and show that a meniscus osculation transition occurs when the Laplace and geometrical radii of curvature of locally parabolic regions coincide. Macroscopically, the osculation transition is of fractional, 7/2, order and separates regimes in which the adsorption is microscopic, containing only a thin wetting layer, and mesoscopic, in which a meniscus exists. We develop a scaling theory for the rounding of the transition due to thin wetting layers and derive critical exponent relations that determine how the interfacial height scales with the geometrical radius of curvature. Connection with the general geometric construction proposed by Rascón and Parry is made. Our predictions are supported by a microscopic model density functional theory for drying at a sinusoidally shaped hard wall where we confirm the order of the transition and also an exact sum rule for the generalized contact theorem due to Upton. We show that as bulk coexistence is approached the adsorption isotherm separates into three regimes: A preosculation regime where it is microscopic, containing only a thin wetting layer; a mesoscopic regime, in which a meniscus sits within the troughs; and finally another microscopic regime where the liquid-gas interface unbinds from the crests of the substrate.

Capillary condensation and depinning transitions in open slits

We study the low-temperature phase equilibria of a fluid confined in an open capillary slit formed by two parallel walls separated by a distance L which are in contact with a reservoir of gas. The top wall of the capillary is of finite length H while the bottom wall is considered of macroscopic extent. This system shows rich phase equilibria arising from the competition between two different types of capillary condensation, corner filling, and meniscus depinning transitions depending on the value of the aspect ratio a=L/H and divides into three regimes: For long capillaries, with a<2/π, the condensation is of type I involving menisci which are pinned at the top edges at the ends of the capillary. For intermediate capillaries, with 2/π<a<1, depending on the value of the contact angle the condensation may be of type I or of type II, in which the menisci overspill into the reservoir and there is no pinning. For short capillaries, with a>1, condensation is always of type II. In all regimes, capillary condensation is completely suppressed for sufficiently large contact angles which is determined explicitly. For long and intermediate capillaries, we show that there is an additional continuous phase transition in the condensed liquid-like phase, associated with the depinning of each meniscus as they round the upper open edges of the slit. Meniscus depinning is third-order for complete wetting and second-order for partial wetting. Detailed scaling theories are developed for these transitions and phase boundaries which connect with the theories of wedge (corner) filling and wetting encompassing interfacial fluctuation effects and the direct influence of intermolecular forces. We test several of our predictions using a fully microscopic density functional theory which allows us to study the two types of capillary condensation and its suppression at the molecular level for different aspect ratios and contact angles.

Bridging of liquid drops at chemically structured walls

Using mesoscopic interfacial models and microscopic density functional theory we study fluid adsorption at a dry wall decorated with three completely wet stripes of width L separated by distances D_{1} and D_{2}. The stripes interact with the fluid with long-range forces inducing a large finite-size contribution to the surface free energy. We show that this nonextensive free-energy contribution scales with lnL and drives different types of bridging transition corresponding to the merging of liquid drops adsorbed at neighboring wetting stripes when the separation between them is molecularly small. We determine the surface phase diagram and show that this exhibits two triple points, where isolated drops, double drops, and triple drops coexist. For the symmetric case, D_{1}=D_{2}≡D, our results also confirm that the equilibrium droplet configuration always has the symmetry of the substrate corresponding to either three isolated drops when D is large or a single triple drop when D is small; however, symmetry-broken configurations do occur in a metastable part of the phase diagram which lies very close to the equilibrium-bridging phase boundary. Implications for phase transitions on other types of patterned surface are considered.

Influence of intermolecular forces at critical-point wedge filling

We use microscopic density functional theory to study filling transitions in systems with long-ranged wall-fluid and short-ranged fluid-fluid forces occurring in a right-angle wedge. By changing the strength of the wall-fluid interaction we can induce both wetting and filling transitions over a wide range of temperatures and study the order of these transitions. At low temperatures we find that both wetting and filling transitions are first order in keeping with predictions of simple local effective Hamiltonian models. However close to the bulk critical point the filling transition is observed to be continuous even though the wetting transition remains first order and the wetting binding potential still exhibits a small activation barrier. The critical singularities for adsorption for the continuous filling transitions depend on whether retarded or nonretarded wall-fluid forces are present and are in excellent agreement with predictions of effective Hamiltonian theory even though the change in the order of the transition was not anticipated.

Filling transitions in acute and open wedges

We present numerical studies of first-order and continuous filling transitions in wedges of arbitrary opening angle ψ, using a microscopic fundamental measure density functional model with short-ranged fluid-fluid forces and long-ranged wall-fluid forces. In this system the wetting transition characteristic of the planar wall-fluid interface is always first order regardless of the strength of the wall-fluid potential ɛ(w). In the wedge geometry, however, the order of the filling transition depends not only on ɛ(w) but also on the opening angle ψ. In particular we show that even if the wetting transition is strongly first order the filling transition is continuous for sufficient acute wedges. We show further that the change in the order of the transition occurs via a tricritical point as opposed to a critical end point. These results extend previous effective Hamiltonian predictions which were limited only to shallow wedges.

Filling and wetting transitions on sinusoidal substrates: a mean-field study of the Landau-Ginzburg model

We study the interfacial phenomenology of a fluid in contact with a one-dimensional array of infinitely long grooves of sinusoidal section, characterized by the periodicity length L and amplitude A. The system is modelled by the Landau-Ginzburg-Wilson functional, with fluid-substrate couplings which control the wettability of the substrate. We investigate the filling and wetting phenomena within the mean-field approximation, and compare with the predictions of the macroscopic and interfacial Hamiltonian theories. For large values of L and under bulk coexistence conditions, we observe first-order filling transitions between dry (D) and partially filled (F) interfacial states, and wetting transitions between partially filled F and completely wet (W) interfacial states of the same order as for the flat substrate. Depending on the order of the wetting transition, the transition temperature is either shifted towards lower temperatures for first-order wetting or it coincides with the wetting temperature on the flat substrate for continuous wetting. On the other hand, if the groove height is of order of the correlation length, only wetting transitions between D and W states are observed under bulk coexistence conditions. For this case, the transition temperature shift obeys approximately Wenzel's phenomenological law if the substrate favors first-order wetting, but it remains unshifted for continuous wetting. The borderline between the small and large L regimes correspond to a D - F - W triple point if wetting is first-order, and a D - F critical point for continuous wetting. Beyond bulk coexistence conditions, filling and first-order wetting transitions continue into off-coexistence filling and prewetting lines, which end up at critical points. Our findings show that the macroscopic theory only describes accurately the filling transition close to bulk coexistence and large L, while microscopic structure of the fluid is essential to understand wetting and filling away from bulk coexistence.

Filling and wetting transitions at grooved substrates

The wetting and filling properties of a fluid adsorbed on a solid grooved substrate are studied by means of a microscopic density functional theory. The grooved substrates are modelled using a solid slab, interacting with the fluid particles via long-range dispersion forces, to which a one-dimensional array of infinitely long rectangular grooves is sculpted. By investigating the effect of the groove periodicity and the width of the grooves and the ridges, a rich variety of different wetting morphologies is found. In particular, we show that for a saturated ambient gas, the adsorbent can occur in one of four wetting states characterized by (i) empty grooves, (ii) filled grooves, (iii) a formation of mesoscopic hemispherical caps (iv) a macroscopically wet surface. The character of the transition between particular regimes, that also extend off-coexistence, sensitively depends on the model geometry. The temperature at which the system becomes completely wet is considerably higher than that for a flat wall.

Critical behaviour of the Ising ferromagnet confined in quasi-cylindrical pores: a Monte Carlo study

The critical behaviour of the Ising ferromagnet confined in pores of radius R and length L is studied by means of Monte Carlo computer simulations. Quasi-cylindrical pores are obtained by replicating n-times a triangular lattice disc of radius R, where L = na and a is the spacing between consecutive replications. So, spins placed at the surface of the pores have less nearest-neighbours (NN) as compared to 8 NN for spins in the bulk. These "missing neighbour" effects undergone by surface spins cause a strong suppression of surface ordering, leading to an ordinary surface transition. Also, the effect propagates into the bulk for small tubes (R ≤ 12) and the effective critical temperature of the pores is shifted towards lower values than in the bulk case. By applying the standard finite-size scaling theory, subsequently supported by numerical data, we concluded that data collapse of relevant observables, e.g., magnetization (m), susceptibility, specific heat, etc., can only be observed by comparing simulation results obtained by keeping the aspect ratio C ≡ R∕L constant. Also, by extrapolating "effective" R-dependent critical temperatures to the thermodynamic limit (R → ∞, C fixed), we obtained T(C)(∞) = 6.208(4). As suggested by finite-size scaling arguments, the magnetization is measured at the critical point scales according to [|m|]Tc R(β/ν) is proportional to [R/L](1/2), where β and ν are the standard exponents for the order parameter and the correlation length, respectively. Furthermore, it is shown that close to criticality the axial correlation length decreases exponentially with the distance. That result is the signature of the formation of (randomly distributed) alternating domains of different magnetization, which can be directly observed by means of snapshot configurations, whose typical length (ξ) is given by the characteristic length of the exponential decay of correlations. Moreover, we show that at criticality ξ = 0.43(2)R.

Density functional study of complete, first-order and critical wedge filling transitions

We present numerical studies of complete, first-order and critical wedge filling transitions, at a right angle corner, using a microscopic fundamental measure density functional theory. We consider systems with short-ranged, cut-off Lennard-Jones, fluid-fluid forces and two types of wall-fluid potential: a purely repulsive hard wall and also a long-ranged potential with three different strengths. For each of these systems we first determine the wetting properties occurring at a planar wall, including any wetting transition and the dependence of the contact angle on temperature. The hard wall corner is completely filled by vapour on approaching bulk coexistence and the numerical results for the growth of the meniscus thickness are in excellent agreement with effective Hamiltonian predictions for the critical exponents and amplitudes, at leading and next-to-leading order. In the presence of the attractive wall-fluid interaction, the corresponding planar wall-fluid interface exhibits a first-order wetting transition for each of the interaction strengths considered. In the right angle wedge geometry the two strongest interactions produce first-order filling transitions while for the weakest interaction strength, for which wetting and filling occur closest to the bulk critical point, the filling transition is second-order. For this continuous transition the critical exponent describing the divergence of the meniscus thickness is found to be in good agreement with effective Hamiltonian predictions.

Critical point wedge filling

We present results of a microscopic density functional theory study of wedge filling transitions, at a right-angle wedge, in the presence of dispersionlike wall-fluid forces. Far from the corner the walls of the wedge show a first-order wetting transition at a temperature T(w) which is progressively closer to the bulk critical temperature T(c) as the strength of the wall forces is reduced. In addition, the meniscus formed near the corner undergoes a filling transition at a temperature T(f)<T(w), the value of which is found to be in excellent agreement with macroscopic predictions. We show that the filling transition is first order if it occurs far from the critical point but is continuous if T(f) is close to T(c) even though the walls still show first-order wetting behavior. For this continuous transition the distance of the meniscus from the apex grows as ℓ(w)≈(T(f)-T)(-β(w)) with the critical exponent β(w)≈0.46±0.05 in good agreement with the phenomenological effective Hamiltonian prediction. Our results suggest that critical filling transitions, with accompanying large scale universal interfacial fluctuation effects, are more generic than thought previously, and are experimentally accessible.

Capillary condensation in a square geometry with surface fields

We study the influence of wetting on capillary condensation for a simple fluid in a square geometry with surface fields, where the reference system is an infinitely long slit. The corner transfer matrix renormalization group method has been extended to study a two-dimensional Ising model confined in an L × L geometry with equal surface fields. Our results have confirmed that in both geometries the coexistence line shift is governed by the same scaling powers, but their prefactors are different.

Wetting transition in the two-dimensional Blume-Capel model: a Monte Carlo study

The wetting transition of the Blume-Capel model is studied by a finite-size scaling analysis of L×M lattices where competing boundary fields ±H_{1} act on the first or last row of the L rows in the strip, respectively. We show that using the appropriate anisotropic version of finite-size scaling, critical wetting in d=2 is equivalent to a "bulk" critical phenomenon with exponents α=-1, β=0, and γ=3. These concepts are also verified for the Ising model. For the Blume-Capel model, it is found that the field strength H_{1c}(T) where critical wetting occurs goes to zero when the bulk second-order transition is approached, while H_{1c}(T) stays nonzero in the region where in the bulk a first-order transition from the ordered phase, with nonzero spontaneous magnetization, to the disordered phase occurs. Interfaces between coexisting phases then show interfacial enrichment of a layer of the disordered phase which exhibits in the second-order case a finite thickness only. A tentative discussion of the scaling behavior of the wetting phase diagram near the tricritical point is also given.

The order of filling transitions in acute wedges

Using a square-gradient density functional model we test the prediction that the filling transition for a fluid in a wedge geometry changes from continuous to first-order as the wedge becomes more acute. Our numerical findings confirm such a change of order, but the value of the tilt angle at which it occurs, α* ≈ 45°, is considerably smaller than the original theoretical prediction. We critically reassess this work, which was based on allowing for the self-interaction of the fluid interface, and argue that the interfacial curvature and effective wavevector dependent surface tension can further lower the predicted value of α*, in keeping with our numerical findings. Interfacial fluctuation effects, occurring beyond mean-field level, are also discussed using effective Hamiltonian theory and are shown to substantially increase the value of α*.

Scaling properties of fluid adsorption near the base of a cylinder

We consider the adsorption of fluid at the foot of a cylinder that protrudes from a flat substrate made of the same material. Provided the contact angle θ is small enough, a drop of liquid condenses near the base, the size of which can be determined using simple macroscopic arguments. The adsorption in this geometry shows scaling behavior related to a number of different interfacial phase transitions and, for systems with short-ranged forces, shows a remarkable property; for small θ, the height of the drop (measured from the base) and the width (measured from the cylinder axis) are near identical to expressions for the thickness and parallel correlation length for microscopic wetting films (at planar walls). The only difference is that the bulk correlation length is replaced by the radius of the cylinder. By taking into account the correct singular behavior of the line tension we show that this geometrical amplification of the microscopic lengths occurs for second-order, first-order, and complete wetting transitions, and is specific to three dimensions. Similar phenomena occurs for long-ranged forces, and shows crossover scaling behavior.

Confinement effects on phase behavior of soft matter systems

When systems that can undergo phase separation between two coexisting phases in the bulk are confined in thin film geometry between parallel walls, the phase behavior can be profoundly modified. These phenomena shall be described and exemplified by computer simulations of the Asakura-Oosawa model for colloid-polymer mixtures, but applications to other soft matter systems (e.g. confined polymer blends) will also be mentioned. Typically a wall will prefer one of the phases, and hence the composition of the system in the direction perpendicular to the walls will not be homogeneous. If both walls are of the same kind, this effect leads to a distortion of the phase diagram of the system in thin film geometry, in comparison with the bulk, analogous to the phenomenon of "capillary condensation" of simple fluids in thin capillaries. In the case of "competing walls", where both walls prefer different phases of the two phases coexisting in the bulk, a state with an interface parallel to the walls gets stabilized. The transition from the disordered phase to this "soft mode phase" is rounded by the finite thickness of the film and is not a sharp phase transition. However, a sharp transition can occur where this interface gets localized at (one of) the walls. The relation of this interface localization transition to wetting phenomena is discussed. Finally, an outlook to related phenomena is given, such as the effects of confinement in cylindrical pores on the phase behavior, and more complicated ordering phenomena (lamellar mesophases of block copolymers or nematic phases of liquid crystals under confinement).

Fluid bridges confined between chemically nanopatterned solid substrates

We discuss equilibrium properties of classical fluids confined to nanoscopic volumes by solid substrates. The substrates themselves are endowed with wettable chemical patterns of variable symmetry. We develop a thermodynamic description suitable for these highly anisotropic systems. Based upon a combination of Monte Carlo simulations in the grand canonical ensemble and lattice density functional theory at mean-field level we analyze the structure and phase behaviour of the confined fluid. Under suitable thermodynamic conditions the fluid may condense partially in regions controlled by the wettable nanopatterns. The resulting fluid bridges are established as thermodynamic phases and exhibit unique rheological features.

Molecular dynamics simulation study on controlling the adsorption behavior of polyethylene by fine tuning the surface nanodecoration of graphite

Molecular dynamics simulations are applied to study the adsorption of polyethylene with different chain lengths on patterned graphite surfaces that contain nanoscale protrusions. The influence of the nanostructure on the strong attractive interaction inherently in the hydrophobic polyethylene and hydrophobic graphite system is investigated by modifying the top surface area and the height and the shape of the protrusions. The results are analyzed in terms of the chain configuration, the adsorption energy, the global orientational order parameter, and the normalized surface-chain contacting pair number in the first adsorption layer. When the size of the protrusion increases, the adsorption energy, the order parameter, and the normalized surface-chain contacting pair number decrease at a fixed chain length. When the size of the protrusion is fixed, the average adsorption energy per monomer and the order parameter decrease with increasing chain length because of the stronger intramolecular interactions between the monomers. Changing the protrusion shape in a suitable way will effectively reduce the strong surface-chain interaction.

Extended wedge covariance for wetting and filling transitions

Fluid adsorption on nonplanar and heterogeneous substrates is studied using a simple interfacial model. For systems with short-ranged forces, we find that, by tuning the local strength of the substrate potential, it is possible to find the exact equilibrium interfacial profile as a functional of the wall shape psi x. The tuning of the local substrate potential takes the form of a gauge condition theta x=+/-psi x, where theta x can be interpreted as a local effective contact angle. For wedgelike geometries with asymptotic tilt angle alpha, the midpoint interfacial height and roughness satisfy the same covariance relations previously found for simple linear wedges. For troughlike geometries satisfying the gauge condition, covariance is also found for the two-point correlation function. Predictions for more microscopic Landau and Ising models are also discussed.

Liquids on topologically nanopatterned surfaces

We report here surface x-ray scattering studies of the adsorption of simple hydrocarbon liquid films on nanostructured surfaces-silicon patterned by an array of nanocavities. Two different regimes, filling and growing, are observed for the wetting film evolution as a function of the chemical potential offset from the bulk liquid-vapor coexistence. The strong influence of geometrical effects is manifested by a dependence of liquid adsorption in the nanocavities that is stronger than the van der Waals behavior for flat surfaces. The observed dependence is, however, much weaker than predicted for the infinitely deep parabolic cavities, suggesting that the finite-size effects contribute significantly to the observed adsorption behavior.

Histogram analysis as a method for determining the line tension of a three-phase contact region by Monte Carlo simulations

A method is proposed for determining the line tension, which is the main physical characteristic of a three-phase contact region, by Monte Carlo (MC) simulations. The key idea of the proposed method is that if a three-phase equilibrium involves a three-phase contact region, the probability distribution of states of a system as a function of two order parameters depends not only on the surface tension, but also on the line tension. This probability distribution can be obtained as a normalized histogram by appropriate MC simulations, so one can use the combination of histogram analysis and finite-size scaling to study the properties of a three phase contact region. Every histogram and results extracted therefrom will depend on the size of the simulated system. Carrying out MC simulations for a series of system sizes and extrapolating the results, obtained from the corresponding series of histograms, to infinite size, one can determine the line tension of the three phase contact region and the interfacial tensions of all three interfaces (and hence the contact angles) in an infinite system. To illustrate the proposed method, it is applied to the three-dimensional ternary fluid mixture, in which molecular pairs of like species do not interact whereas those of unlike species interact as hard spheres. The simulated results are in agreement with expectations.

Effect of fluid-substrate attraction and pore geometry on fluid adsorption

We employ grand canonical ensemble Monte Carlo simulations to investigate the impact of substrate curvature on the phase behavior of an adjacent fluid. The substrates consist of a periodic sequence of grooves in the x direction; the grooves are infinitely long in the y direction. The shape of the grooves is controlled by a parameter eta. For eta = 0 the substrates are planar. If eta = 1, the grooves are wedge shaped. If eta > 1 the grooves become concave and in the limit eta = infinity rectangular. The fluid-substrate potential representing a groove consists of two contributions, namely, that of the homogeneous substrate base corresponding to a semi-infinite solid and that of a finite piece of solid with nonplanar surfaces. Whereas the former contribution can be calculated analytically, the latter needs to be evaluated numerically. For very large values of eta, that is in (almost) rectangular grooves, we observe capillary condensation of that portion of fluid located inside the grooves. As eta decreases capillary condensation gives way to continuous filling. In all cases, a nearly planar film-gas interface eventually forms in the direction normal to the surface of the substrate base and outside the grooves if one increases the chemical potential sufficiently.

Phase transitions in nanosystems caused by interface motion: the Ising bipyramid with competing surface fields

The phase behavior of a large but finite Ising ferromagnet in the presence of competing surface magnetic fields +/-H(s) is studied by Monte Carlo simulations and by phenomenological theory. Specifically, the geometry of a double pyramid of height 2L is considered, such that the surface field is positive on the four upper triangular surfaces of the bipyramid and negative on the lower ones. It is shown that the total spontaneous magnetization vanishes (for L --> infinity) at the temperature Tf(H), related to the "filling transition" of a semi-infinite pyramid, which can be well below the critical temperature of the bulk. The discontinuous vanishing of the magnetization is accompanied by a susceptibility that diverges with a Curie-Weiss power law, when the transition is approached from either side. A Landau theory with size-dependent critical amplitudes is proposed to explain these observations, and confirmed by finite size scaling analysis of the simulation results. The extension of these results to other nanosystems (gas-liquid systems, binary mixtures, etc.) is briefly discussed.