24小时热门版块排行榜    

查看: 1454  |  回复: 22
当前主题已经存档。
本帖产生 1 个 1ST强帖 ,点击这里进行查看

qzhaosdu

金虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
qli0218(金币+1):谢谢 2010-04-01 21:02
应该是你拿实验上测出来的晶胞做一下几何优化,如果收敛了并且和实验上的几何参数相差不大就行了吧,不过这种可能性很小~~
Anewday,anewhour,anewminute,anewpeople.
11楼2010-04-01 19:26:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shelay

木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
没有total energy吧,优化完只给出final enthalpy啊?怎么回事啊?
12楼2010-04-01 19:29:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
qli0218(金币+3):非常感谢! 2010-04-02 08:53
gavinliu7390(金币+2, 1ST强帖+1):回答的很好,很形象! 2010-04-02 19:31
引用回帖:
Originally posted by enola at 2010-04-01 19:12:42:
我看到ms的优化结果里会出现total energy,max stress,max place。。。。
哪位大侠可以讲一下吗?怎么从这里得到优化的进度,或者收敛没?

看来有同学对优化的基本步骤不是很清楚啊,那我给大家稍微科普一下吧。只是我个人的见解而已,更深的理论知识还是看文献吧
几何优化的基本过程就是一系列的scf单点能量计算。对于每步scf计算,除了得到体系总能量之外,还得到了能量的梯度方向。之后根据梯度找能量更低的构形,再重复新的scf计算,直到找到梯度为0的构形。但这个不一定是最小值点。形象点说就象在山上迷路了,想下山,那么只能根据当前所在位置的坡度方向往更低的地方走,直到遇到坡度为0的位置,但这里不一定是海拔最低的山脚下,也有可能是卡在高海拔的山凹中。所以找的是极值点。
在计算中,能量梯度显得有些抽象,那么我把它形象一些:体系总能量对某原子坐标的梯度可以看成该原子的受力(F=dE/dr)。因此几何优化的过程实际就是:根据能量梯度(原子受力)方向搜索能量更低的状态,直到找到极小值点,梯度(原子受力)为0的构形。
优化晶格也是类似的,只是这里的梯度关系是:应变能vs应力。

到这里我们就可以总结total energy,max stress,max displacement, max force这些东西了。
total energy就是总能量,优化的目的就是找它的极值点。在接近极值点,即优化快结束时,能量的起伏已经很小,它接近一条水平线。
max force/stress 就是体系残余的最大力/应力,也就是当前最大的梯度,在接近极值点时,它们接近消失
max displacement 是最大位移,简单点说就是知道某个受力方向了,为了降低能量,需要沿着这个方向移动原子,这个就是displacement,顾名思义,max displacement就是最大的移动幅度。显然,在接近极小值点时,这个值也应该趋于消失的。

以楼主的情况为例,优化的是镍的晶胞,原子是处在正确的点阵位置,原子受力是0,那么位移自然也是0,所以这两条线是没有变化的。而能量变化主要是应变能的贡献,而到了第5步,max stress也已经达到了收敛标准。所以说这个优化是没有问题的
13楼2010-04-01 23:14:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liang7035

金虫 (正式写手)

引用回帖:
Originally posted by acridine at 2010-04-01 23:14:03:


看来有同学对优化的基本步骤不是很清楚啊,那我给大家稍微科普一下吧。只是我个人的见解而已,更深的理论知识还是看文献吧
几何优化的基本过程就是一系列的scf单点能量计算。对于每步scf计算,除了得到体系 ...

谢谢高手啊,真是顿开了。。。谢谢。。。
一个一心向上的人整个世界都会给他让路!
14楼2010-04-02 07:24:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小木虫503

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
高手就是猛啊……谢高手!
15楼2010-04-02 08:55:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

carlaty

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by shelay at 2010-04-01 19:29:50:
没有total energy吧,优化完只给出final enthalpy啊?怎么回事啊?

同问
16楼2010-04-02 09:10:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
qli0218(金币+1):谢谢 2010-04-02 14:48
qli0218(金币+1):谢谢。但是能否能判定是局域极小值还是全局极小值 2010-04-02 08:52



上面已经有人科普过了。
结构优化找的其实就是你所给结构附近的极小值。就好比让一块石头顺着山坡可以滑到山沟里面,而不是世界最低点的***海沟。
要找全局最小值的话,其实是个很麻烦的事情尤其是结构比较复杂的时候。好像有一个软件专门用遗传算法做这个的。

[ Last edited by dxcharlary on 2010-4-2 at 13:46 ]
17楼2010-04-02 13:45:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qli0218

木虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by dxcharlary at 2010-04-02 13:45:38:
qli0218(金币+1):谢谢。但是能否能判定是局域极小值还是全局极小值 2010-04-02 08:52



上面已经有人科普过了。
结构优化找的其实就是你所给结构附近的极小值。就好比让一块石头顺着山坡可以滑到山沟里面, ...

谢谢,那么应该就需要先用动力学计算一下,再来找局域最小值了。
18楼2010-04-02 14:48:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1003307944

银虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
free energy或者total energy之差在1meV至10meV,取最大的那个能量
19楼2010-04-02 15:03:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

enola

捐助贵宾 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by acridine at 2010-04-01 23:14:03:


看来有同学对优化的基本步骤不是很清楚啊,那我给大家稍微科普一下吧。只是我个人的见解而已,更深的理论知识还是看文献吧
几何优化的基本过程就是一系列的scf单点能量计算。对于每步scf计算,除了得到体系 ...

谢谢楼主详细的解释
timeflies..
20楼2010-04-02 15:27:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liang7035 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见