Strength of Alkane–Fluid Attraction Determines the Interfacial Orientation of Liquid Alkanes and Their Crystallization through Heterogeneous or Homogeneous Mechanisms

Origin of the central pit in hemimicelles of semifluorinated alkanes: How molecular dipoles and substrate deformation can determine supra-molecular morphology

HYPOTHESIS

Semifluorinated alkanes amphiphiles spontaneously form highly monodispersed hemimicelles at the surface of water. The origin of the formation and complex structure of these surprising supramolecular aggregates were only recently clarified using molecular dynamics simulations (MD). The existence of a pit at the center of these aggregates made up of almost 3000 molecules was indeed reproduced by the MD simulations, but not understood.

METHOD

A careful strategy of atomistic MD simulations comparing non-electrostatic force fields with force fields that include electrostatic forces, thus bearing an implicit or explicit dipole, allowed demonstrating the roles of dipolar interactions and interactions with the liquid subphase on the morphology of the aggregates.

FINDINGS

The simulation results clearly show that within the hemimicelles the strong molecular dipoles located at the CH2-CF2 junctions tend to align, leading to a collective shift of the PFAA molecules relatively to each other. This shift is responsible for the curvature of the hemimicelles and originates the central pit, provided the possibility of deforming the surface of the water sub-phase. Comparisons with non-electrostatic force field results further contribute to understand the origin of the self-assembling process. The results directly connect for the first time a molecular property with a mesoscopic structural feature. Given the molecular simplicity of these "primitive" amphiphiles compared to the common hydrophilic/hydrophobic surfactants, the results contribute to understand self-assembly in general.

Effects of molecular size and orientation on the interfacial properties and wetting behavior of water/n-alkane systems: a molecular-dynamics study

Molecular dynamics simulations (MD) are performed to study the interfacial structure/tension and wetting behavior of water/n-alkane systems (water/nC5 to water/nC16 where nCx = CxH(2x + 2)). In particular, we study complete-to-partial wetting transitions by changing the n-alkane chain length (N_C) at a constant temperature, T = 295 K. Simulations are carried out with a united-atom TraPPE model for n-alkanes and the TIP4P-2005 model of water. Simulation results are in excellent agreement with the initial spreading coefficients and contact angles calculated using experimental values of the surface and interfacial tensions. In addition, it has been determined that water/(nC5-nC7) and water/(nC8-nC16), respectively, exhibit complete and partial initial wetting modes. Simulations show that the interfacial structures of water/(nC5-nC7) are different from water/(nC8-nC16) systems. In the latter, water preferentially orients near the interface to increase the number of hydrogen bonds and the charge and mass densities. Moreover, the orientation of n-alkane molecules at water/(nC8-nC16) interfaces has a long-range persistence, resulting in layered structures that increase with N_C. In addition, simulation results of the orientational order parameter S_z show alignment behavior of the n-alkane molecules with respect to the interfaces. Simulations predict that the central segments of n-alkane are strongly packed in the interfaces while the end segments (methyl groups) form smaller peaks in the outer edge of the layer. This observation confirms the "horseshoe" or "C-shaped" structure of n-alkane molecules in the water/n-alkane interfaces. At constant temperature, the interface widths of both water and the n-alkanes decrease with increasing n-alkane molecular length. These results suggest that increasing the n-alkane chain length affects the water/n-alkane interfacial properties in a manner similar to that of cooling.

A gel lattice alters the phase state of a solvent

Some low-molecular-weight substances are able to self-assemble into fiber-like structures (strings) to form gels. One of the examples of such substances is trifluoroacetylated alpha-aminoalcohols (TFAAAs) able to gelate in many organic solvents. Here we report the formation and describe the properties of a layer of an altered solvent covering the strings' surface. The altered solvent layer has a different refractive index and melts at a temperature about several degrees lower than that of the bulk solvent. Moreover, the bulk solvent's melting temperature was also decreased by values far beyond the one expected according to Raoult's law. Based on the Gibbs-Thomson equation it is possible to derive the thickness of the special layer as well as the average gel lattice parameters, which were very stable across the variety of systems investigated.

