<DIV>
<DIV>&nbsp;Hi, everyone.</DIV>
<DIV>
<DIV class=FoxDiv20120516095839375739>
<DIV>
<DIV class=FoxDiv20120516095726296235>
<DIV>&nbsp;</DIV>
<DIV>It is my first use of <A href="http://lammps.sandia.gov/">Gromacs</A> and I am looking for a numerical scheme for one specific constrained SDE, the constrain <SPAN class=subtitle><SPAN dir=ltr>is a macroscopic one, i.e., overdamped Langevin(Brownian dynamics) equations with an equality constraint which is </SPAN></SPAN><SPAN class=subtitle><SPAN dir=ltr>expressed in form of expectation(or moment of n-th order). </SPAN></SPAN><SPAN class=subtitle><SPAN dir=ltr>This is unlike the common 'micro' constraint(simply </SPAN></SPAN><SPAN class=subtitle><SPAN dir=ltr>function in terms of variable in SDE) . I&nbsp;can not find such a subrountine </SPAN></SPAN><SPAN class=subtitle><SPAN dir=ltr>for solving constrained SDE in </SPAN></SPAN><SPAN class=subtitle><SPAN dir=ltr>MD codes in <A href="http://lammps.sandia.gov/">Gromacs</A>.&nbsp; </SPAN></SPAN><SPAN class=subtitle><SPAN dir=ltr>I see there is an algorithm called "SHAKE" or "RATTLE "which and it seems&nbsp;they could </SPAN></SPAN><SPAN class=subtitle><SPAN dir=ltr>implement Langevin(or Brownian) dynamics with constraints. If convenient, could anyone help to point out which subroutine </SPAN></SPAN><SPAN class=subtitle><SPAN dir=ltr>is specific for implementing SHAKE/RATTLE? Thanks in advance.</SPAN></SPAN></DIV>
<DIV>
<DIV><SPAN class=subtitle><SPAN dir=ltr></SPAN></SPAN>&nbsp;</DIV>
<DIV><SPAN class=subtitle><SPAN dir=ltr></SPAN></SPAN>&nbsp;</DIV>
<DIV><SPAN class=subtitle><SPAN dir=ltr></SPAN></SPAN>&nbsp;</DIV>
<DIV>
<HR style="WIDTH: 210px; HEIGHT: 1px" align=left color=#b5c4df SIZE=1>
</DIV></DIV>
<DIV><SPAN>
<DIV><SPAN style="FONT-SIZE: 10.5pt; COLOR: #000000; FONT-FAMILY: ËÎÌå">&nbsp; 
<HR style="WIDTH: 100px" align=left color=#b5c4df SIZE=1>

<DIV></DIV>
<DIV></DIV>
<DIV></DIV>
<DIV></DIV>
<DIV>Kevin Len</DIV>
<DIV>Department of Materials Science</DIV>
<DIV>Fudan University</DIV>
<DIV>220 Handan Road Shanghai, CHINA, 200433</DIV></SPAN></DIV></SPAN></DIV></DIV></DIV></DIV></DIV></DIV>