24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2834  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0817化学工程与技术求调剂,一志愿中海洋319 +14 lv945 2026-04-04 14/700 2026-04-06 10:20 by 蓝云思雨
[考研] 一志愿南京航空航天大学 材料与化工329分求调剂 +8 Mr. Z 2026-04-05 8/400 2026-04-06 09:24 by dongzh2009
[考研] 294求调剂 +4 Grey_Ey 2026-04-01 5/250 2026-04-05 23:05 by Grey_Ey
[考研] 材料工程302分求调剂 +9 zyx上岸! 2026-04-04 9/450 2026-04-05 22:08 by 醉翁wl
[考研] 275求调剂 +16 waltzh 2026-04-01 16/800 2026-04-05 17:14 by Hdyxbekcb
[考研] 数一英一274机械调剂 +5 星陨流霞 2026-04-04 6/300 2026-04-05 11:38 by arrow8852
[考研] 求调剂 +3 小沢 2026-04-03 3/150 2026-04-05 09:10 by sihailian3
[考研] 081700化学工程与技术 一志愿中海洋 323 求调剂学校 +16 披星河 2026-04-03 16/800 2026-04-05 09:00 by dick_runner
[考研] 专硕310求调剂 +5 捞捞我…. 2026-04-04 6/300 2026-04-04 23:33 by barlinike
[考研] 290求调剂 +7 luoziheng 2026-04-04 7/350 2026-04-04 23:17 by lqwchd
[考研] 材料调剂 +11 一样YWY 2026-04-02 13/650 2026-04-04 23:10 by 无际的草原
[考研] 材料调剂 +12 一样YWY 2026-04-02 13/650 2026-04-04 20:49 by 蓝云思雨
[考研] 调剂 +4 是可乐不是可乐 2026-04-04 4/200 2026-04-04 19:41 by 唐沐儿
[考研] 317分 一志愿江南大学 化学工程学硕 求调剂 +6 YinTai 2026-04-03 6/300 2026-04-03 22:30 by 无际的草原
[考研] 材料科学与工程339求调剂 +12 hyz0119 2026-03-31 13/650 2026-04-03 18:33 by ls刘帅
[考研] 农学考研求调剂 +3 dkdkxm 2026-04-01 3/150 2026-04-02 16:04 by wangjagri
[考研] 江苏科技大学招材料研究生 +4 Su032713. 2026-04-01 5/250 2026-04-01 22:03 by cccchenso
[考研] 省双一流重点一本大学招收调剂 +4 wwwwffffff 2026-03-31 7/350 2026-04-01 15:23 by wwwwffffff
[考研] 一志愿北交材料工程总分358 +5 cs0106 2026-04-01 7/350 2026-04-01 11:45 by wangjy2002
[考研] 083000环境科学与工程调剂,总分281 +4 橙子(胜意) 2026-03-30 4/200 2026-03-31 00:44 by Linzejun
信息提示
请填处理意见