Surface-induced crystallization in supercooled tetrahedral liquids

Ideal Conductor Model: An Analytical Finite-Size Correction for Nonequilibrium Molecular Dynamics Simulations of Ion Transport through Nanoporous Membranes

Modulating ion transport through nanoporous membranes is critical to many important chemical and biological separation processes. The corresponding transport timescales, however, are often too long to capture accurately using conventional molecular dynamics (MD). Recently, path sampling techniques, such as forward-flux sampling (FFS), have emerged as attractive alternatives for efficiently and accurately estimating arbitrarily long ionic passage times. Here, we use non-equilibrium MD and FFS to explore how the kinetics and mechanisms of pressure-driven chloride transport through a nanoporous graphitic membrane are affected by its lateral dimensions. We not only find ionic passage times and free energy barriers to decrease dramatically upon increasing the membrane surface area but also observe an abrupt and discontinuous change in the locus of the transition state. These strong finite size effects arise due to the cumulative effect of the periodic images of the leading ion entering the pore on the distribution of the induced excess charge at the membrane surface in the feed. By assuming that the feed is an ideal conductor, we analytically derive a finite size correction term that can be computed from the information obtained from a single simulation and successfully use it to obtain corrected free energy profiles with no dependence on the system size. We then estimate ionic passage times in the thermodynamic limit by assuming an Eyring-type dependence of rates on barriers with a size-independent prefactor. This approach constitutes a universal framework for removing finite size artifacts in molecular simulations of ion transport through nanoporous membranes and biological channel proteins.

Anisotropic Molecular Organization at a Liquid/Vapor Interface Promotes Crystal Nucleation with Polymorph Selection

The molecules at the surface of a liquid have different organization and dynamics from those in the bulk, potentially altering the rate of crystal nucleation and polymorphic selection, but this effect remains poorly understood. Here we demonstrate that nucleation at the surface of a pure liquid, d-arabitol, is vastly enhanced, by 12 orders of magnitude, and selects a different polymorph. The surface effect intensifies with cooling and can be inhibited by a dilute, surface-active second component. This phenomenon arises from the anisotropic molecular packing at the interface and its similarity to the surface-nucleating polymorph. Our finding is relevant for controlling the crystallization and polymorphism in any system with a significant interface such as nanodroplets and atmospheric water.

Inducing thermodynamically blocked atomic ordering via strongly driven nonequilibrium kinetics

Ultrafast light-matter interactions enable inducing exotic material phases by promoting access to kinetic processes blocked in equilibrium. Despite potential opportunities, actively using nonequilibrium kinetics for material discovery is limited by the poor understanding on intermediate states of driven systems. Here, using single-pulse time-resolved imaging with x-ray free-electron lasers, we found intermediate states of photoexcited bismuth nanoparticles that showed kinetically reversed surface ordering during ultrafast melting. This entropy-lowering reaction was further investigated by molecular dynamics simulations to reveal that observed kinetics were thermodynamically buried in equilibrium, which emphasized the critical role of electron-mediated ultrafast free-energy modification in inducing exotic material phases. This study demonstrated that ultrafast photoexcitations of electrons provide an efficient strategy to induce hidden material phases by overcoming thermodynamic barriers via nonequilibrium reaction pathways.

Weakly Anisotropic Dielectric Properties of Water Droplets at the Nanoscale

The slab-confined water at the nanoscale exhibits anomalous dielectric properties compared to bulk water, for example, significantly low dielectric constant. In this work, we study the dielectric properties of nanoscale water droplets at room temperature using molecular dynamics simulations. We find that the nanoscale water droplets feature weakly anisotropic dielectric constant: the radial component of dielectric constants is distinctly smaller than the tangential component although they both decrease with reducing droplet size in a similar way. Such dielectric behavior is closely related to the orientational preference of water molecules near the convex surface. The molecular dipole prefers to slightly orientate toward the interior of droplets in contrast to the out-of-plane preference for free-standing water films and slab-confined water, which suppresses the fluctuation of dipole moments in the radial direction. Meanwhile, it facilitates the formation of the open hydrogen-bond network in the surface layer and ultimately leads to the relatively weak suppression of tangential fluctuations. The differential suppression is responsible for the anisotropic dielectric constant of water droplets. This anisotropic characteristic is also found in dielectric relaxation: both the radial and the tangential relaxation are consistently slowed down upon approaching surface but the latter is universally slower.

Identification of critical nuclei in the rapid solidification via configuration heredity

