the summation is written such that the sum runs from n=1, however, inside the sum, there is no n, but (r1aF1b + r2aF2b) instead, I am not quite sure how to interpret this…
Is there a paper/link on how the per atom stress calculation is performed in LAMMPS?
Sorry, what I wrote about kinetic energy term was wrong. Should not answer before thinking It has such form because it is the transfer of a-component of momentum along b direction. Look e.g. Heinz, Mol. Sim. 33 (2007).
the summation is written such that the sum runs from n=1, however, inside the sum, there is no n, but (r1aF1b + r2aF2b) > > instead, I am not quite sure how to interpret this…
The doc page says:
The second term is a pairwise energy contribution where n loops over the Np neighbors of atom I, r1 and r2 are the positions of the 2 atoms in the pairwise interaction, and F1 and F2 are the forces on the 2 atoms resulting from the pairwise interaction.
Is there a paper/link on how the per atom stress calculation is performed in LAMMPS?
yes, see the paper by Aidan Thompson (lead author) on the LAMMPS citation page, in J Chem Phys.
Thank you, so if I understand correctly, to keep the n inside the sum formulation what one has to say is r_n_a*F_n_b + r_(n+1)_a*F_(n+1)_b? The n has to be in the sum somehow, because otherwise it does not seem to make sense mathematically, as the sum goes from n=1, but n is not in the summation itself...
It’s not meant to be a precise mathematical formula. As the text states, it’s just a sum
over neighbors of atom I and R1 and R2 are the coords of I and J, the 2 atoms
in the pairwise interaction.