| 查看: 2527 | 回复: 4 | |||
[求助]
SIESTA几何优化不收敛是不是就说明了这个结构是不稳定的?
|
|
RT:算得是一个BN sheet 掺一个Ti原子 后 再吸附一个氢气分子的模型 算了两次 没有考虑自旋没有收敛 老师让加入自旋再算一次 还是没有收敛 总分子数:55 个 Ti原子位于BN sheet的 H 位 上脚本: SystemName BNSheetGGA # Descriptive name of the system SystemLabel 1H2TiH # Short name for naming files NumberOfAtoms 53 NumberOfSpecies 4 %block ChemicalSpeciesLabel 1 5 B 2 7 N 3 22 Ti 4 1 H %endblock ChemicalSpeciesLabel XC.functional GGA XC.authors PBE #---------------------------------------------------------------- LatticeConstant 1 Ang # Lattice constant alat %block LatticeVectors 12.55 0.0 0.0 0.0 12.55 0.00 0.0 0.0 15.00 # c=15 A for the vacuum space %endblock LatticeVectors #--------------------------------------------------------------- #--------------------------------------------------------------- %block kgrid_Monkhorst_Pack 1 0 0 0.0 0 1 0 0.0 0 0 1 0.0 %endblock kgrid_Monkhorst_Pack #----------------------------------------------------------------------- #----------------------------------------------------------------------- # SCF options MaxSCFIterations 100 # Maximum number of SCF iter DM.NumberPulay 5 # One Pulay every 5 iter DM.MixingWeight 0.25 # New DM amount for next SCF cycle DM.Tolerance 1.d-4 # Tolerance in maximum difference #DM.UseSaveDM true # to use continuae..ion files #SolutionMethod diagon # OrderN or Diagon MeshCutoff 150. Ry # This data is powered by MS OccupationFunction FD ElectronicTemperature 250 meV # Temp. for Fermi smearing # MD MD.TypeOfRun broyden MD.NumCGsteps 16 #Please Pay Enough ATTENTION on this!!! MD.MaxForceTol 0.05 eV/Ang #Spined SpinPolarized T #------------------------------------------------------------------------ #OUTPUT OPTION LongOutput T AtomicCoordinatesFormat Ang %block AtomicCoordinatesAndAtomicSpecies <1H2TiH.xyz |
回复此楼» 猜你喜欢
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
梦想
已经有7人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 几何优化 收敛问题
已经有6人回复
- 求助：vasp结构优化不能收敛
已经有18人回复
- siesta 表面不收敛
已经有3人回复
- 【求助】结构优化不收敛
已经有9人回复
- 【求助完毕】vasp-slab-边界附近原子顶位吸附小分子的结构优化很难收敛吗？
已经有5人回复
- 【求助】还是收敛测试和几何优化的问题，又来麻烦大家了...
已经有3人回复
- 【求助】结构优化不收敛或是收敛太慢，怎么办？
已经有8人回复
- 【求助】求助siesta的Ordern方法不能收敛的问题
已经有4人回复
- 【求助】siesta中transiesta模块的收敛标准是啥~~
已经有20人回复
- 【求助】siesta中关于几何结构优化的问题【完结】
已经有17人回复
- 【求助】ATK优化计算不收敛怎么办
已经有15人回复
- 【求助】castep几何优化收敛的问题
已经有11人回复
ifmc1234
木虫 (著名写手)
- 应助: 18 (小学生)
- 金币: 2797.9
- 散金: 6001
- 红花: 14
- 帖子: 1761
- 在线: 572.6小时
- 虫号: 825500
- 注册: 2009-08-10
- 专业: 理论和计算化学
【答案】应助回帖
★ ★
youzhizhe(金币+2): 谢谢交流。 2011-05-02 19:59:46
xueke227(金币+1): 2011-05-03 14:39:43
youzhizhe(金币+2): 谢谢交流。 2011-05-02 19:59:46
xueke227(金币+1): 2011-05-03 14:39:43
|
首先不知道你这个体系是否有周期性，假如有，在这块 %block kgrid_Monkhorst_Pack 1 0 0 0.0 0 1 0 0.0 0 0 1 0.0 %endblock kgrid_Monkhorst_Pack 增加K点。 第二，对于吸附一个氢气分子的模型 个人觉得除非你用vwd泛函（siesta有development版本），否则考虑的方法会有问题，方法没选对，结果就没多大的意义。 |
2楼2011-05-02 17:11:47
wangsong1016
铜虫 (初入文坛)
- 1ST强帖: 2
- 应助: 0 (幼儿园)
- 金币: 169
- 散金: 4
- 帖子: 45
- 在线: 26.3小时
- 虫号: 408613
- 注册: 2007-06-21
- 性别: GG
- 专业: 物理
3楼2011-05-03 14:06:28
4楼2011-05-03 14:47:02
ifmc1234
木虫 (著名写手)
- 应助: 18 (小学生)
- 金币: 2797.9
- 散金: 6001
- 红花: 14
- 帖子: 1761
- 在线: 572.6小时
- 虫号: 825500
- 注册: 2009-08-10
- 专业: 理论和计算化学
5楼2011-05-03 17:33:01
20
投票:
