| 查看: 478 | 回复: 0 | |||
[求助]
用espresso计算硫化钨的能带结构，无法有效设置k路径
|
|
大家好，我初学espresso，老师让我计算硫化钨的能带结构。 第一个输入文件是ws2.cry.scf.in，用来完成scf计算。 &control calculation = 'scf' restart_mode='from_scratch', prefix='ws2cry' tstress = .true. tprnfor = .true. pseudo_dir = '.', outdir='.' / &system ibrav=4, celldm(1)= 3.17, celldm(3) =3.9093, nat=6, ntyp= 2, ecutwfc =60.0, / &electrons diagonalization='david' mixing_mode = 'plain' mixing_beta = 0.7 conv_thr = 1.0d-8 / ATOMIC_SPECIES W 183.84 W.pbe-nsp-van.UPF S 32.065 S.pbe-van_bm.UPF ATOMIC_POSITIONS S 0.33333 0.66667 0.6225 S 0.33333 0.66667 -0.1225 S -0.33333 -0.66667 -0.6225 S -0.33333 -0.66667 0.1225 W -0.33333 -0.66667 -0.25 W 0.33333 0.66667 0.25 K_POINTS (automatic) 15 15 3 1 1 1 EOF 第二个输入文件是ws2.cry.band.in，用来计算能带。 &control calculation='bands' pseudo_dir = '.', outdir='.', prefix='ws2cry', / &system ibrav= 4, celldm(1)= 3.17, celldm(3) =3.9093, nat= 6, ntyp= 2, ecutwfc =60.0, ecutrho = 400.0, nbnd = 28, / &electrons diagonalization='david', / ATOMIC_SPECIES W 183.84 W.pbe-nsp-van.UPF S 32.065 S.pbe-van_bm.UPF ATOMIC_POSITIONS S 0.33333 0.66667 0.6225 S 0.33333 0.66667 -0.1225 S -0.33333 -0.66667 -0.6225 S -0.33333 -0.66667 0.1225 W -0.33333 -0.66667 -0.25 W 0.33333 0.66667 0.25 K_POINTS 2 0.0 0.0 0.0 4.0 0.36 0.0 0.0 1.0 EOF 还有bands.x的输入文件ws2.cry.bands.in： &bands prefix = 'ws2cry', outdir = '.', filband = 'ws2cry.dat', lsym=.true, / 先执行pw.x <ws2.cry.scf.in> ws2.cry.scf.out，然后执行pw.x <ws2.cry.band.in> ws2.cry.band.out，最后执行bands.x <ws2.cry.bands.in> ws2.cry.bands.out，其中ws2.cry.band.in的k points被完全无视了，我试过各种不同组合的k points，但输出文件ws2cry.dat中的k points数量非常多，经过比对是ws2.cry.scf.in中自动生成的k points。另外输出文件有26个能带，和我的设定不一样（nbnd=28)。我之前用espresso算过铜的能带结构，输出文件的k points与能带的输入文件完全一致，为什么这一次输出文件直接无视了能带的输入文件，而用的是能带的输入文件呢？还有一个细节就是做scf计算的时间很长，而计算能带结构用的的时间却很短，几乎是秒出结果。 |
回复此楼» 猜你喜欢
上海工程技术大学【激光智能制造】课题组招收硕士
已经有6人回复
-
带资进组求博导收留
已经有11人回复
-
自荐读博
已经有5人回复
-
求个博导看看
已经有16人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
-
求助院士们，这个如何合成呀
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
-
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）
已经有3人回复