24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3508  |  回复: 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的回帖

xujc1983

木虫 (著名写手)

苦逼青椒一枚

引用回帖:
8楼: Originally posted by 团小暖 at 2016-07-12 17:12:16
那请问您用这些参数做出来了吗?最初的EPM2模型中C-C/KB的单位是K,这个是什么意思呢?...

应该就是对照EPM2力场的那篇文献做的
KB不是玻尔兹曼常数嘛,所以能量可以借此用温度K表示
9楼2016-07-13 09:14:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

jiezhou5322

银虫 (正式写手)

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

匿名

用户注销 (小有名气)

本帖仅楼主可见
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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 280求调剂 +12 李rien 2026-04-04 12/600 2026-04-10 21:19 by zhouxiaoyu
[考研] 一志愿京区985,085401,与本科专业一致,电子信息工程, +4 阳光开朗的男孩 2026-04-10 4/200 2026-04-10 18:27 by shenrf
[考研] 080500求调剂 +14 黄宇博 2026-04-06 14/700 2026-04-10 13:09 by vgtyfty
[考研] 材料调剂 +5 hzhahg 2026-04-06 5/250 2026-04-10 10:10 by may_新宇
[考研] 292求调剂 +9 笑笑袁 2026-04-09 9/450 2026-04-10 10:05 by LHGeng
[考研] 297求调剂 +8 Kwgyz 2026-04-09 8/400 2026-04-09 23:22 by may_新宇
[考研] 0703调剂 +18 拾玖壹 2026-04-04 20/1000 2026-04-09 20:53 by zhouxiaoyu
[考研] 计算机11408,286分求调剂 +9 木子念晞 2026-04-05 9/450 2026-04-09 20:04 by vgtyfty
[考研] 一志愿085502,267分求调剂 +10 再忙也要吃饭啊 2026-04-08 11/550 2026-04-09 19:51 by gong120082
[考研] 1U盾记得记得就 +9 sanjin020722 2026-04-08 10/500 2026-04-09 14:11 by 诗与自由
[考研] 一志愿北京2,材料与化工308求调剂 +17 熊二想上岸 2026-04-04 18/900 2026-04-09 09:13 by leyan1127
[考研] 298求调剂 +4 manman511 2026-04-05 4/200 2026-04-08 16:50 by tjzhao
[考研] 263分B区求调剂 +6 李nihao 2026-04-08 6/300 2026-04-08 09:38 by 南开小綦
[考研] 机械工程264学硕求调剂 +3 qiushangxian 2026-04-06 3/150 2026-04-08 01:53 by Linzejun
[考研] 323求调剂 +3 林zlu 2026-04-07 4/200 2026-04-07 23:21 by lbsjt
[考研] 071000生物学,一志愿深圳大学296分,求调剂 +12 TIckLw 2026-04-06 13/650 2026-04-07 20:34 by lijunpoly
[考研] 22408 318分求调剂 +4 勤奋的小笼包 2026-04-06 6/300 2026-04-07 15:05 by 纸鹤555
[考研] 348求调剂 +3 车厘子zzz 2026-04-05 3/150 2026-04-05 20:30 by 啵啵啵0119
[考研] 085602调剂 初试总分335 +12 19123253302 2026-04-04 12/600 2026-04-05 08:08 by 544594351
[考研] 求生物学专业调剂-332分 +5 云朵遛弯指南 2026-04-04 5/250 2026-04-04 10:05 by rzh123456
信息提示
请填处理意见