<DIV>Dear users:</DIV>
<DIV>&nbsp;</DIV>
<DIV>Now I am using Gay-Berne potential to simulate the motion of rodlike molecules of nematic liquid crystals with MD.</DIV>
<DIV>But when I use it, I meet these problems,</DIV>
<DIV>1. If I use leapfrog algrithom to simulate, I am not sure how to deal with the values of half time.</DIV>
<DIV>2. In the Gay-Berne potential when the denomiator (r-delt + 1)^(-12) is near 0, the force will become infinite.&nbsp;How can I control&nbsp;this case.</DIV>
<DIV>3. I am not sure the values of the inertia I.</DIV>
<DIV>Thanks for your help</DIV>
<DIV>Shufang Fu</DIV>
<DIV>shufangfu@yahoo.com&nbsp;</DIV><p>
                <hr size=1>Yahoo! for Good<br> 
<a href="http://store.yahoo.com/redcross-donate3/">Click here to donate</a> to the Hurricane Katrina relief effort.