Parallel Verlet neighbor list algorithm for GPU-optimized MD simulations