24小时热门版块排行榜    

查看: 2083  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

混乱的流动

银虫 (小有名气)

[求助] Ewald summation求助 已有1人参与

请问谁有Ewald summation 的程序段,我想用一下。实现以下目的:在三维计算区域内有1000个球体,带电量为q,其中y-z方向是周期性边界条件,x方向没有周期性。用Ewald 应该能算出库伦静电作用势或者静电作用力吧?
回复此楼

» 猜你喜欢

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

lsloneil

专家顾问 (正式写手)

引用回帖:
3楼: Originally posted by 混乱的流动 at 2015-11-09 00:02:28
也就是正负电荷量差不多才行?为什么有这样的特性?
...

是的。因为如果总电荷不为0的话,Ewald summation在k空间求和那一部分的级数会发散。通常为了避免这种问题,分子模拟的程序在处理非电中性系统的Ewald summation是都会给系统加上一个均匀分布的背景电荷密度来中和这部分电荷。如果你的体系净电荷不大(比如说只有-2 e),那么这个背景电荷的影响是比较小的,反之,如果你的体系净电荷很大,那么Ewald summation带来的误差就会很大。
4楼2015-11-09 23:24:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

lsloneil

专家顾问 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
如果使用lammps,可以使用kspace_modify slab这个命令,详细说明参见
http://lammps.sandia.gov/doc/kspace_modify.html

不过说明一下,Ewald summation适用的是电中性体系,如果你的体系非中性,贸然使用Ewald summation会带来误差,而且体系总电荷距离0越远,误差越大。
2楼2015-11-09 16:09:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

混乱的流动

银虫 (小有名气)

引用回帖:
2楼: Originally posted by lsloneil at 2015-11-09 16:09:41
如果使用lammps,可以使用kspace_modify slab这个命令,详细说明参见
http://lammps.sandia.gov/doc/kspace_modify.html

不过说明一下,Ewald summation适用的是电中性体系,如果你的体系非中性,贸然使用Ewald ...

也就是正负电荷量差不多才行?为什么有这样的特性?

发自小木虫IOS客户端
3楼2015-11-09 20:02:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuntao118

新虫 (著名写手)

引用回帖:
2楼: Originally posted by lsloneil at 2015-11-09 16:09:41
如果使用lammps,可以使用kspace_modify slab这个命令,详细说明参见
http://lammps.sandia.gov/doc/kspace_modify.html

不过说明一下,Ewald summation适用的是电中性体系,如果你的体系非中性,贸然使用Ewald ...

如果是谁分子体系的话需要注意些什么?在这里你说的体系是非中性指的是分子本身?还是系统本身是中性、系统的分子可以是极性分子(比如:水分子体系)?
5楼2015-11-10 08:26:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 323求调剂 +6 洼小桶 2026-03-18 6/300 2026-03-23 00:29 by king123!
[考研] 一志愿北京化工大学 070300 学硕 336分 求调剂 +4 vv迷 2026-03-22 4/200 2026-03-22 23:29 by king123!
[考研] 一志愿东华大学化学070300,求调剂 +7 2117205181 2026-03-21 8/400 2026-03-22 22:55 by chixmc
[考研] 277材料科学与工程080500求调剂 +7 自由煎饼果子 2026-03-16 7/350 2026-03-22 22:40 by ACS Nano——
[考研] 石河子大学(211、双一流)硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 324求调剂 +6 lucky呀呀呀鸭 2026-03-20 6/300 2026-03-22 16:01 by ColorlessPI
[考研] 289求调剂 +7 怀瑾握瑜l 2026-03-20 7/350 2026-03-22 15:57 by ColorlessPI
[考研] 319求调剂 +4 小力气珂珂 2026-03-20 4/200 2026-03-22 15:53 by ColorlessPI
[考研] 311求调剂 +3 26研0 2026-03-20 3/150 2026-03-22 14:46 by ColorlessPI
[考研] 286分人工智能专业请求调剂愿意跨考! +4 lemonzzn 2026-03-17 8/400 2026-03-21 22:49 by lemonzzn
[考研] 工科0856求调剂 +3 沐析汀汀 2026-03-21 3/150 2026-03-21 18:30 by 学员8dgXkO
[考研] 材料工程(专)一志愿985 初试335求调剂 +3 hiloiy 2026-03-17 4/200 2026-03-21 03:04 by JourneyLucky
[考研] 二本跨考郑大材料306英一数二 +3 z1z2z3879 2026-03-17 3/150 2026-03-21 02:29 by JourneyLucky
[考研] 一志愿华中科技大学,080502,354分求调剂 +5 守候夕阳CF 2026-03-18 5/250 2026-03-21 01:06 by JourneyLucky
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +4 @taotao 2026-03-19 4/200 2026-03-20 22:14 by JourneyLucky
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
信息提示
请填处理意见