24小时热门版块排行榜    

查看: 686  |  回复: 8
当前主题已经存档。

imunbg

铜虫 (小有名气)

[交流] 【求助】PW里结构优化问题

我用PW在一个压力点进行结构优化时,发现没到目标压力任务就结束了,对dt 调整以后还是这样,不知何原因,是不是结构有问题?希望得到指点,非常感谢!
回复此楼
学习 学习 再学习
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

watermall

金虫 (小有名气)

把结果文件后面给的信息给出来啊
你这样问,别人没办法帮你
2楼2009-11-24 13:42:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

imunbg

铜虫 (小有名气)

问题(1):  当压力设为1200 Kbar时 达到如下值时任务完成了
          total   stress  (ryd/bohr**3)                  (kbar)     P= 1210.66
   0.00823982   0.00000000   0.00000000       1212.12      0.00      0.00
   0.00000000   0.00823982   0.00000000          0.00   1212.12      0.00
   0.00000000   0.00000000   0.00821013          0.00      0.00   1207.75

问题(2) 另一个压力点任务中断了,出现了如下提示:         
Davidson diagonalization (with overlap)
请指点,谢谢!
学习 学习 再学习
3楼2009-11-24 18:48:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjmaxpayne(金币+3,VIP+0):谢谢 11-25 08:47
问题(2)是报错信息么?
对于所有提问的人来说,
如果你不是超级幸运,那么很少有人遇到过你现在遇到的问题。
所以,希望您尽量给出详细的信息;)当然不是全部的输出文件;)

求解本征值时,如果交叠矩阵式非正定的,那么将无法进行下三角矩阵分解,自然无法得到本征能量:)
你给出的提示中有 (with overlap),我猜,仅仅是猜,可能与此有关;)

[ Last edited by xirainbow on 2009-11-24 at 21:36 ]
4楼2009-11-24 21:34:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

imunbg

铜虫 (小有名气)

首先再次感谢您的指点!我降低精度计算以后任务是完成了,但最终结果是
         total   stress  (ryd/bohr**3)                  (kbar)     P= 1258.73
   0.00859456   0.00000000  -0.00000001       1264.30      0.00      0.00
   0.00000000   0.00859457   0.00000000          0.00   1264.30      0.00
  -0.00000001   0.00000000   0.00848099          0.00      0.00   1247.60
我本来设的压力值是 1300 kbar,还是没达到,
学习 学习 再学习
5楼2009-11-24 22:59:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

imunbg

铜虫 (小有名气)

首先再次感谢您的指点!我降低精度计算以后任务是完成了,但最终结果是
         total   stress  (ryd/bohr**3)                  (kbar)     P= 1258.73
   0.00859456   0.00000000  -0.00000001       1264.30      0.00      0.00
   0.00000000   0.00859457   0.00000000          0.00   1264.30      0.00
  -0.00000001   0.00000000   0.00848099          0.00      0.00   1247.60
我本来设的压力值是 1300 kbar,还是没达到,
学习 学习 再学习
6楼2009-11-24 23:02:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

watermall

金虫 (小有名气)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjmaxpayne(金币+2,VIP+0):谢谢 11-25 08:47
并不是你想加多大压力就能加多大压力的
计算模拟毕竟是模拟,你要在软件能模拟的范围内模拟
你首先要搞清楚你的体系能否加到1300Kbar,在1300Kbar时候,你所用的赝势是否仍然符合要求
7楼2009-11-24 23:20:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

imunbg

铜虫 (小有名气)

同一个体系在Castep里优化时结果还可以。那么也许PW里的赝势有点问题,谢谢了!
学习 学习 再学习
8楼2009-11-25 08:49:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

y1ding

铁杆木虫 (著名写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
yjmaxpayne(金币+3,VIP+0):谢谢 11-25 21:25
引用回帖:
Originally posted by imunbg at 2009-11-25 08:49:
同一个体系在Castep里优化时结果还可以。那么也许PW里的赝势有点问题,谢谢了!

PW是超软,Castep则是膜守恒,两者赝势是不同的.
在高压计算中,膜守恒赝势比超软有更好的表现.
你的现象正好说明了这一点.
如果有条件的话,高压情况用VASP的PAW,或者Wien2k的全电子是最为合适的.
9楼2009-11-25 18:36:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 imunbg 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见