 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 2756 | 回復(fù): 5 | |||||||||
| 【獎勵】 本帖被評價4次，作者sobereva增加金幣 3.2 個 | |||||||||
sobereva至尊木蟲 (著名寫手)
|
[資源]
利用Multiwfn計算Gaussian輸出的激發(fā)態(tài)之間的躍遷偶極矩 已有1人參與
|
||||||||
|
利用Multiwfn計算Gaussian輸出的激發(fā)態(tài)之間的躍遷偶極矩 Gaussian的CIS、TDHF、TDDFT可以直接給出基態(tài)到激發(fā)態(tài)的躍遷偶極矩。CIS也提供了alltransitiondensities關(guān)鍵詞可以把激發(fā)態(tài)之間的躍遷偶極矩輸出出來，但是對于TD任務(wù)則沒有這個關(guān)鍵詞。雖然Gaussian有個density=transition=(n,m)關(guān)鍵詞可以產(chǎn)生n到m激發(fā)態(tài)的躍遷密度，然而這個關(guān)鍵詞在目前的版本中貌似功能不正常，而且光是有躍遷密度也沒用，還得利用偶極矩積分才能得到躍遷偶極矩，但Gaussian也并不給出這樣的信息。 筆者多次被問及怎么基于Gaussian得到激發(fā)態(tài)之間的躍遷偶極矩，對于TD任務(wù)確實(shí)沒有可行的辦法。于是筆者在Multiwfn多功能波函數(shù)分析程序的電子激發(fā)分析模塊中加入了這一功能，這里就介紹下怎么使用。這里用的是Multiwfn 3.3版。如果對Multiwfn一無所知的話，建議先閱讀《Multiwfn入門tips》（http://hi.baidu.com/sobereva/item/896ee3a19f6d7d3c020a4d76） Multiwfn計算激發(fā)態(tài)間的躍遷偶極矩需要兩個文件：(1)Gaussian的CIS或TDHF或TDDFT的輸出文件 (2)相應(yīng)任務(wù)的.fch文件。 Gaussian的輸入文件對于CIS、TDHF、TDDFT，分別寫成類似這樣 # B3LYP/6-31+G(d) TD(nstates=10) IOp(9/40=4) # CIS(nstates=10)/6-31+G(d) IOp(9/40=4) # HF/6-31+G(d) TD(nstates=10) IOp(9/40=4) 這里假設(shè)算10個激發(fā)態(tài)。Multiwfn計算時需要利用Gaussian輸出的激發(fā)態(tài)的組態(tài)系數(shù)，默認(rèn)情況下只有大于0.1的系數(shù)會被輸出出來，而較小的都不輸出，這樣的話Multiwfn算出的結(jié)果將會不太準(zhǔn)確。IOp(9/40=x)的含義是將系數(shù)大于10^-x的組態(tài)都輸出出來，因此IOp(9/40=4)會把系數(shù)絕對值大于0.0001的組態(tài)都輸出。x設(shè)的越大，輸出的組態(tài)越多，結(jié)果越準(zhǔn)確，但是x太大的話Multiwfn的耗時也會非常大，通常x=3就夠了，精度足夠，計算速度也比較快。x=4就已經(jīng)很精確了。 在計算之后得到了.out/.log文件，還同時得到了.chk文件。用formchk將.chk轉(zhuǎn)換為.fch，這個文件里記錄了基函數(shù)的定義以及基態(tài)軌道信息，這是Multiwfn要利用的。 下面以苯酚為例進(jìn)行說明怎么利用Multiwfn計算激發(fā)態(tài)之間的躍遷偶極矩，輸入文件如下，計算后得到phenol.out以及phenol.fch。 %chk=C:\gtest\phenol.chk # b3lyp/6-31+G* TD(nstates=5) IOp(9/40=3) b3lyp/6-311G** opted 0 1 C 0.01810200 -1.86802400 0.00000000 C 1.23175000 -1.16732400 0.00000000 C 1.23175000 0.23407600 0.00000000 C 0.01810200 0.93477600 0.00000000 C -1.19554600 0.23407600 0.00000000 C -1.19554600 -1.16732400 0.00000000 H 0.01810200 -2.93802400 0.00000000 H 2.15839700 -1.70232400 0.00000000 H 2.15839700 0.76907600 0.00000000 H -2.12219300 0.76907600 0.00000000 H -2.12219300 -1.70232400 0.00000000 O 0.01810200 2.36477600 0.00000000 H -0.88699500 2.68477600 0.00000000 啟動Multiwfn，依次輸入以下命令，//后的是注釋。 c:\gtest\phenol.fch //先載入電子激發(fā)任務(wù)產(chǎn)生的fch文件 18 //電子激發(fā)分析功能，包含多個子功能。這些功能無與倫比的強(qiáng)大，建議參看Multiwfn手冊4.18節(jié)的例子以及3.21節(jié)的原理和細(xì)節(jié)的介紹來學(xué)習(xí)使用 5 //計算激發(fā)態(tài)間的所有躍遷偶極矩 c:\gtest\phenol.out //電子激發(fā)任務(wù)的Gaussian輸出文件 此時屏幕上首先輸出了激發(fā)態(tài)的匯總信息 Exc.state# Exc.energy(eV) Multi. N_pairs Sum coeff.^2 1 5.06970 1 449 0.500004 2 5.47020 1 102 0.499997 3 5.94190 1 82 0.499975 4 6.26640 1 80 0.499979 5 6.38330 1 83 0.499978 其中N_pairs代表這個激發(fā)態(tài)在輸出文件中通過多少組態(tài)來表示，IOp(9/40=x)的x越大顯然N_pairs也就越大。Multi.是激發(fā)態(tài)的自旋多重度。Sum coeff.^2越接近理想值說明結(jié)果精度越高，對于閉殼層和開殼層情況理想值分別是0.5和1.0。如果偏離理想值比較大，則應(yīng)該加大x來保證結(jié)果精度。 然后選1，結(jié)果就輸出到了屏幕上。也可以選2輸出到文本文件里。結(jié)果如下： Transition dipole moment between excited states (a.u.): i j X Y Z Ene.diff.(eV) Oscil.str 1 1 -0.510929 0.541413 0.000000 0.00000 0.00000 1 2 0.000000 0.000000 0.180185 0.40050 0.00032 1 3 0.000000 0.000000 -0.052359 0.87220 0.00006 1 4 0.000000 0.000000 0.033679 1.19670 0.00003 1 5 0.000000 0.000000 -0.035481 1.31360 0.00004 2 2 2.353643 -3.851667 0.000000 0.00000 0.00000 2 3 0.110576 1.848714 0.000000 0.47170 0.03964 2 4 0.029210 -0.159245 0.000000 0.79620 0.00051 2 5 1.825329 0.672088 0.000000 0.91310 0.08464 3 3 -0.434476 2.840406 0.000000 0.00000 0.00000 3 4 -0.437467 0.064725 0.000000 0.32450 0.00155 3 5 3.994461 -0.441695 0.000000 0.44140 0.17466 4 4 2.503675 -4.112800 0.000000 0.00000 0.00000 4 5 0.489590 -0.576211 0.000000 0.11690 0.00164 5 5 -0.695139 0.745470 0.000000 0.00000 0.00000 每一對兒激發(fā)態(tài)間的躍遷偶極矩的三個分量，彼此間的能量差，以及相應(yīng)的振子強(qiáng)度都輸出出來了。 值得一提的是，對于編號相同的輸出，得到的躍遷偶極矩<m|r|n>顯然就等價于這個激發(fā)態(tài)自身的偶極矩<n|r|n>。比如3 3，如上可見Multiwfn的結(jié)果為-0.434476 2.840406 0.000000，換算成Debye為單位就是-1.104329 7.219591 0.000000。Gaussian直接輸出的3號激發(fā)態(tài)的偶極矩為 X= -1.1033 Y= 7.2217 Z= 0.0000 可見相符得很好。同時這也表明IOp(9/40=3)的精度足夠了，沒太大必要用到IOp(9/40=4)。 PS：在Gaussian中輸出激發(fā)態(tài)偶極矩是使用例如下面這樣的關(guān)鍵詞，在輸出文件末尾就會由L601模塊輸出第3個激發(fā)態(tài)的偶極矩 #P b3lyp/6-31+G* TD(root=3,nstates=5) density=rhoci 注意這里用了rhoci關(guān)鍵詞，如果只寫density的話，那么Gaussian傳遞給L601模塊的激發(fā)態(tài)的密度是弛豫的密度。而當(dāng)用rhoci時，傳遞的是非弛豫的密度，這才是和通過組態(tài)系數(shù)和基態(tài)軌道直接產(chǎn)生的密度直接對應(yīng)的。 [ Last edited by sobereva on 2014-4-14 at 19:08 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 SOB 集錦 | 理論化學(xué)計算知識 | 計算 | 科研之寶，所向披靡 |
| 高斯專輯 | QC Learning | 計算化學(xué) |
» 猜你喜歡
全獎博士 英國利物浦大學(xué) × 臺灣清華大學(xué) 聯(lián)合培養(yǎng)
已經(jīng)有0人回復(fù)
-
尋博導(dǎo)
已經(jīng)有13人回復(fù)
-
物理化學(xué)論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有132人回復(fù)
-
基金正文30頁指的是報告正文還是整個申請書
已經(jīng)有9人回復(fù)
-
哈爾濱工業(yè)大學(xué)（深圳）綠色化學(xué)團(tuán)隊(duì)招收2026年秋季入學(xué)博士生
已經(jīng)有11人回復(fù)
-
澳大利亞昆士蘭科技大學(xué)PhD博士獎學(xué)金
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 激發(fā)態(tài)的躍遷密度（TD）和電荷差異密度(CDD)
已經(jīng)有15人回復(fù)
- multiwfn繪制的ESP等值面圖里靜電勢的默認(rèn)單位是什么？
已經(jīng)有8人回復(fù)
- 求助：如何用Gaussian計算激發(fā)態(tài)偶極矩？
已經(jīng)有6人回復(fù)
- gauss 計算電偶極矩
已經(jīng)有22人回復(fù)
- 使用Multiwfn結(jié)合VMD分析和繪制分子表面靜電勢分布
已經(jīng)有14人回復(fù)
- 一起北京看美女吧：量子化學(xué)波函數(shù)分析暨M(jìn)ultiwfn程序使用培訓(xùn)班
已經(jīng)有18人回復(fù)
- Multiwfn讀取NBO的.40文件是對應(yīng)MO嗎？
已經(jīng)有10人回復(fù)
- 關(guān)于G09計算激發(fā)態(tài)偶極矩，基態(tài)偶極矩，躍遷偶極矩的討論
已經(jīng)有16人回復(fù)
- 建議multiwfn開發(fā)GUI界面 文件打開什么的效率太低了
已經(jīng)有10人回復(fù)
- 如何用Multiwfn軟件做出P DOS?
已經(jīng)有4人回復(fù)
- 使用Multiwfn繪制原子軌道圖形、研究原子殼層結(jié)構(gòu)及相對論效應(yīng)的影響
已經(jīng)有17人回復(fù)
- 使用Multiwfn繪制NBO及相關(guān)軌道
已經(jīng)有30人回復(fù)
- 引用了multiwfn的論文�，F(xiàn)在multiwfn軟件已經(jīng)有19次引用率了呢
已經(jīng)有13人回復(fù)
- 使用Multiwfn做拓?fù)浞治鲆约坝嬎愎聦﹄娮咏嵌?/a>
已經(jīng)有18人回復(fù)
- 回答一些關(guān)于Multiwfn的疑問以及未來Multiwfn的發(fā)展打算
已經(jīng)有37人回復(fù)
2樓2014-03-28 22:12:40
sob姐的帖子，支持！ 3樓2014-03-29 16:19:58
ggdh
專家顧問 (正式寫手)
-
專家經(jīng)驗(yàn): +69 - QC強(qiáng)帖: 2
- 應(yīng)助: 52 (初中生)
- 金幣: 1747.4
- 帖子: 697
- 在線: 206.1小時
- 蟲號: 96968
4樓2014-03-30 13:17:03
簡單回復(fù)
chrinide5樓
2014-03-31 07:54
回復(fù)
五星好評 
blueybz6樓
2014-04-01 09:16
回復(fù)
五星好評
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 sobereva 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 268求調(diào)劑 +7 | 好運(yùn)連綿不絕 2026-03-12 | 8/400 |
|
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
|
[考研] 344求調(diào)劑 +4 | knight344 2026-03-16 | 4/200 |
|
|
|
[考研] 312求調(diào)劑 +4 | 陌宸希 2026-03-16 | 5/250 |
|
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎兩個實(shí)用型專利一項(xiàng)省級立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +6 | SUICHILD 2026-03-12 | 6/300 |
|
|
|
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時間點(diǎn) 2026-03-13 | 5/250 |
|
|
|
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +5 | Liwangman 2026-03-15 | 5/250 |
|
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 快樂的香蕉 2026-03-11 | 4/200 |
|
|
|
[考研] 310求調(diào)劑 +3 | 【上上簽】 2026-03-11 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
|
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3 | 土木小天才 2026-03-12 | 3/150 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-03-12 | 4/200 |
|