24小时热门版块排行榜    

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

jane8556

金虫 (正式写手)

[求助] 困扰很久的难题-结构变形严重-SiO2 已有4人参与

大家好,这个问题困扰了我很长时间,请大家帮我看看是哪里的问题吧,多谢。
问题是这样的,我想模拟水溶液环境下纳米物质在SiO2表面的吸附。但是程序运行的时候总是丢失原子,而且我发现是SiO2变形较大,因而我首先判断是SiO2的问题。所以我将SiO2单独拿出来模拟。首先怀疑是两方面原因,一个是初始模型不合理,因而我试用了MS中的alpha-quartz模型及VMD中的SiO2模型,同时也将模型的大小做了一些调整,例如试用了15*15*4,14*14*4及10*10*4大小的模型。另一个是势函数不合理,因而我试用了BKS,tersoff及comb势,其中BKS也采用了table。但是结果仍然不合理。我将其中一组模拟结果黏贴在附件里,有两张图片,第一张是初始构型,第二张是运行50000步之后的结构。请大家帮忙,谢谢!图片直接上传不上去,给大家造成的麻烦请谅解
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 1.jpg
  • 2015-07-17 16:13:45, 688.83 K
  • 附件 2 : 2.jpg
  • 2015-07-17 16:13:53, 688.83 K

» 猜你喜欢

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

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

oooofish

新虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
jane8556: 金币+10, ★★★很有帮助 2015-09-17 09:59:04
Can't type Chinese in this computer, not on purpose. I also think something must be wrong in the input file. I simulated SiO2 by BKS and Charmm force field, and the structure was ok. And also, can you visualize the structure by using VDW? It's hard to see what happens from the snapshot you pasted.
7楼2015-07-21 02:50:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

chuntao118

新虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
jane8556: 金币+20, ★★★很有帮助 2015-09-17 09:46:50
模拟原子丢失应该是边界条件处理的问题吧!个人认为和势函数以及模拟区域没有关系。因为在势函数的选取可以说别人给定的,模拟区域的大小只影响计算的模拟时间,怎么都不会影响到模拟体系的原子数目的变化的。道理上就讲不通,建议:不要意味的用软件,好好的把道理想明白。目的搞清楚。软件中的细节弄明白。一知半解自然错误多多。

» 本帖已获得的红花(最新10朵)

2楼2015-07-18 08:39:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jane8556

金虫 (正式写手)

送红花一朵
引用回帖:
2楼: Originally posted by chuntao118 at 2015-07-18 08:39:15
模拟原子丢失应该是边界条件处理的问题吧!个人认为和势函数以及模拟区域没有关系。因为在势函数的选取可以说别人给定的,模拟区域的大小只影响计算的模拟时间,怎么都不会影响到模拟体系的原子数目的变化的。道理上 ...

谢谢你,你好专业,上面我可能没说清楚,刚开始是原子丢失,最后将盒子大小变大一点就不出现这个问题了,但是SiO2结构的变形非常大,所以我猜测是势函数的问题,你觉得呢?
3楼2015-07-18 09:46:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

声梦奇缘001

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
jane8556: 金币+10, ★★★很有帮助 2015-09-17 09:47:10
你可以在模拟的时候对SiO2施加一个约束(即施加一个简谐势能),因为你模拟的是在表面的吸附过程,只要表面是不受约束的即可
只做不说,明知没地位,坚信有机会
4楼2015-07-18 12:57:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见