crystals Article Nucleation and Post-Nucleation Growth in Diffusion- Controlled and Hydrodynamic Theory of Solidification Frigyes Podmaniczky 1, * and László Gránásy 1,2, *   Citation: Podmaniczky, F.; Gránásy, L. Nucleation and Post-Nucleation Growth in Diffusion-Controlled and Hydrodynamic Theory of Solidification. Crystals 2021, 11, 437. https://doi.org/10.3390/cryst11040437 Academic Editors: Wolfram Miller and Koichi Kakimoto Received: 13 March 2021 Accepted: 13 April 2021 Published: 17 April 2021 Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affil- iations. Copyright: © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/). 1 Wigner Research Centre for Physics, P.O. Box 49, H-1525 Budapest, Hungary 2 BCAST, Brunel University, Uxbridge, Middlesex UB8 3PH, UK * Correspondence: podmaniczky.frigyes@wigner.hu (F.P.); granasy.laszlo@wigner.hu (L.G.) Abstract: Two-step nucleation and subsequent growth processes were investigated in the framework of the single mode phase-field crystal model combined with diffusive dynamics (corresponding to colloid suspensions) and hydrodynamical density relaxation (simple liquids). It is found that independently of dynamics, nucleation starts with the formation of solid precursor clusters that consist of domains with noncrystalline ordering (ringlike projections are seen from certain angles), and regions that have amorphous structure. Using the average bond order parameter q 6 , we dis- tinguished amorphous, medium range crystallike order (MRCO), and crystalline local orders. We show that crystallization to the stable body-centered cubic phase is preceded by the formation of a mixture of amorphous and MRCO structures. We have determined the time dependence of the phase composition of the forming solid state. We also investigated the time/size dependence of the growth rate for solidification. The bond order analysis indicates similar structural transitions during solidification in the case of diffusive and hydrodynamic density relaxation. Keywords: classical density functional theory; molecular modelling; two-step nucleation; growth kinetics; hydrodynamic theory of freezing 1. Introduction Colloid suspensions are considered as model systems for simple molecular liquids [1,2]. A wealth of information is available for the crystallization of colloids including nucleation, which was obtained via tracing the trajectories of colloid particles via optical methods [24]. For example, in colloids, two-step nucleation, assisted by amorphous/liquid precursor, appears to be a general phenomenon [510]. Analogous experimental information is not accessible for highly undercooled simple liquids. This raises the question, to what extent the observations made on colloids are applicable to metallic liquids. There appear to be differences: e.g., crystal growth is usually diffusion controlled for colloid suspensions, whereas it may be interface controlled for a hypercooled pure liquid. The velocity of a flat crystal–liquid interface depends on the growth mechanism. In the case of diffusionless pro- cesses a steady state velocity is established, whereas in the diffusion-controlled processes the front velocity, v, is proportional to t 1/2 , where t is time [11,12]. The growth rate also depends on the curvature of the interface, as predicted by the classic kinetic model based on monomer attachment and detachment [13,14]. For spherical particles: v = 16D λ 2 3v at 4π 1/3 sinh v at 2k B T Δg 2γ r  (1) where D is the self-diffusion coefficient, λ the jump distance of atoms, v at the atomic volume, k B Boltzmann’s constant, T the temperature, Δg the volumetric free energy difference between the liquid and the solid, γ the solid-liquid interfacial free energy, and r the radius of the crystallite. Note that for the critical size, r*=2γ/Δg, the deterministic growth rate is zero, for larger particles it is positive, whereas for smaller particles it is negative. This Crystals 2021, 11, 437. https://doi.org/10.3390/cryst11040437 https://www.mdpi.com/journal/crystals