<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns="http://www.w3.org/TR/REC-html40">

<head>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=us-ascii">
<meta name=Generator content="Microsoft Word 11 (filtered medium)">
<style>
<!--
 /* Style Definitions */
 p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0in;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman";}
a:link, span.MsoHyperlink
        {color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {color:purple;
        text-decoration:underline;}
span.EmailStyle17
        {mso-style-type:personal-compose;
        font-family:Arial;
        color:windowtext;}
@page Section1
        {size:595.3pt 841.9pt;
        margin:70.85pt 70.85pt 70.85pt 70.85pt;}
div.Section1
        {page:Section1;}
-->
</style>

</head>

<body lang=EN-US link=blue vlink=purple>

<div class=Section1>

<p class=MsoNormal><font size=2 face=Arial><span style='font-size:10.0pt;
font-family:Arial'>Dear GROMACS users<br>
I&nbsp; want to run simulations of polyethyleneoxide melt. In literature I have
found united atoms force fields that performed well for this kind of
simulation. For dihedral potential I need to extend the Ryckaert-Bellemans
potential up to the order of six. I have looked at the mailing list and found
some information which give advice to change the code for this kind of
potential. But I could not figure out where to do this changes and how to do
them. I would need some detail to do it. I do not have experience with
something like this. <br>
<br>
My second question also concerns PEO. I tried to use GROMOS force field as well
for simulations. I used the parameters from ffG53a5 files. In ffG53a5bon.itp I
am not sure which code of dihedrals sould be used for C-C-O-C and O-C-C-O and
which one for the bond angles. In the case of bond angle I found type ga_8 as
suitable, but it seems me from the description that it sould be the same for
all three bonds in PEO. This seems to be not OK, since the bending of C-O-C
angle is characterized by larger bending constant. <br>
<br>
I would appreciate any help. Thank you in advance<br>
<br>
Zuzana.<o:p></o:p></span></font></p>

<p class=MsoNormal><font size=2 face=Arial><span style='font-size:10.0pt;
font-family:Arial'><o:p>&nbsp;</o:p></span></font></p>

</div>

</body>

</html>