| 查看: 480 | 回复: 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计算的时间很长，而计算能带结构用的的时间却很短，几乎是秒出结果。 |
回复此楼» 猜你喜欢
2025年遐想
已经有4人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
-
自然科学基金委宣布启动申请书“瘦身提质”行动
已经有4人回复
-
求个博导看看
已经有18人回复