应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 其他 » 【求助】分子动力学计算文章大修，关于审稿意见的回复
91/1返回列表
查看: 2244  |  回复: 8
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

xh512

金虫 (正式写手)

[交流] 【求助】分子动力学计算文章大修，关于审稿意见的回复已有2人参与

对第一性原理分子动力学的基础知识有点缺乏。请各位大大帮忙解答下我的疑问。
我采用VASP进行了水分子在物质表面的第一性原理分子动力学计算，具体做法如下：

平衡计算： 水分子放在表面上，采用 SMASS=-1 ，初始离子的速度设置为0，步长1fs， 平衡计算了2ps

production run： 把平衡计算得到的最终构型的CONTCAR 改为POSCAR， 采用SMASS=0, 计算了3ps。

有几个不自信的地方，请问：
1.  我做平衡计算时，把初始速度设置为0可不肯以呢，平衡后的构型以及速度与初始构型以及离子速度应该没有联系吧？那我这么做应该是可行的吗？
2. 我用平衡后的构型 开始后面的NVT计算，可行吗？

另外，审稿人提出以下几个质疑，然后被拒，需要大修，请各位帮忙看下。

1:  Are 2 ps  long enough to get the surface in equilibrium?
我看到不少文献上采用2 ps来进行平衡，请问各位，我怎么验证我的计算到达了平衡？
2: time step of 1 fs is too large for the motion of H2O molecular.
有不少采用 1 fs 和 0.5 fs步长计算水分子的， 我怎样说服审稿人，1 fs是合理的呢？ 引用文献，还是计算0.5 fs的然后，然后比较两者的什么差别呢？

由于体系比较多，重新采用更小的步长，更长的平衡时间计算，需要很多时间。请大家帮帮忙，我该怎么做。选出一个体系，进行小步长，延长平衡时间来进行对比说明就可以了吗？主要对比哪些方面呢？
非常感谢大家。
回复此楼

» 收录本帖的淘贴专辑推荐

第一原理资料汇编

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2010-08-14 19:41:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★ ★
xh512(金币+8):谢谢您 2010-08-14 21:28:04
zzy870720z(金币+2):谢谢指导 2010-08-14 21:42:00
回答你的问题就可以写一篇文章了。
SMASS=-1是指速度标定吗？
做平衡计算用什么温度控制方法其实并不重要，关键是电子结构要算的对。如果你的温度不高，平衡后的构型和初始构型有很大的关系。初始速度最好还是按Boltzmann分布来给定的好，毕竟你是做的NVT计算。
ps是分子动力学模拟的最小时长，一般而言，ps可以让体系达到平衡。而我想审稿人对H2O的疑问主要是在H键上，H键的震动周期是很短的。你可以从这方面加以考虑。
赞一下(1人) 回复此楼
2楼2010-08-14 20:29:20
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

condensed

木虫 (著名写手)
★ ★
xh512(金币+8):谢谢您 2010-08-14 21:28:19
zzy870720z(金币+2):谢谢指导 2010-08-14 21:42:22
SMASS不是速度标定,这个是选择系宗和热浴方式的.

1.到没到平衡其实很简单.不管是NVE还是NVT,你可以看P,E,T随着时间的变化,如果震荡幅度很小了基本上就能说达到平衡状态了.

2.这个没经验...应该是对于较轻的元素来说,步长应该短一些,特别是H.

3.拿平衡后的位型来重新跑NVT是可以的.但我不确定是否需要保留原来的速度.
赞一下(1人) 回复此楼
3楼2010-08-14 20:55:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xh512

金虫 (正式写手)
引用回帖:
Originally posted by goldenfisher at 2010-08-14 20:29:20:
回答你的问题就可以写一篇文章了。
SMASS=-1是指速度标定吗？
做平衡计算用什么温度控制方法其实并不重要，关键是电子结构要算的对。如果你的温度不高，平衡后的构型和初始构型有很大的关系。初始速度最好还是按 ...

谢谢goldenfisher专家， 看了下动力学相关的理论，决定还是严谨点，调小步长为0.5 fs， 平衡时间延长至5 ps ，虽然时间要花费很多。

还一个比较担心的问题： 我做平衡和NVT计算的时候，都是把物质表面和水分子放一起，一起计算的。   我现在在想，是先单独把表面平衡到所需温度，然后再在平衡后的表面上加入水分子进行NVT模拟好， 还是像我之前那样，水分子和表面一起模拟好呢。 现在有点觉得，先单独平衡表面，再引入水分子可能比较符合实际情况，您觉得怎么做好一点呢？  非常期盼您的回复。谢谢了。
赞一下 回复此楼
4楼2010-08-14 21:35:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xh512

金虫 (正式写手)
引用回帖:
Originally posted by condensed at 2010-08-14 20:55:31:
SMASS不是速度标定,这个是选择系宗和热浴方式的.

1.到没到平衡其实很简单.不管是NVE还是NVT,你可以看P,E,T随着时间的变化,如果震荡幅度很小了基本上就能说达到平衡状态了.

2.这个没经验...应该是对于较轻的 ...

谢谢condensed 的建议。 明白检验平衡状态的了。 依您的经验， 像我的这种体系（ 水分子 在表面的 吸附行为模拟） ， 是1 还是2方案好一点呢？ 谢谢您了。
1。 水分子和表面放一起，先平衡， 再拿平衡后的构型一起NVT。
2。 表面先平衡。 再在平衡后的表面上，引入水分子，然后再以任意速度，进行NVT计算。但是水分子何时引入，以什么样的速度引入，好像也是好麻烦。
赞一下 回复此楼
5楼2010-08-14 21:42:11
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

