| 查看: 1242 | 回复: 6 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
请教，为什么会出现这样的情况？谢谢大家
|
||
|
我在测试mesh cutoff的时候，一直很顺利，mesh cutoff从200Ry一直增加到750ry，200到700ry的测试中都很顺利完成了scf运算，但是当mesh cutoff增加到750ry时，出现了这样的输出情况 InitMesh: MESH = 80 x 80 x 80 = 512000 InitMesh: Mesh cutoff (required, used) = 750.000 870.761 Ry * Maximum dynamic memory allocated = 168 MB stepf: Fermi-Dirac step function siesta: Program's energy decomposition (eV): siesta: Ebs = -48733249.309373 siesta: Eions = 10539.985266 siesta: Ena =***************** siesta: Ekin = 5039.821056 siesta: Enl = -2320.791946 siesta: DEna = -0.148201 siesta: DUscf = 0.000000 siesta: DUext = 0.000000 siesta: Exc = -1407.884017 siesta: eta*DQ = 0.000000 siesta: Emadel = 0.000000 siesta: Emeta = 0.000000 siesta: Emolmec = 0.000000 siesta: Ekinion = 0.000000 siesta: Eharris =***************** siesta: Etot =***************** siesta: FreeEng =***************** siesta: iscf Eharris(eV) E_KS(eV) FreeEng(eV) dDmax Ef(eV) siesta: 1******************************************390.3300******** timer: Routine,Calls,Time,% = IterSCF 1 2189.282 98.32 elaps: Routine,Calls,Wall,% = IterSCF 1 327.087 98.11 siesta: 2****************************************** 0.0553******** siesta: 3****************************************** 0.0453******** siesta: 4****************************************** 0.0027******** siesta: 5****************************************** 0.0030******** siesta: 6****************************************** 0.0009******** siesta: 7****************************************** 0.0007******** siesta: 8****************************************** 0.0003******** siesta: 9****************************************** 0.0002******** siesta: 10****************************************** 0.0002******** siesta: 11****************************************** 0.0002******** siesta: 12****************************************** 0.0002******** siesta: 13****************************************** 0.0002******** siesta: 14****************************************** 0.0001******** siesta: 15****************************************** 0.0001******** siesta: 16****************************************** 0.0002******** siesta: 17****************************************** 0.0002******** siesta: 18****************************************** 0.0002******** siesta: 19****************************************** 0.0002******** siesta: 20****************************************** 0.0002******** siesta: 21****************************************** 0.0001******** 请问是为什么呢？ 我的输入文件如下： SystemName SystemLabel NumberOfAtoms 10 NumberOfSpecies 4 %block ChemicalSpeciesLabel 1 * * # Species index, atomic number, species label 2 * * 3 * * 4 * * %endblock ChemicalSpeciesLabel LatticeConstant 1.889727 Bohr %block LatticeVectors * * * * * * * * * %endblock LatticeVectors AtomicCoordinatesFormat ScaledByLatticeVectors %block AtomicCoordinatesAndAtomicSpecies * * * 1 * * * 1 * * * 2 * * * 3 * * * 4 * * * 4 * * * 0 4 * * * 4 * * * 4 * * * 4 %endblock AtomicCoordinatesAndAtomicSpecies %block AtomicCoordinatesOrigin 0.000 0.000 0.000 %endblock AtomicCoordinatesOrigin PAO.BasisType split PAO.SplitNorm 0.15 PAO.EnergyShift 0.01 Ry %block PAO.Basis # Define Basis set * 5 2.23376 n=5 0 1 E 75.06727 3.90515 4.47365 1.00000 n=6 0 2 E 71.19116 3.69937 5.35635 4.45480 1.00000 1.00000 n=5 1 1 E 9.59098 0.01349 4.68524 1.00000 n=6 1 1 E 87.28316 3.65766 5.49527 1.00000 n=5 2 1 E 14.20169 4.25960 5.28396 1.00000 * 5 2.25014 n=3 0 1 E 0.64399 2.65247 3.03150 1.00000 n=4 0 2 E 18.52214 3.13305 4.15184 3.42190 1.00000 1.00000 n=3 1 1 E 8.12894 2.27894 3.02759 1.00000 n=4 1 1 E 16.26911 2.00889 4.26349 1.00000 n=3 2 2 E 14.23042 3.38062 3.78964 2.62819 1.00000 1.00000 * 3 0.09 n=2 0 2 E 24.13 2.59 4.601 2.082 1.000 1.000 n=2 1 2 E 47.15 2.61 5.813 2.127 1.000 1.000 n=3 2 1 E 41.49 4.31 7.798 1.000 %endblock PAO.Basis %block ProjectedDensityOfStates -20.00 10.00 0.050 2000 eV %endblock ProjectedDensityOfStates MeshCutoff 750. Ry FilterCutoff 750. Ry ElectronicTemperature 0.001 eV DM.MixingWeight 0.1 DM.NumberPulay 4 MaxSCFiterations 400 SCFMustConverge .true. WriteMullikenPop 1 COOP.Write .true. SpinPolarized yes NonCollinearSpin no %block kgrid_Monkhorst_Pack 6 0 0 0.5 0 6 0 0.5 0 0 6 0.5 %endblock kgrid_Monkhorst_Pack DM.UseSaveDM yes %block DM.InitSpin * +* 0. 0. # atom,spin,theta,phi (deg) * +* 0. 0. %endblock DM.InitSpin 之前的计算中并无问题，请问怎么回事呢？谢谢 |
回复此楼» 猜你喜欢
机械工程264学硕求调剂
已经有3人回复
-
264求调剂
已经有9人回复
-
求调剂
已经有6人回复
-
22408 266求调剂
已经有7人回复
-
求调剂
已经有14人回复
-
273求调剂
已经有43人回复
-
307求调剂
已经有13人回复
-
327求调剂
已经有7人回复
-
338求调剂
已经有6人回复
-
326分，一志愿沪9，求生物学调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教各位前辈，我这样的条件可以申请推免吗？有没有老师会要啊？谢谢~~~
已经有11人回复
- 【求助】有关标准曲线的疑问，请大家帮帮忙，谢谢！
已经有4人回复
- 【请教】簇合物晶体
已经有14人回复
- 【请教】问个空心金属球的TEM表征问题
已经有21人回复
6楼2011-09-06 20:19:36
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
2楼2011-09-06 19:10:07
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1973.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
3楼2011-09-06 19:13:34
4楼2011-09-06 19:20:23
