[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

From: |
Abolfazl Musavi <ab.musavi1990@...24...> |

Date: |
Tue, 5 Dec 2017 13:24:29 +0330 |

Hi

Sorry for that. I mean the manual section for compute pressure command (http://lammps.sandia.gov/doc/compute_pressure.html ) where it is written:

P=NkBT/V+sum_i_N'(r_i.f_i)

"where r_i and f_i are the position and **force vector of atom i**". It seems that fi is the "total force" on atom i.

while in the reference paper (Thompson, Plimpton, Mattson, J Chem Phys, 131, 154107 (2009)) it is told:

P=NkBT/V+<W>/3V (1)

P=NkBT/V+<W>/3V (1)

and

W(rN)=sum_n(sum_i_1_N(r_in.F'in)) (25)

"**F'in is not the total force** on an atom but rather the partial force ".

Is
the lammps method of calculating pressure according to eq.25? If yes,
why the forces are total in the manual and partial in the paper?

thanks for any help to solve my confusion.

thanks for any help to solve my confusion.

Regards

On 5 December 2017 at 05:48, Axel Kohlmeyer <akohlmey@...24...> wrote:

On Mon, Dec 4, 2017 at 9:08 PM, Abolfazl Musavi <ab.musavi1990@...24...> wrote:Hello everyone

According to the reference of compute pressure command (Thompson et al. 2009), the virial term could be calculated with eq.25 (the atom form:W=sum_n sum_i(r_in.F'_in)[i=1:N]in which the "F'in is not the total force on an atom but rather the partial force". But, I could not see in the corresponding section of manual mentioning that the forces are "partial". (W=sum_i(r_i.f_i) [i=1:N')]the LAMMPS manual has over 2000 pages. would you mind being a little more precise than "the corresponding section"?axel.could anyone explain this discrepancy for me please?thank you in advance for your timebest regards

------------------------------------------------------------ ------------------

Check out the vibrant tech community on one of the world's most

engaging tech sites, Slashdot.org! http://sdm.link/slashdot

_______________________________________________

lammps-users mailing list

lammps-users@...655....net

https://lists.sourceforge.net/lists/listinfo/lammps-users

--Dr. Axel Kohlmeyer akohlmey@...24... http://goo.gl/1wk0

College of Science & Technology, Temple University, Philadelphia PA, USA

International Centre for Theoretical Physics, Trieste. Italy.

**Follow-Ups**:**Re: [lammps-users] partial force in virial***From:*Axel Kohlmeyer <akohlmey@...24...>

**References**:**[lammps-users] partial force in virial***From:*Abolfazl Musavi <ab.musavi1990@...24...>

**Re: [lammps-users] partial force in virial***From:*Axel Kohlmeyer <akohlmey@...24...>

- Prev by Date:
**[lammps-users] fix deform alternate** - Next by Date:
**[lammps-users] 12-bit tables for long-range coulomb** - Previous by thread:
**Re: [lammps-users] partial force in virial** - Next by thread:
**Re: [lammps-users] partial force in virial** - Index(es):