<div dir="ltr">Hello,<div><br></div><div>Does that mean the MM radii are the same for all MM_KIND's for the same element? I'm concerned because based on JCTC 1, 1176 (2005), it seems that the appropriate radius can vary somewhat substantially for the same element in different bonding environments. That paper lists 1.2 A (for O) and 0.44 A (for H) as the appropriate radii for water (which add up to substantially more than the OH bond length). Meanwhile, a table of covalent radii (e.g. from the Cambridge Structural Database, http://en.wikipedia.org/wiki/Covalent_radius) lists 0.66 A for O and 0.31 A for H, almost a factor of 2 different for O. Moreover, the radii ought to vary for different bonding environments. Although JCP 116, 6941 (2002) indicates that there is some tolerance for a range of values, it is certainly possible to mess up the calculation by using bad values.</div><div><br></div><div>The best I can think to do now is use 1.2 A and 0.44 A for water, and use values from the Cambridge Structural Database for all the protein atoms. Is that reasonable? I feel this is probably not that big a deal, but I'm new to QM/MM and am trying to ensure I don't make an error that will be difficult to control for.</div><div><br></div><div>Thanks,</div><div><br></div><div>Josh</div><div><br></div><div>On Wednesday, March 25, 2015 at 10:48:44 AM UTC-5, Teo wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div style="word-wrap:break-word"><div>MM radii are just covalent radii (it’s very well explained in few papers).</div><div>Find a table with covalent radii and you’ll have solved your problem (that’s what the mentioned script was doing)...</div><div><br></div><div>Teo</div><div><br></div><div><blockquote type="cite"><div>On 25 Mar 2015, at 15:29, Joshua Carr <<a href="javascript:" target="_blank" gdf-obfuscated-mailto="s3IVn9qJM9kJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">jken...@gmail.com</a>> wrote:</div><br><div><div dir="ltr">I am having the same problem! I need to define radii for some new atom types, and I have no idea how to do so. Someone help, please!<div><br></div><div>Josh Carr<br><br>On Thursday, March 12, 2015 at 8:34:33 AM UTC-5, Kristof wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi all, <br>I am trying to set up a QMMM simulation, for which I will need the covalent radii for the electrostatic coupling.<br>I found some old mails in this list that mentioned a "<a href="http://covalent_radii.pl/" rel="nofollow" target="_blank" onmousedown="this.href='http://www.google.com/url?q\75http%3A%2F%2Fcovalent_radii.pl%2F\46sa\75D\46sntz\0751\46usg\75AFQjCNHZl2l-dTTwkXDzWLxFusy8ySeHag';return true;" onclick="this.href='http://www.google.com/url?q\75http%3A%2F%2Fcovalent_radii.pl%2F\46sa\75D\46sntz\0751\46usg\75AFQjCNHZl2l-dTTwkXDzWLxFusy8ySeHag';return true;">covalent_radii.pl</a>" script that is supposed to yield those. It is supposed to be in the files section of this group.<br>However, I cannot find it (I can't find the files section for that matter), can someone please help me out?<br><br>Best,<br>Kristof<br></div></blockquote></div></div><div><br></div>

-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="s3IVn9qJM9kJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp2k+uns...@googlegroups.<wbr>com</a>.<br>
To post to this group, send email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="s3IVn9qJM9kJ" rel="nofollow" onmousedown="this.href='javascript:';return true;" onclick="this.href='javascript:';return true;">cp...@googlegroups.com</a>.<br>
Visit this group at <a href="http://groups.google.com/group/cp2k" target="_blank" rel="nofollow" onmousedown="this.href='http://groups.google.com/group/cp2k';return true;" onclick="this.href='http://groups.google.com/group/cp2k';return true;">http://groups.google.com/<wbr>group/cp2k</a>.<br>
For more options, visit <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href='https://groups.google.com/d/optout';return true;" onclick="this.href='https://groups.google.com/d/optout';return true;">https://groups.google.com/d/<wbr>optout</a>.<br>
</div></blockquote></div><br></div></blockquote></div></div>