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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】运行出错
111/2返回列表12下一页
查看: 1853  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

gleerat

木虫 (正式写手)

[交流] 【求助】运行出错

我在运行PWscf时，出现如下对称性错误，请各位大牛们帮帮忙哈，谢谢各位了！谢谢！

Subspace diagonalization in iterative solution of the eigenvalue problem:
     a serial algorithm will be used

     Found symmetry operation: I + (  0.0000 -0.5000 -0.5000)
     This is a supercell, fractional translation are disabled
     warning: symmetry operation #  2 not compatible with FFT grid.
   0  -1   1
   0  -1   0
   1  -1   0
     warning: symmetry operation #  3 not compatible with FFT grid.
  -1   0   0
  -1   0   1
  -1   1   0
     warning: symmetry operation #  4 not compatible with FFT grid.
   0   1  -1
   1   0  -1
   0   0  -1
     warning: symmetry operation #  6 not compatible with FFT grid.
   0   1  -1
   0   1   0
  -1   1   0
     warning: symmetry operation #  7 not compatible with FFT grid.
   1   0   0
   1   0  -1
   1  -1   0
     warning: symmetry operation #  8 not compatible with FFT grid.
   0  -1   1
  -1   0   1
   0   0   1
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
明察、慎思、笃行
1楼 2010-09-17 10:55:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★
163.com(金币+1):谢谢参与 2010-09-17 23:02:10
gleerat(金币+1): 2010-10-07 11:09:10
贴出你的输入文件吧，否则不知道你的问题在哪里。
赞一下(1人) 回复此楼
2楼2010-09-17 21:25:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)
引用回帖:
Originally posted by goldenfisher at 2010-09-17 21:25:31:
贴出你的输入文件吧，否则不知道你的问题在哪里。

&CONTROL

                 calculation = 'scf' ,

                restart_mode = 'from_scratch' ,

                      outdir = './' ,

                  pseudo_dir = './' ,

                      prefix = 'CaFe2As2' ,

                     disk_io = 'default' ,

                   verbosity = 'default' ,

                     tprnfor = .true. ,

/

&SYSTEM

                       ibrav = 10,

                   celldm(1) = 5.5062,

                   celldm(2) = 5.4502

                   celldm(3) = 11.664,

                         nat = 5,

                        ntyp = 3,

                     ecutwfc = 30.0 ,

                     ecutrho = 320.0 ,

                       nosym = .false. ,

                        nbnd = 80,

                 occupations = 'smearing' ,

                     degauss = 0.05 ,

                    smearing = 'marzari-vanderbilt' ,

                       nspin = 1 ,

/

&ELECTRONS

                 mixing_beta = 0.3 ,

                 mixing_ndim = 10,

             diagonalization = 'david' ,

/

ATOMIC_SPECIES

   Ca   40.078      Ca.pbe-nsp-van.UPF

   Fe   55.845      Fe.pbe-nd-rrkjus.UPF

   As   74.92159    As.pbe-n-van.UPF

ATOMIC_POSITIONS crystal

  Ca 0.0000000000000000  0.0000000000000000  0.0000000000000000 1 1 1

  Fe 0.2500000000000000  0.2500000000000001  0.2500000000000001 1 1 1

  Fe 0.2500000000000000 -0.2500000000000001  0.2500000000000001 1 1 1

  As 0.0000000000000000  0.0000000000000000  0.3546994010087879 1 1 1

  As 0.0000000000000000  0.5000000000000000  0.1453005989912119 1 1 1









K_POINTS automatic

  5 5 5   0 0 0
赞一下 回复此楼
明察、慎思、笃行
3楼2010-10-07 11:12:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)
引用回帖:
Originally posted by goldenfisher at 2010-09-17 21:25:31:
贴出你的输入文件吧，否则不知道你的问题在哪里。

赞一下 回复此楼
明察、慎思、笃行
4楼2010-10-07 11:21:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)
引用回帖:
Originally posted by goldenfisher at 2010-09-17 21:25:31:
贴出你的输入文件吧，否则不知道你的问题在哪里。

这是我算的体系的晶格结构，输入文件跟出错都在 上面了，麻烦各位高手指点一下，不胜感激。
PW好像是以指定格子类型后找基元的方式建模的，可能这方面做错了，请各位指点。谢谢。。。
赞一下 回复此楼
明察、慎思、笃行
5楼2010-10-07 11:24:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lono75

铁杆木虫 (著名写手)
★
zxzj05(金币+1):谢谢回帖交流 2010-10-07 22:44:36
gleerat(金币+1): 2010-10-08 10:53:59
用xcrysden看一下你的结构吧,
你的celldm(1)(2)(3)是原子单位吗?
建议
ibrav=14
a=
b=
c=
cosbc=
cosac=
cosab=
这样输入结构不易出错.
上面单位埃
赞一下(1人) 回复此楼
由天命,成人事！
6楼2010-10-07 19:03:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★ ★
zxzj05(金币+2):谢谢回帖交流 2010-10-07 22:44:46
gleerat(金币+1): 2010-10-08 10:54:06
ibrav=10是Orthorhombic face-centered晶体构型。你要确认你的晶体确实是这种构型。我认为你写出来的晶体原子坐标和你画出来的图不一致。你的是正交系，我想用ibrav=8,然后把每个不等价的原子坐标都写出来就可以了。
另外，楼上的同仁，我想pwscf中的单位都是原子单位，celldm(1)=a的单位应该是bohr，celld(2)(3)则是相对celldm(1)的值，也就是celldm(2)=b/a，celldm(3)=c/a。
赞一下(1人) 回复此楼
7楼2010-10-07 22:04:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)
引用回帖:
Originally posted by goldenfisher at 2010-10-07 22:04:07:
ibrav=10是Orthorhombic face-centered晶体构型。你要确认你的晶体确实是这种构型。我认为你写出来的晶体原子坐标和你画出来的图不一致。你的是正交系，我想用ibrav=8,然后把每个不等价的原子坐标都写出来就可以了 ...

