Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

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

Open
vemaparna opened this issue Sep 26, 2022 · 0 comments

Comments

@vemaparna
Copy link

vemaparna commented Sep 26, 2022

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>>
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant