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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » BRION应该选几呢？
111/2返回列表12下一页
查看: 1573  |  回复: 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，我:8O.551.O54,科目全,可十急 +4 HFw0lei2R37i 2026-08-14 6/300 2026-08-15 17:34 by 7huP8NPKZCp9
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +7 Ldrop2023 2026-08-13 7/350 2026-08-15 15:57 by mzhh2000
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 HFw0lei2R37i 2026-08-14 8/400 2026-08-15 14:41 by QXtWNz7PI7wZ
[博后之家] 售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
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 HFw0lei2R37i 2026-08-14 9/450 2026-08-15 11:33 by KxMI1BYBxWX1
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 07:09 by 4wMiSEwB6436
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 04:45 by 4wMiSEwB6436
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 7lpolszZVXgi 2026-08-14 5/250 2026-08-15 04:41 by 4wMiSEwB6436
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[基金申请] 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 要乐观耀哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
关闭
请填处理意见
关闭