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

» 猜你喜欢

职称评审没过，求安慰已经有54人回复
垃圾破二本职称评审标准已经有19人回复
毕业后当辅导员了，天天各种学生超烦已经有5人回复
26申博自荐已经有3人回复
A期刊撤稿已经有4人回复

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

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

湘粤cxy

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 118.3
帖子: 36
在线: 13.5小时
虫号: 4800129

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

引用回帖:

4楼: Originally posted by gemucai at 2013-12-18 23:56:55
也许我的经历可以帮帮你。将所有的a2Fq2r.*文件全部删掉，然后再算。这个文件有记忆性，它只保存你第一次计算的网格值，后面的计算跟第一次不一样就会报错。这个问题我跟paolo反映了，他不理我，我也只有在小木虫上 ...

把a2Fq2r.*全删掉也还是会出问题呀