<div dir="ltr">Hi,<div><br></div><div>Please ask questions related to GROMACS use on the users mailing list. This channel is for code development.</div><div><br></div><div>Mark</div></div><br><div class="gmail_quote"><div dir="ltr">On Wed, Oct 14, 2015 at 2:05 PM toannt &lt;<a href="mailto:toannt@vnu.edu.vn">toannt@vnu.edu.vn</a>&gt; wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br>
<br>
Thanks. I&#39;m defining two molecules types, each type have only one<br>
molecule in the simulation. When I said custom bonds, I was referring to<br>
the hydrogen bonds between the two strands. I&#39;d like to understand how<br>
these bonds are handled by Gromacs in an atomistic simulation? If I have<br>
two strands at some distance apart, do they naturally come to each other<br>
to form the dsDNA?<br>
<br>
By the way, I am doing a grand canonical simulation, but I&#39;d like to<br>
reuse libgromacs to avoid reinventing the wheel with some functions.<br>
<br>
Toan<br>
<br>
Vào ngày 14-10-2015 17:56, Berk Hess viết:<br>
&gt; Hi,<br>
&gt;<br>
&gt; There is no direct link between &quot;real&quot; molecules and molecule types.<br>
&gt; The only restriction is that you can only define chemical bonds<br>
&gt; between atoms within the same molecule type. Thus you can have two<br>
&gt; DNA<br>
&gt; strands in one or in two separate molecule types.<br>
&gt; I don&#39;t understand what custom-bonds you are referring to. Usually<br>
&gt; DNA strands are not chemically bonded to each other.<br>
&gt;<br>
&gt; Cheers,<br>
&gt;<br>
&gt; Berk<br>
&gt;<br>
&gt; On 10/14/2015 11:49 AM, toannt wrote:<br>
&gt;&gt;<br>
&gt;&gt; Hi all<br>
&gt;&gt;<br>
&gt;&gt; I&#39;ve just started to look at libgromacs for a custom simulation. I&#39;m<br>
&gt;&gt; studying the gmx data structures to initialize them directly. I wonder<br>
&gt;&gt; if molecules such as DNA is considered by gromacs to be two molecule<br>
&gt;&gt; types (one for each strand) or is only one molecule type with some<br>
&gt;&gt; custom-bonds between the two strands?<br>
&gt;&gt;<br>
&gt;&gt; Thanks<br>
&gt;&gt;<br>
&gt;&gt; Toan<br>
--<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&#39;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>