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 resulting .xtc file from unigbsa-md had not water molecules and ions. #58

Open
Tangbbmc opened this issue Oct 19, 2024 · 4 comments

Comments

@Tangbbmc
Copy link

The resulting .xtc file from unigbsa-md had not water molecules and ions. How to include water molecules and ions in .xtc file. Moreover, when processing pbc with unigbsa-pbc command, I received the following error:

(gbsa) dell@dell-OptiPlex-Tower-Plus-7010:~/Downloads/Uni-GBSA-master/example/1ceb$ unigbsa-pbc -s md.tpr -f md-pbc.xtc -o md-nopbc.xtc -n index.ndx
gmx trjconv -s md.tpr -n index.ndx -f md-pbc.xtc -o md-pbc-1.xtc -pbc mol -ur compact >> /dev/null 2>&1 << EOF
output
EOF
Traceback (most recent call last):
File "/home/dell/anaconda3/envs/gbsa/bin/unigbsa-pbc", line 8, in
sys.exit(PBC_main())
^^^^^^^^^^
File "/home/dell/anaconda3/envs/gbsa/lib/python3.12/site-packages/unigbsa/CLI.py", line 31, in PBC_main
process_pbc(trajfile, tprfile, indexfile=indexfile, outfile=outfile)
File "/home/dell/anaconda3/envs/gbsa/lib/python3.12/site-packages/unigbsa/utils.py", line 117, in process_pbc
raise Exception('Process pbc error: see the %s file for details'%logfile)
Exception: Process pbc error: see the /dev/null file for details

@Aunity
Copy link
Collaborator

Aunity commented Oct 25, 2024

You can save all files during the MD process by adding the -verbose parameter. And the outdir/complex.GMX/md.xtc is the trajectory contains the water molecules and ions. As for the error you mentioned, could you please provide all the files you run(md.tpr, md-pbc.xtc and index.ndx)?

@Tangbbmc
Copy link
Author

Tangbbmc commented Oct 25, 2024 via email

@Aunity
Copy link
Collaborator

Aunity commented Oct 25, 2024

The md.xtc file is the original output trajectory file of the md production. The md-pbc.xtc and the traj_comx.xtc file are the same file in different floder and name. The md-pbc.xtc file is the md.xtc file with PBC removed, and the traj_comx.xtc is a copy of md-pbc.xtc file.

@Tangbbmc
Copy link
Author

Tangbbmc commented Oct 25, 2024 via email

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