Confinement enhanced viscosity vs shear thinning in lubricated ice friction.

The ice surface is known for presenting a very small kinetic friction coefficient, but the origin of this property remains highly controversial to date. In this work, we revisit recent computer simulations of ice sliding on atomically smooth substrates, using newly calculated bulk viscosities for the TIP4P/ice water model. The results show that spontaneously formed premelting films in static conditions exhibit an effective viscosity that is about twice the bulk viscosity. However, upon approaching sliding speeds in the order of m/s, the shear rate becomes very large, and the viscosities decrease by several orders of magnitude. This shows that premelting films can act as an efficient lubrication layer despite their small thickness and illustrates an interesting interplay between confinement enhanced viscosities and shear thinning. Our results suggest that the strongly thinned viscosities that operate under the high speed skating regime could largely reduce the amount of frictional heating.

Surface tension of bulky colloids, capillarity under gravity, and the microscopic origin of the Kardar-Parisi-Zhang equation.

Experimental measurements of the surface tension of colloidal interfaces have long been in conflict with computer simulations. In this Letter we show that the surface tension of colloids as measured by surface fluctuations picks up a gravity-dependent contribution which removes the discrepancy. The presence of this term puts a strong constraint on the structure of the interface which allows one to identify corrections to the fundamental equation of equilibrium capillarity and deduce bottom up the microscopic origin of a growth model with close relation to the Kardar-Parisi-Zhang equation.

Self-diffusion and shear viscosity for the TIP4P/Ice water model.

With an ever-increasing interest in water properties, many intermolecular force fields have been proposed to describe the behavior of water. Unfortunately, good models for liquid water usually cannot provide simultaneously an accurate melting point for ice. For this reason, the TIP4P/Ice model was developed for targeting the melting point and has become the preferred choice for simulating ice at coexistence. Unfortunately, available data for its dynamic properties in the liquid state are scarce. Therefore, we demonstrate a series of simulations aimed at the calculation of transport coefficients for the TIP4P/Ice model over a large range of thermodynamic conditions, ranging from T = 245 K to T = 350 K, for the temperature, and from p = 0 to p = 500 MPa, for the pressure. We have found that the self-diffusion (shear viscosity) exhibits smaller (increased) values than TIP4P/2005 and experiments. However, rescaling the temperature with respect to the triple point temperature, as in a corresponding states plot, we find that TIP4P/Ice compares very well with TIP4P/2005 and experiment. Such observations allow us to infer that despite the different original purposes of these two models examined here, one can benefit from a vast number of reports regarding the behavior of transport coefficients for the TIP4P/2005 model and utilize them following the routine described in this paper.

Ice friction at the nanoscale.

The origin of ice slipperiness has been a matter of great controversy for more than a century, but an atomistic understanding of ice friction is still lacking. Here, we perform computer simulations of an atomically smooth substrate sliding on ice. In a large temperature range between 230 and 266 K, hydrophobic sliders exhibit a premelting layer similar to that found at the ice/air interface. On the contrary, hydrophilic sliders show larger premelting and a strong increase of the first adsorption layer. The nonequilibrium simulations show that premelting films of barely one-nanometer thickness are sufficient to provide a lubricating quasi-liquid layer with rheological properties similar to bulk undercooled water. Upon shearing, the films display a pattern consistent with lubricating Couette flow, but the boundary conditions at the wall vary strongly with the substrate's interactions. Hydrophobic walls exhibit large slip, while hydrophilic walls obey stick boundary conditions with small negative slip. By compressing ice above atmospheric pressure, the lubricating layer grows continuously, and the rheological properties approach bulk-like behavior. Below 260 K, the equilibrium premelting films decrease significantly. However, a very large slip persists on the hydrophobic walls, while the increased friction on hydrophilic walls is sufficient to melt ice and create a lubrication layer in a few nanoseconds. Our results show that the atomic-scale frictional behavior of ice is a combination of spontaneous premelting, pressure melting, and frictional heating.

Intermolecular forces at ice and water interfaces: Premelting, surface freezing, and regelation.

