lammps案例：堆垛层错能的计算

最新推荐文章于 2024-07-07 23:44:50 发布
最新推荐文章于 2024-07-07 23:44:50 发布
阅读量4k 收藏 18
点赞数 2
分类专栏： lammps教程 文章标签： lammps 分子动力学模拟
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/lammps_jiayou/article/details/119241263
版权

大家好，我是小马老师。
本文介绍如何使用lammps计算晶体的堆垛层错能。
正常晶体中的原子按照一定的规律排列，原子面按ABCABC……的顺序排列，如下图所示：
在这里插入图片描述

当某个原子面被抽出后会形成堆垛层错，如图中的蓝色原子层被抽出，原子面排列变为ABCBCA……：
在这里插入图片描述

可通过计算堆垛层错能形成前后的能量差计算堆垛层错能，计算公式为：
在这里插入图片描述

E0为初始能量
Ef为删除一个原子面后的能量
A为界面面积
下面分享一个lammps计算堆垛层错能的in代码，本代码由Richard Glaze完成，本文稍作修改并进行了注释，以供初学者参考，对原作者表示感谢。

# Richard Glaze,
了解本专栏 订阅专栏 解锁全文 超级会员免费看
lammps加油站
关注
  • 2
    点赞
  • 踩
  • 18
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
订阅专栏
defect_lammps_lammps计算缺陷形成能_缺陷能_
10-02
本程序为LAMMPS程序，用于计算缺陷形成能
lammps实例
12-08
关于一些lammps资源的实例，可以对lammps有一个了解。
lammps建模案例分享及详解
最新发布
qq_44963682的博客
07-07 455
本文介绍了如何使用lammps提供的命令建立简单体系的模型，分享了一些lammps的建模案例。如遇到较复杂的模型，可选用Material Stdiuo等建模软件。
3_层错与孪晶_lammps代码_孪晶_lammps_lammps层错能_计算孪晶能
09-10
lammps代码，用于计算层错能和孪晶形成能
lammps软件_LAMMPS学习系列（11）
weixin_39612540的博客
12-03 1772
实例一 在用LAMMPS进行分子动力学模拟之前，最重要的一件事就是找到合适的相互作用势，那么如何判断所得到的经验势是否合适呢？今天就结合一个简单的实例介绍一下如何去判断：直观判断：做个简单的动态弛豫模拟，看看体系是否稳定，结构是否变混乱，具体弛豫模拟可见我之前的推文：LAMMPS学习系列（1）；能量判断（推荐）：计算相关的缺陷（空位，间隙子。。）形成能，层错能，迁移能等能量值或者弹性常数等力...
lammps的具体实例
05-11
lammps的具体实例
lammps经典例子
04-04
lammps经典例子
lammps代码实例分享
qq_44963682的博客
07-19 3710
分享在分子动力学模拟中常用的网站，以及简单的lammps代码实例。
3_层错与孪晶,孪晶与层错能的关系,Unix_Linux
09-10
本主题聚焦于"3_层错与孪晶"，以及它们与"孪晶与层错能"的关系，并涉及到使用LAMMPS来计算这些能量。以下将详细阐述这些概念及其关联。 层错与孪晶是晶体结构中的两种重要缺陷。层错，也称为面缺陷，是晶体中相邻...
LAMMPS 19个典型案例包括in文件后处理
09-23
LAMMPS 的的19个典型案例，包括in文件、分析原理和matlab后处理等，适用于新手的入门教程。
LAMMPS教程
08-14
介绍了分子动力学模拟的基础知识，lammps的指令和实例。
W、CO对Ni合金层错能影响的热力学计算 (2004年)
05-12
运用置换原子计算层错能的热力学模型，计算了Ni-W、Ni-CO合金的层错能．结果表明：随温度升高合金的层错能提高．随着W、CO含量的提高，Ni-W、Ni-CO合金的层错能提高；随W、CO质量分数的增多，Ni合金层错能先增加后减小，存在极大值，但W使Ni合金层错能增加的幅度较小．
基于LAMMPS分子动力学MSD自编程计算程序
03-21
**基于LAMMPS分子动力学MSD自编程计算程序** LAMMPS（Large-scale Atomic/Molecular Massively Parallel Simulator）是一种广泛使用的开源分子模拟软件，适用于研究各种尺度的物质系统，包括原子、分子以及生物大...
lammps案例：扭转
lammps_jiayou的博客
01-10 780
扭转有两种方式，一是固定一端旋转另外一端，二是两端沿相反的方向旋转。案例为两侧旋转方式，扭转后的效果如图所示。
lammps案例分享：Al/聚氨酯界面模拟
lammps_jiayou的博客
11-23 2857
大家好，我是小马老师。 本文分享一个Al和聚氨酯界面的模拟案例。 案例由Nuwan Dewapriya编写，原文可到以下网站下载： https://github.com/nuwan-d/LAMMPS_tutorials_for_short_courses/tree/main/polyurethane_aluminum_interface 本文对案例做一些注解，以供初学者参考。 金属和聚合物界面的研究相对较多，但是对于没有经验的lammps初学者来说，可能还是有些困难。 （1）建模问题 聚合物原子之间有键
LAMMPS案例分析——水滴的蒸发
E_Magic_Cube的博客
06-27 3312
LAMMPS案例分析——水滴的蒸发
lammps案例分享：流体通过多孔膜渗透模拟
qq_44963682的博客
10-22 2392
分享一个简单的lammps模拟案例，模拟流体在刚性多孔膜的渗透流动，本文只模拟了50%膜对流体流动的影响，后续可以改变膜原子的比例来观察其对流体流动的影响。
写一个lammps代码计算112方向的广义层错能
05-12
LAMMPS是一个分子动力学模拟软件，用于模拟原子和分子的相互作用。要计算112方向的广义层错能，可以使用LAMMPS中的dislocation工具来模拟。 以下是一个简单的LAMMPS输入文件： ``` # LAMMPS input file for calculating generalized stacking fault energy # along the 112 direction # define units units metal # define atom style atom_style atomic # define simulation box lattice fcc 3.52 region box block 0 20 0 20 0 20 create_box 1 box create_atoms 1 box # define potentials pair_style eam/alloy pair_coeff * * Ni_u3.eam.alloy Ni # define settings for dislocation simulation lattice orient x 1 -1 0 orient y 1 1 1 orient z -1 -1 2 dislocation 1 line 0.5 0.5 0.0 0.0 loop 20 20 20 start 0 0 0 change_box triclinic displace_atoms all move 0.0 0.0 20.0 units box # define settings for GSF calculation variable dx equal 0.01 variable dy equal 0.01 variable dz equal 0.0 # loop over dx and dy values label loop variable x loop 0 1 ${dx} units box variable y loop 0 1 ${dy} units box displace_atoms all move {x} {y} {dz} units box minimize 1.0e-10 1.0e-10 1000 10000 compute gsf all dislocation/stack/energy 1 1 3.0 variable energy equal c_gsf[1] variable x delete variable y delete print "{dx} {dy} {energy}" jump in.gsf loop # output final results thermo_style custom step pe run 0 ``` 这个输入文件首先定义了模拟所需的单位、原子样式和盒子大小。然后定义了用于模拟的势能函数，并设置了LAMMPS的dislocation工具来模拟位错。在这里，我们使用一个线性位错，其位错线方向为[0.5 0.5 0.0]，沿z轴方向移动20个原子距离。 接下来，我们使用一个循环来遍历dx和dy值，并使用displace_atoms命令在每个循环迭代中将晶格稍微扰动一下。然后我们使用minimize命令进行能量最小化，并使用compute命令计算广义层错能。最后，我们输出广义层错能的值。 要运行这个LAMMPS输入文件，您需要将其保存为一个名为“in.gsf”的文件，并使用以下命令运行LAMMPS： ``` lmp_serial < in.gsf ``` 这将生成一个输出文件，其中包含每个dx和dy值的广义层错能。您可以使用任何适合您的工具来分析这些数据并计算112方向的广义层错能。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

lammps加油站

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值