| 查看: 5323 | 回复: 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 晶体模型计算分析 | 分子动力学 |
» 猜你喜欢
什么时候开奖？
已经有11人回复
-
范进中举一文的中心思想
已经有6人回复
-
建议基金发布提前给出明确的时间点
已经有14人回复
-
2026国自然函评费到账
已经有17人回复
-
让我中一个面上吧！
已经有15人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有5人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
今天放榜吗？
已经有16人回复
-
93BebMhtakh前后11位开头都是大写
已经有8人回复
-
只有每年这种时候来逛逛小木虫
已经有26人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教各位，分子动力学模拟NVE系统，不对速度进行标度可以吗？
已经有7人回复
- 分子动力学模拟长程库仑势的处理问题
已经有8人回复
- 【求助】VASP做分子动力学模拟NVT系综下为啥内能竟然是升高的？
已经有8人回复
- 请问各位熟悉分子动力学模拟的朋友，熟悉精通MS的分子动力学模拟需要多久？
已经有17人回复
- 【讨论】关于汽液两相分子动力学模拟中粒子在汽液两相盒子交换的编程问题
已经有3人回复
- 【讨论】大家好，加水后跑完分子动力学，在log中发现了问题
已经有15人回复
- 【求助】分子动力学模拟尺度问题？
已经有11人回复
- 【原创】分子对接，同源模建，分子动力学软件介绍帖
已经有35人回复
- 【求助】新手请教分子模拟问题
已经有27人回复
倚天舰
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1487.1
- 散金: 4
- 红花: 3
- 帖子: 82
- 在线: 106小时
- 虫号: 858328
- 注册: 2009-09-27
- 专业: 生物化工与食品化工
9楼2011-09-30 18:22:43
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
043114076
木虫 (正式写手)
- 应助: 49 (小学生)
- 金币: 3364.9
- 散金: 52
- 红花: 23
- 帖子: 788
- 在线: 369.5小时
- 虫号: 570592
- 注册: 2008-06-07
- 性别: GG
- 专业: 理论和计算化学
4楼2010-06-03 09:46:33
40