24小时热门版块排行榜    

查看: 4908  |  回复: 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势函数”并且使用“非周期性边界条件”?
谢谢
回复此楼

» 猜你喜欢

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

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

老虎大王

木虫 (著名写手)

★ ★ ★
小木虫(金币+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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

老虎大王

木虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjcmwgk(金币+1,VIP+0): 9-27 21:15
计算库仑势的EWald方法不适用于非周期性边界条件,详见说明书里对kspace_style的说明。
2楼2009-09-27 20:35:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bs340248

铜虫 (小有名气)

我感觉貌似要用非周期性边界条件一定要用kspace_style none
但是要用Born-Mayer-Huggins势函数一定不能用kspace_style none
如何办啊?我哪里理解错了?
谢谢你的回复
3楼2009-09-27 23:41:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Camilleyuxia

铜虫 (初入文坛)


小木虫(金币+0.5):给个红包,谢谢回帖交流
先收了,以后我肯定用得着
好好学习,天天向上
5楼2011-03-21 21:38:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学调剂 +8 妮妮ninicgb 2026-03-15 12/600 2026-03-19 05:55 by anny19840123
[考研] 271材料工程求调剂 +5 .6lL 2026-03-18 5/250 2026-03-19 03:07 by 无懈可击111
[考研] 一志愿华中科技大学,080502,354分求调剂 +4 守候夕阳CF 2026-03-18 4/200 2026-03-18 22:16 by li123456789.
[考研] 085410人工智能专硕317求调剂(0854都可以) +3 xbxudjdn 2026-03-18 3/150 2026-03-18 22:14 by zhq0425
[考研] 298-一志愿中国农业大学-求调剂 +7 手机用户 2026-03-17 7/350 2026-03-18 14:34 by vgtyfty
[考研] 一志愿西南交大,求调剂 +4 材化逐梦人 2026-03-18 4/200 2026-03-18 14:22 by 007_lilei
[考研] 280求调剂 +6 咕噜晓晓 2026-03-18 7/350 2026-03-18 11:25 by 无际的草原
[考研] 303求调剂 +4 睿08 2026-03-17 6/300 2026-03-18 11:01 by Iveryant
[基金申请] 被我言中:新模板不强调格式了,假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 268求调剂 +8 一定有学上- 2026-03-14 9/450 2026-03-17 17:47 by laoshidan
[考研] 梁成伟老师课题组欢迎你的加入 +8 一鸭鸭哟 2026-03-14 10/500 2026-03-17 15:07 by 一鸭鸭哟
[考研] 290求调剂 +6 孔志浩 2026-03-12 11/550 2026-03-17 14:41 by 周舟舟77
[考研] 211本,11408一志愿中科院277分,曾在中科院自动化所实习 +6 Losir 2026-03-12 7/350 2026-03-17 12:09 by danranxie
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 工科278分求调剂 +5 周慢热啊 2026-03-12 7/350 2026-03-13 15:49 by JourneyLucky
[考研] 一志愿211化学学硕310分求调剂 +8 努力奋斗112 2026-03-12 9/450 2026-03-13 15:41 by JourneyLucky
信息提示
请填处理意见