From: Beaven, Andrew (a607b333_at_ku.edu)
Date: Fri Aug 02 2013 - 16:29:05 CDT
Hi - is there a way I can continually update the target position for a harmonic restraint? I'm trying to restrain two protein monomers at the maximal distance away from each other. I see that I could simply use the numerical value of the initial 0.25*$boxa (which is ~18 Å), but I'd like the target value to fluctuate with the box size during the simulation (since I'm starting from a non-equilibrated system). When I use the numerical value, the job runs and appears to keep the monomers in the desired locations, but as mentioned, this isn't exactly what I want…
In my NAMD input file (along with the other simulation variables), I use:
colvars on
colvarsConfig /path/restraint.col
The /path/restraint.col file is shown below.
Thank you,
Andrew
################################################
Colvarstrajfrequency 100
Colvarsrestartfrequency 100
colvar {
name proa_rest
distanceXY {
ref {
dummyAtom ( -0.25*$boxa, -0.25*$boxa, 0.000 )
}
main {
atomsFile /path/proa.ref
atomsCol B
atomsColValue 1.0
}
}
}
colvar {
name prob_rest
distanceXY {
ref {
dummyAtom ( 0.25*$boxa, 0.25*$boxa, 0.000 )
}
main {
atomsFile /path/prob.ref
atomsCol B
atomsColValue 1.0
}
}
}
harmonic {
colvars proa_rest
centers 0.0
forceConstant $fc
}
harmonic {
colvars prob_rest
centers 0.0
forceConstant $fc
}
################################################
This archive was generated by hypermail 2.1.6 : Wed Dec 31 2014 - 23:21:29 CST