24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2836  |  回复: 5

sdwawa

铜虫 (初入文坛)

[求助] LAMMPS中描述C-Si-O势能 已有1人参与

如何在LAMMPS中描述C-Si-O三元相互作用的势能。我问过一篇文章的作者,说直接用SiC和SiO的tersoff势,C-O项用L-J势描述固态CO2。
我自己写了个可以运行,但是不知道这么写对不对,总觉得有问题。
mass         1 12.0   #type 1 = C
mass         2 28.4   #type 2 = Si
mass         3 16.0   #type 3 = O

pair_style   lj/cut 2.5
pair_coeff   *  *  1 1
pair_coeff   1  3  1 1.1 2.8

pair_style  tersoff/zbl  
pair_coeff   *  *  SiC.tersoff.zbl C Si C

pair_style   tersoff
pair_coeff   *  *  SiO.tersoff Si O Si
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsloneil

专家顾问 (正式写手)

建议使用pair_style hybrid/overlay
2楼2014-09-11 10:35:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sdwawa

铜虫 (初入文坛)

引用回帖:
2楼: Originally posted by lsloneil at 2014-09-11 10:35:31
建议使用pair_style hybrid/overlay

谢谢 我之前试过用hybrid但是报错了
pair_style   hybrid/overlay tersoff/zbl tersoff lj/cut
pair_coeff   1*2  1*2  SiC.tersoff.zbl  C Si C
pair_coeff   2*3  2*3  SiO.tersoff  Si O Si
pair_coeff   1  3  1 1.1 2.8

***************************************
pair_style   hybrid/overlay  tersoff/zbl tersoff lj/cut
ERROR: Illegal pair_style command (../pair_lj_cut.cpp:430)
3楼2014-09-11 14:16:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

oxox6085

专家顾问 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
pair_style   hybrid/overlay tersoff/zbl tersoff lj/cut 7.0
pair_coeff   *  *  tersoff/zbl SiC.tersoff.zbl  C Si NULL
pair_coeff   *  *  tersoff SiO.tersoff  NULL Si O
pair_coeff   1  3  lj/cut 1.1 2.8

再试试
好好看买妞,自己进行调试!
有些事情可能只有计算机模拟能告诉我们答案
4楼2014-09-11 17:53:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

oxox6085

专家顾问 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
sdwawa: 金币+5, ★★★★★最佳答案, 非常感谢。还是命令看的不够仔细。谢谢了 2014-09-12 10:27:49
引用回帖:
4楼: Originally posted by oxox6085 at 2014-09-11 10:53:18
pair_style   hybrid/overlay tersoff/zbl tersoff lj/cut 7.0
pair_coeff   *  *  tersoff/zbl SiC.tersoff.zbl  C Si NULL
pair_coeff   *  *  tersoff SiO.tersoff  NULL Si O
pair_coeff   1  3  lj/cut 1.1 ...

lj/cut 我写错了,应该是
pair_style   hybrid/overlay tersoff/zbl tersoff lj/cut 2.8
...
pair_coeff   1  3  lj/cut 1.0 1.1
有些事情可能只有计算机模拟能告诉我们答案
5楼2014-09-11 17:55:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

巴丁

铜虫 (小有名气)

引用回帖:
5楼: Originally posted by oxox6085 at 2014-09-11 17:55:46
lj/cut 我写错了,应该是
pair_style   hybrid/overlay tersoff/zbl tersoff lj/cut 2.8
...
pair_coeff   1  3  lj/cut 1.0 1.1...

请问lj/cut后面epsilon和sigma是如何确定的?
6楼2017-02-21 20:28:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sdwawa 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿郑州大学材料与化工085600,求调剂 +26 吃的不少 2026-04-02 26/1300 2026-04-06 13:40 by BruceLiu320
[考研] 297分083200求助 +8 aekx 2026-04-05 8/400 2026-04-06 10:35 by imissbao
[考研] 301求调剂 +3 XYPLR 2026-04-05 4/200 2026-04-05 19:07 by XYPLR
[考研] 282求调剂 +7 aaa车辆 2026-04-02 11/550 2026-04-05 17:24 by yulian1987
[考研] 275求调剂 +16 waltzh 2026-04-01 16/800 2026-04-05 17:14 by Hdyxbekcb
[考研] 一志愿北京化工大学,初试成绩350求调剂 +9 沿岸?贝壳 2026-04-04 14/700 2026-04-05 01:09 by 沿岸?贝壳
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 14:46 by 湘农储能材料
[考研] 22408求调剂 354分 可跨专业 +3 hannnnnnn 2026-04-04 3/150 2026-04-04 14:35 by 土木硕士招生
[考研] 本科985,专业0812分336求调剂 +4 莫莫很行 2026-04-03 4/200 2026-04-03 21:31 by zhq0425
[考研] 322求调剂 +4 FZAC123 2026-04-03 4/200 2026-04-03 20:55 by zhq0425
[考研] 289-求调剂 +4 这里是_ 2026-04-03 4/200 2026-04-03 14:23 by 1753564080
[考研] 326分求调剂 +3 于是乎呢 2026-04-01 5/250 2026-04-03 14:23 by 于是乎呢
[考研] 生物学硕341求调剂 +4 你笑起来像云朵 2026-04-03 4/200 2026-04-03 10:32 by macy2011
[考研] 材料调剂 +4 一样YWY 2026-04-03 4/200 2026-04-03 09:48 by 蓝云思雨
[考研] 一志愿a区211,085601-307分求调剂 +13 党嘉豪 2026-03-31 26/1300 2026-04-03 08:33 by 495374996
[考研] 279求调剂 +5 傅文秋 2026-04-02 5/250 2026-04-02 18:10 by 笔落锦州
[考研] 一志愿北京科技大学材料学硕328分求调剂 +6 1段时间 2026-03-31 7/350 2026-04-02 13:57 by 3041
[考研] 化学工程专硕324分,一志愿中国矿业大学求调剂 +7 耿耿1314 2026-04-01 7/350 2026-04-02 07:40 by 尚水阁主
[考研] 一志愿北交材料工程总分358 +5 cs0106 2026-04-01 7/350 2026-04-01 11:45 by wangjy2002
[考研] 一志愿食品科学与工程083200求调剂 +4 XQTJZ 2026-03-30 4/200 2026-03-31 04:10 by fmesaito
信息提示
请填处理意见