Improving efficiency of large time‐scale molecular dynamics simulations of hydrogen‐rich systems