| 查看: 4128 | 回复: 8 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
liuxianlv金虫 (正式写手)
|
[交流]
【求助】请教分子动力学模拟已有5人参与
|
||||
|
本人初学MS，本实验室我是第一个搞这个的，纯粹是摸索前进，对文献上的很多不懂特请教各位大哥大姐，谢谢各位。 分子动力学模拟 选取Fe 晶体(001)晶面为吸附表面[18], 表面体系厚度为11 层共计1584 个Fe 原子, 体系大小为3.20 nm×3.20 nm×6.43 nm. 计算过程中“冻结”表面体系中的Fe 原子, 而吸附分子保持与金属表面自由相互作用. 模拟中使用Materials Studio 4.0 软件包中的COMPASS 力场对体系进行优化, 正则系综(NVT)[19]的MD 模拟通过Discover[20] 模块完成. 模温度为298 K, 采用Andersen[21]恒温器进行温度控制. 各分子起始速度由Maxwell-Boltzmann 分布随机产生. 利用周期性边界条件和时间平均等效于系综平均等基本假设, 运用Velocity Verlet 算法[22]求解牛顿运动方程. 范德瓦尔斯和库仑相互作用采用Charge Group 方法[23]计算. 截断半径为1.5 nm, 截断距离之外的分子间相互作用能按平均密度近似方法进行校正. 模拟中时间步长为1 fs, 模拟时间为500 ps, 每1000 步记录一次体系的轨迹信息. 所有计算均在Pentium IV 计算机上完成 分子在金属表面的吸附能由下列公式计算得出[10]adsorption molecule surface E ＝E ＋E －Etotal (1) 其中, Eadsorption 为吸附能; Emolecule 是自由分子的能量;Esurface 是未吸附分子时金属表面的能量; Etotal 是包含一个分子和金属表面体系的总能量. 通过分析体系能量和温度随时间演化轨迹发现: 100 ps 时体系就开始达到平衡, 后250 ps 时体系已充分平衡, 因此选择后250 ps进行取样平均, 求得各分子吸附能和吸附构型结构参数的统计平均值.（文献来源《1-(2-羟乙基)-2-烷基-咪唑啉缓蚀剂缓蚀机理的理论研究》化 学 学 报2008 年第66 卷，张 军） 我对这段模拟过程和设置细节有些不清楚，现摸索如下：各位帮我看看是不是这样的步骤？ 第一步建了表面按要求。 第二步加入表面分子。 第三步  选discover 里的setup选compass直接再点calculation（是不是不用点 和 以及其他设置）。 第四步选discover 里的Dynimics，   （原文中 各分子起始速度由Maxwell-Boltzmann 分布随机产生. 利用周期性边界条件和时间平均等效于系综平均等基本假设, 运用Velocity Verlet 算法[22]求解牛顿运动方程. 范德瓦尔斯和库仑相互作用采用Charge Group 方法[23]计算. 截断半径为1.5 nm, 截断距离之外的分子间相互作用能按平均密度近似方法进行校正.是怎么设的？） 第五步点NVT下的run. Maxwell-Boltzmann 没找到这东西，在哪设？Velocity Verlet 也没有， 范德瓦尔斯和库仑相互作用采用Charge Group 方法[23]计算. 截断半径为1.5 nm好像在Setup里有不知是不是？vdw&coulomb是指范德瓦尔斯和库仑相互作用吗？以及 截断半径为1.5 nm是不是如图设的：  Charge Group 方法是不是指group based |
回复此楼» 收录本帖的淘贴专辑推荐
 晶体模型计算分析 | 分子动力学 |
» 猜你喜欢
对象没有，还非常想读博，难以抉择
已经有40人回复
-
回家两天，不想呆了
已经有8人回复
-
今年什么时候会评啊
已经有7人回复
-
基你太美
已经有6人回复
-
求助析晶问题
已经有4人回复
-
高校辞职，要求赔偿，这到底合不合理
已经有26人回复
-
博后面上今天下午会公布吗？有无消息？
已经有22人回复
-
评审专家会不会很在意申请人的单位啊
已经有7人回复
-
自然基金变化？？
已经有23人回复
-
博后面上今天下午会公布吗？大家有无消息？
已经有5人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教各位，分子动力学模拟NVE系统，不对速度进行标度可以吗？
已经有7人回复
- 分子动力学模拟长程库仑势的处理问题
已经有8人回复
- 【求助】VASP做分子动力学模拟NVT系综下为啥内能竟然是升高的？
已经有8人回复
- 【转】vasp的分子动力学模拟
已经有38人回复
- 请问各位熟悉分子动力学模拟的朋友，熟悉精通MS的分子动力学模拟需要多久？
已经有17人回复
- 【讨论】关于汽液两相分子动力学模拟中粒子在汽液两相盒子交换的编程问题
已经有3人回复
- 【讨论】大家好，加水后跑完分子动力学，在log中发现了问题
已经有15人回复
- 【求助】分子动力学模拟尺度问题？
已经有11人回复
- 【原创】分子对接，同源模建，分子动力学软件介绍帖
已经有35人回复
- 【求助】gulp 分子动力学
已经有26人回复
- 【求助】新手请教分子模拟问题
已经有27人回复
043114076
木虫 (正式写手)
- 应助: 49 (小学生)
- 金币: 3364.9
- 散金: 52
- 红花: 23
- 帖子: 788
- 在线: 369.5小时
- 虫号: 570592
- 注册: 2008-06-07
- 性别: GG
- 专业: 理论和计算化学
4楼2010-06-03 09:46:33
2楼2010-06-02 21:23:40
liuxianlv
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 868
- 散金: 42
- 帖子: 883
- 在线: 61.4小时
- 虫号: 839985
- 注册: 2009-09-04
- 性别: GG
- 专业: 理论和计算化学
|
可以再回答点吗？（原文中 各分子起始速度由Maxwell-Boltzmann 分布随机产生. 利用周期性边界条件和时间平均等效于系综平均等基本假设, 运用Velocity Verlet 算法[22]求解牛顿运动方程. 范德瓦尔斯和库仑相互作用采用Charge Group 方法[23]计算. 截断半径为1.5 nm, 截断距离之外的分子间相互作用能按平均密度近似方法进行校正.是怎么设的？） 第五步点NVT下的run. Maxwell-Boltzmann 没找到这东西，在哪设？Velocity Verlet 也没有， 范德瓦尔斯和库仑相互作用采用Charge Group 方法[23]计算. 截断半径为1.5 nm好像在Setup里有不知是不是？vdw&coulomb是指范德瓦尔斯和库仑相互作用吗？以及 截断半径为1.5 nm是不是如图设 |
3楼2010-06-02 21:43:27
liuxianlv
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 868
- 散金: 42
- 帖子: 883
- 在线: 61.4小时
- 虫号: 839985
- 注册: 2009-09-04
- 性别: GG
- 专业: 理论和计算化学
5楼2010-06-03 10:44:17
5