The identification and characterization of critical nuclei is a long-standing issue in the rapid solidification of metals and alloys. An ambiguous description for their sizes and shapes used to lead to an overestimation or underestimation of homogeneous nucleation ratesITin the framework of classical nucleation theory (CNT). In this paper, a unique method able to distinguish the critical nucleus from numerous embryos is put forward on the basis of configuration heredities of clusters during rapid solidifications. As this technique is applied to analyze the formation and evolution of various fcc-Al single crystal clusters in a large-scale molecular dynamics simulation system, it is found that the sizen_cand geometrical configuration of critical nuclei as well as their liquid-solid interfacial structure can be determined directly. For the present deep super-cooled system with an undercooling ofTm=0.42Tmcal, the average size of critical nuclei is demonstrated to benc̄≈26, but most of which are non-spherical lamellae. Also, their liquid-solid interfaces are revealed to be not an fcc-liquid duplex-phase interface but an fcc/hcp-liquid multi-phase structure. These findings shed some lights on the CNT, and a good agreement with previous simulations and experiments inITindicates this technique can be used to explore the early-stage of nucleation from atomistic levels.

Role of Nanoscale Interfacial Proximity in Contact Freezing in Water

Contact freezing is a mode of atmospheric ice nucleation in which a collision between a dry ice nucleating particle (INP) and a water droplet results in considerably faster heterogeneous nucleation. The molecular mechanism of such an enhancement is, however, still a mystery. While earlier studies had attributed it to collision-induced transient perturbations, recent experiments point to the pivotal role of nanoscale proximity of the INP and the free interface. By simulating the heterogeneous nucleation of ice within INP-supported nanofilms of two model water-like tetrahedral liquids, we demonstrate that such nanoscale proximity is sufficient for inducing rate increases commensurate with those observed in contact freezing experiments, but only if the free interface has a tendency to enhance homogeneous nucleation. Water is suspected of possessing this latter property, known as surface freezing propensity. Our findings therefore establish a connection between the surface freezing propensity and kinetic enhancement during contact nucleation. We also observe that faster nucleation proceeds through a mechanism markedly distinct from classical heterogeneous nucleation, involving the formation of hourglass-shaped crystalline nuclei that conceive at either interface and that have a lower free energy of formation due to the nanoscale proximity of the interfaces and the modulation of the free interfacial structure by the INP. In addition to providing valuable insights into the physics of contact nucleation, our findings can assist in improving the accuracy of heterogeneous nucleation rate measurements in experiments and in advancing our understanding of ice nucleation on nonuniform surfaces such as organic, polymeric, and biological materials.

Mechanisms of ice formation and propagation on superhydrophobic surfaces: A review

Icephobic surfaces, used as passive anti-icing materials, are in high demand due to the costs, damage, and loss of equipment and lives related to ice formation on outdoor surfaces. The proper design of icephobic surfaces is intertwined with the need for a profound understanding of ice formation processes and how ice propagates over a surface. Ice formation (ice nucleation) and interdroplet freezing propagation are processes that determine the onset of freezing and complete ice coverage on a surface, respectively. Evaluating the nature of these phenomena, along with their interactions with substrate and environmental factors, can offer a step toward designing surfaces having an improved icephobic performance. This review paper is organized to discuss ice nucleation and rate, preferable locations of nucleation, and favorable pathways of freezing (desublimation and condensation-freezing) on superhydrophobic surfaces. Furthermore, as the propagation of ice over a substrate plays a more deterministic role for the complete freezing coverage of a surface than that of ice formation, this review also elucidates possible mechanisms of ice propagation, theoretical backgrounds, and strategies to control this propagation using surface characteristics.

Thermal Effects on the Crystallization Kinetics, and Interfacial Adhesion of Single-Crystal Phase-Change Gallium

Although substrates play an important role upon crystallization of supercooled liquids, the influences of surface temperature and thermal property have remained elusive. Here, the crystallization of supercooled phase-change gallium (Ga) on substrates with different thermal conductivity is studied. The effect of interfacial temperature on the crystallization kinetics, which dictates thermo-mechanical stresses between the substrate and the crystallized Ga, is investigated. At an elevated surface temperature, close to the melting point of Ga, an extended single-crystal growth of Ga on dielectric substrates due to layering effect and annealing is realized without the application of external fields. Adhesive strength at the interfaces depends on the thermal conductivity and initial surface temperature of the substrates. This insight can be applicable to other liquid metals for industrial applications, and sheds more light on phase-change memory crystallization.

Liquid infused surfaces with anti-icing properties

