• The effect of Sm doping on the fcc-Al nucleation was investigated in Al-Sm liquids with low Sm concentrations (xSm) with molecular dynamics simulations. The nucleation in the moderately undercooled liquid is achieved by the recently developed persistent-embryo method. Systematically computing the nucleation rate with different xSm (xSm=0%, 1%, 2%, 3%, 5%) at 700 K, we found Sm dopant reduces the nucleation rate by up to 25 orders of magnitudes with only 5% doping concentration. This effect is mostly associated with the increase in the free energy barrier with a minor contribution from suppression of the attachment to the nucleus caused by Sm doping.
  • A metastable {\epsilon}-Al60Sm11 phase appears during the initial devitrification of as-quenched Al-10.2 at.% Sm glasses. The {\epsilon} phase is nonstoichiometric in nature since Al occupation is observed on the 16f Sm lattice sites. Scanning transmission electron microscopic images reveal profound spatial correlation of Sm content on these sites, which cannot be explained by the "average crystal" description from Rietveld analysis of diffraction data. Thermodynamically favorable configurations, established by Monte Carlo (MC) simulations based on a cluster-expansion model, also give qualitatively different correlation functions from experimental observations. On the other hand, molecular dynamics simulations of the growth of {\epsilon}-Al60Sm11 in undercooled liquid show that when the diffusion range of Sm is limited to ~ 4 {\AA}, the correlation function of the as-grown crystal structure agrees well with that of the STEM images. Our results show that kinetic effects, especially the limited diffusivity of Sm atoms plays the fundamental role in determining the nonstoichiometric site occupancies of the {\epsilon}-Al60Sm11 phase during the crystallization process.
  • The crystal nucleation from liquid in most cases is too rare to be accessed within the limited timescales of the conventional molecular dynamics (MD) simulation. Here, we developed a "persistent embryo" method to facilitate crystal nucleation in MD simulations by preventing small crystal embryos from melting using external spring forces. We applied this method to the pure Ni case for a moderate undercooling where no nucleation can be observed in the conventional MD simulation, and obtained nucleation rate in good agreement with the experimental data. Moreover, the method is applied to simulate an even more sluggish event: the nucleation of the B2 phase in a strong glass-forming Cu-Zr alloy. The nucleation rate was found to be 8 orders of magnitude smaller than Ni at the same undercooling, which well explains the good glass formability of the alloy. Thus, our work opens a new avenue to study solidification under realistic experimental conditions via atomistic computer simulation.
  • Phase selection in deeply undercooled liquids and devitrified glasses during heating involves complex interplay between the barriers to nucleation and the ability for these nuclei to grow. During the devitrification of glassy alloys, complicated metastable structures often precipitate instead of simpler, more stable compounds. Here, we access this unconventional type of phase selections by investigating an Al-10%Sm system, where a complicated cubic structure first precipitates with a large lattice parameter of 1.4 nm. We not only solve the structure of this "big cubic" phase containing ~140 atoms but establish an explicit interconnection between the structural orderings of the amorphous alloy and the cubic phase, which provides a low-barrier nucleation pathway at low temperatures. The surprising rapid growth of the crystal is attributed to its high tolerance to point defects, which minimize the short-scale atomic rearrangements to form the crystal. Our study suggests a new scenario of devitrification, where phase transformation proceeds initially without partitioning through a complex intermediate crystal phase.
  • The underlying structural order that transcends the liquid, glass and crystalline states is identified using an efficient genetic algorithm (GA). GA identifies the most common energetically favorable packing motif in crystalline structures close to the alloy's Al-10 at.% Sm composition. These motifs are in turn compared to the observed packing motifs in the actual liquid structures using a cluster-alignment method which reveals the average topology. Conventional descriptions of the short-range order, such as Voronoi tessellation, are too rigid in their analysis of the configurational poly-types when describing the chemical and topological ordering during transition from undercooled metallic liquids to crystalline phases or glass. Our approach here brings new insight into describing mesoscopic order-disorder transitions in condensed matter physics.