Li+與EC(碳酸乙烯酯)結(jié)合能計(jì)算
本人開(kāi)始學(xué)習(xí)計(jì)算，想計(jì)算Li+與不同溶劑分子的結(jié)合能
比如要計(jì)算Li+與EC的結(jié)合能，直接可以從gaussian view讀出charge（1） 和spin（single）值，使用什么method和 basis set比較好，新手還望大佬多賜教，最后能給出簡(jiǎn)單的input文件，萬(wàn)分感謝，求大佬救小弟與水火之中
返回小木蟲(chóng)查看更多
今日熱帖
京公網(wǎng)安備 11010802022153號(hào)
結(jié)合能本身沒(méi)有什么特別的。就是用結(jié)合之后整體的能量 分別減去 Li+ 和 EC它們單獨(dú)的能量即可。這里的能量是看你自己的需要需要哪種，就計(jì)算哪種。電子能（純理論分析） 普通的結(jié)構(gòu)優(yōu)化就可以得到， 跟實(shí)驗(yàn)的鍵”能”對(duì)應(yīng)的一般是 要計(jì)算分子的焓 ，這就需要用頻率計(jì)算得到熱力學(xué)校正量。 有些情況也會(huì)用到結(jié)合自由能，這個(gè)也是熱力學(xué)量量。也是由頻率計(jì)算給出。
至于用什么計(jì)算的方法。 幾何優(yōu)化和頻率計(jì)算 可以就用B3LYP/6-311G(d) EM=GD3BJ 至于后續(xù)要不要用更高級(jí)別的計(jì)算方法去計(jì)算更精確的電子能就看你個(gè)人需求。
但是從你的描述來(lái)看，你高斯都還沒(méi)有入門，還沒(méi)有到需要選擇計(jì)算方法和基組的地步。我上面說(shuō)的如果也都還看不懂在說(shuō)什么。那就趕緊入門一下高斯。（入門高斯，幾乎是“入門量子化學(xué)計(jì)算”的必經(jīng)之路）
怎么入門可以看 http://gaoyang168.com/t-14242735-1 我這個(gè)教程的前言，當(dāng)然教程內(nèi)容本身也是入門的一些內(nèi)容，但是完全都還不會(huì)用高斯的時(shí)候看這個(gè)教程本身用處不大。
官方的教程 我這個(gè)帖子 http://gaoyang168.com/t-14151335-1 有分享。 建議是跟著第二版 堅(jiān)持看上前幾章 跟著里面的練習(xí)和給的例子文件，入門高斯肯定沒(méi)啥問(wèn)題。 如果傾向于看中文，民間縮減版的中文的第二版也可以。入門也還夠。第三版雖然例子更多也有官方中文版不過(guò)里面紕漏比第二版多得多。新手容易被帶偏了。
至于輸入文件可以給一個(gè) 不過(guò)在真正入門前用處也不大的。入門之后這種類型的計(jì)算自己也能信手拈來(lái)的
比如一個(gè)分子的幾何優(yōu)化和頻率計(jì)算
%nproc=16
%mem=20GB
%chk=output_explain_opt_freq.chk
#p opt freq b3lyp/6-311G(D) EM=GD3BJ
h2o opt freq
0 1
O -5.61694298 0.21178637 0.00000000
H -4.65694298 0.21178637 0.00000000
H -5.93739757 1.11672220 0.00000000，
我以前有計(jì)算過(guò)分子的LUMO和HOMO值，麻煩再問(wèn)一下，問(wèn)什么加 EM=GD3BJ，這個(gè)是必須的嗎？
EM=GD3BJ 是給當(dāng)前的泛函加入色散校正。 彌補(bǔ)一些比較老的泛函（比如 B3LYP）沒(méi)有考慮色散作用的缺點(diǎn)。 一般來(lái)說(shuō)這個(gè)校正加上算出來(lái)結(jié)果會(huì)好一些。而且也不增加計(jì)算量。所以可以常規(guī)帶著這個(gè)關(guān)鍵詞。 尤其是最近這一兩年，帶著色散校正成為主流了。 不過(guò)倒是加或者不加對(duì)于HOMO LUMO來(lái)說(shuō)影響應(yīng)該是不大。
祝福