Computational Materials Science Group

Frigyes Podmaniczky¹, Gyula Tóth², György Tegze¹, László Gránásy^1,3

¹Institute for Solid State Physics and Optics, Wigner Research Centre for Physics, P.O. Box 49, Budapest H-1525, Hungary
²Department of Mathematical Sciences, Loughborough University, Loughborough, Leicestershire, LE11 3TU, U.K.
³BCAST, Brunel University, Uxbridge, Middlesex, UB8 3PH, United Kingdom

Crystallization of supersaturated liquids usually starts by epitaxial growth or by heterogeneous nucleation on foreign surfaces. Herein, we review recent advances made in modeling heteroepitaxy and heterogeneous nucleation on flat/modulated surfaces and nanoparticles within the framework of a simple dynamical density functional theory, known as the phase-field crystal model. It will be shown that the contact angle and the nucleation barrier are nonmonotonous functions of the lattice mismatch between the substrate and the crystalline phase. In continuous cooling studies for substrates with lattice mismatch, we recover qualitatively the Matthews–Blakeslee mechanism of stress release via the misfit dislocations. The simulations performed for particle-induced freezing will be confronted with recent analytical results, exploring thus the validity range of the latter. It will be demonstrated that time-dependent studies are essential, as investigations based on equilibrium properties often cannot identify the preferred nucleation pathways. Modeling of these phenomena is essential for designing materials on the basis of controlled nucleation and/or nano-patterning.

László Rátkai¹, Attila Szállás¹, Tamás Pusztai¹, Tetsuo Mohri, László Gránásy^1,2

¹Institute for Solid State Physics and Optics, Wigner Research Centre for Physics, P.O. Box 49, Budapest H-1525, Hungary
²BCAST, Brunel University, Uxbridge, Middlesex, UB8 3PH, United Kingdom

Extending previous work [Pusztai et al., Phys. Rev. E 87, 032401 (2013)], we have studied the formation of eutectic dendrites in a model ternary system within the framework of the phase-field theory. We have mapped out the domain in which two-phase dendritic structures grow. With increasing pulling velocity, the following sequence of growth morphologies is observed: flat front lamellae → eutectic colonies → eutectic dendrites → dendrites with target pattern → partitionless dendrites → partitionless flat front. We confirm that the two-phase and one-phase dendrites have similar forms and display a similar scaling of the dendrite tip radius with the interface free energy. It is also found that the possible eutectic patterns include the target pattern, and single- and multiarm spirals, of which the thermal fluctuations choose. The most probable number of spiral arms increases with increasing tip radius and with decreasing kinetic anisotropy. Our numerical simulations confirm that in agreement with the assumptions of a recent analysis of two-phase dendrites [Akamatsu et al., Phys. Rev. Lett. 112, 105502 (2014)], the Jackson-Hunt scaling of the eutectic wavelength with pulling velocity is obeyed in the parameter domain explored, and that the natural eutectic wavelength is proportional to the tip radius of the two-phase dendrites. Finally, we find that it is very difficult/virtually impossible to form spiraling two-phase dendrites without anisotropy, an observation that seems to contradict the expectations of Akamatsu et al. Yet, it cannot be excluded that in isotropic systems, two-phase dendrites are rare events difficult to observe in simulations.

Tamás Pusztai¹, László Rátkai¹, Attila Szállás¹, László Gránásy^1,2

¹Institute for Solid State Physics and Optics, Wigner Research Centre for Physics, P.O. Box 49, Budapest H-1525, Hungary
²BCAST, Brunel University, Uxbridge, Middlesex, UB8 3PH, United Kingdom

The quantitative phase-field approach has been adapted to model solidification in the presence of Metal Matrix Nanocomposites (MMNCs) in a single-component liquid. Nanoparticles of fixedsize and shape are represented by additional fields. The corresponding equations of motion are assumed to ensure relaxation dynamics, and can be supplemented by random forces (realizing Brownian motion) or external fields. The nanoparticles are characterized by two model parameters: their mobility and the contact angle they realize with the solid-liquid interface. We investigate the question how grain size distribution can be influenced by heterogeneous nucleation on the nanoparticles and by the front-particle interaction. We explore, furthermore, how materials and process parameters, such as temperature, density and size/shape distribution of the nanoparticles, influence microstructure evolution.

László Gránásy^1,2, László Rátkai¹, Attila Szállás¹, Bálint Korbuly¹, Gyula Tóth³, László Környei⁴, Tamás Pusztai¹

¹Institute for Solid State Physics and Optics, Wigner Research Centre for Physics, P.O. Box 49, Budapest H-1525, Hungary
²BCAST, Brunel University, Uxbridge, Middlesex, UB8 3PH, United Kingdom
³Department of Mathematical Sciences, Loughborough University, Loughborough, Leicestershire, LE11 3TU, U.K.
⁴Department of Mathematics and Computational Sciences, Széchenyi István University, Győr 9026, Hungary

Advances in the orientation-field-based phase-field (PF) models made in the past are reviewed. The models applied incorporate homogeneous and heterogeneous nucleation of growth centers and several mechanisms to form new grains at the perimeter of growing crystals, a phenomenon termed growth front nucleation. Examples for PF modeling of such complex polycrystalline structures are shown as impinging symmetric dendrites, polycrystalline growth forms (ranging from disordered dendrites to spherulitic patterns), and various eutectic structures, including spiraling two-phase dendrites. Simulations exploring possible control of solidification patterns in thin films via external fields, confined geometry, particle additives, scratching/piercing the films, etc. are also displayed. Advantages, problems, and possible solutions associated with quantitative PF simulations are discussed briefly.