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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】请教PWSCF晶体建模和自洽计算，不收敛
222/3返回列表上一页123下一页
查看: 4458  |  回复: 21
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

identation

金虫 (著名写手)
引用回帖:
Originally posted by enola at 2010-04-08 11:26:29:

楼主，这个结构从ms导出后，不用先做relax吗？
直接计算？

我是先试验一下，看看和abinit的结果相差多少，嘿嘿，abinit貌似计算铁磁态有问题，于是想试试pwscf，打算先看看abinit的和pwscf的顺磁结果是否一致
赞一下 回复此楼
11楼2010-04-08 11:43:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)
引用回帖:
Originally posted by goldenfisher at 2010-04-07 23:27:19:


这是说你的hamiltonian的本征值计算不收敛。这不影响你的结果，只要你最终得到收敛结果就可以。

那请问一下，如果计算中出现这个提示“     Davidson diagonalization with overlap
     c_bands:  1 eigenvalues not converged
     c_bands:  1 eigenvalues not converged
     c_bands:  2 eigenvalues not converged
     c_bands:  1 eigenvalues not converged
     c_bands:  1 eigenvalues not converged
     c_bands:  1 eigenvalues not converged”

那么计算结果还可靠吗？有没有什么方法可以消除这个呢？
另外，计算的时候如果对于半导体选用occupations=“smearing”，这样可以吗？smearing是金属的，但我看到例子里的NiO也是用smearing的；
另外，LSDA计算时，出错，系统提示要输入那个nelup的数目，自旋向上的电子数目怎么确定呢？
赞一下 回复此楼
12楼2010-04-08 12:08:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★ ★ ★ ★ ★
gavinliu7390:谢谢解答！ 2010-04-08 15:24
aylayl08(金币+5):谢谢解答 2010-04-08 18:31
identation(金币+10): 2010-04-11 20:31
引用回帖:
Originally posted by enola at 2010-04-08 11:41:33:

但是我有个问题啊？这个结构由应力的，计算出来的结构会不会有影响呢？
特别是超胞？
每次计算，如果对自己的结构不确定结构有没有到最优，我们有没有必要都先做relax呢？
谢谢

最好是做一下优化。因为程序的处理总会出现一些精度的问题，而且原子数目越多，这种对结构的预测的误差就越大，每个程序、每个参数对应的最优结构只有一个。有时候一点小小的误差也会导致结果的不同。
对于只包含原胞的体系，优化出来的结果差不多，应该做的是体积优化，也就是说找到最优的晶格常数。
赞一下(2人) 回复此楼
13楼2010-04-08 15:10:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
identation(金币+20):谢谢，嘿嘿，多谢指教哈 2010-04-08 15:30
引用回帖:
Originally posted by identation at 2010-04-08 12:08:04:


那请问一下，如果计算中出现这个提示“     Davidson diagonalization with overlap
     c_bands:  1 eigenvalues not converged
     c_bands:  1 eigenvalues not converged
     c_bands:  2 eigenval ...

这个同样不影响，程序会自动重算这个值。
另外，occupations=“smearing”对于半导体也可以用，不影响计算值。
LSDA计算时并不会要求提供自旋向上向下的值啊。只要你设置好了自旋参量，也就是nspin和statrt_magnetization这两个值，应该就没问题了。
我觉得你的体系可以考虑用LDA+U。
赞一下 回复此楼
14楼2010-04-08 15:14:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)
引用回帖:
Originally posted by goldenfisher at 2010-04-08 15:14:57:

这个同样不影响，程序会自动重算这个值。
另外，occupations=“smearing”对于半导体也可以用，不影响计算值。
LSDA计算时并不会要求提供自旋向上向下的值啊。只要你设置好了自旋参量，也就是nspin和statrt_m ...

谢谢，那+U的话在几何优化的时候也要进行的吧，d电子的强关联作用。但是，PW的几何优化过程，如果晶格常数和原子坐标同时优化，计算时间方面考虑，可行否？abinit试过了，根本没法算，一个星期BFGS才计算了2步，那PW这方面的计算怎样呢？我无经验，请教一下，嘿嘿
赞一下 回复此楼
15楼2010-04-08 15:33:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★
identation(金币+5): 2010-04-08 16:38
qli0218(金币+1):谢谢 2010-04-08 21:20
引用回帖:
Originally posted by identation at 2010-04-08 15:33:14:



谢谢，那+U的话在几何优化的时候也要进行的吧，d电子的强关联作用。但是，PW的几何优化过程，如果晶格常数和原子坐标同时优化，计算时间方面考虑，可行否？abinit试过了，根本没法算，一个星期BFGS才计算了 ...

