Primary and secondary modes of deformation twinning in HCP Mg based on atomistic simulations