LAMMPS WWW Site - LAMMPS Documentation - LAMMPS Mailing List Archives
Re: [lammps-users] LAMMPS moltemplate membrane+protein example
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [lammps-users] LAMMPS moltemplate membrane+protein example

From: Andrew Jewett <jewett@...1937...>
Date: Thu, 22 Jun 2017 14:00:30 -0700

It's true. It's all in the README file.  Thanks Axel.
You can download the code you need for that example at:

I created a new pair style for the silly reason that there are many different particle types in the example, and many of them are purely repulsive.  Right now LAMMPS does not really have an official pair style which supports automatic mixing of Lennard-Jones parameters for purely repulsive particles.  Rather than introduce a new pair style, I'm happy to modify the existing "lj/cut" pair_style to support this, if Steve will accept the change.  (About 8 lines of new code.)

Here are several more comments about this example:

1) The lipid mixture does not work very well. This is my fault.  The original paper by Brannigan&Brown only describes one type of lipid (DPPC), and parameters were chosen accordingly.  I deleberately added a second lipid type to demonstrate how to mix different molecules together using the moltemplate "new random" command.  However this causes strange defects in the membrane over long timescales.  If you don't need a mixture of lipids, just use DPPC and ignore DLPC.

2) the protein model is still a little too long to fit comfortably in the membrane, and it may tilt slightly during the simulation.  Getting rid of the  "DLPC" lipids will probably help, since they are shorter and they make the membrane thinner.

3) David Stelter and Saeed Momeni Bashusqeh have added the Cooke and Desrno lipid model into moltemplate (as well as well as many lipids from the MARTINI and SDK force fields into moltemplate), so there are many other membrane models you can choose from.  (Assuming you don't want to invent your own...  I need to remember to add the new lipids examples to the GitHub page.  If I forget, please remind me.)

Glad to see people are using moltemplate!(especially for coarse grained biomolecules)



On Jun 22, 2017 3:05 AM, "Axel Kohlmeyer" <akohlmey@...24...> wrote:

On Thu, Jun 22, 2017 at 5:28 AM, Bhargav Patel <bhargavpatel091097@...24...> wrote:

So i've tried to run the membrane +protein example of moltemplate and it gives the following error:

ERROR: Illegal pair_style command (../pair_table.cpp:238)
Last command:     pair_style  hybrid table linear 1130                 lj/charmm/coul/charmm/inter es4k4l 14.5 15

I have downloaded the additional LAMMPS requirements from the moltemplate site, added it to the src folder and recompiled LAMMPS. I have also added the USER-MISC package, using the following commands:

make clean-all
make yes-USER-MISC
make serial

Yet the error persists. Please help.

​pair style ​  lj/charmm/coul/charmm/inter is not part of the LAMMPS distribution.
please read the README files. they are called "README" for a reason!

I am using the 17th November 2016 version of LAMMPS in Ubuntu

Thanking You,
Bhargav Patel

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.