Hi,<div><br></div><div>I asked this question before. </div><div>Thanks to Justin, he did reply as follows. </div><div><br></div><div><pre style="margin:0em">Kiwoong Kim wrote:
</pre><blockquote style="margin:0em;padding-left:0.85em;border-left-color:rgb(85,85,238);border-left-width:0.2em;border-left-style:solid"><pre style="margin:0em">Hi,
I'm getting better to use Gromacs owing to many posts on this sites :) .
I have several questions about adding the counter-ions.
</pre><tt>My system has a number of N2 molecules which has charge, -0.40484(for </tt><tt>single N) X 2.</tt><pre style="margin:0em">Thus, I have to add some counter-ions to make the system neutral.
</pre><tt>However, because I'm newbie on Gromacs, I thought of several clumsy ways </tt><tt>myself.</tt><pre style="margin:0em">#1. Add virtual sites (virtual atoms) which has counter-ions like below.
[ atomtypes ]
; name mass charge ptype sigma epsilon
DUM 0 0.80968 V 0.0 0.0
</pre><tt>I set the coordinates of each virtual DUM atoms to the center of N2 </tt><tt>molecules.</tt><pre style="margin:0em">#2. using genion in Gromacs.
</pre><tt>But, I have no idea on this. What molecules do I have to designate to </tt><tt>charge plus ion using genion ??</tt><pre style="margin:0em">I typed below line.
</pre><tt>genion -s md_0_1.tpr -n index.ndx -o ton_genion.pdb -g md_0_1_genion.log </tt><tt>-p topol_genion.top -np 322 -pname dum -pq 0.80968</tt><pre style="margin:0em">and selected N2 molecules which is diffusing particle.
</pre><tt>It results that the name of half of N2 molecules is changed as DUM. </tt><tt>(maybe the system become neutral)</tt><pre style="margin:0em">The number of N2 molecules should be fixed.
</pre><tt>Do I have additional N2 molecule for charging using genion?? Then, what </tt><tt>are the initial coordinates???</tt><pre style="margin:0em">please help me....
any advises would be helpful.
</pre><tt>How can I do that??? </tt><tt></tt></blockquote><pre style="margin:0em"></pre><tt>I see no reason why you should do either. For dinitrogen, which has no net </tt><tt>dipole, it seems intuitive to me that both N atoms should have zero charge.</tt><pre style="margin:0em">
-Justin
</pre><pre style="margin:0em"><span style="font-family:arial;white-space:normal">I made linear rigid model (CO2) by introducing virtual interaction. However, I have troubles in modelling N2 molecules.</span></pre><pre style="margin:0em">
<span style="font-family:arial;white-space:normal">I want to model N2 with quadrupole moment according to the paper (Adsorption (2007) 13: 469-476).</span></pre><pre style="margin:0em"><span style="font-family:arial;white-space:normal">In the paper, the charge N2 model (originally developed by Murthy et al. (1980) can be described as follows.</span></pre>
<pre style="margin:0em"><span style="font-family:arial;white-space:normal">N (-0.40484), Dummy (N2, +0.80968).</span></pre><pre style="margin:0em"><span style="font-family:arial;white-space:normal"><br></span></pre><pre style="margin:0em">
<font face="arial"><span style="white-space:normal">Actually, I don't know how to model N2 molecules with dummy. Please let me know any advises. </span></font></pre><pre style="margin:0em"><font face="arial"><span style="white-space:normal"><br>
</span></font></pre><pre style="margin:0em"><font face="arial"><span style="white-space:normal">Thx. Have a great day. :)</span></font></pre><pre style="margin:0em"></pre><pre style="margin:0em"><br></pre></div>