版塊導航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

返回列表

hanyanli0475

銅蟲 (小有名氣)

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

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

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

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

回復此樓

» 猜你喜歡

求調(diào)劑已經(jīng)有3人回復
265求調(diào)劑已經(jīng)有4人回復
085700資源與環(huán)境308求調(diào)劑已經(jīng)有6人回復
一志愿吉林大學材料學碩321求調(diào)劑已經(jīng)有12人回復
286分人工智能專業(yè)請求調(diào)劑愿意跨考！已經(jīng)有3人回復
329求調(diào)劑已經(jīng)有5人回復
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復已經(jīng)有4人回復
材料學碩318求調(diào)劑已經(jīng)有5人回復
一志愿中國海洋大學，生物學，301分，求調(diào)劑已經(jīng)有6人回復
081700化工學碩調(diào)劑已經(jīng)有3人回復

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

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

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

什么時候能算完！�。�！

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

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

lfhuang

木蟲 (著名寫手)

應助: 39 (小學生)
金幣: 3571.1
散金: 2141
紅花: 37
帖子: 1432
在線: 681.4小時
蟲號: 505824
注冊: 2008-02-17
性別: GG
專業(yè): 金屬材料表面科學與工程

【答案】應助回帖

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

兩個文件里的數(shù)據(jù)相減，你先要把兩個文件中電荷密度之外的信息（比如原子數(shù)，原子坐標...應該在文件頭）弄成完全一致了。
不過我很懷疑樓主算這樣的差分電荷的意義。對于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

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

hanyanli0475

銅蟲 (小有名氣)

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

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

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

贊一下

回復此樓

什么時候能算完�。。。�

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

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

souledge

專家顧問 (著名寫手)

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

【答案】應助回帖

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

引用回帖:

from chdens : error # 1 incompatible celldm

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

贊一下

回復此樓

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

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

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

hanyanli0475

銅蟲 (小有名氣)

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

引用回帖:

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

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

贊一下

回復此樓

什么時候能算完�。。�！

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

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

hanyanli0475

銅蟲 (小有名氣)

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

引用回帖:

謝謝回答，問題解決，celldm（2）中，就差了一點點，呵呵。謝謝您。

贊一下

回復此樓

什么時候能算完！�。�！

6樓2013-02-01 15:57:03

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

chenkunkun

銅蟲 (小有名氣)

應助: 14 (小學生)
金幣: 269.6
帖子: 145
在線: 286.7小時
蟲號: 2370184
注冊: 2013-03-23
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

請問樓主你這加氫的石墨帶在scf算之前有沒有進行relax呢

贊一下

回復此樓

努力

7樓2013-04-18 15:19:59

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

18810950811

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 17
帖子: 3
在線: 20.4小時
蟲號: 3725589
注冊: 2015-03-10
專業(yè): 金屬基復合材料

引用回帖:

計算車紛電荷密度的時候要設(shè)置nx，ny 。這個設(shè)置的根據(jù)是什么呢希望能給將截一下謝謝啊

贊一下

回復此樓

8樓2015-10-24 10:08:46

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 hanyanli0475 的主題更新

返回列表

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

[考研] 085700資源與環(huán)境308求調(diào)劑 +6	墨墨漠 2026-03-18	6/300	2026-03-20 00:02 by 23Postgrad
[考研] 266求調(diào)劑 +5	陽陽哇塞 2026-03-14	10/500	2026-03-19 15:08 by 陽陽哇塞
[考研] 本人考085602 化學工程專碩 +17	不知道叫什么！ 2026-03-15	19/950	2026-03-19 15:06 by 盡舜堯1
[考研] 281求調(diào)劑（0805） +9	煙汐憶海 2026-03-16	19/950	2026-03-19 11:42 by laoshidan
[考研] 0703化學 305求調(diào)劑 +4	FY_yy 2026-03-14	4/200	2026-03-19 05:54 by anny19840123
[考研] 一志愿武理材料305分求調(diào)劑 +5	想上岸的鯉魚 2026-03-18	6/300	2026-03-18 17:53 by 無際的草原
[考研] 311求調(diào)劑 +11	冬十三 2026-03-15	12/600	2026-03-18 14:36 by 星空星月
[考研] 0703化學336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考研] 301求調(diào)劑 +9	yy要上岸呀 2026-03-17	9/450	2026-03-18 08:58 by 無際的草原
[考研] 278求調(diào)劑 +5	煙火先于春 2026-03-17	5/250	2026-03-18 08:43 by 星空星月
[考研] 334求調(diào)劑 +3	志存高遠意在機?/a> 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[考研] 材料專碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學員h26Tkc
[考研] 283求調(diào)劑 +3	聽風就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] [導師推薦]西南科技大學國防/材料導師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[基金申請] 今年的國基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤色
[考研] 0703 物理化學調(diào)劑 +3	我可以上岸的對?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對?/a>
[考研] 085601材料工程315分求調(diào)劑 +3	yang_0104 2026-03-15	3/150	2026-03-15 10:58 by peike
[考研] 080500，材料學碩302分求調(diào)劑學校 +4	初識可樂 2026-03-14	5/250	2026-03-14 21:08 by peike