24小时热门版块排行榜    

CyRhmU.jpeg
查看: 885  |  回复: 8
当前主题已经存档。

freshgirl

木虫 (正式写手)

[交流] 【求助】越优化越离谱

最近用ATK算过渡金属氧化物,结构优化的时候碰到两个问题:
1,能量收敛,但是force tolerance越算越大,从第一步的几个eV/A到最后几百个eV/A,导致失败
2,能量都不收敛。越算差越大。。。

这些离谱情况的共同规律是,前几次SCF算出来的都挺好,突然有个转折点,越算越离谱。

有没有人碰到类似的情况?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

redskywei

木虫 (著名写手)


wuchenwf(金币+1,VIP+0):谢谢 10-1 22:37
没碰见过,但是atk的优化是有点不好说,castep优化好的结构在他那里还是应力很大,我相信castep
坠入无边的网~~
2楼2009-09-29 20:38:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

linggang87

银虫 (小有名气)


wuchenwf(金币+1,VIP+0):谢谢 10-1 22:37
我用VASP计算时也出现楼主的类似情况,一开始力越来越大,不过后来还是收敛了。原因可能是:初始结构偏离基态结构太远,建议修改下结构,或者增大优化时的步长试试。
3楼2009-09-29 20:41:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zdhlover

荣誉版主 (职业作家)

小木虫之莲花佛

★ ★ ★ ★
wuchenwf(金币+3,VIP+0):谢谢 10-1 22:38
freshgirl(金币+1,VIP+0):我严重怀疑ATK的优化能力。另外,pia pia是什么意思? 10-2 12:19
晕倒了,你这个原因有很多,能否告知你过渡元素所用的基失以及截断能和k点的选取以及你的是什么过渡元素

btw:用别人博士论文的一句原话:另外ATK也采用非常有效和稳定的算
法来精确地计算原子所受的力并优化体系的几何结构。



另外castep优化100个以上的原子那速度要多慢有多慢,但是ATK那是pia pia的
4楼2009-09-30 09:09:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fuzp

铁杆木虫 (正式写手)

我最近用siesta优化graphene ribbon,stress超大,达到10万千巴,而且优化着原子就散架了。是不是siesta只能优化离平衡结构比较近的体系啊。
5楼2009-10-01 23:12:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zdhlover

荣誉版主 (职业作家)

小木虫之莲花佛


qasd(金币+1,VIP+0):谢谢~ 10-4 16:42
哈哈,piapia的是小沈阳的原话,就是那速度还是不错的,我觉得ATK优化还好了,关键你要设置好一些参数,另外对于几个特殊过渡组元素还要额外修改一些东西(譬如Fe,你截断能至少要取到250Ry到300Ry等其他事项)
6楼2009-10-02 15:25:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wei.ji

木虫 (正式写手)

★ ★
freshgirl(金币+1,VIP+0):你是怎么定义“实空间程序”?为什么说ATK是“实空间程序”?如果不用ATK,还能用什么程序优化双极系统? 10-3 12:55
qasd(金币+1,VIP+0):谢谢~ 10-4 16:42
实空间程序优化本来就不是强项,atk做优化。。。
引用回帖:
Originally posted by zdhlover at 2009-9-30 09:09:
晕倒了,你这个原因有很多,能否告知你过渡元素所用的基失以及截断能和k点的选取以及你的是什么过渡元素

btw:用别人博士论文的一句原话:另外ATK也采用非常有效和稳定的算
法来精确地计算原子所受的力并优化 ...

http://sim.phys.ruc.edu.cn
7楼2009-10-02 23:38:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wei.ji

木虫 (正式写手)

★ ★ ★
freshgirl(金币+1,VIP+0):我还是不明白。所有DFT算法不都是在倒空间计算的吗?因为引入了k矢量,已经k空间采样。你这个“实空间”的说法,有出处吗?另外,你所提到的VASP+ATK,效果怎样? 10-4 11:38
qasd(金币+2,VIP+0):谢谢~ 10-4 16:42
atk前身是transiesta-c,再前身是transiesta,transiesta基本上就是mcdcal和siesta拼出来的。siesta是在实空间算总能的。这个总能不太容易算得准,自然,力就不太容易算得准,于是优化起来真得很费劲。我们最近用vasp做优化,优化好的结构拿来算输运性质,当然我们也不用atk,呵呵~
http://sim.phys.ruc.edu.cn
8楼2009-10-04 05:58:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wei.ji

木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
freshgirl(金币+1,VIP+0):那什么算法是到空间算总能呢?VASP和ATK都是基于DFT,赝势法呀。你这个概念有没有出处? 10-7 11:52
原子轨道基组怎么用到k空间呢?siesta 2000年前后的文章里面写了他们怎么算总能的,可以看看是不是在实空间算的。我们用自己的程序算输运,基组可以自己做,所以貌似用vasp优化还算靠谱。不过目前还没人仔细确认过是不是靠谱,只是看上去貌似靠谱。
http://sim.phys.ruc.edu.cn
9楼2009-10-06 22:51:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 freshgirl 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见