Dynamics of stress propagation in anharmonic crystals: MD simulations

Zbigniew Kozioł
2024-04-27
Abstract:Anharmonic inter-atomic potential x to power of n, n>1, has been used in molecular dynamics (MD) simulations of stress dynamics of FCC oriented crystal. The model of the chain of masses and springs is found as a convenient and accurate description of simulation results, with masses representing the crystallographic planes. The dynamics of oscillations of two planes is found analytically to be given by Euler beta functions, and its scaling with non-linearity parameter and amplitude of oscillations, or applied shear pressure is discussed on examples of time dependencies of displacements, velocities, and forces acting on masses (planes). The dynamics of stress penetration from the surface of the sample with multiply-planes (an anharmonic crystal) towards its interior is confirmed to be given exactly as a series of Bessel functions, when n=2 (Schrodinger and Pater solutions). When n<>2 the stress dynamics (wave propagation) in bulk material remains qualitatively of the same nature as in the linear case. In particular, results suggest that the quasi-linear relationship between frequency and wave number is preserved. Speed of the transverse sound component, dependent on sound wave amplitude, is found to be a strongly decreasing function of n. The results are useful in the analysis of any MD simulations under pressure, as they help to understand the dynamics of pressure retarded effects, as well help design the proper methodology of performing MD simulations, such as studies of dynamics of dislocations.
Mesoscale and Nanoscale Physics
What problem does this paper attempt to address?