Parameter passing between molecular dynamics and continuum models for droplets on solid substrates: The static case

Applicability of the thermodynamic and mechanical route to Young's equation for rigid and flexible solids: A molecular dynamics simulations study of a Lennard-Jones system model

The wetting properties of a liquid in contact with a solid are commonly described by Young's equation, which defines the relationship between the angle made by a fluid droplet onto the solid surface and the interfacial properties of the different interfaces involved. When modeling such interfacial systems, several assumptions are usually made to determine this angle of contact, such as a completely rigid solid or the use of the tension at the interface instead of the surface free energy. In this work, we perform molecular dynamics simulations of a Lennard-Jones liquid in contact with a Lennard-Jones crystal and compare the contact angles measured from a droplet simulation with those calculated using Young's equation based on surface free energy or surface stress. We analyze cases where the solid atoms are kept frozen in their positions and where they are allowed to relax and simulate surfaces with different wettability and degrees of softness. Our results show that using either surface free energy or surface stress in Young's equation leads to similar contact angles but different interfacial properties. We find that the approximation of keeping the solid atoms frozen must be done carefully, especially if the liquid can efficiently pack at the interface. Finally, we show that to correctly reproduce the measured contact angles when the solid becomes soft, the quantity to be used in Young's equation is the surface free energy only and that the error committed in using the surface stress becomes larger as the softness of the solid increases.

Exploration of contact angle hysteresis mechanisms: From microscopic to macroscopic

Variations from equilibrium Young's angle, known as contact angle hysteresis (CAH), are frequently observed upon droplet deposition on a solid surface. This ubiquitous phenomenon indicates the presence of multiple local surface energy minima for the sessile droplet. Previous research primarily explains CAH via considering macroscopic roughness, such as topographical defects, which alter the effective interfacial energy between the fluid phase and the solid phase, thereby shifting the global surface energy minimum. One typical example is the classic Cassie-Baxter-Wenzel theory. Here, we propose an alternative microscopic mechanism that emphasizes the complexity of molecular rearrangements at the fluid-solid interface, treating their interfacial tensions as variables, which results in multiple local surface energy minima. Our theoretical framework demonstrates that CAH can occur even on chemically homogeneous and mechanically smooth-flat substrates, aligning with previously unexplained experimental observations. In addition, we explore the interplay between macroscopic and microscopic roughness in influencing CAH and clarify the contrasting wetting behaviors-the lotus effect and the rose petal effect-on hierarchical roughness from a thermodynamic perspective. This work provides valuable insights into surface tension determination by restoring the natural physical properties of interfaces and illuminates the multifaceted mechanisms underlying the everyday occurrences of CAH.

Microscopic derivation of the thin film equation using the Mori-Zwanzig formalism

The hydrodynamics of thin films is typically described using macroscopic models whose connection to the microscopic particle dynamics is a subject of ongoing research. Existing methods based on density functional theory provide a good description of static thin films but are not sufficient for understanding nonequilibrium dynamics. In this work, we present a microscopic derivation of the thin film equation using the Mori-Zwanzig projection operator formalism. This method allows to directly obtain the correct gradient dynamics structure along with microscopic expressions for mobility and free energy. Our results are verified against molecular dynamics simulations for both simple fluids and polymers.

Role of Interaction Range and Buoyancy on the Adhesion of Vesicles

