<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:'times new roman', 'new york', times, serif;font-size:12pt"><div>Hi,</div><div>&nbsp;&nbsp;I was interested in calculating the diffusion constant of the center of mass of entire lipid-bilayer ( not individual lipid molecules). &nbsp; Regarding this, I had two doubts I wanted to clarify:</div><div>&nbsp;&nbsp; 1. Since I am interested in calculating the diffusion constant of the bilayer iteself, I guess I should&nbsp;&nbsp;allow the drift of the bilayer and hence &nbsp;not then remove the center of mass motion of bilayer (and solvent ) separately. Is that right ? &nbsp;If so, is it still OK to remove the center of mass of whole system ( bilayer + water together ) ?&nbsp;</div><div><br></div><div>2. &nbsp; How to calculate the mean square displacement of &nbsp;center of mass &nbsp;of the entire bilayer. I guess, g_msd program &nbsp;by default,
 calculate &nbsp;the diffusion constant of the individual atoms or molecules. &nbsp; &nbsp;But, is there a way to get the diffusion constant of the entire bilayer center of mass ? &nbsp;If I specify the index group which consist of all the atoms of the bilayer, will g_msd provide the diffusion constant of its center of mass or it will give the diffusion constant of individual atoms ? &nbsp;I guess, in that case, one should not use -rmcomm option . Is that right ?</div><div><br></div><div><br></div><div><br></div><div>Thanks</div><div>Jagannath</div><div style="position:fixed"></div>


</div><br>







      </body></html>