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.

Effect of Planar Interfaces on Nucleation in Melting and Crystallization

The effect of planar interfaces on nucleation (namely, on the work of critical cluster formation and their shape) is studied both for crystallization and melting. Advancing an approach formulated about 150 years ago by J. W. Gibbs for liquid phase formation at planar liquid−liquid interfaces, we show that nucleation of liquids in the crystal at crystal−vapor planar interfaces proceeds as a rule with a much higher rate compared to nucleation in the bulk of the crystal. Provided the surface tensions crystal−liquid (σcl), liquid−vapor (σlv), and crystal−vapor (σcv) obey the condition σcv=σcl+σlv, the work of critical cluster formation tends to zero; in the range σcv<σcl+σlv, it is less than one half of the work of critical cluster formation for bulk nucleation. The existence of a liquid−vapor planar interface modifies the work of critical cluster formation in crystal nucleation in liquids to a much less significant degree. The work of critical crystal cluster formation is larger than one half of the bulk value of the work of critical cluster formation, reaching this limit at σcv=σcl+σlv. The shape of the critical clusters can be described in both cases by spherical caps with a radius, R, and a width parameter, h. This parameter, h, is the distance from the cutting plane (coinciding with the crystal−vapor and liquid−vapor planar interface, respectively) to the top of the spherical cap. It varies for nucleation of a liquid in a crystal in the range (h/R)≤1 and for crystal nucleation in a liquid in the range 2≥(h/R)≥1. At σcv=σcl+σlv, the ratio (h/R) of the critical cluster for nucleation in melting tends to zero ((h/R)→0). At the same condition, the critical crystallite has the shape of a sphere located tangentially to the liquid−vapor interface inside the liquid ((h/R)≅2). We present experimental data which confirm the results of the theoretical analysis, and potential further developments of the theoretical approach developed here are anticipated.

Freezing of Dilute Aqueous-Alcohol Nanodroplets: The Effect of Molecular Structure

We investigate vapor-liquid nucleation and subsequent freezing of aqueous-alcohol nanodroplets containing 1-pentanol, 1-hexanol, and their 3-isomers. The aerosols are produced in a supersonic nozzle, where condensation and freezing are characterized by static pressure and Fourier transform Infrared (FTIR) spectroscopy measurements. At fixed water concentrations, the presence of alcohol enables particle formation at higher temperatures since both the equilibrium vapor pressure above the critical clusters and the cluster interfacial free energy are decreased relative to the pure water case. The disappearance of a small free OH peak, observed for pure water droplets, when alcohols are added and shifts in the CH peaks as a function of alcohol chain length reveal varying surface partitioning preferences of the alcohols. Changes in the FTIR spectra during freezing, as well as changes in the ice component derived from self-modeling curve resolution analysis, show that 1-hexanol and 1-pentanol perturb freezing less than their branched isomers do. This behavior may reflect the molecular footprints of the alcohols, the available surface area of the droplets, and not only alcohol solubility. The presence of alcohols also lowers the freezing temperature relative to that of pure water, but when there is clear evidence for the formation of ice, the ice nucleation rates change by less than a factor of ∼2-3 for all cases studied.

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.

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.

Localization model description of the interfacial dynamics of crystalline Cu and Cu64Zr36 metallic glass films

Recent studies of structural relaxation in Cu-Zr metallic glass materials having a range of compositions and over a wide range of temperatures and in crystalline UO₂ under superionic conditions have indicated that the localization model (LM) can predict the structural relaxation time τ_α of these materials from the intermediate scattering function without any free parameters from the particle mean square displacement ⟨r²⟩ at a caging time on the order of ps, i.e., the "Debye-Waller factor" (DWF). In the present work, we test whether this remarkable relation between the "fast" picosecond dynamics and the rate of structural relaxation τ_α in these model amorphous and crystalline materials can be extended to the prediction of the local interfacial dynamics of model amorphous and crystalline films. Specifically, we simulate the free-standing amorphous Cu₆₄Zr₃₆ and crystalline Cu films and find that the LM provides an excellent parameter-free prediction for τ_α of the interfacial region. We also show that the Tammann temperature, defining the initial formation of a mobile interfacial layer, can be estimated precisely for both crystalline and glass-forming solid materials from the condition that the DWFs of the interfacial region and the material interior coincide.

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.

Vapor-phase nucleation of n-pentane, n-hexane, and n-heptane: Critical cluster properties

The first and second nucleation theorems provide a way to determine the molecular content and excess internal energies of critical clusters, which rely solely on experimental nucleation rates measured at constant temperatures and supersaturations, respectively. Here, we report the size n^* and excess internal energy E_x(n^*) of n-pentane, n-hexane, and n-heptane critical clusters when particles form under the highly supersaturated conditions present in supersonic expansions. In summary, critical clusters contain from ∼2 to ∼11 molecules and exhibit the expected increase in the critical cluster size with increasing temperature and decreasing supersaturation. Surprisingly, the n^* values for all three alkanes appear to lie along a single line when plotted as a function of supersaturation. Within the framework of the capillarity approximation, the excess internal energies determined for the n-heptane critical clusters formed under the low temperature (∼150 K) conditions in our supersonic nozzle are reasonably consistent with those determined under higher temperature (∼250 K) conditions in the thermal diffusion cloud chamber by Rudek et al. [J. Chem. Phys. 105, 4707 (1996)].

An effective experimental method and apparatus for unsteady water vapor condensation investigation in high speed expansion flow

