| 查看: 4027 | 回复: 12 | ||||
[求助]
原子受力分析比较 已有2人参与
|
||||
如下图，是在优化和静态计算完得到的OUTCAR里的原子受力情况，每一个原子在x,y,z方向的受力都得到了，但是文献中说：convergence of Hellmann- Feynman residual forces less than 0.01 eV/Å per atom was achieved。要求每一个原子的受力都小于0.01ev/A, 这个是理解为x,y,z方向的力都小于0.01，还是什么？如果没有达到要求，继续优化又该怎么进行呢？我的图片中有些原子受力是不是没有达到要求？ 屏幕快照 2017-05-01 上午8.36.32.png |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP计算资料收集byZF |
» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！
已经有9人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- lammps经验心得-2c0q1u4z1n204
已经有105人回复
- 山东高考化学研究
已经有4人回复
- 分子中电荷变化种种
已经有82人回复
- 静电势为什么有正负值？
已经有102人回复
- 常见物理色彩知识
已经有8人回复
- 关于分子间范德华作用的随想
已经有68人回复
- 【转帖】强烈推荐使用phonopy计算声子谱
已经有44人回复
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3578.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
|
我认为可能是按照任意一个atom上force的norm考虑的，（比如说第i个原子）也就是说|f_i|=sqrt(gx_i^2+gy_i^2+gz+i^2)，手册上没说这件事，我的依据是main.F里面： 。。。。。（省略） !-----forces on ions REAL(q) ,ALLOCATABLE:: EWIFOR(:,  , TIFOR(:, !TIFOR应该是每个原子的力，3xN的矩阵。。。。。。（省略） !----------------------------------------------------------------------- ! set DYN%D2C to forces in cartesian coordinates multiplied by FACT ! FACT is determined from timestep in a way, that a stable timestep ! gives a good trial step !----------------------------------------------------------------------- FACT=0 IF (DYN%ISIF<5) FACT=10*DYN%POTIM*EVTOJ/AMTOKG *1E-10_q LSTOP2=.TRUE. NI=1 DO NT=1,T_INFO%NTYP DO NI=NI,T_INFO%NITYP(NT)+NI-1 。。。。。（省略） IF (SQRT(TIFOR(1,NI)**2+TIFOR(2,NI)**2+TIFOR(3,NI)**2) & !每个原子的|f|作为判断，只要有任意原子的|f|>|EDIFFG|，就有LSTOP2=.FALSE & >ABS(DYN%EDIFFG)) LSTOP2=.FALSE. ENDDO ENDDO 。。。。。（省略） IF (DYN%EDIFFG<0) INFO%LSTOP=LSTOP2 !EDIFFF<0的时候以力收敛为标准，此时由逻辑变量LSTOP2赋值给INFO%LSTOP 。。。。。（省略） ! use forces as stopping criterion if EDIFFG<0 IF (INFO%LSTOP) THEN !当INFO%LSTOP为true的时候，打印“reached required accuracy - stopping” 。。。。（省略） WRITE(TIU6,*) 'reached required accuracy - stopping ', & 'structural energy minimisation' ENDIF 虚线下面应该对应的是整个系统的total force。 对于分步优化，这种情况分两种讨论： 1、相同精度分步优化。原则上，这样的情形和一次做到底把NSW设置大点没什么区别，如果在考虑进入程序还要程序初始化，读取一大堆东西其实不见得合算。 2、先精度低、后精度高。(1)实际有可能精度低的时候算出来的结构不一定可靠，甚至有可能更偏离高精度下的平衡结构，此时反而比较慢；(2)低精度下有可能积分误差、SCF的收敛不足等都会会导致势能面毛躁，可能导致构型点在周围几个虚假的极小点来回跳而结构优化难以收敛。 所以我认为大部分情况这样的做法没有特别大的价值。 |
12楼2017-05-01 19:18:19
tsglss
铁杆木虫 (著名写手)
- 应助: 41 (小学生)
- 金币: 6014.2
- 红花: 19
- 帖子: 1875
- 在线: 940小时
- 虫号: 321641
- 注册: 2007-03-11
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2017-05-01 14:21:12
2楼2017-05-01 11:17:05
tsglss
铁杆木虫 (著名写手)
- 应助: 41 (小学生)
- 金币: 6014.2
- 红花: 19
- 帖子: 1875
- 在线: 940小时
- 虫号: 321641
- 注册: 2007-03-11
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2017-05-01 11:58:26
4楼2017-05-01 13:23:41
6楼2017-05-01 14:55:32
7楼2017-05-01 16:28:45
tsglss
铁杆木虫 (著名写手)
- 应助: 41 (小学生)
- 金币: 6014.2
- 红花: 19
- 帖子: 1875
- 在线: 940小时
- 虫号: 321641
- 注册: 2007-03-11
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
8楼2017-05-01 16:35:26
9楼2017-05-01 18:04:11
tsglss
铁杆木虫 (著名写手)
- 应助: 41 (小学生)
- 金币: 6014.2
- 红花: 19
- 帖子: 1875
- 在线: 940小时
- 虫号: 321641
- 注册: 2007-03-11
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
10楼2017-05-01 18:10:12
