LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
Re: [lammps-users] Fixing the position of graphene sheet
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [lammps-users] Fixing the position of graphene sheet

From: Andrew Jewett <jewett@...1937...>
Date: Mon, 2 Oct 2017 15:45:01 -0700

On Thu, Sep 28, 2017 at 3:24 PM, Jatin Kashyap <jk435@...6426...> wrote:
> I’m trying to fix the position of the atoms, throughout whole simulation, of
> a graphene sheet while moving something that is sitting above it while they
> aren’t contacting each other. For that, I’m using following commands to
> constrain the motion of the atoms belonging to “fx_gr” group:
> velocity            fx_gr zero linear
> velocity            fx_gr zero angular

This fails because it does not set the velocity of every atom in the
"fx_gr" group to zero.  It sets the sum of the velocities to zero (and
the sum of the angular velocities).

If you're goal is to completely immobilize these atoms, then set each
atom's velocity to zero using something like this:

velocity fx_gr create 0.0 1
fix 3 fx_gr setforce 0 0 0

For more details:

Alternately you can omit the "fx_gr" group from the group of atoms
which is sent to "fix nvt" (or, equivalently, "fix nve"+"fix
langevin").  This is what I often recommend doing.  Either way should
work, but with "fix setforce", you must be careful to make sure the
velocities are initially zero and no other fixes are acting on the
atoms you want to freeze.

In your case, that would look something like:

group mobile_gr subtract all fx_gr
fix 3 mobile_gr nvt temp 300.0 300.0 100.0

This approach is used here:
For more details, see:

   --- NPT ---
At constant pressure, you can try using the alternate approach below
(although it might not be correct at high pressure):

If your goal is to allow the atoms some freedom of movement, then you
can either use:

1) fix spring   (to tether the atoms to their original locations)

2) fix recenter   (to keep the center of mass from changing)

The second option (fix recenter) only works if the atoms you want to
immobilize exert force on each other.  (If they are bonded to each
other or interacting via some many-body pair_style.)



> But after running simulation, carbons atoms belonging to “fx_gr” group
> appears to be moving right from 0th frame of the dump files. The simulation
> contains water and graphene interaction.
> Thanks to all for your time.
> ——
> Jatin Kashyap (
> Ph.D. Student
> DDLab Group
> Department of Mechanical and Industrial Engineering
> New Jersey Institute of Technology (NJIT)
> Newark, NJ 07102
> ------------------------------------------------------------------------------
> Check out the vibrant tech community on one of the world's most
> engaging tech sites,!
> _______________________________________________
> lammps-users mailing list