| 查看: 2745 | 回复: 5 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
伊点若天金虫 (小有名气)
|
[求助]
Dmol3中参数的设置问题。
|
||
|
在对过渡金属氧化物进行结构优化时，成功优化了其中一个表面，但是对于另外一个表面，总是出现以下错误： Message: SCF not converging. Choose "Use Smearing" on DMol3 SCF panel or set "Occupation Thermal" in the input file You may also need to change spin or use symmetry Resubmit DMol3 Message: DMol3 job failed Error: DMol3 exiting Use Smearing 选项也有勾选，Occupation Thermal值设置为0.005，use symmetry也有勾选，请教下Spin值的调整问题，怎样才能进行合理的参数设置来进行优化？或者是Spin值设定外的其他问题。 （第一个表面spin设定为0.5成功优化，第二个表面spin值也是0.5，出现上述问题，但不确定是否为spin值的原因）。 |
回复此楼» 猜你喜欢
本人考085602 化学工程 专硕
已经有15人回复
-
283求调剂
已经有3人回复
-
274求调剂
已经有5人回复
-
302求调剂
已经有6人回复
-
11408 一志愿西电，277分求调剂
已经有3人回复
-
化学工程321分求调剂
已经有13人回复
-
341求调剂
已经有5人回复
-
材料与化工304求B区调剂
已经有7人回复
-
0703化学336分求调剂
已经有4人回复
-
[复试调剂]西南科技大学国防/材料导师推荐
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- MS软件有机分子三维尺寸计算
已经有13人回复
- 求助origin 中组合图参数设置问题？
已经有9人回复
- MS模拟中能不能进行一些简单的化学变化？
已经有7人回复
- MS计算纳米线收敛问题
已经有6人回复
- 几何优化 收敛问题
已经有6人回复
- 关于使用Dmol3进行过渡态搜索的问题！
已经有18人回复
- 关于MS中通过 Brillouin Zone Path 设置和修改的问题
已经有31人回复
- ms中吸附能的计算
已经有40人回复
- 请教有关DMOL3的设置问题
已经有10人回复
- DMol3计算频率不收敛
已经有4人回复
- 掺杂的周期性碳纳米管， 这样优化设置对吗？
已经有12人回复
- 【求助】关于dmol中部分原子固定的几何优化问题
已经有10人回复
- 【求助】请教画能带的问题,有大红包赠送!
已经有20人回复
- 【求助】请教下 Dmol3 下关于在溶剂里催化问题！
已经有3人回复
- 【求助】DMol3过渡态搜索的相关问题
已经有6人回复
- 【求助】DMol3优化参数的选择
已经有17人回复
vasp001
金虫 (正式写手)
- 应助: 69 (初中生)
- 金币: 1513.7
- 散金: 134
- 红花: 24
- 帖子: 702
- 在线: 326小时
- 虫号: 727265
- 注册: 2009-03-20
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★ ★ ★
franch: 金币+5, 谢谢回帖交流，， 2013-03-25 22:48:44
franch: 金币+5, 谢谢回帖交流，， 2013-03-25 22:48:44
|
密度泛函理论发展到现在，已经能够相当好的解决很多实际的问题。特别是“B3LYP”由于其在有机化学中的巨大成功几乎已经成为DFT的代名词。DFT的核心就是“交换相关泛函”。如果这个东西的精确形式被确定，那么，也许从此以后量子化学家就会失业了。当然，短时间内这种可能性似乎不存在。这里通过一个有趣的现象：“分数占据数”来讨论一些交换相关泛函的一种系统误差。 我们知道，基于波函数的从头算方法是在3N维Hilbert空间中进行的，而DFT则是在3维的Euclid（真实）空间中进行的，在真实空间中，我们可以考察密度的变化对能量的影响。对于真实的N电子体系，我们能够确定其能量。那么，如果是“分数”电子呢？比如，H2+在两个H无穷远的时候，考虑对称性，每个核实际上是带1/2电子的，那么这时的能量如何计算呢？ 利用量子统计的方法，Perdew等人证明，对于电子数N+delta （0 因此，对于精确的密度泛函理论，或者说，精确的交换相关泛函，能量E对电子数N的曲线应为下图1的红色线。很明显，在整数电子处，能量对电子数的一阶微分的不连续的（左，右导数不相等）！这是一个很奇异的性质。也就是说，下面这个关系成立： 并且，由于线性插值性： 其中I和A分别为电离能和电子亲和能。而且，我们发现，E对N是凸的，用化学语言就是说，第一电离能小于第二电离能，第二电离能小于第三电离能，等等。这个结论到目前为止无法从数学上证明，但是在实验和数值计算上没有发现例外。 那么近似的泛函满足这个条件吗？考虑一个模型，假设H2电子数为n（不一定是整数N），则相应的电子密度为： 用这个密度做X-alpha计算，可以得到n belongsto [0,2]的能量曲线，即图1中蓝色线： 图1 很明显，Xα泛函失败了，在N=1处，它的一阶微分是连续的，不存在斜率的跳跃。可以证明，对于包括B3LYP在内的所有现在的泛函，没有一个泛函能满足不连续微分条件！这意味着，对于H2+，有 这是不可能的。这个结果是不符合实际的，这是现代所有交换相关泛函的系统误差之一：它们缺乏在整数电子数处的不可微性。这就是分数占据数（fractional occupation）误差。这使得DFT容易高估离域分子的稳定性。TDDFT常常在计算大共轭体系时出现严重的偏差，原因可能也在于此。 http://prb.aps.org/abstract/PRB/v57/i3/p1427_1 |
5楼2013-03-25 00:42:15
vasp001
金虫 (正式写手)
- 应助: 69 (初中生)
- 金币: 1513.7
- 散金: 134
- 红花: 24
- 帖子: 702
- 在线: 326小时
- 虫号: 727265
- 注册: 2009-03-20
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2013-03-23 13:43:13
伊点若天
金虫 (小有名气)
- 应助: 16 (小学生)
- 金币: 1119.1
- 散金: 1
- 红花: 1
- 帖子: 210
- 在线: 165.6小时
- 虫号: 944459
- 注册: 2010-01-18
- 专业: 理论和计算化学
3楼2013-03-24 21:32:11
vasp001
金虫 (正式写手)
- 应助: 69 (初中生)
- 金币: 1513.7
- 散金: 134
- 红花: 24
- 帖子: 702
- 在线: 326小时
- 虫号: 727265
- 注册: 2009-03-20
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2013-03-25 00:40:58
