| 查看: 4526 | 回复: 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 晶体模型计算分析 | 分子动力学 |
» 猜你喜欢
孩子确诊有中度注意力缺陷
已经有12人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
2026博士申请-功能高分子，水凝胶方向
已经有6人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
-
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复
-
同一篇文章，用不同账号投稿对编辑决定是否送审有没有影响？
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
请教各位，分子动力学模拟NVE系统，不对速度进行标度可以吗？
已经有7人回复- 分子动力学模拟长程库仑势的处理问题
已经有8人回复
- 【求助】VASP做分子动力学模拟NVT系综下为啥内能竟然是升高的？
已经有8人回复
- 请问各位熟悉分子动力学模拟的朋友，熟悉精通MS的分子动力学模拟需要多久？
已经有17人回复
- 【讨论】关于汽液两相分子动力学模拟中粒子在汽液两相盒子交换的编程问题
已经有3人回复
- 【讨论】大家好，加水后跑完分子动力学，在log中发现了问题
已经有15人回复
- 【求助】分子动力学模拟尺度问题？
已经有11人回复
- 【原创】分子对接，同源模建，分子动力学软件介绍帖
已经有35人回复
- 【求助】新手请教分子模拟问题
已经有27人回复
8楼2010-09-21 09:12:03
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