Vesicles on substrates play a fundamental role in many biological processes, ranging from neurotransmitter release at the synapse on small scales to the nutrient intake of trees by large vesicles. For these processes, the adsorption or desorption of vesicles to biological substrates is crucial. Consequently, it is important to understand the factors determining whether and for how long a vesicle adsorbs to a substrate and what shape it will adopt. Here, we systematically study the adsorption of a vesicle to planar substrates with short- and long-range interactions, with and without buoyancy. We assume an axially symmetric system throughout our simulations. Previous studies often considered a contact potential of zero range and neutral buoyancy. The interaction range alters the location and order of the adsorption transition and is particularly important for small vesicles, e.g., in the synapse. Whereas even small density differences between the inside and the outside of the vesicle give rise to strong buoyancy effects for large vesicles, e.g., giant unilamellar vesicles, as buoyancy effects scale with the fourth power of the vesicle size. We find that (i) an attractive membrane-substrate potential with nonzero spatial extension leads to a pinned state, where the vesicle benefits from the attractive membrane-substrate interaction without significant deformation. The adsorption transition is of first order and occurs when the substrate switches from repulsive to attractive. (ii) Buoyancy shifts the transversality condition, which relates the maximal curvature in the contact zone to the adhesion strength and bending rigidity, up/downward, depending on the direction of the buoyancy force. The magnitude of the shift is influenced by the range of the potential. For upward buoyancy, adsorbed vesicles are at most metastable. We determine the stability limit and the desorption mechanisms and compile the thermodynamic data into an adsorption diagram. Our findings reveal that buoyancy, as well as spatially extended interactions, are essential when quantitatively comparing experiments to theory.

Binding potential and wetting behavior of binary liquid mixtures on surfaces

We present a theory for the interfacial wetting phase behavior of binary liquid mixtures on rigid solid substrates, applicable to both miscible and immiscible mixtures. In particular, we calculate the binding potential as a function of the adsorptions, i.e., the excess amounts of each of the two liquids at the substrate. The binding potential fully describes the corresponding interfacial thermodynamics. Our approach is based on classical density functional theory. Binary liquid mixtures can exhibit complex bulk phase behavior, including both liquid-liquid and vapor-liquid phase separation, depending on the nature of the interactions among all the particles of the two different liquids, the temperature, and the chemical potentials. Here we show that the interplay between the bulk phase behavior of the mixture and the properties of the interactions with the substrate gives rise to a wide variety of interfacial phase behaviors, including mixing and demixing situations. We find situations where the final state is a coexistence of up to three different phases. We determine how the liquid density profiles close to the substrate change as the interaction parameters are varied and how these determine the form of the binding potential, which in certain cases can be a multivalued function of the adsorptions. We also present profiles for sessile droplets of both miscible and immiscible binary liquids.

Stress-Induced Acceleration and Ordering in Solid-State Dewetting

We report on the influence of elastic strain on solid-state dewetting. Using continuum modeling, we first study the consequences of elastic stress on the pinching of the film away from the triple line during dewetting. We find that elastic stress in the solid film decreases both the time and the distance at which the film pinches in such a way that the dewetting front is accelerated. In addition, the spatial organization of islands emerging from the dewetting process is affected by strain. As an example, we demonstrate that ordered arrays of quantum dots can be achieved from solid-state dewetting of a square island in the presence of elastic stress.

Shear Banding in Entangled Polymers: Stress Plateau, Banding Location, and Lever Rule

Using molecular dynamics simulation, we study shear banding of entangled polymer melts under a steady shear. The steady shear stress vs shear rate curve exhibits a plateau spanning nearly two decades of shear rates in which shear banding is observed, and the steady shear stress remains unchanged after switching the shear rates halfway in the range of shear rates within the plateau region. In addition, we find strong correlation in the location of the shear bands between different shear rates starting from the same microstate configurations at equilibrium, which suggests the importance of the inherent structural heterogeneity in the entangled polymer network for shear banding. Furthermore, for the steady shear bands persisting to the longest simulated time of 9.0τ_d0 (disengagement time), the shear rate in the slow band and the relative proportion of the bands do not change very much with the increase of imposed shear rate, but the shear rate in the fast band increases approximately in proportion to the imposed shear rates, in contradiction to the lever rule.

Non-locality of the contact line in dynamic wetting phenomena

HYPOTHESIS

