Molecular dynamics simulation of nanodroplet spreading enhanced by linear surfactants

We utilize molecular dynamics simulations to probe the surfactant-mediated spreading of a Lennard-Jones liquid droplet on a solid surface. The surfactants are linear hexamers that are insoluble in the liquid and reduce the surface tension of the liquid-vapor interface. We study how the interaction of the surfactant hexamers with the solid substrate influences spreading, as well as the dependence of spreading on surfactant concentration. We find that the spreading speed is strongly influenced by the attraction of the hydrophobic surfactant tail to the solid surface. When this attraction is sufficiently strong, surfactant molecules partition to the liquid-solid interface and facilitate spreading. This partitioning can lead to an inhomogeneous distribution of surfactant over the liquid-vapor interface, which could drive the Marangoni convection. We also observe that the surfactant molecules can assemble into micelles on the solid surface. The repulsion between micelles at the liquid-solid interface can lead to break-off and migration of the micelles from the liquid-solid to the gas-solid interface and spreading is facilitated in this way. Our model system contains features that are believed to underlie superspreading in experimental studies of droplet spreading.     

Capillary force required to detach micron-sized particles from solid surfaces—Validation with bubbles circulating in water and 2 μm-diameter latex spheres

The adhesion forces holding micron-sized particles to solid surfaces can be studied through the detachment forces developed by the transit of an air–liquid interface in a capillary. Two key variables affect the direction and magnitude of the capillary detachment force: (i) the thickness of the liquid film between the bubble and the capillary walls, and (ii) the effective angle of the triple phase contact between the particles and the interface. Variations in film thickness were calculated using a two-phase flow model. Film thickness was used to determine the time-variation of the capillary force during transit of the bubble. The curve for particle detachment was predicted from the calculated force. This curve proved to be non-linear and gave in situ information on the effective contact angle developing at the particle–bubble interface during detachment. This approach allowed an accurate determination of the detachment force. This theoretical approach was validated using latex particles 2 μm in diameter.     

Anomalous Capillary Pressure, Stress, and Stability of Solids-Coated Bubbles

A two-dimensional theoretical model for solids-coated, or "armored," bubbles shows how the armor can support a liquid-vapor interface of reduced or reversed curvature between the particles, giving the bubble zero or even negative capillary pressure. The inward capillary force pulling the particles into the center of the bubble are balanced by large contact forces between the particles in the armor. Thus the bubble is stabilized against dissolution of gas into surrounding liquid, which otherwise would rapidly collapse the bubble. The stresses between particles in such cases are large and could drive sintering of the particles into a rigid framework. Earlier work on solids-coated bubbles assumed that solids can freely enter or leave the bubble surface as the bubble shrinks or expands. In such a case, armored bubbles would not be stable to gas dissolution into surrounding liquid. A new free-energy analysis, however, suggests that a shrunken bubble would not spontaneously expel a solid particle from its armor to relieve stress and allow the bubble to shrink further. Implications and limitations of the theory are discussed. Copyright 1999 Academic Press.     

Capillary waves at the liquid-vapor interface and the surface tension of water

Capillary waves occurring at the liquid-vapor interface of water are studied using molecular dynamics simulations. In addition, the surface tension, determined thermodynamically from the difference in the normal and tangential pressure at the liquid-vapor interface, is compared for a number of standard three- and four-point water models. We study four three-point models (SPC/E, TIP3P, TIP3P-CHARMM, and TIP3P-Ew) and two four-point models (TIP4P and TIP4P-Ew). All of the models examined underestimate the surface tension; the TIP4P-Ew model comes closest to reproducing the experimental data. The surface tension can also be determined from the amplitude of capillary waves at the liquid-vapor interface by varying the surface area of the interface. The surface tensions determined from the amplitude of the logarithmic divergence of the capillary interfacial width and from the traditional thermodynamic method agree only if the density profile is fitted to an error function instead of a hyperbolic tangent function.     

Universal Two-dimensional forces that Act on Particles at interfaces

This paper starts with a short tribute to the scientific legacy of Peter Kralchevsky and his role in the collaboration with the research group of Kuniaki Nagayama in Japan. Next it presents an overview of the lateral capillary forces, studied jointly by both authors and their groups. Analogies with some other forces in nature are highlighted where relevant. The lateral capillary forces emerge between objects dispersed at a liquid interface or in liquid film and act in direction parallel to the liquid interface, being caused by the liquid surface tension. The lateral capillary forces depend on the size of the objects and the wetting of the particle by the liquid and can be either attractive or repulsive. For fine particles the force–distance relationship follows an inverse-linear law. Their effect is seen in our everyday experience that dust particles and bubbles tend to aggregate on the water surface. These forces are universally acting between any objects regardless of scale and are described by equations similar to those describing the gravitational and electric forces in the two-dimensional (2D) world. The current KN & PK review article provides a didactic overview that is easy to digest as a 2D model of gravity in general relativity, because the origin of the latter force, namely the space distortion, can be clearly visualized.     