Water/oil interfacial tension reduction - an interfacial entropy driven process

The interfacial tension (IFT) of a fluid-fluid interface plays an important role in a wide range of applications and processes. When low IFT is desired, surface active compounds (e.g. surfactants) can be added to the system. Numerous attempts have been made to relate changes in IFT arising from such compounds to the specific nature of the interface. However, the IFT is controlled by an interplay of factors such as temperature and molecular structure of surface-active compounds, which make it difficult to predict IFT as those conditions change. In this study, we present the results from molecular dynamics simulations revealing the specific role surfactants play in IFT. We find that, in addition to reducing direct contact between the two fluids, surfactants serve to increase the disorder at the interface (related to interfacial entropy) and consequently reduce the water/oil IFT, especially when surfactants are present at high surface density. Our results suggest that surfactants that yield more disordered interfacial films (e.g. with flexible and/or unsaturated tails) reduce the water/oil IFT more effectively than surfactants which yield highly ordered interfacial films. Our results shed light on some of the factors that control IFT and could have important practical implications in industrial applications such as the design of cosmetics, food products, and detergents.

Molecular dynamics study of wetting of alkanes on water: from high temperature to the supercooled region and the influence of second inflection points of interfacial tensions

To explore the wetting behavior of alkanes on bulk water interfaces, molecular dynamics (MD) simulations were carried out for united-atom PYS alkane models, and for SPC/E and TIP4P/2005 water models over a wide temperature range. The MD results at each temperature were used to find (1) the surface tension of the alkanes (octane, nonane) and water, and (2) the interfacial tensions of the alkane-water systems. These quantities were then used to calculate the spreading coefficient (S) and contact angle (θc) for each alkane on water. At higher temperatures, the contact angle of octane and nonane on water is found to behave in accord with conventional expectations, i.e., it decreases with increasing temperature for both water models as each system approaches the usual high-temperature transition to perfect wetting. At lower temperatures, we found an unusual temperature dependence of S and θc for each PYS alkane on SPC/E water. In contrast to conventional expectations, θc decreases with a decrease in the temperature. For octane-SPC/E water, this unusual behavior of θc occurs due to the presence of second inflection points (SIP) in the vapor-water and the octane-water interfacial tensions, whereas the SIP effect is much less important for the nonane-water system. The unusual temperature dependence of θc observed for nonane on SPC/E water is also found for nonane on TIP4P/2005 water. On the other hand, such unusual wetting behavior has not been observed in the PYS octane-TIP4P/2005 water system, except possibly for the two lowest temperatures studied.

Coalescence dynamics in oil-in-water emulsions at elevated temperatures

Emulsion stability in a flow field is an extremely important issue relevant for many daily-life applications such as separation processes, food manufacturing, oil recovery etc. Microfluidic studies can provide micro-scale insight of the emulsion behavior but have primarily focussed on droplet breakup rather than on droplet coalescence. The crucial impact of certain conditions such as increased pressure or elevated temperature frequently used in industrial processes is completely overlooked in such micro-scale studies. In this work, we investigate droplet coalescence in flowing oil-in-water emulsions subjected to higher than room temperatures namely between 20 to 70 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{\circ }$$\end{document}∘C. We use a specifically designed lab-on-a-chip application for this purpose. Coalescence frequency is observed to increase with increasing temperature. We associate with this observation the change in viscosity at higher temperatures triggering a stronger perturbation in the thin aqueous film separating the droplets. Using the scaling law for rupture time of such a thin film, we establish a mechanism leading to a higher coalescence frequency at elevated temperatures.

Is Ice Nucleation by Organic Crystals Nonclassical? An Assessment of the Monolayer Hypothesis of Ice Nucleation

