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

Box and atom coordinate scaling not consistent when converting from GROMACS to LAMMPS #2669

Open
chw90 opened this issue Jan 31, 2024 · 1 comment

Comments

@chw90
Copy link

chw90 commented Jan 31, 2024

Environment Information

Open Babel version: 3.1.1
Operating system and version: Ubuntu 22.04

Expected Behavior

The coordinates of atoms and the coordinates of the simulation box bounds are either unchanged or at least scaled with the same factor.

Actual Behavior

Atom coordinates are scaled with a factor of 10. Box bounds in the LAMMPS data file seem uncorrelated with the box vectors given in the GROMACS data file.

Steps to Reproduce

  • An example GROMACS input file is attached (file extension renamed to .txt so GitHub accepts it).
$ mv ice1htip3p.txt ice1htip3p.gro
$ obabel ice1htip3p.gro -olmpdat > ice1htip3p.lmpdat

Now compare the GROMACS box vectors (last line of the `.gro file) and the LAMMPS box sizes:

$ tail -n 1 ice1htip3p.gro
    2.34685163 2.20607179 2.71095555
$ cat ice1htip3p.lmpdat | grep -e ".lo" -e ".hi"
  -0.59000   27.21000 xlo xhi
  -0.59000   27.21000 ylo yhi
  -0.59000   27.21000 zlo zhi

Expected result would be:

$ cat ice1htip3p.lmpdat | grep -e ".lo" -e ".hi"
   0.00000   23.4685163 xlo xhi
   0.00000   22.0607179 ylo yhi
   0.00000   27.1095555 zlo zhi

Is this expected behavior? Does Open Babel compute the simulation box itself based on the atoms instead of taking the definition of the GROMACS input?

ice1htip3p.txt

Copy link

welcome bot commented Jan 31, 2024

Thanks for opening your first issue here! Be sure to follow the issue template!

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