Effect of flexibility on surface tension and coexisting densities of water

Molecular dynamics simulations of pure water at the liquid-vapor interface are performed using direct simulation of interfaces in a liquid slab geometry. The effect of intramolecular flexibility on coexisting densities and surface tension is analyzed. The dipole moment profile across the liquid-vapor interface shows different values for the liquid and vapor phases. The flexible model is a polarizable model. This effect is minor for liquid densities and is large for surface tension. The liquid densities increase from 2% at 300 K to 9% at 550 K when the force field is changed from a fully rigid simple point charge extended (SPCE) model to that of a fully flexible model with the same intermolecular interaction parameters. The increases in surface tension at both temperatures are around 11% and 36%, respectively. The calculated properties of the flexible models are closer to the experimental data than those of the rigid SPCE. The effect of the maximum number of reciprocal vectors (h(z) (max)) and the surface area on the calculated properties at 300 K is also analyzed. The coexiting densities are not sensitive to those variables. The surface tension fluctuates with h(z) (max) with an amplitude larger than 10 mN m(-1). The effect of using small interfacial areas is slightly larger than the error in the simulations.     

Anchoring energy for the nematic liquid crystal-Langmuir Blodgett film interface

We analyse the effect of a Langmuir-Blodgett (LB) multilayer on the surface properties of a nematic liquid crystal (NLC). We show that the easy axis of the LB-NLC interface coincides with that of the LB-solid surface interface. On the contrary, the effective anchoring energy of the LB-NLC interface is lower than that associated with the LB-solid substrate interface. We show in a first approximation that the anchoring energy characterizing the NLC may be separated into three contributions: one connected with the interaction between the LB film and the solid substrate, one due to the direct LB-NLC interaction and the other one having an elastic origin. Nevertheless, to be more precise, one has to consider also the term associated with the interaction energy between the NLC and the substrate, which is screened by the LB film. The elastic contribution is of the order of the elastic constant of the LB film over the thickness of the multilayer. This quantity is estimated to be of the order of 10^-2-10^-1 erg cm^-2, as experimentally observed. Possible extensions of our model are also discussed.     

Interactions between particles with an undulated contact line at a fluid interface: capillary multipoles of arbitrary order

A colloidal particle adsorbed at a fluid interface could have an undulated, or irregular contact line in the presence of surface roughness and/or chemical inhomogeneity. The contact-line undulations produce distortions in the surrounding liquid interface, whose overlap engenders capillary interaction between the particles. The convex and concave local deviations of the meniscus shape from planarity can be formally treated as positive and negative "capillary charges," which form "capillary multipoles." Here, we derive theoretical expressions for the interaction between two capillary multipoles of arbitrary order. Depending on the angle of mutual orientation, the interaction energy could exhibit a minimum, or it could represent a monotonic attraction. For undulation amplitudes larger than 5 nm, the interaction energy is typically much greater than the thermal energy kT. As a consequence, a monolayer from capillary multipoles exhibits considerable shear elasticity, and such monolayer is expected to behave as a two-dimensional elastic solid. These theoretical results could be helpful for the understanding of phenomena related to aggregation and ordering of particles adsorbed at a fluid interface, and for the interpretation of rheological properties of particulate monolayers. Related research fields are the particle-stabilized (Pickering) emulsions and the two-dimensional self-assembly of microscopic particles.     

Direct measurement of the viscous force between two spherical particles trapped in a thin wetting film

Here we present the first direct measurement of the viscous drag force between two spherical particles of millimeter size trapped in a thin wetting film. Each particle is constrained by the liquid/air interface and the solid substrate. The viscous force is counterbalanced by another known force, the attractive capillary immersion force between identical particles protruding from the film surface. The results of the measurements provide evidence for an increased hydrodynamic force due to a non-Stokesian resistance to the particle motion. Our findings can be applied to the self-assembly of colloidal particles in a two-dimensional array for coating and to the friction between small species and a solid. Received: 19 March 1999 Accepted in revised form: 11 May 1999     

Anomalous pH dependent stability behavior of surfactant-free nonpolar oil drops in aqueous electrolyte solutions

