24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3485  |  回复: 8
【悬赏金币】回答本帖问题,作者HJieY将赠送您 15 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

HJieY

新虫 (初入文坛)

[求助] MD simulation之后处理ovito中遇到的问题

前辈们,朋友们:
       你们好!
       我的研究方向是采用分子动力学模拟硬脆性材料在纳米加工中的材料去除机理及变形行为。最近在采用ovito尝试对模拟结果进行后处理时,出现了一些问题不知如何解决,特别希望热心的虫友可以为我释惑。问题如下:
       1、在ovito中进行配位数分析,在视觉效果上感觉配位数设置无误后(截断半径我设置为稍大于键长), 进行成键时却发现配位数与中心原子的键数十分不一致,如5配位原子的中心原子键数应为5才对,但键数却表现的极为任意(0~5个),原子间没有键居多,成5键的原子几乎没有,真的很奇怪;
       2、在SiC材料中,在Si原子和C原子成键时,C原子和C原子、Si原子和Si原子之间是否会成键呢?
       先行感谢致力于研究的你们,盼解惑!
                                                                                                                                                                                         一个在分子动力学上处于迷茫的人儿
回复此楼
人在塔在
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

HJieY

新虫 (初入文坛)

送红花一朵
引用回帖:
6楼: Originally posted by get-it at 2021-07-29 19:46:57
文献的rdf很奇怪。对于周期性体系,截断半径趋于无穷时,g(r)趋于1(你可以找些典型的rdf看看)。但是这里峰值都没到0.5,你可能要看看文献中rdf的公式才行。

前辈,
      您好:
      受到你的回复,我很高兴!我也感觉RDF挺奇怪的,但不只这一篇论文是这样,我看到的另一篇论文也是这样的,RDF峰值最大为1左右(如第二、三图所示),但就像您说的那样,截断半径不可能无穷大(此篇论文中截断半径为5埃,且是在单相聚集区内生成径向分布函数的,如第一个图所示,我就是受到这个划分单相聚集区的方法的启发,来生成RDF的,但结果...);此外,论文中我并没有发现RDF的公式,RDF是用ovito后处理得到的,我想应该也用不到公式吧。
      看起来还是我的做法有误,虽然说模拟数据不同,但绝不会相差很大,再说都是采用球形压头进行单晶硅的纳米压痕加工,应呈现出相似的规律才是。前辈,您是怎么生成RDF,又是怎么进行配位数分析的呢?好困惑啊
      期待您的回复!
MD simulation之后处理ovito中遇到的问题
5.jpg


MD simulation之后处理ovito中遇到的问题-1
6.jpg


MD simulation之后处理ovito中遇到的问题-2
7.jpg

人在塔在
7楼2021-07-29 22:02:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

get-it

新虫 (小有名气)

你用的是说明modification,量过原子间距么

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

2楼2021-07-27 19:12:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

HJieY

新虫 (初入文坛)

送红花一朵
引用回帖:
2楼: Originally posted by get-it at 2021-07-27 19:12:15
你用的是说明modification,量过原子间距么

前辈,
      你好!
      我量过原子间距,未成键的原子基本都是由于原子间距大于截断半径;我在生成键之前,将显示效果设置为只显示5配位原子,哦,怪不得,原来是这样,5配位原子不一定非要和5配位原子才能成键。
      前辈我提的第二个问题其实是针对SiC的径向分布函数来说的,RDF中的横坐标(原子间距)应该是哪种键的长度呢?是Si-C,Si-Si还是C-C呢?RDF图我是这么做的,首先进行配位数分析(根据键长设置为2.5埃),根据配位数着色(配位数范围3-7),然后通过表达式选取单个配位数(如配位数为5)的原子群,再次进行配位数分析,设置截断半径为5埃,得到RDF;但是我所生成的单配位数径向分布函数峰值与我所读论文中的相比,直接打了3个数量级,我有点犯傻了,难道是我生成RDF的方法不正确吗?
      十分期待您的回复!
                                                                                                                                                                                                                                                        迷茫虫友
人在塔在
3楼2021-07-28 09:02:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

get-it

新虫 (小有名气)

为什么要单独选择配位数为5的原子画rdf?文献也是这个流程吗?你自己都没搞懂问题在哪,你得把文献看明白了。
4楼2021-07-28 17:10:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
信息提示
请填处理意见