Numerical Modelling of Domain Structure Annihilation in a Nonlinear Hyperbolic Phase-Field Model

Context and relevance. Modelling the dynamics of the order parameter in nonlinear media with inertial effects is of interest for condensed matter physics and the theory of phase transitions. Hyperbolic generalisations of the Landau–Khalatnikov equation make it possible to account for wave mechanisms of energy transport, which fundamentally distinguishes them from classical relaxation models. Their numerical investigation requires a stable and consistent second-order accurate algorithm capable of correctly reproducing both nonlinear relaxation and inertially induced wave processes. Objective. To investigate the dynamics of the order parameter in a nonlinear medium described by an inertial equation of the Landau–Khalatnikov–Tani type, with application to the process of spontaneous annihilation of a narrow-localised domain structure in a two-dimensional formulation. Hypothesis. It is assumed that inclusion of the inertial term leads to a qualitative modification of the domain boundary evolution, namely to redistribution of gradient energy into dynamical oscillations of the order parameter field, and that the use of an implicit finite-difference scheme with iterative treatment of the nonlinearity ensures correct reproduction of all stages of annihilation. Methods and materials. An implicit Crank–Nicolson finite-difference scheme with second-order temporal accuracy is employed for the numerical solution. The nonlinear term is approximated using the Newton method with iterative refinement at each time step. Correctness of the algorithm implementation is confirmed by comparison of numerical and analytical solutions for a test problem with a known exact solution. Results. Numerical experiments demonstrate successive stages of domain structure evolution: convergence of the interfaces, their disappearance, and the formation of propagating wave disturbances. It is shown that the presence of the inertial term leads to redistribution of part of the gradient energy into dynamical oscillations of the order parameter field. Conclusions. The obtained results confirm the correctness of the implemented numerical approach and demonstrate the significant influence of inertial effects on the dynamics of localised structures. The hyperbolic formulation allows description of wave processes absent in parabolic models, which is important for analysing the energy balance and the evolution of domain boundaries.

