Hi Arun,

Thanks for the information. But I’m still not clear about this whole stress thingy. Probably I’m not framing my question properly. Here is my situation in a more detailed way.

I am more interested in the triangular components of the stress tensor matrix Sxy, Sxz, and Syz. So, What I am actually looking for are these tensor components for the entire system at each time step. Here is the sample data which I want the stress data to look like

(desired Output)

time step Sxy Sxz Syz (for the entire system)

1 **** **** ****

2 **** **** ****

3 **** **** ****

and so on…

But here is what LAMMPS data looks like

time step

1

atom id atom type Sxy Sxz Syz (each individual atom)

1 1 **** **** ****

2 1 **** **** ****

…

and so on…N values

time step

2

atom id atom type Sxy Sxz Syz

1 1 **** **** ****

2 1 **** **** ****

…

and so on…N values

(N is the total number of atoms in the system.

LAMMPS outputs the individual tensor components at every time step. My question is how do I convert them to look like the above desired result. I do not want the tensor components of all the individual atoms. Instead, I want the total value of each component for the system.

This is what I did

SXY = Sigma (Sxy) … added up the contributions of all the atoms

SXZ = Sigma (Syz)

SYZ = Sigma (Syz)

This is not giving me the same result as another software that was used to run a simulation using the same inital configeration. The trajectories match but the stress tensor components do not.

I would like to know wether what I did (summation of all the contributions) is the right way to get the stress tensor components for the system or amI missing on something. If not, could you please take some time to help me figure out this.

Thank you very much

Venkat