 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 6138 | 回復(fù): 22 | ||||
loovfnd至尊木蟲 (著名寫手)
|
[交流]
【求助】Molecular electrostatic potential (MEP) maps 分子靜電勢怎樣計(jì)算？ 已有5人參與
|
|||
|
J. AM. CHEM. SOC. 2006, 128, 1402-1403。 如下圖所示： [ Last edited by yjcmwgk on 2009-11-14 at 13:41 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 自己感覺好的東西 |
» 猜你喜歡
誠邀有志之士加入江蘇大學(xué)環(huán)境學(xué)院合成生物學(xué)課題組！(長期接收學(xué)碩專碩、碩博連讀)
已經(jīng)有16人回復(fù)
-
高效微納米儲(chǔ)能材料與產(chǎn)業(yè)化推進(jìn)實(shí)驗(yàn)室 蘇州 招收碩士研究生（含調(diào)劑）
已經(jīng)有2人回復(fù)
-
物理化學(xué)論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有130人回復(fù)
-
化學(xué)化工學(xué)院 招收化工、化學(xué)、材料等相關(guān)方向研究生（學(xué)碩、專碩都有調(diào)劑名額）
已經(jīng)有1人回復(fù)
-
一志愿南大理論與計(jì)算化學(xué)方向（070300）342分，求調(diào)劑
已經(jīng)有1人回復(fù)
-
歡迎化學(xué)、化工、物理、材料、或計(jì)算機(jī)背景的考研咨詢
已經(jīng)有0人回復(fù)
-
高效微納米儲(chǔ)能材料與產(chǎn)業(yè)化推進(jìn)實(shí)驗(yàn)室 蘇州 招收碩士研究生（含調(diào)劑）
已經(jīng)有0人回復(fù)
-
溫州大學(xué)新能源材料微觀溯源課題組招收2026級(jí)博士研究生（材料、化學(xué)方向）
已經(jīng)有6人回復(fù)
-
材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流
已經(jīng)有5人回復(fù)
-
26申博
已經(jīng)有2人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 【求助】怎么將siesta中的靜電勢ElectrostaticPotential的數(shù)據(jù)導(dǎo)出啊
已經(jīng)有4人回復(fù)
zhou2009
版主 (著名寫手)
- QC強(qiáng)帖: 15
- 應(yīng)助: 337 (大學(xué)生)
- 貴賓: 4.053
- 金幣: 42654.6
- 紅花: 153
- 帖子: 1729
- 在線: 1274.7小時(shí)
- 蟲號(hào): 676800
- 注冊: 2008-12-18
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
- 管轄: 量子化學(xué)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lei0736(金幣+8,VIP+0):謝謝 7-10 16:37
loovfnd(金幣+5,VIP+0):那計(jì)算命令呢？輸入行要怎么寫呢。 7-10 22:11
lei0736(金幣+8,VIP+0):謝謝 7-10 16:37
loovfnd(金幣+5,VIP+0):那計(jì)算命令呢？輸入行要怎么寫呢。 7-10 22:11
|
你可能是問如何作這樣的靜電勢圖吧？ 我沒有作過這樣的圖。如果要我作這樣的圖，我會(huì)這樣去作： 1、 Windows版cubegen的提示信息與回答： 進(jìn)行g(shù)03對(duì)圖示分子進(jìn)行計(jì)算，得到.chk文件，改寫成.fch文件 運(yùn)行cubegen， Property [Density]? Potential …… 這樣就得到圖示分子的分子靜電勢空間格點(diǎn)文件name.cube 2、 運(yùn)行g(shù)sgrid程序，選擇圖示分子有靜電勢等值線圖的那個(gè)截面：如果圖面是xy，即選其為截面，z則選生成該靜電勢等值線圖的那些原子的，它們都在一個(gè)平面，只有一個(gè)z值，且這個(gè)截面之外的原子都不會(huì)截取到，沒有值。形成gsgrid.txt文件，即時(shí)改為自己需要的名字name.txt 3、 進(jìn)入sigmaplot中作圖。 圖板的左邊小方框圖標(biāo)展示了作圖的十多種類型選項(xiàng)，選“Contour Plot”， 出現(xiàn)兩個(gè)選擇：“Contour”，即線條式等值線； “Filled Contour ”， 即色彩填充式等值線。選擇：“Contour”，經(jīng)過幾步選擇后即形成一個(gè)等值線表示的靜電勢等值線圖。進(jìn)一步選擇等值線線型的表達(dá)形式，可以找到上面分子圖的等值線的色帶樣式。得色彩等值線圖。（這個(gè)圖看似Filled Contour 的色彩填充式等值線，實(shí)應(yīng)是等值線的彩色線型。） 4、 就在得到的圖上面另行拷入一幅分子骨架圖，調(diào)整好比例，覆蓋在上面，二圖疊加，并令這上面圖層為“透明”，以便能看見下面的等值線圖。拖動(dòng)這骨架圖，使之如上面分子圖的位置。 這時(shí)會(huì)又調(diào)出等值線圖，選擇等值線的值，使之符合給定的范疇，這點(diǎn)還要推敲、實(shí)踐，或者需要剪裁使之適合范疇，上面這個(gè)圖有明顯的剪裁痕跡。因?yàn)樽髡咧幌肟催@一圈原子對(duì)內(nèi)的清晰靜電勢，這也不為過。圖即完成了，輸出作為插圖即可。 當(dāng)然，我這還是紙上談兵。關(guān)鍵是這樣的圖，是作了一個(gè)靜電勢截面等值線圖，鑲在分子骨架之中。 [ Last edited by zhou2009 on 2009-7-10 at 19:22 ] |
» 本帖已獲得的紅花（最新10朵）
小范范19893樓2009-07-10 15:55:26
zhaoj614
鐵桿木蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 6411.6
- 散金: 177
- 帖子: 281
- 在線: 214.9小時(shí)
- 蟲號(hào): 758060
- 注冊: 2009-04-26
- 性別: GG
- 專業(yè): 有機(jī)分子功能材料化學(xué)
10樓2009-08-11 09:17:27
CKX
木蟲 (著名寫手)
- 應(yīng)助: 45 (小學(xué)生)
- 金幣: 2261.1
- 散金: 133
- 紅花: 8
- 帖子: 1413
- 在線: 186.5小時(shí)
- 蟲號(hào): 602271
- 注冊: 2008-09-14
- 專業(yè): 理論和計(jì)算化學(xué)
★ ★ ★ ★ ★ ★ ★ ★
lei0736(金幣+3,VIP+0):謝謝 7-10 16:36
loovfnd(金幣+5,VIP+0):謝謝你的熱心幫助。 7-10 22:11
lei0736(金幣+3,VIP+0):謝謝 7-10 16:36
loovfnd(金幣+5,VIP+0):謝謝你的熱心幫助。 7-10 22:11
|
可以參閱《化學(xué)學(xué)報(bào)》 1992年05期 加入收藏 獲取最新 靜電勢方法研究分子間的電子轉(zhuǎn)移 Ⅰ.理論背景與方法曹軒 廖沐真 顧憲章 李林蜂 吳國是 【摘要】：光電轉(zhuǎn)換過程中起主導(dǎo)作用的分子間電子轉(zhuǎn)移可視為電子在分子間勢壘穿透的過程。我們在從頭算SCF-MO波函數(shù)基礎(chǔ)上將靜電勢方法用于確定分子間的勢壘,并由此出發(fā)從理論上得到了各種動(dòng)力學(xué)參數(shù)如電子穿透幾率和弛豫時(shí)間等。文中給出了詳細(xì)的數(shù)學(xué)公式和理論模型。 【作者單位】： 清華大學(xué)化學(xué)系 清華大學(xué)化學(xué)系 清華大學(xué)化學(xué)系 清華大學(xué)化學(xué)系 清華大學(xué)化學(xué)系 【關(guān)鍵詞】： 分子靜電勢 分子間電子轉(zhuǎn)移 勢壘高度 弛豫時(shí)間 電子給體 穿透幾率 電子穿透 波函數(shù) 光電轉(zhuǎn)換 電荷分離 【基金】：863基金資助的課題 【DOI】：CNKI:SUN:HXXB.0.1992-05-004 |
2樓2009-07-10 14:48:39
zhou2009
版主 (著名寫手)
- QC強(qiáng)帖: 15
- 應(yīng)助: 337 (大學(xué)生)
- 貴賓: 4.053
- 金幣: 42654.6
- 紅花: 153
- 帖子: 1729
- 在線: 1274.7小時(shí)
- 蟲號(hào): 676800
- 注冊: 2008-12-18
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
- 管轄: 量子化學(xué)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lei0736(金幣+6,VIP+0):謝謝 7-11 10:13
loovfnd(金幣+6,VIP+0):太感謝了！ 7-11 12:37
lei0736(金幣+6,VIP+0):謝謝 7-11 10:13
loovfnd(金幣+6,VIP+0):太感謝了！ 7-11 12:37
|
用G03作單點(diǎn)能計(jì)算，此時(shí)不需要用作CUBE的命令，用%chk=jobname得到CHK文件即成。再調(diào)用cubegen命令作分子的靜電勢空間格點(diǎn)文件。 用G03作單點(diǎn)能計(jì)算，也可用關(guān)鍵字cube=(potential,cards)，自己輸入格點(diǎn)數(shù)、起始點(diǎn)、格點(diǎn)間隔。即得分子的靜電勢空間格點(diǎn)文件。 用cubegen命令 1、在計(jì)算中，用%chk=jobname----------將計(jì)算存為：jobname.chk；默認(rèn)存于：scratch子文件夾。 2、調(diào)用Formchk 將二進(jìn)制的checkpoint 文件(chk)轉(zhuǎn)換成ASCII 碼形式(fchk)；得到j(luò)obname.fch 3、調(diào)用cubegen。 在一個(gè)三維網(wǎng)格（立方）點(diǎn)上求解：Orbitals分子軌道（ MO=n），potential靜電勢，Density電子密度，Gradient密度梯度，NormGradient密度梯度的范數(shù)和Laplacian電子密度的拉普拉斯算符(∇2ρ)。 默認(rèn)計(jì)算電子密度（對(duì)應(yīng)選項(xiàng)Density）。使用的密度由Density 關(guān)鍵字控制；若使用Density=Current 求解來自相關(guān)的或CI-Singles 波函的立方區(qū)域內(nèi)電子密度，而不是默認(rèn)的Hartree-Fock 密度。 Windows版cubegen的提示信息與回答： Property [Density]? MO=n ----------指定要選取的MO數(shù)。一旦回答了MO=n，就不是作Density而是作Orbitals。分子軌道數(shù)n、關(guān)鍵字Homo，Lumo，All，OccA（全部是alpha占據(jù)），OccB（全部是beta 占據(jù)），Valence（所有價(jià)軌道） 和 Virtuals（所有虛軌道）也可以用來代替特定的軌道數(shù)。n 值沒有默認(rèn)值，省略的話會(huì)造成錯(cuò)誤。 Formatted Checkpoint file? jobname ----------輸入的fchk文件的名字，寫路徑，但不寫后綴， 后綴默認(rèn)為.fch。 Cube file [ ]? Jobname ----------輸出的cube文件的名字，寫路徑，但不寫后綴， 后綴默認(rèn)為.cube。 Approximate points per side [0]? 回車 ------點(diǎn)陣每邊的點(diǎn)數(shù)。回車取0為默認(rèn)的80×80×80。 Header in cube file [H]? 回車 ---------輸出檔案的格式：H表示包含標(biāo)題（這是預(yù)設(shè)選擇）； GsGrid要求H�；卮餹，表示不要包含標(biāo)題。若要產(chǎn)生無格式方塊文件，這選項(xiàng)會(huì)被忽略。 如果是： Property [Density]? Potential 得到的即是分子的靜電勢空間格點(diǎn)文件。 |
4樓2009-07-11 10:02:54
5樓2009-07-11 19:38:10
ustc
金蟲 (著名寫手)
- 應(yīng)助: 9 (幼兒園)
- 金幣: 9239.4
- 散金: 207
- 紅花: 7
- 帖子: 1218
- 在線: 480.4小時(shí)
- 蟲號(hào): 28859
- 注冊: 2003-11-13
- 性別: GG
- 專業(yè): 催化化學(xué)
6樓2009-07-12 20:35:02
yjcmwgk
禁蟲 (文壇精英)
密度泛函·小卒
- QC強(qiáng)帖: 38
- 應(yīng)助: 165 (高中生)
- 貴賓: 16.999
- 金幣: 554.6
- 散金: 39880
- 紅花: 395
- 沙發(fā): 2
- 帖子: 10406
- 在線: 1528.1小時(shí)
- 蟲號(hào): 448368
- 注冊: 2007-11-01
- 專業(yè): 理論和計(jì)算化學(xué)
- 管轄: 科研工具資源
★ ★ ★ ★ ★ ★ ★
loovfnd(金幣+4,VIP+0):都是好方法，謝謝。 7-16 14:25
lei0736(金幣+3,VIP+0):謝謝 7-16 15:27
loovfnd(金幣+4,VIP+0):都是好方法，謝謝。 7-16 14:25
lei0736(金幣+3,VIP+0):謝謝 7-16 15:27
|
zhou2009大俠的跟帖挺好，呵呵，但是稍微有些麻煩 生成cube文件的方法，其實(shí)可以用一個(gè)命令就能達(dá)到目的。 zhou2009大俠生成MO和電子密度用了多步輸入，其實(shí)，實(shí)現(xiàn)cubegen的功能，只需要如下1個(gè)命令 cubegen 0 MO=n jobname.fchk Jobname.cube 0 h 其他cube文件，可以如下生成 cubegen 0 MO=n jobname.fchk Jobname.cube 0 h(計(jì)算某個(gè)分子軌道) cubegen 0 density jobname.fchk Jobname.cube 0 h(計(jì)算電子密度) cubegen 0 potential jobname.fchk jobname.cube 0 h(計(jì)算靜電勢) 等等各種cube文件 …… 這些命令均可在Linux界面下直接敲，也可以在dos界面下直接敲 [ Last edited by yjcmwgk on 2009-7-16 at 14:14 ] |
7樓2009-07-16 14:13:20
zhou2009
版主 (著名寫手)
- QC強(qiáng)帖: 15
- 應(yīng)助: 337 (大學(xué)生)
- 貴賓: 4.053
- 金幣: 42654.6
- 紅花: 153
- 帖子: 1729
- 在線: 1274.7小時(shí)
- 蟲號(hào): 676800
- 注冊: 2008-12-18
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
- 管轄: 量子化學(xué)
8樓2009-07-17 09:00:08
yjcmwgk
禁蟲 (文壇精英)
密度泛函·小卒
- QC強(qiáng)帖: 38
- 應(yīng)助: 165 (高中生)
- 貴賓: 16.999
- 金幣: 554.6
- 散金: 39880
- 紅花: 395
- 沙發(fā): 2
- 帖子: 10406
- 在線: 1528.1小時(shí)
- 蟲號(hào): 448368
- 注冊: 2007-11-01
- 專業(yè): 理論和計(jì)算化學(xué)
- 管轄: 科研工具資源
9樓2009-07-18 16:19:25
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 loovfnd 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 08工科 320總分 求調(diào)劑 +4 | 梨花珞晚風(fēng) 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 271求調(diào)劑 +12 | 生如夏花… 2026-03-11 | 14/700 |
|
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研]
 11408 一志愿西電，277分求調(diào)劑
