 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4162 | 回復(fù): 9 | |||||||||
| 【獎(jiǎng)勵(lì)】 本帖被評(píng)價(jià)6次，作者sobereva增加金幣 4.8 個(gè) | |||||||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||||||||
sobereva版主
|
[資源]
談?wù)勚C振頻率校正因子 已有2人參與
|
||||||||
|
談?wù)勚C振頻率校正因子 通常理論計(jì)算的頻率是諧振頻率，由于忽視了非諧振效應(yīng)，而且理論方法本身也有誤差，所以直接基于它計(jì)算基頻、零點(diǎn)能、焓、熵都可能有較大誤差。而將頻率校正因子（后文簡(jiǎn)稱校正因子）乘在理論計(jì)算出的諧振頻率上，基頻以及熱力學(xué)校正量就變得較為準(zhǔn)確了，往往能減小2、3倍誤差。例如一般雜化泛函算出的基頻誤差能到100多cm-1，而校正后一般只有幾十cm-1的誤差。由于忽略了非諧振效應(yīng)會(huì)高估頻率，因此基頻的校正因子一般都小于1.0。 校正因子分為多種，通常有 (1)用于獲得準(zhǔn)確基頻的校正因子。通常對(duì)于所有頻率都是相同的校正值，但有些人為了校正得更準(zhǔn)確，把頻率劃分成不同范圍（比如高頻和低頻范圍），分別給出不同的校正因子 (2)用于獲得準(zhǔn)確諧振頻率的校正因子。注意諧振頻率不是實(shí)驗(yàn)可觀測(cè)的而只是個(gè)概念，但可以通過(guò)一定手段將實(shí)驗(yàn)數(shù)據(jù)轉(zhuǎn)化得到 (3)用于獲得準(zhǔn)確的ZPE的校正因子 (4)用于獲得準(zhǔn)確的焓變?chǔ)=H(T)-H(0)的校正因子。注：H(0)=U(0)=電子能量+ZPE (5)用于獲得準(zhǔn)確的熵的校正因子 這5種校正因子各不相同，不要隨意混用。不過(guò)(1)/(3)以及(2)/(3)的數(shù)值對(duì)各種情況倒是基本總是固定的，分別為0.974和1.014。對(duì)HF或DFT，(4)和(5)的數(shù)值比較接近，混用也無(wú)妨。ΔH和熵的數(shù)值是依賴于計(jì)算時(shí)設(shè)定的溫度的，其校正因子也因此受溫度影響。 由于諧振子模型下ZPE是頻率的線性函數(shù)，即ZPE=∑hν/2，所以第(3)類校正因子乘到原始頻率上還是乘到基于原始頻率算出的ZPE上效果都一樣。但是對(duì)于(4)和(5)類校正因子，要分別乘到原始頻率上再算ΔH和熵，而切不要用原始頻率算完了之后再對(duì)所得的ΔH和S乘上校正因子。 還有一點(diǎn)要說(shuō)明，由于ZPE=∑hν/2，表面上看仿佛只要對(duì)振動(dòng)頻率ν乘上校正因子成為準(zhǔn)確基頻后ZPE就應(yīng)該準(zhǔn)確了，因此(1)和(2)類校正因子仿佛應(yīng)當(dāng)是相同的，這是誤解！因?yàn)檫@種ZPE的計(jì)算公式是對(duì)于諧振模型而言的，既然校正成了實(shí)際的非諧振頻率，那么當(dāng)然ZPE也得按照非諧振子來(lái)算才行。嚴(yán)格的非諧振的ZPE必須直接做非諧振計(jì)算才能得到。ZPE的校正因子則實(shí)際上是先乘到諧振頻率上將之轉(zhuǎn)化成一種特殊頻率，然后再按照Z(yǔ)PE=∑hν/2來(lái)算的時(shí)候，得到的結(jié)果就能正好和實(shí)際的ZPE很接近。 校正因子一般都是用先人擬合好的，通過(guò)向相應(yīng)類型實(shí)驗(yàn)數(shù)據(jù)擬合來(lái)得到。用于擬合的分子少則幾種，多則上百種。常用的校正因子可以在這些地方得到： 1 http://comp.chem.umn.edu/freqscale/version3b2.htm。這是Truhlar組擬合的一大堆校正因子，相當(dāng)全面，包括前述的(1)、(2)、(3)類校正因子。不過(guò)很遺憾的是沒有給出ΔH和熵的校正因子。 2 http://cccbdb.nist.gov/vibscalejust.asp。這是Computational Chemistry Comparison and Benchmark DataBase(CCCBDB)數(shù)據(jù)庫(kù)中的一頁(yè)，給出的是基頻的校正因子。 3 J. Phys. Chem., 100, 16502-16513 (1996)中的表10。給出的是前述(1)、(3)、(4)、(5)類校正因子，列舉的計(jì)算級(jí)別較少。注意表10的ΔH和S的校正因子是298.15K下的，其它溫度下的見文中圖2。 雖然這三種來(lái)源擬合校正因子過(guò)程所用的測(cè)試集和方法有異，但對(duì)于同種級(jí)別，結(jié)果是接近的，因此它們是有可比性的。還有一些校正因子零零碎碎地在一些文獻(xiàn)中報(bào)道，用得較少，可靠性難說(shuō)。 對(duì)于HF和各種GGA泛函，基頻的校正因子分別接近0.9和1.0。這即是說(shuō)GGA泛函由于理論方法誤差和非諧振效應(yīng)抵消，不做校正算出來(lái)的基頻也挺準(zhǔn)確。而HF由于顯著高估力常數(shù)，所以算出來(lái)的頻率太高，算出來(lái)的3000cm-1實(shí)際上只有約2700cm-1而已。雜化泛函的基頻校正因子則介于0.9~1.0之間，而且往往HF成分越高數(shù)值越小。 注意振動(dòng)頻率和熱力學(xué)校正量的計(jì)算絕對(duì)不是用的理論級(jí)別越高越好。CCSD(T)/cc-pVTZ可能還不如B3LYP/6-31G*算出來(lái)的頻率準(zhǔn)確，因?yàn)檎`差既來(lái)自于理論方法給出的力常數(shù)的誤差，也來(lái)自于忽略了非諧振效應(yīng)的誤差。在某些低水平下，可能恰好這兩部分抵消得好，導(dǎo)致結(jié)果不錯(cuò)，GGA泛函算的基頻幾乎不需要校正就是這個(gè)原因。而如果盲目地用很高的計(jì)算級(jí)別，雖然基本消除了理論方法本身的誤差，但非諧振效應(yīng)缺乏較好的抵消，那么所得的基頻以及熱力學(xué)校正量很可能比很爛的級(jí)別算出來(lái)的還糟。另外還有關(guān)鍵的一點(diǎn)，就是在使用了校正因子后，不同方法間的差距會(huì)很大程度上拉平，原本差的方法在校正后結(jié)果可能很好。這正是為什么熱力學(xué)組合方法G3，雖然目標(biāo)是QCISD(T)的精度，卻采用HF/6-31G*（經(jīng)過(guò)校正因子校正）這樣表面看上去如此爛的級(jí)別來(lái)算ZPE。所以說(shuō)，用很高級(jí)別的方法去做振動(dòng)分析是⑨的行為，何況計(jì)算二階導(dǎo)數(shù)還那么費(fèi)時(shí)，純屬白耽誤功夫。半經(jīng)驗(yàn)方法也有對(duì)應(yīng)的校正因子，但是即便是經(jīng)過(guò)校正后，誤差還是較大，大于從頭算方法很多，所以除非迫不得已不宜用半經(jīng)驗(yàn)方法算頻率。基于分子力場(chǎng)也可以算頻率和熱力學(xué)校正量，如果分子力場(chǎng)選得合適，那么結(jié)果可能很不錯(cuò)，且不需要做校正，例如MMFF94和MM3對(duì)于典型有機(jī)分子基頻平均誤差在60cm-1左右。 通常來(lái)說(shuō)，使用B3LYP/6-31G*結(jié)合校正因子就能得到很好的計(jì)算結(jié)果，計(jì)算量小，而且比起其它多數(shù)方法經(jīng)過(guò)校正后的結(jié)果更好。雖然Truhlar的文章表明VSXC經(jīng)過(guò)校正后更準(zhǔn)確一點(diǎn)，但畢竟是相對(duì)小眾向泛函，寡人不推薦。BMC-CCSD經(jīng)過(guò)校正后特別準(zhǔn)確，但畢竟比較貴，除非是要做高精度計(jì)算否則就不推薦了。 各類校正因子對(duì)于所用基組都不敏感（除非是和很小的基組相比），所以如果你用的基組沒有對(duì)應(yīng)的校正因子，可以直接用此方法在相近的基組下的校正因子。對(duì)于同一種方法，無(wú)論是校正前還是校正后，都不代表越大的基組下結(jié)果越準(zhǔn)確，而通常是差不多（除非基組爛到3-21G這種層次），因此根本沒必要用雙zeta以上的基組、帶彌散函數(shù)做振動(dòng)分析。 要提醒的是，振動(dòng)分析所用的級(jí)別和優(yōu)化幾何結(jié)構(gòu)的級(jí)別必須嚴(yán)格一致，而得到的熱力學(xué)校正量是可以加到其它更高水平方法算的電子能量上的。比如說(shuō)，要算一個(gè)體系某溫度下的氣相下的自由能，那么可以用很好的方法諸如PWPB95/def2-QZVP去優(yōu)化氣相結(jié)構(gòu)并獲得電子能量E，然后再用B3LYP/6-31G*去優(yōu)化體系然后在考慮校正因子的情況下做振動(dòng)分析，得到ZPE、ΔH和熵。最后體系的自由能G=E+ZPE+ΔH-T*S。 雖然校正后的B3LYP/6-31G*是通常推薦的做振動(dòng)分析的級(jí)別，但是有些體系，可能在B3LYP/6-31G*下面優(yōu)化的結(jié)果很糟，諸如弱相互作用顯著的情況，進(jìn)而可能由此導(dǎo)致算出的頻率和熱力學(xué)校正量都很差，盡管已經(jīng)經(jīng)過(guò)校正。此時(shí)就應(yīng)該選擇更合適的方法做振動(dòng)分析以及對(duì)應(yīng)的優(yōu)化，例如可以用M06-2X、wB97XD。一個(gè)有點(diǎn)麻煩的問(wèn)題是諸如M06-2X、wB97XD這樣的較新的泛函缺乏ΔH和熵的校正因子，不過(guò)從JPC,100,16502的數(shù)據(jù)來(lái)看，雜化泛函的這兩類校正因子都比較接近于1.0，所以不做校正也沒什么問(wèn)題�；蛘吒鶕�(jù)已有的DFT的校正因子去猜也成，例如Truhlar給出的wB97XD/ma-TZVP和JPC,100,16502給出的B3P86/6-31G*的基頻和ZPE校正因子都差不多，因此可以假設(shè)wB97XD/ma-TZVP的ΔH和熵的校正因子也與之相仿佛，可以挪來(lái)使用。 注意校正因子并非總能很好解決振動(dòng)問(wèn)題，盡管通�？梢浴Ｓ眯Ｕ蜃雍笾皇钦w統(tǒng)計(jì)結(jié)果變好了，并不代表對(duì)每種體系都能很好處理，比如很多級(jí)別下對(duì)于氟分子的振動(dòng)頻率算得都不咋樣。并且對(duì)于非諧振效應(yīng)很強(qiáng)、很低頻率的情況，用了校正因子后結(jié)果肯定也還是不好，只能直接去做非諧振計(jì)算才行。 下面簡(jiǎn)單說(shuō)一下在Gaussian中使用校正因子的方式。在freq任務(wù)過(guò)程中使用scale=X就代表校正因子為X，例如 #P HF/6-31G* freq scale=0.9135 temperature=323 注意輸出信息中只是在計(jì)算熱力學(xué)校正量的時(shí)候?qū)λ妙l率乘上了校正因子，而輸出的頻率還是原始的頻率！所以要想得到校正后的頻率，得自己手動(dòng)把頻率乘上校正因子。Gaussian顯然是不會(huì)自動(dòng)用校正因子的，即默認(rèn)scale=1.0。 如果已經(jīng)做過(guò)一次頻率計(jì)算并且保存了chk文件，想查看另一個(gè)校正因子下的熱力學(xué)數(shù)據(jù)，最好的方法是用Gaussian自帶的freqchk工具，它會(huì)從chk文件中讀取儲(chǔ)存的Hessian矩陣并作分析。啟動(dòng)它之后，輸入chk文件路徑，然后依次輸入 N //不輸出Hyperchem文件 323 //溫度。如果輸入0則為默認(rèn)的298.15K 1 //壓力。如果輸入0則為默認(rèn)的1atm 0.8905 //校正因子 [回車] [回車] 然后立刻就看到了原始頻率以及等同于scale=0.8905 temperature=323時(shí)freq任務(wù)給出的熱力學(xué)數(shù)據(jù)，和Gaussian自身輸出的一樣。 利用freqchk，就可以反復(fù)利用不同校正因子得到較準(zhǔn)確的ZPE、ΔH和熵，并獲得自由能G=H-T*S，過(guò)程在下面說(shuō)一下。以下是freq任務(wù)會(huì)輸出的四個(gè)量，為了省事用字母表示，使用不同校正因子時(shí)它們都會(huì)有變化。 Zero-point correction= A Thermal correction to Energy= B Thermal correction to Enthalpy= C 注意這一項(xiàng)不是ΔH，而是ZPE+ΔH Thermal correction to Gibbs Free Energy= D 首先，在ZPE校正因子下得到ZPE（即A）。然后在熵校正因子下得到-T*S（即D-C）。之后在ΔH校正因子下得到ΔH（即C-A）。最終將這三個(gè)量加在一起，再加上高精度方法算的電子能量E，就得到了自由能G=E+ZPE+ΔH-T*S。 在作IR或Raman光譜圖時(shí)，要先把頻率乘上校正因子后再做圖。手動(dòng)一一去改Gaussian輸出文件里的頻率值顯然太費(fèi)事了。Multiwfn(http://multiwfn.codeplex.com)的繪制光譜圖的功能可以方便地設(shè)定校正因子。啟動(dòng)Multiwfn后，載入g03或g09的freq任務(wù)的輸出文件，選11進(jìn)入光譜繪制功能，選IR或Raman，之后選0就可以顯示光譜圖。如果選14 Multiply the vibrational frequencies by a factor，就可以令不同編號(hào)范圍的頻率乘以指定的校正因子，重新作圖后就能看到效果了。用戶因此不僅可以對(duì)所有頻率乘上統(tǒng)一的校正因子，還可以根據(jù)一些文章的思想，對(duì)于不同范圍的頻率乘上不同的校正因子以得到更準(zhǔn)確結(jié)果。 最后順帶一提的是，除了對(duì)頻率進(jìn)行校正，還有人提出了直接對(duì)力常數(shù)進(jìn)行校正的方法，不過(guò)這樣的方法涉及到太多參數(shù)，難搞，所以沒流行起來(lái)。對(duì)于NMR、UV/Vis等其它譜的計(jì)算也有人提出過(guò)不同級(jí)別下的各種校正方法，不過(guò)都遠(yuǎn)沒有像頻率校正因子這樣被廣泛使用。 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 量化軟件學(xué)習(xí) | 量化學(xué)習(xí) | SOB 集錦 | Gaussian |
| 氣態(tài)分子化學(xué)動(dòng)力學(xué) | raman | 計(jì)算 |
» 猜你喜歡
PbS量子點(diǎn)如何確定濃度
已經(jīng)有11人回復(fù)
-
200年來(lái)人們一直無(wú)法破解的難題： 氣液臨界點(diǎn)狀態(tài)方程
已經(jīng)有7人回復(fù)
-
物理化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有53人回復(fù)
-
關(guān)于引力延遲下雙星系統(tǒng)角動(dòng)量不守恒的疑惑---求高手解答！
已經(jīng)有10人回復(fù)
-
大灣區(qū)大學(xué)劉天輝課題組招收2026級(jí)博士生（與中山大學(xué)物理學(xué)院聯(lián)合培養(yǎng)）
已經(jīng)有22人回復(fù)
-
PbS量子點(diǎn)紫外吸收
已經(jīng)有0人回復(fù)
-
比利時(shí)魯汶大學(xué)與國(guó)家留學(xué)基金委共同資助博士研究生CSC-KU Leuven PhD Scholarship
已經(jīng)有0人回復(fù)
-
半透性容器進(jìn)行高濕或者低濕研究
已經(jīng)有0人回復(fù)
-
倫敦大學(xué)學(xué)院Benjamin Moss博士招收人工智能、原位光譜、催化相關(guān)26年博士生
已經(jīng)有11人回復(fù)
-
中國(guó)科學(xué)院成都山地災(zāi)害與環(huán)境研究所2026年特別研究助理（博士后）招聘啟事
已經(jīng)有0人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 雜質(zhì)的校正因子是咋回事��？
已經(jīng)有11人回復(fù)
- 校正因子和原研不一致，怎么辦
已經(jīng)有16人回復(fù)
- 校正因子的影響因素
已經(jīng)有5人回復(fù)
- HPLC有關(guān)物質(zhì)方法學(xué)驗(yàn)證校正因子問(wèn)題
已經(jīng)有19人回復(fù)
- 氣相色譜內(nèi)標(biāo)法測(cè)校正因子，內(nèi)標(biāo)物的選取有什么要求？校正因子可能會(huì)小于1嗎？
已經(jīng)有5人回復(fù)
- 外標(biāo)法和相對(duì)校正因子加自身對(duì)照法對(duì)比計(jì)算雜質(zhì)含量結(jié)果偏差比較大-已解決
已經(jīng)有12人回復(fù)
- 待測(cè)組分對(duì)于內(nèi)標(biāo)物的相對(duì)質(zhì)量校正因子
已經(jīng)有3人回復(fù)
- 誰(shuí)了解新藥報(bào)批中用校正因子方法測(cè)定有關(guān)物質(zhì)
已經(jīng)有11人回復(fù)
- DMM 的FID校正因子？
已經(jīng)有4人回復(fù)
- 菜鳥求助色譜標(biāo)準(zhǔn)曲線與校正因子的關(guān)系
已經(jīng)有6人回復(fù)
- 【求助】GC測(cè)氣體校正因子
已經(jīng)有9人回復(fù)
- 【求助】用色譜校正因子求轉(zhuǎn)化率的公式
已經(jīng)有15人回復(fù)
- 【求助】求詢B3LYP/6-311++G(3df,3pd)的校正因子
已經(jīng)有6人回復(fù)
- 【求助】新人，請(qǐng)教怎樣做零點(diǎn)能校正？
已經(jīng)有17人回復(fù)
- 【討論】光譜還可以校正？
已經(jīng)有8人回復(fù)
7樓2014-03-23 01:03:39
2樓2014-01-19 13:50:58
3樓2014-01-19 14:15:26
6樓2014-01-20 08:33:09
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +14 | 丹青奶蓋 2026-03-26 | 16/800 |
|
|---|---|---|---|---|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-04-01 | 5/250 |
|
|
|
[考研] 英一數(shù)二生物信息學(xué)287分，本科生物科學(xué)，求調(diào)劑 +3 | 碧水xyz 2026-03-29 | 4/200 |
|
|
|
[考研] 材料調(diào)劑 +6 | 一樣YWY 2026-04-01 | 6/300 |
|
|
|
[考研] 085600，材料與化工321分調(diào)劑 +5 | 大饞小子 2026-03-27 | 7/350 |
|
|
|
[考研] 362求調(diào)劑 +11 | 西南交材料專碩3 2026-03-31 | 11/550 |
|
|
|
[考研] 322求調(diào)劑 +8 | 三水sss 2026-04-01 | 8/400 |
|
|
|
[考研] 考研生物與醫(yī)藥調(diào)劑 +7 | 鐵憨憨123425 2026-03-31 | 7/350 |
|
|
|
[考研] 一志愿a區(qū)211，085601-307分求調(diào)劑 +10 | 黨嘉豪 2026-03-31 | 23/1150 |
|
|
|
[考研] 0856 335分 | 封裝or激光加工 老師看看我��！ +9 | cccchenso 2026-03-29 | 9/450 |
|
|
|
[基金申請(qǐng)] 面上5B能上會(huì)嗎？ +8 | redcom 2026-03-29 | 8/400 |
|
|
|
[考研] 286求調(diào)劑 +6 | Faune 2026-03-30 | 6/300 |
|
|
|
[考研] 0703化學(xué) +20 | 妮妮ninicgb 2026-03-27 | 20/1000 |
|
|
|
[考研] 313求調(diào)劑 +6 | 賣個(gè)關(guān)子吧 2026-03-31 | 6/300 |
|
|
|
[考研] 085600，專業(yè)課化工原理，320分求調(diào)劑 +6 | 大饞小子 2026-03-29 | 6/300 |
|
|
|
[考研] 327求調(diào)劑 +5 | 小卡不卡. 2026-03-29 | 5/250 |
|
|
|
[考研] 0703 化學(xué) 求調(diào)劑，一志愿山東大學(xué) 342 分 +7 | Shern—- 2026-03-28 | 7/350 |
|
|
|
[考研] 環(huán)境科學(xué)與工程334分求調(diào)劑 +6 | 王一一依依 2026-03-30 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +10 | 張zz111 2026-03-27 | 11/550 |
|
|
|
[考研] 295求調(diào)劑 +5 | wei-5 2026-03-26 | 5/250 |
|