RE: question about Distance XY

From: Song, Hyundeok (songhk) (songhk_at_mail.uc.edu)
Date: Tue Sep 13 2011 - 09:10:07 CDT

Dear Jerome

I removed the harmonic parts, and it works. That's what I want.
I realized I want to keep an ion in cylinder of radius 10A, not inside a circle.
Thanks for your reply with clear explanation.

Hyun

________________________________________
From: heninj_at_gmail.com [heninj_at_gmail.com] on behalf of Jérôme Hénin [jhenin_at_ifr88.cnrs-mrs.fr]
Sent: Tuesday, September 13, 2011 3:50 AM
To: Song, Hyundeok (songhk)
Cc: namd-l_at_ks.uiuc.edu
Subject: Re: namd-l: question about Distance XY

Dear Hyun,

Your input activates two separate functions:
1) the half-harmonic restraint ("wall") included in the colvar
2) a generic harmonic restraint that keeps the colvar around the value
given by 'centers'

For your purpose, you only need the wall potential, so you can just
remove the 'harmonic' block. The ProjectionXY colvar should do the job
and restrain itself.

Finally, this restraint will not keep the ion in a horizontal plane,
it will keep it within a cylinder of radius 10 A. If you want to
restrain it to a disc-shaped region, you'll need to bias the distanceZ
coordinate for the same atoms.

Cheers,
Jerome

On 13 September 2011 06:27, Song, Hyundeok (songhk) <songhk_at_mail.uc.edu> wrote:
> Dear NAMD users
>
> First I red many postings about Distance XY from mailing list, but I didn't understand some parts.
>
> I want to restrain an ion on XY plane such as the circle of radius 10A
> ex: reference - atom groups from 1 to 32.
> main - one ion, number 50
> xydistane between a reference group and an ion.
>
> I included my colvar script. Please read my colvar script. If my script has a problem, please let me know.
> I have two questions.
> 1. Does the harmonic force will not apply as long as an ion is inside the boundary(10A from cm of references)?
>
> 2. Which value do I choose for 'centers', initial center of restraint ?
> Can I set 'centers' as '0' (A) ?
>
> Initially an ion is inside the circle, and if the harmonic force is applied only if an ion try to escape the circle of radius 10A, setting up 'centers' less than 10A is fine.
> I didn't specify the lowerboundary because it's zero.
>
> Thanks.
>
> Sincerely
> Hyun.
>
> ************************colvar script ****************************************
> colvar {
> name ProjectionXY
> upperboundary 10.0
> upperwallconstant 10.0
>
> distanceXY {
> ref {
> atomNumbersRange {1 - 32 }
> }
> main {
> atomnumbers { 50 }
> }
> }
> }
>
>
> harmonic {
> colvars ProjectionXY
> centers 0.0
> forceConstant 10
> }
>
>

This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:57:44 CST