| 查看: 1212 | 回复: 4 | |||
yangxuezhang金虫 (正式写手)
|
[求助]
关于siesta计算铁的结合能的问题
|
|
因为是新手，想从简单的问题算起，不过算过之后发现算出来的结合能比实验数字大了蛮多，实际是4.28左右，算出来是5.6左右，不知道问题出在哪，附上我的输入文件，希望高手指教一下 SystemName bcc Fe ferro GGA # Descriptive name of the system SystemLabel Fe # Short name for naming files # Output options WriteCoorXmol WriteMullikenPop 1 # Species and atoms NumberOfSpecies 1 NumberOfAtoms 1 %block ChemicalSpeciesLabel 1 26 Fe %endblock ChemicalSpeciesLabel # Basis PAO.EnergyShift 50 meV PAO.BasisSize DZP %block PAO.Basis Fe 2 0 2 P 6. 0. 2 2 0. 0. %endblock PAO.Basis LatticeConstant 2.87 Ang %block LatticeVectors 0.50000 0.500000 0.500000 0.50000 -0.500000 0.500000 0.50000 0.500000 -0.500000 %endblock LatticeVectors KgridCutoff 15. Ang %block BandLines 1 0.00000 0.000000 0.000000 \Gamma 40 2.00000 0.000000 0.000000 H 28 1.00000 1.000000 0.000000 N 28 0.00000 0.000000 0.000000 \Gamma 34 1.00000 1.000000 1.000000 P %endblock BandLines xc.functional GGA # Exchange-correlation functional xc.authors PBE # Exchange-correlation version SpinPolarized true # Logical parameters are: yes or no MeshCutoff 150. Ry # Mesh cutoff. real space mesh # SCF options MaxSCFIterations 40 # Maximum number of SCF iter DM.MixingWeight 0.1 # New DM amount for next SCF cycle DM.Tolerance 1.d-3 # Tolerance in maximum difference # between input and output DM DM.UseSaveDM true # to use continuation files DM.NumberPulay 3 SolutionMethod diagon # OrderN or Diagon ElectronicTemperature 25 meV # Temp. for Fermi smearing # MD options MD.TypeOfRun cg # Type of dynamics: MD.NumCGsteps 0 # Number of CG steps for # coordinate optimization MD.MaxCGDispl 0.1 Ang # Maximum atomic displacement # in one CG step (Bohr) MD.MaxForceTol 0.04 eV/Ang # Tolerance in the maximum # atomic force (Ry/Bohr) # Atomic coordinates AtomicCoordinatesFormat Fractional %block AtomicCoordinatesAndAtomicSpecies 0.000000000000 0.000000000000 0.000000000000 1 %endblock AtomicCoordinatesAndAtomicSpecies 然后把晶格常数改成10A，算出单个原子的能量，当然晶格常数我取了很多的点2.85.2.88之类的 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有20人回复
-
今天务委会开完了，明天出结果吗
已经有22人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有9人回复
-
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
放榜前的不淡定
已经有19人回复
» 本主题相关价值贴推荐，对您同样有帮助:
yangxuezhang
金虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 2250.4
- 散金: 216
- 红花: 5
- 帖子: 564
- 在线: 211.1小时
- 虫号: 1967414
- 注册: 2012-09-01
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2012-10-04 18:37:13
3楼2012-10-04 18:39:29
yangxuezhang
金虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 2250.4
- 散金: 216
- 红花: 5
- 帖子: 564
- 在线: 211.1小时
- 虫号: 1967414
- 注册: 2012-09-01
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
4楼2012-10-05 15:46:56
yangxuezhang
金虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 2250.4
- 散金: 216
- 红花: 5
- 帖子: 564
- 在线: 211.1小时
- 虫号: 1967414
- 注册: 2012-09-01
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
5楼2012-10-10 19:33:34
28