版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

返回列表

當(dāng)前只顯示滿(mǎn)足指定條件的回帖，點(diǎn)擊這里查看本話(huà)題的所有回帖

hanyanli0475

銅蟲(chóng) (小有名氣)

應(yīng)助: 11 (小學(xué)生)
金幣: 173.1
散金: 32
紅花: 3
帖子: 296
在線(xiàn): 130.5小時(shí)
蟲(chóng)號(hào): 1160645
注冊(cè): 2010-12-02
性別: GG
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] 石墨烯納米帶差分電荷密度

各位朋友大家好，我現(xiàn)在想計(jì)算一下，石墨烯納米帶用氫修飾邊緣后的差分電荷密度，我的思路是：1.計(jì)算含氫修飾的石墨烯納米帶的自洽，然后再做其對(duì)應(yīng)的電荷密度pp1，2.計(jì)算沒(méi)有被氫修飾的石墨烯納米帶的自洽，然后再做其對(duì)應(yīng)的電荷密度pp2，3.用pp1減去pp2，就得出差分電荷密度了。
但是，在運(yùn)行最后一步相減時(shí)，出現(xiàn)如下錯(cuò)誤：
from chdens : error #       1    incompatible celldm

大家看看這是怎么回事啊？？？我猜測(cè)，由于氫修飾石墨烯納米帶有18個(gè)原子，而沒(méi)有修飾的納米帶有14個(gè)原子，所以報(bào)錯(cuò)了……，對(duì)嗎？大家?guī)蛶兔�，這兩天著急看結(jié)果。謝謝大家了。
我貼出輸入文件，大家看看。
各位高手幫幫忙啊，現(xiàn)在很著急想想看看這個(gè)結(jié)果，已經(jīng)好幾天了。謝謝大家了。
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.沒(méi)有修飾
&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',
/

回復(fù)此樓

» 猜你喜歡

材料080500調(diào)劑求收留已經(jīng)有7人回復(fù)
一志愿南昌大學(xué)，327分，材料與化工085600 已經(jīng)有6人回復(fù)
環(huán)境工程調(diào)劑已經(jīng)有9人回復(fù)
工科材料085601 279求調(diào)劑已經(jīng)有9人回復(fù)
287求調(diào)劑已經(jīng)有5人回復(fù)
22408 344分求調(diào)劑一志愿華電計(jì)算機(jī)技術(shù) 已經(jīng)有3人回復(fù)
289求調(diào)劑已經(jīng)有3人回復(fù)
317求調(diào)劑已經(jīng)有10人回復(fù)
廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計(jì)算機(jī)專(zhuān)業(yè)調(diào)劑）已經(jīng)有3人回復(fù)
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？已經(jīng)有8人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

石墨烯的電荷性為什么為負(fù)？已經(jīng)有10人回復(fù)
修改RWIGS 維格納半徑之后，差分電荷密度會(huì)不會(huì)變化已經(jīng)有4人回復(fù)
怎么利用PWSCF計(jì)算差分電荷密度已經(jīng)有25人回復(fù)
電荷差分密度已經(jīng)有9人回復(fù)
VESTA做的差分電荷密度圖分析已經(jīng)有20人回復(fù)
再次求問(wèn)，差分電荷密度圖的定性分析已經(jīng)有17人回復(fù)
雙層石墨烯的能帶、pdos、以及差分電荷密度已經(jīng)有5人回復(fù)
【求教】電子密度，差分密度，原子電荷這三個(gè)之間有什么關(guān)系？已經(jīng)有7人回復(fù)
vasp如何計(jì)算差分電荷密度已經(jīng)有7人回復(fù)
差分電荷密度的計(jì)算已經(jīng)有11人回復(fù)
【求助】差分電荷密度已經(jīng)有19人回復(fù)
【請(qǐng)教】石墨烯的密度是多少啊，已經(jīng)有8人回復(fù)

什么時(shí)候能算完！�。�！

1樓 2013-01-31 22:38:10

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hanyanli0475

銅蟲(chóng) (小有名氣)

應(yīng)助: 11 (小學(xué)生)
金幣: 173.1
散金: 32
紅花: 3
帖子: 296
在線(xiàn): 130.5小時(shí)
蟲(chóng)號(hào): 1160645
注冊(cè): 2010-12-02
性別: GG
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

您好，謝謝您的回答，這個(gè)計(jì)算電荷密度的思路，也是在小木蟲(chóng)學(xué)習(xí)的，文獻(xiàn)上也是這樣給的。我也會(huì)好好看一下您提供的文獻(xiàn)，學(xué)習(xí)一下。今天早上我把原子數(shù)改相同了，但是還是不行，您能在幫我看看嗎?謝謝哈

[ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ]

贊一下

回復(fù)此樓

什么時(shí)候能算完！�。�！

3樓2013-02-01 11:07:31

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 8 個(gè)回答

lfhuang

木蟲(chóng) (著名寫(xiě)手)

應(yīng)助: 39 (小學(xué)生)
金幣: 3571.1
散金: 2141
紅花: 37
帖子: 1432
在線(xiàn): 681.4小時(shí)
蟲(chóng)號(hào): 505824
注冊(cè): 2008-02-17
性別: GG
專(zhuān)業(yè): 金屬材料表面科學(xué)與工程

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

