LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
Re: [lammps-users] Fix rigid and minimization
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [lammps-users] Fix rigid and minimization

From: Stefan Paquay <stefanpaquay@...24...>
Date: Wed, 2 May 2018 13:43:19 -0400

There's pretty much two ways to go about this:
1. Run normal dynamics but slowly ramp down your temperature and apply some viscous force (like fix langevin or fix viscous) to prevent your simulation from blowing up. The ramping down of the temperature is "annealing".
2. Replace the rigid body by a body that is held together by very stiff springs, and after minimizing, activate fix rigid to perform your normal dynamics.

The complexity your input and rigid bodies determine which is the most convenient/straightforward option.

On Wed, May 2, 2018 at 1:33 PM, Mohammad Rasool Vazirisereshk <mvazirisereshk@...24...> wrote:
Hi everyone,

I am trying to perform a molecular static simulation for a group of atoms residing on a substrate (both atom types are Copper). I am forced to keep the top group of the atoms as rigid (they need to have rigid body movement). I would like to use fix rigid however this fix is not invoked during the minimization. I cannot use fix setforce 0 0 0 instead of fix rigid because these atoms are coupled with spring to another atom and therefore the force applied to the rigid group cannot be ignored. 

Is there any workaround to have fix rigid and minimize working together? I know this might be a long shot but has anyone written/modified the LAMMPS minimize to include fix rigid? Looking at the mailing archives, I saw "simulated annealing MD run" might work instead of minimization. Can someone explain to me what annealing MD is?

Thank you in advance.


Check out the vibrant tech community on one of the world's most
engaging tech sites,!
lammps-users mailing list