Ice accretion on solid surfaces, a ubiquitous phenomenon that occurs in winter, brings much inconvenience to daily life and can even cause serious catastrophes. Icephobic surfaces, a passive way of processing surfaces to prevent surface destruction from ice accumulation, have attracted much attention from scientists because of their special ice-repellent properties, and many efforts have been made to rationally design durable icephobic coatings. This review is aimed at providing a brief and crucial overview of ice formation processes and feasible de-icing strategies. Here, the excellent anti-icing performance of liquid infused surfaces (LIS) inspired from Nepenthes is emphatically introduced. After a short introduction, the recent progresses in ice nucleation theory and ice adhesion decrease mechanism are comprehensively reviewed to gain a general understanding of the long freeze process and low ice adhesion on LIS. Subsequently, the anti-icing performance of LIS is systematically evaluated from four aspects regarding water repellence, condensation-frosting, long freeze process, and low ice adhesion. Finally, this review focuses on discussing the advantages and disadvantages of LIS and the potential measures to eliminate and alleviate these drawbacks.

Entropic colloidal crystallization pathways via fluid-fluid transitions and multidimensional prenucleation motifs

Complex crystallization pathways are common in protein crystallization, tetrahedrally coordinated systems, and biomineralization, where single or multiple precursors temporarily appear before the formation of the crystal. The emergence of precursors is often explained by a unique property of the system, such as short-range attraction, directional bonding, or ion association. But, structural characteristics of the prenucleation phases found in multistep crystallization remain unclear, and models are needed for testing and expanding the understanding of fluid-to-solid ordering pathways. Here, we report 3 instances of 2-step crystallization of hard-particle fluids. Crystallization in these systems proceeds via a high-density precursor fluid phase with prenucleation motifs in the form of clusters, fibers and layers, and networks, respectively. The density and diffusivity change across the fluid-fluid phase transition increases with motif dimension. We observe crystal nucleation to be catalyzed by the interface between the 2 fluid phases. The crystals that form are complex, including, notably, a crystal with 432 particles in the cubic unit cell. Our results establish the existence of complex crystallization pathways in entropic systems and reveal prenucleation motifs of various dimensions.

Thermodynamically driven assemblies and liquid-liquid phase separations in biology

The sustenance of life depends on the high degree of organization that prevails through different levels of living organisms, from subcellular structures such as biomolecular complexes and organelles to tissues and organs. The physical origin of such organization is not fully understood, and even though it is clear that cells and organisms cannot maintain their integrity without consuming energy, there is growing evidence that individual assembly processes can be thermodynamically driven and occur spontaneously due to changes in thermodynamic variables such as intermolecular interactions and concentration. Understanding the phase separation in vivo requires a multidisciplinary approach, integrating the theory and physics of phase separation with experimental and computational techniques. This paper aims at providing a brief overview of the physics of phase separation and its biological implications, with a particular focus on the assembly of membraneless organelles. We discuss the underlying physical principles of phase separation from its thermodynamics to its kinetics. We also overview the wide range of methods utilized for experimental verification and characterization of phase separation of membraneless organelles, as well as the utility of molecular simulations rooted in thermodynamics and statistical physics in understanding the governing principles of thermodynamically driven biological self-assembly processes.

Free energy landscape and molecular pathways of gas hydrate nucleation

Despite the significance of gas hydrates in diverse areas, a quantitative knowledge of hydrate formation at a molecular level is missing. The impediment to acquiring this understanding is primarily attributed to the stochastic nature and ultra-fine scales of nucleation events, posing a great challenge for both experiment and simulation to explore hydrate nucleation. Here we employ advanced molecular simulation methods, including forward flux sampling (FFS), p_B histogram analysis, and backward flux sampling, to overcome the limit of direct molecular simulation for exploring both the free energy landscape and molecular pathways of hydrate nucleation. First we test the half-cage order parameter (H-COP) which we developed for driving FFS, through conducting the p_B histogram analysis. Our results indeed show that H-COP describes well the reaction coordinates of hydrate nucleation. Through the verified order parameter, we then directly compute the free energy landscape for hydrate nucleation by combining both forward and backward flux sampling. The calculated stationary distribution density, which is obtained independently of nucleation theory, is found to fit well against the classical nucleation theory (CNT). Subsequent analysis of the obtained large ensemble of hydrate nucleation trajectories show that although on average, hydrate formation is facilitated by a two-step like mechanism involving a gradual transition from an amorphous to a crystalline structure, there also exist nucleation pathways where hydrate crystallizes directly, without going through the amorphous stage. The CNT-like free energy profile and the structural diversity suggest the existence of multiple active transition pathways for hydrate nucleation, and possibly also imply the near degeneracy in their free energy profiles among different pathways. Our results thus bring a new perspective to the long standing question of how hydrates crystallize.

Chemical ordering and crystal nucleation at the liquid surface: A comparison of Cu50Zr50 and Ni50Al50 alloys

