| 查看: 2222 | 回复: 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 ] |
回复此楼» 猜你喜欢
FileCode能看出啥？
已经有23人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
关于Filecode分析方法
已经有9人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有10人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp自旋极化计算时，初始设置自旋能否不被优化
已经有6人回复
- 关于siesta计算铁的结合能的问题
已经有4人回复
- VASP 算不准氧气键能
已经有16人回复
wtianxing
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 296.6
- 帖子: 86
- 在线: 72.6小时
- 虫号: 438680
- 注册: 2007-09-27
- 专业: 凝聚态物性I:结构、力学和
10楼2013-12-21 09:36:04
wtianxing
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 296.6
- 帖子: 86
- 在线: 72.6小时
- 虫号: 438680
- 注册: 2007-09-27
- 专业: 凝聚态物性I:结构、力学和
2楼2013-12-17 08:29:53
3楼2013-12-17 18:19:00
wtianxing
银虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 296.6
- 帖子: 86
- 在线: 72.6小时
- 虫号: 438680
- 注册: 2007-09-27
- 专业: 凝聚态物性I:结构、力学和
4楼2013-12-18 08:36:07
