atomsk建模——任意取向模型的注意点

atomsk建模需要注意晶格对称性

Atomsk 会利用 fcc和bcc 晶格的对称性来减小单元的尺寸

在 fcc 晶格中,当一个米勒指数为奇数而另外两个为偶数时,沿给定 [hkl] 方向的最短晶格向量实际上是 1/2[hkl]。

在 bcc 晶格中,当所有三个米勒指数都是奇数时,沿 [hkl] 方向的周期减少到 1/2[hkl]。

atomsk官网原文:

​​​​​​​Atomsk - 教程 - 晶体取向 (univ-lille.fr)

实例

例如用Atomsk建立[1-10] [11-2] [111]取向的NiTi模型,

atomsk --create bcc 3.01 Ti Ni orient [1-10] [11-2] [111] -wrap NiTi.xsf

​​​​​​​[111]方向得到的尺寸并不是预期的3√3=5.2,而是2.6,得到的模型中Ti:Ni=4:2,明显不是我们想要的

但我们可以根据晶格的对称性推导出最小的晶胞

取巧的,可以用以下命令大致观察晶胞的原子分布

atomsk --create bcc 3.01 Ti Ni orient [1-10] [11-2] [111] -wrap Ti.xsf

atomsk --create bcc 3.01 Ni Ti orient [1-10] [11-2] [111] -wrap Ni.xsf

atomsk --merge z 2 Ti.xsf Ni.xsf  NiTi.xsf lmp

需要注意的是,得到的原子坐标不能马上使用,因为使用merge命令的时候原子会比实际上靠的更近,所以此步骤目的是为了观察晶胞的原子分布,具体的原子坐标需要用晶格的对称性推导

在使用Atomsk进行建模替换时,可以使用命令"atomsk -remove-atom <index>"或者"atomsk -remove-atom <species>"来删除给定索引或者种类的原子。另外,也可以使用"atomsk -remove-atom select"命令来删除所选择的原子。\[1\] 例如,如果要替换原子生成合金结构,可以使用以下命令: ``` units metal dimension 3 boundary p p p atom_style atomic timestep 0.001 neighbor 0.2 bin read_data final.lmp set type 1 type/ratio 2 0.33 8793 set type 1 type/ratio 3 0.5 56332 #替换原子,1:1:1 write_data FeNiCr.data ``` 这个命令将会替换原子,生成合金结构,并将结果写入FeNiCr.data文件中。\[2\] 另外,如果想使用Atomsk建立管道类模型,可以使用命令"atomsk -roll <direction> <θ> <axis>"。这个命令可以根据指定的方向、角度和轴来建立管道类模型。\[2\] 希望以上信息对您有所帮助。如果还有其他问题,请随时提问。 #### 引用[.reference_title] - *1* *2* [lammps和atomsk中常用的建模方法](https://blog.csdn.net/qq_44963682/article/details/125883094)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] - *3* [Atomsk空位缺陷建模技巧](https://blog.csdn.net/lammps_jiayou/article/details/119208769)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v91^control_2,239^v3^insert_chatgpt"}} ] [.reference_item] [ .reference_list ]
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值