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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 石墨烯纳米带差分电荷密度
51/1返回列表
查看: 2200  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

hanyanli0475

铜虫 (小有名气)

[求助] 石墨烯纳米带差分电荷密度

各位朋友大家好，我现在想计算一下，石墨烯纳米带用氢修饰边缘后的差分电荷密度，我的思路是：1.计算含氢修饰的石墨烯纳米带的自洽，然后再做其对应的电荷密度pp1，2.计算没有被氢修饰的石墨烯纳米带的自洽，然后再做其对应的电荷密度pp2，3.用pp1减去pp2，就得出差分电荷密度了。
但是，在运行最后一步相减时，出现如下错误：
from chdens : error #         1      incompatible celldm

大家看看这是怎么回事啊？？？我猜测，由于氢修饰石墨烯纳米带有18个原子，而没有修饰的纳米带有14个原子，所以报错了……，对吗？大家帮帮忙，这两天着急看结果。谢谢大家了。
我贴出输入文件，大家看看。
各位高手帮帮忙啊，现在很着急想想看看这个结果，已经好几天了。谢谢大家了。
1.氢修饰
&CONTROL
calculation = 'scf' ,
restart_mode = 'from_scratch' ,
outdir = '/tmp/' ,
pseudo_dir = '/pseudo/' ,
prefix = 'GNRs1' ,
etot_conv_thr = 1.0D-5 ,
forc_conv_thr = 1.0D-4 ,
tstress = .true. ,
tprnfor = .true. ,
/
&SYSTEM
ibrav = 8,
celldm(1) = 12.85444,celldm(2) =2.647059,celldm(3) =0.621529,
nat = 18,ntyp = 2,
ecutwfc = 125.0 ,ecutrho = 500 ,nosym = .false.,
occupations = 'smearing' ,degauss = 0.02 ,smearing = 'mp' ,
/
&ELECTRONS
conv_thr = 1.0d-8 ,
mixing_beta = 0.7 ,
mixing_mode = 'plain' ,
diagonalization = 'david',
/
ATOMIC_SPECIES
C    12.00000  C.pz-vbc.UPF
H    1.00000   H.pz-vbc.UPF   
ATOMIC_POSITIONS crystal
C        0.866030823   0.704519630   0.694313986
C        0.866030561   0.300457518   0.694314669
C        0.866027893   0.434500388   0.701171789
C        0.866027929   0.570477950   0.701171509
C        0.866028609   0.638094818   0.201893797
C        0.866028570   0.366882713   0.201894300
C        0.866027532   0.502489004   0.199704828
C        0.866028510   0.366882792   0.866296944
C        0.866027540   0.502488999   0.868484013
C        0.866028591   0.638094715   0.866297098
C        0.866030847   0.704519700   0.373879718
C        0.866030637   0.300457490   0.373879469
C        0.866027908   0.434500383   0.367017742
C        0.866027925   0.570477984   0.367017815
H        0.866034074   0.757566662   0.818374200
H        0.866034006   0.247411263   0.818376532
H        0.866034046   0.757566942   0.249821844
H        0.866034000   0.247411048   0.249819748
K_POINTS automatic
1 1 21     0 0 0
电荷密度（氢修饰）
&inputpp
    prefix  = 'GNRs1',
    outdir = '/tmp/',
    filplot = 'GNRscharge1',
    plot_num= 0
/
&plot
    nfile = 1,
    filepp(1) = 'GNRscharge1',
    weight(1) = 1.0,
    iflag = 3,
    output_format = 5,
    fileout = 'GNRs1.rho.dat',
    x0(1) = 0.0, x0(2) =0.0, x0(3)= 0.0,
    e1(1) =1.0, e1(2)= 0.0, e1(3) = 0.0,
    e2(1) =0.0, e2(2)= 1.0, e2(3) = 0.0,
    e3(1) =0.0, e3(2)= 0.0, e3(3) = 1.0,
    nx=60, ny=70,nz=100,
/
2.没有修饰
&CONTROL
calculation = 'scf' ,
restart_mode = 'from_scratch' ,
outdir = '/tmp/' ,
pseudo_dir = '/pseudo/' ,
prefix = 'GNRs' ,
etot_conv_thr = 1.0D-5 ,
forc_conv_thr = 1.0D-4 ,
tstress = .true. ,
tprnfor = .true. ,
/
&SYSTEM
ibrav = 8,
celldm(1) = 12.85444,celldm(2) = 2.64706,celldm(3) =0.621529,
nat = 14,ntyp = 1,
ecutwfc = 125.0 ,ecutrho = 500 ,nosym = .false.,
occupations = 'smearing' ,degauss = 0.02 ,smearing = 'mp' ,
/
&ELECTRONS
conv_thr = 1.0d-8 ,
mixing_beta = 0.7 ,
mixing_mode = 'plain' ,
diagonalization = 'david',
/
ATOMIC_SPECIES
C 12.00000  C.pz-vbc.UPF  
ATOMIC_POSITIONS crystal
C        0.866030145   0.700155656   0.679081001
C        0.866029909   0.304821192   0.679081287
C        0.866030105   0.434673629   0.700325500
C        0.866030142   0.570303990   0.700325660
C        0.866029908   0.638126091   0.197297690
C        0.866029779   0.366850839   0.197297344
C        0.866029841   0.502488587   0.199817032
C        0.866029844   0.366850921   0.870892872
C        0.866029877   0.502488591   0.868370667
C        0.866029971   0.638126029   0.870892498
C        0.866030295   0.700155402   0.389111299
C        0.866030010   0.304821478   0.389111239
C        0.866030061   0.434673655   0.367863033
C        0.866030114   0.570303941   0.367862876
K_POINTS automatic
1 1 21     0 0 0
电荷密度
&inputpp
    prefix  = 'GNRs2',
    outdir = '/tmp/',
    filplot = 'GNRscharge2',
    plot_num= 0
