yank icon indicating copy to clipboard operation
yank copied to clipboard

The periodic box size has decreased to less than twice the nonbonded cutoff

Open vemaparna opened this issue 2 years ago • 0 comments

I just changed the input files in the yank default yaml script. After running yank using amber prmtop and inpcrd files for 25 iterations , i got the following error. Unable to understand how the periodic box size decreased? Can you help me to figure out the problem

2022-09-25 21:03:01,339: Iteration 25/1000
2022-09-25 21:03:01,339: ********************************************************************************
2022-09-25 21:03:01,339: Single node: executing <function ReplicaExchangeSampler._mix_replicas at 0x7fd0af2cc0d0>
2022-09-25 21:03:01,339: Mixing replicas...
2022-09-25 21:03:01,347: Mixing of replicas took    0.008s
2022-09-25 21:03:01,347: Accepted 32498/388962 attempted swaps (8.4%)
2022-09-25 21:03:01,347: Propagating all replicas...
2022-09-25 21:03:01,347: Running _propagate_replica serially.
2022-09-25 21:03:07,902: Running _get_replica_move_statistics serially.
2022-09-25 21:03:07,903: Propagating all replicas took    6.556s
2022-09-25 21:03:07,903: Running _compute_replica_energies serially.
Traceback (most recent call last):
  File "/home/aparna/anaconda3/envs/yank1/bin/yank", line 10, in <module>
    sys.exit(main())
  File "/home/aparna/anaconda3/envs/yank1/lib/python3.9/site-packages/yank/cli.py", line 73, in main
    dispatched = getattr(commands, command).dispatch(command_args)
  File "/home/aparna/anaconda3/envs/yank1/lib/python3.9/site-packages/yank/commands/script.py", line 155, in dispatch
    yaml_builder.run_experiments(write_status=write_status)
  File "/home/aparna/anaconda3/envs/yank1/lib/python3.9/site-packages/yank/experiment.py", line 751, in run_experiments
    completed[exp_index] = self._run_experiment(exp, write_status=write_status)
  File "/home/aparna/anaconda3/envs/yank1/lib/python3.9/site-packages/yank/experiment.py", line 3251, in _run_experiment
    built_experiment.run(n_iterations=switch_experiment_interval)
  File "/home/aparna/anaconda3/envs/yank1/lib/python3.9/site-packages/yank/experiment.py", line 427, in run
    alchemical_phase.run(n_iterations=iterations_to_run)
  File "/home/aparna/anaconda3/envs/yank1/lib/python3.9/site-packages/yank/yank.py", line 1207, in run
    self._sampler.run(n_iterations=n_iterations)
  File "/home/aparna/anaconda3/envs/yank1/lib/python3.9/site-packages/openmmtools/multistate/multistatesampler.py", line 684, in run
    self._compute_energies()
  File "/home/aparna/anaconda3/envs/yank1/lib/python3.9/site-packages/openmmtools/utils.py", line 87, in _wrapper
    return func(*args, **kwargs)
  File "/home/aparna/anaconda3/envs/yank1/lib/python3.9/site-packages/openmmtools/multistate/multistatesampler.py", line 1318, in _compute_energies
    new_energies, replica_ids = mpiplus.distribute(self._compute_replica_energies, range(self.n_replicas),
  File "/home/aparna/anaconda3/envs/yank1/lib/python3.9/site-packages/mpiplus/mpiplus.py", line 512, in distribute
    all_results = [task(job_args, *other_args, **kwargs) for job_args in distributed_args]
  File "/home/aparna/anaconda3/envs/yank1/lib/python3.9/site-packages/mpiplus/mpiplus.py", line 512, in <listcomp>
    all_results = [task(job_args, *other_args, **kwargs) for job_args in distributed_args]
  File "/home/aparna/anaconda3/envs/yank1/lib/python3.9/site-packages/openmmtools/multistate/multistatesampler.py", line 1358, in _compute_replica_energies
    compatible_energies = states.ThermodynamicState.reduced_potential_at_states(
  File "/home/aparna/anaconda3/envs/yank1/lib/python3.9/site-packages/openmmtools/states.py", line 937, in reduced_potential_at_states
    return [thermodynamic_states[0].reduced_potential(context)]
  File "/home/aparna/anaconda3/envs/yank1/lib/python3.9/site-packages/openmmtools/states.py", line 891, in reduced_potential
    openmm_state = context_state.getState(getEnergy=True)
  File "/home/aparna/anaconda3/envs/yank1/lib/python3.9/site-packages/openmm/openmm.py", line 12318, in getState
    state = _openmm.Context_getState(self, types, enforcePeriodicBox, groups_mask)
openmm.OpenMMException: The periodic box size has decreased to less than twice the nonbonded cutoff.
2022-09-25 21:03:08,386: Single node: executing <bound method MultiStateReporter.close of <openmmtools.multistate.multistatereporter.MultiStateReporter object at 0x7fd0abe1b580>>

vemaparna avatar Sep 26 '22 05:09 vemaparna