您好，我想问一下，按你说的用ibrav=8，那么请问，您觉得不等价原子应该是哪几个呢？我先前也用过ibrav＝8，但是感觉不等价原子是不是太多了，PW好像不会自动去找更小的原胞的吧？在算声子的时候原子太多是不是不好呢？谢谢！
赞一下 回复此楼
明察、慎思、笃行
8楼2010-10-08 11:10:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)
引用回帖:
Originally posted by goldenfisher at 2010-10-07 22:04:07:
ibrav=10是Orthorhombic face-centered晶体构型。你要确认你的晶体确实是这种构型。我认为你写出来的晶体原子坐标和你画出来的图不一致。你的是正交系，我想用ibrav=8,然后把每个不等价的原子坐标都写出来就可以了 ...

warning: symmetry operation #  3 not allowed.   fractional translation:
       0.5000000  0.0000000  0.5000000  in crystal coordinates
     warning: symmetry operation #  4 not allowed.   fractional translation:
       0.0000000  0.5000000  0.5000000  in crystal coordinates
     warning: symmetry operation #  7 not allowed.   fractional translation:
       0.0000000  0.5000000 -0.5000000  in crystal coordinates
     warning: symmetry operation #  8 not allowed.   fractional translation:
       0.5000000  0.0000000 -0.5000000  in crystal coordinates
我改了一下输入文件，出现以上的错误，一下是输入文件，请指教，谢谢你。


&CONTROL

                 calculation = 'scf' ,

                restart_mode = 'from_scratch' ,

                      outdir = './' ,

                  pseudo_dir = './' ,

                      prefix = 'CaFe2As2' ,

                     disk_io = 'default' ,

                   verbosity = 'default' ,

                     tprnfor = .true. ,

/

&SYSTEM

                       ibrav = 8,

                   celldm(1) = 5.5062,

                   celldm(2) = 5.4502

                   celldm(3) = 11.664,

                         nat = 14,

                        ntyp = 3,

                     ecutwfc = 30.0 ,

                     ecutrho = 320.0 ,

                       nosym = .false. ,

                        nbnd = 80,

                 occupations = 'smearing' ,

                     degauss = 0.05 ,

                    smearing = 'marzari-vanderbilt' ,

                       nspin = 1 ,

/

&ELECTRONS

                 mixing_beta = 0.3 ,

                 mixing_ndim = 10,

             diagonalization = 'david' ,

/

ATOMIC_SPECIES

   Ca   40.078      Ca.pbe-nsp-van.UPF

   Fe   55.845      Fe.pbe-nd-rrkjus.UPF

   As   74.92159    As.pbe-n-van.UPF

ATOMIC_POSITIONS crystal

  Ca 0.0000000000000000  0.0000000000000000  0.0000000000000000 1 1 1

  Ca 0.5000000000000000  0.5000000000000000  0.0000000000000000 1 1 1

  Ca 0.5000000000000000  0.0000000000000000  0.5000000000000000 1 1 1

  Ca 0.0000000000000000  0.5000000000000000  0.5000000000000000 1 1 1

  Fe 0.2500000000000000  0.2500000000000001  0.2500000000000001 1 1 1

  Fe 0.2500000000000000  0.2500000000000001  0.7500000000000001 1 1 1

  As 0.0000000000000000  0.0000000000000000  0.3546994010087879 1 1 1

  As 0.5000000000000000  0.5000000000000000  0.3546994010087879 1 1 1

  As 0.0000000000000000  0.5000000000000000  0.1453005989912119 1 1 1

  As 0.5000000000000000  0.0000000000000000  0.1453005989912119 1 1 1

  As 0.0000000000000000  0.0000000000000000  0.6453005990000000 1 1 1

  As 0.5000000000000000  0.5000000000000000  0.6453005990000000 1 1 1

  As 0.0000000000000000  0.5000000000000000  0.8546994010000000 1 1 1

  As 0.5000000000000000  0.0000000000000000  0.8546994010000000 1 1 1

  









K_POINTS automatic

  5 5 5   0 0 0
赞一下 回复此楼
明察、慎思、笃行
9楼2010-10-08 11:56:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)
★
cenwanglai(金币+1):谢谢回复！ 2010-10-10 09:11:48
你的第二个输入文件和第一个得到的晶体不一样
其次，这个问题貌似在于你的晶体构型有点奇怪，你确定a只有3埃不到？你可以修改看看。我没有遇到过这个问题，可以上mailist看看

[ Last edited by goldenfisher on 2010-10-9 at 11:46 ]
赞一下(1人) 回复此楼
10楼2010-10-09 11:44:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gleerat 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +55 医学老男孩 2026-08-13 123/6150 2026-08-19 23:44 by stan_et
[基金申请] 今天基金会出结果吗？20260819 +13 kkkl_v 2026-08-19 13/650 2026-08-19 22:37 by 宝贝虫子
[基金申请] 时间戳变了，能看出什么问题？ +12 基诺咪客 2026-08-17 16/800 2026-08-19 22:04 by kk 的小木虫
[基金申请] filecode，4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低？主要原因有哪些？ +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜吗？ +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +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
[基金申请] 什么时候开奖？ +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 重要消息，中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
信息提示
关闭
请填处理意见
关闭