Arxiv:1707.03044V1 [Cond-Mat.Mtrl-Sci] 10 Jul 2017 Employing Body-Centered Cubic (Bcc) [12–16] and Fcc 1 2 [14, 17–20] Crystal Structures
Total Page:16
File Type:pdf, Size:1020Kb
Simulating complex crystal structures using the phase-field crystal model Eli Alster,1, 2 David Montiel,3 Katsuyo Thornton,3 and Peter W. Voorhees2, ∗ 1Department of Chemical and Biological Engineering, Northwestern University, Evanston, Illinois 60208, USA 2Department of Materials Science and Engineering, Northwestern University, Evanston, Illinois 60208, USA 3Department of Materials Science and Engineering, University of Michigan, Ann Arbor, MI 48109, USA (Dated: July 8, 2018) We introduce a phase-field crystal model that creates an array of complex three- and two- dimensional crystal structures via a numerically tractable three-point correlation function. The three-point correlation function is designed in order to energetically favor the principal interplanar angles of a target crystal structure. This is achieved via an analysis performed by examining the crystal’s structure factor. This approach successfully yields energetically stable simple cubic, dia- mond cubic, simple hexagonal, graphene layers, and CaF2 crystals. To illustrate the ability of the method to yield a particularly complex and technologically important crystal structure, we show how this three-point correlation function method can be used to generate perovskite crystals. Multiscale phenomena in materials are notoriously dif- design interaction potentials between discrete particles ficult to model because they elude conventional tech- such that a given structure is a global energy minimum. niques like molecular dynamics, continuum mechanics, They have also found this task to be non-trivial [21–24]. and phase-field theory. This is a problem because most Traditionally, the single-component PFC free energy interesting phenomena span multiple orders of magnitude functional is expressed as a combination of one-body and in both time and length scales. One promising method two-body interactions. Namely, for multiscale simulation of crystalline materials is the Z 1 2 1 3 1 4 phase-field crystal (PFC) method. The PFC method uti- F [n] = F1[n] + F2[n] = n − n + n dr lizes a free energy functional that is averaged over the V 2 6 12 Z time scale of atomic vibrations but retains patterns as- 1 − n(r)C2 ∗ n dr, (1) sociated with lattice structures as equilibrium states [1]. 2 V It reproduces Read-Shockley grain boundary energies [2] where F is a nondimensionalized free energy, n is a nondi- and has been used to examine graphene grain boundary mensionalized density, F is the ideal free energy term, structure [3], to model step-flow growth from a supersat- 1 F is the energy from two-point interactions, V is the urated vapor [4], and to study electromigration in metal 2 system volume, C2 is an isotropic two-point correlation interconnects [5], among other applications [6]. R 0 0 0 function, and C2 ∗ n ≡ V 0 C2(|r − r |)n(r )dr [1, 7, 8]. A major shortcoming of the PFC method is the limited Since this free energy functional is rotationally and trans- number of crystal structures it can simulate. Progress lationally invariant, it is not a trivial task to design a C2 in the PFC community has largely been made by con- function that produces the desired crystal structure as sidering various phenomenological forms for the Fourier an energy minimum. In fact, symmetry considerations ˆ transform of the two-point correlation function, C2 [7, 8]. suggest that most forms of C2 will only result in a bcc or ˆ For a C2 containing up to three peaks, there exist two- lower-dimensional structure [25]. dimensional (2D) PFC models for all five Bravais lattices In an effort to derive an improved model for graphene, [9] and various chiral phases [10]. In three-dimensions, Seymour et al. added the energy due to three-point cor- ˆ PFC models with this C2 can form simple cubic [8], face- relations [26], centered cubic (fcc) [8], and diamond cubic structures [11]. Consequently, applications of three-dimensional F [n] = F1 + F2 + F3, (2) (3D) PFC models have been dominated by simulations where F and F are the same as in Eq. 1 and arXiv:1707.03044v1 [cond-mat.mtrl-sci] 10 Jul 2017 employing body-centered cubic (bcc) [12–16] and fcc 1 2 [14, 17–20] crystal structures. Obviously, there exist 1 ZZZ F [n] = − n(r)C (r−r0, r−r00)n(r0)n(r00)drdr0dr00. many more crystal structures, and no attempts have been 3 6 3 made to describe a crystal as complex as perovskite, a (3) deficit that we address in this Letter. Although in general the calculation of F3 is of compu- 3 The challenge of producing complex crystal structures tational complexity O(N ), if the C3 function is of the in PFC models is similar to that faced by the self- form assembly community. Their goal is to solve the so- X (i) (i) C3(r1, r2) = C3 (r1)C3 (r2), (4) called “inverse” statistical mechanics problem: how to i 2 0 00 where r1 ≡ r − r and r2 ≡ r − r , the free energy and boundary energies as a function of misorientation, which evolution equations simplify into a number of convolu- would be impossible otherwise. tions of O(N log N) computational complexity that are This ansatz was chosen not only because it is rotation- easily computed via the fast-Fourier transform. Seymour ally invariant but also because the Legendre polynomials et al. proposed a single length scale, 2D real space form are both separable and form a complete orthogonal set. (i) for C3 that yields equilibrium states with a specified More explicitly, the separability of the Legendre polyno- bond angle when the bond angle, θ, satisfies the relation mials means that l ◦ X 4π 360 mod θ = 0, (5) P (kˆ · kˆ ) = Y (kˆ )Y (kˆ ), (7) l 1 2 2l + 1 lm 1 lm 2 for example 60◦, 90◦, and 120◦ [26]. Although this m=−l yielded an improved 2D graphene model [3], the model where Ylm are the normalized real spherical harmonics was not flexible enough to generate any new crystal struc- [34]. Thus, Cˆ3 can be written as a sum of products of tures in either two or three dimensions [26, 27]. two-point correlation functions, i.e., In this Letter, we introduce a form for C3 that sta- l l bilizes angles between specified crystallographic planes Xmax X Cˆ (k , k ) = α (−1)l Cˆ(lm)(k , kˆ )Cˆ(lm)(k , kˆ ), and can include multiple length scales, multiple preferred 3 1 2 l 1 1 2 2 l=0 m=−l angles, and angles not restricted by Eq. 5. Then, we (8) will discuss how to choose parameters in order to yield where energy-minimizing single-component crystal structures. r We show that this method can produce a wide array of 4π Cˆ(lm)(k, kˆ) ≡ (−i)l βR(k)Y (kˆ). (9) energy-minimizing crystal structures, from simple cubic 2l + 1 lm (e.g., Po [28]) and diamond cubic (e.g., C-diamond, Si, l α-tin [29]), to graphene layers and disordered CaF (e.g., The factor of (−1) in Eq. 8 was introduced to cancel the 2 l (lm) 0 phase factors, (−i) , in Eq. 9 so that C is real (see the structure of the θc phase of Al2Cu, which is com- mercially very important for strengthening in aluminum Supplementary Materials for details). alloys [30]). Additionally, it produces an unnamed crystal Substituting the inverse Fourier transform of Eq. 8 into Eq. 3 results in structure corresponding to the X atoms in ABX3 per- ovskite (a structure we will call X3), which is necessary l l 2 1 Xmax X Z for modeling perovskite. F = − α (−1)l n(r) C(lm) ∗ n dr (10) 3 6 l Finally, as a capstone demonstration of the method, l=0 m=−l we combine the X3 and simple cubic models to generate a perovskite crystal structure. Such compounds are of and great interest due to applications ranging from high ef- lmax l l δF3 X αl(−1) X ficiency solar cells [31, 32] to light-emitting diodes [33]. = − (C(lm) ∗ n)2 δn 6 However, the development of perovskite microstructure is l=0 m=−l governed by phenomena that occur on a diffusional time l (lm) (lm) scale, not the nanoseconds afforded by molecular dynam- + 2(−1) C ∗ [n(C ∗ n)] , (11) ics. Consequently, this model will provide a new avenue to investigate multiscale phenomena in these important since C(lm)(−r) = (−1)lC(lm)(r) by the parity property materials. of real spherical harmonics (note the similarity between The ansatz we use for C3 is these expressions and Eq. 10 and Eq. 42 respectively from [26]). l Xmax Since the Legendre polynomials form a complete or- Cˆ (k , k ) = β2R(k )R(k ) α P (kˆ · kˆ ), (6) 3 1 2 1 2 l l 1 2 thogonal set, if l=0 lmax ˆ X where C3 is the Fourier transform of the three-point cor- B(x) ≡ αlPl(x), (12) relation function, ki = |ki|, kˆi = ki/|ki|, β is an inter- l=0 action strength parameter, R(k) is a real radial function, then each αl is given by Pl are the Legendre polynomials, and αl are constant co- efficients. As is explained later, αl should be interpreted 2l + 1 Z 1 as determining the preferred interplanar angles. Because αl = B(x)Pl(x)dx. (13) 2 −1 Cˆ3(Rk1, Rk2) = Cˆ3(k1, k2) for any rotation matrix R, Cˆ3 is rotationally invariant and so is C3. By keeping This is convenient because it implies that the angular the free energy rotationally invariant, it is possible to portion of Eq. 6 can represent any function through a study phenomena such as solid-liquid interfaces and grain series of Legendre polynomials.