<div dir="ltr">Hi,<div><br></div><div>Thanks for your interest! On this list we try to discuss issues about how to develop with GROMACS or for GROMACS. If there's a known bug, as it looks like you have, then opening at issue at <a href="https://redmine.gromacs.org">https://redmine.gromacs.org</a> would be a very welcome start (particularly mentioning details of which GROMACS version and attaching tarballs input files).</div><div><br></div><div>Generating initial coordinate files is a better question for the gmx-users mailing list. What tools you might want depend on lot on what you're trying to do. While there's a partial list at <a href="http://www.gromacs.org/Downloads/Related_Software">http://www.gromacs.org/Downloads/Related_Software</a>, others doing projects similar to yours will tend to have better suggestions than the developers!</div><div><br></div><div>Mark</div></div><br><div class="gmail_quote"><div dir="ltr">On Fri, May 25, 2018 at 4:09 PM Marc Barbry <<a href="mailto:marc.barbry@mailoo.org">marc.barbry@mailoo.org</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div text="#000000" bgcolor="#FFFFFF">
Dear Gromax developers,<br>
<br>
First of all, I would like to indicate that I'm a new user of
Gromacs, therefore some of my comments/questions may be basics.
Furthermore, it seems that I'm not using the right method to do what
is explained below, if you could indicate me a better technique to
realize this it would greatly help me.<br>
<br>
<br>
Running the gmx solvant program, I detected a bug related to the
input. When using a basic pdb file for water containing only the
coordinates of the atoms,<br>
<br>
MODEL 1<br>
ATOM 0 O MOL 1 0.000 0.000 0.119 1.00
0.00 O <br>
ATOM 1 H MOL 1 0.000 0.763 -0.477 1.00
0.00 H <br>
ATOM 2 H MOL 1 0.000 -0.763 -0.477 1.00
0.00 H <br>
ENDMDL<br>
<br>
Then the program was stuck in a while loop, because the box in
solvate.cpp was null.<br>
<br>
<br>
In more details, I'm trying to do the tutorial present in this page<br>
<a href="http://www.gromacs.org/@api/deki/files/198/=gmx-tutorial.pdf" target="_blank">http://www.gromacs.org/@api/deki/files/198/=gmx-tutorial.pdf</a><br>
<br>
You can download the input files with the following link<br>
<a href="https://mbarbry.website.fr.to/nextcloud/index.php/s/Wts5ejbyp76ymDp" target="_blank">https://mbarbry.website.fr.to/nextcloud/index.php/s/Wts5ejbyp76ymDp</a><br>
<br>
Thus I have a Lysozyme protein, and I want to add a solvent around
it, therefore I used the command.<br>
<br>
gmx solvate -cs H2O-nocell.pdb -cp centered.gro -p topol.top
-o solvated.gro<br>
<br>
Where H2O-nocell.pdb is the file described above.<br>
In this case the program get stuck in the while loop line ~250 of
solvant.cpp because box[d][d] = 0.0. Obviously, a check of the box
is missing here in order to avoid such kind of issue, which are
difficult to debug.<br>
<br>
Second problem, I then added a cell to the H2O.pdb file,<br>
<br>
MODEL 1<br>
CRYST1 8.0 10.0 8.0 90.00 90.00 90.00<br>
ORIGX1 1.000000 0.000000 0.000000 0.00000
<br>
ORIGX2 0.000000 1.000000 0.000000 0.00000
<br>
ORIGX3 0.000000 0.000000 1.000000 0.00000<br>
ATOM 0 O MOL 1 0.000 0.000 0.119 1.00
0.00 O <br>
ATOM 1 H MOL 1 0.000 0.763 -0.477 1.00
0.00 H <br>
ATOM 2 H MOL 1 0.000 -0.763 -0.477 1.00
0.00 H <br>
ENDMDL<br>
<br>
<br>
The program goes a bit further but then return a memory error<br>
<br>
Removed 0 solvent atoms due to solvent-solvent overlap<br>
Removed 120 solvent atoms due to solute-solvent overlap<br>
<br>
Program received signal SIGSEGV, Segmentation fault.<br>
0x00007ffff7027dd7 in
gmx::internal::AnalysisNeighborhoodSearchImpl::~AnalysisNeighborhoodSearchImpl()<br>
() from /home/marc/.local/gromacs/lib/libgromacs.so.3<br>
(gdb) quit<br>
<br>
<br>
I guess my input file is not correctly constructed, I'm a newbie
with pdb files and MD in general. Do you know a good way to generate
pdb file? At the moment I use <a href="http://avogadro.cc/" target="_blank">Avogadro</a> and <a href="https://wiki.fysik.dtu.dk/ase/index.html" target="_blank">ASE</a>, but
obviously they were not think to be use with Gromacs.<br>
<br>
Thanks for your helps,<br>
<br>
Best regards,<br>
Marc Barbry<br>
</div>
-- <br>
Gromacs Developers mailing list<br>
<br>
* Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/GMX-developers_List" rel="noreferrer" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/GMX-developers_List</a> before posting!<br>
<br>
* Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" rel="noreferrer" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br>
<br>
* For (un)subscribe requests visit<br>
<a href="https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-developers" rel="noreferrer" target="_blank">https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-developers</a> or send a mail to <a href="mailto:gmx-developers-request@gromacs.org" target="_blank">gmx-developers-request@gromacs.org</a>.</blockquote></div>