求助Gaussian09計(jì)算HOMO-LUMO 和 三重態(tài),單重態(tài) △E ST
各位蟲友，本人gaussian計(jì)算小白，最近嘗試用gaussian09計(jì)算一個(gè)化合物單體的homo-lumo，單線態(tài)能級分布，三線態(tài)能級分布，以及系間竄越的能差Δes╟t值。在此基礎(chǔ)上繼續(xù)進(jìn)行這個(gè)化合物單體聚集形成二聚體、三聚體以及四聚體的homo-lumo，單線態(tài)能級分布，三線態(tài)能級分布，以及系間竄越的能差Δes╟t值�；衔飭误w含有一個(gè)碘代季銨鹽，請問大家該如何進(jìn)行計(jì)算？大概的步驟是什么？有沒有大神能幫個(gè)忙��？有償?shù)模?br />
另外，我的化合物單體能夠形成70個(gè)有效相干長度的聚集體，請問各位能不能計(jì)算70個(gè)聚集體的homo-lumo，單線態(tài)能級分布，三線態(tài)能級分布，以及系間竄越的能差Δes╟t值呢？
[ Last edited by paramecium86 on 2020-8-10 at 16:26 ]
返回小木蟲查看更多
今日熱帖
京公網(wǎng)安備 11010802022153號(hào)
單體計(jì)算的話，你自己的機(jī)器還跑得動(dòng)。
聚集體的話，理論來說我覺得也可以跑哈哈，但這七十個(gè)單體聚合，自己的電腦得跑爆炸。
我想嘗試自己跑一下，但是對于高斯還不會(huì)用，所以想請教一下大家?guī)兔χ更c(diǎn)一下
各位大神幫幫忙
首先高斯一般情況下（20-40核心左右的服務(wù)器，linux版的高斯）用普通dft搭配最最保底的基組對于基態(tài)結(jié)構(gòu)的優(yōu)化，最多可以算300-400個(gè)原子數(shù)目的體系。所以看你體系的原子數(shù)目可以大概判斷出能不能算得了。而且因?yàn)榭赡芤婕暗搅思ぐl(fā)態(tài)的優(yōu)化，這樣的話200個(gè)原子以上就基本算不動(dòng)了。所以70個(gè)分子聚在一起是肯定算不動(dòng)的。
單體如果比較大可能連四聚體都不那么好算。畢竟聚在一起多了可能還有不同的構(gòu)型。
以下是計(jì)算的簡單過程，我默認(rèn)樓主知道如何用高斯進(jìn)行，幾何優(yōu)化計(jì)算， 并且知道dft， tddft都是什么。知道激發(fā)態(tài)是什么。如果這些都不知道的話，最好先去入門一下。（具體如何入門該看什么資料。見我這個(gè)教程前三頁http://www.gaoyang168.com/t-14242735-1）
就單體而言，計(jì)算homo-lumo軌道其實(shí)就是對這個(gè)分子的基態(tài)進(jìn)行結(jié)構(gòu)優(yōu)化之后，就可以查看所有分子軌道的信息了，包括各個(gè)軌道的能量或者要作圖也可以看分子軌道的等值面圖�；蛘呖赡苓�涉及到homo-lumo軌道重合度分析。就都是這一步就可以完成。
至于系間竄越就涉及到了激發(fā)態(tài)計(jì)算。現(xiàn)在比較流行的方法就是tddft了。（簡單的介紹可看http://www.gaoyang168.com/t-14242735-1教程第1.5.1節(jié)和1.5.4 節(jié)推薦的幾篇博文。）
而涉及到系間竄越的速率k isc 一般用Marcus公式來計(jì)算竄越速率什么的。這個(gè)公式里的 △E ST 一般要用到垂直能量差（如果是TADF之類的要計(jì)算反系間竄越 k risc 有時(shí)候要用到絕熱能量差）
垂直能差計(jì)算比較容易
大體步驟就是 優(yōu)化好分子的基態(tài)S0結(jié)構(gòu)之后。用優(yōu)化好的結(jié)構(gòu)計(jì)算tddft(50-50) 關(guān)鍵詞來計(jì)算基態(tài)結(jié)構(gòu)下的各個(gè)單重和三重激發(fā)態(tài)的能量。把S1 T1能量作差就是 垂直 △E ST 。在這個(gè)計(jì)算中因?yàn)闊o論哪個(gè)激發(fā)態(tài)，都是基于基態(tài)的結(jié)構(gòu)，它們結(jié)構(gòu)相同，所以叫做垂直能量差。
而絕熱能差，既然叫做絕熱，就需要 對S1 和 T1 分別都進(jìn)行幾何優(yōu)化（顯然S1 和 T1的最穩(wěn)定結(jié)構(gòu)肯定不一樣）。這樣S1 和T1 各自都處于不同結(jié)構(gòu)下，這樣計(jì)算出來的能量差就叫絕熱能量差。
對S1 的優(yōu)化可以使用 tddft opt 關(guān)鍵詞進(jìn)行優(yōu)化。計(jì)算完成后可以得到s1的電子能。（用tddft計(jì)算后得到的激發(fā)態(tài)的電子能如何尋找請看我這個(gè)教程http://www.gaoyang168.com/t-14242735-1最后一頁）
對T1的優(yōu)化同樣可以使用tddft(triplet) opt來進(jìn)行優(yōu)化。計(jì)算完成后也可以得到t1的電子能。二者相減就是 絕熱能量差。
對于二聚體 三聚體 四聚體 計(jì)算方法基本一樣。只是聚在一起可能會(huì)有不同構(gòu)型。只要順利找到優(yōu)化出最穩(wěn)定的構(gòu)型。剩下的步驟就都一樣，
非常有幫助，謝謝！另外，我想問一個(gè)問題，就是我的二聚體、三聚體、四聚體都是非共價(jià)相互作用聚集到一起的，我想知道的是怎么在高斯里導(dǎo)入這種非共價(jià)相互作用結(jié)合的分子，并且可以拿來計(jì)算優(yōu)化？或者說是怎么在高斯里模擬這種結(jié)構(gòu)進(jìn)行計(jì)算優(yōu)化？
高斯里也不區(qū)分共價(jià)還是非共價(jià)（相互作用比較弱就把各個(gè)單體相互擺的遠(yuǎn)一點(diǎn)）。你把初始結(jié)構(gòu)設(shè)定好 算出來該是什么樣子就就是什么樣子。算出來之后 想知道它們之間是什么作用，可以進(jìn)行進(jìn)一步的波函數(shù)分析。但是算之前不需要考慮它們之間是什么相互作用。