24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

wilsun7

银虫 (初入文坛)

应助: 2 (幼儿园)
金币: 233.6
帖子: 45
在线: 45.4小时
虫号: 1863760
注册: 2012-06-19
专业: 凝聚态物性 II ：电子结构

[求助] 计算费米面时遇到的问题

大家好，本人新手，在计算费米面时遇到一个问题：在运行 bands_FS.x > bands_FS.out 时没有任何反应，查看系统监视器，bands_FS.x 根本没有运行。希望大神帮我看看问题处在那。计算是参考的这个网站： http://valenhou.blog.edu.cn/2006/133295.html

以下附上本人的计算过程，希望和大家一起研究

第一步：自洽计算（pw.x < scf-FeSe.in > scf-FeSe.out）
输入文件：
&CONTROL
               calculation = 'scf' ,
            restart_mode = 'from_scratch' ,
               wf_collect = .true. ,
                  outdir = './tmp' ,
               pseudo_dir = '/home/wilson/Downloads/espresso-4.3.1/Psudoptential_files',
                  prefix = 'FeSe' ,
               verbosity = 'high'
/
&SYSTEM
                     ibrav = 0,
               celldm(1) = 7.1148,
                     nat = 4,
                     ntyp = 2,
                     nspin = 1,
                  ecutwfc = 30.0 ,
                  ecutrho = 300.0 ,
               occupations = 'smearing' ,
                  degauss = 0.01 ,
                  smearing = 'methfessel-paxton' ,
/
&ELECTRONS
electron_maxstep = 500,
                  conv_thr = 1.0e-10 ,
               mixing_beta = 0.1 ,
         diagonalization = 'david' ,
/
CELL_PARAMETERS
   1.000000000 0.000000000 0.000000000
   0.000000000 1.000000000 0.000000000
   0.000000000 0.000000000 1.465600000
ATOMIC_SPECIES
Fe  55.85 Fe.pbe-sp-van.UPF
Se  78.96 Se.pbe-van.UPF
ATOMIC_POSITIONS crystal
Fe 0.250000 0.75000 0.00000
Fe 0.750000 0.25000 0.00000
Se 0.250000 0.25000 0.26000
Se 0.75000 0.75000 0.74000
K_POINTS automatic
8  8  6 1 1 1

第二步：用kvecs_FS.x产生要计算费米面时的密集网格k点（kvecs_FS.x < kvecs_FS.in > kvecs_FS.out）
kvecs_FS.in 文件如下：
kvecs_FS.in
1.000000 0.000000 0.000000
0.000000 1.000000 0.000000
0.000000 0.000000 0.682314
16 16 16
FeSe

然后 kvecs_FS.x < kvecs_FS.in > kvecs_FS.out
这时会生成文件kvecs_fese 包含了所有k点

第三步：进行一次非自洽计算，计算这些k点的本征值(cat scf.in kvecs_FeSe > nscf.in)
把 scf 的输入文件和生成的k点附加到nscf中。

输入文件：
&CONTROL
               calculation = 'nscf' ,
            restart_mode = 'from_scratch' ,
               wf_collect = .true. ,
                  outdir = './tmp' ,
               pseudo_dir = '/home/wilson/Downloads/espresso-4.3.1/Psudoptential_files',
                  prefix = 'FeSe' ,
               verbosity = 'high'
/
&SYSTEM
                     ibrav = 0,
               celldm(1) = 7.1148,
                     nat = 4,
                     ntyp = 2,
                     nspin = 1,
                  ecutwfc = 30.0 ,
                  ecutrho = 300.0 ,
               occupations = 'smearing' ,
                  degauss = 0.01 ,
                  smearing = 'methfessel-paxton' ,
/
&ELECTRONS
electron_maxstep = 500,
                  conv_thr = 1.0e-10 ,
               mixing_beta = 0.1 ,
         diagonalization = 'david' ,
/
CELL_PARAMETERS
   1.000000000 0.000000000 0.000000000
   0.000000000 1.000000000 0.000000000
   0.000000000 0.000000000 1.465600000
ATOMIC_SPECIES
Fe  55.85 Fe.pbe-sp-van.UPF
Se  78.96 Se.pbe-van.UPF
ATOMIC_POSITIONS crystal
Fe 0.250000 0.75000 0.00000
Fe 0.750000 0.25000 0.00000
Se 0.250000 0.25000 0.26000
Se 0.75000 0.75000 0.74000
K_POINTS
  4913
0.000000 0.000000 0.000000  1.00
0.000000 0.000000 0.042645  1.00
0.000000 0.000000 0.085289  1.00
   .....             .....          ......          ...
再run pw.x < nscf-FeSe.in > nscf-FeSe.out

第四步：用bands_FS.x将计算的k点以及相应的本征值，转换成xcrysden软件的bxsf格式。
首先把 nscf-FeSe.out 复制到 Bands.out
然后新建文件input_FS ：
20 24
8.8449
FeSe
16 16 16
1.000000 0.000000 0.000000
0.000000 1.000000 0.000000
0.000000 0.000000 0.682314

再run bands_FS.x > bands_FS.out
然后应该就可以了。我的问题出现在最后一步，前面3步一切正常。最后一步 bands_FS.x > bands_FS.out 后就没反应了。。。希望大家帮我看看问题出在哪儿。谢谢

回复此楼

1楼2012-06-27 06:44:38

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

1ST强帖: 13
应助: 241 (大学生)
金币: 25701.7
散金: 602
红花: 166
帖子: 3782
在线: 873.8小时
虫号: 1007127
注册: 2010-04-27
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

在pw.x的计算中设置合理的nbnd的值，比如体系的总的价电子数/2 + 10或更大的数（整数）

在bands_FS.x的输入文件 input_FS中设置n_start， n_last的值要合理，要在Fermi energy所穿过的能带数的左右。
Fermi energy所穿过能带数可通过体系的总价电子数来估计。

bands_FS.x除了读入input_FS之外，还要读入前一步的pw.x进行本征值计算的out文件，运行的命令行：
bands_FS.x < 前一步的pw.x进行本征值计算的out文件 > log.out
或
bands_FS.x < 前一步的pw.x进行本征值计算的out文件