We study the influence of the liquid-vapor surface on the crystallization kinetics of supercooled metal alloys. While a good glass former, Cu₅₀Zr₅₀, shows no evidence of surface enhancement of crystallization, Ni₅₀Al₅₀ exhibits an increased rate of crystallization due to heterogeneous nucleation at the free liquid surface. The difference in the compositional fluctuations at the interface is proposed as the explanation of the distinction between the two alloys. Specifically, we observe compositional ordering at the surface of Ni₅₀Al₅₀, while the Cu₅₀Zr₅₀ alloy only exhibits a diffuse adsorption of the Cu at the interface. We argue that the general difference in composition susceptibilities at planar surfaces represents an important factor in understanding the difference in the glass forming ability of the two alloys.

Subsurface nucleation of supercooled acetaminophen

Nucleation of acetaminophen glass at room temperature originates from ∼50 μm below the surface.

Perspective on Structural Evolution and Relations with Thermophysical Properties of Metallic Liquids

The relationship between the structural evolution and properties of metallic liquids is a long-standing hot issue in condensed-matter physics and materials science. Here, recent progress is reviewed in several fundamental aspects of metallic liquids, including the methods to study their atomic structures, liquid-liquid transition, physical properties, fragility, and their correlations with local structures, together with potential applications of liquid metals at room temperature. Involved with more experimentally and theoretically advanced techniques, these studies provide more in-depth understanding of the structure-property relationship of metallic liquids and promote the design of new metallic materials with superior properties.

Reversible crystalline-amorphous phase transformation in Si nanosheets with lithi-/delithiation

Silicon (Si) has a large theoretical capacity of 4200 mAhg^-1 and has great potential as a high-performance anode material for Li ion batteries (LIBs). Meanwhile, nanostructures can exploit the potential of Si and, accordingly, many zero-dimensional (0D) and one-dimensional (1D) Si nanostructures have been studied. Herein, we report on two-dimensional (2D) Si nanostructures, Si nanosheets (SiNSs), as anodes for LIBs. These 2D Si nanostructures, with a thickness as low 5 nm and widths of several micrometers, show reversible crystalline-amorphous phase transformations with the lithi-/delithiation by the dimensionality of morphology and large surface area. The reversible crystalline-amorphous phase transformation provides a structural stability of Li⁺ insertions and makes SiNSs promising candidates for reliable high-performance LIBs anode materials.

Effect of antifreeze protein on heterogeneous ice nucleation based on a two-dimensional random-field Ising model

Antifreeze proteins (AFPs) are the key biomolecules that protect many species from suffering the extreme conditions. Their unique properties of antifreezing provide the potential of a wide range of applications. Inspired by the present experimental approaches of creating an antifreeze surface by coating AFPs, here we present a two-dimensional random-field lattice Ising model to study the effect of AFPs on heterogeneous ice nucleation. The model shows that both the size and the free-energy effect of individual AFPs and their surface coverage dominate the antifreeze capacity of an AFP-coated surface. The simulation results are consistent with the recent experiments qualitatively, revealing the origin of the surprisingly low antifreeze capacity of an AFP-coated surface when the coverage is not particularly high as shown in experiment. These results will hopefully deepen our understanding of the antifreeze effects and thus be potentially useful for designing novel antifreeze coating materials based on biomolecules.

Computational investigation of surface freezing in a molecular model of water

Water freezes in a wide variety of low-temperature environments, from meteors and atmospheric clouds to soil and biological cells. In nature, ice usually nucleates at or near interfaces, because homogenous nucleation in the bulk can only be observed at deep supercoolings. Although the effect of proximal surfaces on freezing has been extensively studied, major gaps in understanding remain regarding freezing near vapor-liquid interfaces, with earlier experimental studies being mostly inconclusive. The question of how a vapor-liquid interface affects freezing in its vicinity is therefore still a major open question in ice physics. Here, we address this question computationally by using the forward-flux sampling algorithm to compute the nucleation rate in a freestanding nanofilm of supercooled water. We use the TIP4P/ice force field, one of the best existing molecular models of water, and observe that the nucleation rate in the film increases by seven orders of magnitude with respect to bulk at the same temperature. By analyzing the nucleation pathway, we conclude that freezing in the film initiates not at the surface, but within an interior region where the formation of double-diamond cages (DDCs) is favored in comparison with the bulk. This, in turn, facilitates freezing by favoring the formation of nuclei rich in cubic ice, which, as demonstrated by us earlier, are more likely to grow and overcome the nucleation barrier. The films considered here are ultrathin because their interior regions are not truly bulk-like, due to their subtle structural differences with the bulk.

