| 查看: 6058 | 回复: 7 | ||
| 【奖励】 本帖被评价4次，作者漫天飘雪增加金币 2.8 个 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[资源]
Quantum expresso计算计算费米面
|
||
|
本教程转发自https://blog.sciencenet.cn/home.php?mod=space&uid=3352196&do=blog&quickforward=1&id=1070348 1.普通自洽计算(pw.x<scf.in>scf.out) &CONTROL calculation = 'scf' prefix = 'silicon' restart_mode = 'from_scratch' wf_collect = .false. pseudo_dir = './', outdir = './' disk_io = 'low' tstress = .true. tprnfor = .true. etot_conv_thr = 1.0d-5 forc_conv_thr = 1.0d-4 / &SYSTEM ibrav = 2, celldm(1) = 10.2, nat = 2, ntyp = 1, ecutwfc = 50 ecutrho = 800 occupations = 'smearing' smearing = 'mp' degauss = 0.02 nbnd = 28 noncolin = .true. lspinorb = .true. / &electrons diagonalization = 'david' mixing_mode = 'plain' mixing_beta = 0.7 conv_thr = 1.0d-9 / ATOMIC_SPECIES Si 28.08600 Si.pbe-n-kjpaw_psl.0.1.UPF ATOMIC_POSITIONS {crystal} Si 0.000000000 0.000000000 0.000000000 Si 0.250000000 0.250000000 0.250000000 K_POINTS {automatic} 6 6 6 0 0 0 2.用密集网格k点进行一次非自洽计算(pw.x<nscf.in>nscf.out) &CONTROL calculation = 'nscf' prefix = 'silicon' restart_mode = 'from_scratch' wf_collect = .false. pseudo_dir = './', outdir = './' disk_io = 'low' tstress = .true. tprnfor = .true. etot_conv_thr = 1.0d-5 forc_conv_thr = 1.0d-4 verbosity = 'high' / &SYSTEM ibrav = 2, celldm(1) = 10.2, nat = 2, ntyp = 1, ecutwfc = 50 ecutrho = 738 occupations = 'smearing' smearing = 'mp' degauss = 0.02 nbnd = 28 noncolin = .true. lspinorb = .true. / &electrons diagonalization = 'david' mixing_mode = 'plain' mixing_beta = 0.7 conv_thr = 1.0d-9 / ATOMIC_SPECIES Si 28.08600 Si.pbe-n-kjpaw_psl.0.1.UPF ATOMIC_POSITIONS {crystal} Si 0.000000000 0.000000000 0.000000000 Si 0.250000000 0.250000000 0.250000000 K_POINTS {automatic} 40 40 40 0 0 0 3.将计算的k点以及相应的本征值，转换成xcrysden软件的bxsf格式(fs.x<fs.in>fs.out) &fermi outdir='./' prefix='silicon' / 4.采用xcrysden来画图,输入文件为silicon_fs.bxsf |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : silicon.zip
2017-09-27 08:48:11, 1.92 M
» 猜你喜欢
建议基金发布提前给出明确的时间点
已经有13人回复
-
2026国自然函评费到账
已经有7人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有20人回复
-
科研孤儿太难了
已经有18人回复
-
看来今天不会放榜了？
已经有11人回复
-
时间戳又变了
已经有19人回复
-
filecode，4个jtjc了
已经有14人回复
-
让我中一个面上吧！
已经有13人回复
-
人气不行了
已经有7人回复
-
感觉是下周放榜了
已经有12人回复
5楼2017-10-14 12:27:10
dft_family
禁虫 (小有名气)
|
本帖内容被屏蔽 |
4楼2017-10-14 09:38:47
6楼2017-12-02 19:50:33
7楼2021-06-28 23:49:05
简单回复
李蕊1993878楼
2022-04-27 17:06
回复
五星好评 顶一下，感谢分享！
40