|From:||Axel Kohlmeyer <akohlmey@...24...>|
|Date:||Mon, 28 Aug 2017 14:38:00 -0400|
On Mon, Aug 28, 2017 at 12:14 AM, Trung Nguyen <ndactrung@...24...> wrote:Hi Wes,it's possible to support triclinic boxes for fix rigid/npt and rigid/nph either by following the reference cited the fix rigid doc (Kamberaj et al), or by borrowing what is implemented in fix nh for updating the full box dimension tensor.However, please note that fix rigid/npt and rigid/nph with triclinic boxes would be mostly useful for the cases where you only have rigid bodies in the system. For systems of polymer-grafted rigid bodies like yours, I think the workaround along the line Axel suggested would work to the similar effect:# assuming you're having a triclinic box setupgroup nanoparticles type 1group tethers type 2fix 1 nanoparticles rigid/nvt temp ...fix 2 tethers npt temp 1.0 1.0 1.0 xy 1.0 1.0 10.0 ... dilate allfix_modify 2 press thermo_pressThe reason that it would work is that the simulation triclinic box is to be barostatted by the fix npt, which uses the system pressure tensor as feedback.Thanks Axel and Trung, this helps. I now see that this method is one of the options listed on the "fix rigid" page for using pressure coupling with a mixture of rigid and non-rigid molecules.Is the thermo_press compute modification (the last line) always necessary in a system like this?
--That's one step I did not see on the "fix rigid" page's section on how to use pressure coupling with a system containing both rigid and non-rigid bodies.