The notion of the contact line is fundamental to capillary science, where in a large category of wetting phenomena, it was always regarded as a one-dimensional object involving only microscopic length scales. This prevailing opinion had a strong impact and repercussions on the developing theories and methodologies used to interpret experimental data. It is hypothesised that this is not the case under certain conditions leading to non-local effects and requiring the development of a modified force balance at the contact line.

THEORY AND SIMULATIONS

Using the first principles of molecular dynamic simulations and a unique combination of steady state conditions and observables, the microscopic structure of the contact region and its connections with macroscopic quantities of capillary flows was revealed for the first time.

FINDINGS

The contact line is shown to become a non-local, macroscopic object involving rather complex interplay between microscopic distributions of density, velocity and friction force. It was established that the non-locality effects, which cannot be in principle captured by localised methodologies, kick off at a universal tipping point and lead to a modified force balance. The developed framework is applicable to a wide range of capillary flows to identify and analyse this regime in applications.

Disjoining-pressure-induced acceleration of mass shedding in solid-state dewetting

Surface-diffusion mediated solid-state dewetting has been observed and studied in a number of different systems during the past two decades. This process can be accompanied by the pinching of the film at a finite distance from the retracting triple line. The repetition of this pinching is often referred to as periodic mass shedding. We show that the disjoining pressure of the film can accelerate mass shedding by orders of magnitude in ultrathin films with nanometric thickness. In the presence of power-law disjoining pressures induced by van der Waals forces, the mass shedding time exhibits an approximate power-law dependence on film thickness t_{ms}∼h[over ¯]^{ν}, with ν≈6. Exponentially decaying disjoining forces also give rise to a strong acceleration of mass shedding. However, due to the finite range of the exponential potential, the mass shedding time does not exhibit a simple power-law dependence on the thickness, and is controlled by a cutoff thickness. In addition, two-dimensional simulations indicate that, within the range of thicknesses that we have studied and for isotropic dynamics, the transversal instability of a straight front does not lead to fingering, and mass shedding is the dominant instability of the dewetting front. Finally, we also show that no significant difference is observed in the dewetting dynamics between simulations based on a model with a wetting potential integrated over the film surface area, or over the projected substrate area.

Binding potentials for vapour nanobubbles on surfaces using density functional theory

We calculate density profiles of a simple model fluid in contact with a planar surface using density functional theory (DFT), in particular for the case where there is a vapour layer intruding between the wall and the bulk liquid. We apply the method of Hughes et al (2015 J. Chem. Phys. 142 074702) to calculate the density profiles for varying (specified) amounts of the vapour adsorbed at the wall. This is equivalent to varying the thickness h of the vapour at the surface. From the resulting sequence of density profiles we calculate the thermodynamic grand potential as h is varied and thereby determine the binding potential as a function of h. The binding potential obtained via this coarse-graining approach allows us to determine the disjoining pressure in the film and also to predict the shape of vapour nano-bubbles on the surface. Our microscopic DFT based approach captures information from length scales much smaller than some commonly used models in continuum mechanics.

Equilibrium Contact Angle and Adsorption Layer Properties with Surfactants

The three-phase contact line of a droplet on a smooth surface can be characterized by the Young equation. It relates the interfacial energies to the macroscopic contact angle θ_e. On the mesoscale, wettability is modeled by a film-height-dependent wetting energy f( h). Macro- and mesoscale descriptions are consistent if γ cos θ_e = γ + f( h_a), where γ and h_a are the liquid-gas interface energy and the thickness of the equilibrium liquid adsorption layer, respectively. Here, we derive a similar consistency condition for the case of a liquid covered by an insoluble surfactant. At equilibrium, the surfactant is spatially inhomogeneously distributed, implying a nontrivial dependence of θ_e on surfactant concentration. We derive macroscopic and mesoscopic descriptions of a contact line at equilibrium and show that they are consistent only if a particular dependence of the wetting energy on the surfactant concentration is imposed. This is illustrated by a simple example of dilute surfactants, for which we show excellent agreement between theory and time-dependent numerical simulations.

