We investigate several important issues regarding the Random Batch Method (RBM) for second order interacting particle systems. We first show the uniform-in-time strong convergence for second order systems under suitable contraction conditions. Secondly, we propose the application of RBM for singular interaction kernels via kernel splitting strategy, and investigate numerically the application to molecular dynamics.