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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » BRION应该选几呢？
北京石油化工学院2026年研究生招生接收调剂公告
111/2返回列表12下一页
查看: 1293  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 081700，311，求调剂 +13 冬十三 2026-04-04 14/700 2026-04-05 22:23 by 醉翁wl
[考研] 本科211，293分请求调剂 +8 莲菜就是藕吧 2026-04-03 9/450 2026-04-05 19:12 by 蓝云思雨
[考研] 0703化学321分求调剂 +17 三dd. 2026-03-30 18/900 2026-04-05 18:07 by 蓝云思雨
[考研] 080200学硕，机械工程专业277分，求带走！ +7 瓶子PZ 2026-03-31 7/350 2026-04-05 17:49 by liucky
[考研] 298分 070300求调剂 +15 zwen03 2026-04-02 15/750 2026-04-05 12:52 by Hdyxbekcb
[考研] 323求调剂 +8 李佳乐1 2026-04-04 8/400 2026-04-04 22:26 by hemengdong
[考研] 调剂 +9 19945159693 2026-04-03 10/500 2026-04-04 20:16 by dongzh2009
[考研] 调剂 +4 是可乐不是可乐 2026-04-04 4/200 2026-04-04 19:41 by 唐沐儿
[考研] 材料295 +13 小英11 2026-04-03 14/700 2026-04-04 09:02 by 来看流星雨10
[考研] 400分求调剂 +3 尴尬且挠头 2026-04-04 3/150 2026-04-04 08:41 by jp9609
[考研] 305求调剂 +3 77Qi 2026-04-03 3/150 2026-04-03 23:01 by qzxyhcsy
[考研] 调剂 +5 asdasdassda 2026-04-03 6/300 2026-04-03 20:27 by 岸上的一条鱼
[考研] 336求调剂 +8 kiyy 2026-04-01 8/400 2026-04-03 19:41 by lijunpoly
[考研] 338求调剂 +4 zzz，，r 2026-04-03 4/200 2026-04-03 16:39 by lijunpoly
[考研] 一志愿北京科技材料科学与工程288分，求调剂 +14 是辰啊 2026-04-02 14/700 2026-04-02 21:10 by dongzh2009
[考研] 一志愿北交大材料工程总分358 +8 cs0106 2026-04-01 9/450 2026-04-02 10:36 by 不吃魚的貓
[考研] 270调剂 +7 maxjxbsk 2026-04-02 7/350 2026-04-02 09:50 by yulian1987
[考研] 085410 一志愿211 22408分数359求调剂 +3 123456789qw 2026-03-31 4/200 2026-04-02 00:06 by 义文wang
[考研] 材料与化工（0856）304求B区调剂 +8 邱gl 2026-03-30 16/800 2026-04-01 17:58 by 邱gl
[考研] 土木304求调剂 +5 顶级擦擦 2026-03-31 5/250 2026-04-01 08:15 by fdcxdystjk￥
信息提示
关闭
请填处理意见
关闭