Re: [ESPResSo-users] multiple potentials on a polymer

[Olaf Lenz]

Hi!

Again, please write to the mailing list, not an individual person.

On 01/10/2011 12:06 PM, Edward Kim wrote:
> Thank you very much!!
> 
> FENE + harmonic was just an example. I am just trying to produce a
> rigid structure that was implemented with FENE +
> "V(theta)=K(1+cos(theta))". (it looks like bond angle cosine)

This means that besides the bond-length potential, you also want to have
a bond-angle potential. This is a very common case. Note, that a
bond-angle potential is a three-body potential as opposed to a
bond-length two body potential!
To set this up, you will have to use the part command for the various
particles, as you can read in the User's Guide, sec. 5.3.5.

> By the way, In the case of using syntax "polymer", can I just use
> another polymer command with harmonic bond? (please see (1) below) or
> use part command (2)?

You have to use the "part" command. "polymer" will always generate new
particles.

Olaf
-- 
Dr. rer. nat. Olaf Lenz
Institut für Computerphysik, Pfaffenwaldring 27, D-70569 Stuttgart
Phone: +49-711-685-63607

olenz.vcf
Description: Vcard