计算化学公社

 找回密码 Forget password
 注册 Register
Views: 1833|回复 Reply: 5
打印 Print 上一主题 Last thread 下一主题 Next thread

[Lammps] 求助:ovito分析多次级联碰撞产生的缺陷异常增多的问题

[复制链接 Copy URL]

25

帖子

0

威望

131

eV
积分
156

Level 3 能力者

本帖最后由 墨茶茶 于 2024-2-14 00:19 编辑

各位老师好!我最近在用lammps做多次级联碰撞的研究。我用lammps建立了一个模型,首先使用displace_atoms将盒子中所有的原子向任意方向平移一个矢量,之后在模型中随机选择一个中间位置的原子作为PKA原子进行级联碰撞模拟(每次级联碰撞选择的原子种类和赋予的能量都相同),级联碰撞结束后将盒子中的原子整体反向平移相同的矢量,输出模型,之后再重复上述步骤进行下一次级联碰撞。之后我将每次级联碰撞结束后的模型放入ovtio中用ovito自带的Wigner-Seitz defect analysis功能进行点缺陷分析。分析发现,点缺陷的数量会在某些时刻发生极其异常的暴增,之后又会很快的回到暴增前的水平。比如第n次级联后模型中残余的缺陷为10^2个左右,但到第n+1次级联碰撞后模型中的点缺陷又会暴增10^5个左右(就像凭空出现一般),并且到第n+3次左右又会回到10^2个左右,之后会反复出现这种情况。我翻遍国内外论坛,有人认为是ovito自带的软件识别错误,有人认为是模型发生了漂移(我只在级联碰撞前后使用了displace_atoms命令),但没有人提出解决方案。因此我想请问一下各位老师有没有碰到和我相同的情况,请问是如何解决的。

第n次级联碰撞时的缺陷数量

第n+1次时级联碰撞缺陷数量

第n+2次时级联碰撞缺陷数量

2

帖子

0

威望

19

eV
积分
21

Level 1 能力者

2#
发表于 Post on 2024-5-18 10:26:37 | 只看该作者 Only view this author
您好,我也是做辐照相关的,最近遇到类似的问题,可以交流下吗?

2

帖子

0

威望

19

eV
积分
21

Level 1 能力者

3#
发表于 Post on 2024-5-18 10:28:17 | 只看该作者 Only view this author
小石 发表于 2024-5-18 10:26
您好,我也是做辐照相关的,最近遇到类似的问题,可以交流下吗?

联系方式qq:2622442820

25

帖子

0

威望

131

eV
积分
156

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2024-10-7 16:59:19 | 只看该作者 Only view this author
小石 发表于 2024-5-18 10:26
您好,我也是做辐照相关的,最近遇到类似的问题,可以交流下吗?

你好,原因是在级联模拟过程中盒子中的原子整体发生了平移,因此在级联过程中要加一个fix   4 all recenter INIT INIT INIT,以固定盒子内的原子不要发生除碰撞外的平移。

1

帖子

0

威望

21

eV
积分
22

Level 1 能力者

5#
发表于 Post on 2024-12-30 23:17:51 | 只看该作者 Only view this author
墨茶茶 发表于 2024-10-7 16:59
你好,原因是在级联模拟过程中盒子中的原子整体发生了平移,因此在级联过程中要加一个fix   4 all recent ...

你好老师,我用了您的这个命令后还是会出现这种异常现象,是什么原因呢,我CNA和位错线分析都是没问题的,但就是在看frankelpair这里出了问题

25

帖子

0

威望

131

eV
积分
156

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2025-1-6 16:01:16 | 只看该作者 Only view this author
Haozhou 发表于 2024-12-30 23:17
你好老师,我用了您的这个命令后还是会出现这种异常现象,是什么原因呢,我CNA和位错线分析都是没问题的 ...

我也不太清楚了,你可以附上几个图看看

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-22 05:43 , Processed in 0.207609 second(s), 29 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list