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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】运行出错
北京石油化工学院2026年研究生招生接收调剂公告
111/2返回列表12下一页
查看: 1480  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 能动调剂326专硕 +4 wan112233 2026-04-04 4/200 2026-04-04 22:47 by yu221
[考研] 315求调剂 +8 欣喜777 2026-04-04 9/450 2026-04-04 22:45 by 理想在实践之中
[考研] 本科211 分数293请求调剂 +4 莲菜就是藕吧 2026-04-01 4/200 2026-04-04 22:32 by hemengdong
[考研] 278求调剂 +14 范婷娜 2026-04-04 15/750 2026-04-04 22:15 by lqwchd
[考研] 282求调剂 +6 aaa车辆 2026-04-02 9/450 2026-04-04 22:15 by arrow8852
[考研] 333求调剂 +12 wfh030413@ 2026-04-03 13/650 2026-04-04 21:02 by jj987
[考研] 调剂 +9 19945159693 2026-04-03 10/500 2026-04-04 20:16 by dongzh2009
[考研] 一志愿安徽某211 0703化学总分339求调剂 +6 晚风不晚 2026-04-04 6/300 2026-04-04 20:11 by dongzh2009
[考研] 考研求调剂 +3 木心想继续深造 2026-04-03 3/150 2026-04-03 21:56 by 啵啵啵0119
[考研] 320求调剂 +5 振—TZ 2026-04-02 5/250 2026-04-03 14:42 by fxue1114
[考研] 复试调剂 +3 bvzz 2026-04-01 3/150 2026-04-03 09:47 by 蓝云思雨
[考研] 0703化学 +7 goldtt 2026-04-02 9/450 2026-04-03 09:36 by 蓝云思雨
[考研] 生物学308分求调剂（一志愿华东师大） +6 相信必会光芒万� 2026-03-31 7/350 2026-04-02 23:16 by JourneyLucky
[考研] 环境工程297分求调剂一志愿杭高院 +15 GENJIOW 2026-03-31 16/800 2026-04-02 17:56 by cyh—315
[考研] 318求调剂，计算材料方向 +10 吸喵有害笙命 2026-04-01 11/550 2026-04-02 16:29 by oooqiao
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-04-02 3/150 2026-04-02 15:06 by cal0306
[考研] 275学硕081000服从调剂到其他专业，保不住本专业了 +7 一只小小水牛 2026-04-02 8/400 2026-04-02 14:23 by alice-2022
[考研] 考研调剂 +12 Amber00 2026-03-31 12/600 2026-04-02 09:04 by sanrepian
[考研] 311求调剂 +10 李芷新1 2026-03-31 10/500 2026-04-01 14:38 by chenqifeng666
[考研] 086000生物与医药 初试274求调剂 +5 小叮当来了 2026-03-30 6/300 2026-04-01 14:25 by mxj1203
信息提示
关闭
请填处理意见
关闭