| 查看: 3728 | 回复: 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  EPM2-m.jpg |
回复此楼» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- MS里的forcite 模块只能用自带的那三个力场吗？想用trappe力场是不是要用lammps?
已经有6人回复
- LAMMPS中ReaxFF力场是怎么回事？
已经有17人回复
- 用MS建模导入LAMMPS计算，是不是不能用compass力场？
已经有9人回复
匿名
用户注销 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1014.3
- 散金: 482
- 红花: 13
- 帖子: 249
- 在线: 118.8小时
- 虫号: 0
- 注册: 2014-10-29
- 性别: MM
- 专业: 工程热物理相关交叉领域
4楼2016-07-07 15:37:37
jiezhou5322
银虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 2072.4
- 散金: 5
- 红花: 1
- 帖子: 390
- 在线: 332.6小时
- 虫号: 3093434
- 注册: 2014-03-27
- 专业: 岩土力学与岩土工程
3楼2016-03-22 16:07:04
xujc1983
木虫 (著名写手)
苦逼青椒一枚
- 应助: 71 (初中生)
- 金币: 2065.3
- 散金: 2373
- 红花: 18
- 帖子: 1096
- 在线: 575.3小时
- 虫号: 285496
- 注册: 2006-10-14
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
5楼2016-07-07 16:42:10
匿名
用户注销 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1014.3
- 散金: 482
- 红花: 13
- 帖子: 249
- 在线: 118.8小时
- 虫号: 0
- 注册: 2014-10-29
- 性别: MM
- 专业: 工程热物理相关交叉领域
6楼2016-07-08 01:11:47
20