Optimizing and improving the fidelity of reactive, polarizable molecular dynamics simulations on modern high performance computing architectures