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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】请教PWSCF晶体建模和自洽计算，不收敛
222/3返回列表上一页123下一页
查看: 4409  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 工科材料085601 279求调剂 +6 困于星晨 2026-03-17 8/400 2026-03-20 11:24 by kkcoco25
[考研] 274求调剂 +7 S.H1 2026-03-18 7/350 2026-03-20 11:11 by Delta2012
[考研] 286分人工智能专业请求调剂愿意跨考！ +3 lemonzzn 2026-03-17 4/200 2026-03-20 11:04 by lemonzzn
[考研] 08工学调剂 +4 用户573181 2026-03-20 4/200 2026-03-20 10:32 by 朗月清风1
[考研] 298-一志愿中国农业大学-求调剂 +8 手机用户 2026-03-17 8/400 2026-03-20 09:32 by 陆伶越
[考研] 271材料工程求调剂 +7 .6lL 2026-03-18 7/350 2026-03-20 09:10 by xingguangj
[考研] 材料专硕英一数二306 +6 z1z2z3879 2026-03-18 6/300 2026-03-20 08:49 by xingguangj
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 324分 085600材料化工求调剂 +3 llllkkkhh 2026-03-18 3/150 2026-03-19 14:22 by houyaoxu
[考研] 328求调剂，英语六级551，有科研经历 +4 生物工程调剂 2026-03-16 12/600 2026-03-19 11:10 by 生物工程调剂
[考研] 311求调剂 +4 冬十三 2026-03-18 4/200 2026-03-18 21:47 by 尽舜尧1
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 08工科 320总分 求调剂 +5 梨花珞晚风 2026-03-17 5/250 2026-03-18 14:49 by haxia
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
信息提示
关闭
请填处理意见
关闭