Thank you for your response, Dr. Moore. I am using real units for my simulation and not sure if the prefactor will be 1.0 in this case. In the case of coulombic pair style, a prefactor C is defined which is the same as 1/(4piepsilon_0). See:
The value of this prefactor C for real units is defined in update.cpp as:
force->qqr2e = 332.06371; // see also force->qqr2d_lammps_real