Fast Surface Dynamics of Metallic Glass Enable Superlatticelike Nanostructure Growth

Contrary to the formation of complicated polycrystals induced by general crystallization, a modulated superlatticelike nanostructure, which grows layer by layer from the surface to the interior of a Pd_{40}Ni_{10}Cu_{30}P_{20} metallic glass, is observed via isothermal annealing below the glass transition temperature. The generation of the modulated nanostructure can be solely controlled by the annealing temperature, and it can be understood based on the fast dynamic and liquidlike behavior of the glass surface. The observations have implications for understanding the glassy surface dynamics and pave a way for the controllable fabrication of a unique and sophisticated nanostructure on a glass surface to realize the properties' modification.

Silicon-wall interfacial free energy via thermodynamics integration

We compute the interfacial free energy of a silicon system in contact with flat and structured walls by molecular dynamics simulation. The thermodynamics integration method, previously applied to Lennard-Jones potentials [R. Benjamin and J. Horbach, J. Chem. Phys. 137, 044707 (2012)], has been extended and implemented in Tersoff potentials with two-body and three-body interactions taken into consideration. The thermodynamic integration scheme includes two steps. In the first step, the bulk Tersoff system is reversibly transformed to a state where it interacts with a structureless flat wall, and in a second step, the flat structureless wall is reversibly transformed into an atomistic SiO₂ wall. Interfacial energies for liquid silicon-wall interfaces and crystal silicon-wall interfaces have been calculated. The calculated interfacial energies have been employed to predict the nucleation mechanisms in a slab of liquid silicon confined by two walls and compared with MD simulation results.

Polymorph selection and nucleation pathway in the crystallization of Hertzian spheres

The crystallization process of Hertzian spheres is studied by means of molecular dynamics simulations in an NPT ensemble where the total number of particles N, the pressure P, and the temperature T are kept constant. It has been observed that the bond orientational ordering rather than the translational ordering (density) plays a primary role. The crystal polymorphs are determined by the state points. Under the conditions of small supercooling, the system is likely to be nucleated into crystals that have a preference for the metastable bcc structure, which can be regarded as a manifestation of the Alexander-McTague mechanism. In contrast, small nuclei are found to have a preference for fcc symmetry under conditions of a high degree of supercooling. Prestructured precursors that act as seeds and wet on the nuclei during nucleation always have a high degree of bcc-like ordering, despite different state points. The results above may provide a clue to the understanding of the crystallization process in core-softened particles.

Hydrogen-plasma-induced Rapid, Low-Temperature Crystallization of μm-thick a-Si:H Films

Being a low-cost, mass-production-compatible route to attain crystalline silicon, post-deposition crystallization of amorphous silicon has received intensive research interest. Here we report a low-temperature (300 °C), rapid (crystallization rate of ~17 nm/min) means of a-Si:H crystallization based on high-density hydrogen plasma. A model integrating the three processes of hydrogen insertion, etching, and diffusion, which jointly determined the hydrogenation depth of the excess hydrogen into the treated micrometer thick a-Si:H, is proposed to elucidate the hydrogenation depth evolution and the crystallization mechanism. The effective temperature deduced from the hydrogen diffusion coefficient is far beyond the substrate temperature of 300 °C, which implies additional driving forces for crystallization, i.e., the chemical annealing/plasma heating and the high plasma sheath electric field. The features of LFICP (low-frequency inductively coupled plasma) and LFICP-grown a-Si:H are also briefly discussed to reveal the underlying mechanism of rapid crystallization at low temperatures.

Computational investigation of structure, dynamics and nucleation kinetics of a family of modified Stillinger-Weber model fluids in bulk and free-standing thin films

In recent years, computer simulations have found increasingly widespread use as powerful tools for studying phase transitions in wide variety of systems. In the particular and very important case of aqueous systems, the commonly used force-fields tend to offer quite different predictions with respect to a wide range of thermodynamic and kinetic properties, including the ease of ice nucleation, the propensity to freeze at a vapor-liquid interface, and the existence of a liquid-liquid phase transition. It is thus of fundamental and practical interest to understand how different features of a given water model affect its thermodynamic and kinetic properties. In this work, we use the forward-flux sampling technique to study the crystallization kinetics of a family of modified Stillinger-Weber (SW) potentials with energy (ε) and length (σ) scales taken from the monoatomic water (mW) model, but with different tetrahedrality parameters (λ). By increasing λ from 21 to 24, we observe the nucleation rate increases by 48 orders of magnitude at a supercooling of ζ = T/Tm = 0.845. Using classical nucleation theory, we are able to demonstrate that this change can largely be accounted for by the increase in |Δμ|, the thermodynamic driving force. We also perform rate calculations in freestanding thin films of the supercooled liquid, and observe a crossover from surface-enhanced crystallization at λ = 21 to bulk-dominated crystallization for λ ≥ 22.

Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations

The nucleation of crystals in liquids is one of nature's most ubiquitous phenomena, playing an important role in areas such as climate change and the production of drugs. As the early stages of nucleation involve exceedingly small time and length scales, atomistic computer simulations can provide unique insights into the microscopic aspects of crystallization. In this review, we take stock of the numerous molecular dynamics simulations that, in the past few decades, have unraveled crucial aspects of crystal nucleation in liquids. We put into context the theoretical framework of classical nucleation theory and the state-of-the-art computational methods by reviewing simulations of such processes as ice nucleation and the crystallization of molecules in solutions. We shall see that molecular dynamics simulations have provided key insights into diverse nucleation scenarios, ranging from colloidal particles to natural gas hydrates, and that, as a result, the general applicability of classical nucleation theory has been repeatedly called into question. We have attempted to identify the most pressing open questions in the field. We believe that, by improving (i) existing interatomic potentials and (ii) currently available enhanced sampling methods, the community can move toward accurate investigations of realistic systems of practical interest, thus bringing simulations a step closer to experiments.

Ice nucleation on nanotextured surfaces: the influence of surface fraction, pillar height and wetting states

Ice nucleation and growth on nanostructured surfaces.

Thermotropic interface and core relaxation dynamics of liquid crystals in silica glass nanochannels: a dielectric spectroscopy study

We report dielectric relaxation spectroscopy experiments on two rod-like liquid crystals of the cyanobiphenyl family (5CB and 6CB) confined in tubular nanochannels with 7 nm radius and 340 micrometer length in a monolithic, mesoporous silica membrane. The measurements were performed on composites for two distinct regimes of fractional filling: monolayer coverage at the pore walls and complete filling of the pores. For the layer coverage a slow surface relaxation dominates the dielectric properties. For the entirely filled channels the dielectric spectra are governed by two thermally-activated relaxation processes with considerably different relaxation rates: a slow relaxation in the interface layer next to the channel walls and a fast relaxation in the core region of the channel filling. The strengths and characteristic frequencies of both relaxation processes have been extracted and analysed as a function of temperature. Whereas the temperature dependence of the static capacitance reflects the effective (average) molecular ordering over the pore volume and is well described within a Landau-de Gennes theory, the extracted relaxation strengths of the slow and fast relaxation processes provide an access to distinct local molecular ordering mechanisms. The order parameter in the core region exhibits a bulk-like behaviour with a strong increase in the nematic ordering just below the paranematic-to-nematic transition temperature TPN and subsequent saturation during cooling. By contrast, the surface ordering evolves continuously with a kink near TPN. A comparison of the thermotropic behaviour of the monolayer with the complete filling reveals that the molecular order in the core region of the pore filling affects the order of the peripheral molecular layers at the wall.

Nanoscale size effects in crystallization of metallic glass nanorods

Atomistic understanding of crystallization in solids is incomplete due to the lack of appropriate materials and direct experimental tools. Metallic glasses possess simple metallic bonds and slow crystallization kinetics, making them suitable to study crystallization. Here, we investigate crystallization of metallic glass-forming liquids by in-situ heating metallic glass nanorods inside a transmission electron microscope. We unveil that the crystallization kinetics is affected by the nanorod diameter. With decreasing diameters, crystallization temperature decreases initially, exhibiting a minimum at a certain diameter, and then rapidly increases below that. This unusual crystallization kinetics is a consequence of multiple competing factors: increase in apparent viscosity, reduced nucleation probability and enhanced heterogeneous nucleation. The first two are verified by slowed grain growth and scatter in crystallization temperature with decreasing diameters. Our findings provide insight into relevant length scales in crystallization of supercooled metallic glasses, thus offering accurate processing conditions for predictable metallic glass nanomolding.

Homogeneous nucleation and microstructure evolution in million-atom molecular dynamics simulation

Homogeneous nucleation from an undercooled iron melt is investigated by the statistical sampling of million-atom molecular dynamics (MD) simulations performed on a graphics processing unit (GPU). Fifty independent instances of isothermal MD calculations with one million atoms in a quasi-two-dimensional cell over a nanosecond reveal that the nucleation rate and the incubation time of nucleation as functions of temperature have characteristic shapes with a nose at the critical temperature. This indicates that thermally activated homogeneous nucleation occurs spontaneously in MD simulations without any inducing factor, whereas most previous studies have employed factors such as pressure, surface effect, and continuous cooling to induce nucleation. Moreover, further calculations over ten nanoseconds capture the microstructure evolution on the order of tens of nanometers from the atomistic viewpoint and the grain growth exponent is directly estimated. Our novel approach based on the concept of “melting pots in a supercomputer” is opening a new phase in computational metallurgy with the aid of rapid advances in computational environments.