Dynamical Density Functional Theory for the Evaporation of Droplets of Nanoparticle Suspension

We develop a lattice gas model for the drying of droplets of a nanoparticle suspension on a planar surface, using dynamical density functional theory (DDFT) to describe the time evolution of the solvent and nanoparticle density profiles. The DDFT assumes a diffusive dynamics but does not include the advective hydrodynamics of the solvent, so the model is relevant to highly viscous or near to equilibrium systems. Nonetheless, we see an equivalent of the coffee-ring stain effect, but in the present model it occurs for thermodynamic rather the fluid-mechanical reasons. The model incorporates the effect of phase separation and vertical density variations within the droplet and the consequence of these on the nanoparticle deposition pattern on the surface. We show how to include the effect of slip or no-slip at the surface and how this is related to the receding contact angle. We also determine how the equilibrium contact angle depends on the microscopic interaction parameters.

Analytical solutions for the profile of two-dimensional droplets with finite-length precursor films

By means of the lubrication approximation we obtain the full family of static bidimensional profiles of a liquid resting on a substrate under partial-wetting conditions imposed by a disjoining-conjoining pressure. We show that for a set of quite general disjoining-conjoining pressure potentials, the free surface can adopt only five nontrivial static patterns; in particular, we find solutions when the height goes to zero which describe satisfactorily the complete free surface for a finite amount of fluid deposited on a substrate. To test the extension of the applicability of our solutions, we compare them with those obtained when the lubrication approximations are not employed and under conditions where the lubrication hypothesis are not strictly valid, and also with axisymmetric solutions. For a given disjoining-conjoining potential, we report a new analytical solution that accounts for all the five possible solutions.

Multiscale Modeling of the Three-Dimensional Meniscus Shape of a Wetting Liquid Film on Micro-/Nanostructured Surfaces

The design of structured surfaces for increasing the heat flux dissipated during boiling and evaporation processes via enhanced liquid rewetting requires prediction of the liquid meniscus shape on these surfaces. In this study, a general continuum model is developed to predict the three-dimensional meniscus shape of liquid films on micro/nanostructured surfaces based on a minimization of the system free energy that includes solid-liquid van der Waals interaction energy, surface energy, and gravitational potential. The continuum model is validated at the nanoscale against molecular dynamics simulations of water films on gold surfaces with pyramidal indentations, and against experimental measurements of water films on silicon V-groove channels at the microscale. The validated model is used to investigate the effect of film thickness and surface structure depth on the meniscus shape. The meniscus is shown to become more conformal with the surface structure as the film thickness decreases and the structure depth increases. Assuming small interface slope and small variation in film thickness, the continuum model can be linearized to obtain an explicit expression for the meniscus shape. The error of this linearized model is quantitatively assessed and shown to increase with increasing structure depth and decreasing structure pitch. The model developed can be used for accurate prediction of three-dimensional meniscus shape on structured surfaces with micro/nano-scale features, which is necessary for determining the liquid delivery rate and heat flux dissipated during thin-film evaporation. The linearized model is useful for rapid prediction of meniscus shape when the structure depth is smaller than or comparable to the liquid film thickness.

Modelling the evaporation of nanoparticle suspensions from heterogeneous surfaces

We present a Monte Carlo (MC) grid-based model for the drying of drops of a nanoparticle suspension upon a heterogeneous surface. The model consists of a generalised lattice-gas in which the interaction parameters in the Hamiltonian can be varied to model different properties of the materials involved. We show how to correctly choose the interactions, to minimise the effects of the underlying grid so that hemispherical droplets form. We also include the effects of surface roughness to examine the effects of contact-line pinning on the dynamics. When there is a 'lid' above the system, which prevents evaporation, equilibrium drops form on the surface, which we use to determine the contact angle and how it varies as the parameters of the model are changed. This enables us to relate the interaction parameters to the materials used in applications. The model has also been applied to drying on heterogeneous surfaces, in particular to the case where the suspension is deposited on a surface consisting of a pair of hydrophilic conducting metal surfaces that are either side of a band of hydrophobic insulating polymer. This situation occurs when using inkjet printing to manufacture electrical connections between the metallic parts of the surface. The process is not always without problems, since the liquid can dewet from the hydrophobic part of the surface, breaking the bridge before the drying process is complete. The MC model reproduces the observed dewetting, allowing the parameters to be varied so that the conditions for the best connection can be established. We show that if the hydrophobic portion of the surface is located at a step below the height of the neighbouring metal, the chance of dewetting of the liquid during the drying process is significantly reduced.

