Abstract:The phase-field crystal model (PFC) describes crystal structures at diffusive timescales through a periodic order parameter representing the atomic density. One of its main features is that it naturally incorporates elastic and plastic deformation. To correctly interpret numerical simulation results or devise extensions related to the elasticity description, it is important to have direct access to the elastic field. In this work, we discuss its evaluation in classical PFC models based on the Swift-Hohenberg energy functional. We consider approaches where the stress field can be derived from the microscopic density field (i.e., the order parameter), and a simple novel numerical routine is proposed. By numerical simulations, we demonstrate that it overcomes some limitations of currently used methods. Moreover, we shed light on the elasticity description conveyed by classical PFC models, characterizing a residual stress effect present at equilibrium. We show explicitly and discuss the evaluation of the elastic fields in prototypical representative cases involving an elastic inclusion, a grain boundary, and dislocations.
What problem does this paper attempt to address?
This paper aims to solve the problem of how to accurately evaluate the elastic field in the Phase - Field Crystal (PFC) model. Specifically, the paper focuses on how to extract the elastic field from the microscopic density field (i.e., the order parameter in the PFC model) and proposes a new numerical method to achieve this. This method can overcome some limitations of existing methods, especially when dealing with complex structures such as inclusions, grain boundaries and dislocations.
### Main research contents:
1. **Background introduction**:
- **PFC model**: Describes the behavior of crystal structures on the diffusion time scale and represents the atomic density through a periodic order parameter.
- **Elastic and plastic deformation**: The PFC model naturally includes the description of elastic and plastic deformation.
- **Swift - Hohenberg energy functional**: The basis of the PFC model, used to describe the free energy.
2. **Elastic field evaluation**:
- **Derivation of stress field**: A method for deriving the stress field from the microscopic density field.
- **New numerical method**: A new numerical method is proposed to more accurately evaluate the elastic field by filtering high - order harmonics in the density field.
3. **Benchmark tests**:
- **Eshelby inclusion problem**: To verify the effectiveness of the new method in dealing with inclusion problems.
- **Symmetric tilted Σ5 grain boundary**: To study the stress distribution at the grain boundary.
- **Rotating spherical inclusion**: To simulate the dislocation arrangement between the inclusion and the matrix and its evolution.
### Key findings:
- **Residual stress**: In the classical PFC model, residual stress may occur due to the influence of periodic boundary conditions. The paper effectively eliminates this residual stress by filtering high - order harmonics.
- **Stress field matching**: The stress field obtained by the new method is highly consistent with the analytical solution and the results of other methods, especially when dealing with complex structures.
- **Applicability**: The proposed numerical method is not only applicable to the PFC model, but can also be applied to other models based on the microscopic density field.
### Conclusion:
By proposing a new numerical method, the paper successfully solves the problem of the accuracy of elastic field evaluation in the PFC model. This method performs well in dealing with complex structures such as inclusions, grain boundaries and dislocations, providing a powerful tool for future materials science research.
### Formula examples:
- **Free energy functional**:
\[
F[\psi]=\int_{\Omega} f(\psi) \, dr, \quad f(\psi)=\lambda-\frac{\kappa}{2} \psi^{2}-\frac{\delta}{6} \psi^{3}+\frac{\psi^{4}}{12}+\frac{\kappa}{2} \psi L^{2} \psi
\]
- **Stress tensor**:
\[
\sigma_{ij}=\frac{\partial \delta F}{\partial B_{i} u_{j}}
\]
For triangular lattices:
\[
\sigma_{ij}=\kappa\left(\frac{2\beta^{2} \psi_{,j}+\psi_{,k} k_{i}+\psi \beta^{2}_{,i}}{\psi_{,j}-\psi_{,k} k \psi_{,ij}}\right)
\]
For tetragonal lattices:
\[
\sigma_{ij}=-2\kappa(2\psi + 3\psi_{,kk}+\psi_{,kk}(p)(3\psi_{,ij}+2\psi_{,rr})_{ij})
\]
These formulas and methods provide theoretical basis and technical support for understanding and evaluating the elastic field in the PFC model.