| 查看: 1037 | 回复: 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之类的 |
回复此楼» 猜你喜欢
08工学调剂
已经有5人回复
-
281求调剂（0805）
已经有25人回复
-
085600材料与化工
已经有6人回复
-
265求调剂
已经有8人回复
-
0817 化学工程 299分求调剂 有科研经历 有二区文章
已经有21人回复
-
295材料求调剂，一志愿武汉理工085601专硕
已经有4人回复
-
0856调剂，是学校就去
已经有9人回复
-
298-一志愿中国农业大学-求调剂
已经有9人回复
-
求调剂
已经有4人回复
-
学校已经提交到NSFC，还能修改吗？
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
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
10