Nudged elastic band calculation of the binding potential for liquids at interfaces

The wetting behavior of a liquid on solid substrates is governed by the nature of the effective interaction between the liquid-gas and the solid-liquid interfaces, which is described by the binding or wetting potential g(h) which is an excess free energy per unit area that depends on the liquid film height h. Given a microscopic theory for the liquid, to determine g(h), one must calculate the free energy for liquid films of any given value of h, i.e., one needs to create and analyze out-of-equilibrium states, since at equilibrium there is a unique value of h, specified by the temperature and chemical potential of the surrounding gas. Here we introduce a Nudged Elastic Band (NEB) approach to calculate g(h) and illustrate the method by applying it in conjunction with a microscopic lattice density functional theory for the liquid. We also show that the NEB results are identical to those obtained with an established method based on using a fictitious additional potential to stabilize the non-equilibrium states. The advantages of the NEB approach are discussed.

Films, layers, and droplets: The effect of near-wall fluid structure on spreading dynamics

We present a study of the spreading of liquid droplets on a solid substrate at very small scales. We focus on the regime where effective wetting energy (binding potential) and surface tension effects significantly influence steady and spreading droplets. In particular, we focus on strong packing and layering effects in the liquid near the substrate due to underlying density oscillations in the fluid caused by attractive substrate-liquid interactions. We show that such phenomena can be described by a thin-film (or long-wave or lubrication) model including an oscillatory Derjaguin (or disjoining or conjoining) pressure and explore the effects it has on steady droplet shapes and the spreading dynamics of droplets on both an adsorption (or precursor) layer and completely dry substrates. At the molecular scale, commonly used two-term binding potentials with a single preferred minimum controlling the adsorption layer height are inadequate to capture the rich behavior caused by the near-wall layered molecular packing. The adsorption layer is often submonolayer in thickness, i.e., the dynamics along the layer consists of single-particle hopping, leading to a diffusive dynamics, rather than the collective hydrodynamic motion implicit in standard thin-film models. We therefore modify the model in such a way that for thicker films the standard hydrodynamic theory is realized, but for very thin layers a diffusion equation is recovered.

Thermal fluctuations of an interface near a contact line

The effect of thermal fluctuations near a contact line of a liquid interface partially wetting an impenetrable substrate is studied analytically and numerically. Promoting both the interface profile and the contact line position to random variables, we explore the equilibrium properties of the corresponding fluctuating contact line problem based on an interfacial Hamiltonian involving a "contact" binding potential. To facilitate an analytical treatment, we consider the case of a one-dimensional interface. The effective boundary condition at the contact line is determined by a dimensionless parameter that encodes the relative importance of thermal energy and substrate energy at the microscopic scale. We find that this parameter controls the transition from a partial wetting to a pseudopartial wetting state, the latter being characterized by a thin prewetting film of fixed thickness. In the partial wetting regime, instead, the profile typically approaches the substrate via an exponentially thinning prewetting film. We show that, independently of the physics at the microscopic scale, Young's angle is recovered sufficiently far from the substrate. The fluctuations of the interface and of the contact line give rise to an effective disjoining pressure, exponentially decreasing with height. Fluctuations therefore provide a regularization of the singular contact forces occurring in the corresponding deterministic problem.

Crystallization of soft matter under confinement at interfaces and in wedges

