24小时热门版块排行榜    

查看: 2670  |  回复: 43

cj4566

木虫 (正式写手)

[求助] 紧急求助:运行pwscf出现对称性操作警告

我在运行pwscf进行结构优化时,出现下面的错误,麻烦大家帮忙指点一下啊。
     Subspace diagonalization in iterative solution of the eigenvalue problem:
     scalapack distributed-memory algorithm (size of sub-group:  2*  2 procs)

     warning: symmetry operation #  3 not compatible with FFT grid.
  -1   0   0
   1   1   0
   0   0  -1
     warning: symmetry operation #  4 not compatible with FFT grid.
   1   0   0
  -1  -1   0
   0   0  -1
     warning: symmetry operation #  7 not compatible with FFT grid.
   1   0   0
  -1  -1   0
   0   0   1
     warning: symmetry operation #  8 not compatible with FFT grid.
  -1   0   0
   1   1   0
   0   0   1

输入文件:
&CONTROL
                 calculation  = 'relax',
                 restart_mode = 'from_scratch',
                 prefix       = "101-2",
                 pseudo_dir   = "/pkg/suse11/quantum/5.0/pseudo",
                 outdir       = "/home/n7388683/QE/tmp",  
                      tstress = .true
                      tprnfor = .true                  
/
&SYSTEM
                       ibrav = 12,
                           a = 3.287,
                           b = 6.237385,
                           c = 23.868112,
                       cosab = -0.263491891704172,
                         nat = 8,
                        ntyp = 2,
                 occupations = 'smearing',
                     degauss = 0.02D0,
                     ecutwfc = 30,
                     ecutrho = 300,
/
&ELECTRONS
                    conv_thr = 1.D-8,
                 mixing_beta = 0.4,
/
&IONS
                ion_dynamics = 'bfgs',           
/
ATOMIC_SPECIES
   Zn   65.3799972534  Zn.pbe-van.UPF
    O   15.9989995956  O.pbe-van_ak.UPF
ATOMIC_POSITIONS crystal
  O   0.2199834158725568   0.3844112496986048   0.0175122228436065
  O   0.8318951398867058   0.6082346977269033   0.1225855599052453
  O   0.5632431858843476   0.0709307897221868   0.1050733370616389
  O   0.9513314618701982   0.8471073416938884   0.0000000000000000
Zn   0.0741775052813360   0.0927994285161632   0.0569885856639280
Zn   0.6860892292954852   0.3166228765444618   0.1620619227255668
Zn   0.8055255512789776   0.5554955205114469   0.0394763628203215
Zn   0.4174372752931270   0.7793189685397455   0.1445496998819603
K_POINTS automatic
  8 8 1   0 0 0

计算的是氧化锌的一个极性面,会不会因为这个面不稳定出现上面情况啊,计算其他晶面的时候都没出现过,谢谢了!!!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
cj4566: 金币+20, ★★★★★最佳答案 2012-09-25 09:58:41
franch: 金币+2, 谢谢回帖交流 2012-09-25 10:03:07
QE5.0改成这样的警告了?挺有意思的
是这样的,你所使用的实空间FFT网格的格点并没有刚好取到平移对称性所在的位置,所以它只能忽略你体系的平移对称性。
解决方法也好办:
1.如果你对对称群没要求,或者不用进行声子模式分析,可忽略
2.结构优化时不用管,真正自洽计算时如果出现该警告再考虑消除。
3.如果想消除本警告,可设置nr1,nr2,nr3的值,这三个值决定了你FFT网格在三个方向上有多少格点。如果不清楚要取多少,请参考输出文件里的“FFT dimensions”一行,然后酌情稍作修改,使得你取的FFT网格中包含平移对称性。
2楼2012-09-25 08:52:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

其实如果嫌麻烦且体系对称性简单,直接改ecutrho就行了。
3楼2012-09-25 08:55:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

可以参考这个http://muchong.com/html/201010/2409072.html

虽然里面没有提出来如何解决,但是如果你能像帖子里那样改变一下晶格类型,将警告
warning: symmetry operation #  3 not compatible with FFT grid.
  -1   0   0
   1   1   0
   0   0  -1
变成如下形式:
warning: symmetry operation #  3 not allowed.   fractional translation:
       a  b  c  in crystal coordinates

出现平移坐标(a,b,c)之后,具体解决方案就是如下:
设置的nr1一定要能整除a, nr2要能整除b, nr3要能整除c,然后就没有警告了。
4楼2012-09-25 09:06:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

