       vars->gpu_mForcex,
       vars->gpu_mForcey,
       vars->gpu_mForcez,
-      vars->gpu_molIndex,
-      vars->gpu_kindIndex,
-      vars->gpu_lambdaVDW,
-      vars->gpu_lambdaCoulomb,
-      vars->gpu_isFraction,
       sc_coul,
       sc_sigma_6,
       sc_alpha,

  double distSq;
  double virX, virY, virZ;
  double pRF = 0.0, qi_qj, pVF = 0.0;
  //tensors for VDW and real part of electrostatic
  gpu_vT11[threadID] = 0.0, gpu_vT22[threadID] = 0.0, gpu_vT33[threadID] = 0.0;
  gpu_rT11[threadID] = 0.0, gpu_rT22[threadID] = 0.0, gpu_rT33[threadID] = 0.0;

                                gpu_rCutCoulomb[box], gpu_diElectric_1[0],
                                gpu_sigmaSq, sc_coul, sc_sigma_6, sc_alpha,
                                sc_power, lambdaCoulomb, gpu_count[0],
                                kA, kB);

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