The surface freezing and surface melting transitions that are exhibited by a model two-dimensional soft matter system are studied. The behaviour when confined within a wedge is also considered. The system consists of particles interacting via a soft purely repulsive pair potential. Density functional theory (DFT) is used to calculate density profiles and thermodynamic quantities. The external potential due to the confining walls is modelled via a hard wall with an additional repulsive Yukawa potential. The surface phase behaviour depends on the range and strength of this repulsion: when the repulsion is weak, the wall promotes freezing at the surface of the wall. The thickness of this frozen layer grows logarithmically as the bulk liquid-solid phase coexistence is approached. Our mean-field DFT predicts that this crystalline layer at the wall must be nucleated (i.e. there is a free energy barrier) and its formation is necessarily a first-order transition, referred to as 'prefreezing', by analogy with the prewetting transition. However, in contrast to the latter, prefreezing cannot terminate in a critical point, since the phase transition involves a change in symmetry. If the wall-fluid interaction is sufficiently long ranged and the repulsion is strong enough, surface melting can occur instead. Then the interface between the wall and the bulk crystalline solid is wetted by the liquid phase as the chemical potential is decreased towards the value at liquid-solid coexistence. It is observed that the finite thickness fluid film at the wall has a broken translational symmetry due to its proximity to the bulk crystal, and so the nucleation of the wetting film can be either first order or continuous. Our mean-field theory predicts that for certain wall potentials there is a premelting critical point analogous to the surface critical point for the prewetting transition. When the fluid is confined within a linear wedge, this can strongly promote freezing when the opening angle of the wedge is commensurate with the crystal lattice.

Wetting and evaporation of salt-water nanodroplets: A molecular dynamics investigation

We employ molecular dynamics simulations to study the wetting and evaporation of salt-water nanodroplets on platinum surfaces. Our results show that the contact angle of the droplets increases with the salt concentration. To verify this, a second simulation system of a thin salt-water film on a platinum surface is used to calculate the various surface tensions. We find that both the solid-liquid and liquid-vapor surface tensions increase with salt concentration and as a result these cause an increase in the contact angle. However, the evaporation rate of salt-water droplets decreases as the salt concentration increases, due to the hydration of salt ions. When the water molecules have all evaporated from the droplet, two forms of salt crystals are deposited, clump and ringlike, depending on the solid-liquid interaction strength and the evaporation rate. To form salt crystals in a ring, it is crucial that there is a pinned stage in the evaporation process, during which salt ions can move from the center to the rim of the droplets. With a stronger solid-liquid interaction strength, a slower evaporation rate, and a higher salt concentration, a complete salt crystal ring can be deposited on the surface.

Directed transport of polymer drops on vibrating superhydrophobic substrates: a molecular dynamics study

Using molecular dynamics simulations of a coarse-grained polymer liquid we investigate the transport of droplets on asymmetrically structured (saw-tooth shaped), vibrating substrates. Due to a continuous supply of power by substrate vibrations and the asymmetry of its topography, the droplets are driven in a preferred direction. We study this directed motion as a function of the size of the droplets, the linear dimensions of the substrate corrugation, and the period of vibrations. Two mechanisms of driven transport are identified: (i) one that relies on the droplet's contact lines and (ii) in a range of vibration periods, the entire contact area contributes to the driving. In this latter regime, the set-up may be used in experiments for sorting droplets according to their size. Additionally, we show that the linear dimension of the substrate corrugation affects the flux inside the droplet. While on a substrate with a fine corrugation droplets mostly slide, on a more coarsely corrugated substrate the flux may exhibit an additional rotation pattern.

Force-Field Parameters from the SAFT-γ Equation of State for Use in Coarse-Grained Molecular Simulations

