[lammps-users] Computation of total stress?

I am performing a uniaxial compression test and I am using buckingham, stillinger and web potential, and coul/long potentials. I am interested in computing the total stress in the system after applying a constant strainrate.

Since per-atom stress does not include contributions due to long-range Coulombic interactions and the virial term from the pressure computation (compute pressure command) is computed for all pairwise as well as 2-body, 3-body, and 4-body, and long-range interactions. My question is would it be more accurate to use compute pressure command to compute the total pressure in the system? I am not sure this is the way stress is computed in md when long range forces are considered

Compute stress/atom does not include long-range Coulombics. Compute
pressure does, in a global sense. Pressure and stress are the same thing,
modulo a volume factor.

Steve