24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

xiaoqiu007

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 356.5
散金: 7
帖子: 100
在线: 20.8小时
虫号: 707167
注册: 2009-02-23
专业: 金属功能材料

[求助] fcc格子，执行q2r.x时出错，提示，q not allowed

这是q2r.out中的信息
reading grid info from file c225.dyn0
  reading force constants from file c225.dyn1
  Dielectric Tensor not found
  nqs=          1
q= 0.00000000  0.00000000  0.00000000
  reading force constants from file c225.dyn2
  nqs=          3
q= 0.50000000 -0.50000000  0.50000000
q= 0.50000000  0.50000000 -0.50000000
q=  -0.50000000 -0.50000000 -0.50000000
  reading force constants from file c225.dyn3
  nqs=          6
q= 0.25000000  0.25000000  0.25000000
q=  -0.25000000 -0.25000000  0.25000000
q=  -0.25000000 -0.25000000 -0.25000000
q=  -0.25000000  0.25000000 -0.25000000

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
   from init : error #       1
   q not allowed
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

   stopping ...

回复此楼

» 猜你喜欢

谈谈两天一夜的“延安行” 已经有3人回复
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？已经有11人回复
氨基封端PDMS和HDI反应快速固化已经有10人回复
博士读完未来一定会好吗已经有37人回复
论文投稿求助已经有4人回复
Applied Surface Science 这个期刊。有哪位虫友投过的能把word模板发给我参考一下嘛已经有3人回复
投稿精细化工已经有6人回复
博士申请都是内定的吗？已经有9人回复
心脉受损已经有8人回复

» 本主题相关价值贴推荐，对您同样有帮助:

PH计使用，校正时数值显示出错，什么问题已经有6人回复
为什么要引入倒格子已经有5人回复
CP2K 跑 MD 为什么分子跑到了周期性格子外面已经有5人回复
四方格子的电导公式已经有4人回复
FCC2010,上海，至今未被检索已经有7人回复
fcc晶体的衍射花样已经有4人回复
炼厂各装置（FCC、焦化）生产的液化气中会不会含有氧气已经有3人回复
【求助】碳纳米管的倒格子空间已经有8人回复
【交流】一般情况下为什么fcc结构比bcc结构的塑性更好已经有24人回复
【讨论】哪位知道为什么FCC，BCC的枝晶的生长方向都是100？？已经有9人回复
【原创】endnote提示文件不相联即this URL could not be lauched的解决办法已经有16人回复
【求助】高斯03出错已经有9人回复

1楼2011-05-01 21:02:35

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

léon_WANG

木虫 (小有名气)

应助: 5 (幼儿园)
金币: 2577.5
散金: 20
红花: 4
帖子: 156
在线: 226.1小时
虫号: 2458360
注册: 2013-05-10
专业: 凝聚态物性I:结构、力学和

DO nq = 1,nqs WRITE(6,'(a,3f12.8)') ' q= ',(q(i,nq),i=1,3) lq = .TRUE. DO ipol=1,3 xq = 0.0d0 DO icar=1,3 xq = xq + at(icar,ipol) * q(icar,nq) * nr(ipol) END DO lq = lq .AND. (ABS(NINT(xq) - xq) .LT. eps) iq = NINT(xq) ! m(ipol)= MOD(iq,nr(ipol)) + 1 IF (m(ipol) .LT. 1) m(ipol) = m(ipol) + nr(ipol) END DO IF (.NOT.lq) CALL errore('init','q not allowed',1)
这怎么解？