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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】运行出错
111/2返回列表12下一页
查看: 1443  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 伙伴们，祝我生日快乐吧 +17 myrtle 2026-03-10 26/1300 2026-03-16 18:32 by 青橙Ln
[基金申请] 今年的国基金是打分制吗？ 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 0703化学调剂 +6 妮妮ninicgb 2026-03-15 9/450 2026-03-16 16:40 by houyaoxu
[考研] 321求调剂 +5 大米饭！ 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 290求调剂 +5 孔志浩 2026-03-12 10/500 2026-03-16 09:01 by 余晖&
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 289求调剂 +5 步川酷紫123 2026-03-11 5/250 2026-03-15 00:45 by kruisytel
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 327求调剂 +4 Ffff03 2026-03-10 4/200 2026-03-14 00:17 by JourneyLucky
[考研] 【考研调剂求收留】 +3 Ceciilia 2026-03-11 3/150 2026-03-13 20:18 by JourneyLucky
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 328化工专硕求调剂 +4 。，。，。，。i 2026-03-12 4/200 2026-03-13 14:44 by JourneyLucky
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考研] 420求调剂 +4 莫向外求11 2026-03-10 6/300 2026-03-12 14:41 by ruiyingmiao
[考研] 279求调剂 +3 莫xiao 2026-03-10 4/200 2026-03-11 08:06 by 斩魂滴兔子！
[考研] 调剂 +5 呵唔哦豁 2026-03-10 5/250 2026-03-10 22:00 by 28375m
信息提示
关闭
请填处理意见
关闭