1. Landau, L.D., Lifshitz, E.M. (1980). Statistical Physics. Part 1. Oxford: Pergamon Press. 584 p. (In Eng.).
2. Lines, M.E., Glass, A.M. (1977). Principles and Applications of Ferroelectrics and Related Materials. Oxford: Clarendon Press. 610 p. (In Eng.).
3. Onuki, A. (2002). Phase Transition Dynamics. Cambridge: Cambridge University Press. 308 p. (In Eng.). https://doi.org/10.1017/CBO9780511534874
4. Chandra, P., Littlewood, P.B. (2007). Physics of Ferroelectrics. In: K.M. Rabe, C.H. Ahn, J.-F. Scott (Eds.), Physics of Ferroelectrics: A Modern Perspective (pp. 23—56). Berlin: Springer. (In Eng.). https://doi.org/10.1007/978-3-540-34591-6_3
5. Khalatnikov, I.M. (1965). An Introduction to the Theory of Superfluidity. New York: W.A. Benjamin. 266 p. (In Eng.).
6. Cahn, J.W., Hilliard, J.E. (1958). Free energy of a nonuniform system. I. Interfacial free energy. Journal of Chemical Physics, 28, 258—267. (In Eng.). https://doi.org/10.1063/1.1744102
7. Bray, A.J. (1994). Theory of phase-ordering kinetics. Advances in Physics, 43, 357—459. (In Eng.). https://doi.org/10.1080/00018739400101505
8. Chen, L.-Q. (2002). Phase-field models for microstructure evolution. Annual Review of Materials Research, 32, 113—140. (In Eng.). https://doi.org/10.1146/annurev.matsci.32.112001.132041
9. Frey, R., et al. (2023). Accelerated search for new ferroelectric materials using machine learning and high throughput DFT. Physical Review Research, 5, 023122. (In Eng.). https://doi.org/10.1103/PhysRevResearch.5.023122
10. Allen, S.M., Cahn, J.W. (1979). A microscopic theory for antiphase boundary motion and its application to antiphase domain coarsening. Acta Metallurgica, 27, 1085—1095. (In Eng.). https://doi.org/10.1016/0001-6160(79)90196-2
11. Chen, L.-Q. (2008). Phase-field method for microstructure evolution. JournaloftheAmericanCeramicSociety, 91(6), 1835—1844. (InEng.). https://doi.org/10.1111/j.1551-2916.2008.02413.x
12. Wang, B., Xia, R., Fan, H., Woo, C.H. (2003). Dynamic process of domain switching in ferroelectric films. Journal of Applied Physics, 94, 3384—3391. (In Eng.). https://doi.org/10.1063/1.1599958
13. Wang, Y.L., Chen, L.-Q., Nagarajan, V. (2013). Phase-field modeling of domain wall dynamics in ferroelectrics with inertial effects. Physical Review B, 88, 054404. (In Eng.). https://doi.org/10.1103/PhysRevB.88.054404
14. Tani, S. (1972). Relaxation of order parameter in ferroelectrics. Progress of Theoretical Physics, 47, 1740—1752. (In Eng.).
15. Cross, L.E. (1987). Ferroelectric relaxors. Ferroelectrics, 69, 1—10. (In Eng.). https://doi.org/10.1080/00150198708016945
16. Zhong, W., Vanderbilt, D., Rabe, K.M. (1994). Phase transitions in BaTiO₃ from first principles. Physical Review Letters, 73, 1861—1864. (In Eng.). https://doi.org/10.1103/PhysRevLett.73.1861
17. Bhattacharya, K. (2003). Microstructure of Martensite: Why it Forms and How it Gives Rise to the Shape-Memory Effect. Oxford: Oxford University Press. 268 p. (In Eng.).
18. Scott, J.F. (2000). Ferroelectric Memories. Heidelberg: Springer Science & Business Media. 248 p. (In Eng.). https://doi.org/10.1007/978-3-662-04307-3
19. Temam, R. (1997). Infinite-Dimensional Dynamical Systems in Mechanics and Physics. New York: Springer. 428 p. (In Eng.). https://doi.org/10.1007/978-1-4612-0645-3
20. Rabe, K.M., Ahn, C.H., Scott, J.-F. (Eds.). (2007). Physics of Ferroelectrics: A Modern Perspective. Berlin: Springer. (In Eng.).
21. Moroz, L., Maslovskaya, A. (2021). Fractional differential model of domain boundary kinetics in ferroelectrics: A computational approach. AIP Conference Proceedings, 2328(1), 020011. (In Eng.). https://doi.org/10.1063/5.0042140
22. Zhang, J., Yang, X. (2022). Efficient Fully Discrete Finite-Element Numerical Scheme with Second-Order Temporal Accuracy for the Phase-Field Crystal Model. Mathematics, 10(1), 155. (In Eng.). https://doi.org/10.3390/math10010155
23. Formalev, V.F., Reyznikov, D.L. (2006). Numerical Methods. Moscow: FIZMATLIT, 400 p. (In Russ.)
24. Alikin, Yu.M., Turygin, A.P., Kosobokov, M.S., Alikin, D.O., Shur, V.Ya. (2022). Decay of domains created by local switching on non-polar cut of MgO doped LiNbO₃ single crystals. Ferroelectrics, 592(1), 12—18. (In Eng.). https://doi.org/10.1080/00150193.2022.2052240
25. Su, Y., Landis, C.M. (2007). Continuum thermodynamics of ferroelectric domain evolution: Theory, finite element implementation, and application to domain wall pinning. Journal of the Mechanics and Physics of Solids, 55(2), 280—305. (In Eng.). https://doi.org/10.1016/j.jmps.2006.07.006
26. Richman, M.S., Rulis, P., Caruso, A.N. (2017). Ferroelectric system dynamics simulated by a second-order Landau model. Journal of Applied Physics, 122(9), 094101. (In Eng.). https://doi.org/10.1063/1.5000139