Using Lifshitz theory, we assess the role of van der Waals forces at interfaces of ice and water. The results are combined with measured structural forces from computer simulations to develop a quantitative model of the surface free energy of premelting films. This input is employed within the framework of wetting theory and allows us to predict qualitatively the behavior of quasi-liquid layer thickness as a function of ambient conditions. Our results emphasize the significance of vapor pressure. The ice-vapor interface is shown to exhibit only incomplete premelting, but the situation can shift to a state of complete surface melting above water saturation. The results obtained serve also to assess the role of subsurface freezing at the water-vapor interface, and we show that intermolecular forces favor subsurface ice nucleation only in conditions of water undersaturation. We show that ice regelation at ambient pressure may be explained as a process of capillary freezing, without the need to invoke the action of bulk pressure melting. Our results for van der Waals forces are exploited in order to gauge dispersion interactions in empirical point charge models of water.

Analytical theory for the crossover from retarded to non-retarded interactions between metal plates.

The van der Waals force established between two surfaces plays a central role in many phenomena, such as adhesion or friction. However, the dependence of this forces on the distance of separation between plates is very complex. Two widely different non-retarded and retarded regimes are well known, but these have been traditionally studied separately. Much less is known about the important experimentally accessible cross-over regime. In this study, we provide analytical approximations for the van der Waals forces between two plates that interpolates exactly between the short distance and long distance behavior, and provides new insight into the crossover from London to Casimir forces at finite temperature. At short distance, where the behavior is dominated by non-retarded interactions, we work out a very accurate simplified approximation for the Hamaker constant which adopts analytical form for both the Drude and Lorentz models of dielectric response. We apply our analytical expressions for the study of forces between metallic plates, and observe very good agreement with exact results from numerical calculations. Our results show that contributions of interband transitions remain important in the experimentally accessible regime of decades nm for several metals, including gold.

Lifshitz theory of wetting films at three phase coexistence: The case of ice nucleation on Silver Iodide (AgI).

HYPOTHESIS: As a fluid approaches three phase coexistence, adsorption may take place by the successive formation of two intervening wetting films. The equilibrium thickness of these wetting layers is the result of a delicate balance of intermolecular forces, as dictated by an underlying interface potential. The van der Waals forces for the two variable adsorption layers may be formulated exactly from Dzyaloshinskii-Lifshitz-Pitaevskii theory, and analytical approximations may be derived that extent well beyond the validity of conventional Hamaker theory. CALCULATIONS: We consider the adsorption equilibrium of water vapor on Silver Iodide where both ice and a water layers can form simultaneously and compete for the vapor as the triple point is approached. We perform numerical calculations of Lifshitz theory for this complex system and work out analytical approximations which provide quantitative agreement with the numerical results. FINDINGS: At the three phase contact line between AgI/water/air, surface forces promote growth of ice both on the AgI/air and the water/vapor interfaces, lending support to a contact nucleation mode of AgI in the atmosphere. Our approach provides a framework for the description of adsorption at three phase coexistence, and allows for the study of ice nucleation efficiency on atmospheric aerosols.

How ice grows from premelting films and water droplets.

Close to the triple point, the surface of ice is covered by a thin liquid layer (so-called quasi-liquid layer) which crucially impacts growth and melting rates. Experimental probes cannot observe the growth processes below this layer, and classical models of growth by vapor deposition do not account for the formation of premelting films. Here, we develop a mesoscopic model of liquid-film mediated ice growth, and identify the various resulting growth regimes. At low saturation, freezing proceeds by terrace spreading, but the motion of the buried solid is conveyed through the liquid to the outer liquid-vapor interface. At higher saturations water droplets condense, a large crater forms below, and freezing proceeds undetectably beneath the droplet. Our approach is a general framework that naturally models freezing close to three phase coexistence and provides a first principle theory of ice growth and melting which may prove useful in the geosciences.

Surface phase transitions and crystal habits of ice in the atmosphere.

