½ð³æ (ÖøÃûдÊÖ)
ÁõËÕÖÝ ![]()
|
¡ï ¡ï ¡ï ¡ï Сľ³æ(½ð±Ò+0.5):¸ø¸öºì°ü£¬Ð»Ð»»ØÌû½»Á÷ lei0736(½ð±Ò+3):лл 2010-04-14 15:56
Á¦³¡ÖеĵçºÉ²¢·ÇÔ×ÓÕæÊµµÄµçºÉ, ¶øÊǸù¾Ýijһ¼ÆËã·½·¨µÃ³öµÄÊýÖµ
·Åµ½Õû¸öÁ¦³¡²ÎÊýÖнøÐÐÄâºÏµÃµ½µÄ. ²»Í¬Á¦³¡¼ÆËãµçºÉµÄ·½·¨²»Í¬, ÊýÖµ²î±ðÒ²ºÜ´ó
±ÈÈçcharmm ʹÓÃhf/6-31g* »ù×éϵÄMullikin µçºÉ, ¶øamber ÓõÄÔòÊÇESP.
Methodologies for Developing the Partial Charges
Now is a good time to discuss the determination of the partial charges appearing
in the topology file. Two widely used force fields are CHARMM and AMBER
¨C and they both use di?erent protocols to develop the partial atomic charges.
CHARMM starts with the Mullikin charges and then strategically places a
water molecule near any polar group. CHARMM¡¯s charges are refined so that
they best reproduce the interaction energy between the water and the polar
group. To recreate these values, the partial charges of the polar group are then
altered until this interaction energy is accurately replicated with the MD sim-
ulations. Several concrete examples of this charge-fitting strategy are worked
out in Ref. [7].
AMBER, on the other hand, begins with the ESP charges and then uses a
program they distribute called RESP (for Restricted ESP). This charge-fitting
scheme replicates the electrostatic potential at a distance out from the nucleus
(much like the surface you built in the previous exercise) with only the atom
centers having charge. Of course, the electrostatic surface generated using this
methodisoftendegenerate; meaningtherecouldbemorethanonesetofcharges
that would reproduce the potential surface. To resolve this issue, it is best to
further restrict your potential surface by supplying it an initial set of charges for
the program¡¯s iterations, or further restrictions derived from chemical insight
However they are generated, the importance of a good set of partial charges
cannot be overstated. Much of the energy (and indirectly then, the behavior)
of the system will be calculated using these values. |
|