Arxiv:1207.0257V1 [Cond-Mat.Soft] 2 Jul 2012
Total Page:16
File Type:pdf, Size:1020Kb
November 2, 2018 17:7 Advances in Physics aip Advances in Physics Vol. 00, No. 00, XXX 2012, 1{95 RESEARCH ARTICLE Phase-field-crystal models for condensed matter dynamics on atomic length and diffusive time scales: an overview Heike Emmericha∗, Hartmut L¨owenb∗, Raphael Wittkowskib, Thomas Gruhna, Gyula I. T´othc, Gy¨orgyTegzec, and L´aszl´oGr´an´asyc;d∗ aLehrstuhl f¨urMaterial- und Prozesssimulation, Universit¨atBayreuth, D-95440 Bayreuth, Germany, bInstitut f¨urTheoretische Physik II, Weiche Materie, Heinrich-Heine-Universit¨at D¨usseldorf,D-40225 D¨usseldorf,Germany, cInstitute for Solid State Physics and Optics, Wigner Research Centre for Physics, P.O. Box 49, H-1525 Budapest, Hungary, dBCAST, Brunel University, Uxbridge, Middlesex, UB8 3PH, U.K. (Received 00 Month 200x; final version received 00 Month 200x) Here, we review the basic concepts and applications of the phase-field-crystal (PFC) method, which is one of the latest simulation methodologies in materials science for problems, where atomic- and microscales are tightly coupled. The PFC method operates on atomic length and diffusive time scales, and thus constitutes a computationally efficient alternative to molecu- lar simulation methods. Its intense development in materials science started fairly recently following the work by Elder et al. [Phys. Rev. Lett. 88 (2002), p. 245701]. Since these initial studies, dynamical density functional theory and thermodynamic concepts have been linked to the PFC approach to serve as further theoretical fundaments for the latter. In this review, we summarize these methodological development steps as well as the most important appli- cations of the PFC method with a special focus on the interaction of development steps taken in hard and soft matter physics, respectively. Doing so, we hope to present today's state of the art in PFC modelling as well as the potential, which might still arise from this method in physics and materials science in the nearby future. Keywords: phase-field-crystal (PFC) models, static and dynamical density functional theory (DFT and DDFT), condensed matter dynamics of liquid crystals and polymers, nucleation and pattern formation, simulations in materials science, colloidal crystal growth and growth anisotropy arXiv:1207.0257v1 [cond-mat.soft] 2 Jul 2012 ∗Corresponding authors. Emails: [email protected], [email protected], and [email protected] ISSN: 0001-8732 print/ISSN 1460-6976 online c 2012 Taylor & Francis DOI: 10.1080/0001873YYxxxxxxxx http://www.informaworld.com November 2, 2018 17:7 Advances in Physics aip 2 Taylor & Francis and I. T. Consultant Table of contents 4. Phase-field-crystal models applied to nucleation and pattern formation 1. Introduction in metals 2. From density functional theory 4.1. Properties of nuclei from to phase-field-crystal models extremum principles 2.1. Density functional theory 4.1.1. Homogeneous nucleation 2.2. Dynamical density functional theory 4.1.2. Heterogeneous nucleation 2.2.1. Basic equations 4.2. Pattern formation 2.2.2. Brownian dynamics: Langevin 4.2.1. PFC modelling of surface and Smoluchowski picture patterns 2.2.3. Derivation of DDFT 4.2.2. Pattern formation in binary 2.2.4. Application of DDFT to solidification colloidal crystal growth 4.3. Phenomena in the solid state 2.3. Derivation of the PFC model 4.3.1. Dislocation dynamics and grain for isotropic particles from DFT boundary melting 2.3.1. Free-energy functional 4.3.2. Crack formation and 2.3.2. Dynamical equations propagation 2.3.3. Colloidal crystal growth: 4.3.3. Strain-induced morphologies DDFT versus PFC modelling 4.3.4. Kirkendall effect 3. Phase-field-crystal modelling 4.3.5. Density/solute trapping in condensed matter physics 4.3.6. Vacancy/atom transport in the 3.1. The original PFC model and its VPFC model generalisations 5. Phase-field-crystal modelling 3.1.1. Single-component PFC models in soft matter physics 3.1.2. Binary PFC models 5.1. Application to colloids 3.1.3. PFC models for liquid crystals 5.1.1. Nucleation in colloidal 3.1.4. Numerical methods crystal aggregation 3.1.5. Coarse-graining the PFC 5.1.2. Pattern formation in colloidal models crystal aggregation 3.2. Phase diagrams the PFC models 5.1.3. Colloid patterning realize 5.2. Application to polymers 3.2.1. Phase diagram of single-compo- 5.2.1. Self-consistent field theory nent and binary systems for polymer systems 3.2.2. Phase diagram of two- 5.2.2. Phase-field methods dimensional liquid crystals for polymer systems 3.3. Anisotropies in the PFC models 5.2.3. Comparison of the methods 3.3.1. Free energy of the liquid-solid 5.3. Application to liquid crystals interface 6. Summary and outlook 3.3.2. Growth anisotropy Appendix 3.4. Glass formation Acknowledgements 3.5. Phase-field-crystal modelling of foams References 3.6. Coupling to hydrodynamics List of abbreviations November 2, 2018 17:7 Advances in Physics aip Advances in Physics 3 Figure 1. Pattern formation on microscopic to cosmic length scales. From left to right: multiple spiralling nanoscale terraces starting from a central heterogeneity. (Reproduced with permission from C. Klemenz, Hollow cores and step-bunching effects on (001) YBCO surfaces grown by liquid-phase epitaxy, J. Cryst. Growth 187 (1998), 221-227, no. 2, DOI: 10.1016/S0022-0248(97)00866-X c 1998 by Elsevier.) Cellular slime mould self-organized into a five-arm spiral structure. (Reproduced with permission from B. Vasiev, F. Siegert, and C. Weijer, Multiarmed spirals in excitable media, Phys. Rev. Lett. 78 (1997), 2489-2492, no. 12, DOI: 10.1103/PhysRevLett.78.2489 c 1997 by the American Physical Society.) Messier 100, a multi- arm spiral galaxy in the Virgo Supercluster, 60 million light-years from earth. (Credit: ESO/IDA/Danish 1.5m/R. Gendler, J.-E. Ovaldsen, C. C. Th¨one,and C. F´eron.) 1. Introduction Pattern formation has been observed in complex systems from microscopic to cos- mic scales (for examples, see figure 1), a phenomenon that has been exciting the fantasy of humanity for a long time. Nonequilibrium systems in physics, chem- istry, biology, mathematics, cosmology, and other fields produce an amazingly rich and visually fascinating variety of spatiotemporal behaviour. Experiments and sim- ulations show that many of such systems { reacting chemicals, bacteria colonies, granular matter, plasmas { often display analogous dynamical behaviour. The wish to find the origin of the common behaviour has been driving the efforts for find- ing unifying schemes that allow the assigning of many of these processes into a few universality classes. Pattern formation and the associated nonlinear dynamics have received a continuous attention of the statistical physics community over the past decades. Reviews of the advances made in different directions are available in the literature and range from early works on critical dynamics [1] via phase- separation [2] and pattern formation in nonequilibrium systems [3, 4] to recent detailed treatments of the field in books (e. g., references [5{7]). In particular, Seul and Andelman [4] described pattern formation on the mesoscale as manifestation of modulated structures. Within this approach, the modulated phases are stabilized by competing attractive and repulsive interac- tions, which favour inhomogeneities characterized by a certain modulation length scale. The modulations are described by a single scalar order parameter. As out- lined in reference [4], the idea of Seul and Andelman can be applied to a large variety of systems ranging from Langmuir films over semiconductor surfaces and magnet garnets to polyelectrolyte solutions. Furthermore, the pioneering theories of spontaneous domain formation in magnetic materials and in the intermediate state of type I superconductors has been reinterpreted within this framework. In the past decade, special attention has been paid to a similar model, whose mathematical formulation has been laid down decades earlier to address hydro- dynamic instabilities [10] and to describe the transition to the antiferromagnetic state in liquid 3He or to a non-uniform state in cholesteric liquid crystals [11], whereas recently it has been employed for the modelling of crystallization in un- dercooled liquids on the atomic scale [12]. This approach is known to the materials science community as the phase-field-crystal (PFC) model [12], and has proven to be an amazingly efficient tool for addressing crystalline self-organization/pattern formation on the atomistic scale. November 2, 2018 17:7 Advances in Physics aip 4 Taylor & Francis and I. T. Consultant The PFC approach attracts attention owing to a unique situation: the crystal- lization of liquids is traditionally addressed on this scale by the density functional theory (DFT) [13{15], whose best developed non-perturbative version, known as the fundamental-measure theory (FMT) [16], leads to unprecedented accuracy for such properties as the liquid-solid interfacial free energy [17, 18] or the nucleation barrier [17]. However, handling of large systems is hampered by the complexity of such models. In turn, the PFC model, being a simplistic DFT itself, incor- porates most of the essential physics required to handle freezing: it is atomistic, anisotropies and elasticity are automatically there, the system may choose from a variety of periodic states [such as body-centred cubic (bcc), face-centred cubic (fcc), and hexagonal close packed (hcp)] besides the homogeneous fluid, etc. The free-energy functional is fairly simple having the well-known Swift-Hohenberg (SH) form Z ~ ~4 F~ = d~r − β + k2 + r22 ~ + ; (1) 2 0 ~r 4 where ~ is the reduced particle density and β a reduced temperature, while k0 is the absolute value of the wave number vector the system prefers. (In equation (1) all quantities are dimensionless.) This together with the assumption of overdamped conservative (diffusive) dynamics (a major deviation from the non-conservative dy- namics of the SH model) leads to a relatively simple equation of motion (EOM) that, in turn, allows the handling of a few times 107 atoms on the diffusive time scale.