24小时热门版块排行榜

返回列表

Materialman

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1.5
帖子: 21
在线: 31小时
虫号: 2765757
注册: 2013-10-30
专业: 高分子材料与环境

[求助] 跪求各位大神，帮忙看看小弟的多体势命令编写是否合适，求指点！！

跪求各位大神，谁在计算中用到了多体势的命令，小弟初次使用不知道编写的有没有问题，希望各位大神不吝赐教，多多指导：
   研究的体系包括聚乙烯和石墨烯，在data文件中出现了三种原子类型，原子类型1为H,原子类型2为C(聚乙烯中的碳原子)，第三中类型为C(石墨烯中的碳原子)。
具体的要求如下：1.COMPASS力场应用于聚乙烯中；
                        2. AIREBO应用于石墨烯中；
                        3.聚乙烯与石墨烯之间采用L-J势【Lorentz–Berthelot 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