24小时热门版块排行榜    

查看: 1567  |  回复: 10

Gina88

木虫 (正式写手)

[求助] BRION应该选几呢?

我在优化的时候,BRION=0, 1, 2
BRION=2为分子动力学优化,优化结果和BRION=1或者BRION=2可以理解。

但是初始结构一样,我的优化用BRION=1牛顿法和BRION=2共轭梯度法,得到的最终结构却完全不同。我的理解,BRION=1和BRION=2原理及理论应该是一样的,仅仅是算法不同,所以得到的最终结果应该是一致的,顶多精度不同。

可是事实是最终结果完全不一样啊。
大侠们也遇到过这样的情况么,如果你们遇到这个情况,应该怎么选择呢?
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuli8

荣誉版主 (知名作家)

2楼2011-08-31 23:07:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gina88

木虫 (正式写手)

引用回帖:
2楼: Originally posted by wuli8 at 2011-08-31 23:07:22:

那您在优化的时候,BRION=1和BRION=2得到的结构一样吗?
3楼2011-09-01 00:22:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuli8

荣誉版主 (知名作家)

…………

优秀版主优秀版主

【答案】应助回帖


Gina88(金币+1): 2011-09-01 09:51:07
zzy870720z(金币+1): 鼓励交流 2011-09-01 10:51:11
引用回帖:
3楼: Originally posted by Gina88 at 2011-09-01 00:22:23:
那您在优化的时候,BRION=1和BRION=2得到的结构一样吗?

不一样。
磁性也不一样。
IBRION=2的时候有时候不能算
…………
4楼2011-09-01 08:48:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gina88

木虫 (正式写手)

引用回帖:
4楼: Originally posted by wuli8 at 2011-09-01 08:48:52:
不一样。
磁性也不一样。
IBRION=2的时候有时候不能算

他们两个的计算结果为什么不一样呢?
IBRION=2为什么不能算呢?
IBRION=2是优化不能收敛?还是优化后的结果是错误的?
我怎么看见有人说IBRION=2是更常用的啊?
5楼2011-09-01 09:53:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)

【答案】应助回帖

★ ★
Gina88(金币+1): 2011-09-01 17:01:53
youzhizhe(金币+2): 谢谢交流。 2011-09-01 19:07:45
"BRION=1和BRION=2原理及理论应该是一样的"
这个理解本身就是错的
不同的寻找最小值的优化算法不能保证每个都能找到“真正”的最小值,每个算法各有千秋
6楼2011-09-01 12:24:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gina88

木虫 (正式写手)

引用回帖:
6楼: Originally posted by bingmou at 2011-09-01 12:24:44:
"BRION=1和BRION=2原理及理论应该是一样的"
这个理解本身就是错的
不同的寻找最小值的优化算法不能保证每个都能找到“真正”的最小值,每个算法各有千秋

谢谢您了,
我只是觉得这两个只不过是一种算法而已,在物理意义上没有什么区别。
所以想当然的觉得得到的结果应该是一样的,差别顶多来自于精度不同。
那这两个优化算法各有什么千秋呢,我是指具体的?
7楼2011-09-01 17:30:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)

【答案】应助回帖

★ ★
gzqdyouxia(金币+2): 谢谢指教 2011-09-02 15:21:31
Gina88(金币+1): 2011-09-02 15:37:48
引用回帖:
7楼: Originally posted by Gina88 at 2011-09-01 17:30:42:
谢谢您了,
我只是觉得这两个只不过是一种算法而已,在物理意义上没有什么区别。
所以想当然的觉得得到的结果应该是一样的,差别顶多来自于精度不同。
那这两个优化算法各有什么千秋呢,我是指具体的?

是的,他们只是一种算法而已,在物理意义上没有区别,但是效率和效用上就是不同的,就像高斯消元法和那个用行列式计算的算法(名字忘记了),都是可以用来求解方程的,但是前者是后者效率的至少100倍,精算精度也不同
这两种算法的优劣在VASP手册中有,对应这个IBRION的部分
8楼2011-09-02 14:24:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gina88

木虫 (正式写手)

引用回帖:
8楼: Originally posted by bingmou at 2011-09-02 14:24:44:
是的,他们只是一种算法而已,在物理意义上没有区别,但是效率和效用上就是不同的,就像高斯消元法和那个用行列式计算的算法(名字忘记了),都是可以用来求解方程的,但是前者是后者效率的至少100倍,精算精度也 ...

bingmou (站内联系TA)
QUOTE:
Originally posted by cenwanglai at 2010-12-14 09:10:04:


是周期性边界条件。你完全明白了我的问题所在。

我的POTIM调到0.1了。是不是可以再调小一点?
我觉得你可以仔细研究是一下是哪些原子的力无法收敛,以及他们为什么在那里无法收敛。POTIM已经很小了,我觉得没有必要修改,而且VASP在弛豫的过程中是修改POTIM的,如果发现不对会缩小POTIM,如果发现正确会加大POTIM

在百度上看见您回答的别人的帖子。
想弱弱的问一下,上哪里去看每个原子的力呢?怎么看每个原子的力有没有收敛呢?如果没有收敛,应该采取什么措施呢?谢谢了!
9楼2011-09-02 18:24:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bingmou

金虫 (著名写手)

【答案】应助回帖


Gina88(金币+2): 2011-09-05 16:44:56
zzy870720z(金币+1): 鼓励交流 2011-09-05 17:34:57
引用回帖:
9楼: Originally posted by Gina88 at 2011-09-02 18:24:34:
bingmou (站内联系TA)
QUOTE:
Originally posted by cenwanglai at 2010-12-14 09:10:04:


是周期性边界条件。你完全明白了我的问题所在。

我的POTIM调到0.1了。是不是可以再调小一点?
我觉得你可以 ...

在VASP中的OUTCAR文件,最后部分会给出每个原子受到的力
10楼2011-09-05 13:11:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Gina88 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +27 医学老男孩 2026-08-13 58/2900 2026-08-15 15:57 by GZUZM
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +7 Ldrop2023 2026-08-13 7/350 2026-08-15 15:57 by mzhh2000
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 P4JhWrbJGjSe 2026-08-14 5/250 2026-08-15 15:46 by QXtWNz7PI7wZ
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 15:10 by QXtWNz7PI7wZ
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 HFw0lei2R37i 2026-08-14 8/400 2026-08-15 14:41 by QXtWNz7PI7wZ
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +7 majunge000 2026-08-11 9/450 2026-08-15 13:34 by kiz2011
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 k0dTPqJtl0jt 2026-08-14 3/150 2026-08-15 12:33 by KxMI1BYBxWX1
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 7lpolszZVXgi 2026-08-14 6/300 2026-08-15 09:28 by sunzitan
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[论文投稿] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-15 00:16 by 4wMiSEwB6436
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见