Potent ice nucleating organic crystals display an increase in nucleation efficiency with pressure and memory effect after pressurization that set them apart from inorganic nucleants. These characteristics were proposed to arise from an ordered water monolayer at the organic-water interface. It was interpreted that ordering of the monolayer is the limiting step for ice nucleation on organic crystals, rendering their mechanism of nucleation nonclassical. Despite the importance of organics in atmospheric ice nucleation, that explanation has never been investigated. Here we elucidate the structure of interfacial water and its role in ice nucleation at ambient pressure on phloroglucinol dihydrate, the paradigmatic example of outstanding ice nucleating organic crystal, using molecular simulations. The simulations confirm the existence of an interfacial monolayer that orders on cooling and becomes fully ordered upon ice formation. The monolayer does not resemble any ice face but seamlessly connects the distinct hydrogen-bonding orders of ice and the organic surface. Although large ordered patches develop in the monolayer before ice nucleates, we find that the critical step is the formation of the ice crystallite, indicating that the mechanism is classical. We predict that the fully ordered, crystalline monolayer nucleates ice above -2 °C and could be responsible for the exceptional ice nucleation by the organic crystal at high pressures. The lifetime of the fully ordered monolayer around 0 °C, however, is too short to account for the memory effect reported in the experiments. The latter could arise from an increase in the melting temperature of ice confined by strongly ice-binding surfaces.

Mechanism of surface freezing of alkanes

Using molecular dynamics simulation of octane (C₈) and nonadecane (C₁₉), we probe the mechanism of n-alkane surface freezing, the appearance of a crystalline monolayer above the liquid at a temperature T_sf above the bulk freezing point T_f. Formation of a crystalline monolayer occurs robustly in these systems. When T_f > T_sf, the surface frozen phase is metastable with respect to the solid but persists for long periods for study in simulations. Surface freezing of both C₈ and C₁₉ is driven by significant energy-lowering when alkane chains become ordered along the surface normal, and we elucidate the origins of this phenomenon. The degree of configurational disorder in the surface frozen layer relative to the solid is much larger for C₈ compared to C₁₉. From the Gibbsian viewpoint, we extract the excess energy and entropy of the liquid and surface frozen phases. We also consider the surface frozen layer as an intervening third phase, the viewpoint taken in previous theoretical analyses. Here, we find significantly increased entropy of the surface frozen phase of C₈ associated with configurational disorder, while the energy and entropy of the surface frozen phase of C₁₉ are marginally different from the bulk solid. Finally, by combining our previously determined solid-vapor surface free energies of C₈ and C₁₉ with liquid-vapor surface tensions from this work, we eliminate wetting as a possible mechanism for C₈ surface freezing, but it remains a possibility for C₁₉. We analyze the molecular structure of the liquid, surface frozen, and solid surfaces and discuss its relevance to thermodynamic properties.

Pathway of orientational symmetry breaking in crystallization of short n-alkane droplets: A molecular dynamics study

We carry out molecular dynamics simulations by using an all-atom model to study the nucleation and crystallization of n-alkane droplets under three-dimensional and quasi-two-dimensional conditions. We focus on the development of orientational order of chains from a random state to a neatly ordered one. Two new methods, the map of symmetry breaking and the information entropy of chain orientations, are introduced to characterize the emerge and remelting phenomena of a primary nucleus at the early stage of crystallization. Stepwise nucleation, as well as the surface induced nucleation, of large droplets is observed. We elucidate the kinetic process of the formation of a primary nucleus and the rearrangement of every single molecule involved in a primary nucleus. We found that density fluctuation and orientational preordering are coupled together and occur simultaneously in nucleation. Our results show the pathway of orientational symmetry breaking in the crystallization of n-alkane droplets that are heuristic for the deeper understanding of the crystallization in more complex molecules such as polymers.

Impact of water on the lubricating properties of hexadecane at the nanoscale

