<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.7655.11">
<TITLE>CMAP Dihedrals and free energy perturbations</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/plain format -->
<BR>

<P><FONT SIZE=2>Hi<BR>
I was wondering if anyone could help me.&nbsp; I am trying to use Gromacs 4.5.4 to do some free energy calculations.&nbsp; I am using the CHARMM forcefield with CMAP and keep getting the error message:<BR>
<BR>
Fatal error:<BR>
Function type CMAP Dih. not implemented in ip_pert<BR>
<BR>
I noticed that someone had posted a while ago about a similar message for Improper Dihs. and that the solution to this was an extra line added to topsort.c.&nbsp; I think I need to add a line to topsort.c too for F_CMAP but am slightly unsure of what to do.<BR>
Any help would be appreciated.&nbsp; Is there an easy way to carry to carry out free energy perturbation calculations in Gromacs when using CMAP in the potential was well?<BR>
<BR>
Thanks in advance<BR>
Louise<BR>
<BR>
<BR>
<BR>
</FONT>
</P>

</BODY>
</HTML>