24小时热门版块排行榜    

查看: 2323  |  回复: 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

木虫 (正式写手)

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

我给一个式子就明白了



里面的Ecut是ecutrho,即能量截断,Ni是i方向上的点的数目,即输入设置里面的nri,ai是i方向上的晶格常数,它们之间有这样的关系。所以你改变Ecut的目的还是为了改变Ni,其实是一样的。
14楼2012-09-26 09:31:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 44 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本人考085602 化学工程 专硕 +7 不知道叫什么! 2026-03-15 8/400 2026-03-15 20:11 by 棒棒球手
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 东南大学364求调剂 +3 JasonYuiui 2026-03-15 3/150 2026-03-15 18:57 by 无际的草原
[考研] 化学工程321分求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-15 18:49 by a不易
[考研] 290求调剂 +3 孔志浩 2026-03-12 8/400 2026-03-15 15:30 by 孔志浩
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 化学工程321分求调剂(南京工业,浙江工业) +3 大米饭! 2026-03-09 4/200 2026-03-14 02:34 by JourneyLucky
[考研] 301求调剂 +3 归零lbm 2026-03-09 3/150 2026-03-14 02:20 by JourneyLucky
[考研] 332分材料工程调剂 +3 莓好时光海苔 2026-03-09 3/150 2026-03-14 02:03 by JourneyLucky
[考研] 求调剂 +6 yfihxh 2026-03-09 6/300 2026-03-14 01:18 by JourneyLucky
[考研] 26考研调剂 +3 ying123. 2026-03-10 3/150 2026-03-14 00:18 by JourneyLucky
[考研] 285 求调剂 资源与环境 一志愿北京化工大学 +3 未名考生 2026-03-10 3/150 2026-03-13 23:04 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 一志愿西南交大,材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 7/350 2026-03-13 16:06 by ruiyingmiao
[考研] 274求调剂 +3 S.H1 2026-03-12 3/150 2026-03-13 15:15 by JourneyLucky
[考研] 085600 材料与化工 295 求调剂 +10 dream…… 2026-03-10 12/600 2026-03-12 13:46 by dream……
[基金申请] 提交后的基金本子,已让学校撤回了,可否换口子提交 +3 dut_pfx 2026-03-10 3/150 2026-03-11 08:38 by kudofaye
[考研] 085602化工求调剂 +7 董boxing 2026-03-10 7/350 2026-03-10 17:07 by BruceLiu320
信息提示
请填处理意见