<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Thanks for the links, I will have a look.<br>
<br>
Marc<br>
<br>
<div class="moz-cite-prefix">On 25-05-18 16:15, Mark Abraham wrote:<br>
</div>
<blockquote type="cite"
cite="mid:CAMNuMAT1pbtHHxHX2z6TOfKrnXu9BCs_NprbRBQqkwA_vPNHbg@mail.gmail.com">
<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"
moz-do-not-send="true">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"
moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">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"
moz-do-not-send="true">Avogadro</a> and <a
href="https://wiki.fysik.dtu.dk/ase/index.html"
target="_blank" moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">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" moz-do-not-send="true">gmx-developers-request@gromacs.org</a>.</blockquote>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
</blockquote>
<br>
</body>
</html>