 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 6629 | 回復(fù): 14 | |||||||||
Gina88木蟲(chóng) (正式寫手)
|
[求助]
bader電荷布局分析怎么運(yùn)行��？
|
||||||||
|
bader電荷布局分析怎么運(yùn)行啊？ 先謝謝大俠們了！ |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 第一性原理相關(guān)文檔 | vasp計(jì)算資料 | 第一性 | 計(jì)算化學(xué) |
| 第一性原理 | 科研啊科研 | soft study |
» 猜你喜歡
265求調(diào)劑
已經(jīng)有8人回復(fù)
-
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有21人回復(fù)
-
295材料求調(diào)劑，一志愿武漢理工085601專碩
已經(jīng)有4人回復(fù)
-
0856調(diào)劑，是學(xué)校就去
已經(jīng)有9人回復(fù)
-
298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有9人回復(fù)
-
求調(diào)劑
已經(jīng)有4人回復(fù)
-
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？
已經(jīng)有7人回復(fù)
-
材料學(xué)碩297已過(guò)四六級(jí)求調(diào)劑推薦
已經(jīng)有6人回復(fù)
-
295復(fù)試調(diào)劑
已經(jīng)有7人回復(fù)
-
本人考085602 化學(xué)工程 專碩
已經(jīng)有20人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- Vasp能計(jì)算電荷轉(zhuǎn)移的具體數(shù)值么？
已經(jīng)有6人回復(fù)
- 求助vasp5.2能用的chsum.pl和bader小程序
已經(jīng)有16人回復(fù)
- Mulliken 布局分析任何分析電荷分布、轉(zhuǎn)移和化學(xué)性質(zhì)
已經(jīng)有10人回復(fù)
- 怎么使用bader電荷分析出原子軌道上的電子數(shù)，或者其它能夠得出軌道電子的軟件？
已經(jīng)有29人回復(fù)
- 缺陷(摻雜)化合價(jià)的確定
已經(jīng)有3人回復(fù)
- siesta計(jì)算dipole moment（電偶極矩）
已經(jīng)有7人回復(fù)
- 用DMol3計(jì)算離子型化合物電荷的布局分布
已經(jīng)有9人回復(fù)
- OUTCAR中的 total charge 怎么來(lái)的�。�
已經(jīng)有16人回復(fù)
- 【任務(wù)征集】VASP跑分子動(dòng)力學(xué)數(shù)據(jù)處理方法
已經(jīng)有20人回復(fù)
- 【分享】SIESTA中Bader analysis
已經(jīng)有43人回復(fù)
- 【其他】學(xué)習(xí)bader的交流
已經(jīng)有9人回復(fù)
- 【求助】vasp如何來(lái)查看電荷轉(zhuǎn)移情況
已經(jīng)有15人回復(fù)
- 【求助】castep計(jì)算電荷布局?jǐn)?shù)分析
已經(jīng)有9人回復(fù)
- 【求助】怎么看電荷布局?jǐn)?shù)
已經(jīng)有4人回復(fù)
- 【求助】bader電荷分析的幾個(gè)問(wèn)題
已經(jīng)有11人回復(fù)
- 【求助】如何解釋氫鍵的鍵電荷布局是負(fù)數(shù)
已經(jīng)有3人回復(fù)
Gina88
木蟲(chóng) (正式寫手)
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 9 (幼兒園)
- 金幣: 3328.8
- 散金: 124
- 紅花: 8
- 帖子: 395
- 在線: 113小時(shí)
- 蟲(chóng)號(hào): 1047684
- 注冊(cè): 2010-06-26
- 專業(yè): 理論和計(jì)算化學(xué)
|
Mulliken電荷是一種計(jì)算原子partial charge（局部電荷）的方法。本來(lái)，電子是在空間分布的一團(tuán)電子云，但是廣大玩計(jì)算的童鞋們希望能有一種更簡(jiǎn)潔的方法來(lái)描述這種混沌不清的分布，看看誰(shuí)多誰(shuí)少，誰(shuí)剝削了誰(shuí)，因此計(jì)算原子電荷就非常流行，稱為布局分析(population analysis)，通俗說(shuō)，如何在原子間分割大家共同的家產(chǎn)。大致說(shuō)先后提出了三種分割家產(chǎn)的思路。 第一類是以原子軌道對(duì)分子軌道的貢獻(xiàn)劃分。其中最著名的是Mulliken童鞋1955年提出的Mulliken population analuysis. Mullken童鞋說(shuō)，如果兩個(gè)原子共用了一些電子對(duì)，那么電子可以分成3部分，屬于A的(凈原子布居)， 屬于B的， 和共同的(重疊布居)。最公平的方法是，把大家共同擁有的一分為二，各拿一半，加上自己原有的，就是你的局部電荷。怎么計(jì)算凈原子布居和重疊布居？根據(jù)分子軌道中的原子軌道系數(shù)。 因?yàn)榉肿榆壍朗窃榆壍赖木�性組合，假設(shè)兩原子核外各有一條電子軌道A, B組成了一個(gè)電子對(duì)分子軌道C=xA+yB，總電荷等于波函數(shù)的平方：|C|^2=2個(gè)電子； 仔細(xì)看這個(gè)公式，誰(shuí)是誰(shuí)的不用說(shuō)了吧。。。類似依軌道劃分的還有所謂NBO電荷，Born電荷。 有人問(wèn)平面波基組并非原子中心基組，沒(méi)有原子軌道系數(shù)，怎么求Mulliken電荷呢？一般的平面波方法都是另外用一套原子中心基函數(shù)去擬合平面波計(jì)算得到的電子云分布密度，然后再套用Mulliken分析方法。 第二種分割思路是按波函數(shù)節(jié)點(diǎn)劃分，這就是nb哄哄的Bader charge， 依據(jù)所謂的AIM理論，atom in molecule。Bader童鞋非常憤青，經(jīng)常用很犀利的語(yǔ)言駁斥各種電荷劃分。他說(shuō)依軌道劃分類似于在兩原子間砍了一刀，各拿一半，這是沒(méi)有任何物理意義的。那么什么叫有物理意義？他說(shuō)你這一刀砍下去，要順著波函數(shù)的節(jié)點(diǎn)砍下去。我們知道波函數(shù)的節(jié)點(diǎn)是波函數(shù)變號(hào)的地方，對(duì)應(yīng)的電荷密度為0. 他說(shuō)這是電子云空間分割最理想最天然的分割點(diǎn)，比你那軌道劃分好5倍，呵呵。不可否認(rèn)，分子軌道的空間節(jié)點(diǎn)可以顯示成鍵的中心位置等信息，的確是一種比較好的劃分方法。 第三種思路，都是搞力場(chǎng)的童鞋感興趣的，叫ESP電荷，包括原裝ESP, 改進(jìn)版的ChelpG, Hershifiled等方法。搞力場(chǎng)的童鞋說(shuō)，你們搞切割的，雖然思路清晰，但是我拿到我的力場(chǎng)中一用，計(jì)算靜電作用，發(fā)現(xiàn)不行，連分子偶極距都是錯(cuò)的。怎么辦呢？俺們也不搞切割了，俺們拿一個(gè)探測(cè)電荷放到分子外頭探，探探受到了多大的靜電作用；把空間均勻劃分成幾千個(gè)格點(diǎn)，探幾千次；然后呢，我就在原子核的位置點(diǎn)嘗試放置一些電荷，來(lái)擬合探測(cè)電荷受到的靜電作用。不斷嘗試，最后最能夠重現(xiàn)分子靜電場(chǎng)的那一套原子電荷就是俺要的電荷。這就是原始的ESP電荷分析。后來(lái)發(fā)展的ChelpG什么，就是在擬合中加各種約束，比如增加偶極距的權(quán)重。這種電荷分析方法基本上是經(jīng)驗(yàn)擬合，但是能夠正確重現(xiàn)偶極距，分子間作用等性質(zhì)，搞力場(chǎng)的最喜歡。對(duì)他們來(lái)說(shuō)，只要結(jié)果正確即可，何須在意經(jīng)由何種途徑？唯一不爽的是屏蔽效應(yīng)很嚴(yán)重，即包裹在分子里面的原子的電荷不容易探到它們的靜電作用，擬合時(shí)你電荷給大點(diǎn)給小點(diǎn)，擬合結(jié)果都差不多，郁悶。。。 不全的大家補(bǔ)充。 |
3樓2013-02-04 16:31:36
Gina88
木蟲(chóng) (正式寫手)
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 9 (幼兒園)
- 金幣: 3328.8
- 散金: 124
- 紅花: 8
- 帖子: 395
- 在線: 113小時(shí)
- 蟲(chóng)號(hào): 1047684
- 注冊(cè): 2010-06-26
- 專業(yè): 理論和計(jì)算化學(xué)
|
幾個(gè)月前，當(dāng)我第一次打開(kāi)makefile的時(shí)候，我很暈，不知道該看什么，也不知道該干什么。如果能讀懂這個(gè)簡(jiǎn)單的makefile，那就基本上可以解決問(wèn)題了。 一個(gè)簡(jiǎn)單的程序，如果只有不到10行的主程序project.F90，那直接用編譯命令 ifort project.F90 -o project 就可以編譯成可執(zhí)行文件project了。 make文件是一個(gè)腳本。如果一個(gè)復(fù)雜的程序有多個(gè)模塊需要編譯，那可能需要檢查文件依賴關(guān)系，硬件和操作系統(tǒng)，設(shè)置環(huán)境變量，調(diào)用鏈接和數(shù)學(xué)庫(kù)之類的。一項(xiàng)一項(xiàng)的手動(dòng)編譯就很繁瑣了。所以就有了makefile。它是把所有編譯命令和參數(shù)都集中在一個(gè)命名為Makefile或者makefile的文件中，稱之為make文件。然后在存放makefile的目錄下運(yùn)行make命令就好了。make會(huì)自動(dòng)尋找當(dāng)前目錄下的Makefile或者makefile. 也可以make makefile，這時(shí)可以指定makefile文件名，比如make yindu_a_san bader程序是對(duì)vasp電荷做布局分析的. 我的平臺(tái)是linux,32位；編譯器是ifort. 所以進(jìn)入makefiles文件夾，選擇Makefile.lnx_ifort_32這個(gè)make文件。 Makefile.lnx_ifort_32內(nèi)容如下： .SUFFIXES: .f90 #源文件名后綴為f90。如果是fortran和C源代碼都有的呢？.SUFFIXES: .f90 .c ? FC = ifort #這里指出編譯器。如果不能在搜索路徑$PATH中找到，那就需要給出絕對(duì)路徑。如/opt/intel/fc/9.1.036/bin/ifort FFLAGS = -O2 #定義優(yōu)化級(jí)別，不知道是什么意思�？梢赃xO1,O2,O3. LINK = -static #定義編譯鏈接方式。靜態(tài)鏈接。編譯時(shí)把鏈接的庫(kù)文件編譯進(jìn)程序，而不是動(dòng)態(tài)調(diào)用。不清楚。 OBJS = kind_mod.o \ #定義變量。把 多個(gè)*.o文件定義為OBJS. matrix_mod.o \ #"\"表示繼續(xù)上一行？也就是按照一行輸入？ ions_mod.o \ options_mod.o \ charge_mod.o \ chgcar_mod.o \ cube_mod.o \ io_mod.o \ bader_mod.o \ voronoi_mod.o \ command_arg.o %.o %.mod : %.f90 #把*.f90文件編譯成什么？ 不明白。 $(FC) $(FFLAGS) -c $*.f90 #編譯命令。啥意思？ bader: $(OBJS) main.o #為什么不把main.o也寫進(jìn)OBJS呢？ rm -f bader #如果有老版本的bader可執(zhí)行文件，刪除之。為新版本的bader留位子。 $(FC) $(LINK) main.o -o $@ $(OBJS) #啥意思？ @？ 跟前面一條編譯命令有什么關(guān)系？ dist: bader #這個(gè)什么意思？ 我看刪掉這一行以及下面兩行都沒(méi)有關(guān)系。 tar -cf bader_lnx_32.tar bader gzip -9 bader_lnx_32.tar clean: #告訴make clean的時(shí)候刪除那些文件。有則刪之。 rm -f *.o *.mod bader bader_lnx_32.tar.gz #上面一行說(shuō)要?jiǎng)h除bader，呵呵，如果make后接著就make clean，那就得不到可執(zhí)行文件bader了。所以make后，及時(shí)把bader轉(zhuǎn)移到搜索路徑，如/user/bin下。或者直接修改為 rm -f *.o *.mod 很多我也不明白。不過(guò)已經(jīng)不怕它了。baidu google一下基本也可以解決一下小問(wèn)題。但是要自己從空白來(lái)寫makefile，我就無(wú)能為力了。 寫過(guò)Fortran的人，是否也會(huì)自己寫makefile呢？希望來(lái)分享解釋一下。 |
4樓2013-02-04 17:11:07
Gina88
木蟲(chóng) (正式寫手)
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 9 (幼兒園)
- 金幣: 3328.8
- 散金: 124
- 紅花: 8
- 帖子: 395
- 在線: 113小時(shí)
- 蟲(chóng)號(hào): 1047684
- 注冊(cè): 2010-06-26
- 專業(yè): 理論和計(jì)算化學(xué)
|
請(qǐng)高人科普一下bader 電荷分析一些細(xì)節(jié)。。 先發(fā)一個(gè)布局?jǐn)?shù)分析的帖子吧。需要到bader analysis 下載兩個(gè)腳本，chsum.pl, bader , 1。在靜態(tài)計(jì)算時(shí)，在INCAR 中加入LAECHARG 參數(shù)，并設(shè)置LAECHG =.TRUE. 運(yùn)算完畢會(huì)生成AECCAR0 AECCAR1 AECCER2, 然后運(yùn)行chsum.pl AECCAR0 AECCAR2, 會(huì)生成CHARGCAR_sum 文件， 2.然后運(yùn)行 bader CHGCAR -ref CHGCAR_sum 其中ACF.dat 就包含價(jià)電子的電荷信息。 有幾點(diǎn)要注意： A. Bader 只支持PAW 的方法， B。需要vasp4.6.31 以及以上版本。 C.計(jì)算電荷布局分析時(shí)，需要比較密的FFT 網(wǎng)格，可以通過(guò)NG(X,Y,Z) F來(lái)調(diào)整，增大這些量，直到總的電荷數(shù)目正確。 D。Bader 電荷布居分析基于電荷密度。 E和Mulliken ,Hirhfeld 分析比較，原子電荷：Hirshfeld |
博士PHD2樓2013-02-04 15:37:16
5樓2013-10-10 11:35:32
6樓2013-10-10 11:35:40
博士PHD
新蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 107.4
- 帖子: 111
- 在線: 82小時(shí)
- 蟲(chóng)號(hào): 1372524
- 注冊(cè): 2011-08-19
- 專業(yè): 理論和計(jì)算化學(xué)
7樓2015-07-08 21:10:12
博士PHD
新蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 107.4
- 帖子: 111
- 在線: 82小時(shí)
- 蟲(chóng)號(hào): 1372524
- 注冊(cè): 2011-08-19
- 專業(yè): 理論和計(jì)算化學(xué)
8樓2015-07-08 21:10:20
微6涼
鐵蟲(chóng) (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 202.8
- 帖子: 89
- 在線: 35.6小時(shí)
- 蟲(chóng)號(hào): 2724445
- 注冊(cè): 2013-10-14
- 專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
9樓2015-08-01 18:26:57
微6涼
鐵蟲(chóng) (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 202.8
- 帖子: 89
- 在線: 35.6小時(shí)
- 蟲(chóng)號(hào): 2724445
- 注冊(cè): 2013-10-14
- 專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
10樓2015-08-01 19:14:07
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 Gina88 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +21 | rare12345 2026-03-18 | 21/1050 |
|
|---|---|---|---|---|
|
[考研] 0856調(diào)劑，是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 8/400 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +6 | z1z2z3879 2026-03-18 | 6/300 |
|
|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分，有科研經(jīng)歷，求調(diào)劑 +3 | 一瓶苯甲酸 2026-03-14 | 3/150 |
|
|
|
[考研] 324分 085600材料化工求調(diào)劑 +3 | llllkkkhh 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-19 | 3/150 |
|
|
|
[考研] 346求調(diào)劑[0856] +3 | WayneLim327 2026-03-16 | 6/300 |
|
|
|
[考研] 330求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
|
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭(zhēng) 2026-03-15 | 3/150 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|