/
&plot
    nfile = 1,
    filepp(1) = 'GNRscharge2',
    weight(1) = 1.0,
    iflag = 3,
    output_format = 5,
    fileout = 'GNRs2.rho.dat',
    x0(1) = 0.0, x0(2) =0.0, x0(3)= 0.0,
    e1(1) =1.0, e1(2)= 0.0, e1(3) = 0.0,
    e2(1) =0.0, e2(2)= 1.0, e2(3) = 0.0,
    e3(1) =0.0, e3(2)= 0.0, e3(3) = 1.0,
    nx=60, ny=70,nz=100,
/

3.相减
&inputpp
/
&plot
    nfile = 2,
    filepp(1) = 'GNRscharge1',
    filepp(2) = 'GNRscharge2',
    weight(1) = 1.0,
    weight(2) = -1.0,
    iflag = 3,
    output_format = 5,
    fileout = 'GNRsdiff.rho.dat',
/
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
什么时候能算完！！！！
1楼 2013-01-31 22:38:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

18810950811

新虫 (初入文坛)
引用回帖:
4楼: Originally posted by souledge at 2013-02-01 12:07:58
出错信息写得再明白不过了，incompatible celldm，不兼容的celldm。什么是celldm？那不就是&system中的celldm(i)参数么？说白了，就是两个结构的尺寸不同。尺寸不同的东西你让电脑怎么减？
我们知道，差分电荷密 ...

计算车纷电荷密度的时候要设置nx，ny  。 这个设置的根据是什么呢  希望能给将截一下 谢谢啊
赞一下 回复此楼
8楼2015-10-24 10:08:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

lfhuang

木虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
两个文件里的数据相减，你先要把两个文件中电荷密度之外的信息（比如原子数，原子坐标...应该在文件头）弄成完全一致了。
不过我很怀疑楼主算这样的差分电荷的意义。对于graphene nanoribbon体系的差分电荷处理，类似信息，可以参考下面参考中的公式13和图6：
http://iopscience.iop.org/0953-8984/25/5/055304
赞一下 回复此楼
先进使役材料计算（https://aidme.nimte.ac.cn）
2楼2013-02-01 01:53:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanyanli0475

铜虫 (小有名气)
您好，谢谢您的回答，这个计算电荷密度的思路，也是在小木虫学习的，文献上也是这样给的。我也会好好看一下您提供的文献，学习一下。今天早上我把原子数改相同了，但是还是不行，您能在帮我看看吗?谢谢哈

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
什么时候能算完！！！！
3楼2013-02-01 11:07:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
hanyanli0475: 金币+20, ★★★★★最佳答案 2013-02-01 14:13:53
引用回帖:
from chdens : error #         1      incompatible celldm

出错信息写得再明白不过了，incompatible celldm，不兼容的celldm。什么是celldm？那不就是&system中的celldm(i)参数么？说白了，就是两个结构的尺寸不同。尺寸不同的东西你让电脑怎么减？
我们知道，差分电荷密度是，只有每个xyz能对应完全相同才有意义。
另外，补充一下二楼的疑问。二楼的疑问是这么做差分有没有什么物理意义，而不是这么做究竟对不对，是不是正统做差分的方法。事实上，差分这个东西也可以根据自己的需要来自己定义差分的方法。
赞一下 回复此楼
思想重于技巧，内涵重于表象
4楼2013-02-01 12:07:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 能源材料化学课题组招收硕士研究生8-10名 +3 脱颖而出 2026-03-16 6/300 2026-03-17 21:19 by z1z2z3879
[考研] 材料与化工求调剂 +6 为学666 2026-03-16 6/300 2026-03-17 20:15 by peike
[考研] 304求调剂 +8 小熊joy 2026-03-14 8/400 2026-03-17 17:29 by ruiyingmiao
[考研] 271求调剂 +12 生如夏花… 2026-03-11 14/700 2026-03-17 10:56 by lovewei0727
[考研] 一志愿，福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 283求调剂 +10 小楼。 2026-03-12 14/700 2026-03-16 16:08 by 13811244083
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 255求调剂 +3 李嘉慧， 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[考研] 招收0805（材料）调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 337一志愿华南理工0805材料求调剂 +7 mysdl 2026-03-11 9/450 2026-03-13 22:43 by JourneyLucky
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] 一志愿西南交大，材料专硕317求调剂 +5 lx8568 2026-03-11 5/250 2026-03-13 21:43 by peike
[考研] 工科，求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] （081700）化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
信息提示
关闭
请填处理意见
关闭