应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Abinit » 【求助】nband和occ 设置
51/1返回列表
查看: 1999  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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 关于这怎么 觉的 不一样呐
回复此楼

» 猜你喜欢

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

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

查看全部散金贴
1楼 2010-12-23 20:42:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzx2008

荣誉版主 (著名写手)
★
gzqdyouxia(金币+1):谢谢指点 2010-12-24 11:03:18
若occopt=2，需手动设置占据occ，且此时nband（nkpt*nsppol），是nkpt*nsppol维数组，occ（nband），先按k点排列自旋朝上，然后自旋朝下
赞一下(1人) 回复此楼
4楼2010-12-24 09:22:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

identation

金虫 (著名写手)
★ ★
youzhizhe(金币+2):谢谢指导。 2010-12-23 23:37:24
nband在自洽计算的时候，如果是不考虑自旋绝缘体一般可指定为总电子数的一半，如果是金属一般要多余电子数目一半的百分之二十；
occ是电子占有数，考虑自旋的系统，你可以分别指定自旋上和自旋下的电子，两者总数和是电子数目；其实是自旋上和自旋下的occ分别是两个数组，应该是nband维的数组
赞一下(1人) 回复此楼
2楼2010-12-23 22:55:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mazuju028

木虫 (正式写手)
★ ★
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.
赞一下(1人) 回复此楼
3楼2010-12-24 09:06:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