Simulation of Harmonic Oscillators on the Lattice

[EN] This work deals with the simulation of a two¿dimensional ideal lattice having simple tetragonal geometry. The harmonic character of the oscillators give rise to a system of second¿order linear differential equations, which can be recast into matrix form. The explicit solutions which govern the dynamics of this system can be expressed in terms of matrix trigonometric functions. For the derivation we employ the Lagrangian formalism to determine the correct solutions, which extremize the underlying action of the system. In the numerical evaluation we develop diverse state¿of¿the¿art algorithms which efficiently tackle equations with matrix sine and cosine functions. For this purpose, we introduce two special series related to trigonometric functions. They provide approximate solutions of the system through a suitable combination. For the final computation an algorithm based on Taylor expansion with forward and backward error analysis for computing those series had to be devised. We also implement several MATLAB programs which simulate and visualize the two¿dimensional lattice and check its energy conservation. This work has been supported by the Spanish Ministerio de Economia y Competitividad, the European Regional Development Fund (ERDF) under grant TIN2017-89314-P, and the Programa de Apoyo a la Investigacion y Desarrollo 2018 (PAID-06-18) of the Universitat Politecnica de Valencia under grant SP20180016. Tung, MM.; Ibáñez González, JJ.; Defez Candel, E.; Sastre, J. (2020). Simulation of Harmonic Oscillators on the Lattice. Mathematical Methods in the Applied Sciences. 43(14):8237-8252. https://doi.org/10.1002/mma.6510 Dehghan, M., & Hajarian, M. (2009). Determination of a matrix function using the divided difference method of Newton and the interpolation technique of Hermite. Journal of Computational and Applied Mathematics, 231(1), 67-81. doi:10.1016/j.cam.2009.01.021 Dehghan, M., & Hajarian, M. (2010). Computing matrix functions using mixed interpolation methods. Mathematical and Computer Modelling, 52(5-6), 826-836. doi:10.1016/j.mcm.2010.05.013 Kazem, S., & Dehghan, M. (2017). Application of finite difference method of lines on the heat equation. Numerical Methods for Partial Differential Equations, 34(2), 626-660. doi:10.1002/num.22218 Kazem, S., & Dehghan, M. (2018). Semi-analytical solution for time-fractional diffusion equation based on finite difference method of lines (MOL). Engineering with Computers, 35(1), 229-241. doi:10.1007/s00366-018-