| 查看: 1193 | 回复: 2 | |||
xuebiao金虫 (小有名气)
沧州金丝小枣专卖
|
[交流]
【求助】求助一下：关于dmol3的收敛问题 已有1人参与
|
|
最近在学用dmol3，老师给了一个题目：优化Co13，并求得其磁矩和自旋密度图。 然后偶就查文献，然后在一篇文章里找到了Co13的基本构型，此文章里也给出了Co13的平均磁矩是2.38，平均键长也给出了：2.499 （此文献是用siesta优化的） 经过计算可知此Co13的总磁矩是31，多重度是32 有了这些资料那在我自己的计算中就可以简单一些了 Co13是一个二十面体的外壳，中心处还有一个Co原子 用dmol3计算： 主要更改： 点选：spin unrestricted项 并在input里面自己加了:Spin 31 Basis Set  ND Core treatment:All electron然后还点选了：Use smearing （default） 然后进行优化，很快就收敛了。 但是这样的话加了温度结果不是很正确。 所以我一步一步的减少温度 从thermal 0.005 一直减到thermal 0.0002 都能收敛并达到fine的优化结果 但是当我想要继续降温，并到fermi时 ，确实迟迟不能收敛。 所以想问下高手，有什么解决办法吗？ [ Last edited by xuebiao on 2010-7-31 at 19:19 ] |
回复此楼» 猜你喜欢
梦想
已经有7人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
面上合作单位盖章
已经有8人回复
-
面上函评意见出来了，像什么等级？
已经有15人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- dmol3吸附优化不收敛
已经有3人回复
- 关于dmol3的收敛性测试
已经有12人回复
- DMol3计算频率不收敛
已经有4人回复
- 【求助】关于DMol3 收敛阈值
已经有4人回复
- 【求助】关于MS5.0的计算设置问题
已经有12人回复
- 【求助】有人计算过铁卟啉吗？
已经有9人回复
xuebiao
金虫 (小有名气)
沧州金丝小枣专卖
- 应助: 0 (幼儿园)
- 金币: 1066.6
- 散金: 189
- 红花: 1
- 帖子: 203
- 在线: 153.6小时
- 虫号: 719681
- 注册: 2009-03-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
下面是inputfile，在这种情况下是能够收敛的
|
# Task parameters Calculate optimize Opt_energy_convergence 1.0000e-005 Opt_gradient_convergence 2.0000e-003 A Opt_displacement_convergence 5.0000e-003 A Opt_iterations 500 Opt_max_displacement 0.3000 A Symmetry on # Electronic parameters Spin_polarization unrestricted Charge 0 Spin 31 Basis dnp Pseudopotential none Functional pbe Aux_density octupole Integration_grid fine Occupation thermal 0.0002Cutoff_Global 4.5000 angstrom Scf_density_convergence 1.0000e-006 Scf_charge_mixing 0.2000 Scf_spin_mixing 0.5000 Scf_iterations 1000 Scf_diis 6 pulay # Print options Print eigval_last_it # Calculated properties Plot homo Plot lumo Plot density Mulliken_analysis charge Hirshfeld_analysis charge Grid msbox 3 0.2500 0.2500 0.2500 3.0000 |
2楼2010-07-31 19:40:15
xuebiao
金虫 (小有名气)
沧州金丝小枣专卖
- 应助: 0 (幼儿园)
- 金币: 1066.6
- 散金: 189
- 红花: 1
- 帖子: 203
- 在线: 153.6小时
- 虫号: 719681
- 注册: 2009-03-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2010-07-31 19:41:58
20