24小时热门版块排行榜    

查看: 2203  |  回复: 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的回帖

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的回帖

cj4566

木虫 (正式写手)

引用回帖:
9楼: Originally posted by gemucai at 2012-09-25 20:04:50
结果出来了,并没有出现任何的警告。然后我用优化好的结构去做scf计算也没有出现警告。你用的QE是什么版本?

另外,nr1, nr2, nr3的值是Dense  grid,不是smooth grid。
比方说你的Dense  grid是36,  72, 250, ...

这奇怪了,我对这个晶面建了好几个模型计算完都有这个警告。我用的是QE 5.0 并行版本啊,你呢?
     Program PWSCF v.5.0        starts on 24Sep2012 at 15:55: 0
     This program is part of the open-source Quantum ESPRESSO suite for quantum simulation of materials;
。。。
     Parallel version (MPI), running on     8 processors
     R & G space division:  proc/pool =    8
。。。
11楼2012-09-26 06:32:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cj4566

木虫 (正式写手)

引用回帖:
10楼: Originally posted by tangosnow at 2012-09-26 03:20:16
加大ecutwfc,QE的QFA里有这个问题

谢谢啊,我加大ecutwfc试试,不太会改nr1,2,3
12楼2012-09-26 06:36:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cj4566

木虫 (正式写手)

引用回帖:
15楼: Originally posted by gemucai at 2012-09-26 09:38:29
别忘了改ecutwfc的同时要改ecutrho,能量截断才是影响FFT网格的因素。...

好的,谢谢啊,这两个要一起改是吧。我刚才没改ecutwfc,把ecutrho增加到500还是出现警告,现在同时改ecutwfc和ecutrho试试。
16楼2012-09-26 10:48:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cj4566

木虫 (正式写手)

引用回帖:
17楼: Originally posted by gemucai at 2012-09-26 11:27:16
这警告对能量没有影响。需要看声子模式吗?如果不需要,这警告真就可以忽略的。
因为你是单斜体系,用的又是正交的FFT网格,你算的可能就是那个倾斜的面,两套格子本身就不匹配,问题恐怕很难弄好。
刚刚帮你找了 ...

好的,谢谢了,可能的确是晶格问题,我计算的是斜切的面。我不计算声子的,我就不管它了。
18楼2012-09-26 11:41:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cj4566

木虫 (正式写手)

引用回帖:
17楼: Originally posted by gemucai at 2012-09-26 11:27:16
这警告对能量没有影响。需要看声子模式吗?如果不需要,这警告真就可以忽略的。
因为你是单斜体系,用的又是正交的FFT网格,你算的可能就是那个倾斜的面,两套格子本身就不匹配,问题恐怕很难弄好。
刚刚帮你找了 ...

这个问题搞定了,把单斜晶系改成正交晶系就好了。
19楼2012-09-27 05:48:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cj4566

木虫 (正式写手)

引用回帖:
34楼: Originally posted by 刘仕晨 at 2013-03-29 13:26:13
不是的,当ecutwfc=150时,总能为-167.74382003 Ry,当ecutwfc=160时,总能量为-167.78218062 Ry,相差0.52ev,我觉得是不是还需要继续加大ecutwfc?正常应该收敛到相差10mev么,就是奇怪,居然要这么大的ecut...

是不是结构有问题
35楼2013-03-29 14:41:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cj4566 的主题更新
信息提示
请填处理意见