24小时热门版块排行榜    

查看: 1597  |  回复: 15

CMA17

新虫 (初入文坛)

[求助] 高分子模板剂诱导晶体生长的模拟

目前有实验数据证明高分子模板剂可以诱导晶体按照一个方向生长,现在想要用MS来模拟一下这个过程,请问模拟过程要怎样设计一下?或者有相关文献推荐也非常感谢。还请前辈指导
PS:目前金币较少,我马上跟别人借点金币过来补充。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

clickme

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
CMA17: 金币+20, ★★★★★最佳答案 2017-06-19 08:41:29
计算模板剂在不同晶面上的吸附能差异。

[ 发自手机版 http://muchong.com/3g ]
2楼2017-06-18 18:06:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CMA17

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by clickme at 2017-06-18 18:06:41
计算模板剂在不同晶面上的吸附能差异。

谢谢您的回复。我确实也是这么做了,但是老觉得写在文章里,这个过程跟实验数据分析是反着的。您那边有符合得比较好的文献吗?我可能关键词搜错了并没有找到关于模板剂诱导晶体生长的文献,只是找到了一些与我的晶体相同的,不同高分子吸附的文献。很期待您的指导
3楼2017-06-19 08:50:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

clickme

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
CMA17: 金币+10, ★★★★★最佳答案 2017-06-20 10:30:21
如果可以的话做一下频率分析,得到零点能和室温下的热矫正,可能你的结论就能正过来了。

[ 发自手机版 http://muchong.com/3g ]
4楼2017-06-19 12:11:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CMA17

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by clickme at 2017-06-19 12:11:30
如果可以的话做一下频率分析,得到零点能和室温下的热矫正,可能你的结论就能正过来了。

有这两个能量了然后怎样分析呢?您能不能详细地讲讲,例如我现在算出来各晶面吸附能的大小是E(001)>E(100)>E(111)。还有用哪个模块计算可以得到高分子在金属氧化物表面吸附能并且能做频率分析?先谢谢您的耐心指导
5楼2017-06-20 10:36:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

clickme

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
CMA17: 金币+10, ★★★★★最佳答案 2017-06-20 15:03:51
你结构优化得到的能量只是绝对零度下的电子能而已,还应该加上零点震动能和键的伸缩摇摆振动在室温下蕴含的能量,这些都加在一起再比较吸附能,就有可能把你之前算的能量顺序完全推翻。你只要在计算的时候进行频率分析计算就能得到这两个值。

[ 发自手机版 http://muchong.com/3g ]
6楼2017-06-20 14:27:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CMA17

新虫 (初入文坛)

引用回帖:
6楼: Originally posted by clickme at 2017-06-20 14:27:05
你结构优化得到的能量只是绝对零度下的电子能而已,还应该加上零点震动能和键的伸缩摇摆振动在室温下蕴含的能量,这些都加在一起再比较吸附能,就有可能把你之前算的能量顺序完全推翻。你只要在计算的时候进行频率分 ...

谢谢您的回答,我用的是MS中的Discover模块计算的,在Dynamics的过程中已经设置了温度(为实验的反应温度),并且根据吸附能公式“吸附能=体系总能量-(高分子能量+晶体表面能量)”计算得到了吸附能,这种情况是否还需要加上零点振动能和键的能量?
另外,根据我目前得到的计算结果,在(001)的吸附能最大,正好符合我的实验数据中晶体沿[001]方向生长的结果。我之前说的计算过程与实验过程是反的是指,实验过程中是晶体与高分子发生晶格匹配,然后开始沿[001]方向生长,而计算过程是让高分子去吸附晶面,我是觉得这两个过程是反的,并不是说计算与实验结果相反,您能帮我解释一下这边的疑惑吗?期待您的指导
7楼2017-06-20 16:27:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

clickme

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
CMA17: 金币+10, ★★★★★最佳答案 2017-06-21 09:36:47
引用回帖:
7楼: Originally posted by CMA17 at 2017-06-20 16:27:16
谢谢您的回答,我用的是MS中的Discover模块计算的,在Dynamics的过程中已经设置了温度(为实验的反应温度),并且根据吸附能公式“吸附能=体系总能量-(高分子能量+晶体表面能量)”计算得到了吸附能,这种情况是否 ...

好的,我明白你的意思了。
1. 我之前说的零点能校正是对应与量子化学计算方法需要考虑的问题,你用MD的话就可以不用考虑了。
2. 我用的是MS2017, 我看了下Forcite模块(新版中Discover与Forcite合并和),没有声字或者频率分析的选项,所有温度校正也不好再加入了。
3. 你是否用Dynamics后的能量相减得到吸附能?如果是这样的话有不妥,因为此时的结构还有温度热涨落导致的形变能。应该把dynamics的结构再做一此结构优化。
4. 我对晶体生长不大熟悉。如果高分子与晶体某个面结合力强,就能让这个面在生长过程中保持的久一些,最后在表征的时候这个面就多,所以看起来像是沿着这个面生长?我猜的。
8楼2017-06-20 16:46:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CMA17

新虫 (初入文坛)

引用回帖:
8楼: Originally posted by clickme at 2017-06-20 16:46:08
好的,我明白你的意思了。
1. 我之前说的零点能校正是对应与量子化学计算方法需要考虑的问题,你用MD的话就可以不用考虑了。
2. 我用的是MS2017, 我看了下Forcite模块(新版中Discover与Forcite合并和),没有声 ...

谢谢您的指导!

目前我做模拟计算的过程是:1、晶体切面,在Discover中加力场;2、固定下层原子,在Discover中进行Minimizer,构建超胞并加真空层;3、构建高分子,在Forcite中进行结构优化;4、将高分子放在晶体表面,固定晶体所有原子,进行动力学计算;5、最后能量计算。

是不是要改成:1、晶体切面,在Discover中加力场;2、固定下层原子,在Discover中进行Minimizer,构建超胞并加真空层;3、构建高分子,在Forcite中进行结构优化;4、将高分子放在晶体表面,固定晶体所有原子,将整个结构在Forcite中进行结构优化;5、将整个结构进行动力学计算;6、最后能量计算。

修改后的计算过程我有2个问题:1、将整个结构进行结构优化的过程并没有加入温度场,这样是否能满足温度涨落导致的形变能的计算。在目前的计算过程中,动力学计算过程未固定高分子,且计算过程中有温度参与,是否已经包括温度造成的分子形变能;2、动力学计算之前将整个体系进行结构优化,当不同晶面吸附同一高分子时,由于晶面结构不同,结构优化后的高分子在不同晶面上的结构一定不同,这样的两个体系算出的吸附能是否有可比性。就是说,进行动力学计算的高分子组成不变但构型不同,这样得到的吸附能能不能算同高分子对晶体不同晶面的吸附能?

问题比较多,期待您的再次耐心指导,谢谢
9楼2017-06-21 10:12:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

clickme

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
CMA17: 金币+10, ★★★★★最佳答案 2017-06-21 15:11:24
原则是这样,用带温度的分子动力学模拟让高分子链找到大致的全局最构型,在后优再化一下就得到了全局最优结构的精确能量。所以你只要把你之前得到的构型再优化下即可。到要保证之前动力学足够久,高分子得到了充分驰预,找到了自己最好的大致吸附构型。所以动力学模拟阶段温度不宜太高,五六百k即可。

[ 发自手机版 http://muchong.com/3g ]
10楼2017-06-21 11:28:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 CMA17 的主题更新
信息提示
请填处理意见