Direct calculation of ice homogeneous nucleation rate for a molecular model of water

Ice formation is ubiquitous in nature, with important consequences in a variety of environments, including biological cells, soil, aircraft, transportation infrastructure, and atmospheric clouds. However, its intrinsic kinetics and microscopic mechanism are difficult to discern with current experiments. Molecular simulations of ice nucleation are also challenging, and direct rate calculations have only been performed for coarse-grained models of water. For molecular models, only indirect estimates have been obtained, e.g., by assuming the validity of classical nucleation theory. We use a path sampling approach to perform, to our knowledge, the first direct rate calculation of homogeneous nucleation of ice in a molecular model of water. We use TIP4P/Ice, the most accurate among existing molecular models for studying ice polymorphs. By using a novel topological approach to distinguish different polymorphs, we are able to identify a freezing mechanism that involves a competition between cubic and hexagonal ice in the early stages of nucleation. In this competition, the cubic polymorph takes over because the addition of new topological structural motifs consistent with cubic ice leads to the formation of more compact crystallites. This is not true for topological hexagonal motifs, which give rise to elongated crystallites that are not able to grow. This leads to transition states that are rich in cubic ice, and not the thermodynamically stable hexagonal polymorph. This mechanism provides a molecular explanation for the earlier experimental and computational observations of the preference for cubic ice in the literature.

Physics of icing and rational design of surfaces with extraordinary icephobicity

Icing of surfaces is commonplace in nature and technology, affecting everyday life and sometimes causing catastrophic events. Understanding (and counteracting) surface icing brings with it significant scientific challenges that requires interdisciplinary knowledge from diverse scientific fields such as nucleation thermodynamics and heat transfer, fluid dynamics, surface chemistry, and surface nanoengineering. Here we discuss key aspects and findings related to the physics of ice formation on surfaces and show how such knowledge could be employed to rationally develop surfaces with extreme resistance to icing (extraordinary icephobicity). Although superhydrophobic surfaces with micro-, nano-, or (often biomimetic) hierarchical roughnesses have shown in laboratory settings (under certain conditions) excellent repellency and low adhesion to water down to temperatures near or below the freezing point, extreme icephobicity necessitates additional important functionalities. Other approaches, such as lubricant-impregnated surfaces, exhibit both advantages and serious limitations with respect to icing. In all, a clear path toward passive surfaces with extreme resistance to ice formation remains a challenge, but it is one well worth undertaking. Equally important to potential applications is scalable surface manufacturing and the ability of icephobic surfaces to perform reliably and sustainably outside the laboratory under adverse conditions. Surfaces should possess mechanical and chemical stability, and they should be thermally resilient. Such issues and related research directions are also addressed in this article.

Inhomogeneous relaxation dynamics and phase behaviour of a liquid crystal confined in a nanoporous solid

We report filling-fraction dependent dielectric spectroscopy measurements on the relaxation dynamics of the rod-like nematogen 7CB condensed in 13 nm silica nanochannels. In the film-condensed regime, a slow interface relaxation dominates the dielectric spectra, whereas from the capillary-condensed state up to complete filling an additional, fast relaxation in the core of the channels is found. The temperature-dependence of the static capacitance, representative of the averaged, collective molecular orientational ordering, indicates a continuous, paranematic-to-nematic (P-N) transition, in contrast to the discontinuous bulk behaviour. It is well described by a Landau-de-Gennes free energy model for a phase transition in cylindrical confinement. The large tensile pressure of 10 MPa in the capillary-condensed state, resulting from the Young-Laplace pressure at highly curved liquid menisci, quantitatively accounts for a downward-shift of the P-N transition and an increased molecular mobility in comparison to the unstretched liquid state of the complete filling. The strengths of the slow and fast relaxations provide local information on the orientational order: the thermotropic behaviour in the core region is bulk-like, i.e. it is characterized by an abrupt onset of the nematic order at the P-N transition. By contrast, the interface ordering exhibits a continuous evolution at the P-N transition. Thus, the phase behaviour of the entirely filled liquid crystal-silica nanocomposite can be quantitatively described by a linear superposition of these distinct nematic order contributions.

Nucleation at the contact line observed on nanotextured surfaces

