当前位置: 首页 > 分子模拟 >MS Discover计算纯表面的问题

MS Discover计算纯表面的问题

作者 月只蓝
来源: 小木虫 750 15 举报帖子
+关注

一直用MS Discover模块,没有注意过这个问题,直到今天有个同学问起。
问题是这样的,MS Discover自带的教程,该教程目录如下
Help --> Tutorials --> Discover tutorials --> Polymer interactions with a metal oxide surface
跟着这个教程,一步步地做,一直做到附图1中的红框部分,此时教程的意思是,优化Al2O3的 (0 0 -1)晶面。按教程里面的操作,是在晶体导入之后,直接切面得(0 0 -1)晶面,然后直接地对该晶面进行优化(没有做诸如添加真空层之类的操作)。
优化完后,得到如附图2的结果,从附图2中的红框可以看到,优化前后,该晶面确实被优化了,因为总能量从初始的-456左右,下降到了-8532左右。
问题是这样出现的,既然上面说到该晶面被优化了,而且优化能量变化得很大,那么应该能从优化前后的结构中观察到区别,观察结果如附图3所示,附图3中,左边的晶面是优化之后的,右边的是优化之前的,目测来看,看不见差异,一点差异都看不出来。考虑到目测的误差,于是,计算优化后的晶面的能量,如附图4中所示,附图4中,左图是优化前后的能量对比,右图优化后的晶面的能量,可以看到,优化后的晶面的能量,与优化前的一模一样!
也就是说,对 (0 0 -1)晶面作优化后,输出的晶面结构和初始的是一样的。这个时候,我们会想到,是不是优化后的结构没有及时地更新呢?同时考虑到,优化计算的过程时间很短,所以,在discover的setup中,作了如附图5的设置。

结果还是和上述的一样:Discover对纯晶面优化的时候,尽管输出文件.out中,会给出优化前后能量的差异,但是输出的结构不会更新。

为了排除教程所涉及体系的特殊性,又对其他的晶体进行了类似计算,发现Discover 对不添加真空层的纯表面,全部出现上述情况。
而且,我使用MS6.0 MS6.1两个版本重复了上述计算,出现一模一样的情况。

但是一旦给这个晶面加上真空层,使其从二维变到三维后,结算后输出的结构就会更新。

综上所述,现在的问题是,Discover对纯晶面优化后,输出的结构不会更新,今天被这个问题困扰了。

请使用Discover的朋友给予解答,金币奉上!
MS Discover计算纯表面的问题
1.png


MS Discover计算纯表面的问题-1
2.png


MS Discover计算纯表面的问题-2
3.png


MS Discover计算纯表面的问题-3
4.png


MS Discover计算纯表面的问题-4
5.png

[ Last edited by 月只蓝 on 2013-10-20 at 18:51 ] 返回小木虫查看更多

今日热帖
  • 精华评论
  • pigrass

    引用回帖:
    8楼: Originally posted by pigrass at 2014-03-03 22:57:59
    初步试下来可能是fix cartesian position的原因
    教程中minimize操作解除了之前部分原子cartesian position的fix,minimize前后结构的能量表达式改变,所以虽然结构看起来一样,但是能量计算方式不同,结果不同

    ...

    郁闷图片传不上来,用渣浪贴图了
    图片地址
    https://photo.weibo.com/22702667 ... 2y7b7gdxj211x0k9q6v
    可以看到模拟的两个结构能量不同

  • pigrass

    引用回帖:
    8楼: Originally posted by pigrass at 2014-03-03 22:57:59
    初步试下来可能是fix cartesian position的原因
    教程中minimize操作解除了之前部分原子cartesian position的fix,minimize前后结构的能量表达式改变,所以虽然结构看起来一样,但是能量计算方式不同,结果不同

    ...

    烦屎了烦屎了
    图片不好传
    9L的渣浪图退出后就看不到
    最后试一次,不行就算了
    https://ww2.sinaimg.cn/large/87518556gw1ee2y7b7gdxj211x0k9q6v.jpg

  • 月只蓝

    引用回帖:
    8楼: Originally posted by pigrass at 2014-03-03 22:57:59
    初步试下来可能是fix cartesian position的原因
    教程中minimize操作解除了之前部分原子cartesian position的fix,minimize前后结构的能量表达式改变,所以虽然结构看起来一样,但是能量计算方式不同,结果不同

    ...

    非常感谢您时间的付出。
    如果按照教程来,我也发现“最后结构没有变化,但是out文件的能量不一样”,那么是不是可以认为,尽管out文件中显示能量有差异,但实际上原结构并没有被优化呢?

  • pigrass

    引用回帖:
    11楼: Originally posted by 月只蓝 at 2014-03-04 08:27:01
    非常感谢您时间的付出。
    如果按照教程来,我也发现“最后结构没有变化,但是out文件的能量不一样”,那么是不是可以认为,尽管out文件中显示能量有差异,但实际上原结构并没有被优化呢?...

    我觉得是这么理解吧
    结构是被优化了,只是正好优化后的结构和初始结构一致
    因为原子位置限制的原因,所以两个结构能量有较大差异

  • pigrass

    dis里初始结构和最终结果原子的限定状态不一样
    在forcite里,初始结构fix,最终结构也限定了
    这点上感觉dis不是很靠谱……

    forcite数据
    1、不限定原子
    几何优化前后能量是差异不大
    ---- Initial structure ----
    Total enthalpy                  : -21161.553598 kcal/mol
      External pressure term        : 0.000000 kcal/mol
    Total energy                    :      -21161.553598 kcal/mol
    ---- Final structure ----
    Total enthalpy                  : -21218.565110 kcal/mol
      External pressure term        : 0.000000 kcal/mol
    Total energy                    :      -21218.565110 kcal/mol
    2、限定原子
    几何优化前后能量是差异不大
    ---- Initial structure ----
    Total enthalpy                  : -4791.154883 kcal/mol
      External pressure term        : 0.000000 kcal/mol
    Total energy                    :       -4791.154883 kcal/mol
    ---- Final structure ----
    Total enthalpy                  : -4819.258852 kcal/mol
      External pressure term        : 0.000000 kcal/mol
    Total energy                    :       -4819.258852 kcal/mol,

  • 月只蓝

    引用回帖:
    13楼: Originally posted by pigrass at 2014-03-04 09:53:35
    dis里初始结构和最终结果原子的限定状态不一样
    在forcite里,初始结构fix,最终结构也限定了
    这点上感觉dis不是很靠谱……

    forcite数据
    1、不限定原子
    几何优化前后能量是差异不大
    ---- Initial structure  ...

    https://muchong.com/bbs/viewthread.php?tid=7004165&target=self&page=1
    帮忙给这个帖子看看吧。

  • pigrass

    昨天我也不知道为啥没有应助选项,就一个“只回复”和“送红花”

  • 月只蓝

    引用回帖:
    15楼: Originally posted by pigrass at 2014-03-04 13:04:22
    昨天我也不知道为啥没有应助选项,就一个“只回复”和“送红花”

    哦,那是因为金币已经送完,主题已经完结,就只有以上选项。

猜你喜欢
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