| 5 | 1/1 | 返回列表 |
| 查看: 2509 | 回復: 8 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
zzyuss鐵桿木蟲 (知名作家)
抗旱救災委員會名譽會長
|
[求助]
求助Gaussian09計算HOMO-LUMO 和 三重態(tài),單重態(tài) △E ST 已有1人參與
|
||
|
各位蟲友,本人gaussian計算小白,最近嘗試用gaussian09計算一個化合物單體的homo-lumo,單線態(tài)能級分布,三線態(tài)能級分布,以及系間竄越的能差Δes╟t值。在此基礎(chǔ)上繼續(xù)進行這個化合物單體聚集形成二聚體、三聚體以及四聚體的homo-lumo,單線態(tài)能級分布,三線態(tài)能級分布,以及系間竄越的能差Δes╟t值。化合物單體含有一個碘代季銨鹽,請問大家該如何進行計算?大概的步驟是什么?有沒有大神能幫個忙?有償?shù)模?br />
另外,我的化合物單體能夠形成70個有效相干長度的聚集體,請問各位能不能計算70個聚集體的homo-lumo,單線態(tài)能級分布,三線態(tài)能級分布,以及系間竄越的能差Δes╟t值呢? [ Last edited by paramecium86 on 2020-8-10 at 16:26 ] |

版主 (著名寫手)
|
首先高斯一般情況下(20-40核心左右的服務(wù)器,linux版的高斯)用普通dft搭配最最保底的基組對于基態(tài)結(jié)構(gòu)的優(yōu)化,最多可以算300-400個原子數(shù)目的體系。所以看你體系的原子數(shù)目可以大概判斷出能不能算得了。而且因為可能要涉及到了激發(fā)態(tài)的優(yōu)化,這樣的話200個原子以上就基本算不動了。所以70個分子聚在一起是肯定算不動的。 單體如果比較大可能連四聚體都不那么好算。畢竟聚在一起多了可能還有不同的構(gòu)型。 以下是計算的簡單過程,我默認樓主知道如何用高斯進行,幾何優(yōu)化計算, 并且知道dft, tddft都是什么。知道激發(fā)態(tài)是什么。如果這些都不知道的話,最好先去入門一下。(具體如何入門該看什么資料。見我這個教程前三頁http://m.gaoyang168.com/t-14242735-1) 就單體而言,計算homo-lumo軌道其實就是對這個分子的基態(tài)進行結(jié)構(gòu)優(yōu)化之后,就可以查看所有分子軌道的信息了,包括各個軌道的能量或者要作圖也可以看分子軌道的等值面圖;蛘呖赡苓涉及到homo-lumo軌道重合度分析。就都是這一步就可以完成。 至于系間竄越就涉及到了激發(fā)態(tài)計算,F(xiàn)在比較流行的方法就是tddft了。(簡單的介紹可看http://m.gaoyang168.com/t-14242735-1教程第1.5.1節(jié)和1.5.4 節(jié)推薦的幾篇博文。) 而涉及到系間竄越的速率k isc 一般用Marcus公式來計算竄越速率什么的。這個公式里的 △E ST 一般要用到垂直能量差(如果是TADF之類的要計算反系間竄越 k risc 有時候要用到絕熱能量差) 垂直能差計算比較容易 大體步驟就是 優(yōu)化好分子的基態(tài)S0結(jié)構(gòu)之后。用優(yōu)化好的結(jié)構(gòu)計算tddft(50-50) 關(guān)鍵詞來計算基態(tài)結(jié)構(gòu)下的各個單重和三重激發(fā)態(tài)的能量。把S1 T1能量作差就是 垂直 △E ST 。在這個計算中因為無論哪個激發(fā)態(tài),都是基于基態(tài)的結(jié)構(gòu),它們結(jié)構(gòu)相同,所以叫做垂直能量差。 而絕熱能差,既然叫做絕熱,就需要 對S1 和 T1 分別都進行幾何優(yōu)化(顯然S1 和 T1的最穩(wěn)定結(jié)構(gòu)肯定不一樣)。這樣S1 和T1 各自都處于不同結(jié)構(gòu)下,這樣計算出來的能量差就叫絕熱能量差。 對S1 的優(yōu)化可以使用 tddft opt 關(guān)鍵詞進行優(yōu)化。計算完成后可以得到s1的電子能。(用tddft計算后得到的激發(fā)態(tài)的電子能如何尋找請看我這個教程http://m.gaoyang168.com/t-14242735-1最后一頁) 對T1的優(yōu)化同樣可以使用tddft(triplet) opt來進行優(yōu)化。計算完成后也可以得到t1的電子能。二者相減就是 絕熱能量差。 對于二聚體 三聚體 四聚體 計算方法基本一樣。只是聚在一起可能會有不同構(gòu)型。只要順利找到優(yōu)化出最穩(wěn)定的構(gòu)型。剩下的步驟就都一樣。 |
鐵桿木蟲 (知名作家)
抗旱救災委員會名譽會長

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料085600 303求調(diào)劑 +5 | 1bygone 2026-03-04 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 0703化學 一志愿東北大學 279分求調(diào)劑 +4 | 聶樂平 2026-03-03 | 4/200 |
|
|
[考研] 0703 東華大學 理學化學均過 a 線,總分 281 求調(diào)劑 +3 | 香香 Lu 2026-03-03 | 5/250 |
|
|
[考研] 本科太原理工采礦工程,求調(diào)劑 +3 | onlx 2026-03-01 | 3/150 |
|
|
[考研] 接收調(diào)劑 +15 | 津萌津萌 2026-03-02 | 23/1150 |
|
|
[考研] 0703化學 學碩 理工科均可 不區(qū)分研究方向 總分279求調(diào)劑 +5 | 1一11 2026-03-03 | 5/250 |
|
|
[考博] 26申博 求博導 +3 | 愛讀書的小帥 2026-02-28 | 5/250 |
|
|
[考研] 298求調(diào)劑一志愿中海洋 +3 | lour. 2026-03-03 | 3/150 |
|
|
[碩博家園] 2025屆雙非化工碩士畢業(yè),申博 +4 | 更多的是 2026-02-27 | 5/250 |
|
|
[基金申請] 沒有青基直接申請面上,感覺自己瘋了 +5 | kevin63t 2026-03-02 | 6/300 |
|
|
[考研] 求調(diào)劑 +4 | Guo_yuxuan 2026-03-02 | 5/250 |
|
|
[考研] 材料學碩318求調(diào)劑 +11 | February_Feb 2026-03-01 | 11/550 |
|
|
[考研] 321求調(diào)劑一志愿東北林業(yè)大學材料與化工英二數(shù)二 +5 | 蟲蟲蟲蟲蟲7 2026-03-01 | 9/450 |
|
|
[考研] 295求調(diào)劑。一志愿報考鄭州大學化學工藝學碩,總分295分 +8 | yl1 2026-03-02 | 9/450 |
|
|
[考研] 0856材料調(diào)劑 +5 | 沿岸有貝殼OUC 2026-03-02 | 5/250 |
|
|
[考研] 0856材料求調(diào)劑 +12 | hyf hyf hyf 2026-02-28 | 13/650 |
|
|
[考研] 一志愿華南理工大學材料與化工326分,求調(diào)劑 +3 | wujinrui1 2026-02-28 | 3/150 |
|
|
[碩博家園] 博士自薦 +7 | 科研狗111 2026-02-26 | 11/550 |
|
|
[考研]
|
LYidhsjabdj 2026-02-28 | 4/200 |
|
|
[論文投稿]
Optics letters投稿被拒求助
30+3
|
luckyry 2026-02-26 | 4/200 |
|