PULSEDYN: Open source code for dynamical simulations of strongly nonlinear systems

We present PULSEDYN which is an open source C++ code written with the goal of performing highly accurate and late time simulations on strongly and weakly nonlinear chains while requiring minimal effort on the user’s part. PULSEDYN offers a suite of potentials and solvers and simulations and may be set up using a parameter file and running the code provided as an executable. PULSEDYN is also written with an emphasis on modularity, allowing for easy changes to be made to the code. PULSEDYN is available to download from Github and is accompanied by documentation and a detailed user manual [1]. We show standard results reproduced from literature using PULSEDYN for benchmarking purposes. We then use PULSEDYN to demonstrate how a strongly nonlinear Fermi-Pasta-Ulam-Tsingou system transitions from quasi-equilibrium to the equipartitioned state and verify our results by correctly predicting the exact equilibrium specific heat of the system from the simulated equipartitioned state.
[1] R. Kashyap and S. Sen, Comput. Phys. Commun. (submitted).