★ ★
franch: 金币+2, 谢谢回帖交流 2012-09-25 10:03:22
引用回帖:
4楼: Originally posted by gemucai at 2012-09-25 09:06:52
可以参考这个http://muchong.com/html/201010/2409072.html

虽然里面没有提出来如何解决,但是如果你能像帖子里那样改变一下晶格类型,将警告
warning: symmetry operation #  3 not compatible with FFT grid.
...

最后一句错了,设置的nr1乘以a一定要是整数,nr2乘以b是整数,nr3乘以c是整数而不是分数。
5楼2012-09-25 09:09:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cj4566

木虫 (正式写手)

引用回帖:
2楼: Originally posted by gemucai at 2012-09-25 08:52:10
QE5.0改成这样的警告了?挺有意思的
是这样的,你所使用的实空间FFT网格的格点并没有刚好取到平移对称性所在的位置,所以它只能忽略你体系的平移对称性。
解决方法也好办:
1.如果你对对称群没要求,或者不用进行 ...

非常感谢,这就试试
6楼2012-09-25 09:58:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cj4566

木虫 (正式写手)

引用回帖:
2楼: Originally posted by gemucai at 2012-09-25 08:52:10
QE5.0改成这样的警告了?挺有意思的
是这样的,你所使用的实空间FFT网格的格点并没有刚好取到平移对称性所在的位置,所以它只能忽略你体系的平移对称性。
解决方法也好办:
1.如果你对对称群没要求,或者不用进行 ...

你好,我在输出文件里找到两个FFT dimensions,应该看哪一个啊?
     Dense  grid:   279535 G-vectors     FFT dimensions: (  36,  72, 250)
     Smooth grid:    70809 G-vectors     FFT dimensions: (  24,  45, 160)

另外你说的贴子我之前看过了,怎么改晶格类型,让警告变成下面形式呢?我已经按照那贴子的建议用了a,b,c,cosab的。
warning: symmetry operation #  3 not allowed.   fractional translation:       a  b  c  in crystal coordinates

还有改ecutrho又该怎么改呢,麻烦您拿我这例子讲下啊,谢谢了!
7楼2012-09-25 10:12:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

引用回帖:
7楼: Originally posted by cj4566 at 2012-09-25 10:12:08
你好,我在输出文件里找到两个FFT dimensions,应该看哪一个啊?
     Dense  grid:   279535 G-vectors     FFT dimensions: (  36,  72, 250)
     Smooth grid:    70809 G-vectors     FFT dimensions: (  2 ...

首先你可以放心,不是因为这个极性面不稳定才出现的警告,只是因为你做优化移动原子破坏了对称性而已。

既然是warning而不是error,那就是说后果不严重,不影响电子计算结果。

我正在重复你的计算,有结果了会告诉你的。
8楼2012-09-25 15:37:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)

结果出来了,并没有出现任何的警告。然后我用优化好的结构去做scf计算也没有出现警告。你用的QE是什么版本?

另外,nr1, nr2, nr3的值是Dense  grid,不是smooth grid。
比方说你的Dense  grid是36,  72, 250,你可以设置nr1=40, nr2=80, nr3=300,然后运行pw.x,再看输出文件就会发现Dense  grid变成了40, 80 ,300。

但我重复你的计算并没有任何这个警告啊。。。
9楼2012-09-25 20:04:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tangosnow

木虫 (小有名气)

加大ecutwfc,QE的QFA里有这个问题
10楼2012-09-26 03:20:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cj4566 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 9/450 2026-08-18 20:09 by 奶牛小黑
[基金申请] 朋友圈看到的 +4 wangzilk 2026-08-18 5/250 2026-08-18 19:25 by wangzilk
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +51 医学老男孩 2026-08-13 110/5500 2026-08-18 19:07 by icerianxia
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 时间戳变了,能看出什么问题? +8 基诺咪客 2026-08-17 11/550 2026-08-18 17:13 by 悠然熊猫
[基金申请] 重要消息,中午系统在维护 +9 yuleib84 2026-08-18 10/500 2026-08-18 16:52 by archvillain
[教师之家] 两姐妹,一个辛辛苦苦30岁工作当大学教师,另一个23岁就工作当开车教练 +3 瞬息宇宙 2026-08-11 7/350 2026-08-18 16:42 by LNP@mRNA
[论文投稿] 投稿咨询 +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 淀粉搬运工
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[基金申请] 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
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
信息提示
请填处理意见