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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】运行出错
111/2返回列表12下一页
查看: 1463  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 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
[考研] 298-一志愿中国农业大学-求调剂 +7 手机用户 2026-03-17 7/350 2026-03-18 14:34 by vgtyfty
[考研] 299求调剂 +5 △小透明* 2026-03-17 5/250 2026-03-18 11:49 by 尽舜尧1
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[基金申请] 被我言中：新模板不强调格式了，假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 考研化学学硕调剂，一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 211本，11408一志愿中科院277分，曾在中科院自动化所实习 +6 Losir 2026-03-12 7/350 2026-03-17 12:09 by danranxie
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 283求调剂 +3 听风就是雨； 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 326求调剂 +4 诺贝尔化学奖觊� 2026-03-15 7/350 2026-03-16 17:11 by 诺贝尔化学奖觊�
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 材料与化工 323 英一+数二+物化，一志愿：哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
信息提示
关闭
请填处理意见
关闭