Q1：朱老师，想请教一个问题，对层状结构ZnIn2S4进行结构优化，优化前上层Zn和下层S 键长为3.86埃，优化完后变为2.61D3。另外同一层的Zn-S键长由2.456埃变为3.458埃了。我用castep和vasp试了都是这样的变化，请问一下这种情况该怎么办呢？Castep中我选了Grimme修正。vasp中用了DFT-D3。
A：
①3.458 Å基本已不是正常Zn–S键，说明Zn从原位跑到了别处；
②CASTEP和VASP都这样，说明更像模型本身不稳定；
③先检查晶型、原子占位、层间堆垛是否建对；
④建议先优化体相晶胞验证结构，再做分步优化。
Q2：朱老师，计算hse能带的istart我设置为1然后icharge为11时候发现了报错？
A：先PBE 优化/自洽，保留 WAVECAR，然后做 HSE 时设 LHFCALC=.TRUE.，ISTART=1，ICHARG 不要设 11，用普通自洽方式跑，若要高对称路径能带，就用 0-weight KPOINTS 方法。
一、HSE 不能配 ICHARG=11 用，报错是正常的；
二、HSE 能带正确做法一般有两种：
①直接做 HSE 自洽；
②先做普通网格点的自洽，再在 KPOINTS 里把能带路径点作为 0-weight 点 加进去，让 HSE 继续自洽输出能带。
A：常建一个晶胞里放要研究的那一组分子，最常见就是两个，不要随便放多个：若研究两分子相互作用，就在一个大真空盒子里放这两个分子；
真空要足够大，避免与周期镜像相互作用，一般各方向至少15–20 Å；
放多个分子只适用于研究团簇、浓度效应或协同作用，不是常规做法。
Q4：朱老师，这个VASP报错是怎么回事呢？我该怎么调整？
1、先查结构，是否有原子过近、初始键长不合理;
2、IBRION=3 + ISIF=3 偏激进，先改成 IBRION=2，加 POTIM=0.2~0.3，必要时先用 ISIF=2 只放原子不放晶胞；
3、保留 ALGO=Normal，把 NELM 提到 120，必要时增大 NBANDS，若还不稳，再去换 ALGO、调混合参数。
Q5：朱老师，构建缺陷以后，在 没有缺陷的带隙内多出来的就是缺陷态吗？
A：不一定，但很多时候是，还要区分真正缺陷态和带边扰动态
若缺陷引入的能级落在本征带隙内，且电荷明显局域在缺陷附近，通常可认作缺陷态；
但有些缺陷态不一定在带隙里，也可能与价带或导带混合，所以“带隙内多出来”是强信号，不是唯一标准；
超胞缺陷计算会带来能带折叠，图上多出来的线不一定全是新缺陷态，最好做 band unfolding 辅助判断；
最稳妥的判据是一起看：PDOS、分波电荷密度/partial charge、局域性、自旋分裂，再和无缺陷体系对比。
