On 7/23/12 8:31 AM, tarak karmakar wrote:
Thanks for the reply.
But then what should I give as Kb [ bond stretching constant ] in the
ffbonded.itp file ? should I give 00.000 there or something else ???


Constraints make bonds rigid.  The force constant is irrelevant.

-Justin

On Mon, Jul 23, 2012 at 4:11 PM, Justin Lemkul <jalem...@vt.edu> wrote:


On 7/23/12 2:39 AM, tarak karmakar wrote:

Dear All,


    In my simulation I want to fix the S=O bond length of SO2, but the
angle has to be kept flexible. So I search for the distance restraints
in gromacs mailing list and
accordingly I have incorporated the distant_restraints block in the
topology file. After a short equilibration, when I am seeing the movie
of the trajectory, the bond length is changing. So what is the right
way to keep the distance fixed and the angle flexible ?


Use constraints for the bonds, not distance restraints.

-Justin


     Part of the topology file is as follows


; Include forcefield parameters
#include "./amber94ff/forcefield.itp"

[ moleculetype ]
; Name            nrexcl
Other_chain_A       3

[ atoms ]
;   nr       type  resnr residue  atom   cgnr     charge       mass
typeB    chargeB      massB
; residue 502 SR  rtp CR   q  0.0
       1         S    502     SR      S      1     0.6512      32.07  ;
qtot 0.6512
       2         O    502     SR    O1      2    -0.3256         16   ;
qtot 0.3256
       3         O    502     SR    O2      3    -0.3256         16   ;
qtot 0

[ bonds ]
;  ai    aj funct            c0            c1            c2            c3
      1     2     1
      1     3     1

[ distance_restraints ]
;  ai    aj type index type' low up1 up2 fac
     1     2     1     0      1   1.149  1.490  1.490  1.0
     1     3     1     0      1   1.149  1.490  1.490  1.0

[ angles ]
;  ai    aj    ak funct            c0            c1            c2
c3
      2     1     3     1

; Include Position restraint file
#ifdef POSRES
#include "posre.itp"
#endif

; Include water topology
#include "./amber94.ff/spc.itp"

#ifdef POSRES_WATER
; Position restraint for each water oxygen
[ position_restraints ]
;  i funct       fcx        fcy        fcz
     1    1       1000       1000       1000
#endif

; Include topology for ions
#include "./amber94.ff/ions.itp"

[ system ]
; Name
Protein in water

[ molecules ]
; Compound        #mols
Other_chain_A       1
SOL               622


--
========================================

Justin A. Lemkul, Ph.D.
Research Scientist
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin

========================================


--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Only plain text messages are allowed!
* Please search the archive at
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the www
interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists




--
========================================

Justin A. Lemkul, Ph.D.
Research Scientist
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin

========================================


--
gmx-users mailing list    gmx-users@gromacs.org
http://lists.gromacs.org/mailman/listinfo/gmx-users
* Only plain text messages are allowed!
* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/Search before posting!
* Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org.
* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

Reply via email to