| 查看: 1060 | 回复: 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之类的 |
回复此楼» 猜你喜欢
296求调剂
已经有14人回复
-
生物学调剂
已经有6人回复
-
300分求调剂 （085501机械专硕，本科扬大）
已经有3人回复
-
283求调剂 086004考英二数二
已经有6人回复
-
农学0904 312求调剂
已经有5人回复
-
085600材料与化工329分求调剂
已经有10人回复
-
复试调剂
已经有9人回复
-
材料与化工调剂
已经有35人回复
-
调剂求收留
已经有10人回复
-
一志愿哈工大 085600 277 12材科基求调剂
已经有18人回复
» 本主题相关价值贴推荐，对您同样有帮助:
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
