LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
Re: [lammps-users] Using GROMOS Force Field in lammps
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [lammps-users] Using GROMOS Force Field in lammps

From: Axel Kohlmeyer <akohlmey@...24...>
Date: Mon, 16 Oct 2017 15:12:42 -0400

On Mon, Oct 16, 2017 at 12:25 AM, Faust SHI <shi86@...107...> wrote:
> Hi,
> I want to use GROMOS Force Field in lammps. I noticed that I can use
> some_styles in lammps source field as the energy functions for angle and
> dihedral. And I can also use lj/gromacs for non-bonded interactions. However
> I cannot find that for bond calculation while for GROMOS it should be k(b^2
> - b0^2)^2. I think it should be similar with the harmonic bond style and I
> want to adapt the cpp file for harmonic bond style in src folder so that I
> can use that to calculate bond interaction in GROMOS force field. What
> variables in the cpp file should I change? Just the variable for energy like
> energy = k(b - b0)^2 to energy = k(b^2 - b0^2)^2? Or do I need to change the
> variables for others like force?

you need to include forces as well. since this is straightforward to
implement, i've just created this bond style and submitted it for
inclusion into the next LAMMPS patch release.
you can find the pull request here:


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

Dr. Axel Kohlmeyer  akohlmey@...24...
College of Science & Technology, Temple University, Philadelphia PA, USA
International Centre for Theoretical Physics, Trieste. Italy.