what's the value of energy-conversion constant in lammps

Dear Lammps users/developers,

In manual, I have found that the coul/cut style computes the standard Coulombic interaction potential given by E = Cqiqj/(epsilon*r). where C is an energy-conversion constant, and epsilon is the dielectric constant which can be set by the dielectric command (the constant is unitless). I want to know the value and unit of C.

In the mail list, I have found a mail talking about how lammps deals with C. however, the exact value isn’t shown. I need the value and unit of C for some analysis.

I would appreciate any help and suggestion.

Thanks!

Zhang Chao

Dear Lammps users/developers,

In manual, I have found that the coul/cut style computes the standard
Coulombic interaction potential given by E = C*qi*qj/(epsilon*r). where C is
an energy-conversion constant, and epsilon is the dielectric constant which
can be set by the dielectric command (the constant is unitless). I want to
know the value and unit of C.

In the mail list, I have found a mail talking about how lammps deals with C.
however, the exact value isn’t shown. I need the value and unit of C for
some analysis.

you can look this up in any text book describing Coulomb's law. C
stands for 1.0/(4*pi*epsilon_0), with epsilon_0 being the dielectric
permittivity.
in the unit and value of C in LAMMPS depends on your choice of units,
of course. the value LAMMPS uses for each set of units is assigned in
the function Update::set_units() in the file update.cpp

axel.