24小时热门版块排行榜

返回列表

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

identation

金虫 (著名写手)

1ST强帖: 16
应助: 46 (小学生)
贵宾: 2.516
金币: 253.4
散金: 2603
红花: 44
帖子: 1283
在线: 500.5小时
虫号: 842155
注册: 2009-09-06
性别: GG
专业: 凝聚态物性 II ：电子结构

[交流] 【求助】PWSCF计算介电极化时为何出现如下提示？已有7人参与

==================================================
                        POLARIZATION CALCULATION
                     !!! NOT THOROUGHLY TESTED !!!
            --------------------------------------------------

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   from c_phase : error #       1
   Polarization only for insulators and no empty bands
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

   stopping ...
计算的是铁磁性体系，自洽计算正常，非自洽计算中出现这样的提示，请问是怎么回事呢？谢谢

回复此楼

» 猜你喜欢

依托企业入选了国家启明计划青年人才。有无高校可以引进的。已经有14人回复
表哥与省会女结婚，父母去帮带孩子被省会女气回家生重病了已经有11人回复
江汉大学解明教授课题组招博士研究生/博士后已经有3人回复
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性已经有11人回复
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗已经有10人回复
依托企业入选了国家启明计划青年人才。有无高校可以引进的。已经有11人回复

1楼 2010-06-17 09:02:00

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

identation

金虫 (著名写手)

1ST强帖: 16
应助: 46 (小学生)
贵宾: 2.516
金币: 253.4
散金: 2603
红花: 44
帖子: 1283
在线: 500.5小时
虫号: 842155
注册: 2009-09-06
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

Originally posted by valenhou001 at 2010-08-21 20:48:46:
"另外，计算得到的能带和DOS图对不上呀，能带图的带隙比DOS的宽很多呀，这是怎么回事呢"

这个较正常，态密度中，由于采用展宽的方法来确定态密度值，在价带边或导带边的能量处而出现非0的值，使得带 ...

"IF ((degauss > 0.01) .OR. (nbnd /= nelec/2)) CALL errore('c_phase', &
'Polarization only for insulators and no empty bands',1)
"
这是pwscf源程序的一句话，请问老师是不是一定要把nbnd指定成电子数的一半呢？是价电子的一半还是电子数的一半呢？

&control
calculation = 'nscf'
pseudo_dir = './'
outdir = './tmp'
lberry = .true.
gdir = 3
nppstr = 7
/
&system
ibrav = 1
celldm(1) = 7.3699
nat = 5
ntyp = 3
nbnd = 22
ecutwfc = 30.0
occupations = 'fixed'
degauss = 0.00
/
&electrons
conv_thr = 1e-5
mixing_beta = 0.3
/
ATOMIC_SPECIES
Pb 207.2 Pb.vdb.UPF
Ti 47.867 Ti.vdb.UPF
O 15.9994 O.vdb.UPF
ATOMIC_POSITIONS
Pb 0.000 0.000 0.010
Ti 0.500 0.500 0.500
O 0.000 0.500 0.500
O 0.500 0.500 0.000
O 0.500 0.000 0.500
K_POINTS {automatic}
4 4 7 1 1 1

上面的是计算PbTiO3的程序，这里的nbnd好像不是电子数的一半呀
麻烦老师赐教一下，谢谢！