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

From: |
周永祥 <yxzhou1991@...34...> |

Date: |
Wed, 10 Dec 2014 21:09:44 +0800 (GMT+08:00) |

Hi all:

When I read your codes about computing rigid pressure, I found some problem . FIrst, You compute the constrained force , then compute virial term by equation : pxx = 0.5 * f.x * r.x, f.x means the constrained force in x direction, r.x means the unwrapped coordinate. But I think r.x should be the distance from partices to masscenter of the rigid. Then the computed pressure are almost zero because postive and negative are almost equal.

**Follow-Ups**:**Re: [lammps-users] compute rigid pressure***From:*Axel Kohlmeyer <akohlmey@...24...>

- Prev by Date:
**Re: [lammps-users] dihedral_style/opls parameters for lammps** - Next by Date:
**Re: [lammps-users] compute rigid pressure** - Previous by thread:
**[lammps-users] no output in rdf file with rerun - regd** - Next by thread:
**Re: [lammps-users] compute rigid pressure** - Index(es):