24小时热门版块排行榜    

查看: 229  |  回复: 2
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

pengdou

至尊木虫 (正式写手)

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

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

[ Last edited by pengdou on 2008-12-24 at 15:58 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxzj05

荣誉版主 (著名写手)

★ ★ ★
zdhlover(金币+3,VIP+0):多谢,欢迎常来
zdhlover(金币+0,VIP+0):包含圣诞特别奖励
按照这个教程边学边用,
建模型表面时后实在不会就可以学学MScleave表面,它这方面很强。
然后,把坐标转换成abinit所用的xred格式。

很多时候大家把自己的事情都忙不过来,
所以最好还是学方法,
自己干,这样学的扎实!
储氢家族欢迎储氢研究者!
3楼2008-12-24 21:00:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

zxzj05

荣誉版主 (著名写手)

★ ★ ★ ★ ★ ★ ★
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
储氢家族欢迎储氢研究者!
2楼2008-12-24 20:56:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 08:14 by C79jjtjAKjEn
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 07:58 by C79jjtjAKjEn
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 s3fFTmArrBt6 2026-09-14 6/300 2026-09-15 07:38 by C79jjtjAKjEn
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 s3fFTmArrBt6 2026-09-14 6/300 2026-09-15 07:14 by C79jjtjAKjEn
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:49 by 5BDX0d0WFp7t
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-15 05:49 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-15 05:34 by 5BDX0d0WFp7t
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 05:13 by 5BDX0d0WFp7t
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 4/200 2026-09-14 23:24 by vZfe6xYu34yj
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 22:45 by vZfe6xYu34yj
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:13 by vZfe6xYu34yj
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:36 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 15:09 by Aj1rhIDL5ixY
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见