24小时热门版块排行榜    

查看: 2217  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 能否退出参与的面上项目解除限项 +19 koalala 2026-08-24 21/1050 2026-08-24 16:48 by Wei_ren
[基金申请] 没有任何消息-是不是就凉了 +3 图啦图啦 2026-08-24 3/150 2026-08-24 16:43 by GZUZM
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 我面上完蛋了 +8 且听虎啸 2026-08-20 9/450 2026-08-24 13:52 by zzuzxg
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 18/900 2026-08-24 12:03 by iaeyuan
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见