Extension of classical nucleation theory for uniformly sheared systems

Is icosahedral short-range order presented in supercooled transition metals?

Supercooled transition metals are characterized by the absence of long-range order and the presence of a specific short-range order in the arrangement of atoms. So, the presence of shoulders and broadenings at the second maximum in the experimentally measured quantity-in the static structure factorS(k)-is usually interpreted as a manifestation of the icosahedral (ideal or distorted) short-range order (ISRO-Icosahedral Short-Range Order). Icosahedral short-range order is a structure with fivefold symmetry in the arrangement of atoms, which can lead to the possibility of achieving deep supercooling. In this work, we study the local structural features of equilibrium and supercooled nickel melts under various cooling protocols (γ∈1010,1014 K s^-1) in order to clarify the mechanism of formation of the icosahedral short-range order in pure transition metals. Comprehensive studies of the properties of nickel melts were carried out using experiments on x-ray diffraction, large-scale molecular dynamics (MD) simulations with subsequent structural and cluster analysis. A good agreement was found between the results of x-ray diffraction data and the MD simulations results for an equilibrium nickel melt. It was found that the nickel melt is characterized by a short-range order, formed by fragments of icosahedra and distorted icosahedral clusters. It was revealed that the 'liquid-crystal' phase transition in nickel is accompanied by the transformation of distorted icosahedral clusters into clusters with thefcc/hcp-symmetry. It is shown that, in contrast to the Voronoi tessellation method, the cluster analysis method based on the (q4,q6) rotational invariants does not allow sufficiently correct identification of the distorted icosahedral short-range order in metal melts.

Adsorption Energetics of Amino Acid Analogs on Polymer/Water Interfaces Studied by All-Atom Molecular Dynamics Simulation and a Theory of Solutions

Energetics of adsorption was addressed with all-atom molecular dynamics simulation on the interfaces of poly(2-methoxyethyl acrylate) (PMEA), poly(methyl methacrylate) (PMMA), and poly(butyl acrylate) (PBA) with water. A wide variety of adsorbate solutes were examined, and the free energy of adsorption was computed with the method of energy representation. It was found that the adsorption free energy was favorable (negative) for all the combinations of solute and polymer, and among PMEA, PMMA, and PBA, the strongest adsorption was observed on PMMA for the hydrophobic solutes and on PMEA for the hydrophilic ones. According to the decomposition of the adsorption free energy into the contributions from polymer and water, it was seen that the polymer contribution is larger in magnitude with the solute size. The total free energy of adsorption was correlated well with the solvation free energy in bulk water only for hydrophobic solutes. The roles of the intermolecular interaction components such as electrostatic, van der Waals, and excluded-volume were further studied, and the electrostatic component was influential only in determining the polymer dependences of the adsorption propensities of hydrophilic solutes. The extent of adsorption was shown to be ranked by the van der Waals component in the solute-polymer interaction separately over the hydrophilic and hydrophobic solutes, with the excluded-volume effect from water pointed out to also drive the adsorption.

Water Dissolved in a Variety of Polymers Studied by Molecular Dynamics Simulation and a Theory of Solutions

The performance of a polymer medium as a separation membrane is determined by the dissolution free energy ΔG and diffusion coefficient D of the permeant. In this work, ΔG and D of water are investigated with all-atom molecular dynamics simulation in a wide variety of polymer species in the amorphous state. The computed ΔG is shown to agree well with the experimental value for linear homopolymers, and the degrees of polymerization of the homopolymers do not affect ΔG when they are beyond ∼10. The copolymers of ethylene-vinylidene difluoride, ethylene-vinyl acetate, and ethylene-acrylamide are then examined by changing the repeating patterns of the constituent monomers in both the periodic and graft forms. It is found that ΔG is determined primarily by the overall compositions of the monomers and is not affected by the copolymerization topology (periodic or graft). The hydrophobicity of the copolymer is enhanced, furthermore, when the hydrophobicity and hydrophilicity of the ethylene and non-ethylene parts are well contrasted and those parts are fragmented along the polymer chain. According to the computed D, the diffusivity of water tends to be larger when the (co)polymer is more hydrophobic and ΔG is more positive. D is actually seen to vary by orders of magnitude with the polymer structures, while the effect of the polymer species on the water permeation is stronger for ΔG than for D.

Universal Nucleation Behavior of Sheared Systems

Using molecular simulations and a modified classical nucleation theory, we study the nucleation, under flow, of a variety of liquids: different water models, Lennard-Jones, and hard sphere colloids. Our approach enables us to analyze a wide range of shear rates inaccessible to brute-force simulations. Our results reveal that the variation of the nucleation rate with shear is universal. A simplified version of the theory successfully captures the nonmonotonic temperature dependence of the nucleation behavior, which is shown to originate from the violation of the Stokes-Einstein relation.

Homogeneous nucleation of sheared liquids: advances and insights from simulations and theory

One of the most ubiquitous and technologically important phenomena in nature is the nucleation of homogeneous flowing systems. The microscopic effects of shear on a nucleating system are still imperfectly understood, although in recent years a consistent picture has emerged. The opposing effects of shear can be split into two major contributions for simple atomic and molecular liquids: increase of the energetic cost of nucleation, and enhancement of the kinetics. In this perspective, we describe the latest computational and theoretical techniques which have been developed over the past two decades. We collate and unify the overarching influences of shear, temperature, and supersaturation on the process of homogeneous nucleation. Experimental techniques and capabilities are discussed, against the backdrop of results from simulations and theory. Although we primarily focus on simple systems, we also touch upon the sheared nucleation of more complex systems, including glasses and polymer melts. We speculate on the promising directions and possible advances that could come to fruition in the future.