With climate modeling predicting a raise of at least 2°C by year 2100, the fate of ice has become a serious concern, but we still do not understand how ice grows (or melts). In the atmosphere, crystal growth rates of basal and prism facets exhibit an enigmatic temperature dependence and crossover up to three times in a range between 0° and -40°. Here, we use large-scale computer simulations to characterize the ice surface and identify a sequence of previously unidentified phase transitions on the main facets of ice crystallites. Unexpectedly, we find that as temperature is increased, the crystal surface transforms from a disordered phase with proliferation of steps to a smooth phase with small step density. This causes the anomalous increase of step free energies and provides the long sought explanation for the enigmatic crossover of snow crystal growth rates found in the atmosphere.

Micelle-directed chiral seeded growth on anisotropic gold nanocrystals.

Surfactant-assisted seeded growth of metal nanoparticles (NPs) can be engineered to produce anisotropic gold nanocrystals with high chiroptical activity through the templating effect of chiral micelles formed in the presence of dissymmetric cosurfactants. Mixed micelles adsorb on gold nanorods, forming quasihelical patterns that direct seeded growth into NPs with pronounced morphological and optical handedness. Sharp chiral wrinkles lead to chiral plasmon modes with high dissymmetry factors (~0.20). Through variation of the dimensions of chiral wrinkles, the chiroptical properties can be tuned within the visible and near-infrared electromagnetic spectrum. The micelle-directed mechanism allows extension to other systems, such as the seeded growth of chiral platinum shells on gold nanorods. This approach provides a reproducible, simple, and scalable method toward the fabrication of NPs with high chiral optical activity.

Rounded Layering Transitions on the Surface of Ice.

Understanding the wetting properties of premelting films requires knowledge of the film's equation of state, which is not usually available. Here we calculate the disjoining pressure curve of premelting films and perform a detailed thermodynamic characterization of premelting behavior on ice. Analysis of the density profiles reveals the signature of weak layering phenomena, from one to two and from two to three water molecular layers. However, disjoining pressure curves, which closely follow expectations from a renormalized mean field liquid state theory, show that there are no layering phase transitions in the thermodynamic sense along the sublimation line. Instead, we find that transitions at mean field level are rounded due to capillary wave fluctuations. We see signatures that true first order layering transitions could arise at low temperatures, for pressures between the metastable line of water-vapor coexistence and the sublimation line. The extrapolation of the disjoining pressure curve above water-vapor saturation displays a true first order phase transition from a thin to a thick film consistent with experimental observations.

Structure and water attachment rates of ice in the atmosphere: role of nitrogen.

In this work we perform computer simulations of the ice surface in order to elucidate the role of nitrogen in the crystal growth rates and crystal habits of snow in the atmosphere. In pure water vapor at temperatures typical of ice crystal formation in cirrus clouds, we find that basal and primary prismatic facets exhibit a layer of premelted ice, with thickness in the subnanometer range. For partial pressures of 1 bar, well above the expected values in the troposphere, we find that only small amounts of nitrogen are adsorbed. The adsorption takes place onto the premelted surface, and hardly any nitrogen dissolves within the premelting film. The premelting film thickness does not change either. We quantify the resulting change of the ice/vapor surface tension to be in the hundredth of mN m-1 and find that the structure of the pristine ice surface is not changed in a significant manner. We perform a trajectory analysis of colliding water molecules, and find that the attachment rates from direct ballistic collision are very close to unity irrespective of the nitrogen pressure. Nitrogen is however at sufficient density to deflect a fraction of trajectories with smaller distance than the mean free path. Our results show explicitly that the reported differences in growth rates measured in pure water vapor and a controlled nitrogen atmosphere are not related to a significant disruption of the ice surface due to nitrogen adsorption. On the contrary, we show clearly from our trajectory analysis that nitrogen slows down the crystal growth rates due to collisions between water molecules with bulk nitrogen gas. This clarifies the long standing controversy of the role of inert gases on crystal growth rates and demonstrates their influence is solely related to the diffusion limited flow of water vapor across the gas phase.

Disconnecting Symmetry Breaking from Seeded Growth for the Reproducible Synthesis of High Quality Gold Nanorods.

