<br><br>Hi<br>molecule dipole is 48.0 sum of q_i x_i<br><br>based on the following two websites, <br><i><b>x_i </b></i> is the <a href="http://en.wikipedia.org/wiki/Displacement_%28vector%29" title="Displacement (vector)">displacement vector</a> pointing from the negative charge to the positive charge. <br>
<br>what about the x_i for the salt-molecule, which dissociates into one counter ion and the rest of the molecule in water?<br><br><br><a href="http://en.wikipedia.org/wiki/Bond_dipole_moment">http://en.wikipedia.org/wiki/Bond_dipole_moment</a><br>
<a href="http://en.wikipedia.org/wiki/Electric_dipole_moment">http://en.wikipedia.org/wiki/Electric_dipole_moment</a><br><br><br>Thank you<br>Lin<br><br><br><br>On 2010-10-20 06.06, Chih-Ying Lin wrote:<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; Hi<br>
&gt; molecule dipole is 48.0 sum of q_i x_i<br>
&gt; x_i the bond length for covalent bond.<br>
<br>
No. x_i is the atomic position.<br>
<br>
&gt; but what is &quot;x_i&quot; for salt-molecule?<br>
&gt;<br>
&gt;<br>
&gt; For salt-molecule, the ionic bonds are broken in water solvent and the<br>
&gt; counter ions are spread among the water.<br>
&gt;<br>
&gt; What is the &quot;x_i&quot; of the ionic bond in the dipole moment calculation?<br>
&gt; Is x_i equal to the distance of the two parts of the salt-molecules (the<br>
&gt; counter ion and the rest of the molecule) even though the salt molecule<br>
&gt; has dissolved in the water?<br>
&gt;<br>
&gt; I mean, is x_i equal to the length of simulation box if the counter ion<br>
&gt; and the rest of the molecule are in the two sides of the simulation box?<br>
&gt;<br>
&gt;<br>
&gt; I mean, if Gromacs takes x_i as the length of simulation box if the<br>
&gt; counter ion and the rest of the molecule are in the two sides of the<br>
&gt; simulation box?<br>
&gt;<br>
&gt; Thank you<br>
&gt; Lin<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt;<br>
&gt; Try <a href="http://en.wikipedia.org/wiki/Electric_dipole_moment" target="_blank">http://en.wikipedia.org/wiki/Electric_dipole_moment</a> , you might want<br>
&gt; to read the part about calculating dipole moments for an array of point<br>
&gt; charges, it is not difficult. 33 point charges are doable using pencil<br>
&gt; and calculator in about 10min. Do not worry about the reference point as<br>
&gt; long as your system is neutral, just set it to (0,0,0). Otherwise, take<br>
&gt; any kind of first year physics book it will contain very similar<br>
&gt; information.<br>
&gt;<br>
&gt; On 10/19/2010 05:39 AM, Chih-Ying Lin wrote:<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;  Hi<br>
&gt;&gt;  According to the following website,<br>
&gt;&gt;<br>
&gt;&gt;  <a href="http://en.wikipedia.org/wiki/Bond_dipole_moment" target="_blank">http://en.wikipedia.org/wiki/Bond_dipole_moment</a><br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;  \mu = \delta \, d.<br>
&gt;&gt;  The bond dipole is modeled as +ä -- ä- with a distance /d/ between the<br>
&gt;&gt;  partial charges &lt;<a href="http://en.wikipedia.org/wiki/Partial_charges" target="_blank">http://en.wikipedia.org/wiki/Partial_charges</a>&gt; +ä and ä-.<br>
&gt;&gt;  For a complete molecule the total molecular dipole moment may be<br>
&gt;&gt;  approximated as the vector sum of individual bond dipole moments.<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;  However, for a molecule of multiple atoms,<br>
&gt;&gt;  There may be more than one bond connected on one atom.<br>
&gt;&gt;  E<br>
&gt;&gt;  |<br>
&gt;&gt;  B - A - C<br>
&gt;&gt;  partial charge of atom_A = -0.5<br>
&gt;&gt;  partial charge of atom_B = 0.2<br>
&gt;&gt;  partial charge of atom_C = 0.35<br>
&gt;&gt;  partial charge of atom_E = 0.4<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;  Which partial charges should I use when I calculate bond-dipole-moment<br>
&gt;&gt;  of A-B ?<br>
&gt;&gt;  Which partial charges should I use when I calculate bond-dipole-moment<br>
&gt;&gt;  of A-C ?<br>
&gt;&gt;  Which partial charges should I use when I calculate bond-dipole-moment<br>
&gt;&gt;  of A-E ?<br>
&gt;&gt;<br>
&gt;&gt;  Thank you<br>
&gt;&gt;  Lin<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;<br>
&gt;&gt;  On 2010-10-18 03.30, Chih-Ying Lin wrote:<br>
&gt;&gt; &gt;  HI<br>
&gt;&gt; &gt;  I confined one molecule in the center of box and issue the g_dipole<br>
&gt;&gt;  command.<br>
&gt;&gt; &gt;  The average dipole moment is still around 32.<br>
&gt;&gt; &gt;  It is the molecule with 33 atoms / united atoms of most carbon groups,<br>
&gt;&gt; &gt;  isn&#39;t the dipole moment around 32 too high?<br>
&gt;&gt; &gt;  How can I test next and know that the dipole moment around 32 is<br>
&gt;&gt; &gt;  acceptable?<br>
&gt;&gt;  By calculating on paper: dipole is 48.0 sum of q_i x_i, therefore if you<br>
&gt;&gt;  have large charge separation you will get a large dipole.<br>
&gt;&gt;<br>
&gt;&gt; &gt;  Thank you<br>
&gt;&gt; &gt;  Lin<br>
&gt;&gt; &gt;  On 2010-10-16 21.36, Chih-Ying Lin wrote:<br>
&gt;&gt; &gt; &gt;<br>
&gt;&gt; &gt; &gt; Hi<br>
&gt;&gt; &gt; &gt; I issue the g_dipole command on Gromacs =&gt; And, the following<br>
&gt;&gt; &gt; &gt; information is shown.<br>
&gt;&gt; &gt; &gt; There are 10 molecules in the selection,<br>
&gt;&gt; &gt; &gt; Does the Average =32.1611 refer to the average for a single over the<br>
&gt;&gt; &gt; &gt; simulation time?<br>
&gt;&gt; &gt; &gt; Or, the Average = 32.1611 summing for all the 10 molecules over the<br>
&gt;&gt; &gt; &gt; simulation time?<br>
&gt;&gt; &gt; &gt; If the average = 32.1611 for a single molecule with 33 atoms / united<br>
&gt;&gt; &gt; &gt; atoms of most carbon groups, isn&#39;t the dipole moment too high?<br>
&gt;&gt; &gt;  I think this is the average per molecule. You may need to run trjconv<br>
&gt;&gt; &gt;  -pbc whole, because mdrun may break molecules in two parts, meaning that<br>
&gt;&gt; &gt;  the molecule becomes as big as the box.<br>
&gt;&gt; &gt;<br>
&gt;&gt; &gt; &gt;<br>
&gt;&gt; &gt; &gt;<br>
&gt;&gt; &gt; &gt;<br>
&gt;&gt; &gt; &gt;<br>
&gt;&gt; &gt; &gt; What does &quot;will subtract their charge at their center of mass&quot; this<br>
&gt;&gt; &gt;  mean?<br>
&gt;&gt; &gt; &gt; Why &quot;will subtract their charge at their center of mass&quot; ?<br>
&gt;&gt; &gt; &gt;<br>
&gt;&gt; &gt; &gt;<br>
&gt;&gt; &gt; &gt;<br>
&gt;&gt;<br>
&gt; -------------- next part --------------<br>
&gt; An HTML attachment was scrubbed...<br>
&gt; URL:<br>
&gt; <a href="http://lists.gromacs.org/pipermail/gmx-users/attachments/20101019/4cf32833/attachment.html" target="_blank">http://lists.gromacs.org/pipermail/gmx-users/attachments/20101019/4cf32833/attachment.html</a><br>

&gt;<br>
<br>
<br>
--<br>
David van der Spoel, Ph.D., Professor of Biology<br>
Dept. of Cell &amp; Molec. Biol., Uppsala University.<br>
Box 596, 75124 Uppsala, Sweden. Phone:  +46184714205.<br>
<a href="mailto:spoel@xray.bmc.uu.se">spoel@xray.bmc.uu.se</a>    <a href="http://folding.bmc.uu.se/" target="_blank">http://folding.bmc.uu.se</a><br>
<br>
<br>
------------------------------<br>
<font color="#888888"><br>
--<br>
gmx-users mailing list<br>
<a href="mailto:gmx-users@gromacs.org">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/Search" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/Search</a> before posting!<br>
<br>
End of gmx-users Digest, Vol 78, Issue 142<br>
******************************************</font><br>