24小时热门版块排行榜    

查看: 3734  |  回复: 13
本帖产生 1 个 模拟EPI ,点击这里进行查看

awaken2013

铜虫 (小有名气)

引用回帖:
10楼: Originally posted by zhangmao511 at 2015-01-20 22:12:12
感谢你的指点,深深的表示谢意!!
Gromacs手册我总是需要的时候看一下,对某个概念不懂的时候看一下,还没完整的看过一遍。以后抽时间看。
我看了一下手册,这是我个人的理解,不知道对不?? ...

rlist=rc=rw的时候,你的理解完全正确。
当rlist《rw=rc的时候,已经是双程截断的概念了。此时rlist内的原子计算非键相互作用,如键角力 键伸缩 二面角等,短程力(但是不包括库仑力和范德华力,这两个力被分离出去)此时rlist的力就算完了。
(这里开始可以抛开rlist的概念)
rc和rw是范德华力的截断,在此范围内的范德华力和库仑力是短程范围内的力,由此计算,属于短程力。
rlistlong一般比rw和rc大0.3.。而且在mdp中会被自动设置为-1 自动调整。当两个原子距离等于rlistlong时 短程的力被开关函数调节为0.
在rlistlong之外的力,按照P长程力计算。。
11楼2015-01-20 23:14:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmao511

铁虫 (小有名气)

引用回帖:
11楼: Originally posted by awaken2013 at 2015-01-20 23:14:13
rlist=rc=rw的时候,你的理解完全正确。
当rlist《rw=rc的时候,已经是双程截断的概念了。此时rlist内的原子计算非键相互作用,如键角力 键伸缩 二面角等,短程力(但是不包括库仑力和范德华力,这两个力被分离出 ...

P长程力,指的是……??
谢谢讲解,顺便弱弱的问下,能加你好友不??
12楼2015-01-21 11:31:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmao511

铁虫 (小有名气)

引用回帖:
11楼: Originally posted by awaken2013 at 2015-01-20 23:14:13
rlist=rc=rw的时候,你的理解完全正确。
当rlist《rw=rc的时候,已经是双程截断的概念了。此时rlist内的原子计算非键相互作用,如键角力 键伸缩 二面角等,短程力(但是不包括库仑力和范德华力,这两个力被分离出 ...

你好,想问一下。下面这段话摘自一篇文献的method,想问一下最后一句话是指静电相互作用cutoff 是10A,Lennard-Jones interactions cutoff 是12A吗???Lennard-Jones是指的什么呢?谢谢。。
The MD simulations were conducted with NAMD (28) and performed with periodic boundary conditions at a constant pressure (1 atm) and temperature (300K) using Langevin dynamics with Nosé-Hoover Langevin piston pressure control. The long-range electrostatic interactions were calculated using the particle mesh Ewald algorithm. Both the electrostatic and the Lennard-Jones interactions had a twin-range cutoff of 10-12 A.
13楼2015-01-21 16:18:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

awaken2013

铜虫 (小有名气)

引用回帖:
13楼: Originally posted by zhangmao511 at 2015-01-21 16:18:41
你好,想问一下。下面这段话摘自一篇文献的method,想问一下最后一句话是指静电相互作用cutoff 是10A,Lennard-Jones interactions cutoff 是12A吗???Lennard-Jones是指的什么呢?谢谢。。
The MD simulations ...

我之前的解释和理解有误。以下是一位老师回复我的。

rlist、cutoff这些都是对范德华、静电这两类非键相互作用而言的,bond,angle,dihedral,1-4那些和成键有关的作用与此没有丝毫关系。

当库仑和范德华作用都使用cutoff时,rlist小于rcoulomb和rvdw时是twin-range,即rlist内的作用每步都计算,rlist与rcoulomb和rvdw之间的是每nstlist步更新一次。更远部分的作用完全忽略。

PME计算静电作用的情况与此截然不同。只有一个cutoff,rlist必须等于rcoulomb,在rcoulomb以内的在实空间计算,在rcoulomb以外的在倒易空间计算。都是每一步都计算,而且这样处理静电作用是精确的,不管多远的静电相互作用都会被计算,因此PME是计算周期性体系静电相互作用最理想的方法。而且rcoulomb在理论上只影响计算效率而不影响计算结果。

switching function在一般的研究中没用。PME时静电作用本来就是精确的,对于范德华作用,由于收敛速度远比静电作用快,rvdw设得大一些也就基本计算精确了。

但是对于你的文献引用的情况我也不知了 慢慢思考吧
14楼2015-01-24 22:15:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhangmao511 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 296求调剂 +5 大口吃饭 身体健 2026-03-13 5/250 2026-03-17 21:05 by 不惑可乐
[考研] 材料与化工求调剂 +6 为学666 2026-03-16 6/300 2026-03-17 20:15 by peike
[考研] 341求调剂 +5 捣蛋猪猪 2026-03-11 7/350 2026-03-17 19:09 by 捣蛋猪猪
[考研] 302求调剂 +9 负心者当诛 2026-03-11 9/450 2026-03-17 17:13 by ruiyingmiao
[考研] 085600材料与化工求调剂 +5 绪幸与子 2026-03-17 5/250 2026-03-17 16:40 by laoshidan
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[考研] 一志愿,福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 机械专硕325,寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 环境工程调剂 +6 大可digkids 2026-03-16 6/300 2026-03-16 17:16 by barlinike
[考研] 0703化学调剂,求各位老师收留 +8 秋有木北 2026-03-14 8/400 2026-03-16 15:21 by 哦哦123
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 255求调剂 +3 李嘉慧, 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[考研] 328求调剂 +3 5201314Lsy! 2026-03-13 6/300 2026-03-14 15:31 by hyswxzs
[考研] 0703化学调剂 +4 快乐的香蕉 2026-03-11 4/200 2026-03-13 22:41 by JourneyLucky
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
信息提示
请填处理意见