| 查看: 1349 | 回复: 5 | ||||
[求助]
菜鸟求助关于加压优化的问题。。。
|
|
最近在做关于加压优化的东西， 第一个问题：我在toturial 中没找见具体的加压优化的例子（只有lesson3中有不加压优化的例子），请问有没有具体关于加压优化的例子？有的话在那一节？ 第二个问题：Out 文件中：压力张量Cartesian components of stress tensor (hartree/bohr^3) sigma(1 1)= -4.76514276E-08 sigma(3 2)= 0.00000000E+00 sigma(2 2)= -4.76514276E-08 sigma(3 1)= 0.00000000E+00 sigma(3 3)= -4.76514276E-08 sigma(2 1)= 0.00000000E+00 和梯度收敛标准：At Broyd/MD step 3, gradients are converged : max grad (force/stress) = 4.7651E-06 < tolmxf= 5.0000E-05 ha/bohr (free atoms) 之间为什么呈现sigma(1 1)+22+33整体除以三 再乘以10的2次方的关系？按道理说max grad的单位是 ha/boh，而压力张量单位是hartree/bohr^3。这是不同量纲的东西啊？ 三 请问为什么立方晶系sigma(1 1)=22=33而四方晶系（比如金红石）sigma(1 1)=22不等于33呢，11 22 和33之间的关系式什么呢？如果我压价优化金红石晶胞时应该怎眼设置11 22 33呢？ 最后 跪求高手了。。。。 |
回复此楼» 猜你喜欢
估计是周四
已经有3人回复
-
基金啊基金
已经有3人回复
-
科研孤儿太难了
已经有15人回复
-
让我中一个面上吧！
已经有10人回复
-
言之凿凿，明天官网会通知放榜？
已经有17人回复
-
时间戳又变了
已经有17人回复
-
今天系统多次维护，明天很可能放榜！
已经有11人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
放榜前的不淡定
已经有12人回复
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
fzx2008: 金币+5, 专家考核, 谢谢指教！ 2012-12-04 18:38:56
感谢参与，应助指数 +1
fzx2008: 金币+5, 专家考核, 谢谢指教！ 2012-12-04 18:38:56
|
1，没有加压优化的例子，因为所谓的不加压优化，就是为strtarget = 0的优化，而加压优化，是strtarget为一个finite value的优化，造个例子就属于重复建设了。 2，Hartree是能量单位，可以看作是F*l的结果，Bohr是长度单位，也就是l。于是Hartree/Bohr^3就可以看作是F*l/l^3=F/l^2，正好是压强的单位。 3，理论上也应该相等，但是由于四方晶的c轴与a、b非等效轴，所以难以优化后完全一致，但是该问题显然在对称性极高的立方晶上不存在。其实可以看到，如果优化的收敛条件很高，33方向的压强与11，22的差别很小。至于加压，一般都是加等静压，当然是11=22=33了。 |
» 本帖已获得的红花（最新10朵）
2楼2012-12-04 17:43:08
 |
3楼2012-12-04 18:14:05
送鲜花一朵 |
首先 谢谢了 关于第二个问题：我想问max grad（就是max grad (force/stress) = 4.7651E-06），和 sigma(1 1) 22 33 的关系，这三个是压强，但max grad明显是在和tolmxf作比较，所以应该是力。（11+22+33）除以三等于 -4.76514276E-08，这是压强，但max grad (force/stress) = 4.7651E-06 ，这是力，难道说压强乘以100就是力（二者之间相差100倍），这是为什么呢？ 另外一个问题：在计算金红石0压力优化时，残余压力为： -Cartesian components of stress tensor (GPa) [Pressure= 4.8974E-03 GPa] - sigma(1 1)= -1.70282756E+00 sigma(3 2)= 0.00000000E+00 - sigma(2 2)= -1.70282756E+00 sigma(3 1)= 0.00000000E+00 - sigma(3 3)= 3.39096300E+00 sigma(2 1)= 0.00000000E+00 可以看到虽然整体的 [Pressure= 4.8974E-03 GPa] 相对比较小，但是11 22 和33都超过了1GPa................，这合适吗？（就连我设置1G压力优化时三个分量才设为-1G，而这个是0压力优化！）也就是说我应该看压力张量分量呢还是以整体残余压力为标准？ 最后，我悬赏的10个金币应该怎样支付给你？给你两朵小红花吗？ |
4楼2012-12-05 10:30:33
souledge
专家顾问 (著名写手)
-
专家经验: +174 - 1ST强帖: 19
- 应助: 303 (大学生)
- 金币: 11306.4
- 散金: 1123
- 红花: 108
- 帖子: 2814
- 在线: 916.1小时
- 虫号: 941486
- 注册: 2010-01-12
- 专业: 结构陶瓷
- 管辖: 第一性原理
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
若博特: 金币+5, ★★★很有帮助, 很热心的专家 2012-12-14 09:03:28
若博特: 金币+5, ★★★很有帮助, 很好的写手，谢谢了，希望能及时回复下边的问题。 2012-12-15 16:20:41
若博特: 金币+5, ★★★很有帮助, 很热心的专家 2012-12-14 09:03:28
若博特: 金币+5, ★★★很有帮助, 很好的写手，谢谢了，希望能及时回复下边的问题。 2012-12-15 16:20:41
|
1，根据fzx2008版主的提示，该问题可以由下面几个信息解答： a）At Broyd/MD step 14, gradients are converged : max grad (force/stress) = 5.5774E-06 < tolmxf= 1.0000E-04 ha/bohr (free atoms) 所以max grad表示最大梯度。 b） 所以在Abinit中，直接由stress来预测力。一般系数是100. 综上，该方法是取stress矩阵中最大的值来乘100（当然可以少乘，自己设置strfact即可），然后与tolmxf比较。 2，LZ的优化有点粗糙啊，虽然总Pressure过得去，但是，各个方向的确很大，建议让收敛条件更严苛点，顺便看看截断能之类的参数是否收敛。 3，这个问题应该问版主而不是我…… |
5楼2012-12-05 15:06:45
6楼2012-12-10 17:12:09
40