+3
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
|
[考研] 341求調(diào)劑 +5 | 搗蛋豬豬 2026-03-11 | 5/250 |
|
|
|
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +5 | Liwangman 2026-03-15 | 5/250 |
|
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3 | 自由煎餅果子 2026-03-16 | 3/150 |
|
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-15 | 6/300 |
|
|
|
[考研] 【0703化學(xué)調(diào)劑】-一志愿華中師范大學(xué)-六級(jí)475 +5 | Becho359 2026-03-11 | 5/250 |
|
|
|
[考研] 材料371求調(diào)劑 +9 | 鱷魚? 2026-03-11 | 11/550 |
|
|
|
[考研] 0703，333分求調(diào)劑 一志愿鄭州大學(xué)-物理化學(xué) +3 | 李魔女斗篷 2026-03-11 | 3/150 |
|
|
|
[考研] 336求調(diào)劑 +6 | Iuruoh 2026-03-11 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +7 | 18880831720 2026-03-11 | 7/350 |
|
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
|
[考研] 277求調(diào)劑 +4 | anchor17 2026-03-12 | 4/200 |
|
|
|
[考研] 274求調(diào)劑0856材料化工 +12 | z2839474511 2026-03-11 | 13/650 |
|
|
|
[考博] 2026年博士申請(qǐng) +3 | QwQwQW10 2026-03-11 | 3/150 |
|
|
|
[考研] 420求調(diào)劑 +4 | 莫向外求11 2026-03-10 | 6/300 |
|