From: Jérôme Hénin (jhenin_at_ifr88.cnrs-mrs.fr)
Date: Wed Jun 16 2010 - 11:42:09 CDT
Mike:
Just a small addition to Giacomo's explanation. The "half-harmonic"
boundary potential is the following:
V(x) = 0.5 * k * (x - upperWall)^2 if x > upperWall
V(x) = 0 otherwise
and the default value of upperWall is upperBoundary.
Jerome
On 16 June 2010 18:20, Giacomo Fiorin <giacomo.fiorin_at_temple.edu> wrote:
> Hi Michael, upperBoundary is chosen by you. There is no way to tell which
> is the highest possible RMSD value without knowing the system. You chose 2
> Å for a specific reason (i.e. to keep the protein close to the ref
> structure), but in general it can be much higher.
>
> Anyway, the force will kick in only after RMSD goes above 2 Å: no matter how
> strong the force constant, you will always have some RMSD values a little
> over 2 Å. There is no way to enforce it as a hard limit without using
> infinitely large forces (which is actually what would happen if the code
> were to automatically enforce the hard boundary). However, you can control
> the tolerance by lowering or increasing upperWallConstant.
>
> Also, although you can set upperBoundary wherever you want, its more
> appropriate usage is as you suggested in your mail, as the highest
> theoretically possible above which the RMSD really never goes. Use
> lowerWall, instead, to put, a wall on a smaller RMSD value, and keep in mind
> that you'll always have to allow for some tolerance.
>
> Giacomo
>
> ---- ----
> Dr. Giacomo Fiorin
> ICMS - Institute for Computational Molecular Science - Temple University
> 1900 N 12 th Street, Philadelphia, PA 19122
> giacomo.fiorin_at_temple.edu
> ---- ----
>
>
>
> On Wed, Jun 16, 2010 at 11:24 AM, Michael Zimmermann <michaelz_at_iastate.edu>
> wrote:
>>
>> It looks like dos2unix fixed that initial error.
>>
>> I've never had carriage control interfere with namd, so thanks for the
>> pointer.
>>
>> It is still not clear to me how "upperBoundary" and "upperWallConstant"
>> are interpreted in a colvar that calculates RMSD. Is the upperBoundary the
>> greatest allowable RMSD? If so, how will it accomplish such a bound? Will
>> it increase the spring stiffness sufficient to not allow further RMSD
>> changes? Is the upperWallConstant then the greatest allowable spring
>> constant?
>>
>> The user's guide gives more general descriptions:
>> *lowerWallConstant < (colvar) Lower wall force constant (kcal/mol) >
>> Description: If lowerWall or lowerBoundary is defined, provides the force
>> constant. The
>> energy unit of the constant is kcal/mol, while the spatial unit is that of
>> the colvar.
>> *upperWallConstant < (colvar) Upper wall force constant (kcal/mol) >
>> Description: Similar to lowerWallConstant.
>>
>> Thanks again.
>>
>> Mike Z
>>
>> On Wed, Jun 16, 2010 at 7:24 AM, Jérôme Hénin <jhenin_at_ifr88.cnrs-mrs.fr>
>> wrote:
>>>
>>> Hi Mike,
>>>
>>> I cannot reproduce this error, something in the file seems to confuse
>>> the parser. If the file has been edited under windows at any point,
>>> you could try treating it with dos2unix.
>>>
>>> Another point, though: the last two keywords, "upperBoundary" and
>>> "upperWallConstant" are colvar-related, not colvar-component-related,
>>> so they should be outside the rmsd block (but still within the colvar
>>> block).
>>>
>>> Jerome
>>>
>>> On 15 June 2010 22:59, Michael Zimmermann <michaelz_at_iastate.edu> wrote:
>>> > Thank you all for your replies thus far. I am making progress, but
>>> > still
>>> > running into problems with the RMSD colvar.
>>> >
>>> > I am running NAMD version 2.7b2 Linux-x86_64.
>>> >
>>> > ------------------------------------------------------------
>>> > ---> The contents of my colvar.conf file is:
>>> > ------------------------------------------------------------
>>> > colvar {
>>> > name RMSDColvar1
>>> > rmsd {
>>> > atoms {
>>> > # this file has occupancy marked as 2 for the atoms
>>> > # we want restrained by RMSD
>>> > atomsFile prot.colvar.mark.pdb
>>> > atomsCol O
>>> > atomsColValue 2.0
>>> > }
>>> > refPositionsFile prot.colvar.mark.pdb
>>> > refPositionsCol O
>>> > refPositionsColValue 2.0
>>> > upperBoundary 2 # does this mean the RMSD won't go above 2A?
>>> > upperWallConstant 10 # is this the force constant used on
>>> > springs to
>>> > restrain the atom group?
>>> > }
>>> > }
>>> >
>>> > ------------------------------------------------------------
>>> > ---> I don't know what the error I receive means. All my parameters
>>> > have a
>>> > value as far as I can tell
>>> > ------------------------------------------------------------
>>> > colvars: Initializing a new collective variable.
>>> > colvars: # name = "colvar1" [default]
>>> > colvars: Initializing a new "RMSD" component.
>>> > colvars: # componentCoeff = 1 [default]
>>> > colvars: # componentExp = 1 [default]
>>> > colvars: Initializing atom group "atoms".
>>> > colvars: Error: keyword "" is not supported, or not recognized in
>>> > this
>>> > context.
>>> > colvars: If this error message is unclear, try recompile with
>>> > -DCOLVARS_DEBUG.
>>> > FATAL ERROR: Error in the collective variables module: exiting.
>>> > ------------------------------------------------------------
>>> >
>>> > Any advice on how to continue would be appreciated.
>>> >
>>> > Mike Z
>>> >
>>> >
>>> > On Thu, Jun 10, 2010 at 5:20 PM, Giacomo Fiorin
>>> > <giacomo.fiorin_at_temple.edu>
>>> > wrote:
>>> >>
>>> >> Michael: again, which version of NAMD are you using? There isn't
>>> >> unfortunately a wide set of examples currently available, but the NAMD
>>> >> 2.7b2
>>> >> Users Guide contain exactly the example you're looking for, with two
>>> >> atomic
>>> >> distances and an harmonic restraint applied to both. You just need
>>> >> to fill
>>> >> in the atom numbers and the distance values that you want to restrain.
>>> >>
>>> >> If you have too many interatomic distances to restrain, you may also
>>> >> try
>>> >> the "rmsd" variable and apply a harmonic potential to it, with center
>>> >> on
>>> >> zero.
>>> >>
>>> >> Giacomo
>>> >>
>>> >> ---- ----
>>> >> Dr. Giacomo Fiorin
>>> >> ICMS - Institute for Computational Molecular Science - Temple
>>> >> University
>>> >> 1900 N 12 th Street, Philadelphia, PA 19122
>>> >> giacomo.fiorin_at_temple.edu
>>> >> ---- ----
>>> >>
>>> >>
>>> >>
>>> >> On Thu, Jun 10, 2010 at 4:15 PM, Michael Zimmermann
>>> >> <michaelz_at_iastate.edu>
>>> >> wrote:
>>> >>>
>>> >>> I am having trouble finding an example colvarsConfig file. NAMD Some
>>> >>> of
>>> >>> the Colvar documentation says that the format is similar to a regular
>>> >>> NAMD
>>> >>> configuration file, but what parameters need to be set in it?
>>> >>>
>>> >>> All I want is to restrict the internal motion of a subset of atoms in
>>> >>> my
>>> >>> simulation. I do not want them fixed in space, but I want their
>>> >>> relative
>>> >>> positions to remain roughly unchanged.
>>> >>>
>>> >>> I do not have the colvarsConfig parameter set (no calvar config file)
>>> >>> and
>>> >>> the simulation has not thrown any errors.
>>> >>>
>>> >>> Mike
>>> >>>
>>> >>> On Tue, Jun 8, 2010 at 1:11 PM, Michael Zimmermann
>>> >>> <michaelz_at_iastate.edu>
>>> >>> wrote:
>>> >>>>
>>> >>>> Dear NAMD users,
>>> >>>>
>>> >>>> I am having trouble finding the proper syntax for atom assignments
>>> >>>> within colvars. I want to keep a group of atoms that are in
>>> >>>> different
>>> >>>> segments constrained by RMSD to a reference position. This is what
>>> >>>> I have
>>> >>>> in my config file so far, but there seems to be something wrong...
>>> >>>>
>>> >>>> ================================
>>> >>>> colvars on
>>> >>>> colvar {
>>> >>>> rmsd {
>>> >>>> atoms {
>>> >>>> psfSegID P6
>>> >>>> atomNameResidueRange 231-235
>>> >>>> psfSegID P16
>>> >>>> atomNameResidueRange 232-235
>>> >>>> psfSegID P29
>>> >>>> atomNameResidueRange 232-235
>>> >>>> psfSegID P33
>>> >>>> atomNameResidueRange 25-28
>>> >>>> psfSegID P33
>>> >>>> }
>>> >>>> refPositionsFile markedProt.colvar.pdb
>>> >>>> refPositionsCol B
>>> >>>> }
>>> >>>> }
>>> >>>> ================================
>>> >>>>
>>> >>>> --
>>> >>>> Michael Zimmermann
>>> >>>> Ph.D. student in Bioinformatics and Computational Biology
>>> >>>> Department of Biochemistry, Biophysics and Molecular Biology
>>> >>>> Iowa State University
>>> >>>
>>> >>>
>>> >>>
>>> >>> --
>>> >>> Michael Zimmermann
>>> >>> Ph.D. student in Bioinformatics and Computational Biology
>>> >>> Department of Biochemistry, Biophysics and Molecular Biology
>>> >>> Iowa State University
>>> >>
>>> >
>>> >
>>> >
>>> > --
>>> > Michael Zimmermann
>>> > Ph.D. student in Bioinformatics and Computational Biology
>>> > Department of Biochemistry, Biophysics and Molecular Biology
>>> > Iowa State University
>>> >
>>
>>
>>
>> --
>> Michael Zimmermann
>> Ph.D. student in Bioinformatics and Computational Biology
>> Department of Biochemistry, Biophysics and Molecular Biology
>> Iowa State University
>
>
This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:55:53 CST