手机版
           

MS计算异质结:从晶格匹配到能带对齐的CASTEP实战路径

发布时间:2026-06-22   来源:科研学术网    
字号:

MS计算异质结——在Materials Studio环境下用CASTEP模块做半导体异质界面模拟——是光电材料、催化异质结构和二维材料器件研究的标准工具链之一。但界面建模阶段的一个晶格匹配错误,传到后面能带对齐环节可以产生0.2-0.3 eV的人为偏移,这个量级足以改变界面能带排布的类型判定(type-I vs type-II)。

Build Layers:晶格匹配是费时的,但跳不过去

MS的Build Layers工具是界面建模的入口。选择两个表面模型→指定劈裂面指数→输入晶格失配容忍度→生成界面结构。问题在于,大多数半导体异质结对(如ZnO/GaN、TiO₂/SrTiO₃)的晶格常数天然不匹配——ZnO a=3.25 Å,GaN a=3.19 Å,失配度~1.9%——不算大但也不能直接忽略。

处理失配的策略只有三条路:压缩张力侧、拉伸受压侧、或者做超胞匹配。压缩张力侧(即让晶格常数较小的材料拉伸匹配较大的)会引入面内拉伸应变,改变电子结构——ZnO在2%的双轴拉伸下带隙收缩约0.1 eV(PBE)。拉伸受压侧则反之。哪种方案更好取决于实际异质结的应变分布——如果ZnO外延生长在GaN衬底上,ZnO会承受拉伸应变,那么在计算中拉伸ZnO匹配GaN的晶格常数更接近真实物理。

超胞匹配是精度最高的方案——如用4×4的ZnO (a=4×3.25=13.00 Å)匹配5×5的GaN (5×3.19=15.95 Å),虽然面积匹配仍有~18%的失配,但在一个方向上可以做到完全共格。这种方案的代价是超胞原子数暴涨——240+原子——CASTEP的SCF收敛速度显著下降,对计算资源有硬要求。

应变分配与界面间距优化

晶格匹配完成后,界面间距(两个表面之间的真空区厚度)需要扫描优化。在MS中,固定面内晶格→沿界面法向扫描两层间距→在每个间距点计算总能,总能最低的间距就是最优界面间距。

这个过程不能跳过。以ZnO/GaN界面为例,初始间距设为层间键长的平均值(2.0 Å),扫描间距1.5-3.0 Å(步长0.1 Å),总能曲线呈典型的Morse-like形态——最小值出现在2.1 Å附近,向两侧偏离后总能上升。如果跳过扫描直接用默认间距,0.3 Å的间距偏差对应总能高估~0.5 eV——这个误差在后续能带对齐中会乘法式放大。

另一个容易被忽视的参数是偶极修正。在非对称异质结(如polar ZnO/polar GaN)中,界面的自发极化会在真空区建立一个人工电场——这个电场在周期性边界条件下无法通过标准DFT自洽消除,必须加偶极修正(dipole correction)。CASTEP的偶极修正在Electronic选项卡的SCF设置中打开,在polar/polar界面中这个开关的ON/OFF差距可达0.15-0.2 eV。

能带对齐:静电势参照与带边能量

能带对齐是异质结计算的核心产出——两张能带图排在一起没有意义,必须共用能量零点。CASTEP的静电势参照法是通过平面平均静电势的分布,在两个材料各自的体相区域提取平均值,差值即为带边对齐所需的能量偏移。

实际操作:先分别对两个体相材料做SCF→提取平面平均静电势在真空/体相区域的平台值→计算VBM相对于静电势的位置→对异质结界面模型,从界面区域的静电势确定共同的能量零点。这个流程的精度瓶颈在于体相计算的k点网格必须足够密,保证提取的VBM能量误差<0.02 eV——4×4×1的k点通常不够,6×6×1或更密才能将统计噪声压低。

在ZnO/GaN异质结的项目中,CASTEP算出的价带偏移(VBO)约1.3 eV,导带偏移约0.6 eV,形成type-II能带排布。这些数字跟实验光电子能谱(XPS)的测量偏差~0.15 eV——考虑到DFT-PBE的带隙本身被低估约35-50%,这个VBO的一致性已经相当好。因为价带偏移的计算更多依赖静电势和价带边的相对位置,对带隙低估的敏感度远低于导带偏移。

应变带来的悄然变化

界面应变对带隙的影响在异质结中不容忽视。ZnO在1.9%双轴压缩下带隙从0.76 eV(PBE裸值)升至0.85 eV——带隙变化0.09 eV虽不算大,但在能带对齐的语境中足以让导带偏移移动0.05 eV,改变界面态的充放电行为。

应变效应的另一个维度是改变了带边的轨道成分。ZnO的价带顶主要来自O-2p轨道,在压缩应变下O-O间距缩小,p-p排斥增强,价带展宽——这会改变VBM的精确位置,进而影响VBO的数值。MS环境下处理应变效应的优势在于几何建模的直观性——应变可以在Build Layers时直接指定,比手动改POSCAR坐标高效且不易出错。

MS计算异质结的最终价值体现在一条物理链上:晶格匹配→界面间距优化→静电势对齐→能带偏移→器件性能预测。其中每一环的误差传递都需要量化判断,而MS平台的前处理效率让这条链的参数扫描和误差评估成为实际可行的操作。

更多内容请访问 https://www.keyanxueshu.com/

图说天下

×
cp2k计算
dft计算
Gaussian计算
MS计算
VASP计算