It has been conjectured that roughness plays a role in surface nucleation, the tendency for freezing to begin preferentially at the liquid-gas interface. Using high speed imaging, we sought evidence for freezing at the contact line on catalyst substrates with imposed characteristic length scales (texture). Length scales consistent with the critical nucleus size and with δ∼τ/σ, where τ is a relevant line tension and σ is the surface tension, range from nanometers to micrometers. It is found that nanoscale texture causes a shift in the nucleation of ice in supercooled water to the three-phase contact line, while microscale texture does not.

Suppression of sub-surface freezing in free-standing thin films of a coarse-grained model of water

Freezing in the vicinity of water-vapor interfaces is of considerable interest to a wide range of disciplines, most notably the atmospheric sciences. In this work, we use molecular dynamics and two advanced sampling techniques, forward flux sampling and umbrella sampling, to study homogeneous nucleation of ice in free-standing thin films of supercooled water. We use a coarse-grained monoatomic model of water, known as mW, and we find that in this model a vapor-liquid interface suppresses crystallization in its vicinity. This suppression occurs in the vicinity of flat interfaces where no net Laplace pressure in induced. Our free energy calculations reveal that the pre-critical crystalline nuclei that emerge near the interface are thermodynamically less stable than those that emerge in the bulk. We investigate the origin of this instability by computing the average asphericity of nuclei that form in different regions of the film, and observe that average asphericity increases closer to the interface, which is consistent with an increase in the free energy due to increased surface-to-volume ratios.

Probing methane hydrate nucleation through the forward flux sampling method

Understanding the nucleation of hydrate is the key to developing effective strategies for controlling methane hydrate formation. Here we present a computational study of methane hydrate nucleation, by combining the forward flux sampling (FFS) method and the coarse-grained water model mW. To facilitate the application of FFS in studying the formation of methane hydrate, we developed an effective order parameter λ on the basis of the topological analysis of the tetrahedral network. The order parameter capitalizes the signature of hydrate structure, i.e., polyhedral cages, and is capable of efficiently distinguishing hydrate from ice and liquid water while allowing the formation of different hydrate phases, i.e., sI, sII, and amorphous. Integration of the order parameter λ with FFS allows explicitly computing hydrate nucleation rates and obtaining an ensemble of nucleation trajectories under conditions where spontaneous hydrate nucleation becomes too slow to occur in direct simulation. The convergence of the obtained hydrate nucleation rate was found to depend crucially on the convergence of the spatial distribution for the spontaneously formed hydrate seeds obtained from the initial sampling of FFS. The validity of the approach is also verified by the agreement between the calculated nucleation rate and that inferred from the direct simulation. Analyzing the obtained large ensemble of hydrate nucleation trajectories, we show hydrate formation at 220 K and 500 bar is initiated by the nucleation events occurring in the vicinity of water-methane interface, and facilitated by a gradual transition from amorphous to crystalline structure. The latter provides the direct support to the proposed two-step nucleation mechanism of methane hydrate.

Low activation energy for the crystallization of amorphous silicon nanoparticles

We have experimentally determined the crystallization rate of plasma-produced amorphous silicon powder undergoing in-flight thermal annealing, and have found a significant reduction in the activation energy for crystallization compared to amorphous silicon thin films. This finding allows us to shed light onto the mechanism leading to the formation of high quality nanocrystals in non-thermal plasmas.

Ice nucleation at the nanoscale probes no man's land of water

At a given thermodynamic condition, nucleation events occur at a frequency that scales with the volume of the system. Therefore at the nanoscale, one may expect to obtain supercooled liquids below the bulk homogeneous nucleation temperature. Here we report direct computational evidence that in supercooled water nano-droplets ice nucleation rates are strongly size dependent and at the nanoscale they are several orders of magnitude smaller than in bulk water. Using a thermodynamic model based on classical nucleation theory, we show that the Laplace pressure is partially responsible for the suppression of ice crystallization. Our simulations show that the nucleation rates found for droplets are similar to those of liquid water subject to a pressure of the order of the Laplace pressure within droplets. Our findings aid the interpretation of molecular beam experiments and support the hypothesis of surface crystallization of ice in microscopic water droplets in clouds.

Crystalline nucleation in undercooled liquids: a Bayesian data-analysis approach for a nonhomogeneous Poisson process

A Bayesian data-analysis approach to data sets of maximum undercooling temperatures recorded in repeated melting-cooling cycles of high-purity samples is proposed. The crystallization phenomenon is described in terms of a nonhomogeneous Poisson process driven by a temperature-dependent sample nucleation rate J(T). The method was extensively tested by computer simulations and applied to real data for undercooled liquid Ge. It proved to be particularly useful in the case of scarce data sets where the usage of binned data would degrade the available experimental information.