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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 【求助】优化后成功晶胞达到最佳的标志！
232/3返回列表上一页123下一页
查看: 1904  |  回复: 22
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
本帖产生 1 个 1ST强帖 ，点击这里进行查看

qzhaosdu

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

shelay

木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
没有total energy吧，优化完只给出final enthalpy啊？怎么回事啊？
赞一下(1人) 回复此楼
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也已经达到了收敛标准。所以说这个优化是没有问题的
赞一下(3人) 回复此楼
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):给个红包，谢谢回帖交流
高手就是猛啊……谢高手！
赞一下(1人) 回复此楼
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 ]
赞一下(2人) 回复此楼
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



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

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

1003307944

银虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
free energy或者total energy之差在1meV至10meV，取最大的那个能量
赞一下(1人) 回复此楼
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计算，除了得到体系 ...

谢谢楼主详细的解释
赞一下(1人) 回复此楼
timeflies..
20楼2010-04-02 15:27:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liang7035 的主题更新
232/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] filecode，4个jtjc了 +11 ziyangfang 2026-08-19 14/700 2026-08-21 16:51 by applepetty
[基金申请] 让我中一个面上吧！ +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 只有每年这种时候来逛逛小木虫 +23 yaoyewhu2008 2026-08-20 25/1250 2026-08-21 12:19 by Poppy1104
[基金申请] 什么时候开奖？ +7 CrisMessi 2026-08-18 7/350 2026-08-21 12:15 by weiyin
[基金申请] 时间戳又变了8-15 +15 archvillain 2026-08-15 27/1350 2026-08-21 11:32 by tim76
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 13/650 2026-08-21 08:30 by 北京莱茵编辑
[基金申请] 时间戳今天，20号变了 +4 archvillain 2026-08-20 4/200 2026-08-21 08:10 by gatelove
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[教师之家] 为什么余额宝的年化利率越来越低？主要原因有哪些？ +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜吗？ +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
信息提示
关闭
请填处理意见
关闭