<br> Thanks for Mark's fast reply, and It really works.<br> Then, I have another qusetion about the manual: It seems that the exclusion just excludes the nonbond potential of LJ interaction, but Buckingham. If we want to turn off all VDW interaction between some atoms, the energy monitor group is a nice choice. I have no idea if I got it right.<br>
And, we can't use LJ and Buckingham in a same system. Unfortunately, I met a force-field that most of the nonbond interactions were described by Buckingham, but LJ only for Ow-Ow among water. I chose nbfunct=2 in [defaults] section,and nbfunct=1 in [nonbond_params] for Ow-Ow. As a result, a error message read: <b>Trying to add LJ (SR) while the default nonbond type is Buck.ham (SR). </b>I don't know whether it could be settled in gromacs.<br>
<br> Sincerely<br> Shijun<br> <br><div class="gmail_quote"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
----------------------------------------------------------------------<br></blockquote><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
Message: 1<br>
Date: Mon, 21 Sep 2009 06:24:34 +0000 (GMT)<br>
From: Mark Abraham <<a href="mailto:mark.abraham@anu.edu.au">mark.abraham@anu.edu.au</a>><br>
Subject: Re: [gmx-users] Re: who can help me with the force-field?<br>
To: Discussion list for GROMACS users <<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a>><br>
Message-ID: <<a href="mailto:fbcb8a00405b.4ab71c22@anu.edu.au">fbcb8a00405b.4ab71c22@anu.edu.au</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
On 09/21/09, xiao shijun <<a href="mailto:xshijun@gmail.com">xshijun@gmail.com</a>> wrote:<br>
><br>
> Mark wrote:<br>
> > xiao shijun wrote:<br>
><br>
> >> Hi everyone,<br>
> >> Recently, I met a problem on the force-field. I am doing something<br>
> >> on the calcium carbonate, and want to simulation it with gromacs.<br>
> >> The shell model for oxygen has been emploied in my system. When I<br>
> >> considered the non-boned potential associated with interactions between<br>
> >> atoms<br>
> >> in the same CO3 group, I only need to take O-O buckingham interaction<br>
> >> into account, but any electric interaction. How could I realize that in<br>
> >> my .top file.<br>
> >> First, I tried to use [pairs] to eliminate electric interaction, but<br>
> >> I can't add the buckingham interaction of O-O because there is no<br>
> >> buckingham<br>
> >> options for LJC14. I want to know if anyone of you have such the<br>
> >> experience. Thanks!!<br>
> <br>
> >Usually you would not define a nonbonded interaction for atoms that are<br>
> >so close together, since their interaction should be able to be modelled<br>
> >adequately with the C-O bond and O-C-O angle (and perhaps with CHARMM,<br>
> >an interaction harmonic the O-O distance as a Urey-Bradley term).<br>
> <br>
> Thanks for Mark's suggestion. But I want to reproduce the simulation from the other paper,<br>
> The paper describe the force-field that O-O nonbonded interaction in the same carbonate should be considered.<br>
> In addition to that, I can't use buckingham potential in [pairs], even choose funct = 2. I am wondering<br>
> if I want to describe 1-4 interaction by Buckingham, what could I do?<br>
This is one application for setting nrexcl for moleculetypes. You don't want exclusions for particles 2 bonds apart, so set this to one. See chapter 5.<br>
<br>
If that doesn't work, then if the "bonded" interactions are merely harmonic distance potentials, then GROMACS permits a form of them that doesn't create a bond (bond type 6, IIRC, check chapter 5; or constraint type 2). Thus, the automated exclusion mechanism won't be involved, and you can have Buckingham O-O so long as you have Buckingham for all other LJ interactions. You will need manual exclusions for C-O pairs, I guess. Thus<br>
<br>
[ atomtypes ]<br>
C #nonbonded parameters go here<br>
<br>
O #nonbonded parameters go here<br>
<br>
[moleculetype]<br>
carbonate<br>
<br>
[atoms]<br>
C<br>
O1<br>
O2<br>
O3<br>
<br>
[bonds]<br>
C-O1 # type 6<br>
C-O2 # type 6<br>
C-O3 # type 6<br>
<br>
[exclusions]<br>
1 2<br>
1 3<br>
1 4<br>
<br>
Whether this works may depend on whatever else is in your system, though.<br>
<br>
Mark<br>
<br>
<br>
</blockquote></div><br>