Re: Restart multiple walker ABF

From: Jeff Comer (jeffcomer_at_gmail.com)
Date: Thu May 19 2016 - 15:58:22 CDT

You have 4 replicas that write four different restart files indexed by
[myReplica]. When you restart, each system should restart from its restart
file:
binCoordinates $inputname.[myReplica].coor
binVelocities $inputname.[myReplica].vel
extendedSystem $inputname.[myReplica].xsc

Jeff

–––––––––––––––––––––––––––––––––––———————
Jeffrey Comer, PhD
Assistant Professor
Institute of Computational Comparative Medicine
Nanotechnology Innovation Center of Kansas State
Kansas State University
Office: P-213 Mosier Hall
Phone: 785-532-6311
Website: http://jeffcomer.us

On Thu, May 19, 2016 at 3:32 AM, Atanu Maity <atanuchem48_at_gmail.com> wrote:

> Dear users,
>
> I am facing this problem while trying to restart an REMD (multiple walker)
> run as a part of ABF simulation. I start the simulation using the following
> command line –
>
>
> mpirun -np 32 namd2 +replicas 4 abf-replica.namd +stdout
> replica4/%d/job4.%d.log &
>
>
> Four replicas were running using 8 core each properly.
>
>
> The input contains the following lines –
>
>
>
> structure step5_assembly.xplor_ext.psf
>
> coordinates step5_assembly.pdb
>
> set temp 303.15;
>
> set outputname replica4/abf-22-prod2.[myReplica]
>
> set inputname abf-22-prod1.1
>
> binCoordinates $inputname.coor
>
> binVelocities $inputname.vel
>
> extendedSystem $inputname.xsc
>
> outputname $outputname
>
> restartfreq 10000;
>
> dcdfreq 2000;
>
> dcdUnitCell yes;
>
> xstFreq 10000;
>
> outputEnergies 10000;
>
> outputTiming 10000;
>
>
> # Force-Field Parameters
>
> paraTypeCharmm on;
>
> parameters toppar-namd/par_all36_prot.prm;
>
> parameters toppar-namd/par_all36_lipid.prm;
>
> parameters toppar-namd/toppar_all36_lipid_cholesterol.str;
>
> parameters toppar-namd/toppar_all36_lipid_glycolipid.str;
>
> parameters toppar-namd/toppar_all36_lipid_bacterial.str;
>
> # These are specified by CHARMM
>
> exclude scaled1-4
>
> langevinPistonTemp $temp;
>
>
>
> # Constant Temperature Control
>
> langevin on; # langevin dynamics
>
> langevinDamping 1.0; # damping coefficient of 1/ps (keep
> low)
>
> langevinTemp $temp; # random noise at this level
>
> langevinHydrogen off; # don't couple bath to hydrogens
>
>
>
> colvars on
>
> colvarsConfig abf-22-bb-replica.in
>
> ## Replica Exchange ##
>
> ##SElectionRules will be applied at each 100 steps###
>
> source /apps/NAMD_2.11b1_Source/lib/selectionRules.tcl
>
> source /apps/NAMD_2.11b1_Source/lib/resampleWalkers.tcl
>
> source /apps/NAMD_2.11b1_Source/lib/minExchanges.tcl
>
> firsttimestep 230000
>
> replicaUniformPatchGrids on
>
> set n 12500
>
> set sharedFreq 100
>
> for {set i 0} {$i < $n} {incr i} {
>
> run $sharedFreq
>
> cv bias abf1 share
>
> if {$i % 50 == 49} {
>
> selectionRules
>
> }
>
> }
>
>
>
> The problem with the restart is that I did not understand which restart
> files (.coor, .vel, .state) should I use because I have four sets of such
> files.
>

This archive was generated by hypermail 2.1.6 : Tue Dec 27 2016 - 23:22:09 CST