One of the major difficulties hindering the widespread application of colloidal anisotropic plasmonic nanoparticles is the limited robustness and reproducibility of multistep synthetic methods. We demonstrate herein that the reproducibility and reliability of colloidal gold nanorod (AuNR) synthesis can be greatly improved by disconnecting the symmetry-breaking event from the seeded growth process. We have used a modified silver-assisted seeded growth method in the presence of the surfactant hexadecyltrimethylammonium bromide and n-decanol as a co-surfactant to prepare small AuNRs in high yield, which were then used as seeds for the growth of high quality AuNR colloids. Whereas the use of n-decanol provides a more-rigid micellar system, the growth on anisotropic seeds avoids sources of irreproducibility during the symmetry breaking step, yielding uniform AuNR colloids with narrow plasmon bands, ranging from 600 to 1270 nm, and allowing the fine-tuning of the final dimensions. This method provides a robust route for the preparation of high quality AuNR colloids with tunable morphology, size, and optical response in a reproducible and scalable manner.

Surface van der Waals forces in a nutshell.

Most often in chemical physics, long range van der Waals surface interactions are approximated by the exact asymptotic result at vanishing distance, the well known additive approximation of London dispersion forces due to Hamaker. However, the description of retardation effects that is known since the time of Casimir is completely neglected for the lack of a tractable expression. Here we show that it is possible to describe surface van der Waals forces at arbitrary distances in one single simple equation. The result captures the long sought crossover from non-retarded (London) to retarded (Casimir) interactions, the effect of polarization in condensed media, and the full suppression of retarded interactions at large distance. This is achieved with similar accuracy and the same material properties that are used to approximate the Hamaker constant in conventional applications. The results show that at ambient temperature, retardation effects significantly change the power law exponent of the conventional Hamaker result for distances of just a few nanometers.

Nanocapillarity and Liquid Bridge-Mediated Force between Colloidal Nanoparticles.

In this work, we probe the concept of interface tension for ultrathin adsorbed liquid films on the nanoscale by studying the surface fluctuations of films down to the monolayer. Our results show that the spectrum of film height fluctuations of a liquid-vapor surface may be extended to ultrathin films provided we take into account the interactions of the substrate with the surface. Global fluctuations of the film height are described in terms of disjoining pressure, whereas surface deformations that are proportional to the interface area are accounted for by a film thickness-dependent surface tension. As a proof of concept, we model the capillary forces between colloidal nanoparticles held together by liquid bridges. Our results indicate that the classical equations for capillarity follow very precisely down to the nanoscale provided we account for the film height dependence of the surface tension.

Femtosecond laser reshaping yields gold nanorods with ultranarrow surface plasmon resonances.

The irradiation of gold nanorod colloids with a femtosecond laser can be tuned to induce controlled nanorod reshaping, yielding colloids with exceptionally narrow localized surface plasmon resonance bands. The process relies on a regime characterized by a gentle multishot reduction of the aspect ratio, whereas the rod shape and volume are barely affected. Successful reshaping can only occur within a narrow window of the heat dissipation rate: Low cooling rates lead to drastic morphological changes, and fast cooling has nearly no effect. Hence, a delicate balance must be achieved between irradiation fluence and surface density of the surfactant on the nanorods. This perfection process is appealing because it provides a simple, fast, reproducible, and scalable route toward gold nanorods with an optical response of exceptional quality, near the theoretical limit.

Capillary wave theory of adsorbed liquid films and the structure of the liquid-vapor interface.

