def run_rfe():
pdb = AbFile('systems/ptp1b-l1-l2-complex.pdb', tag='pdb')
top = AbFile('systems/ptp1b-l1-l2-complex.top', tag='topology')
tag = AbFile('systems/ptp1b-l1-l2-tags.pdb', tag='alchemicaltags')
cor = AbFile('systems/ptp1b-l1-l2-complex.inpcrd', tag='coordinate')
system = System(name='ptp1b-l1-l2', files=[pdb, top, tag, cor])
p = Protocol(clone_settings=False)
for step, numsteps in zip(Rfe.steps, [5000, 50000]):
rfe = Simulation()
rfe.system = system
rfe.engine = 'namd_mpi'
rfe.cores = 32
rfe.cutoff = 12.0
rfe.switchdist = 10.0
rfe.pairlistdist = 13.5
rfe.numminsteps = 5000
rfe.numsteps = numsteps
rfe.add_input_file(step, is_executable_argument=True)
rfe.add_ensemble('replica', range(1))
# to increase the number of EnTK tasks: change the lambdawindow parameter
rfe.add_ensemble('lambdawindow', [1.])
p.append(rfe)
ht = Runner('bw_aprun', comm_server=('two.radical-project.org', 33158))
ht.add_protocol(p)
ht.run(walltime=480, queue='high')
def run_afe():
system = System(prefix='systems/dov-wt')
afe = Simulation()
afe.system = system
afe.engine = 'namd'
afe.cores = 2
afe.numminsteps = 10
afe.numsteps = 10
afe.add_input_file(Afe.step0, is_executable_argument=True)
afe.add_input_file('inputs/restraint.in', is_executable_argument=False)
afe.add_ensemble('replica', range(2))
afe.add_ensemble('lambdawindow', [0.0, 1.0])
afe.cutoff = 12.0
afe.switchdist = 10.0
afe.pairlistdist = 13.5
afe.k1 = 10
afe.k2 = 500
ht = Runner(comm_server=('two.radical-project.org', 33146))
ht.add_protocol(afe)
ht.run(walltime=1000)