可以啊，无论是加U还是你说的变胞，都可以。前提是你的参数合理。
变胞就是calculation='vc-relax'。设置目标压强为0即可
赞一下(1人) 回复此楼
16楼2010-04-08 16:37:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

y1ding

铁杆木虫 (著名写手)
★
identation(金币+2):谢谢参与
引用回帖:
Originally posted by identation at 2010-04-08 15:33:14:



谢谢，那+U的话在几何优化的时候也要进行的吧，d电子的强关联作用。但是，PW的几何优化过程，如果晶格常数和原子坐标同时优化，计算时间方面考虑，可行否？abinit试过了，根本没法算，一个星期BFGS才计算了 ...

pw和abinit速度差不多的
赞一下(1人) 回复此楼
17楼2010-04-08 17:07:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

identation

金虫 (著名写手)
引用回帖:
Originally posted by goldenfisher at 2010-04-08 16:37:30:


可以啊，无论是加U还是你说的变胞，都可以。前提是你的参数合理。
变胞就是calculation='vc-relax'。设置目标压强为0即可

那加U是不是只要把Mn的Hubbard_U指定为非零即可?其他元素都不设定，这样可否？
赞一下 回复此楼
18楼2010-04-11 10:42:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★ ★
qasd(金币+2):谢谢 2010-04-12 21:49
是的，设定lsda_u=.true.和你需要设定的原子的U值大小即可。
这个模块的计算需要点技术含量阿。
赞一下(1人) 回复此楼
19楼2010-04-11 23:27:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

风轻云淡7

金虫 (小有名气)
★
identation(金币+2):谢谢参与
我也遇到多你这样的错误，依据我的经验，这样作物的最根本原因是你的结构有问题，特别是原子坐标与晶体结构的饿符合程度……
赞一下(1人) 回复此楼
20楼2010-08-29 21:29:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 identation 的主题更新
222/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +3 archer.. 2026-04-09 3/150 2026-04-09 17:47 by vgtyfty
[考研] 297求调剂 +5 Kwgyz 2026-04-09 5/250 2026-04-09 17:18 by wp06
[考研] 337求调剂 +4 Gky09300550， 2026-04-09 4/200 2026-04-09 17:18 by 帕尔马拉特
[考研] 材料工程调剂 +12 小刘同学吖吖 2026-04-06 13/650 2026-04-09 17:07 by luoyongfeng
[考研] 考研调剂 +13 冰冰，，， 2026-04-07 13/650 2026-04-09 17:01 by Lilly_Li
[考研] 22408 270分 +9 sanjin020722 2026-04-08 10/500 2026-04-09 14:11 by 诗与自由
[考研] 334求调剂 +16 Riot2025 2026-04-08 17/850 2026-04-09 09:28 by wdyheheeh
[基金申请] 有爆料，一个青年教师卖房得400万，然后换了一个四青帽子 +6 babu2015 2026-04-08 6/300 2026-04-09 08:14 by kudofaye
[考研] 302分求调剂 +4 凡语祈愿 2026-04-08 5/250 2026-04-08 22:03 by 土木硕士招生
[考研] 277求调剂 +16 倪建设 2026-04-06 16/800 2026-04-08 21:28 by zhouxiaoyu
[考研] 材料调剂 +4 hzhahg 2026-04-06 4/200 2026-04-08 16:20 by luoyongfeng
[考研] 344求调剂 +11 魏子per 2026-04-07 11/550 2026-04-07 23:01 by JourneyLucky
[考研] 调剂 +18 不逢春 2026-04-05 19/950 2026-04-07 22:04 by lijunpoly
[考研] 312求调剂 +18 gtw1 2026-04-06 20/1000 2026-04-07 18:16 by 蓝云思雨
[考研] 一志愿西电085401求调剂 +4 sunw1306 2026-04-07 4/200 2026-04-07 16:40 by 啵啵啵0119
[考研] 材料调剂 +13 一样YWY 2026-04-03 14/700 2026-04-05 18:20 by 蓝云思雨
[考研] 数一英一274机械调剂 +5 星陨流霞 2026-04-04 6/300 2026-04-05 11:38 by arrow8852
[考研] 调剂 +8 熊二想上岸 2026-04-04 8/400 2026-04-05 05:27 by houyaoxu
[考研] 320调剂 +4 农业工程与信息� 2026-04-03 4/200 2026-04-03 21:40 by lbsjt
[考研] 一志愿安徽大学0817化学工程与技术，求调剂 +14 我不是只因 2026-04-02 15/750 2026-04-03 09:49 by 蓝云思雨
信息提示
关闭
请填处理意见
关闭