Improving the convergence order of binary neutron star merger simulations in the Baumgarte- Shapiro-Shibata-Nakamura formulation