Si能带计算讲义.docx

下载文档 降价啦

373
0
约2.45千字
约 8页
2017-04-02 发布于湖北
举报
版权申诉
保障服务

Si能带计算讲义.docx

1、本文档共8页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

Si的能带及光学特性计算梁自堃115113001184GGA计算能带system=Si ENCUT=300 ISTART=1 ICHARG = 11 PREC=Accurate GGA=PE ISMEAR=-5 EDIFF=1E-5 LOPTICS = TRUE CSHIFT = 0.1 NBANDS = 36 NPAR = 1 NEDOS = 2000 LRPA = .FALSE. LREAL = .FALSE. LWAVE = .FALSE. LCHARG = .FALSE. INCARSystem = fcc Si ICHARG=11 #charge read fileENCUT = 240ISMEAR = 0; SIGMA = 0.1;LORBIT=11KPOINTSk-points for bandstructure L-G-X-U K-G 10lineReciprocal0.50000 0.50000 0.50000 10.00000 0.00000 0.00000 10.00000 0.00000 0.00000 1 0.00000 0.50000 0. 50000 10.00000 0.50000 0.50000 10.25000 0.62500 0.62500 1 0.37500 0.7500 0.37500 10.00000 0.00000 0.00000 1POSCARfcc Si: 3.90.5 0.5 0.00.0 0.5 0.50.5 0.0 0.5 1cartesian0 0 0HSE计算能带首先，用 PBE 函数进行正常的自恰计算，这里需要输出波函数和电荷密度；其次，用 HSE 函数进行自恰计算，这里最好用 ALGO＝damp；第三，用 HSE 函数进行自恰计算，去除 ALGO 参数，这里使用 DAV 算法；最后，计算能带和态密度。有时可以将第三步和第四步合并。第一步，INCAR? SYSTEM=Si? ENCUT=500 ? EDIFF=1e-5 ? IBRION=-1 ? ISIF=2? NSW=0 ISMEAR=0；SIGMA=0.05 ? POTIM=0.1 ? #Wave function and charge ? LWAVE=T ? LCHARG=T第二步，? SYSTEM=Si? ENCUT=500 ? EDIFF=1e-5 ? IBRION=-1 ? ISIF=2? NSW=0ISMEAR=0；SIGMA=0.05POTIM=0.1#Wave function and chargeLWAVE=TLCHARG=TLHFCALG=.TRUE. #采用杂化泛函方法HFSCREEN=0.2ALGO=DampedTIME=0.4第三步，SYSTEM=Si? ENCUT=500 ? EDIFF=1e-5 ? IBRION=-1 ? ISIF=2? NSW=0 ? ISMEAR=0；SIGMA=0.05 ? POTIM=0.1 ? #Wave function and charge ? LWAVE=T ? LCHARG=TLHFCALG=.TRUE. #采用杂化泛函方法HFSCREEN=0.2 第四步，计算能带SYSTEM=Si? ENCUT=500 ? EDIFF=1e-5 ? IBRION=-1 ? ISIF=2? NSW=0 ? ISMEAR=0；SIGMA=0.05 ? POTIM=0.1 ? #Wave function and charge ? LWAVE=T ? LCHARG=TLHFCALG=.TRUE. #采用杂化泛函方法HFSCREEN=0.2 在第二步的基础上，使用 DAV 算法计算，好处是收敛较快，本征矢量没有任何问题。这一步的 DOS 是没有问题的。第四步计算能带和第三步没有差别。主要是 k 点的设置，由于 HSE 方法无法从电荷密度进行非自恰，所以只能用自恰的方法计算能带。通常是先把 IBZKPT拷贝成 KPOINTS。用GGA和HSE方法算出的Si能带光学特性HSE06INCARSYSTEM = Si-DiamondENCUT = 400ISTART = 1?ICHARG = 2PREC = AccurateNSW = 0IBRION = -1ISIF = 2EDIFF = 0.00001EDIFFG = -0.01LRPA = FALSENPAR = 4IALGO = 48LHFCALC = TRUEHFSCREEN = 0.2ALGO = AllTIME = 0.4ENCUTFO