       gpu_rT23[threadID] = pRF * (0.5 * (virY * diff_comz + virZ * diff_comy));
     }
 
-    pVF = CalcEnForceGPU(distSq, gpu_particleKind[gpu_pair1[threadID]],
-                         gpu_particleKind[gpu_pair2[threadID]],
+    pVF = CalcEnForceGPU(distSq, kA, kB,
                          gpu_sigmaSq, gpu_n, gpu_epsilon_Cn, gpu_rCut[0],
                          gpu_rOn[0], gpu_isMartini[0], gpu_VDW_Kind[0],
                          gpu_count[0], lambdaVDW, sc_sigma_6, sc_alpha,

                                gpu_rCutCoulomb[box], gpu_diElectric_1[0],
                                gpu_sigmaSq, sc_coul, sc_sigma_6, sc_alpha,
                                sc_power, lambdaCoulomb, gpu_count[0],
                                gpu_particleKind[gpu_pair1[threadID]],
                                gpu_particleKind[gpu_pair2[threadID]]);

      gpu_rT11[threadID] = pRF * (virX * diff_comx);
      gpu_rT22[threadID] = pRF * (virY * diff_comy);

      vars->gpu_mForcex,
      vars->gpu_mForcey,
      vars->gpu_mForcez,
      sc_coul,
      sc_sigma_6,
      sc_alpha,
