| 查看: 1966 | 回复: 8 | |||
| 【悬赏金币】回答本帖问题，作者HJieY将赠送您 15 个金币 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
MD simulation之后处理ovito中遇到的问题
|
|||
|
前辈们，朋友们： 你们好！ 我的研究方向是采用分子动力学模拟硬脆性材料在纳米加工中的材料去除机理及变形行为。最近在采用ovito尝试对模拟结果进行后处理时，出现了一些问题不知如何解决，特别希望热心的虫友可以为我释惑。问题如下： 1、在ovito中进行配位数分析，在视觉效果上感觉配位数设置无误后（截断半径我设置为稍大于键长）， 进行成键时却发现配位数与中心原子的键数十分不一致，如5配位原子的中心原子键数应为5才对，但键数却表现的极为任意（0~5个），原子间没有键居多，成5键的原子几乎没有，真的很奇怪； 2、在SiC材料中，在Si原子和C原子成键时，C原子和C原子、Si原子和Si原子之间是否会成键呢？ 先行感谢致力于研究的你们，盼解惑！ 一个在分子动力学上处于迷茫的人儿 |
回复此楼» 猜你喜欢
我标书代码变了
已经有74人回复
-
试用期辞职
已经有7人回复
-
豫北虫友互识
已经有5人回复
-
F口信息学部拿面上，大概需要什么样的成果
已经有13人回复
-
累并快乐着
已经有16人回复
-
关于硕博连读的一些疑问？
已经有7人回复
-
2024国社科通讯评审
已经有30人回复
-
F口401需要啥文章水平
已经有7人回复
-
去英国的小伙伴儿都在哪儿租的房子呀？
已经有9人回复
-
跑板能跑开，过柱过不纯怎么办
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
送红花一朵|
前辈， 您好： 受到你的回复，我很高兴！我也感觉RDF挺奇怪的，但不只这一篇论文是这样，我看到的另一篇论文也是这样的，RDF峰值最大为1左右（如第二、三图所示），但就像您说的那样，截断半径不可能无穷大（此篇论文中截断半径为5埃，且是在单相聚集区内生成径向分布函数的，如第一个图所示，我就是受到这个划分单相聚集区的方法的启发，来生成RDF的，但结果...）；此外，论文中我并没有发现RDF的公式，RDF是用ovito后处理得到的，我想应该也用不到公式吧。 看起来还是我的做法有误，虽然说模拟数据不同，但绝不会相差很大，再说都是采用球形压头进行单晶硅的纳米压痕加工，应呈现出相似的规律才是。前辈，您是怎么生成RDF，又是怎么进行配位数分析的呢？好困惑啊 期待您的回复！  5.jpg  6.jpg  7.jpg |
7楼2021-07-29 22:02:46
get-it
新虫 (小有名气)
- 应助: 23 (小学生)
- 金币: 1825.8
- 散金: 26
- 红花: 4
- 帖子: 273
- 在线: 448.3小时
- 虫号: 3212202
- 注册: 2014-05-16
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2021-07-27 19:12:15
送红花一朵|
前辈， 你好！ 我量过原子间距，未成键的原子基本都是由于原子间距大于截断半径；我在生成键之前，将显示效果设置为只显示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
get-it
新虫 (小有名气)
- 应助: 23 (小学生)
- 金币: 1825.8
- 散金: 26
- 红花: 4
- 帖子: 273
- 在线: 448.3小时
- 虫号: 3212202
- 注册: 2014-05-16
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2021-07-28 17:10:42
50