 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (6502)
- >考研 (3212)
- >導(dǎo)師招生 (2970)
- >文獻(xiàn)求助 (530)
- >蟲(chóng)友互識(shí) (446)
- >碩博家園 (296)
- >考博 (213)
- >招聘信息布告欄 (185)
- >博后之家 (123)
- >休閑灌水 (112)
- >論文投稿 (92)
- >綠色求助(高懸賞) (73)
- >基金申請(qǐng) (68)
- >找工作 (49)
- >SciFinder/Reaxys (39)
- >外文書(shū)籍求助 (37)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 823 | 回復(fù): 1 | |||
feynman1965金蟲(chóng) (正式寫(xiě)手)
|
[求助]
關(guān)于計(jì)算改變核間距時(shí)兩原子分子HOMO 已有1人參與
|
|
各位蟲(chóng)友好， 我是做實(shí)驗(yàn)的，現(xiàn)在想用高斯做一個(gè)這樣的計(jì)算，兩原子分子比如H2的某一個(gè)排斥態(tài)的HOMO。 這還沒(méi)完，分子從基態(tài)躍遷到這個(gè)排斥態(tài)肯定核間距不斷增加，直到解離。 怎么計(jì)算在這個(gè)解離過(guò)程中，不同核間距的HOMO? 請(qǐng)大家給個(gè)方向性和個(gè)別關(guān)鍵細(xì)節(jié)指導(dǎo)，詳細(xì)的我可以自己學(xué)習(xí)。 謝謝大家。 |
回復(fù)此樓» 猜你喜歡
深圳信息職業(yè)技術(shù)學(xué)院2026年教師招聘公告（最新）
已經(jīng)有0人回復(fù)
-
誠(chéng)邀有志之士加入江蘇大學(xué)環(huán)境學(xué)院合成生物學(xué)課題組！(長(zhǎng)期接收學(xué)碩專(zhuān)碩、碩博連讀)
已經(jīng)有16人回復(fù)
-
物理化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有145人回復(fù)
-
研究生招生
已經(jīng)有0人回復(fù)
-
一志愿南大理論與計(jì)算化學(xué)方向（070300）342分，求調(diào)劑
已經(jīng)有1人回復(fù)
-
歡迎化學(xué)、化工、物理、材料、或計(jì)算機(jī)背景的考研咨詢(xún)
已經(jīng)有0人回復(fù)
-
山西大學(xué)資源與環(huán)境工程研究所化工分離/環(huán)境功能材料小組招收2026年博士生
已經(jīng)有3人回復(fù)
-
一志愿南師大0703化學(xué) 275求調(diào)劑
已經(jīng)有5人回復(fù)
-
A 區(qū)過(guò)線(xiàn)就可申！重慶雙非公辦高校 0856調(diào)劑，接受跨專(zhuān)業(yè)！納米材料 / 潤(rùn)滑熱門(mén)方向
已經(jīng)有0人回復(fù)
-
歡迎調(diào)劑長(zhǎng)江師范學(xué)院湯老師課題組的碩士研究生
已經(jīng)有0人回復(fù)
-
調(diào)劑招生
已經(jīng)有3人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
雪狼乖乖 
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 18 (小學(xué)生)
- 貴賓: 0.005
- 金幣: 2805.8
- 散金: 567
- 紅花: 34
- 帖子: 901
- 在線(xiàn): 402.2小時(shí)
- 蟲(chóng)號(hào): 532524
- 注冊(cè): 2008-03-25
- 性別: GG
- 專(zhuān)業(yè): 理論和計(jì)算化學(xué)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
feynman1965: 金幣+20, ★★★很有幫助 2014-06-23 11:37:35
感謝參與，應(yīng)助指數(shù) +1
feynman1965: 金幣+20, ★★★很有幫助 2014-06-23 11:37:35
|
單重態(tài)H2拉伸過(guò)程的電子密度差等值線(xiàn)動(dòng)畫(huà) 轉(zhuǎn)載自Sobereva http://hi.baidu.com/sobereva/item/c69893149c743a6e70d5e8c3 做電子密度差圖必須先獲得涉及到的各個(gè)元素原子的波函數(shù)文件。每次執(zhí)行Multiwfn時(shí)程序會(huì)從當(dāng)前目錄下atomwfn目錄里尋找是否已有原子的波函數(shù)文件，沒(méi)有則需要調(diào)用Gaussian重算。為了方便、省時(shí)，這里先事先算好氫原子的波函數(shù)文件，放在atomwfn目錄下，并命名為H .wfn，注意空格必須有，否則Multiwfn不認(rèn)�；M建議與計(jì)算分子時(shí)一致（若基組級(jí)別已較高，不一致時(shí)問(wèn)題也不大）。其實(shí)計(jì)算密度差圖還有個(gè)既重要又麻煩的球?qū)ΨQ(chēng)化原子密度問(wèn)題，由于氫的密度本來(lái)就是球?qū)ΨQ(chēng)的，這里就不提了。 執(zhí)行壓縮包內(nèi)的H2sin.gjf，這是單重態(tài)H2的scan任務(wù)，從0.2埃開(kāi)始24步，每步0.1埃。注意盡管是單重態(tài)，但是拉遠(yuǎn)后對(duì)稱(chēng)破缺態(tài)是基態(tài)，所以要用非限制性方法，且需要加上guess(always,mix)打破對(duì)稱(chēng)性。計(jì)算完畢后將儲(chǔ)存波函數(shù)文件的目錄H2sin拷到Multiwfn所在目錄下。 編寫(xiě)一個(gè)文本文件batchcontourfix.txt放在Multiwfn目錄下，內(nèi)容是： 4 -2 1 2 200,200 6 //自行輸入格點(diǎn)設(shè)定 0.00000 -4.50000 -6.95664 //原點(diǎn) 0.00000 0.04500 0.00000 //平移向量1 0.00000 0.00000 0.06957 //平移向量2 0 [空行] [空行] 其中-2代表要做變形屬性，由于函數(shù)選了1，計(jì)算的就是變形密度（密度差）。第一節(jié)的例子中坐標(biāo)軸尺寸和刻度在不斷變動(dòng)，這是為了使當(dāng)前幀的結(jié)構(gòu)能充滿(mǎn)圖像，然而有些人不喜歡這樣，因此此例中自行設(shè)定原點(diǎn)和平移向量，坐標(biāo)軸就固定住了。顯然很難直觀判斷數(shù)值應(yīng)該設(shè)多少，我建議用體系涉及空間尺度最廣的一幀作參考獲取這些信息。比如此例H2_0024.wfn對(duì)應(yīng)拉得最遠(yuǎn)的一幀，就先手動(dòng)對(duì)這個(gè)波函數(shù)的任意一個(gè)函數(shù)做圖（X=0的YZ平面，延展距離用默認(rèn)的4.5 bohr），Multiwfn會(huì)輸出作圖時(shí)所用的格點(diǎn)信息： Origin is 0.00000 -4.50000 -6.95664 bohr The x/y/z of transition vector 1 is 0.00000 0.04500 0.00000 bohr The x/y/z of transition vector 2 is 0.00000 0.00000 0.06957 bohr 令所有幀作圖時(shí)都使用這個(gè)格點(diǎn)設(shè)定，就能保證每一幀感興趣的區(qū)域都不出界。 然后，將前面batchrun.bat里的H2tri替換為H2sin，batchcolor.txt替換為batchcontourfix.txt，執(zhí)行之，就得到了每幀的.PNG圖像。將圖像復(fù)制到/sob/H2sinpic目錄下，并進(jìn)入/sob目錄，執(zhí)行: convert -delay 20 -resize 500 -shave '100x0' -colors 128 -monitor -reverse -fill red -annotate +5+395 "Created by Sobereva" H2sinpic/DISLIN.PNG H2sinpic/DISLIN_?.PNG H2sinpic/DISLI_??.PNG H2sinpic/H2sin_all.gif 和第一節(jié)的命令沒(méi)什么區(qū)別，只是目錄換了，并且加入了一個(gè)新的設(shè)定-shave '100x0'，這是因?yàn)檫@次的圖像的坐標(biāo)軸橫軸固定了，左右的白邊長(zhǎng)度也固定了，將每幀在resize為寬度為500后，左右都截掉100像素可以讓動(dòng)畫(huà)的畫(huà)布更緊湊。最終所得動(dòng)畫(huà)H2sin_all.gif如圖所示： 隨距離拉近后，電子密度迅速在原子間聚集，表明共價(jià)鍵正在加強(qiáng)。聚集的區(qū)域一開(kāi)始逐漸變大并逐漸變得扁長(zhǎng)，但后來(lái)卻開(kāi)始收縮并變得接近于球型，最后一幀電子密度減少區(qū)域和增加區(qū)域都幾乎變成了球?qū)ΨQ(chēng)狀，前者包著后者。這是因?yàn)樵有驍?shù)越大，核吸引勢(shì)越強(qiáng)，電子密度分布范圍收縮得越厲害，基態(tài)He的1s電子密度延展范圍小于基態(tài)H的，所以當(dāng)兩個(gè)氫核距離接近到氦核特征明顯出現(xiàn)時(shí)內(nèi)部電子密度就增加了，相應(yīng)地外側(cè)減少。 e824b899a9014c08b18022600a7b02087af4f4b8.jpg.gif |
2樓2014-06-23 11:36:20
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 feynman1965 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[電化學(xué)] 070300化學(xué)調(diào)劑 +7 | 山頂見(jiàn)α 2026-03-25 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 一志愿西安交大材料學(xué)碩（英一數(shù)二）347，求調(diào)劑到高分子/材料相關(guān)專(zhuān)業(yè) +4 | zju51 2026-03-31 | 6/300 |
|
|
|
[考研] 080200專(zhuān)業(yè)277分，求帶走！ +4 | 瓶子PZ 2026-03-31 | 4/200 |
|
|
|
[考研] 一志愿南師大0703化學(xué) 275求調(diào)劑 +5 | Ripcord上岸 2026-03-27 | 5/250 |
|
|
|
[考研] 材料科學(xué)與工程求調(diào)劑 +13 | 深V宿舍吧 2026-03-29 | 13/650 |
|
|
|
[考研] 070300化學(xué)專(zhuān)業(yè)279調(diào)劑 +9 | 哈哈哈^_^ 2026-03-31 | 9/450 |
|
|
|
[考研] 070300一志愿211，312分求調(diào)劑院校 +12 | 小黃鴨寶 2026-03-30 | 12/600 |
|
|
|
[考研] 080500-315分復(fù)試調(diào)劑 +9 | 上岸3821 2026-03-31 | 9/450 |
|
|
|
[考研] 282求調(diào)劑 不挑專(zhuān)業(yè) 求收留 +4 | Yam. 2026-03-30 | 5/250 |
|
|
|
[考研] 生物學(xué) 296 求調(diào)劑 +7 | 朵朵- 2026-03-26 | 9/450 |
|
|
|
[考研] 085600 295分求調(diào)劑 +13 | W55j 2026-03-30 | 15/750 |
|
|
|
[考研] 一志愿中海洋320化學(xué)工程與技術(shù)學(xué)碩求調(diào)劑 +8 | 披星河 2026-03-30 | 8/400 |
|
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +11 | 多久上課 2026-03-27 | 12/600 |
|
|
|
[考研] 材料化工340求調(diào)劑 +3 | jhx777 2026-03-30 | 3/150 |
|
|
|
[考研] 298求調(diào)劑 +3 | 什么是胖頭魚(yú) 2026-03-30 | 5/250 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
|
[考研] 085701環(huán)境工程，267求調(diào)劑 +16 | minht 2026-03-26 | 16/800 |
|
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲(chóng)不懶了 2026-03-27 | 5/250 |
|
|
|
[考研] 機(jī)械學(xué)碩310分，數(shù)一英一，一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
|
[考研] 環(huán)境專(zhuān)碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|