Revisiting the Rheology of Neutron Star Crusts with Molecular Dynamics

Explosive events from magnetars are likely due to the catastrophic release of stress in their crusts, but the behavior of crustal matter beyond linear elasticity is poorly understood. We argue here that seminal results from molecular dynamics informing crust breaking calculations are non-converged, and must be revisited. We estimate the criteria for quasi-static, rate-independent flow by comparing imposed deformation timescales to grain boundary diffusion in polycrystals. We argue that convergence in this regime should be observed at strain rates slower than $10^{-5}\,\omega_p$ (plasma frequency $\omega_p$) in simulations of $N\approx10^5$ particles across order 10 grains at a quarter of the melting temperature. Though computationally expensive, this is tractable with modern methods and GPU supercomputers.