This paper presents an effective experimental method and system for mechanism study of unsteady water vapor condensation encountered in high speed expansion flow. We proposed an experimental method and designed a simplified expansion system to fulfill the study. Tunable diode laser absorption spectroscopy (TDLAS) and light sheet technique are integrated in the system to monitor the unsteady condensation process in a high speed expansion flow generated by the expansion system. Two near infrared water vapor absorption transitions (1395.0 nm and 1409.27 nm) and one near infrared methane absorption transition (1653.73 nm) are applied in the TDLAS measurement to measure the transient flow parameters during the condensation process. Using the experimental method, time dependent condensation processes are monitored with different expansion time scales. The light sheet results visually reveal the condensation phenomena during the expansion process, while TDLAS results quantitatively follow the condensation process. The experimental results are compared with computational fluid dynamics simulations and a good agreement between them is observed, which indicates that the presented experimental method and system is effective in investigating unsteady water vapor condensation in high speed expansion flow.

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.

How Cubic Can Ice Be?

Using an X-ray laser, we investigated the crystal structure of ice formed by homogeneous ice nucleation in deeply supercooled water nanodrops (r ≈ 10 nm) at ∼225 K. The nanodrops were formed by condensation of vapor in a supersonic nozzle, and the ice was probed within 100 μs of freezing using femtosecond wide-angle X-ray scattering at the Linac Coherent Light Source free-electron X-ray laser. The X-ray diffraction spectra indicate that this ice has a metastable, predominantly cubic structure; the shape of the first ice diffraction peak suggests stacking-disordered ice with a cubicity value, χ, in the range of 0.78 ± 0.05. The cubicity value determined here is higher than those determined in experiments with micron-sized drops but comparable to those found in molecular dynamics simulations. The high cubicity is most likely caused by the extremely low freezing temperatures and by the rapid freezing, which occurs on a ∼1 μs time scale in single nanodroplets.

Freezing of supercooledn-decane nanodroplets: from surface driven to frustrated crystallization

Decane condenses to form nanodroplets that freeze from the outside in, until it's so cold that fractal like frustrated crystals form instead.

On the determination of the crystal-vapor surface free energy, and why a Gaussian expression can be accurate for a system far from Gaussian

The crystal-vapor surface free energy γ is an important physical parameter governing physical processes, such as wetting and adhesion. We explore exact and approximate routes to calculate γ based on cleaving an intact crystal into non-interacting sub-systems with crystal-vapor interfaces. We do this by turning off the interactions, ΔV, between the sub-systems. Using the soft-core scheme for turning off ΔV, we find that the free energy varies smoothly with the coupling parameter λ, and a single thermodynamic integration yields the exact γ. We generate another exact method, and a cumulant expansion for γ by expressing the surface free energy in terms of an average of e(-βΔV) in the intact crystal. The second cumulant, or Gaussian approximation for γ is surprisingly accurate in most situations, even though we find that the underlying probability distribution for ΔV is clearly not Gaussian. We account for this fact by developing a non-Gaussian theory for γ and find that the difference between the non-Gaussian and Gaussian expressions for γ consist of terms that are negligible in many situations. Exact and approximate methods are applied to the (111) surface of a Lennard-Jones crystal and are also tested for more complex molecular solids, the surface of octane and nonadecane. Alkane surfaces were chosen for study because their crystal-vapor surface free energy has been of particular interest for understanding surface freezing in these systems.

Low-temperature Bessel beam trap for single submicrometer aerosol particle studies

We report on a new instrument for single aerosol particle studies at low temperatures that combines an optical trap consisting of two counter-propagating Bessel beams (CPBBs) and temperature control down to 223 K (-50 °C). The apparatus is capable of capturing and stably trapping individual submicrometer- to micrometer-sized aerosol particles for up to several hours. First results from studies of hexadecane, dodecane, and water aerosols reveal that we can trap and freeze supercooled droplets ranging in size from ~450 nm to 5500 nm (radius). We have conducted homogeneous and heterogeneous freezing experiments, freezing-melting cycles, and evaporation studies. To our knowledge, this is the first reported observation of the freezing process for levitated single submicrometer-sized droplets in air using optical trapping techniques. These results show that a temperature-controlled CPBB trap is an attractive new method for studying phase transitions of individual submicrometer aerosol particles.

Stretching and folding of 2-nanometer hydrocarbon rods

Linear alkanes CnH2n+2 in vacuum isolation are finite models for an infinite polyethylene chain. Using spontaneous Raman scattering in supersonic jet expansions for n = 13-21 in different spectral ranges, we determine the minimal chain length nh for the cohesion-driven folding of the preferred extended all-trans conformation into a hairpin structure. We treat fully stretched all-trans alkanes as molecular "nanorods" and derive Young's modulus E for the stretching of an isolated single-strand polyethylene fibre by extrapolating the longitudinal acoustic mode to infinite chain length. Two key quality parameters for accurate intra- and intermolecular force fields of hydrocarbons (nh = 18 ± 1, E = 305 ± 5 GPa) are thus derived with high accuracy from experimental spectroscopy.

Freezing of stressed bilayers and vesicles

We show using a minimalistic theoretical framework that phase transition decoupling in lipid bilayers is caused by a nonuniform stress profile due to an asymmetric distribution of lipids between the two leaflets. Applying this framework to vesicles, we demonstrate that their anomalous freezing is also caused by a stress asymmetry, but that this is due to lipid tail extension on freezing. Finally, we predict a previously unknown dependence of surface tension on temperature, find the phenomenon also in coarse grained molecular dynamics simulations, and suggest that it might have relevance in thermosensitive protein gating.