24小时热门版块排行榜    

查看: 2018  |  回复: 3

xsw841020

铜虫 (初入文坛)


[交流] 【求助】nband和occ 设置

ixc 7
nsppol 2
nbanddsto -1

# Crystalline aluminum : optimization of the lattice parameter
# at fixed number of k points and broadening.

#Definition of occupation numbers
occopt 2
nband 3 3 3 3 3 3 3 3
occ 1 1 0 1 1 0 1 1 0 1 1 0 1 0 0 1 0 0 1 0 0 1 0 0
tsmear 0.05

#Definition of the unit cell
acell 3*7.60           # This is equivalent to   7.60 7.60 7.60
rprim  0.0  0.5  0.5   # FCC primitive vectors (to be scaled by acell)
       0.5  0.0  0.5
       0.5  0.5  0.0

#Optimization of the lattice parameters
#optcell 1
#ionmov  3
ntime  10
dilatmx 1.05
#ecutsm  0.5

#Definition of the atom types
ntype 1           # There is only one type of atom
zatnum 13         # The keyword "zatnum" refers to the atomic number of the
                  # possible type(s) of atom. The pseudopotential(s)
                  # mentioned in the "files" file must correspond
                  # to the type(s) of atom. Here, the only type is Aluminum
                        

#Definition of the atoms
natom 1           # There is only one atom per cell
type 1            # This atom is of type 1, that is, Aluminum
xred              # This keyword indicate that the location of the atoms
                  # will follow, one triplet of number for each atom
   0.0  0.0  0.0     # Triplet giving the REDUCED coordinate of atom 1.

#Definition of the planewave basis set
ecut  6.0         # Maximal kinetic energy cut-off, in Hartree

#Definition of the k-point grid
kptopt 2          # Option for the automatic generation of k points, taking
                  # into account the symmetry
ngkpt 2 2 2       # This is a 2x2x2 FCC grid, based on the primitive vectors
nshiftk 1         # of the reciprocal space. For a FCC real space lattice,
                  # like the present one, it actually corresponds to the
                  # so-called 4x4x4 Monkhorst-Pack grid, if the following shifts
                  # are used :
shiftk 0.5 0.5 0.5
#       0.5 0.0 0.0
#       0.0 0.5 0.0
#       0.0 0.0 0.5

#Definition of the SCF procedure
nstep 10          # Maximal number of SCF cycles
toldfe 1.0d-6     # Will stop when, twice in a row, the difference
                  # between two consecutive evaluations of total energy
                  # differ by less than toldfe (in Hartree)
新手  我就想问下 怎么 得出 nband 和 occ 是这些呐  我看 tutorial 关于这怎么 觉的 不一样呐
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
youzhizhe(金币+2):谢谢指导。 2010-12-23 23:37:24
nband在自洽计算的时候,如果是不考虑自旋绝缘体一般可指定为总电子数的一半,如果是金属一般要多余电子数目一半的百分之二十;
occ是电子占有数,考虑自旋的系统,你可以分别指定自旋上和自旋下的电子,两者总数和是电子数目;其实是自旋上和自旋下的occ分别是两个数组,应该是nband维的数组
2楼2010-12-23 22:55:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★
gzqdyouxia(金币+2):谢谢指点 2010-12-24 11:03:07
你这个例子中的nband,ocd的形式是由occopt=2决定的,当occopt=2时,表示可以标定每个k点对应的带数(nband),还可以标定每个能级上电子占据数(occ)。
所以你的nband 3 3 3 3 3 3 3 3,表示的是4个K点,每个K点都有自旋向上和向下电子,且占据三个能级。
occ  1 1 0 1 1 0 1 1 0 1 1 0 1 0 0 1 0 0 1 0 0 1 0 0 表示4个K点上,自旋向上电子排布为1 1 0,自旋向下电子排布为1 0 0.
3楼2010-12-24 09:06:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gzqdyouxia(金币+1):谢谢指点 2010-12-24 11:03:18
若occopt=2,需手动设置占据occ,且此时nband(nkpt*nsppol),是nkpt*nsppol维数组,occ(nband),先按k点排列自旋朝上,然后自旋朝下
4楼2010-12-24 09:22:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xsw841020 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +5 Liwangman 2026-03-15 5/250 2026-03-16 17:10 by 我的船我的海
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 311求调剂 +5 26研0 2026-03-15 5/250 2026-03-16 16:21 by a不易
[考研] 085600材料与化工 求调剂 +13 enenenhui 2026-03-13 14/700 2026-03-16 15:19 by 了了了了。。
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 331求调剂(0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 311求调剂 +8 zchqwer 2026-03-10 8/400 2026-03-14 00:01 by JourneyLucky
[考研] 一志愿西南交大,材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 0703化学求调剂 +7 绿豆芹菜汤 2026-03-12 7/350 2026-03-13 17:25 by njzyff
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 26考研求调剂 +5 丶宏Sir 2026-03-13 5/250 2026-03-13 13:05 by JourneyLucky
[考研] 279求调剂 +3 莫xiao 2026-03-10 4/200 2026-03-11 08:06 by 斩魂滴兔子!
[考研] 0856材料与化工353分求调剂 +11 NIFFFfff 2026-03-09 11/550 2026-03-10 18:36 by suyuanhai
信息提示
请填处理意见