| 查看: 6370 | 回复: 29 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
[求助]
LAMMPS中金属氧化物的势能 已有1人参与
|
|||||
| 楼主一直用LAMMPS里EAM形式的势能，最近遇到一UO2的势能构成为：Buckingham+Morse+Coloumb type potential， 相关文献：Journal of Alloys and Compounds 360 (2003) 210–216。从势能形式上看，貌似是阴阳离子相互作用，楼主不清楚这应该用LAMMPS中的那种形式的势能及对应的命令，求大神指导。 |
回复此楼» 猜你喜欢
现代”学阀”该如何界定
已经有12人回复
-
国社科系统bug了，是不是要放榜了？
已经有10人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有4人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有3人回复
-
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告
已经有11人回复
-
我的奶奶
已经有3人回复
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
月只蓝: 金币+10, 感谢专家指导！ 2015-01-13 14:20:28
月只蓝: 金币+10, 感谢专家指导！ 2015-01-13 14:20:28
|
很难说，具体取决于你的需要。后处理的话其实用lammps自带的pizza.py结合VMD，有时再用一下lammps的rerun和各种compute命令，已经能够处理各种文件类型转化，计算相互作用能量，部分热力学性质和rdf了。如果你的分析仅局限于这些，那么把这些掌握熟就足够了（话说把这些掌握熟也不那么容易，pizza.py和VMD我就用得不好）。 编程的好处是灵活度高，尤其是用C或fortran编程，处理数据速度比脚本语言或者matlab都快不少。像我有时要计算精度比较高的MSD（lammps自带的计算MSD功能没有充分利用所有数据，精度较低），这时编程的优势就体现出来了。另外如果你要做一些pizza.py和VMD不能做的分析，比如各种时间关联函数以及一些自由能计算，这时也需要自己编程。当然编程一开始都比较费时，几百行的程序有时都能debug几天（可能我debug能力比较弱）。 |
20楼2015-01-13 12:05:52
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
gulubaozi: 金币+5, ★★★很有帮助 2015-01-05 16:17:24
abinitio: 金币+1, 鼓励交流 2015-01-06 07:50:37
gulubaozi: 金币+5, ★★★★★最佳答案 2015-01-08 10:21:18
感谢参与，应助指数 +1
gulubaozi: 金币+5, ★★★很有帮助 2015-01-05 16:17:24
abinitio: 金币+1, 鼓励交流 2015-01-06 07:50:37
gulubaozi: 金币+5, ★★★★★最佳答案 2015-01-08 10:21:18
|
应该用pair_style hybrid/overlay命令，具体用法参照manual. 需要包含以下两种pair_style的相互作用系数。 pair_style buck/coul/long pair_style morse Good luck! |
2楼2015-01-05 11:53:42
3楼2015-01-05 16:18:02
4楼2015-01-05 17:26:27