A description of fluid systems with molecular-based algebraic equations of state (EoSs) and by direct molecular simulation is common practice in chemical engineering and the physical sciences, but the two approaches are rarely closely coupled. The key for an integrated representation is through a well-defined force field and Hamiltonian at the molecular level. In developing coarse-grained intermolecular potential functions for the fluid state, one typically starts with a detailed, bottom-up quantum-mechanical or atomic-level description and then integrates out the unwanted degrees of freedom using a variety of techniques; an iterative heuristic simulation procedure is then used to refine the parameters of the model. By contrast, with a top-down technique, one can use an accurate EoS to link the macroscopic properties of the fluid and the force-field parameters. We discuss the latest developments in a top-down representation of fluids, with a particular focus on a group-contribution formulation of the statistical associating fluid theory (SAFT-γ). The accurate SAFT-γ EoS is used to estimate the parameters of the Mie force field, which can then be used with confidence in direct molecular simulations to obtain thermodynamic, structural, interfacial, and dynamical properties that are otherwise inaccessible from the EoS. This is exemplified for several prototypical fluids and mixtures, including carbon dioxide, hydrocarbons, perfluorohydrocarbons, and aqueous surfactants.

Disjoining pressure and the film-height-dependent surface tension of thin liquid films: new insight from capillary wave fluctuations

In this paper we review simulation and experimental studies of thermal capillary wave fluctuations as an ideal means for probing the underlying disjoining pressure and surface tensions, and more generally, fine details of the Interfacial Hamiltonian Model. We discuss recent simulation results that reveal a film-height-dependent surface tension not accounted for in the classical Interfacial Hamiltonian Model. We show how this observation may be explained bottom-up from sound principles of statistical thermodynamics and discuss some of its implications.

Rupture mechanism of liquid crystal thin films realized by large-scale molecular simulations

The ability of liquid crystal (LC) molecules to respond to changes in their environment makes them an interesting candidate for thin film applications, particularly in bio-sensing, bio-mimicking devices, and optics. Yet the understanding of the (in)stability of this family of thin films has been limited by the inherent challenges encountered by experiment and continuum models. Using unprecedented large-scale molecular dynamics (MD) simulations, we address the rupture origin of LC thin films wetting a solid substrate at length scales similar to those in experiment. Our simulations show the key signatures of spinodal instability in isotropic and nematic films on top of thermal nucleation, and importantly, for the first time, evidence of a common rupture mechanism independent of initial thickness and LC orientational ordering. We further demonstrate that the primary driving force for rupture is closely related to the tendency of the LC mesogens to recover their local environment in the bulk state. Our study not only provides new insights into the rupture mechanism of liquid crystal films, but also sets the stage for future investigations of thin film systems using peta-scale molecular dynamics simulations.

Coexistence of spinodal instability and thermal nucleation in thin-film rupture: insights from molecular levels

Despite extensive investigation using hydrodynamic models and experiments over the past decades, there remain open questions regarding the origin of the initial rupture of thin liquid films. One of the reasons that makes it difficult to identify the rupture origin is the coexistence of two dewetting mechanisms, namely, thermal nucleation and spinodal instability, as observed in many experimental studies. Using a coarse-grained model and large-scale molecular dynamics simulations, we are able to characterize the very early stage of dewetting in nanometer-thick liquid-metal films wetting a solid substrate. We observe the features characteristic of both spinodal instability and thermal nucleation in the spontaneously dewetting films and show that these two macroscopic mechanisms share a common origin at molecular levels.

Stability of thin liquid films and sessile droplets under confinement

The stability of nonvolatile thin liquid films and of sessile droplets is strongly affected by finite size effects. We analyze their stability within the framework of density functional theory using the sharp kink approximation, i.e., on the basis of an effective interface Hamiltonian. We show that finite size effects suppress spinodal dewetting of films because it is driven by a long-wavelength instability. Therefore nonvolatile films are stable if the substrate area is too small. Similarly, nonvolatile droplets connected to a wetting film become unstable if the substrate area is too large. This instability of a nonvolatile sessile droplet turns out to be equivalent to the instability of a volatile drop which can attain chemical equilibrium with its vapor.