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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】请教PWSCF晶体建模和自洽计算，不收敛
222/3返回列表上一页123下一页
查看: 4162  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