Relaxation of surface tension in the free-surface boundary layer of simple Lennard-Jones liquids

Dynamic surface stress field of the pure liquid-vapor interface subjected to the cyclic loads

We demonstrate a methodology for computationally investigating the mechanical response of a pure molten lead surface system to the lateral mechanical cyclic loads and try to answer the following question: how does the dynamically driven liquid surface system follow the classical physics of the elastic-driven oscillation? The steady-state oscillation of the dynamic surface tension (or excess stress) under cyclic load, including the excitation of high-frequency vibration mode at different driving frequencies and amplitudes, was compared with the classical theory of a single-body driven damped oscillator. Under the highest studied frequency (50 GHz) and amplitude (5%) of the load, the increase of in (mean value) dynamic surface tension could reach ∼5%. The peak and trough values of the instantaneous dynamic surface tension could reach (up to) 40% increase and (up to) 20% decrease compared to the equilibrium surface tension, respectively. The extracted generalized natural frequencies seem to be intimately related to the intrinsic timescales of the atomic temporal-spatial correlation functions of the liquids both in the bulk region and in the outermost surface layers. These insights uncovered could be helpful for quantitative manipulation of the liquid surface using ultrafast shockwaves or laser pulses.

Heat flow through a liquid-vapor interface in a nano-channel: the effect of end-grafting polymers on a wall

Heat transfer through a liquid-vapor interface is a complex phenomenon and crucially relevant in heat-removal and cryogenic applications. The physical coupling among confining walls, liquid and vapor phases is very important for controlling or improving cooling rates or condensation efficiency. Surface modification is a promising route, which has been explored to taylor the heat transfer through confined two-phase systems. We use coarse-grained molecular-dynamics simulations to study the heat transfer through a nano-confined liquid-vapor interface as a function of fluid filling. We set up a stationary heat flow through a liquid-vapor interface, stabilized with the liquid in contact with a colder wall and a vapor in contact with a hotter wall. For these physical conditions, we perform extensive simulations by progressively increasing the number of fluid particles, i.e. the channel filling, and measure the fluid distribution in the channel, density, pressure and temperature profiles We also compare the heat flux behavior between a bare-surfaces nano-channel and others where the hot surface was coated with end-grafted polymers, with different wetting affinities and bending properties. We take extreme cases of polymer properties to obtain a general picture of the polymer effect on the heat transfer, as compared with the bare surfaces. We find that walls covered by end-grafted solvophylic polymers change the heat flux by a factor of 6, as compared with bare walls, if the liquid phase is in contact with the polymers. Once the liquid wets the coated wall, the improve on heat flux is smaller and dominated by the grafting density. We also find that for a wall coated with stiff polymers, the jump in heat flux takes place at a significantly lower channel filling, when the polymers' free ends interact with the liquid surface. Interestingly, the morphology of the polymers induces a 'liquid bridge' between the liquid phase and the hot wall, through which heat is transported with high (liquid-like) thermal conductivity.

Entropic surface segregation from athermal polymer blends: Polymer flexibility vs bulkiness

We examine athermal binary blends composed of conformationally asymmetric polymers of equal molecular volume next to a surface of width ξ. The self-consistent field theory (SCFT) of Gaussian chains predicts that the more compact polymer with the shorter average end-to-end length, R0, is entropically favored at the surface. Here, we extend the SCFT to worm-like chains with small persistence lengths, ℓ p, relative to their contour lengths, ℓ c, for which [Formula: see text]. In the limit of ℓ p ≪ ξ, we recover the Gaussian-chain prediction where the segregation depends only on the product ℓ p ℓ c, but for realistic polymer/air surfaces with ξ ∼ ℓ p, the segregation depends separately on the two quantities. Although the surface continues to favor flexible polymers with smaller ℓ p and bulky polymers with shorter ℓ c, the effect of bulkiness is more pronounced. This imbalance can, under specific conditions, lead to anomalous surface segregation of the more extended polymer. For this to happen, the polymer must be bulkier and stiffer, with a stiffness that is sufficient to produce a larger R0 yet not so rigid as to reverse the surface affinity that favors bulky polymers.

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.

Effect of chain stiffness on the entropic segregation of chain ends to the surface of a polymer melt

