A python implementation of Phase Field Crystal using pseudospectral method
PythonGPL-3.0
PyPFC
A python implementation of Phase Field Crystal solved with semi-implicit pseudospectral method
For the PFC, the conserved order parameter is a relative density $n(\boldsymbol{r},t)$ which evolves following a conservation equation as
$$\frac{\partial n}{\partial t} = \nabla \cdot \left[ M \nabla \left( \frac{\delta F[n]}{\delta n} \right) \right]$$
with M being a mobility, $\delta F[n]/\delta n$ is the functional derivative of free-energy. The free-energy functional 1 given by
$$F[n(\boldsymbol{r},t)] = \int_V \left[ \frac{1}{2} n \left(1+\nabla^2\right)^2 n + \frac{1}{4} n^2(2r+n^2) \right]\ d{\boldsymbol{r}}$$
with $r$ being a constant proportional to the temperature deviation from the melting point 2.
Example
The following figure is a result for the system with $M=1.0$, $r=-0.25$, $L_x=L_y=16\pi$, $N_x=N_y=2^8$, and $dt=0.1$. The initial condition is given by a normal distribution described by $n(x,y,t=0) = -0.285(1.0 - 0.02\mathcal{N}(0,1))$. The system evolves until $T=1500$.
