<p dir="ltr">Hi Alex,</p>
<p dir="ltr">I'm one of the authors of the article mentioned by Gerrit. As you also observe, there is no easy way of combining graphene and biological compounds, each community has its own tools which only covers a part of the possible applications. F.e. we were interested in mechanical properties, so our FF tries to reproduce existing experimental data, but ignores electrostatics... although electrostatics are of great importance in biological molecular systems which often carry partial charges. To my knowledge (which is likely already obsolete due to the high pace of research related to graphene and CNTs), there isn't a FF which comprehensively reproduces all kinds of physical and chemical properties of these compounds. We mention specifically in the article that our FF makes compromises in order to keep a simple functional form and be fast, which is a must for any FF that hopes to be used for simulations on biological time scales. In particular when used in GROMACS which is fast by default :-)</p>
<p dir="ltr">I beg to disagree with your statement on the Morse potential: it is present in GROMACS for a long time, although no FF for biological compounds use it at the moment. We also chose it because it is a pairwise potential, making it simple to integrate in other (bio oriented) simulation packages like CHARMM.</p>
<p dir="ltr">In my opinion, the integration of a graphene/CNT FF in GROMACS is a multidimensional problem, having to consider the speed, the ease of implementation, the compatibility to existing (bio) FFs, as well as the ability of reproducing experimental data. Sadly the great majority of the (many) graphene/CNT FFs ignores the first two or three aspects, making them a poor choice for this task... How do the proposed FFs fare in these directions?</p>
<p dir="ltr">Cheers,<br>
Bogdan</p>
<div class="gmail_quote">On Nov 19, 2014 9:23 AM, "Alex" <<a href="mailto:nedomacho@gmail.com">nedomacho@gmail.com</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi there,<div><br></div><div>Sorry, I do not have institutional access to this journal. </div><div>As a general comment, I am not at all doubting that there is work on the subject of developing Morse-like, or any other popular potential forms. For instance, the standard bond-order potential is also of the Morse type with bond coordination terms. Of course,the work you mentioned could be used as basis (additional to what I cited) for parameterizing the FFs already implemented within GROMACS). </div><div><br></div><div>As far as I know, no bonding (intramolecular) interactions in GROMACS are based on Morse-like expressions. As of right now, people who want to do large-scale simulations of graphene or BN just go to LAMMPS, and those who want to do biomolecular simulations involving those membranes within GROMACS, are sort of out of luck. I am just generally hoping there's some interest on your end in actually having something built-in. I am not really here to promote my paper, I just like GROMACS. :)</div><div><br></div><div>Also, could someone rename the thread subject? I really meant intramolecular (bonding) interactions.</div><div><br></div><div>Thanks,</div><div><br></div><div>Sasha</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Nov 19, 2014 at 1:07 AM, Gerrit Groenhof <span dir="ltr"><<a href="mailto:ggroenh@gwdg.de" target="_blank">ggroenh@gwdg.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Hi, <div><br></div><div>THere has been some work done on FF parameterization recently by the Grater group, see: <div><br></div><div><a href="http://pubs.rsc.org/en/content/articlelanding/2014/cp/c3cp55340j#!divAbstract" target="_blank">http://pubs.rsc.org/en/content/articlelanding/2014/cp/c3cp55340j#!divAbstract</a><div><div><div><br></div><div>Gerrit</div><div><br></div><div><br></div><div>On Nov 19, 2014, at 9:41 AM, Alex wrote:</div><br><blockquote type="cite"><div dir="ltr"><font face="arial, helvetica, sans-serif">Dear developers,</font><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">I really hope to get your attention.</font></div><div><font face="arial, helvetica, sans-serif">I may be wrong, but GROMACS still doesn't have a reasonable built-in parameterization for graphene as part of any FF. There are some tutorials online, but they are all quite flaky, based on the parameterizations of organic molecules (as expected). Given the few rounds of chatting with the reviewers I've had for my last paper on the subject, the parameterization is sorely desired. So, maybe it could be a good idea for you to look into it, as people do want a proper model of graphene to simulate it in water in the presence of all the other things that GROMACS simulates so well.</font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">For a DREIDING-like model, here are possibly useful papers containing pretty much complete descriptions of the parameter sets:</font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">1. J. H. Los et al<i> Scaling properties of flexible membranes from atomistic simulations: Application to graphene</i>, Phys Rev B. 80, 121405 R 2009 . </font></div><div><font face="arial, helvetica, sans-serif">(this model is parameterized by differentiating the LCBOPII bond-order potential, which fails at describing the phonon spectrum of graphene. No dihedral constraints, see below Fig. 1)</font></div><div><font face="arial, helvetica, sans-serif"><br></font></div><div><font face="arial, helvetica, sans-serif">2. <span style="color:rgb(0,0,0);font-size:12px;line-height:16.2000007629395px">A Shakouri, T Y Ng</span><span style="color:rgb(0,0,0);font-size:12px;line-height:16.2000007629395px"> and R M Lin </span></font><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:12px;line-height:16.2000007629395px"><i>A new REBO potential based atomistic structural model for graphene sheets </i>Nanotechnology 22(29), 295711, 2011.</span></font></div><div><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:12px;line-height:16.2000007629395px">(no dihedrals here either, based on the second-generation bond-order potential, which fails in the harmonic spectrum)</span></font></div><div><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:12px;line-height:16.2000007629395px"><br></span></font></div><div><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:12px;line-height:16.2000007629395px">3. A. Smolyanitsky Molecular dynamics simulation of thermal ripples in graphene with bond-order-informed harmonic constraints, Nanotechnology 25(48), 485701, 2014.</span></font></div><div><font color="#000000" face="arial, helvetica, sans-serif"><span style="font-size:12px;line-height:16.2000007629395px">(this is my work, the parameterization does have dihedrals, and it is based on the phonon-optimized second-gen bond-order potential, which is described in ref #6 of this work).</span></font></div><div><br></div><div>For the most part, all of these parameterizations do not directly correspond to the "parent" potentials (except around 0 K), as spring-like FFs will generally present a relatively stiff model, but I think this info could be a small step forward. I hope this is useful. If you have any questions on the numbers, please do let me know. </div><div><br></div><div>Thanks,</div><div><br></div><div>Sasha</div><span><font color="#888888"><div><br></div></font></span></div><span><font color="#888888">
-- <br>Gromacs Developers mailing list<br><br>* Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/GMX-developers_List" target="_blank">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" 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" 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>.</font></span></blockquote></div><br></div></div></div></div><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" target="_blank">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" 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" 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>.<br></blockquote></div><br></div>
<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" target="_blank">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" 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" 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">gmx-developers-request@gromacs.org</a>.<br></blockquote></div>