Entropic segregation of chain ends to the surface of a monodisperse polymer melt and its effect on surface tension are examined using self-consistent field theory (SCFT). In order to assess the dependence on chain stiffness, the SCFT is solved for worm-like chains. Our focus is still on relatively flexible polymers, where the persistence length of the polymer, ℓ _p , is comparable to the width of the surface profile, ξ, but still much smaller than the total contour length of the polymer, ℓ _c . Even this small degree of rigidity causes a substantial increase in the level of segregation, relative to that of totally flexible Gaussian chains. Nevertheless, the long-range depletion that balances the surface excess still exhibits the same universal shape derived for Gaussian chains. Furthermore, the excess continues to reduce the surface tension by one unit of k _B T per chain end, which results in the usual N ^-1 reduction in surface tension observed by experiments. This enhanced segregation will also extend to polydisperse melts, causing the molecular-weight distribution at the surface to shift towards smaller N _n relative to the bulk. This provides a partial explanation for recent quantitative differences between experiments and SCFT calculations for flexible polymers.

Unexpected finite size effects in interfacial systems: Why bigger is not always better-Increase in uncertainty of surface tension with bulk phase width

We present an unexpected finite size effect affecting interfacial molecular simulations that is proportional to the width-to-surface-area ratio of the bulk phase L_l/A. This finite size effect has a significant impact on the variance of surface tension values calculated using the virial summation method. A theoretical derivation of the origin of the effect is proposed, giving a new insight into the importance of optimising system dimensions in interfacial simulations. We demonstrate the consequences of this finite size effect via a new way to estimate the surface energetic and entropic properties of simulated air-liquid interfaces. Our method is based on macroscopic thermodynamic theory and involves comparing the internal energies of systems with varying dimensions. We present the testing of these methods using simulations of the TIP4P/2005 water forcefield and a Lennard-Jones fluid model of argon. Finally, we provide suggestions of additional situations, in which this finite size effect is expected to be significant, as well as possible ways to avoid its impact.

Hydrodynamics of Moving Contact Lines: Macroscopic versus Microscopic

The fluid-mechanics community is currently divided in assessing the boundaries of applicability of the macroscopic approach to fluid mechanical problems. Can the dynamics of nanodroplets be described by the same macroscopic equations as are used for macrodroplets? To the greatest degree, this question should be addressed to the moving-contact-line problem. The problem is naturally multiscale, where even using slip boundary conditions results in spurious numerical solutions and transcendental stagnation regions in modeling in the vicinity of the contact line. In this article, it is demonstrated through mutual comparisons between macroscopic modeling and molecular dynamics simulations that a small, albeit natural, change in the boundary conditions is all that is necessary to completely regularize the problem and eliminate these nonphysical effects. The limits of the macroscopic approach applied to the moving-contact-line problem have been tested and validated on the basis of microscopic first-principles molecular dynamics simulations.

Young's Equation for a Two-Liquid System on the Nanometer Scale

We use large-scale molecular dynamics simulations to study the Lennard-Jones forces acting at the various interfaces of a liquid bridge (liquid 1) between two realistic solid plates on the scale of few nanometers when the two free surfaces are in contact with a second immiscible liquid (liquid 2) with an interfacial tension of γ₁₂. Each plate comprises a regular square planar lattice of atoms arranged in three atomic layers. To maintain rigidity while allowing momentum exchange with the liquid, solid atoms are allowed to vibrate thermally around their initial positions by a strong harmonic potential. By varying the solid-liquid coupling, we investigate a range of nonzero contact angles between the liquid-liquid interface and the solid. We first compute the forces when the plates are stationary (equilibrium case), from the perspectives of both the liquid and the solid. Our results confirm that the normal and tangential components of the computed interfacial forces at each contact line are consistent with Young's equation on this small scale. In particular, we show that the tangential force exerted by the liquid-liquid interface on the plates is given by the difference in the individual works of adhesion of the two liquids and equal to γ₁₂ cos θ_1,2⁰, where θ_1,2⁰ is the equilibrium contact angle measured through liquid 1. This result, which differs from that expected for a single liquid, is relevant to the interactions and behavior of two liquid-solid systems in nanotechnology. We then study the forces when the plates are translated at equal speeds in opposite directions over a range of steady velocities (dynamic case) and repeat the measurements of the force exerted by the liquid-liquid interface on the solid. We find that the normal and tangential components of this force are still correctly predicted by the normal and tangential components of the interfacial tension, provided only that the equilibrium contact angle is replaced by its dynamic analogue θ_1,2^D. Usually assumed without proof, this result is significant for our proper understanding of dynamic wetting at all scales.

