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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】请教PWSCF晶体建模和自洽计算，不收敛
222/3返回列表上一页123下一页
查看: 4839  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 让我中一个面上吧！ +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 人气不行了 +7 fansofjerry 2026-08-21 7/350 2026-08-21 14:13 by Vivilian
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 时间戳又变了8-15 +15 archvillain 2026-08-15 27/1350 2026-08-21 11:32 by tim76
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 时间戳今天，20号变了 +4 archvillain 2026-08-20 4/200 2026-08-21 08:10 by gatelove
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 今天系统多次维护，明天很可能放榜！ +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 时间戳变了，能看出什么问题？ +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天放榜吗？ +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