Fluid lubricants are routinely used to reduce friction in a wide range of applications, from car engines to machinery and hard-disk drives. However, their efficiency can be significantly influenced by the ambient conditions they are exposed to, in particular humidity. Our understanding of the molecular mechanisms responsible for the well-documented impact of water on lubrication remains limited, hindering the improvement of tribological formulations. Here, we use Atomic Force Microscopy (AFM) and shear force spectroscopy to investigate the structural and dynamical behaviour of a model lubricant, hexadecane, confined between an AFM probe and a hydrophilic mica surface at different temperatures and humidities. We show that both the nanoscale structure and the tribological behaviour of the system are dominated by the nucleation of water nanodroplets at the interface. The process is favoured at higher temperature and can be explained with classical nucleation theory whereby the droplets become stable when larger than 20 nm to 50 nm size, depending on the ambient conditions. Below this threshold, a molecularly thin film of water molecules coats the surface uniformly. Highly localised shear measurements demonstrate a detrimental impact of the nanodroplets on shear with a twofold increase in the lubricated friction force. However, this can be mitigated by the adjunction of an amphiphilic additive, here oleic acid.

How Do Surfactants Control the Agglomeration of Clathrate Hydrates?

Clathrate hydrates can spontaneously form under typical conditions found in oil and gas pipelines. The agglomeration of clathrates into large solid masses plugs the pipelines, posing adverse safety, economic, and environmental threats. Surfactants are customarily used to prevent the aggregation of clathrate particles and their coalescence with water droplets. It is generally assumed that a large contact angle between the surfactant-covered clathrate and water is a key predictor of the antiagglomerant performance of the surfactant. Here we use molecular dynamic simulations to investigate the structure and dynamics of surfactant films at the clathrate-oil interface, and their impact on the contact angle and coalescence between water droplets and hydrate particles. In agreement with the experiments, the simulations predict that surfactant-covered clathrate-oil interfaces are oil wet but super-hydrophobic to water. Although the water contact angle determines the driving force for coalescence, we find that a large contact angle is not sufficient to predict good antiagglomerant performance of a surfactant. We conclude that the length of the surfactant molecules, the density of the interfacial film, and the strength of binding of its molecules to the clathrate surface are the main factors in preventing the coalescence and agglomeration of clathrate particles with water droplets in oil. Our analysis provides a molecular foundation to guide the molecular design of effective clathrate antiagglomerants.

CO2 condensation onto alkanes: unconventional cases of heterogeneous nucleation

The classical picture invoked for heterogeneous nucleation is frequently that of a liquid condensing onto an immiscible solid particle. Here, we examine heterogeneous nucleation of CO2 onto particles comprised of n-pentane or n-hexane under conditions where CO2 should be a solid and the seed particles may be liquid or solid. Although CO2 condensed under all but one of the six conditions investigated, these experiments do not easily fit into the framework of standard heterogeneous nucleation experiments. Rather they explore unconventional regimes of heterogeneous nucleation in which the state of the seed particle may both affect whether deposition can proceed, and, in turn, be influenced by the presence of the condensing species. The work complements the earlier work of Tanimura et al. [RSC Adv., 2015, 5, 105537-105550] that investigated CO2 condensation onto ice nanoparticles, by using seed particles comprised of non-polar compounds that form and freeze under conditions where CO2 is already supersaturated with respect to the solid ice. In some cases, the conditions for seed formation approach the limit of homogeneous CO2 nucleation. Vibrational spectroscopy measurements help pinpoint where CO2 starts to condense. Furthermore, these IR measurements suggest that the n-alkanes never freeze in the presence of CO2, even if the temperatures are well below those required for them to freeze when CO2 is absent. Over the temperature range 65 < T/K < 140, the conditions corresponding to the onset of CO2 heterogeneous nucleation on pre-existing seed particle almost all lie very close to the extrapolated vapor-liquid equilibrium line of CO2 for a broad range of seed materials.

Long-term deep-supercooling of large-volume water and red cell suspensions via surface sealing with immiscible liquids

Supercooling of aqueous solutions is a fundamentally and practically important physical phenomenon with numerous applications in biopreservation and beyond. Under normal conditions, heterogeneous nucleation mechanisms critically prohibit the simultaneous long-term (> 1 week), large volume (> 1 ml), and low temperatures (< -10 °C) supercooling of aqueous solutions. Here, we report on the use of surface sealing of water by an oil phase to significantly diminish the primary heterogeneous nucleation at the water/air interface. We achieve deep supercooling (down to -20 °C) of large volumes of water (up to 100 ml) for long periods (up to 100 days) simultaneously via this approach. Since oils are mixtures of various hydrocarbons we also report on the use of pure alkanes and primary alcohols of various lengths to achieve the same. Furthermore, we demonstrate the utility of deep supercooling via preliminary studies on extended (100 days) preservation of human red blood cells.

