| 查看: 181 | 回复: 2 | |||
| 当前主题已经存档。 | |||
pengdou至尊木虫 (正式写手)
|
[交流]
【求助】求Si(111)表面7x7重构的坐标
|
||
|
求Si(111)表面7x7重构的坐标,包括元胞基矢 [ Last edited by pengdou on 2008-12-24 at 15:58 ] |
回复此楼» 猜你喜欢
国自然面上和省基金B类撒花
已经有16人回复
-
有没有学校收留
已经有3人回复
-
312求调剂
已经有3人回复
-
华师大读博
已经有5人回复
-
又一批高校组建人工智能学院 师资行吗 不是骗人吗
已经有5人回复
-
急需审稿人！！！
已经有3人回复
-
申博/考博
已经有8人回复
-
295分求调剂
已经有6人回复
-
085600材料与化工调剂
已经有6人回复
-
有没有接收比较快的sci期刊呀，最好在一个月之内的，研三孩子求毕业
已经有7人回复
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):非常之感谢！！！
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
zxzj05
荣誉版主 (著名写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 贵宾: 1.242
- 金币: 5912.5
- 红花: 2
- 帖子: 1528
- 在线: 144.5小时
- 虫号: 406147
- 注册: 2007-06-18
- 性别: GG
- 专业: 金属功能材料
- 管辖: 第一性原理
3楼2008-12-24 21:00:21
5