| 查看: 5226 | 回复: 4 | |||
bs340248铜虫 (小有名气)
|
[交流]
【求助】LAMMPS中非周期性边界条件和Born-Mayer-Huggins势函数不会设置 已有1人参与
|
|
LAMMPS下 希望使用“Born-Mayer-Huggins势函数”和“非周期性边界条件”，但是 当设置 boundary f f f pair_style born/coul/long 6.872 提示ERROR: Pair style is incompatible with KSpace style 当设置 kspace_style eward 1.0e-4 时 提示ERROR: Cannot use nonperiodic boundaries with Ewald 我如何才能实现“Born-Mayer-Huggins势函数”并且使用“非周期性边界条件”？ 谢谢 |
回复此楼» 猜你喜欢
分子筛催化方向 一篇cej共一第一 可以申请上海什么学校呢
已经有4人回复
-
求助共晶分子筛的电镜图分析，有偿
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有166人回复
-
中南大学化学化工学院招收电池储能、理论计算方向博士生
已经有16人回复
-
请问哪里可以代烧带压力的炉子？
已经有6人回复
-
关于合成MIL-125的求助
已经有0人回复
-
双非本硕光催化方向想读博，请各位老师看看我！
已经有2人回复
-
瑞典林雪平大学范立洲课题组现招聘博士后、岗位制博士生
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
寻求物理化学学报最新的投稿模板！
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- word中，设置字体为多倍行距后，选中时字体总是靠上，请问如何调整为垂直方向的居中
已经有17人回复
- 关于科研经验版资源帖设置的说明
已经有4人回复
- fluent中怎样定义边界层厚度和动量层厚度
已经有5人回复
- 周期性边界条件相关问题
已经有6人回复
- BoltzTraP载流子浓度设置
已经有6人回复
- 在Lammps中安装xmovie工具时，Makefile文件怎么修改X库路径
已经有3人回复
- ansys 13.0 网格划分和边界定义怎么设置，拜求高手
已经有13人回复
- dynaform 软件模拟问题，三通管成形的冲头怎么设置
已经有5人回复
- Fluent 入口边界条件设置：是质量流量还是速度入口？
已经有11人回复
- Tecplot 怎么可以得到FLUENT计算的结果的 热边界层厚度
已经有10人回复
- 计算结果和初始边界设定不一样？
已经有8人回复
- 优化时候温度设置
已经有10人回复
- 系统重装要怎么样设置BIOS？
已经有5人回复
老虎大王
木虫 (著名写手)
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
2楼2009-09-27 20:35:34
bs340248
铜虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 2435.3
- 散金: 12
- 帖子: 108
- 在线: 61.4小时
- 虫号: 706708
- 注册: 2009-02-23
- 性别: GG
- 专业: 理论和计算化学
3楼2009-09-27 23:41:19
老虎大王
木虫 (著名写手)
- 应助: 26 (小学生)
- 贵宾: 0.17
- 金币: 4774.1
- 散金: 8
- 红花: 42
- 帖子: 1361
- 在线: 215.2小时
- 虫号: 659094
- 注册: 2008-11-21
- 专业: 金属结构材料
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+2,VIP+0): 9-30 10:34
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+2,VIP+0): 9-30 10:34
|
Kspace这个东西是为了处理库仑势。因为库仑势是1/r的函数，而且是长程势，计算1/r这样的函数很费时间，所以要用一个比较省时间的办法。Lammp目前共提供五种选择，除了none以外，就是r ewald or pppm or pppm/tip4p or ewald/n ，其中经典和基本的算法是Ewald，还有一种需要用Fourier变换的方法，叫PPPM。最后两种是它们的变体，实际上就是两类算法。但是，这两类算法都必须用到周期性边界条件（见说明书中Kspace_style的说明），所以，不适于你的体系。 那么你怎么办？我认为你只能用kspace_style none。也就是直接计算库仑势。 我的理解是这样的：kspace_style none指的是不用任何技巧，硬算库仑势1/r项。你可能理解成了kspace_style none指的是不计算库仑势。呵呵。谁理解得对呢？你不妨做个实验：你用kspace_style none来做两次计算，先把原子电荷全设成0，再把原子电荷设成真实电荷，运行个几步，比较一下计算的能量结果有多大区别，你就清楚了。 |
4楼2009-09-30 10:20:16
Camilleyuxia
铜虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 81.8
- 红花: 1
- 帖子: 21
- 在线: 8.7小时
- 虫号: 1158571
- 注册: 2010-11-29
- 性别: MM
- 专业: 流体力学
5楼2011-03-21 21:38:08
