24小时热门版块排行榜    

查看: 1946  |  回复: 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的回帖

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的回帖

HJieY

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by get-it at 2021-07-28 17:10:42
为什么要单独选择配位数为5的原子画rdf?文献也是这个流程吗?你自己都没搞懂问题在哪,你得把文献看明白了。

前辈,
      您好:
      我看过的论文中确实展示出了单个配位数的RDF和ADF,但是没有提示是怎么做出来的,其实我自己做的径向分布函数是自己想出来的。
      前两个图是我读到的论文中的插图,最后一个是我自己按自己的想法做出的图,但是,与论文中的RDF的峰值相差了好几个数量级,不知是怎么回事?
      附件中的图片分别就是4配位原子和6配位原子的RDF和ADF,而Si-I、Si-II、、Si-III、Si-XII、Si-XIII是单晶Si在切削过程中所演变出的不同晶相结构,图中虚线分别是各种标准晶相所对应的RDF位置。
      十分期待您的回复!
MD simulation之后处理ovito中遇到的问题
1.jpg


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


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

人在塔在
5楼2021-07-29 09:11:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

get-it

新虫 (小有名气)

文献的rdf很奇怪。对于周期性体系,截断半径趋于无穷时,g(r)趋于1(你可以找些典型的rdf看看)。但是这里峰值都没到0.5,你可能要看看文献中rdf的公式才行。

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

6楼2021-07-29 19:46:57
已阅   回复此楼   关注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中遇到的问题-3
5.jpg


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


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

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

get-it

新虫 (小有名气)

看了文献我大概理解了,文献算的是有限的体系。你用ovito分析时可能需要把周期性关掉。
至于纳米压痕,我也没做过,建议你直接写邮件给作者询问。
8楼2021-07-29 23:20:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

HJieY

新虫 (初入文坛)

引用回帖:
8楼: Originally posted by get-it at 2021-07-29 23:20:40
看了文献我大概理解了,文献算的是有限的体系。你用ovito分析时可能需要把周期性关掉。
至于纳米压痕,我也没做过,建议你直接写邮件给作者询问。

前辈,谢谢您的建议!  哎...
人在塔在
9楼2021-07-30 09:53:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 HJieY 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] mitsunobu反应求助 15+3 cuicuicb 2024-06-03 4/200 2024-06-03 17:28 by w87437985
[基金申请] 为什么美国要封锁华为而不封锁丘C同呢? +14 zju2000 2024-06-02 18/900 2024-06-03 16:57 by ducan21
[教师之家] 统计一下:硕士毕业答辩后的谢师宴是学生出钱,还是老师出钱? +23 苏东坡二世 2024-06-02 28/1400 2024-06-03 16:40 by jurkat.1640
[基金申请] 前些天开会有个人见到人就搞关系,一查此人全是MDPI/Hindawi论文,鄙视! +19 zju2000 2024-06-02 20/1000 2024-06-03 16:25 by lkhjiayou
[基金申请] 入职高校3年发表10+SCI,尽人事听天命 +34 kaoyan250 2024-05-27 52/2600 2024-06-03 13:49 by kaoyan250
[论文投稿] 投稿的时候编辑发来邮件说作者机构信息没填写怎么办 5+3 飞飞0421 2024-06-02 7/350 2024-06-03 13:39 by topedit
[硕博家园] 材料博士毕业 +3 112233ssg 2024-06-02 3/150 2024-06-03 09:59 by Lwyy1106
[论文投稿] Scientific reports 投稿 +4 chencome12 2024-06-02 6/300 2024-06-03 09:13 by chencome12
[电化学] 电化学传感器 +3 哈哈嘿哟 2024-06-01 5/250 2024-06-03 08:46 by mecart2021
[论文投稿] 求助大神,Fe和Al离子对MOF都有淬灭,当两种离子共存时,怎么区分两种离子? 10+3 maoxiao 2024-06-02 3/150 2024-06-03 01:43 by nono2009
[基金申请] +5 河马の史诗 2024-06-02 5/250 2024-06-02 19:05 by 328838485
[考博] 申请2024或2025年博士研究生 +5 嘟噜嘟1 2024-05-29 12/600 2024-06-01 22:36 by 嘟噜嘟1
[考博] 24or25材料专业申博 +4 农夫三拳有点痛 2024-05-30 11/550 2024-06-01 14:45 by Napoleonsky
[基金申请] 博后特别资助状态变化 +24 随梦而飞2017 2024-05-30 35/1750 2024-06-01 10:10 by 青岛阳仔
[基金申请] 数理的人才答辩通知发了吗 +7 yzy3327 2024-05-30 7/350 2024-06-01 09:15 by 撒哈拉西部
[硕博家园] 哈工大硕博招生!博士每月入学! +4 nailooo 2024-05-30 5/250 2024-06-01 06:47 by anevay
[教师之家] 中年 (金币+3) +18 459582015 2024-05-28 19/950 2024-06-01 00:41 by 沈婉婷.Girl
[论文投稿] 审稿专家比较坚定的让补充实验,但实在没法补充实验,修回还有希望吗? (EPI+1) 3+3 qweasd12345 2024-05-29 6/300 2024-05-30 08:11 by qweasd12345
[论文投稿] 高手朋友推荐比较容易投稿和录用的SCI期刊,不投稿收费SC,对分区没有要求 5+3 xintangren 2024-05-28 4/200 2024-05-29 10:46 by xintangren
[基金申请] E10开始送了,希望有好运 +5 sail 2024-05-27 5/250 2024-05-28 18:36 by 芝小芝
信息提示
请填处理意见