[gmx-developers] Fluctuating charge code, Q-Chem interface, force matching.

Gerrit Groenhof ggroenh at gwdg.de
Mon Nov 5 08:05:00 CET 2012

Dear Lee-Ping,

On 1 and 2 we already had some exchange in the past, and I am pleased to see that you got it done. 

Concerning the interface with Q-chem, if it is documented, and works without too much coding on the Q-chem part (in contrast to gaussian!!), I'd be happy to have it. There is an initiative to restructure the complete qmmm setup, as we now are getting too many program-specific interfaces.

On 1, I think it will be interesting to have, but I think it also depends on how well it is integrated (and integratable). I suppose it would also have to work on GPU now?

Concerning the forcematching. Is this the same procedure as that of Voth and co-workers. We used that, but as a seperate tool. Unless it is something else, I'd think such code would be more suitable as a gmxtool?


On Nov 4, 2012, at 9:43 PM, Lee-Ping Wang wrote:

> Dear developers,
> I'm a postdoctoral fellow in the Stanford Chemistry department working with Vijay Pande and Todd Martinez.  I earned my Ph.D. from the MIT Chemistry department under the supervision of Troy Van Voorhis.  I've been using MD software as a central part of my research for over five years, and I am definitely a big fan of GROMACS. :)
> Over the course of my Ph.D. I did a lot of method development in Gromacs.  I'm writing to find out whether there's interest in integrating any of my contributed code into the Gromacs code base.  
> My contributions are in three main categories:
> 1) A fluctuating charge code following the work of Steve Rick and Bruce Berne, with further improvements following the work of Jiahao Chen and Todd Martinez.
> 2) A QM/MM interface with Q-Chem, mainly following the existing QM/MM interface with Gaussian.
> 3) Energy and force matching code for force field optimization.  This consists mainly of adding analytic first and second parametric derivatives for the energy and force terms into the low-level subroutines that compute the force field contributions.
> Items (1) and (3) are incorporated into an upcoming publication which was just accepted to JCTC.  
> These methods have been implemented into version 4.0.7.  Most of the code (perhaps 95%) is in separate files with as little modification to the original code as possible.  All of the functionality is controlled through the input files and command line arguments.  There is no performance impact for normal MD runs because the parametric derivatives are only switched on when requested.
> If there is interest on your side, I'd be happy to do all of the work.  I think the work involved is mainly (1) updating the code to be compatible with the latest development branch, (2) writing unit tests, and (3) creating documentation.
> Please let me know and I'll do my best to follow the directions on the Gromacs website and stay in touch.
> Thanks,
> - Lee-Ping--
> gmx-developers mailing list
> gmx-developers at gromacs.org
> http://lists.gromacs.org/mailman/listinfo/gmx-developers
> Please don't post (un)subscribe requests to the list. Use the
> www interface or send it to gmx-developers-request at gromacs.org.

More information about the gromacs.org_gmx-developers mailing list