Capillary instabilities by fluctuation induced forces

Critical Thickness of Free-Standing Nanothin Films Made of Melted Polyethylene Chains via Molecular Dynamics

The mechanical stability of nanothin free-standing films made of melted polyethylene chains was predicted via molecular dynamics simulations in the range of 373.15-673.15 K. The predicted critical thickness, tc, increased with the square of the temperature, T, with additional chains needed as T increased. From T = 373.15 K up to the thermal limit of stability for polyethylene, tc values were in the range of nanothin thicknesses (3.42-5.63 nm), which approximately corresponds to 44-55 chains per 100 nm2. The density at the center of the layer and the interfacial properties studied (density profiles, interfacial thickness, and radius of gyration) showed independence from the film thickness at the same T. The polyethylene layer at its tc showed a lower melting T (<373.15 K) than bulk polyethylene.

Field-induced formation and growth of pillars on films of bisphenol-A-polycarbonate

AFM image and FFT of the surface topologies of the BPAPC films of thickness 60 nm under the action of an electric voltage of 30 V at 200 °C for 70 min.

Kinetics of Field-Induced Surface Patterns on PMMA

A simple model was developed to analyze the growth of a liquid pillar under the action of an electric field between two parallel electrodes. A quadratic relationship between time and the diameter of the pillar was obtained. The diameter of the pillar increases with time. Large electric field assists the growth of the liquid pillar, while a liquid with a large viscosity hinders the growth of the liquid pillar. The field-induced formation and growth of PMMA pillars on PMMA films were observed using the configuration of a parallel capacitor. Pillars of larger sizes and smaller densities were formed on thicker PMMA films than on thinner PMMA films. The root-mean-square ( https://en.wikipedia.org/wiki/Root_mean_square ) diameter of the pillars increases with the increase of the annealing time and annealing temperature. The growth behavior of the pillars can be described by an Arrhenius relation with an activation energy of 24.4 kJ/mol, suggesting that the growth of the pillars is controlled by a thermal activation process.

Hydrodynamic fluctuation-induced forces in confined fluids

We study thermal, fluctuation-induced hydrodynamic interaction forces in a classical, compressible, viscous fluid confined between two rigid, planar walls with no-slip boundary conditions. We calculate hydrodynamic fluctuations using the linearized, stochastic Navier-Stokes formalism of Landau and Lifshitz. The mean fluctuation-induced force acting on the fluid boundaries vanishes in this system, so we evaluate the two-point, time-dependent force correlations. The equal-time correlation function of the forces acting on a single wall gives the force variance, which we show to be finite and independent of the plate separation at large inter-plate distances. The equal-time, cross-plate force correlation, on the other hand, decays with the inverse inter-plate distance and is independent of the fluid viscosity at large distances; it turns out to be negative over the whole range of plate separations, indicating that the two bounding plates are subjected to counter-phase correlations. We show that the time-dependent force correlations exhibit damped temporal oscillations for small plate separations and a more irregular oscillatory behavior at large separations. The long-range hydrodynamic correlations reported here represent a "secondary Casimir effect", because the mean fluctuation-induced force, which represents the primary Casimir effect, is absent.

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.

Accelerating dewetting on deformable substrates by adding a liquid underlayer

We investigated the dependence of the dewetting velocity of a thin, low-viscosity polystyrene (PS) top film on a poly(methyl methacrylate) (PMMA) double layer consisting of a low-viscosity underlayer of thickness h(L) coated with a high-viscosity middle layer of thickness h(M). The addition of the liquid underlayer generated complex nonmonotonic behavior of the dewetting velocity as a function of increasing h(M). In particular, we observed an acceleration of dewetting for an intermediate range of h(M). This phenomenon has been interpreted by a combination deformation of the middle elastic layer and a concurrent change in the contact angle. On one hand, deformation led to the formation of a trench that dissipated energy during its movement through the liquid underlayer and thus caused a slowing down of dewetting. However, with an increase in the thickness of the elastic middle layer, the size of the trench decreased and its influence on the dewetting velocity also decreased. On the other hand, the deformation of the elastic layer also led to an increase in the contact angle. This increase in the driving capillary forces caused an increase in the dewetting velocity.

Dewetting induced by complete versus nonretarded van der Waals forces

Spontaneous rupture of some polymer films upon heating is commonplace. The very criterion for this instability is the system free energy, G(L), possessing a negative curvature. In films that are apolar with h < or = 100 nm van der Waals (vdW) interactions usually constitute a major contribution to G(L) for which the approximate form G(L) = -A/12piL(2) (where A is the Hamaker constant), ignoring retardation, has been widely used. In this work, we investigate the limits to this approximation by calculating the complete vdW interactions for popular polymer film systems in dewetting experiments including air-polystyrene-SiO2-Si, air-polystyrene-poly(methyl methacrylate)-Si, and air-poly(methyl methacrylate)-polystyrene-Si based on the theory of Dzyaloshinskii, Lifshitz, and Pitaevskii (DLP). We found that retardation effects could produce significant modifications to G(L) even when the thickness of the polymer and/or the interlayer is only 1-2 nm, contrary to conventional presumption.

Morphology changes in the evolution of liquid two-layer films

We consider a thin film consisting of two layers of immiscible liquids on a solid horizontal (heated) substrate. Both the free liquid-liquid and the liquid-gas interface of such a bilayer liquid film may be unstable due to effective molecular interactions relevant for ultrathin layers below 100-nm thickness, or due to temperature-gradient-caused Marangoni flows in the heated case. Using a long-wave approximation, we derive coupled evolution equations for the interface profiles for the general nonisothermal situation allowing for slip at the substrate. Linear and nonlinear analyses of the short- and long-time film evolution are performed for isothermal ultrathin layers, taking into account destabilizing long-range and stabilizing short-range molecular interactions. It is shown that the initial instability can be of a varicose, zigzag, or mixed type. However, in the nonlinear stage of the evolution the mode type, and therefore the pattern morphology, can change via switching between two different branches of stationary solutions or via coarsening along a single branch.

Alternative pathways of dewetting for a thin liquid two-layer film

We consider two stacked ultrathin layers of different liquids on a solid substrate. Using long-wave theory, we derive coupled evolution equations for the free liquid-liquid and liquid-gas interfaces. Depending on the long-range van der Waals forces and the ratio of the layer thicknesses, the system follows different pathways of dewetting. The instability may be driven by varicose or zigzag modes and leads to film rupture either at the liquid-gas interface or at the substrate. We predict that the faster layer drives the evolution and may accelerate the rupture of the slower layer by orders of magnitude, thereby promoting the rupture of rather thick films.

Molecular forces caused by the confinement of thermal noise

We have investigated spontaneous surface instabilities of very thin polymer films. Film stability and the wavelength of the dominating unstable mode were found to depend sensitively on the media adjacent to the film. Our experimental results cannot be explained by van der Waals interactions alone. To account for the presence of an additional destabilizing force, we propose that the geometrical confinement of thermally excited acoustic waves gives rise to a force that is strong enough to destabilize thin films. This thermoacoustic effect is of similar magnitude as van der Waals forces.