condensed

木虫 (著名写手)
xh512(金币+2):谢谢您 2010-08-14 21:52:29
介个我没做过...抱歉不能告诉你那个方法更好一些.
赞一下 回复此楼
6楼2010-08-14 21:44:23
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wei.ji

木虫 (正式写手)
★ ★
xh512(金币+1):谢谢参与
xh512(金币+6):谢谢了，意识到这个问题的严重性了，已经在重算了。 2010-08-15 18:22:24
zzy870720z(金币+1):谢谢指教 2010-08-15 20:39:19
怎么说呢。。。
1fs对于H来说的确太大了，2ps对于表面来说的确有点短，审稿人说的很靠谱呀。
引用回帖:
Originally posted by xh512 at 2010-08-14 19:41:10:
对第一性原理分子动力学的基础知识有点缺乏。请各位大大帮忙解答下我的疑问。
我采用VASP进行了水分子在物质表面的第一性原理分子动力学计算，具体做法如下：

平衡计算： 水分子放在表面上，采用 SMASS=-1 ， ...

赞一下(2人) 回复此楼
http://sim.phys.ruc.edu.cn
7楼2010-08-15 18:14:08
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
xh512(金币+6):谢谢您了，我知道怎么做了。 2010-08-15 21:21:44
zxzj05(金币+2):奖励一下 2010-08-16 13:13:50
引用回帖:
Originally posted by xh512 at 2010-08-14 21:35:40:


谢谢goldenfisher专家， 看了下动力学相关的理论，决定还是严谨点，调小步长为0.5 fs， 平衡时间延长至5 ps ，虽然时间要花费很多。

还一个比较担心的问题： 我做平衡和NVT计算的时候，都是把物质表面和水 ...

你说的两种方式都是可行的，平衡后无论你用哪种方式，得到的结果都应该类似。我觉得这要根据你研究的对象来定，如果你的对象是想加温再考虑吸附，就先分别做MD。如果是先吸附再加温，就应该一起做MD.
赞一下(2人) 回复此楼
8楼2010-08-15 21:04:47
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dq20081935

铜虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
问一下，速标是什么？MS软件中的discover模块中有要求设置速标的吗？
赞一下 回复此楼
9楼2012-04-30 11:05:04
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xh512 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件) 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面青地会评时间 +5 tanjydd 2024-06-15 5/250 2024-06-15 23:08 by andywei1028
[教师之家] 关于2023的收入 +33 小龙虾2008 2024-06-10 34/1700 2024-06-15 23:01 by zeolitess
[文学芳草园] 累并快乐着 +13 MYHLD521 2024-06-14 13/650 2024-06-15 22:59 by zeolitess
[教师之家] 每次骚扰女学生的都是院系领导，而不是普通教师，小编们要注意措辞正确 +7 zju2000 2024-06-15 9/450 2024-06-15 22:54 by zeolitess
[基金申请] Nature 11日发文，中国著名学者们称造假迫不得已 +5 babu2015 2024-06-14 5/250 2024-06-15 20:40 by lc231001
[论文投稿] 审稿人含糊拒稿，还需要回复吗？如何回复？ 20+3 BruceChum 2024-06-15 17/850 2024-06-15 20:19 by arthas_007
[找工作] 江西双非一本和四川双一流高校如何选择？ 5+7 寒山敲钟 2024-06-12 22/1100 2024-06-15 20:00 by LittleBush
[基金申请] 博后基金，博管会会提前知道消息吗？ +3 yuyiang 2024-06-13 3/150 2024-06-15 19:52 by Lynn212
[硕博家园] 博士毕业高校和就业的相关问题 +4 SCITOPPP 2024-06-14 6/300 2024-06-15 18:54 by SCITOPPP
[考博] 2025考博 +4 自强不息a?a 2024-06-15 6/300 2024-06-15 12:36 by 自强不息a?a
[食品] 食品博士导师 +6 小李醒yy 2024-06-11 9/450 2024-06-14 23:37 by 小李醒yy
[找工作] 应届大龄土博教职选择，石河子大学和广西大学 10+4 茅山老道士 2024-06-10 31/1550 2024-06-14 10:40 by ChmistryCHEN
[基金申请] 工材E10口函评结束了吗 10+3 我1的飞翔 2024-06-13 5/250 2024-06-14 06:35 by nono2009
[考博] 还能读上博士吗 +14 小伙走起 2024-06-09 24/1200 2024-06-14 06:32 by 18787065176
[教师之家] 请各位老师谈谈身边非升即走的例子 +10 Lucas1121 2024-06-09 25/1250 2024-06-13 21:52 by Lucas1121
[硕博家园] 科研求助 +5 杲www 2024-06-12 6/300 2024-06-13 16:16 by 姓李名明
[有机交流] 原料反应完了，怎么知道是产物还是中间体 +6 小胡在努力 2024-06-11 8/400 2024-06-13 13:33 by 091602
[基金申请] C口食品 +6 好好不过分 2024-06-09 19/950 2024-06-13 12:26 by zydemon525
[论文投稿] with editor日期变更 +3 慎独的小花卷 2024-06-12 8/400 2024-06-13 11:00 by 慎独的小花卷
[找工作] 成都产品质量检测研究院 200+3 鲸鱼663 2024-06-11 8/400 2024-06-13 08:10 by 加纳居士
信息提示
关闭
请填处理意见
关闭