24小时热门版块排行榜    

查看: 2565  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

xiaoxiao136

木虫 (正式写手)

[交流] 【求助】结构优化不收敛 已有4人参与

我用Gaussian优化好的结构用到adf,但是还是不收敛,该怎么调节呢,大家说一说啊,手册上写的看得有点迷茫
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiechong

金虫 (著名写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
ADF好贵呀, 买不起, 倾家荡产也只能用一年

咋这么贵呢
9楼2011-03-22 10:18:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

九仙居士

捐助贵宾 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
以下内容是我以前在一个坛子里看到的,抄到这里,你看看是否有用:

对于含有单电子的开壳层体系,如过渡金属或者稀土化合物,往往容易遇到SCF 不收敛的问题。对于这类体系,ADF 软件的使用手册(ADFUsersGuide.pdf)针对SCF 不收敛问题提供了以下几种建议:
1, New DIIS
使用NEWDIIS 关键字,在SCF 过程中调用新的DIIS 方法,这可以解决一部分的不收敛问题;
2, Keep Orbitals
对于某些SCF 特别难收敛的体系,如开壳层高自旋体系可以使用关键字OCCUPATIONS 并指定keeporbitals,把轨道顺序固定下来避免能量震荡。这种方法几乎一定能让SCF 收敛,但是轨道的顺序是否正确却无法保证,除非事先指定轨道排布。图形界面下ADFinput 可以先列出各个轨道并分别指定其中的电子数,对于较小的分子体系不妨一试,利用这种方法甚至可以设置出一个激发态的电子轨道占据。
3, Mixing
SCF 关键字中有个选项mix,其意义是根据占据轨道算的当前势的相对权重,默认值为0.2。对于需要强阻尼的困难体系,适当减少mix 的值有助于SCF 收敛,一般设置为0.02-0.05 之间比较合适。
4, Level Shift
SCF 关键字中还有个选项能级移动参数lShift,可以人为地将非占据轨道的能量拔高( 单位:hatree)。在SCF 迭代中,当电荷在能量非常接近且都位于Fermi 能级附近的
不同轨道之间前后摆动时,这可帮助解决收敛问题。不过需要注意,涉及到能级跃迁的计算( 如预测UV 谱) 使用此参数将导致结果出现偏差。
5, Smearing
对于非常难收敛的体系,还可以使用OCCUPATIONS关键字的smearq 选项, 这相当于给体系一个温度并按照Boltzman 分布排电子, 电子排布不再是aufbau 的。Smearing 选项可以让SCF 很快收敛,不过这将会引入分数形式的轨道占据数,对此可以先设定一个smearq 的值,然后逐步降低至零,类似于一个退火的过程。早期版本中这需要依次修改smearq 的大小然后多次restart, 现在直接在"smearq=" 之后填一串数字即可。ADFinput 中的设置更简单,分别指定smearq 的初始值和终值,再给出退火的步数,它将自动生成一系列递减数值并且步长也递减,比手动设置的会更精细一些。对于smearq 的每一个值,SCF 都将作一次自洽,因此SCF 的最大循环圈数需要放宽,一般按每个值50 圈估算。smearq 的初始值如何设置,这取决于具体的体系,需要尝试。
2楼2011-03-15 08:49:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaoxiao136

木虫 (正式写手)

引用回帖:
Originally posted by 九仙居士 at 2011-03-15 08:49:50:
以下内容是我以前在一个坛子里看到的,抄到这里,你看看是否有用:

对于含有单电子的开壳层体系,如过渡金属或者稀土化合物,往往容易遇到SCF 不收敛的问题。对于这类体系,ADF 软件的使用手册(ADFUsersGuide ...

好的,我会试试的,谢谢啦
3楼2011-03-15 15:39:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

levyleo

金虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
什么东西这么诡异!
通常的C、H、O、N体系不至于这样啊!
我太渺小了。
4楼2011-03-18 09:16:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见