Recent advances in atomic force microscopy (AFM) force measurement techniques have allowed the direct measurement and theoretical interpretation of the interaction between a liquid droplet and a solid surface or between two liquid droplets. In this study, we investigated the interaction across an aqueous thin film between fluorocarbon (perfluoropentane) droplets, hydrocarbon (tetradecane) droplets, and a droplet and a flat mica surface in the absence of stabilizers. It was found that even at a relatively elevated electrolyte concentration of 0.1 M NaNO3, depending on the solution pH, interactions between two identical droplets or a droplet and a mica surface could be repulsive. A simple theoretical analysis of the magnitude and range of these interactive forces suggests that the DLVO theory cannot explain the observed behavior. The measured force behavior is discussed in the context of ion adsorption, and the arising charging effects, at the bare oil-water interface.     

Molecular simulation of fluid-solid interfaces at nanoscale

The equilibrium states of vapor and liquid coexistent phases in contact with a solid surface are studied at the nanoscale by molecular dynamics simulations for a temperature close to the fluid triple point. The characteristics of the solid-fluid interfaces are determined when the interaction strength between the fluid and the solid varies in order to go from a situation of complete drying to that of complete wetting. From the vapor-liquid density profiles of liquid drops lying on the substrate surface or menisci of liquid films confined in slit pores, the contact angles made by the vapor-liquid interface with the solid are computed. The angle values are similar for the drops and the films. They are also in good qualitative agreement with the estimates obtained through the Young's relation from the surface tensions associated with the vapor-solid, liquid-solid, and vapor-liquid interfaces. However, at this scale, the uncertainties inherent to the angle computation and, to a lesser extent, to size effects seem to preclude that the quantitative agreement between the angle estimates obtained from the interface geometry and calculated from the Young's relation can be better than few degrees.     

Alignment of liquid crystals using Langmuir‒Blodgett films of unsymmetrical bent-core liquid crystals

ABSTRACT

The properties of the thin films of liquid crystal (LC) molecules can be governed easily by external fields. The anisotropic structure of the LC molecules has a large impact on the electrical and optical properties of the film. The Langmuir monolayer (LM) of LC molecules at the air–water interface is known to exhibit a variety of surface phases which can be transferred onto a solid substrate using the Langmuir?Blodgett (LB) technique. Here, we have studied the LM and LB films of asymmetrically substituted bent-core LC molecules. The morphology of LB film of the molecules is found to be a controlling parameter for aligning bulk LC in the nematic phase. It was found that the LB films of the bent-core molecules possessing defects favour the planar orientation of nematic LC, whereas the LB films with fewer defects show homeotropic alignment. The defect in LB films may introduce splay or bend distortions in the nematic near the alignment layer which can govern the planar alignment of the bulk LC. The uniform layer of LB film facilitates the molecules of nematic to anchor vertically due to a strong van der Waals interaction between the aliphatic chains leading to a homeotropic alignment.     

Influences of hydration force and elastic strain energy on the stability of solid film in a very thin solid-on-liquid structure

Since hydration forces become very strong at short range and are particularly important for determining the magnitude of the adhesion between two surfaces or interaction energy, the influences of the hydration force and elastic strain energy due to hydration-induced layering of liquid molecules close to a solid film surface on the stability of a solid film in a solid-on-liquid (SOL) nanostructure are studied in this paper. The liquid of this thin SOL structure is a kind of water solution. Since the surface forces play an important role in the structure, the total free energy change of SOL structures consists of the changes in the bulk elastic energy within the solid film, the surface energy at the solid-liquid interface and the solid-air interface, and highly nonlinear volumetric component associated with interfacial forces. The critical wavelength of one-dimensional undulation, the critical thickness of the solid film, and the critical thickness of the liquid layer are studied, and the stability regions of the solid film have been determined. Emphasis is placed on calculation of critical values, which are the basis of analyzing the stability of the very thin solid film.     

Modeling Pendular Liquid Bridges with a Reducing Solid-Liquid Interface

Liquid bridges formed between particles of dissimilar surface properties are important in many processes involving the handling of powders in mixtures. For instance, growth kinetic models for wet granulation frequently incorporate the evolution and resistance to breakage of individual liquid bridges between particles in a statistical form. These models generally propose a confusing definition of liquid-to-solid contact angles. Taken as a single thermodynamic value, they typically neglect the influence of wetting hysteresis on the liquid bridge. In this paper, a simple model based on the interfacial energies is proposed for the evolution of liquid bridges when one solid-liquid interface reduces. This receding process is well described by a balance between the adhesion energy of the bridge liquid on the particle surface and the capillary energy stored by the liquid free surface. The extent of solid-liquid interfacial area reduction can hence be predicted from the initial liquid bridge configuration. The liquid bridge shape is approximated by a parabolic curve, which is validated from the good agreement between measured and calculated contact angles or liquid-vapor interfacial area. Copyright 2000 Academic Press.     

