24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3518  |  回复: 12
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

xujc1983

木虫 (著名写手)

苦逼青椒一枚

[求助] 如何在LAMMPS中利用EPM2-m力场描述超临界CO2

如题,想利用LAMMPS结合EPM2-m力场模拟超临界CO2
EPM2-m力场如附图
自己试着写的in文件如下,虽说可以运行,但C-O键长总是不对,不知哪位大侠帮忙解惑?
units            real
boundary         p p p
atom_style       full
atom_modify      map array
read_data        data.500CO2

pair_style       lj/charmm/coul/charmm 16.0 18.0 16.0 18.0
pair_coeff       1 1 0.0023266827 2.800
pair_coeff       1 2 0.0039798693 2.925
pair_coeff       2 2 0.0068077011 3.050
neighbor         2.5 bin
neigh_modify     every 100 delay 0 check no

fix                          1 all restrain bond 1 2 2000.0 2000.0 1.1491
thermo                          100
minimize         1.0e-6 1.0e-6 10000 10000
fix                          2 all npt temp 304.1 304.1 100.0 iso 72.8 72.8 2000.0
dump             1 all dcd 100 scCO2.dcd
dump_modify      1 unwrap yes
dump                         2 all xtc 100 scCO2.xtc
dump_modify      2 unwrap yes
thermo                          100
timestep         0.1    # fs
run              1000    # 50 ps

如何在LAMMPS中利用EPM2-m力场描述超临界CO2
EPM2-m.jpg
回复此楼

» 猜你喜欢

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

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

jiezhou5322

银虫 (正式写手)

楼主你的问题解决了吗,lammps下怎样编辑力场呢?
每个人只能陪你走一段路程
3楼2016-03-22 16:07:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

匿名

用户注销 (小有名气)

本帖仅楼主可见
4楼2016-07-07 15:37:37
已阅   申请模拟EPI   回复此楼   编辑   查看我的主页

xujc1983

木虫 (著名写手)

苦逼青椒一枚

引用回帖:
4楼: Originally posted by 团小暖 at 2016-07-07 15:37:37
楼主,您好。EPM2模型里已经显示键长是1.149埃米了啊,你是怎么知道键长不对的?我看到你的文件里没有键能参数,请问bond_style 选harmonic的话,键能参数是多少呢?这个力场的文章上我没找到。

bond_style       harmonic
angle_style      harmonic
pair_style       lj/charmm/coul/charmm 16.0 18.0 16.0 18.0
bond_coeff       1 2000.0 1.1491     # CO2:C-O
angle_coeff      1 147.75 174.2      # CO2:O-C-O
pair_coeff       1 1 0.055852 2.757  # CO2:C-C
pair_coeff       1 2 0.094489 2.892  # CO2:C-O
pair_coeff       2 2 0.159852 3.033  # CO2:O-O
5楼2016-07-07 16:42:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (小有名气)

本帖仅楼主可见
6楼2016-07-08 01:11:47
已阅   申请模拟EPI   回复此楼   编辑   查看我的主页
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工329分求调剂 +8 叶zilin 2026-04-10 8/400 2026-04-10 23:26 by 314126402
[论文投稿] mdpi小修rvr时间四五天了 20+3 哈哈high 2026-04-08 5/250 2026-04-10 16:02 by 北京莱茵润色
[基金申请] 有爆料,一个青年教师卖房得400万,然后换了一个四青帽子 +9 babu2015 2026-04-08 9/450 2026-04-10 11:43 by 苏东坡二世
[考研] 一志愿211 0703化学 346分求调剂 +22 土豆er? 2026-04-09 23/1150 2026-04-10 10:58 by 高维春
[考研] 调剂申请086000一志愿西北农林科技大学生物与医药320分-本科齐鲁工业大学 +3 美美女士 2026-04-09 3/150 2026-04-10 10:31 by liuhuiying09
[考研] 材料工程302分求调剂 +18 zyx上岸! 2026-04-04 18/900 2026-04-10 10:07 by 314126402
[考研] 生物与医药调剂 +5 十七sa 2026-04-05 5/250 2026-04-10 08:14 by kangsm
[考研] 化工求调剂! +34 RichLi_ 2026-04-06 34/1700 2026-04-09 20:39 by zhouxiaoyu
[考研] 材料化工总分334求调剂 +16 Riot2025 2026-04-08 17/850 2026-04-09 20:19 by maddjdld
[考研] 一志愿武理车辆 281 求调剂 +5 上岸研究生. 2026-04-07 5/250 2026-04-09 15:56 by only周
[考研] 一志愿0807 数一英一 313 有没有二轮调剂 +11 emokidd 2026-04-08 12/600 2026-04-09 09:24 by wyf236
[考研] 283求调剂 +19 A child 2026-04-04 19/950 2026-04-08 14:26 by xingguangj
[考研] 调剂 +3 电气300求调剂不 2026-04-08 6/300 2026-04-08 09:39 by 电气300求调剂不
[考研] 263分B区求调剂 +6 李nihao 2026-04-08 6/300 2026-04-08 09:38 by 南开小綦
[考研] 求调剂 +10 chenxrlkx 2026-04-05 10/500 2026-04-06 11:31 by 猪会飞
[考研] 求调剂 +7 张.1 2026-04-05 7/350 2026-04-05 20:40 by 啵啵啵0119
[考研] 调剂 +3 李广火 2026-04-05 3/150 2026-04-05 18:57 by 蓝云思雨
[考研] 考研调剂 +3 mcbbc 2026-04-04 3/150 2026-04-05 10:03 by barlinike
[考研] 0854求调剂 +4 assdll 2026-04-04 4/200 2026-04-05 09:44 by zhq0425
[考研] 316求调剂 +9 墨辰_Orion926 2026-04-04 9/450 2026-04-04 21:35 by lbsjt
信息提示
请填处理意见