兩個(gè)文件里的數(shù)據(jù)相減，你先要把兩個(gè)文件中電荷密度之外的信息（比如原子數(shù)，原子坐標(biāo)...應(yīng)該在文件頭）弄成完全一致了。
不過(guò)我很懷疑樓主算這樣的差分電荷的意義。對(duì)于graphene nanoribbon體系的差分電荷處理，類(lèi)似信息，可以參考下面參考中的公式13和圖6：
http://iopscience.iop.org/0953-8984/25/5/055304

贊一下

回復(fù)此樓

先進(jìn)使役材料計(jì)算（https://aidme.nimte.ac.cn）

2樓2013-02-01 01:53:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

souledge

專(zhuān)家顧問(wèn) (著名寫(xiě)手)

專(zhuān)家經(jīng)驗(yàn): +174
1ST強(qiáng)帖: 19
應(yīng)助: 303 (大學(xué)生)
金幣: 11306.4
散金: 1123
紅花: 108
帖子: 2814
在線(xiàn): 916.1小時(shí)
蟲(chóng)號(hào): 941486
注冊(cè): 2010-01-12
專(zhuān)業(yè): 結(jié)構(gòu)陶瓷
管轄: 第一性原理

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
hanyanli0475: 金幣+20, ★★★★★最佳答案 2013-02-01 14:13:53

引用回帖:

from chdens : error # 1 incompatible celldm

出錯(cuò)信息寫(xiě)得再明白不過(guò)了，incompatible celldm，不兼容的celldm。什么是celldm？那不就是&system中的celldm(i)參數(shù)么？說(shuō)白了，就是兩個(gè)結(jié)構(gòu)的尺寸不同。尺寸不同的東西你讓電腦怎么減？
我們知道，差分電荷密度是 $d \rho (x,y,z)=\rho (x,y,z) - \sum [{\rho}^{i} (x, y, z)]$ ，只有每個(gè)xyz能對(duì)應(yīng)完全相同才有意義。
另外，補(bǔ)充一下二樓的疑問(wèn)。二樓的疑問(wèn)是這么做差分有沒(méi)有什么物理意義，而不是這么做究竟對(duì)不對(duì)，是不是正統(tǒng)做差分的方法。事實(shí)上，差分這個(gè)東西也可以根據(jù)自己的需要來(lái)自己定義差分的方法。

贊一下

回復(fù)此樓

思想重于技巧，內(nèi)涵重于表象

4樓2013-02-01 12:07:58

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hanyanli0475

銅蟲(chóng) (小有名氣)

應(yīng)助: 11 (小學(xué)生)
金幣: 173.1
散金: 32
紅花: 3
帖子: 296
在線(xiàn): 130.5小時(shí)
蟲(chóng)號(hào): 1160645
注冊(cè): 2010-12-02
性別: GG
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

4樓: Originally posted by souledge at 2013-02-01 12:07:58
出錯(cuò)信息寫(xiě)得再明白不過(guò)了，incompatible celldm，不兼容的celldm。什么是celldm？那不就是&system中的celldm(i)參數(shù)么？說(shuō)白了，就是兩個(gè)結(jié)構(gòu)的尺寸不同。尺寸不同的東西你讓電腦怎么減？
我們知道，差分電 ...

非常感謝您的回答，我看了一下我的晶格常數(shù)，設(shè)置的都是一樣的呀，而且在計(jì)算摻雜時(shí)，就可以相減，不知道是為何原因？我應(yīng)該怎樣修改呢？

贊一下

回復(fù)此樓

什么時(shí)候能算完�。。。�

5樓2013-02-01 14:13:47

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 8 個(gè)回答

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 317求調(diào)劑 +5	申子申申 2026-03-19	10/500	2026-03-20 15:58 by 蔡誠(chéng)
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +9	手機(jī)用戶(hù) 2026-03-17	9/450	2026-03-20 14:24 by 無(wú)懈可擊111
[考研] 279分求調(diào)劑一志愿211 +9	chaojifeixia 2026-03-19	10/500	2026-03-20 12:29 by lature00
[考研] 274求調(diào)劑 +8	S.H1 2026-03-18	8/400	2026-03-20 11:53 by 學(xué)員8dgXkO
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11	Ymlll 2026-03-18	14/700	2026-03-20 10:46 by EBSD
[考研] 一志愿武漢理工材料工程專(zhuān)碩調(diào)劑 +7	Doleres 2026-03-19	7/350	2026-03-20 10:39 by EBSD
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5	脫穎而出 2026-03-16	14/700	2026-03-20 09:30 by kkcoco25
[考研] 266求調(diào)劑 +5	陽(yáng)陽(yáng)哇塞 2026-03-14	10/500	2026-03-19 15:08 by 陽(yáng)陽(yáng)哇塞
[考研] 一志愿西安交通大學(xué)材料工程專(zhuān)業(yè) 282分求調(diào)劑 +5	楓橋ZL 2026-03-18	7/350	2026-03-19 14:52 by 功夫瘋狂
[考研] 化學(xué)求調(diào)劑 +3	臨澤境llllll 2026-03-17	4/200	2026-03-19 13:59 by houyaoxu
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 0703化學(xué)336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考博] 26博士申請(qǐng) +3	1042136743 2026-03-17	3/150	2026-03-17 23:30 by 輕松不少隨
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 材料工程專(zhuān)碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 藥學(xué)383 求調(diào)劑 +3	藥學(xué)chy 2026-03-15	4/200	2026-03-16 20:51 by 元子^0^
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 求老師收留調(diào)劑 +4	jiang姜66 2026-03-14	5/250	2026-03-15 20:11 by Winj1e