Probing chitosan and phospholipid interactions using Langmuir and Langmuir-Blodgett films as cell membrane models

The interaction between chitosan and Langmuir and Langmuir-Blodgett (LB) films of dimyristoyl phosphatidic acid (DMPA) is investigated, with the films serving as simplified cell membrane models. At the air-water interface, chitosan modulates the structural properties of DMPA monolayers, causing expansion and decreasing the monolayer elasticity. As the surface pressure increased, some chitosan molecules remained at the interface, but others were expelled. Chitosan could be transferred onto solid supports alongside DMPA using the LB technique, as confirmed by infrared spectroscopy and quartz crystal microbalance measurements. The analysis of sum-frequency vibration spectroscopy data for the LB films combined with surface potential measurements for the monolayers pointed to chitosan inducing the ordering of the DMPA alkyl chains. Furthermore, the morphology of DMPA LB films, studied with atomic force microscopy, was affected significantly by the incorporation of chitosan, with the mixed chitosan-DMPA films displaying considerably higher thickness and roughness, in addition to chitosan aggregates. Because chitosan affected DMPA films even at pressures characteristic of cell membranes, we believe this study may help elucidate the role of chitosan in biological systems.     

Self-organization of a wedge-shaped surfactant in monolayers and multilayers

The self-organization behavior of a wedge-shaped surfactant, disodium-3,4,5-tris(dodecyloxy)phenylmethylphosphonate, was studied in Langmuir monolayers (at the air-water interface), Langmuir-Blodgett (LB) monolayers and multilayers, and films adsorbed spontaneously from isooctane solution onto a mica substrate (self-assembled films). This compound forms an inverted hexagonal lyotropic liquid crystal phase in the bulk and in thick adsorbed films. Surface pressure isotherm and Brewster angle microscope (BAM) studies of Langmuir monolayers revealed three phases: gas (G), liquid expanded (LE), and liquid condensed (LC). The surface pressure-temperature phase diagram was determined in detail; a triple point was found at approximately 10 degrees C. Atomic force microscope (AFM) images of LB monolayers transferred from various regions of the phase diagram were consistent with the BAM images and indicated that the LE regions are approximately 0.5 nm thinner than the LC regions. AFM images were also obtained of self-assembled films after various adsorption times. For short adsorption times, when monolayer self-assembly was incomplete, the film topography indicated the coexistence of two distinct monolayer phases. The height difference between these two phases was again 0.5 nm, suggesting a correspondence with the LE/LC coexistence observed in the Langmuir monolayers. For longer immersion times, adsorbed multilayers assembled into highly organized periodic arrays of inverse cylindrical micelles. Similar periodic structures, with the same repeat distance of 4.5 nm, were also observed in three-layer LB films. However, the regions of organized periodic structure were much smaller and more poorly correlated in the LB multilayers than in the films adsorbed from solution. Collectively, these observations indicate a high degree of similarity between the molecular organization in Langmuir layers/LB films and adsorbed self-assembled films. In both cases, monolayers progress through an LE phase, into LE/LC coexistence, and finally into LC phase as surface density increases. Following the deposition of an additional bilayer, the film reorganizes to form an array of inverted cylindrical micelles.     

十二烷基苯磺酸钠在SiO2表面聚集的分子动力学模拟

采用分子动力学方法研究了阴离子表面活性剂十二烷基苯磺酸钠(SDBS)在无定形SiO2固体表面的吸附. 设置不同的水层厚度, 观察固液界面和气液界面吸附的差异. 模拟发现表面活性剂分子能够在短时间内吸附到SiO2表面, 受碳链和固体表面之间相互作用的影响形成表面活性剂分子层, 并依据吸附量的大小形成不同的聚集结构; 在水层足够厚的情况下, 由于有较多的表面活性剂分子吸附在固体表面,从而形成带有疏水核心的半胶束结构; 计算得到的成对势表明极性头与钠离子或水分子之间的结合或解离与二者之间的能垒有关, 解离能垒远大于结合能垒, 引起更多Na+聚集在极性头周围而只有少数Na+存在于溶液中; 无论气液还是固液界面, 极性头均伸向水相, 与水分子形成不同类型的氢键. 模拟表明, 分子动力学方法可以作为实验的一种补充, 为实验提供必要的微观结构信息.