   double distSq;
   double virX, virY, virZ;
   double pRF = 0.0, qi_qj, pVF = 0.0;
+  double lambdaVDW = 0.0, lambdaCoulomb = 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;

      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,

      vars->gpu_rMin,
      vars->gpu_rMaxSq,
      vars->gpu_expConst,
      vars->gpu_molIndex,
      vars->gpu_kindIndex,
      vars->gpu_lambdaVDW,
      vars->gpu_lambdaCoulomb,
      vars->gpu_isFraction,
      box);

  checkLastErrorCUDA(__FILE__, __LINE__);
