Re: Definition of two colvars for ABF

From: Jérôme Hénin (jhenin_at_ifr88.cnrs-mrs.fr)
Date: Tue Jan 05 2010 - 08:12:47 CST

Keeping the mailing list in CC, since others might have ideas about
what is happening. At this point, I don't, unfortunately.

2010/1/5 Branko <bdrakuli_at_chem.bg.ac.yu>:
> Dear Jerome,
>
> Thank you very much for the answer.For windows  execute abf_integrate, I'll
> must ask colleagues skilled in programming, but I try to reproduce (short)
> simulation using just one colvar for simplicity, and using the low
> colvarsRestartFrequency 100 in place of 10000, but <files>.state still did
> not appears. In namd out file, again there is data that those files should
> be generated (please see bellow, underlined0. If this can be helpful part of
> colvar definition file:
>
> colvarsTrajFrequency    1  (same as in previous simulations)
> colvarsRestartFrequency 100 (lower)
> analysis    on (same as in previous simulations)
>
> from NAMD input:
>
> restartfreq         1000  (same as in previous simulations)
>
> and from NAMD output (part that repeat in every iteration, i.e. recorded
> step):
>
> WRITING EXTENDED SYSTEM TO RESTART FILE AT STEP 24000
> WRITING COORDINATES TO DCD FILE AT STEP 24000
> WRITING COORDINATES TO RESTART FILE AT STEP 24000
> FINISHED WRITING RESTART COORDINATES
> WRITING VELOCITIES TO RESTART FILE AT STEP 24000
> FINISHED WRITING RESTART VELOCITIES
> colvars: Synchronizing trajectory file.
> colvars: Synchronizing trajectory file.
> colvars: Synchronizing trajectory file.
> colvars: Synchronizing trajectory file.
> colvars: Synchronizing trajectory file.
> colvars: Synchronizing trajectory file.
> colvars: Synchronizing trajectory file.
> colvars: Synchronizing trajectory file.
> colvars: Synchronizing trajectory file.
> colvars: Synchronizing trajectory file.
> ENERGY:   25000      1774.3401      2569.9087       590.3621
> 0.0705            340.3763     -2520.7006         0.0000         0.0000
> 5115.0868           7869.4439       305.2329      2754.3571
> 7919.9032       309.9006             66.6905       168.5968
> 120000.0000       -75.6140       -75.9222
>
> This "colvars: Synchronizing trajectory file" didn't appear if
> colvarsRestartFrequency is 10000.
> Hope that this could be helpful.
>
> Sincerely
>
> Branko
>
>
>
>
>
> colvars:
> ----------------------------------------------------------------------
> colvars: Initializing the collective variables module, version 29-09-2008.
> colvars: # colvarsTrajFrequency = 1
> colvars: # colvarsRestartFrequency = 100
> 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: Collective variables initialized, 1 in total.
> colvars:
> ----------------------------------------------------------------------
> colvars:   Initializing a new "abf" instance.
> colvars:   # name = "abf1" [default]
> colvars:   # colvars = { d }
> 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:   Finished ABF setup.
> colvars:
> ----------------------------------------------------------------------
> colvars: Collective variables biases initialized, 1 in total.
> colvars:
> ----------------------------------------------------------------------
> colvars: Collective variables module initialized.
> colvars:
> ----------------------------------------------------------------------
>
>
>
>
> Jérôme Hénin wrote:
>
> 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
>
>
>
>
>
>
>
> ________________________________
>
> No virus found in this incoming message.
> Checked by AVG - www.avg.com
>
> Version: 9.0.725 / Virus Database: 270.14.125/2600 - Release Date: 01/04/10
> 20:35:00
>
>
>

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