应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 【求助】越优化越离谱
91/1返回列表
查看: 1061  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

freshgirl

木虫 (正式写手)

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

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

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

有没有人碰到类似的情况？
回复此楼

» 猜你喜欢
1楼 2009-09-29 20:30:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

redskywei

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

linggang87

银虫 (小有名气)
★
wuchenwf(金币+1,VIP+0):谢谢 10-1 22:37
我用VASP计算时也出现楼主的类似情况，一开始力越来越大，不过后来还是收敛了。原因可能是：初始结构偏离基态结构太远，建议修改下结构，或者增大优化时的步长试试。
赞一下(10人) 回复此楼
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的
赞一下(11人) 回复此楼
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等其他事项)
赞一下(10人) 回复此楼
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也采用非常有效和稳定的算
法来精确地计算原子所受的力并优化 ...

赞一下(10人) 回复此楼
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，呵呵～
赞一下(10人) 回复此楼
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优化还算靠谱。不过目前还没人仔细确认过是不是靠谱，只是看上去貌似靠谱。
赞一下(11人) 回复此楼
http://sim.phys.ruc.edu.cn
9楼2009-10-06 22:51:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 freshgirl 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上提前没消息，有中的吗 +10 archvillain 2026-08-02 12/600 2026-08-02 18:46 by hitlend
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 AS6fX1tr1on4 2026-08-02 8/400 2026-08-02 18:44 by ca3DVAx4OB19
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 6/300 2026-08-02 18:34 by ca3DVAx4OB19
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 AS6fX1tr1on4 2026-08-02 8/400 2026-08-02 18:24 by ca3DVAx4OB19
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 AS6fX1tr1on4 2026-08-02 10/500 2026-08-02 18:18 by ca3DVAx4OB19
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 AS6fX1tr1on4 2026-08-02 13/650 2026-08-02 18:14 by ca3DVAx4OB19
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 cBGhfJFuNMca 2026-08-02 16/800 2026-08-02 18:12 by ca3DVAx4OB19
[博后之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 cBGhfJFuNMca 2026-08-02 16/800 2026-08-02 18:11 by ca3DVAx4OB19
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 cBGhfJFuNMca 2026-08-02 20/1000 2026-08-02 18:04 by ca3DVAx4OB19
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 15/750 2026-08-02 16:23 by jHcrHMhbKCFV
[找工作] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 cBGhfJFuNMca 2026-08-02 12/600 2026-08-02 13:43 by jHcrHMhbKCFV
[教师之家] 基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 10/500 2026-08-02 10:00 by 锦衣卫寒战
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[有机交流] 一个有机合成实验室都需要哪些设备？ 50+3 kf2781974 2026-07-31 7/350 2026-08-02 07:59 by 88817753
[基金申请] 微信指数没变化，科研之友没阅读 +14 wangze12014 2026-07-28 18/900 2026-07-31 18:35 by lzchaihuo
[高分子] HXDI做水性聚氨酯乳液，是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了，估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 同事接到电话了，我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
信息提示
关闭
请填处理意见
关闭