24小时热门版块排行榜

返回列表

当前主题已经存档。

pengdou

至尊木虫 (正式写手)

应助: 0 (幼儿园)
金币: 22671.7
帖子: 957
在线: 182小时
虫号: 461689
注册: 2007-11-18
专业: 凝聚态物性 II ：电子结构

[交流] 【求助】求Si(111)表面7x7重构的坐标

求Si(111)表面7x7重构的坐标,包括元胞基矢

[ Last edited by pengdou on 2008-12-24 at 15:58 ]

回复此楼

» 猜你喜欢

酰胺脱乙酰基已经有11人回复
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗已经有4人回复
有时候真觉得大城市人没有县城人甚至个体户幸福已经有10人回复
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）已经有5人回复
有院领导为了换新车，用横向课题经费买了俩车已经有5人回复
CSC & MSCA 博洛尼亚大学能源材料课题组博士/博士后招生｜MSCA经费充足、排名优已经有6人回复
面上项目申报已经有3人回复
博士延得我，科研能力直往上蹿已经有7人回复
面上基金申报没有其他的参与者成吗已经有5人回复
遇见不省心的家人很难过已经有22人回复

1楼 2008-12-24 14:04:25

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

zxzj05

荣誉版主 (著名写手)

1ST强帖: 1
应助: 0 (幼儿园)
贵宾: 1.242
金币: 5912.5
红花: 2
帖子: 1528
在线: 144.5小时
虫号: 406147
注册: 2007-06-18
性别: GG
专业: 金属功能材料
管辖: 第一性原理

★ ★ ★ ★ ★ ★ ★
zdhlover(金币+5,VIP+0):圣诞节特别奖励
pengdou(金币+2,VIP+0):非常之感谢！！！

五.  表面 (Aluminum)
  5-1 总能和晶格常数计算
#Definition of occupation numbers
   occopt 4             # Occupation option , control how input parameters nband, occ and wtk
   tsmear 0.05          # Temperature of smearing (default 0.04)
#Optimization of the lattice parameters
   optcell 1             # Optimize the cell shape and dimesions
   ionmov 3             # Ionic moves (default 0)
   ntime 10          # Number of time steps
   dilatmx 1.05          # Dilatation : maximal value
   ecutsm 0.5          # Energy cutoff smearing
#Definition of the atoms
   kptopt 1             # Kpoint option (default 0.0)
   nshiftk 4             # Number of shifts of k point grids
   shiftk 0.5 0.5 0.5 # These shifts will be the same for all grids
            0.5 0.0 0.0
            0.0 0.5 0.0
            0.0 0.0 0.5
   ngkpt 2  2  2    # Number of grid points for k points generation
   nstep 10          # Maximal number of SCF cycles
      ecut 6.0          # Maximal kinetic energy cut-off, in Hartree
      ixc 1             # LDA Teter Pade parametrization
   ntypat 1             # There is only one type of atom
   znucl 13          # The atomic number of the type(s) of atom
   natom 1          # There is only one atom per cell
   toldfe 1.0d-6    # Tolerance on the difference of total energy
   typat 1             # This atom is of type 1, that is, Aluminum
   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
      xred 0.0  0.0  0.0  # The location of the atoms
Note :
1)这里ecut= 6Ha是一个很粗略的值，在真实计算时必须做ecut的收敛计算以得到正确的ecut
2)不须要用 diemac，因为 diemac 的默认值已设定在金属
计算结束后，在output file里找Fermi 这个关键词可以找到 Fermi energy 为 0.26800 Hatree
      Fermi energy (hartree) = 0.26800
      Eigenvalues (hartree) for nkpt= 2  k points:
      kpt# 1, nband=  3, wtk=  0.75000, kpt= -0.2500  0.5000  0.0000 (reduced coord)
         0.09391 0.25391 0.41846
            occupation numbers for kpt# 1  2.00003 1.33306 0.00014
      kpt# 2, nband=  3, wtk=  0.25000, kpt= -0.2500  0.0000  0.0000 (reduced coord)
      -0.07058 0.41033 0.68787
            occupation numbers for kpt# 2  2.00000 0.00030 0.00000