Dynamic Contact Angle at the Nanoscale: A Unified View

Generation of a dynamic contact angle in the course of wetting is a fundamental phenomenon of nature. Dynamic wetting processes have a direct impact on flows at the nanoscale, and therefore, understanding them is exceptionally important to emerging technologies. Here, we reveal the microscopic mechanism of dynamic contact angle generation. It has been demonstrated using large-scale molecular dynamics simulations of bead-spring model fluids that the main cause of local contact angle variations is the distribution of microscopic force acting at the contact line region. We were able to retrieve this elusive force with high accuracy. It has been directly established that the force distribution can be solely predicted on the basis of a general friction law for liquid flow at solid surfaces by Thompson and Troian. The relationship with the friction law provides both an explanation of the phenomenon of dynamic contact angle and a methodology for future predictions. The mechanism is intrinsically microscopic, universal, and irreducible and is applicable to a wide range of problems associated with wetting phenomena.

Crystal growth from a supersaturated melt: relaxation of the solid-liquid dynamic stiffness

We discuss the growth process of a crystalline phase out of a metastable over-compressed liquid that is brought into contact with a crystalline substrate. The process is modeled by means of molecular dynamics. The particles interact via the Lennard-Jones potential and their motion is locally thermalized by Langevin dynamics. We characterize the relaxation process of the solid-liquid interface, showing that the growth speed is maximal for liquid densities above the solid coexistence density, and that the structural properties of the interface rapidly converge to equilibrium-like properties. In particular, we show that the off-equilibrium dynamic stiffness can be extracted using capillary wave theory arguments, even if the growth front moves fast compared to the typical diffusion time of the compressed liquid, and that the dynamic stiffness converges to the equilibrium stiffness in times much shorter than the diffusion time.

Nonequilibrium interfaces in colloidal fluids

The time-dependent structure, interfacial tension, and evaporation of an oversaturated colloid-rich (liquid) phase in contact with an undersaturated colloid-poor (vapor) phase of a colloidal dispersion is investigated theoretically during the early-stage relaxation, where the interface is relaxing towards a local equilibrium state while the bulk phases are still out of equilibrium. Since systems of this type exhibit a clear separation of colloidal and solvent relaxation time scales with typical times of interfacial tension measurements in between, they can be expected to be suitable for analogous experimental studies, too. The major finding is that, irrespective of how much the bulk phases differ from two-phase coexistence, the interfacial structure and the interfacial tension approach those at two-phase coexistence during the early-stage relaxation process. This is a surprising observation since it implies that the relaxation towards global equilibrium of the interface is not following but preceding that of the bulk phases. Scaling forms for the local chemical potential, the flux, and the dissipation rate exhibit qualitatively different leading order contributions depending on whether an equilibrium or a nonequilibrium system is considered. The degree of nonquilibrium between the bulk phases is found to not influence the qualitative relaxation behavior (i.e., the values of power-law exponents), but to determine the quantitative deviation of the observed quantities from their values at two-phase coexistence. Whereas the underlying dynamics differs between colloidal and molecular fluids, the behavior of quantities such as the interfacial tension approaching the equilibrium values during the early-stage relaxation process, during which nonequilibrium conditions of the bulk phases are not changed, can be expected to occur for both types of systems.

Relaxation of surface tension in the liquid-solid interfaces of Lennard-Jones liquids

We have established the surface tension relaxation time in the liquid-solid interfaces of Lennard-Jones (LJ) liquids by means of direct measurements in molecular dynamics (MD) simulations. The main result is that the relaxation time is found to be almost independent of the molecular structures and viscosity of the liquids (at 70-fold change) used in our study and lies in such a range that in slow hydrodynamic motion the interfaces are expected to be at equilibrium. The implications of our results for the modeling of dynamic wetting processes and interpretation of dynamic contact angle data are discussed.