| 查看: 1043 | 回复: 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之类的 |
回复此楼» 猜你喜欢
26考研一志愿中国石油大学(华东)305分求调剂
已经有4人回复
-
286求调剂
已经有5人回复
-
材料与化工（0856）304求 B区 调剂
已经有3人回复
-
268求调剂
已经有8人回复
-
265求调剂
已经有14人回复
-
350 本科985求调剂，求老登收留
已经有3人回复
-
【考研调剂】化学专业 281分，一志愿四川大学，诚心求调剂
已经有9人回复
-
初试 317
已经有3人回复
-
279分求调剂 一志愿211
已经有15人回复
-
279求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
3楼2012-10-04 18:39:29
yangxuezhang
金虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 2250.4
- 散金: 216
- 红花: 5
- 帖子: 564
- 在线: 211.1小时
- 虫号: 1967414
- 注册: 2012-09-01
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
2楼2012-10-04 18:37:13
yangxuezhang
金虫 (正式写手)
- 应助: 6 (幼儿园)
- 金币: 2250.4
- 散金: 216
- 红花: 5
- 帖子: 564
- 在线: 211.1小时
- 虫号: 1967414
- 注册: 2012-09-01
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
4楼2012-10-05 15:46:56
