|
|
[ÇóÖú]
¹òÇó¸÷λ´óÉñ£¬°ïæ¿´¿´Ð¡µÜµÄ¶àÌåÊÆÃüÁî±àдÊÇ·ñºÏÊÊ£¬ÇóÖ¸µã£¡£¡
¹òÇó¸÷λ´óÉñ£¬ËÔÚ¼ÆËãÖÐÓõ½Á˶àÌåÊÆµÄÃüÁСµÜ³õ´ÎʹÓò»ÖªµÀ±àдµÄÓÐûÓÐÎÊÌ⣬ϣÍû¸÷λ´óÉñ²»Áߴͽ̣¬¶à¶àÖ¸µ¼£º
Ñо¿µÄÌåϵ°üÀ¨¾ÛÒÒÏ©ºÍʯīϩ£¬ÔÚdataÎļþÖгöÏÖÁËÈýÖÖÔ×ÓÀàÐÍ£¬Ô×ÓÀàÐÍ1ΪH,Ô×ÓÀàÐÍ2ΪC(¾ÛÒÒÏ©ÖеÄ̼Ô×Ó)£¬µÚÈýÖÐÀàÐÍΪC(ʯīϩÖеÄ̼Ô×Ó)¡£
¾ßÌåµÄÒªÇóÈçÏ£º1.COMPASSÁ¦³¡Ó¦ÓÃÓÚ¾ÛÒÒÏ©ÖУ»
2. AIREBOÓ¦ÓÃÓÚʯīϩÖУ»
3.¾ÛÒÒÏ©Óëʯīϩ֮¼ä²ÉÓÃL-JÊÆ¡¾Lorentz¨CBerthelot mixing rules (¦Åij = sqrt£¨¦Åi¦Åj£©, wherej ); ¦Òij =(¦Òi +¦Òj )/2,where ¦Å and ¦Ò are the energy and distance constants, respectively. The subscripts i,j refer to different atom species¡¿¡£
4.²ÎÊýÉèÖÃÒªÇóÈ總¼þͼÖÐËùʾ
ÒÔÏÂΪСµÜ±àдµÄ¶àÌåÊÆÔËÐÐÓï¾ä£º
pair_style hybrid lj/class2/coul/cut 10.0 9.5 airebo 3.0 1 1
pair_coeff 1 3 lj/cut/coul/cut 0.001586 3.2035¡¾ sqrt(0.000867*0.002390)=0.001586 ,(2.995+3.412)/2=3.2035 ¡¿
pair_coeff 2 3 lj/cut/coul/cut 0.002366 3.711 ¡¾ sqrt(0.002342*0.002390)=0.002366 ,(4.010+3.412)/2=3.711 ¡¿
pair_coeff 3 3 ../potentials/CH.airebo C
pair_coeff 1 1 lj/class2/coul/cut 0.000867 2.955
pair_coeff 2 2 lj/class2/coul/cut 0.002342 4.010
![¹òÇó¸÷λ´óÉñ£¬°ïæ¿´¿´Ð¡µÜµÄ¶àÌåÊÆÃüÁî±àдÊÇ·ñºÏÊÊ£¬ÇóÖ¸µã£¡£¡]()
1.png
![¹òÇó¸÷λ´óÉñ£¬°ïæ¿´¿´Ð¡µÜµÄ¶àÌåÊÆÃüÁî±àдÊÇ·ñºÏÊÊ£¬ÇóÖ¸µã£¡£¡-1]()
1.png |
|