In this paper we try to work out in detail the implications of a microscopic theory for capillary waves under the assumption that the density is given along lines normal to the interface. Within this approximation, which may be justified in terms of symmetry arguments, the Fisk-Widom scaling of the density profile holds for frozen realizations of the interface profile. Upon thermal averaging of capillary wave fluctuations, the resulting density profile yields results consistent with renormalization group calculations in the one-loop approximation. The thermal average over capillary waves may be expressed in terms of a modified convolution approximation where normals to the interface are Gaussian distributed. In the absence of an external field we show that the phenomenological density profile applied to the square-gradient free energy functional recovers the capillary wave Hamiltonian exactly. We extend the theory to the case of liquid films adsorbed on a substrate. For systems with short-range forces, we recover an effective interface Hamiltonian with a film height dependent surface tension that stems from the distortion of the liquid-vapor interface by the substrate, in agreement with the Fisher-Jin theory of short-range wetting. In the presence of long-range interactions, the surface tension picks up an explicit dependence on the external field and recovers the wave vector dependent logarithmic contribution observed by Napiorkowski and Dietrich. Using an error function for the intrinsic density profile, we obtain closed expressions for the surface tension and the interface width. We show the external field contribution to the surface tension may be given in terms of the film's disjoining pressure. From literature values of the Hamaker constant, it is found that the fluid-substrate forces may be able to double the surface tension for films in the nanometer range. The film height dependence of the surface tension described here is in full agreement with results of the capillary wave spectrum obtained recently in computer simulations, and the predicted translation mode of surface fluctuations reproduces to linear order in field strength an exact solution of the density correlation function for the Landau-Ginzburg-Wilson Hamiltonian in an external field.

Vapour-liquid interfacial properties of square-well chains from density functional theory and Monte Carlo simulation.

The statistical associating fluid theory for attractive potentials of variable range (SAFT-VR) density functional theory (DFT) developed by [Gloor et al., J. Chem. Phys., 2004, 121, 12740-12759] is used to predict the interfacial behaviour of molecules modelled as fully-flexible square-well chains formed from tangentially-bonded monomers of diameter σ and potential range λ = 1.5σ. Four different model systems, comprising 4, 8, 12, and 16 monomers per molecule, are considered. In addition to that, we also compute a number of interfacial properties of molecular chains from direct simulation of the vapour-liquid interface. The simulations are performed in the canonical ensemble, and the vapour-liquid interfacial tension is evaluated using the wandering interface (WIM) method, a technique based on the thermodynamic definition of surface tension. Apart from surface tension, we also obtain density profiles, coexistence densities, vapour pressures, and critical temperature and density, paying particular attention to the effect of the chain length on these properties. According to our results, the main effect of increasing the chain length (at fixed temperature) is to sharpen the vapour-liquid interface and to increase the width of the biphasic coexistence region. As a result, the interfacial thickness decreases and the surface tension increases as the molecular chains get longer. The interfacial thickness and surface tension appear to exhibit an asymptotic limiting behaviour for long chains. A similar behaviour is also observed for the coexistence densities and critical properties. Agreement between theory and simulation results indicates that SAFT-VR DFT is only able to predict qualitatively the interfacial properties of the model. Our results are also compared with simulation data taken from the literature, including the vapour-liquid coexistence densities, vapour pressures, and surface tension.

Premelting-Induced Smoothening of the Ice-Vapor Interface.

We perform computer simulations of the quasiliquid layer of ice formed at the ice-vapor interface close to the ice Ih-liquid-vapor triple point of water. Our study shows that the two distinct surfaces bounding the film behave at small wavelengths as atomically rough and independent ice-water and water-vapor interfaces. For long wavelengths, however, the two surfaces couple, large scale parallel fluctuations are inhibited, and the ice-vapor interface becomes smooth. Our results could help explain the complex morphology of ice crystallites.

Mesoscopic Hamiltonian for the fluctuations of adsorbed Lennard-Jones liquid films.

We use Monte Carlo simulations of a Lennard-Jones fluid adsorbed on a short-range planar wall substrate to study the fluctuations in the thickness of the wetting layer, and we get a quantitative and consistent characterization of their mesoscopic Hamiltonian, H[ξ]. We have observed important finite-size effects, which were hampering the analysis of previous results obtained with smaller systems. The results presented here support an appealing simple functional form for H[ξ], close but not exactly equal to the theoretical nonlocal proposal made on the basis a generic density-functional analysis by Parry and coworkers. We have analyzed systems under different wetting conditions, as a proof of principle for a method that provides a quantitative bridge between the molecular interactions and the phenomenology of wetting films at mesoscopic scales.