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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » BRION应该选几呢？
111/2返回列表12下一页
查看: 1566  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

Gina88

木虫 (正式写手)

[求助] BRION应该选几呢？

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

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

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

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-08-31 16:09:14
已阅   回复此楼   关注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的时候有时候不能算
赞一下(1人) 回复此楼
…………
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原理及理论应该是一样的"
这个理解本身就是错的
不同的寻找最小值的优化算法不能保证每个都能找到“真正”的最小值，每个算法各有千秋
赞一下(1人) 回复此楼
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的部分
赞一下(1人) 回复此楼
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文件，最后部分会给出每个原子受到的力
赞一下(1人) 回复此楼
10楼2011-09-05 13:11:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Gina88 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 HFw0lei2R37i 2026-08-14 8/400 2026-08-15 14:41 by QXtWNz7PI7wZ
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +26 医学老男孩 2026-08-13 57/2850 2026-08-15 14:25 by JChen2024
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 k0dTPqJtl0jt 2026-08-14 3/150 2026-08-15 12:41 by KxMI1BYBxWX1
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 k0dTPqJtl0jt 2026-08-14 3/150 2026-08-15 12:33 by KxMI1BYBxWX1
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-15 11:41 by KxMI1BYBxWX1
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 HFw0lei2R37i 2026-08-14 9/450 2026-08-15 11:33 by KxMI1BYBxWX1
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +6 Ldrop2023 2026-08-13 6/300 2026-08-15 10:09 by tuanggou
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 7lpolszZVXgi 2026-08-14 5/250 2026-08-15 04:41 by 4wMiSEwB6436
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 4/200 2026-08-15 04:17 by 4wMiSEwB6436
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +11 Tide man 2026-08-10 12/600 2026-08-15 02:12 by home3163
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-14 22:32 by 4wMiSEwB6436
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
关闭
请填处理意见
关闭