Re: Definition of two colvars for ABF

From: Jérôme Hénin (jhenin_at_ifr88.cnrs-mrs.fr)
Date: Mon Jan 04 2010 - 16:59:22 CST

Dear Branko,

So far, to use abf_integrate under windows, you need to compile it
yourself (the source files are in the lib/ directory of the NAMD
source distribution). I have no experience in windows programming, so
I can't really help you with that.

The lack of a colvars.state file sounds like a problem. What happens
if you set a fairly low colvarsRestartFrequency?

Best,
Jerome

2010/1/1 Branko <bdrakuli_at_chem.bg.ac.yu>:
> Dear All,
>
> Introducing of boundaries within "gyration colvar" solve the problem and
> simulation proceed to the end. Good choice is that one have some previous
> simulation of the similar type to estimate boundaries. Now I appreciate if
> one  could  advice me how to  made integration. As I use Win version and
> abf_integrate exist only under the linux version, I try to copy whole folder
> within "lib" directory, but "integrate" or "abf_integrate" can not be
> recognized as a command (command prompt with properly chosen path). Also can
> VMD console could be helpful?. Thanks in advance.
>
> Branko
>
> Branko wrote:
>
> Dear All,
>
> I am trying to define two colvars for ABF (2.7b1 Windows version), using
> configuration given below, all atoms belong to the same molecule and atoms
> in first colvar do not overlapped with the atoms of the second one (i.e. are
> orthogonal). Simulation using the  first colvars only proceed smoothly, but
> the  second  one  is likely the problem.  So the question is are boundaries
> for  "gyration" must be defined, if  are requisite how to estimate values .
> As a additional note still I can't find  colvars.state file during
> successfully simulation. Best wishes to developers and users in New 2010.
>
> Kind regards
>
> Branko
>
> Output (problem underlined):
>
> colvars:
> ----------------------------------------------------------------------
> colvars: Initializing the collective variables module, version 29-09-2008.
> colvars: # colvarsTrajFrequency = 1
> colvars: # colvarsRestartFrequency = 10000
> colvars: # trajAppend = off [default]
> colvars: The restart output state file will be
> "C:\namd2\246triEtinthd.restart.colvars.state".
> colvars: The trajectory file will be "C:\namd2\246triEtinthd.colvars.traj".
> colvars: The final output state file will be
> "C:\namd2\246triEtinthd.colvars.state".
> colvars: # analysis = on
> colvars:
> ----------------------------------------------------------------------
> colvars: # readTrajectory = "" [default]
> colvars: # readBegin = 0 [default]
> colvars: # readEnd = 0 [default]
> colvars:
> ----------------------------------------------------------------------
> colvars:   Initializing a new collective variable.
> colvars:   # name = d
> colvars:   Initializing a new "distance" component.
> colvars:     # componentCoeff = 1 [default]
> colvars:     # componentExp = 1 [default]
> colvars:     # oneSiteSystemForce = off [default]
> colvars:       Initializing atom group "group1".
> colvars:       Atom group "group1" defined, 3 atoms, total mass = 36.0321.
> colvars:       Initializing atom group "group2".
> colvars:       Atom group "group2" defined, 3 atoms, total mass = 44.0095.
> colvars:   All components initialized.
> colvars:   # width = 0.01
> colvars:   # lowerBoundary = 3.35
> colvars:   Lower boundary defined.
> colvars:   # upperBoundary = 9.7
> colvars:   Upper boundary defined.
> colvars:   # lowerWallConstant = 0 [default]
> colvars:   # upperWallConstant = 0 [default]
> colvars:   # extendedLagrangian = off [default]
> colvars:   # outputValue = on [default]
> colvars:   # outputVelocity = off [default]
> colvars:   # outputSystemForce = off [default]
> colvars:   # outputAppliedForce = off [default]
> colvars:   # runAve = off [default]
> colvars:   # corrFunc = off [default]
> colvars:
> ----------------------------------------------------------------------
> colvars:   Initializing a new collective variable.
> colvars:   # name = gyration
> colvars:   Initializing a new "radius of gyration" component.
> colvars:     # componentCoeff = 1 [default]
> colvars:     # componentExp = 1 [default]
> colvars:       Initializing atom group "atoms".
> colvars:       Atom group "atoms" defined, 15 atoms, total mass = 212.181.
> colvars:   All components initialized.
> colvars:   # width = 0.01
> colvars:   # lowerBoundary = 0 [default]
> colvars:   Lower boundary was not defined.
> colvars:   # upperBoundary = 0 [default]
> colvars:   Upper boundary was not defined.
> colvars:   # extendedLagrangian = off [default]
> colvars:   # outputValue = on [default]
> colvars:   # outputVelocity = off [default]
> colvars:   # outputSystemForce = off [default]
> colvars:   # outputAppliedForce = off [default]
> colvars:   # runAve = off [default]
> colvars:   # corrFunc = off [default]
> colvars:
> ----------------------------------------------------------------------
> colvars: Collective variables initialized, 2 in total.
> colvars:
> ----------------------------------------------------------------------
> colvars:   Initializing a new "abf" instance.
> colvars:   # name = "abf1" [default]
> colvars:   # colvars = { d, gyration }
> colvars:   # applybias = on
> colvars:   # hidejacobian = off
> colvars:   Jacobian (geometric) forces will be included in reported free
> energy gradients.
> colvars:   # fullsamples = 500
> colvars:   # inputprefix =  [default]
> colvars:   # outputfreq = 1000
> colvars:   Tried to initialize a count grid on variable with undefined
> boundaries.
> colvars:   If this error message is unclear, try recompile with
> -DCOLVARS_DEBUG.
> FATAL ERROR: Error in the collective variables module: exiting.
>
> Input:
>
> colvarsTrajFrequency    1
> colvarsRestartFrequency 10000
> analysis    on
>
>
> colvar {
>
>    name d
>
>    lowerBoundary  3.35
>    upperBoundary 9.70
>    width 0.01
>
>    distance {
>        group1  {
>           atomnumbers 1 3 5
>        }
>        group2  {
>           atomnumbers 16 17 18
>        }
>    }
> }
>
> colvar {
>    name gyration
>    width 0.01
>
>    gyration {
>        atoms  {
>           atomNumbers 7 8 9 10 11 12 13 14 15 19 20 21 22 23 24
>        }
>      }
> }
>
> abf {
>    colvars d gyration
>    fullSamples    500
>    hideJacobian    no
>    outputFreq    1000
>    applyBias    yes
> }
>
> ________________________________
>
> No virus found in this incoming message.
> Checked by AVG - www.avg.com
> Version: 9.0.722 / Virus Database: 270.14.123/2594 - Release Date: 12/30/09
> 08:27:00
>
>
>

This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:53:39 CST