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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » BRION应该选几呢？
111/2返回列表12下一页
查看: 1172  |  回复: 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下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