24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

bs340248

铜虫 (小有名气)

应助: 2 (幼儿园)
金币: 2435.3
散金: 12
帖子: 108
在线: 61.4小时
虫号: 706708
注册: 2009-02-23
性别: GG
专业: 理论和计算化学

[交流] 【求助】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势函数”并且使用“非周期性边界条件”？
谢谢

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

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人回复

1楼 2009-09-27 13:19:02

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

应助: 26 (小学生)
贵宾: 0.17
金币: 4774.1
散金: 8
红花: 42
帖子: 1361
在线: 215.2小时
虫号: 659094
注册: 2008-11-21
专业: 金属结构材料

★ ★ ★
小木虫(金币+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，再把原子电荷设成真实电荷，运行个几步，比较一下计算的能量结果有多大区别，你就清楚了。