| 查看: 1835 | 回复: 9 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
wtianxing银虫 (小有名气)
|
[求助]
用siesta算Li在Si中结合能，是否开自旋？ 已有1人参与
|
|
|
用siesta计算单个Li在64个完整Si晶体中的结合能，公式如下： E（64Si+Li）-E（64Si）-E（Li） 其中E（64Si+Li）是将Li放入含64个Si的晶体中，用周期边界条件算的，E（64Si）是直接计算含64个Si晶体的能量，E（Li）是将一个锂加三边真空层算的。 问题出在单个锂用真空层算能量上，考虑自旋（SpinPolaried .true.）后的能量要比不考虑的低0.32eV,而其他两个能量用不用自旋没有影响。 按说考虑自旋后能量更低，应该是对的，但是算出的结合能与其他文献采用vasp（不知道考没考虑自旋）的对不上，就差0.3左右，而不考虑自旋就对上了。请高手帮忙，非常感谢。 算单个Li的输入文件： Siesta Version: siesta-2.0.1 Architecture : Ifort(10.1.012)+MKL(10.0.5.025)+OpenMPI(1.2.6)@openSUSE-10.3_X86-64 Compiler flags: mpif90 -O2 -i-static SystemLabel work NumberOfAtoms 1 NumberOfSpecies 1 %block ChemicalSpeciesLabel 1 3 Li %endblock ChemicalSpeciesLabel PAO.BasisType split PAO.BasisSize DZP PAO.EnergyShift 0.010000 Ry LatticeConstant 1.0 Ang %block LatticeParameters 21.6804 23.578 15.30271078 90.000000 90.000000 90.000000 %endblock LatticeParameters AtomicCoordinatesFormat NotScaledCartesianAng %block AtomicCoordinatesAndAtomicSpecies 12.605130 5.880100 8.131973 1 %endblock AtomicCoordinatesAndAtomicSpecies #---------------------------------------- # Energy control parameters for SCF cycles MeshCutoff 200.0 Ry # Energy cut-off for the planewave %block kgrid_Monkhorst_Pack # k mesh points for integration 1 0 0 0.500 0 1 0 0.500 0 0 1 0.500 %endblock kgrid_Monkhorst_Pack MaxSCFIterations 300 DM.MixingWeight 0.1 DM.NumberPulay 10 DM.Tolerance 1.d-4 DM.UseSaveDM SolutionMethod diagon ElectronicTemperature 300 K XC.functional GGA XC.authors PBE SpinPolarized .true. 开之后能量低0.3左右 #--------------------------------------- # Control the files printed WriteKpoints .false. WriteForces .true. WriteCoorXmol .true. WriteCoorCerius .true. WriteCoorStep .true. #--------------------------------------- MD.TypeOfRun CG MD.NumCGsteps 300 MD.MaxCGDispl 0.15 Ang MD.MaxForceTol 0.010 eV/Ang MD.VariableCell .false. Diag.DivideAndConquer .false. ************************** End of input data file ***************************** [ Last edited by wtianxing on 2013-12-16 at 15:51 ] |
回复此楼» 猜你喜欢
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？
已经有3人回复
-
自荐读博
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有5人回复
-
论文终于录用啦！满足毕业条件了
已经有22人回复
-
不自信的我
已经有5人回复
-
磺酰氟产物，毕不了业了！
已经有4人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp自旋极化计算时，初始设置自旋能否不被优化
已经有6人回复
- 关于siesta计算铁的结合能的问题
已经有4人回复
- VASP 算不准氧气键能
已经有16人回复
gmmy
木虫 (小有名气)
- 应助: 7 (幼儿园)
- 金币: 3487.4
- 红花: 1
- 帖子: 103
- 在线: 262.1小时
- 虫号: 447465
- 注册: 2007-11-01
- 专业: 凝聚态物性 II ：电子结构
7楼2013-12-18 18:08:28
gmmy
木虫 (小有名气)
- 应助: 7 (幼儿园)
- 金币: 3487.4
- 红花: 1
- 帖子: 103
- 在线: 262.1小时
- 虫号: 447465
- 注册: 2007-11-01
- 专业: 凝聚态物性 II ：电子结构
9楼2013-12-19 21:23:27