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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】请教PWSCF晶体建模和自洽计算，不收敛
222/3返回列表上一页123下一页
查看: 3611  |  回复: 21
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 identation 的 49 个金币 ，回帖就立即获得 2 个金币，每人有 1 次机会

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下一页
提示: 如果您在30分钟内回复过其他散金贴，则可能无法领取此贴金币
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 回家两天，不想呆了 +7 368ghnf 2024-06-22 13/650 2024-06-26 19:07 by ZZZemmm
[基金申请] 国基在研影响申请结果吗 +10 WOWO159357 2024-06-26 10/500 2024-06-26 18:19 by 德尚中行
[有机交流] 大佬们，打的核磁氢谱与chemdraw预测的有偏差 +5 来了个 2024-06-26 5/250 2024-06-26 17:28 by Jasonlin7758
[硕博家园] 联培博士文章第一单位署名问题交流~ +9 橙成成c 2024-06-23 21/1050 2024-06-26 16:45 by 橙成成c
[论文投稿] 投稿求助 +5 平凡的日子 2024-06-19 6/300 2024-06-26 16:38 by 大大熊呀
[找工作] 老博士选高校 5+5 woodhead8386 2024-06-24 8/400 2024-06-26 16:35 by 大大熊呀
[论文投稿] 我在写关于多维空间的论文，希望能与大家交流，准备在science上发表 +8 lgf519 2024-06-24 10/500 2024-06-26 16:25 by 梦渺岚烟
[硕博家园] 怎么带研究生？ +24 豆豆小小2 2024-06-20 42/2100 2024-06-26 16:13 by wangyu1986w
[有机交流] 过柱子，产品，杂质在是 pe:ea=100:1 也一起出来？ +5 w256 2024-06-25 5/250 2024-06-26 09:31 by 小木木cc
[硕博家园] 申博雅思求助 +6 刘嘉研究究 2024-06-24 7/350 2024-06-26 07:08 by 小李醒yy
[考博] 申博好难 +6 自强不息a?a 2024-06-24 6/300 2024-06-25 23:02 by 考研吒儿
[硕博家园] 数据不好 +5 Hetai 2024-06-23 7/350 2024-06-25 12:37 by 1591099
[考博] 申请25博士，可以提前进组做科研助理 +6 逐梦途中w 2024-06-22 8/400 2024-06-24 21:49 by keyaner23
[有机交流] 求助析晶问题 20+4 dengdawang 2024-06-24 5/250 2024-06-24 21:22 by cc116
[金属] EBSD的解析率只有10% +3 wallace6666 2024-06-20 7/350 2024-06-24 16:52 by wallace6666
[金属] 寻找钛合金热压缩代做 +4 liuyang358 2024-06-23 4/200 2024-06-24 13:50 by 搬砖狗不放弃
[催化] 求助原位红外 +3 jfdhj 2024-06-21 7/350 2024-06-24 10:52 by qvhm2609
[基金申请] 请问三个字的中文名字，在文中到底怎么表示？ 5+3 zhongyuan764 2024-06-22 6/300 2024-06-23 17:42 by 鱼翔浅底1
[论文投稿] OSA期刊审稿逾期 +3 Thomas_Squid 2024-06-22 3/150 2024-06-23 15:20 by wspglt
[基金申请] 工材口青年基金上会可能性 +10 今晚推荐22 2024-06-19 13/650 2024-06-23 15:11 by Pickfoot
信息提示
关闭
请填处理意见
关闭