Search

A widely accepted phenomenological rule states that solids with free surfaces cannot be overheated. In this work we discuss this statement critically under the light of the statistical thermodynamics of interfacial roughening transitions. Our results show that the basal face of ice as described by the TIP4P/Ice model can remain mechanically stable for more than one hundred nanoseconds when overheated by 1 K, and for several hundreds of nanoseconds at smaller overheating despite the presence of a significant quasi-liquid layer at the surface. Such time scales, which are often of little experimental significance, can become a concern for the determination of melting points by computer simulations using the direct coexistence method. In the light of this observation, we reinterpret computer simulations of ice premelting and show that current results for the TIP4P/Ice model all imply a scenario of incomplete surface melting. Using a thermodynamic integration path, we reassess our own estimates for the Laplace pressure difference between water and vapor. These calculations are used to measure the disjoining pressure of premelting liquid films and allow us to confirm a minimum of the interfacial free energy at finite premelting thickness of about one nanometer., Comment: Manuscript as version of record except for minor edits

Achieving the formation of target open crystalline lattices from colloidal particles is of paramount importance for their potential application in photonics. Examples of such desired structures are the diamond, tetrastack, and pyrochlore lattices. Here, we demonstrate that the self-assembly of tetravalent patchy particles results in the selective formation of cubic tetrastack crystals, both in the bulk and in the systems subjected to external fields exerted by the solid substrate. It is demonstrated that the presence of an external field allows for the formation of well-defined single crystals with a low density of defects. Moreover, depending on the strength of the applied external field, the mechanism of epitaxial growth changes. For weakly attractive external fields, the crystallization occurs in a similar manner as in the bulk, since the fluid does not wet the substrate. Nonetheless, the formed crystal is considerably better ordered than the crystals formed in bulk, since the surface induces the ordering in the first layer. On the other hand, it is demonstrated that the formation of well-ordered cubic tetrastack crystals is considerably enhanced by the increase in external field strength, and the formation of the thick crystalline film occurs via a series of layering transitions.

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 which 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., Comment: 5 pages, 1 figure

The endeavor to selectively fabricate a cubic diamond is challenging due to the formation of competing phases such as its hexagonal polymorph or others possessing similar free energy. The necessity to achieve this is of paramount importance since the cubic diamond is the only polymorph exhibiting a complete photonic bandgap, making it a promising candidate in view of photonic applications. Herein, we demonstrate that due to the presence of an external field and delicate manipulation of its strength we can attain selectivity in the formation of a cubic diamond in a one-component system comprised of designer tetrahedral patchy particles. The driving force of such a phenomenon is the structure of the first adlayer which is commensurate with the (110) face of the cubic diamond. Moreover, after a successful nucleation event, once the external field is turned off, the structure remains stable, paving an avenue for further post-synthetic treatment.

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 at 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 TIP4P/Ice compares very well with TIP4P/2005 and to 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.

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 non equilibrium 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 substrates 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 the atomic scale frictional behavior of ice is a combination of spontaneous premelting, pressure melting and frictional heating., Comment: 9 pages, 4 figures + supplementary material, 24 pgs, 14 figures, 5 tables