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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】运行出错
111/2返回列表12下一页
查看: 1367  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