作者： Zhang Zhaofu 以下是我阅读文献、结合自己的经验做的一些关于表面界面建模中的要点 / 经验小结，如果有任何不对，欢迎及时联系我做修改。希望转载注明源出处： http://blog.sciencenet.cn/blog-2686986-1175883.html 在涉及到表面或者界面计算的时候，要做的第一件事便是表面或者界面建模工作。因为在网上已经有很多建模经验总结帖子，所以对一些已经广泛流传的经验（比如如何建模根号表面或者界面）就不详细介绍，而是会更系统概括性列举出需要考虑的要点并列举一些经验。 首先要对涉及到的表面和界面的类型进行分析归类：对表面来说，主要包括 a). 2d 材料的表面模型 ; b) 3d 材料的表面模型。对界面来说，主要包括 a). 2d/2d 材料的界面模型； b). 2d/3d 材料的界面模型； c). 3d/3d 材料的界面模型。（注意， 3d 材料模型是指切面前材料就是 3d ，切面并在 c 方向加上真空组成 surface slab 或者 interface slab 后，在概念上是属于 2d ，即只有 ab 方向是无限延伸） 2d 材料的表面界面模型比较简单直接，放到最后面来讲。因为界面是由两种材料的 surface slab 构成，所以可以把表面建模中需要考虑的因素并规到界面建模需要考虑的因素之内。 对于 3d 材料的界面模型 ，首先要分两大类： A. 有真空，单一界面模型 ； B. 无真空，双界面模型 。前者模型从上往下分别为：真空层/材料A/材料B，视情况在B的底端有额外的钝化原子层，即在周期性边界条件下只有一个 材料A/材料B的界面 ；后者从上往下分别为：材料A/材料B，或者 材料A/材料B/材料A，无真空无钝化原子层，即在周期性边界条件下有两个 材料A/材料B的界面 。 对 A. 有真空单一界面模型 ，模型包括原子层和真空层。所以要考虑： 1. 考虑原子层的厚度 。厚度太薄不足以充分体现出材料的表面性质，厚度太厚是会造成额外的计算资源浪费。合理的模型厚度是指既包括体现表面材料 surface 性质的原子层，又包括体现体材料 bulk 性质的原子层。至于如何确定原子层厚度，可以考虑： a. Layer-projected DOS (LDOS)判断, b. 根据 averaged potential 即功函数判断 , c. 根据 quantum confinement 判断， d. 根据论文中取几层判断。这里结合我的论文 介绍下根据 LDOS 如何判断： Fig. 1 The local density of states of each layer in (a) GaN (001) and (b) HfO2 (001) surface model for layer thickness test. The vertical lines are given to identify the VBMs and CBMs. 以上图我发表的论文 HfO2/GaN 界面模型为例 介绍。建立一个很厚的 GaN 表面 slab ，然后直接算 DOS 即可。然后做出每一层的 LDOS 画在一起，会发现最表面几层原子是存在表面态，以及因为 quantum confinement 表面几层的 bandgap 会比较大。随着层数厚度加深会发现 surface state 消失且 VBM CBM 的能量位置趋于稳定。所以就可以以趋于稳定的层数来确定一个适合的层数厚度。 一些论文中也有介绍其他的判断方法，包括：功函数判据（建立不同层数的表面模型计算功函数，看看多少层功函数稳定）和 quantum confinement 判据（ PBE 会低估 Eg ，但是因为 quantum confinement 材料的 Eg 会变大，会存在某个层数使得算出来的 Eg 就等于或者接近实验值），这里不详细讲解了。 2. 考虑真空层的厚度。 由于主流的计算软件 (VASP CASTEP 等 ) 都是基于周期性边界条件 , 所以实际的计算模型在 c 方向是原子层 / 真空层 / 原子层 / 真空层 /… 无限重复模型。所以当真空层较薄时候，原子层会跟它的 image structure 镜像结构存在 interaction造成误差 ；但是真空层太厚又会导致计算量增大。测试真空层厚度常用的方法是 averaged potential 即功函数方法：做出 averaged potential 曲线，看真空区域是否是平的。 经验性的数据是：真空一定不能小于 10A ， 15-20A 为宜， 20-30A 更好。 3. 考虑终端钝化即 bulk termination 。 切出来的 slab 模型通常 ( 注 ) 在最上层和最下层存在悬挂键（ dangling bonds, DB ）。在表面或者界面处的悬挂键是真实存在的，对他们的不同处理可以反应表面界面的不同属性。但是最下层的悬挂键却不应该存在，因为这部分的原子层在模拟体相 bulk properties ，所以需要进行合理的钝化处理。最常见的方式是加氢（或者赝氢）钝化，具体方式为：计算出来悬挂键的电子数目 N ，然后选择 2-N 个电子的赝氢即可。以 III-V 为例，可以参考我的文章 中均有提及：对于 Ga ，每一个 DB 有 0.75e ，所以需要 1.25e 的赝氢钝化；对于 N ，每一个 DB 有 1.25e ，所以需要 0.75e 的赝氢钝化。 对于加赝氢的位置：通常选 DB 切开之前的键的中点位置加赝氢。 另外，对于离子型金属氧化物，除了加 H 钝化外，还可以选择均匀去掉半层 O 来钝化，可以参考 中对 HfO2 的处理。 注：对 non-polar surface 是不存在悬挂键，比如 cubic 材料的 (110) 面； 2d 材料的情况切出来单层也通常没有悬挂键。 4. 切出来晶格匹配 lattice matched 表面 slab 。 晶格匹配即包括 ab 晶格常数匹配，也包括 ab 夹角相等。对于 semiconductor/semiconductor or insulator/semiconductor ，建议 mismatch 不要超过 5% 。对于 metal/semiconductor or metal/insulator 来说，考虑到金属有较好的延展性，形变并不太会影响金属的性质（比如功函数）所以 mismatch tolerance 可以大一些，但是最好不要超过 10% 。 如何确定 matched 的界面呢，可以考虑如下的旋转方案：先对模型做很大的扩胞，只看一层原子，找出某个晶格参数下俩 surface slab 的晶格参数是 matched ，然后就可以根据根号建模方案切割表面即可。 cleave surface 的细节可以网上参考前人分享的根号建模经验帖，这里就不展开讲了。 Fig. 2 (a) Ga monolayer of the GaN (001) surface and (b) Hf monolayer of the HfO2 (001) surface. Red atoms are Ga, and purple atoms are Hf. 另外注意，在拼接界面的时候要确定界面的晶格常数，通常是要选择基底材料（下层材料）的晶格常数做为界面材料的晶格常数，即保持下层材料晶格常数不变，改变上层金属或者介质层的晶格常数，以此来模拟实验中材料生长的细节性质。 5. 初始构型的确定。 由于 界面模型 是人工直接 build layer 产生，界面成键细节不同会很大程度上影响系统总能；对于初始层间距的选择也是如此：合理的初始层间距可以很大的减小计算量。并且， DFT 在做结构优化的时候，有可能陷入 local energy minimum 而不是 global energy minimum 。为了尽量找到一个较为合理的初始构型，可以做 movement test. Fig. 3 The movement test of HfO2 surface in (a) z, (b) y and (c) x direction when GaN surface is fixed in the GaN/HfO2 interface model. 以上图我发表的论文 HfO2/GaN 界面模型为例 介绍。 Firstly, fix HfO2 (001) surface, move GaN (001) surface in z direction step by step, perform self-consistently calculation, and find the lowest energy position along z direction (see Figure S5(a)). Starting from the point with lowest energy, move GaN along x-(see Figure S5(c)) and y-(see Figure S5(b)) direction and locate the lowest energy before a full structure relaxation. 6. 固定优化。 这里包括两个固定：固定晶格常数，即不改变晶胞仅优化原子结构；通常选择固定最底层赝氢（以及靠近赝氢的一层或者几层 bulk layers) 来更好的模拟 bulk properties 。 对 B. 无真空，双界面模型 ，模型无真空，是由两层界面组成。要考虑： 1. 考虑原子层的厚度 。基本同 A1 。因为模型不包含真空，每一个界面两侧都是原子层，受 screening effect 的影响，原子层厚度其实可以比单界面模型稍薄一些。但是依然最好先做 layer thickness test 。 2. 不需要考虑真空层。 由于双界面模型是存在两个完全相同的界面，没有真空层，计算量也就可以略小一些。 3. 不需要考虑终端钝化即 bulk termination 。 由于双界面模型是存在两个完全相同的界面， 一样的 bonding characteristics, 所以不存在不合理的悬挂键，不需要考虑终端钝化 / 悬挂键处理。 4. 切出来晶格匹配 lattice matched 表面 slab 。 基本同 A4 。 5. 初始构型的确定。 基本同 A5 。 6. 固定优化。 无真空模型，但是因为存在 interfacial bonding and interfacial interaction 是会造成键长（层间距）改变，所以可以用固定 ab 方向、开放 c 方向晶格常数方案来做结构优化。并且固定基底材料的中间原子层即可。 通过以上对比 A. 有真空单一界面模型 , B. 无真空双界面模型 的建模要点，可以发现， B. 无真空双界面模型 计算量可以小一些，但是它要求必须是对称的模型，保证两个界面在建模时候完全相同。所以对于一些材料不适合这个建模方案。 注意：上述内容仅围绕表面界面建模的要点进行描述。对于表面材料，通常在表 面存在surface reconstruction ；对于界面材料，需要考虑 electron counting rule 。这些更复杂的考虑是要具体情况具体分析，这里不展开详细描述。关于界面考虑 electron counting rule的例子 ，可以参考我的论文 以及更多其他的论文。 对于 2d 材料的表面界面模型 ，相对来说比较简单，通常来说只需要考虑真空厚度 / 固定优化即可。 在 2d/2d 异质结中，两侧 的surface slab就具有很高的对称性。这样的对称性很高的异质结界面材料，如何确定一个合适的初始结构，以我MoS2/GaN 2d/3d 异质结 为例讲一下。 ( 虽然我的体系是 2d/3d 界面，但是方法是一样的 ) Fig. 4 relaxed MoS2/GaN heterostructures with various high-symmetry atomic arrangements. 界面两侧的材料本身就具有高对称性，且基本上是以 vdW 相互作用连接。可以：首先根据高对称性原则，得到几个不同的初始结构（上图的 6 个 type ）；然后对每一个 type 调整其层间距，算静态能量，找到他们的能量最低层间距，即为这个 type 的最佳初始构型；再从 6 个最佳初始构型中选择能量最低的来做系统研究即可。 以上便是我阅读文献、结合自己的经验做的一些关于表面界面建模中的要点 / 经验小结，如果有任何不对，欢迎及时联系我做修改。另外欢迎引用我的文章，即帖子中提及的 . 希望转载注明源出处： http://blog.sciencenet.cn/blog-2686986-1175883.html Zhaofu Zhang, et al, ACS Appl. Mater. Interfaces 2015, 7, 5141−5149 Zhaofu Zhang, et al, ACS Appl. Mater. Interfaces 2018, 10, 17419−1 7426 Zhaofu Zhang et al, Appl. Phys. Express 11 081003 (2018) Zhaofu Zhang et al, IEEE Trans. Electron Devices 2017 64 (10), 4036-4043 Zhaofu Zhang et al, Appl. Phys. Lett. 114, 161601 (2019)
标签: 界面建模