Homogeneous nucleation in a Poiseuille flow

Nucleation in a dynamical environment plays an important role in the synthesis and manufacturing of quantum dots and nanocrystals. In this work, we investigate the effects of fluid flow (low Reynolds number flow) on the homogeneous nucleation in a circular microchannel in the framework of the classical nucleation theory. The contributions of the configuration entropy from the momentum-phase space and the kinetic energy and strain energy of a microcluster are incorporated in the calculation of the change of the Gibbs free energy from a flow state without a microcluster to a flow state with a microcluster. An analytical equation is derived for the determination of the critical nucleus size. Using this analytical equation, an analytical solution of the critical nucleus size for the formation of a critical liquid nucleus is found. For the formation of a critical solid nucleus, the contributions from both the kinetic energy and the strain energy are generally negligible. We perform numerical analysis of the homogeneous nucleation of a sucrose microcluster in a representative volume element of an aqueous solution, which flows through a circular microchannel. The numerical results reveal the decrease of the critical nucleus size and the corresponding work of formation of a critical nucleus with the increase of the distance to axisymmetric axis for the same numbers of solvent atoms and solute atoms/particles.

Classical nucleation theory for the crystallization kinetics in sheared liquids

While statistical mechanics provides a comprehensive framework for the understanding of equilibrium phase behavior, predicting the kinetics of phase transformations remains a challenge. Classical nucleation theory (CNT) provides a thermodynamic framework to relate the nucleation rate to thermodynamic quantities such as pressure difference and interfacial tension through the nucleation work necessary to spawn critical nuclei. However, it remains unclear whether such an approach can be extended to the crystallization of driven melts that are subjected to mechanical stresses and flows. Here, we demonstrate numerically for hard spheres that the impact of simple shear on the crystallization rate can be rationalized within the CNT framework by an additional elastic work proportional to the droplet volume. We extract the local stress and strain inside solid droplets, which yield size-dependent values for the shear modulus that are about half of the bulk value. Finally, we show that for a complete description one also has to take into account the change of interfacial work between the strained droplet and the sheared liquid. From scaling reasons, we expect this extra contribution to dominate the work formation of small nuclei but become negligible compared to the elastic work for droplets composed of a few hundreds of particles.

Kinetics of crystalline nuclei growth in glassy systems

This work reports results for crystalline nuclei growth in glassy systems. The crystal growth laws rescaled onto the waiting times of critically-sized nuclei follow a unified dependence. The scaled crystal growth rate characteristics as functions of reduced temperature follow unified power-law dependencies.

Effects of shear flow on phase nucleation and crystallization

Classical nucleation theory offers a good framework for understanding the common features of new phase formation processes in metastable homogeneous media at rest. However, nucleation processes in liquids are ubiquitously affected by hydrodynamic flow, and there is no satisfactory understanding of whether shear promotes or slows down the nucleation process. We developed a classical nucleation theory for sheared systems starting from the molecular level of the Becker-Doering master kinetic equation and we analytically derived a closed-form expression for the nucleation rate. The theory accounts for the effect of flow-mediated transport of molecules to the nucleus of the new phase, as well as for the mechanical deformation imparted to the nucleus by the flow field. The competition between flow-induced molecular transport, which accelerates nucleation, and flow-induced nucleus straining, which lowers the nucleation rate by increasing the nucleation energy barrier, gives rise to a marked nonmonotonic dependence of the nucleation rate on the shear rate. The theory predicts an optimal shear rate at which the nucleation rate is one order of magnitude larger than in the absence of flow.

Shear melting at the crystal-liquid interface: Erosion and the asymmetric suppression of interface fluctuations

The influence of an applied shear on the planar crystal-melt interface is modeled by a nonlinear stochastic partial differential equation of the interface fluctuations. A feature of this theory is the asymmetric destruction of interface fluctuations due to advection of the crystal protrusions on the liquid side of the interface only. We show that this model is able to qualitatively reproduce the nonequilibrium coexistence line found in simulations. The impact of shear on spherical clusters is also addressed.

Scaling law for crystal nucleation time in glasses

Due to high viscosity, glassy systems evolve slowly to the ordered state. Results of molecular dynamics simulation reveal that the structural ordering in glasses becomes observable over "experimental" (finite) time-scale for the range of phase diagram with high values of pressure. We show that the structural ordering in glasses at such conditions is initiated through the nucleation mechanism, and the mechanism spreads to the states at extremely deep levels of supercooling. We find that the scaled values of the nucleation time, τ1 (average waiting time of the first nucleus with the critical size), in glassy systems as a function of the reduced temperature, T˜, are collapsed onto a single line reproducible by the power-law dependence. This scaling is supported by the simulation results for the model glassy systems for a wide range of temperatures as well as by the experimental data for the stoichiometric glasses at the temperatures near the glass transition.

Shear-accelerated crystallization in a supercooled atomic liquid

A bulk metallic glass forming alloy is subjected to shear flow in its supercooled state by compression of a short rod to produce a flat disk. The resulting material exhibits enhanced crystallization kinetics during isothermal annealing as reflected in the decrease of the crystallization time relative to the nondeformed case. The transition from quiescent to shear-accelerated crystallization is linked to strain accumulated during shear flow above a critical shear rate γ̇(c)≈0.3 s(-1) which corresponds to Péclet number, Pe∼O(1). The observation of shear-accelerated crystallization in an atomic system at modest shear rates is uncommon. It is made possible here by the substantial viscosity of the supercooled liquid which increases strongly with temperature in the approach to the glass transition. We may therefore anticipate the encounter of nontrivial shear-related effects during thermoplastic deformation of similar systems.