Mechanism of surface freezing of alkanes.

Using molecular dynamics simulation of octane (C8) and nonadecane (C19), we probe the mechanism of n-alkane surface freezing, the appearance of a crystalline monolayer above the liquid at a temperature Tsf above the bulk freezing point Tf. Formation of a crystalline monolayer occurs robustly in these systems. When Tf > Tsf, the surface frozen phase is metastable with respect to the solid but persists for long periods for study in simulations. Surface freezing of both C8 and C19 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 C8 compared to C19. 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 C8 associated with configurational disorder, while the energy and entropy of the surface frozen phase of C19 are marginally different from the bulk solid. Finally, by combining our previously determined solid-vapor surface free energies of C8 and C19 with liquid-vapor surface tensions from this work, we eliminate wetting as a possible mechanism for C8 surface freezing, but it remains a possibility for C19. We analyze the molecular structure of the liquid, surface frozen, and solid surfaces and discuss its relevance to thermodynamic properties.

Freezing of supercooled n-decane nanodroplets: from surface driven to frustrated crystallization.

Whether crystallization starts at the liquid-vapor interface or randomly throughout the bulk has been the subject of intense debate. In our earlier work, we investigated the freezing of supercooled nanodroplets of short chain (C8, C9) n-alkanes formed by homogeneous condensation in a supersonic nozzle. The rate at which the solid appeared suggested freezing starts at the droplet surface well before the rest of the droplet freezes. Experiments were, however, limited to a single condition for each compound and it was not clear whether freezing of n-alkanes always occurs in this two step manner. Here, we expand our work to include freezing of a third n-alkane, n-decane, and, furthermore, we vary the temperatures at which droplets are formed and freeze. The phase transitions are again characterized using three experimental techniques - pressure trace measurements (PTM), Fourier Transform Infrared Spectroscopy (FTIR), and Small Angle X-ray Scattering (SAXS). We also use Wide Angle X-ray Scattering (WAXS) to confirm, for the first time, the crystalline nature of our frozen n-alkane nanodroplets. As the temperature at which the droplets form and freeze decreases, the kinetics of the phase transition changes. At higher temperatures, the phase transition occurs in two steps characterized by different rates, whereas at lower temperatures we observe only a single step. Finally, in the lowest temperature experiment, where droplets start to form and freeze ∼50 K below the bulk melting temperature, we found that the particles develop a fractal structure and appear locked in a "frustrated" crystalline state.

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.

Experimental evidence for surface freezing in supercooled n-alkane nanodroplets.

Intermediate chain length (16 ≤i≤ 50) n-alkanes are known to surface freeze at temperatures that are up to three degrees higher than the equilibrium melting point [B. M. Ocko et al., Phys. Rev. E, 1997, 55, 3164-3182]. Our recent experimental results suggest that highly supercooled nanodroplets of n-octane and n-nonane also surface freeze, and subsequently bulk crystallization occurs. The data yield surface and bulk nucleation rates on the order of ~10(15) cm(-2) s(-1) and ~10(22) cm(-3) s(-1), respectively, at temperatures between 180 K and 200 K. Molecular dynamics simulations at the united atom level were used to follow the freezing of a supercooled n-octane drop and show that an ordered monolayer develops on the surface of the droplet almost immediately, and the rest of the droplet then freezes in a layer-by-layer manner.

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.