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

Conversion DESMOND--> AMBER #365

Open
icamps opened this issue Aug 18, 2020 · 5 comments
Open

Conversion DESMOND--> AMBER #365

icamps opened this issue Aug 18, 2020 · 5 comments

Comments

@icamps
Copy link

icamps commented Aug 18, 2020

I would like to know how to make the conversion DESMOND-->AMBER.

From here, it is tell that the conversion will be DESMOND--> GROMACS --> AMBER.

Using the command:

python convert.py --des_in ~/Documents/01_Sim1-ARB/desmond_md_Sim1-ARB-out.cms --amber --odir ~/Documents/AMBERFormat/01_Sim1-ARB/

I am getting the message:

ERROR 2020-08-18 16:17:04 list index out of range
Traceback (most recent call last):
  File "convert.py", line 663, in _save_amber
    top = pmd.load_file(oname + '.top', xyz=oname + '.gro')
  File "/home/icamps/anaconda3/lib/python3.7/site-packages/parmed/formats/registry.py", line 214, in load_file
    return cls(filename, *args, **kwargs)
  File "/home/icamps/anaconda3/lib/python3.7/site-packages/parmed/gromacs/gromacstop.py", line 251, in __init__
    self.read(fname, defines, parametrize)
  File "/home/icamps/anaconda3/lib/python3.7/site-packages/parmed/gromacs/gromacstop.py", line 374, in read
    atoms = [molecule.atoms[int(w)-1] for w in line.split()]
  File "/home/icamps/anaconda3/lib/python3.7/site-packages/parmed/gromacs/gromacstop.py", line 374, in <listcomp>
    atoms = [molecule.atoms[int(w)-1] for w in line.split()]
IndexError: list index out of range
INFO 2020-08-18 16:17:04 Finished!

and two files are created:

desmond_md_Sim1-ARB-out_converted.top
desmond_md_Sim1-ARB-out_converted.gro

What I am doing wrong? I appreciate any help.

@mrshirts
Copy link
Contributor

So, it's converting Desmond into GROMACS correctly, but not from GROMACS into AMBER. So there is probably something wrong with the GROMACS files. It's impossible to tell from here, though. How many atoms are the .cms files?

@icamps
Copy link
Author

icamps commented Aug 19, 2020

In the .cms file are 89531 atoms.

@icamps
Copy link
Author

icamps commented Aug 31, 2020

Any ideas?

@mrshirts
Copy link
Contributor

mrshirts commented Sep 1, 2020 via email

@icamps
Copy link
Author

icamps commented Sep 2, 2020

The files are too big (2.4GB) to share here. Please, you can download the DESMOND files from https://bit.ly/3lBOYzY.

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

2 participants