| 查看: 6273 | 回复: 29 | ||||
[求助]
LAMMPS中金属氧化物的势能 已有1人参与
|
||||
| 楼主一直用LAMMPS里EAM形式的势能，最近遇到一UO2的势能构成为：Buckingham+Morse+Coloumb type potential， 相关文献：Journal of Alloys and Compounds 360 (2003) 210–216。从势能形式上看，貌似是阴阳离子相互作用，楼主不清楚这应该用LAMMPS中的那种形式的势能及对应的命令，求大神指导。 |
回复此楼» 猜你喜欢
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有7人回复
-
今天基金会出结果吗？20260819
已经有18人回复
-
我面上完蛋了
已经有14人回复
-
范进中举一文的中心思想
已经有8人回复
-
只有每年这种时候来逛逛小木虫
已经有27人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
-
没有任何消息-是不是就凉了
已经有8人回复
-
人气不行了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- LAMMPS中描述C-Si-O势能
已经有5人回复
- lammps中界面的建立
已经有9人回复
- 如何通过分子模拟计算材料的渗透率？
已经有11人回复
- lammps势文件求助
已经有7人回复
- 请问gromacs导入EAM势能函数的问题和lammps导入martini粗粒化力场的问题
已经有4人回复
- 一个薄膜体系中掺杂一个磁性原子，不考虑磁性原子之间相互作用时如何优化结构？？
已经有13人回复
- 请教大家一个关于压力、温度模拟的问题
已经有4人回复
- 关于势函数的求助
已经有23人回复
- 关于化学式写法的问题，Ca4LaO(BO3)3，这里面写在括号外面的O是什么意思？
已经有10人回复
- 寻求lammps的应用 欢迎大家畅所欲言
已经有8人回复
- 【求助】gulp参数如何获得
已经有5人回复
24楼2015-01-15 11:23:49
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
8楼2015-01-08 11:06:35
9楼2015-01-08 11:39:14
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
gulubaozi10楼2015-01-08 12:06:48
|
楼主，经过两天的郁闷弄清了LAMMPS手册中coulomb表达式及其设置使用。 coulomb力的计算公式可以参看wiki ，通过力与势能的关系，可得到coulomb势，与LAMMPS手册中所用公式对比，可得手册中的C （energy-conversion constant） 实际上是wiki 中coulomb力的那个常数项。 虫子cgzhang_gg在他的帖子里计算了LAMMPS中的这个C值（帖子http://muchong.com/html/201306/6016306.html），楼主通过计算coulomb力常数项表达式的值，并与上值对比，确信LAMMPS手册中所用coulomb势公式中的C （energy-conversion constant）就是coulomb力常数项，而LAMMPS公式中的epsilon（其值通过dielectric command设置）仅仅是个比例常数，不具备任何物理意义。 |
12楼2015-01-08 16:56:49
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
18楼2015-01-10 11:55:48
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
|
不需要ewald/disp，直接用ewald就行。 你这里的pair_style是 buck/coul/long，意思是仅对库仑相互作用采取长程处理。但如果你的pair_style是buck/long/coul/long，也就是说对Buckingham也要作长程处理，根据kspace_style部分的说明 http://lammps.sandia.gov/doc/kspace_style.html 这时就要用ewald/disp了。 |
» 本帖已获得的红花（最新10朵）
26楼2015-05-07 14:04:40
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
|
yes |
5楼2015-01-05 23:17:28
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
6楼2015-01-05 23:18:55
7楼2015-01-08 10:20:48
