MS Discover計(jì)算純表面的問題
一直用MS Discover模塊,沒有注意過這個(gè)問題,直到今天有個(gè)同學(xué)問起。
問題是這樣的,MS Discover自帶的教程,該教程目錄如下
Help --> Tutorials --> Discover tutorials --> Polymer interactions with a metal oxide surface
跟著這個(gè)教程,一步步地做,一直做到附圖1中的紅框部分,此時(shí)教程的意思是,優(yōu)化Al2O3的 (0 0 -1)晶面。按教程里面的操作,是在晶體導(dǎo)入之后,直接切面得(0 0 -1)晶面,然后直接地對(duì)該晶面進(jìn)行優(yōu)化(沒有做諸如添加真空層之類的操作)。
優(yōu)化完后,得到如附圖2的結(jié)果,從附圖2中的紅框可以看到,優(yōu)化前后,該晶面確實(shí)被優(yōu)化了,因?yàn)榭偰芰繌某跏嫉?456左右,下降到了-8532左右。
問題是這樣出現(xiàn)的,既然上面說到該晶面被優(yōu)化了,而且優(yōu)化能量變化得很大,那么應(yīng)該能從優(yōu)化前后的結(jié)構(gòu)中觀察到區(qū)別,觀察結(jié)果如附圖3所示,附圖3中,左邊的晶面是優(yōu)化之后的,右邊的是優(yōu)化之前的,目測來看,看不見差異,一點(diǎn)差異都看不出來?紤]到目測的誤差,于是,計(jì)算優(yōu)化后的晶面的能量,如附圖4中所示,附圖4中,左圖是優(yōu)化前后的能量對(duì)比,右圖優(yōu)化后的晶面的能量,可以看到,優(yōu)化后的晶面的能量,與優(yōu)化前的一模一樣!
也就是說,對(duì) (0 0 -1)晶面作優(yōu)化后,輸出的晶面結(jié)構(gòu)和初始的是一樣的。這個(gè)時(shí)候,我們會(huì)想到,是不是優(yōu)化后的結(jié)構(gòu)沒有及時(shí)地更新呢?同時(shí)考慮到,優(yōu)化計(jì)算的過程時(shí)間很短,所以,在discover的setup中,作了如附圖5的設(shè)置。
結(jié)果還是和上述的一樣:Discover對(duì)純晶面優(yōu)化的時(shí)候,盡管輸出文件.out中,會(huì)給出優(yōu)化前后能量的差異,但是輸出的結(jié)構(gòu)不會(huì)更新。
為了排除教程所涉及體系的特殊性,又對(duì)其他的晶體進(jìn)行了類似計(jì)算,發(fā)現(xiàn)Discover 對(duì)不添加真空層的純表面,全部出現(xiàn)上述情況。
而且,我使用MS6.0 MS6.1兩個(gè)版本重復(fù)了上述計(jì)算,出現(xiàn)一模一樣的情況。
但是一旦給這個(gè)晶面加上真空層,使其從二維變到三維后,結(jié)算后輸出的結(jié)構(gòu)就會(huì)更新。
綜上所述,現(xiàn)在的問題是,Discover對(duì)純晶面優(yōu)化后,輸出的結(jié)構(gòu)不會(huì)更新,今天被這個(gè)問題困擾了。
請(qǐng)使用Discover的朋友給予解答,金幣奉上!

1.png

2.png

3.png

4.png

5.png
[ Last edited by 月只藍(lán) on 2013-10-20 at 18:51 ]
返回小木蟲查看更多
京公網(wǎng)安備 11010802022153號(hào)
郁悶圖片傳不上來,用渣浪貼圖了
圖片地址
https://photo.weibo.com/22702667 ... 2y7b7gdxj211x0k9q6v
可以看到模擬的兩個(gè)結(jié)構(gòu)能量不同
煩屎了煩屎了
圖片不好傳
9L的渣浪圖退出后就看不到
最后試一次,不行就算了
https://ww2.sinaimg.cn/large/87518556gw1ee2y7b7gdxj211x0k9q6v.jpg
非常感謝您時(shí)間的付出。
如果按照教程來,我也發(fā)現(xiàn)“最后結(jié)構(gòu)沒有變化,但是out文件的能量不一樣”,那么是不是可以認(rèn)為,盡管out文件中顯示能量有差異,但實(shí)際上原結(jié)構(gòu)并沒有被優(yōu)化呢?
我覺得是這么理解吧
結(jié)構(gòu)是被優(yōu)化了,只是正好優(yōu)化后的結(jié)構(gòu)和初始結(jié)構(gòu)一致
因?yàn)樵游恢孟拗频脑,所以兩個(gè)結(jié)構(gòu)能量有較大差異
dis里初始結(jié)構(gòu)和最終結(jié)果原子的限定狀態(tài)不一樣
在forcite里,初始結(jié)構(gòu)fix,最終結(jié)構(gòu)也限定了
這點(diǎn)上感覺dis不是很靠譜……
forcite數(shù)據(jù)
1、不限定原子
幾何優(yōu)化前后能量是差異不大
---- 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、限定原子
幾何優(yōu)化前后能量是差異不大
---- 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,
http://www.gaoyang168.com/bbs/viewthread.php?tid=7004165&target=self&page=1
幫忙給這個(gè)帖子看看吧。
昨天我也不知道為啥沒有應(yīng)助選項(xiàng),就一個(gè)“只回復(fù)”和“送紅花”
哦,那是因?yàn)榻饚乓呀?jīng)送完,主題已經(jīng)完結(jié),就只有以上選項(xiàng)。