24小时热门版块排行榜    

查看: 1767  |  回复: 14
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zdz103979358

新虫 (初入文坛)

[求助] Fe(001)表面自组装成膜分子力学计算

模拟内容如下:采用分子力学方法研究咪唑啉缓蚀剂分子在金属表面成膜机制,选取Fe(001)表面为模拟的金属表面,假定在Fe(001)表面的覆盖度为0.64nm/个,即咪唑啉在Fe(001)表面达到饱和吸附。模拟体系大小为:3.24nm×3.24nm×6.03nm,包含16个咪唑啉缓蚀剂分子。模拟力场选择CAMPASS,计算方法选择Smart Minimizer ,采用Ultra-fine收敛精度,模拟的步数设为20000步。
    我的步骤是:1.Import/structure/metals/pure-metals/ Fe晶体;
2.切Fe(001)表面,设置深度为6:Build/Surface/Cleave   surface,Cleave  plane(h,k,l)输 入001,Depth输入6;
3.建立超晶胞:Build/Symmetry/Supercell,u=11,v=11.
4.去周期性:Build /Symmetry/non-periodic superstructure.
5.导入已经用高斯优化好的咪唑啉分子,并通过复制粘贴,使16个缓蚀剂分子均匀分布在Fe(001)表面。头接近Fe表面,尾部背离表面,16个缓蚀剂分子直立在Fe表面。
6.固定Fe。选定所有Fe原子,Modify/Constraints/点击Fix Carsesian position。
7.开始优化:Modules/Discover/Minimizer,计算方法选择Smart Minimizer ,采用Ultra-fine收敛精度,模拟的步数设为20000步。
    我的问题是:优化完成后,文献上说缓蚀剂分子吸附在Fe表面,形成吸附角,与Fe表面形成一定的覆盖度。但是我模拟出来的图形是几乎没有吸附角,缓蚀剂仍然成直立,没有吸附角,且分布乱七八招的。
    这东西我弄很难久了,谁懂的教一下,在下感激不尽。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

visualization

金虫 (小有名气)

从最后结果看,感觉没有运行多长时间啊
力场中怎么设置的电荷
我也想知道为什么要去掉周期性
7楼2012-05-08 22:48:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

zdz103979358

新虫 (初入文坛)

谁能帮我看看我哪里错了

文献上的结果,缓蚀剂分子成一定的吸附角,排列规律



我构造的原图



我计算出来的结果,几乎没有吸附角,排列不规律。

2楼2012-05-04 11:56:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

geniusew

金虫 (小有名气)

【答案】应助回帖


感谢参与,应助指数 +1
chaizhm: 金币+1, 谢谢~ 2012-05-04 19:40:53
选择一个分子做下模拟试试。。不要16个
3楼2012-05-04 13:20:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zdz103979358

新虫 (初入文坛)

谁懂的帮忙解决一下,急~~
4楼2012-05-05 12:06:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见