3.1.2 进展过程中将问题简化，暂定只研究单原子铜在铍多晶中的扩散问题2024.04.03

由于在计算刚开始并未使用kspacing等操作，且不同的计算体系能量出现整体偏移情况。因此后续把所有数据重新进行一遍相同设置的计算。有了初次的经验，再次计算时决定把问题简单化。本次计算主要针对缺陷处铜扩散较快问题做研究，且目标定为温度对晶体中的铜扩散影响研究，这样还能简化晶界等构型的影响。进而在后续的计算中，只计算铍单质，单原子铜掺杂和双原子铜掺杂情况，并且不考虑表面构型。

1. 考虑铍单晶的不同构性

单晶的训练集只计算hcp、fcc构型，这两个构型具有相似的能量。将cubic和tetragonal的单晶构型删去，减少初始训练集的准备量。

Be16_hcp
Be54_hcp
Be128_hcp
Be32_fcc
Be108_fcc

2. 考虑铜单原子和双原子掺杂构型

其中以hcp结构为主要构型,因为实验上观测到磁控溅射生长出来的薄膜儿基本呈现六角立方结构。

Be15_Cu1_hcp
Be53_Cu1_hcp
Be127_Cu1_hcp
Be53_Cu1_hcp_sa_plane
Be32_Cu1_fcc
Be47_Cu1_tetragonal
Be53_Cu1_cubic

Be52_Cu2_hcp
Be52_Cu2_hcp_sp_plane

Be31_Cu2_fcc
Be46_Cu2_tetragonal
Be52_Cu2_cubic

同时筛选的时候，基本以50多个原子构型为基础进行初始训练机的训练。即只对hcp结构进行了扩包。其他构型可能不是计算的主要内容，因此初始训练集只是部分涉及相关构型，且把掺杂构型跟单晶构型合并到了一起。

当然到目前为止并未针对特定的问题进行数据集的准备，例如我后续想研究晶界对Cu扩散的影响，那么我可以把晶界构型加入到训练集中，进行势函数的扩充计算。

作者：朱雪刚 邮箱:xuegangzhu@qq.com； 工作单位：石家庄学院 理学院/北京科学智能研究院（AISI）访问学者2023.07-2024.09，访问导师北京大学陈默涵; 徐张满仓 邮箱: xuzhangmancang@dp.tech

截止2024.07.21录制视频教程已上传至Bohrium的课程《DeePMD应用案例讲解：铜原子掺杂铍晶体的机器学习势函数拟合过程演示》网址： https://bohrium.dp.tech/courses/1075495070?tab=courses 后续会在Bohrium平台更新

注意：后续的更新，大部分会在Bohrium课程平台进行，请看教程入门的同学加入课程进行学习；且录制的视频课程会把一些个人观点给加入，从个人观点来看比文字教程的内容更多

写教程内容讨论QQ群:143276924 DPGEN+ABACUS教程准备；

ABACUS软件的QQ群：759914681，群内有专职开发人员，目前ABACUS提问问题以github的issue为主，群内可作为辅助提问。

欢迎大家推广本教程，让更多的dp入门学习者有个参考，目前2024.09.16我联系的微信公众号推广是 lammps加油站的小马老师（我也报名了小马老师的一对一辅导）。当然也看到lammps爱好者在转发推广。感谢大家公众号的推广。