Free energy of the bcc-liquid interface and the Wulff shape as predicted by the Phase-Field Crystal model
The Euler-Lagrange equation of the phase-field crystal (PFC) model has been solved under appropriate boundary conditions to obtain the equilibrium free energy of the body centered cubic crystal-liquid interface for 18 orientations at various reduced temperatures in the range $ε\in\left[0,0.5\right]$. While the maximum free energy corresponds to the $\left\{ 100\right\} $ orientation for all $ε$ values, the minimum is realized by the $\left\{ 111\right\} $ direction for small $ε\,(<0.13)$, and by the $\left\{ 211\right\} $ orientation for higher $ε$. The predicted dependence on the reduced temperature is consistent with the respective mean field critical exponent. The results are fitted with an eight-term Kubic harmonic series, and are used to create stereographic plots displaying the anisotropy of the interface free energy. We have also derived the corresponding Wulff shapes that vary with increasing $ε$ from sphere to a polyhedral form that differs from the rhombo-dodecahedron obtained previously by growing a bcc seed until reaching equilibrium with the remaining liquid.