Title: PULSEDYN - A dynamical simulation tool for studying strongly nonlinear chains,
Abstract: We introduce PULSEDYN, a particle dynamics program in $C++$, to solve
many-body nonlinear systems in one dimension. PULSEDYN is designed to make
computing accessible to non-specialists in the field of nonlinear dynamics of
many-body systems and to ensure transparency and easy benchmarking of numerical
results for an integrable model (Toda chain) and three non-integrable models
(Fermi-Pasta-Ulam-Tsingou, Morse and Lennard-Jones). To achieve the latter, we
have made our code open source and free to distribute. We examine (i) soliton
propagation and two-soliton collision in the Toda system, (ii) the recurrence
phenomenon in the Fermi-Pasta-Ulam-Tsingou system and the decay of a single
localized nonlinear excitation in the same system through quasi-equilibrium to
an equipartitioned state, and SW propagation in chains with (iii) Morse and
(iv) Lennard-Jones potentials. We recover well known results from theory and
other numerical results in the literature. We have obtained these results by
setting up a parameter file interface which allows the code to be used as a
black box. Therefore, we anticipate that the code would prove useful to
students and non-specialists. At the same time, PULSEDYN provides
scientifically accurate simulations thus making the study of rich dynamical
processes broadly accessible.