 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2212 | 回復: 7 | ||||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||||
hanyanli0475銅蟲 (小有名氣)
|
[求助]
石墨烯納米帶差分電荷密度
|
|||
|
各位朋友大家好，我現(xiàn)在想計算一下，石墨烯納米帶用氫修飾邊緣后的差分電荷密度，我的思路是：1.計算含氫修飾的石墨烯納米帶的自洽，然后再做其對應的電荷密度pp1，2.計算沒有被氫修飾的石墨烯納米帶的自洽，然后再做其對應的電荷密度pp2，3.用pp1減去pp2，就得出差分電荷密度了。 但是，在運行最后一步相減時，出現(xiàn)如下錯誤： from chdens : error # 1 incompatible celldm 大家看看這是怎么回事��？？？我猜測，由于氫修飾石墨烯納米帶有18個原子，而沒有修飾的納米帶有14個原子，所以報錯了……，對嗎？大家?guī)蛶兔�，這兩天著急看結果。謝謝大家了。 我貼出輸入文件，大家看看。 各位高手幫幫忙啊，現(xiàn)在很著急想想看看這個結果，已經(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', / |
回復此樓» 猜你喜歡
求調劑
已經(jīng)有3人回復
-
265求調劑
已經(jīng)有4人回復
-
085700資源與環(huán)境308求調劑
已經(jīng)有6人回復
-
一志愿吉林大學材料學碩321求調劑
已經(jīng)有12人回復
-
286分人工智能專業(yè)請求調劑愿意跨考！
已經(jīng)有3人回復
-
329求調劑
已經(jīng)有5人回復
-
申請回稿延期一個月，編輯同意了。但系統(tǒng)上的時間沒變，給編輯又寫郵件了，沒回復
已經(jīng)有4人回復
-
材料學碩318求調劑
已經(jīng)有5人回復
-
一志愿中國海洋大學，生物學，301分，求調劑
已經(jīng)有6人回復
-
081700化工學碩調劑
已經(jīng)有3人回復
» 本主題相關商家推薦: (我也要在這里推廣)
» 本主題相關價值貼推薦，對您同樣有幫助:
- 石墨烯的電荷性為什么為負？
已經(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人回復
- 【求教】電子密度，差分密度，原子電荷這三個之間有什么關系？
已經(jīng)有7人回復
- vasp如何計算差分電荷密度
已經(jīng)有7人回復
- 差分電荷密度的計算
已經(jīng)有11人回復
- 【求助】差分電荷密度
已經(jīng)有19人回復
- 【請教】石墨烯的密度是多少啊，
已經(jīng)有8人回復
chenkunkun
銅蟲 (小有名氣)
- 應助: 14 (小學生)
- 金幣: 269.6
- 帖子: 145
- 在線: 286.7小時
- 蟲號: 2370184
- 注冊: 2013-03-23
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
7樓2013-04-18 15:19:59
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 |
2樓2013-02-01 01:53:42
hanyanli0475
銅蟲 (小有名氣)
- 應助: 11 (小學生)
- 金幣: 173.1
- 散金: 32
- 紅花: 3
- 帖子: 296
- 在線: 130.5小時
- 蟲號: 1160645
- 注冊: 2010-12-02
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
3樓2013-02-01 11:07:31
souledge
專家顧問 (著名寫手)
-
專家經(jīng)驗: +174 - 1ST強帖: 19
- 應助: 303 (大學生)
- 金幣: 11306.4
- 散金: 1123
- 紅花: 108
- 帖子: 2814
- 在線: 916.1小時
- 蟲號: 941486
- 注冊: 2010-01-12
- 專業(yè): 結構陶瓷
- 管轄: 第一性原理
【答案】應助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
hanyanli0475: 金幣+20, ★★★★★最佳答案 2013-02-01 14:13:53
感謝參與，應助指數(shù) +1
hanyanli0475: 金幣+20, ★★★★★最佳答案 2013-02-01 14:13:53
|
出錯信息寫得再明白不過了，incompatible celldm，不兼容的celldm。什么是celldm？那不就是&system中的celldm(i)參數(shù)么？說白了，就是兩個結構的尺寸不同。尺寸不同的東西你讓電腦怎么減？ 我們知道，差分電荷密度是 另外，補充一下二樓的疑問。二樓的疑問是這么做差分有沒有什么物理意義，而不是這么做究竟對不對，是不是正統(tǒng)做差分的方法。事實上，差分這個東西也可以根據(jù)自己的需要來自己定義差分的方法。 |
4樓2013-02-01 12:07:58
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085700資源與環(huán)境308求調劑 +6 | 墨墨漠 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 生物學調劑招人�。�！ +3 | 山海天嵐 2026-03-17 | 4/200 |
|
|
|
[考研] 復試調劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
|
[考研] 266求調劑 +5 | 陽陽哇塞 2026-03-14 | 10/500 |
|
|
|
[考研] 一志愿 西北大學 ，070300化學學碩，總分287，雙非一本，求調劑。 +3 | 晨昏線與星海 2026-03-19 | 3/150 |
|
|
|
[考研] 346求調劑[0856] +3 | WayneLim327 2026-03-16 | 6/300 |
|
|
|
[考研] 0817調劑 +3 | 沒有答案_ 2026-03-14 | 3/150 |
|
|
|
[考研] 本科鄭州大學物理學院，一志愿華科070200學碩，346求調劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
|
[考研] 材料專業(yè)求調劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|
|
[考研] 297求調劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
|
[考研] 304求調劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
|
[考博] 環(huán)境領域全國重點實驗室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
|
[考研] 考研求調劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[考研] 302求調劑 +4 | 小賈同學123 2026-03-15 | 8/400 |
|
|
|
[考研] 283求調劑 +3 | 聽風就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研] 326求調劑 +4 | 諾貝爾化學獎覬?/a> 2026-03-15 | 7/350 |
|
|
|
[考研] 0703 物理化學調劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
|
[考研] 0856求調劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|
|
|
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流 +4 | 自由的_飛翔 2026-03-13 | 5/250 |
|
|
|
[考研] 中科大材料專碩319求調劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|