Computer Simulations: Essential Tools for Crystal Growth Studies

This special issue discusses recent advances in computer simulation studies of crystal growth. Crystal growth is a key to innovation in science and technology. Owing to recent progress in computer performance, computer simulation studies of crystal growth have become increasingly important. This special issue covers a variety of simulation methods, including the Monte Carlo, molecular dynamics, first-principles, multiscale, and continuum simulation methods, which are used for studies on the fundamentals and applications of crystal growth and related phenomena for different materials, such as hard-sphere systems, ice, organic crystals, semiconductors, and graphene.

Ice-Nucleating and Antifreeze Proteins Recognize Ice through a Diversity of Anchored Clathrate and Ice-like Motifs

Cold-adapted organisms produce antifreeze and ice-nucleating proteins to prevent and promote ice formation. The crystal structure of hyperactive bacterial antifreeze protein (AFP) MpAFP suggests that this protein binds ice through an anchored clathrate motif. It is not known whether other hyperactive AFPs and ice-nucleating proteins (INPs) use the same motif to recognize or nucleate ice. Here we use molecular simulations to elucidate the ice-binding motifs of hyperactive insect AFPs and a model INP of Pseudomonas syringae. We find that insect AFPs recognize ice through anchored clathrate motifs distinct from that of MpAFP. By performing simulations of ice nucleation by PsINP, we identify two distinct ice-binding sites on opposite sides of the β-helix. The ice-nucleating sequences identified in the simulations agree with those previously proposed for the closely related INP of Pseudomonas borealis based on the structure of the protein. The simulations indicate that these sites have comparable ice-nucleating efficiency, but distinct binding motifs, controlled by the amino acid sequence: one is an anchored clathrate and the other ice-like. We conclude that anchored clathrate and ice-like motifs can be equally effective for binding proteins to ice and promoting ice nucleation.

Pre-critical fluctuations and what they disclose about heterogeneous crystal nucleation

Heterogeneous crystal nucleation is ubiquitous in nature and at the heart of many industrial applications. At the molecular scale, however, major gaps in understanding this phenomenon persist. Here we investigate through molecular dynamics simulations how the formation of precritical crystalline clusters is connected to the kinetics of nucleation. Considering heterogeneous water freezing as a prototypical scenario of practical relevance, we find that precritical fluctuations connote which crystalline polymorph will form. The emergence of metastable phases can thus be promoted by templating crystal faces characteristic of specific polymorphs. As a consequence, heterogeneous classical nucleation theory cannot describe our simulation results, because the different substrates lead to the formation of different ice polytypes. We discuss how the issue of polymorphism needs to be incorporated into analysis and comparison of heterogeneous and homogeneous nucleation. Our results will help to interpret and analyze the growing number of experiments and simulations dealing with crystal polymorph selection.

Reaction Coordinate for Ice Crystallization on a Soft Surface

The control of assembly and crystallization of molecules is becoming increasingly important in chemistry, engineering, and materials sciences. Crystallization is also central to understand natural processes that include the formation of atmospheric ice and biomineralization. Organic surfaces, biomolecules, and even liquid/vapor interfaces can promote the nucleation of crystals. These soft surfaces present significant structural fluctuations, which have been shown to strongly impact the rate of crystallization. This raises the question of whether degrees of freedom of soft surfaces play a role in the reaction coordinate for crystal nucleation. Here we use molecular simulations to investigate the mechanism of ice nucleation promoted by an alcohol monolayer. Our analysis indicates that while the flexibility of the surface strongly depresses its ice nucleation ability, it does not play a role in the coordinate that controls the transformation from liquid to ice. We find that the variable that drives the transformation is the size of the crystalline cluster, the same as that for the homogeneous crystallization. We argue that this is a general result that arises from the separation of time scales between surface fluctuations and the crossing of the transition state barrier for crystallization.