版块导航: 正在加载中...

客户端APP下载

登录注册

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

24小时热门版块排行榜

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

xijinyang

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 358.4
帖子: 261
在线: 76.6小时
虫号: 886038

[交流] 【求助】运行q2r时候报错！

我想计算graphene的电声相互作用量，在计算到第四步q2r.x（由倒空间的动力学矩阵元变为实空间的）报错，具体错误提示是：
   from init : error #       1
   nc already filled: wrong q grid or wrong nr

我的第一步较密k-mesh的scf的输入文件为：

&control
calculation = 'scf'
restart_mode='from_scratch'
wf_collect=.true.
outdir = './scf'
prefix = 'graphene'
pseudo_dir = './'
/
&system
ibrav = 4
celldm(1) = 4.65088, celldm(3)=3.0119
nat = 2
ntyp = 1
nbnd = 16
ecutwfc = 60
occupations = 'smearing'
degauss = 0.001
smearing = 'gaussian'
la2F = .true.
/
&electrons
conv_thr = 1.D-8
diagonalization = 'david'
/
ATOMIC_SPECIES
C 12.0107 C.pz-vbc.UPF
ATOMIC_POSITIONS {crystal}
C  6.6666765924E-01  3.3333531848E-01  2.6727647101E-51
C  3.3333234076E-01  6.6666468152E-01 -2.6727647101E-51
K_POINTS {automatic}
72 72 1 0 0 0

第二步较疏k-mesh的scf输入文件为：
&control
calculation = 'scf'
restart_mode='from_scratch'
wf_collect=.true.
outdir = './scf'
prefix = 'graphene'
pseudo_dir = './'
/
&system
ibrav = 4
celldm(1) = 4.65088, celldm(3)=3.0119
nat = 2
ntyp = 1
nbnd = 16
ecutwfc = 60
occupations = 'smearing'
degauss = 0.001
smearing = 'gaussian'
/
&electrons
conv_thr = 1.D-8
diagonalization = 'david'
/
ATOMIC_SPECIES
C 12.0107 C.pz-vbc.UPF
ATOMIC_POSITIONS {crystal}
C  6.6666765924E-01  3.3333531848E-01  2.6727647101E-51
C  3.3333234076E-01  6.6666468152E-01 -2.6727647101E-51
K_POINTS {automatic}
36 36 1 0 0 0

第三步较疏q-mesh的ph.x的输入文件为：

Electron-phonon coefficients for graphene
&inputph
recover = .true.
tr2_ph=1.0d-14
prefix='graphene'
fildvscf='graphene_dv'
amass(1)=12.0107
outdir='./scf'
fildyn='graphene.dyn'
elph=.true.
trans=.true.
ldisp=.true.
nq1=24, nq2=24, nq3=1
/

第四步q2r.x的输入文件为：
&input
zasr='simple',
fildyn='graphene.dyn',
flfrc='graphene.fc',
la2F=.true.
/

难道是我q-mesh取得有问题？请大家指教！谢谢！

» 本帖已获得的红花（最新10朵）

qingeraaa

» 猜你喜欢

拟解决的关键科学问题还要不要写已经有10人回复
救命帖已经有5人回复
限项规定已经有5人回复
为什么nbs上溴没有产物点出现呢已经有9人回复
招博士已经有3人回复
存款400万可以在学校里躺平吗已经有35人回复
最失望的一年已经有18人回复
求推荐博导已经有4人回复
求推荐英文EI期刊已经有5人回复
疑惑？已经有5人回复

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

MS5.0安装全部成功了，运行时总是报错，根本用不了已经有13人回复
vasp跨节点运行出错，mpiexec_node-1 (handle_stdin_input 1089) 已经有5人回复
程序没报错中断，但停止运行已经有3人回复
linux下运行python报错，不知原因已经有12人回复
求助matlab-simulink运行时报错问题已经有3人回复
请教，为什么siesta编译中没报错，然而运行的时候跑不起来呢已经有9人回复
【求助】LAMMPS-运行报错已经有9人回复
【求助】vasp的运行时候报错：Error reading item 'IMAGES' from file INCAR. 已经有13人回复
【求助】MS运行报错 The job has failed. Download any results generated so far? 已经有20人回复
【求助】BSSE of O2运行出错已经有17人回复
【求助】vasp运行时mpi报错已经有12人回复
【求助】Orca 2.8的一个计算赝势的例子报错已经有12人回复
【求助】Gambit 中的两个报错问题已经有3人回复
【求助】请教运行报错分析已经有3人回复
【讨论】lammps运行eam/fs出错已经有38人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-01-20 09:11:22

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

souledge

专家顾问 (著名写手)

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

LZ的错误可能是和la2F有关的，虽然在说明中没有写，而且好象PWgui中也没有，不过若是要做el-ph的话，在SCF时，system namelist下也有个类似la2F的参数的，这里没有就会在后面的q2r.x处挂点~
具体的参数可以参考官方给出的la2F的算例。

3楼2013-11-30 20:00:14

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

查看全部 6 个回答

qingeraaa

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 10.7
帖子: 219
在线: 195.9小时
虫号: 1503611

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

送红花一朵

刚开始用Qe，有看到说q点数目需要整除K点网格楼主的问题如何解决的？

2楼2013-11-21 19:42:28

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

gemucai

木虫 (正式写手)

应助: 16 (小学生)
金币: 1774.6
帖子: 693
在线: 201.7小时
虫号: 816478

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

也许我的经历可以帮帮你。将所有的a2Fq2r.*文件全部删掉，然后再算。这个文件有记忆性，它只保存你第一次计算的网格值，后面的计算跟第一次不一样就会报错。这个问题我跟paolo反映了，他不理我，我也只有在小木虫上提醒大家了。

赞一下(1人)

4楼2013-12-19 10:56:55

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

xiangzhe123

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 31
帖子: 2
在线: 4.1小时
虫号: 3217806

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

楼主，我最近在运行q2r时候，出现了你这相同的错误，你这问题后来是怎么解决的啊？能指教一下吗？

5楼2015-01-27 19:12:00

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

查看全部 6 个回答