| 查看: 2233 | 回复: 6 | ||||
[交流]
【求助完毕】VASP结构优化时离子步原子是如何移动的？
|
|
这个没理解。 结构优化实际上可以理解为多目标优化过程。对于box不变，EDIFFG<0的过程，结构优化是求能量和受力最小的原子坐标，写成方程是： S(E；F1，F2，...， Fn)=f(r1，r2， ...， rn). 其中E,F分别代表体系总能量和原子受力，n表示原子个数，ri表示第i个原子的位置。 由于F是能量E对坐标r的导数，因此
这个方程似乎不是显式的，至少对于多体体系没有解析解。对于任一R,对应一个离子步，可以用自洽场方法找到对应最低能量，以及计算出此时的受力。 然后调整一下R, R(n=new)=R(o=old)+deltaR R(n)对应下一个离子步，然后再通过自洽场方法得到对应能量和受力。 问题是，在S'(E)=f(R)方程具体形式未知的情况下，如何确定deltaR呢？ deltaR一般与能量对座标的一阶导数，或者二阶导数有关。对于势能曲面表达式未知的情况（似乎所有多粒子体系都是这样的），如何求导数呢？岂不是要连导数都要试探？即先预取一点，然后在该点领域通过自洽场迭代算一圈能量，再求导数？ 最后是力的问题。VASP的OUTCAR中,TATAL-FORCE后面三列对应力的三个分量,其负号表示什么意思？如前所述，原子受力是能量对位置微扰的一阶导数（负号）。那么。例如在X轴正向移动delatRx>0，有下面的情况： Fx=-deltaE/deltaRx !*************************************************************! ! 势能曲面函数应该是显式的。尽管没有准确形式，但是各种交换关联泛 ! ！ 函的提出，加上各意义较明确的能量分量，可以近似给出势能曲面表达式。 ! !*************************************************************! [ Last edited by cenwanglai on 2010-12-5 at 10:00 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP |
» 猜你喜欢
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有20人回复
-
今天务委会开完了，明天出结果吗
已经有22人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有9人回复
-
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
放榜前的不淡定
已经有19人回复
-
某些机构，以效率低为荣，以效率低作为存在感
已经有10人回复
-
2026国自然函评费到账
已经有23人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- VASP结构优化时出现的POTIM问题
已经有4人回复
- vasp 离子步数不够，没有达到精度就自动停了，如何继续计算？
已经有6人回复
- vasp 结构优化问题
已经有14人回复
- 求助：vasp结构优化不能收敛
已经有18人回复
- vasp进行固定某些原子的某些坐标的结构优化
已经有5人回复
- vasp如何优化指定的几个原子
已经有9人回复
- 【活动or资源】vasp晶体结构优化
已经有142人回复
- 请假：vasp如何优化单层原子的体系呢？本人新手，希望能说的明白点！
已经有4人回复
- 【求助】vasp 高压相结构优化
已经有17人回复
- 【求助完毕】vasp-slab-边界附近原子顶位吸附小分子的结构优化很难收敛吗？
已经有5人回复
- 【求助】VASP怎么做晶格参数的结构优化
已经有17人回复
- 【求助】VASP结构优化
已经有7人回复
- 【求助成功】VASP中，如何固定磁矩做结构优化
已经有14人回复
- [分享]CLNEB计算实例-Li在石墨上的迁移.
已经有23人回复
2楼2010-12-03 21:49:29
3楼2010-12-04 00:03:59
4楼2010-12-04 09:31:04
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
cenwanglai(金币+2):呵呵，谢谢回复！ 2010-12-04 14:11:54
zzy870720z(金币+1):谢谢交流 2010-12-04 21:14:29
小木虫(金币+0.5):给个红包，谢谢回帖交流
cenwanglai(金币+2):呵呵，谢谢回复！ 2010-12-04 14:11:54
zzy870720z(金币+1):谢谢交流 2010-12-04 21:14:29
|
http://baike.baidu.com/view/2565822.htm 这里看看共轭梯度法， 其余算法推荐陈正隆的分子模拟 |
5楼2010-12-04 09:45:52
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
cenwanglai(金币+3):周末愉快！ 2010-12-04 14:12:07
zzy870720z(金币+1):谢谢交流 2010-12-04 21:14:36
小木虫(金币+0.5):给个红包，谢谢回帖交流
cenwanglai(金币+3):周末愉快！ 2010-12-04 14:12:07
zzy870720z(金币+1):谢谢交流 2010-12-04 21:14:36
|
嗯，半牛顿算法我这里有个PDF，文库里有，了大家看看吧 http://wenku.baidu.com/view/f45b4a7f5acfa1c7aa00cc13.html |
6楼2010-12-04 09:51:40
7楼2010-12-07 20:51:59
