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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】PW里结构优化问题
91/1返回列表
查看: 793  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

imunbg

铜虫 (小有名气)

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

我用PW在一个压力点进行结构优化时，发现没到目标压力任务就结束了，对dt 调整以后还是这样，不知何原因，是不是结构有问题？希望得到指点，非常感谢！
回复此楼

» 猜你喜欢
学习 学习 再学习
1楼 2009-11-23 16:30:21
已阅   回复此楼   关注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 ]
赞一下(7人) 回复此楼
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人) 回复此楼
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的全电子是最为合适的.
赞一下(7人) 回复此楼
9楼2009-11-25 18:36:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 imunbg 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 9zw1gfyi7x 2026-08-01 3/150 2026-08-01 18:51 by Kddn1pB6pAzX
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 9zw1gfyi7x 2026-08-01 7/350 2026-08-01 18:43 by Kddn1pB6pAzX
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 9zw1gfyi7x 2026-08-01 6/300 2026-08-01 18:37 by Kddn1pB6pAzX
[考博] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 9zw1gfyi7x 2026-08-01 8/400 2026-08-01 18:31 by Kddn1pB6pAzX
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 9zw1gfyi7x 2026-08-01 8/400 2026-08-01 18:31 by Kddn1pB6pAzX
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 9zw1gfyi7x 2026-08-01 6/300 2026-08-01 18:31 by Kddn1pB6pAzX
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 9zw1gfyi7x 2026-07-31 10/500 2026-08-01 18:20 by Kddn1pB6pAzX
[基金申请] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 11/550 2026-08-01 18:11 by Kddn1pB6pAzX
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 9/450 2026-08-01 18:10 by Kddn1pB6pAzX
[博后之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 9zw1gfyi7x 2026-07-31 9/450 2026-08-01 18:02 by Kddn1pB6pAzX
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 jzg3xk1po8 2026-07-31 12/600 2026-08-01 17:43 by Kddn1pB6pAzX
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 jzg3xk1po8 2026-07-31 9/450 2026-08-01 17:40 by Kddn1pB6pAzX
[基金申请] 2026年国自然面上资助率 +15 布布和一二 2026-07-30 19/950 2026-08-01 15:18 by 40862026
[基金申请] 微信指数没变化，科研之友没阅读 +14 wangze12014 2026-07-28 18/900 2026-07-31 18:35 by lzchaihuo
[高分子] HXDI做水性聚氨酯乳液，是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了，估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[基金申请] 今年的WR进展到哪一步了？ +5 wsgjhwz 2026-07-25 9/450 2026-07-30 12:59 by 可淡不可忘
[基金申请] 时间戳他又来了 +14 晓晓爱翠翠 2026-07-26 17/850 2026-07-29 16:30 by 呵呵二哥
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
信息提示
关闭
请填处理意见
关闭