 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4966 | 回復(fù): 14 | ||||||||||||||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||||||||||||||
honey寶貝金蟲 (初入文壇)
|
[交流]
【轉(zhuǎn)載】Gaussian入門技巧六（關(guān)于虛頻） 已有14人參與
|
|||||||||||||
|
關(guān)于虛頻 首先，什么是頻率。 中學(xué)的時候我們學(xué)過簡諧振動，對應(yīng)的回復(fù)力是f=-kx,對應(yīng)的能量曲線，是一個開口向上的二次函數(shù)E=kx^2/2. 這樣的振動，對應(yīng)的x=0的點是能量極小值點（簡單情況下也就是最小值點）。這時的振動頻率我們也會求：ω=2π sqrt(k/m)。顯然它是一個正的頻率，也就是通常意義下的振動頻率。 那么，一維情況下，如果能量曲線是一個開口向下的二次曲線呢？首先，從能量上看，這是個不穩(wěn)定的點，中學(xué)的物理書上稱為“不穩(wěn)平衡”。用現(xiàn)在的觀點看，就是這一點導(dǎo)數(shù)是零（受力為0），且是能量極大值。如果套用上面的公式，“回復(fù)力”f=-k'x（實際上已經(jīng)不是回復(fù)，而是讓x越來越遠(yuǎn)了），這里k'是個負(fù)數(shù)，ω=2π sqrt(k'/m)顯然就是一個虛數(shù)了，即所謂的虛頻。Gaussian里面給出一個負(fù)的頻率，就是對應(yīng)這個虛頻的。 實際情況下，分子的能量是一個高維的勢能面，構(gòu)型優(yōu)化的時候，有時得到了極小值點，這樣這個點的任意方向上，都可以近似為開口向上的二次函數(shù)，這樣這里對應(yīng)的振動頻率就都是正的。對于極大值點，在每個方向都是開口向下的二次函數(shù)，那么頻率就會都是負(fù)的——當(dāng)然一般優(yōu)化很少會遇到這樣的情況。對于頻率有正有負(fù)的情況，說明找到的點在某些方向上是極大值，有些方向上是極小值。如果要得到穩(wěn)定的能量最低構(gòu)型，顯然需要通過微調(diào)分子的構(gòu)型，消去所有的虛頻。如何微調(diào)？要看虛頻的振動方向。想象著虛頻對應(yīng)的就是開口向下的二次函數(shù)，顯然，把分子坐標(biāo)按照振動的方向移動一點點，分子應(yīng)該就可以順著勢能面找到新的穩(wěn)定點，但是也不能太小。而所謂的過渡態(tài)，則是連接反應(yīng)物和產(chǎn)物之間的最低能量路徑上的能量極大值。好比山谷中的A，B兩點，它們之間的一個小土丘，就是過渡態(tài)，從A到B的反應(yīng)，需要越過的是這個小土丘，而不是兩邊的高山。這樣，過渡態(tài)就是在一個方向上是極大值，而在其它方向上都是極小值的點。因此，過渡態(tài)只有一個虛頻。 優(yōu)化得到虛頻，消虛頻的方法，就是根據(jù)虛頻對應(yīng)的振動位移(輸出文件每一個頻率下面都有一個類似坐標(biāo)的3列數(shù)值，給出的就是每個原子在這個振動頻率時候的移動方式)，把這些位移，乘以一個適當(dāng)大小的因子，直接和原來平衡構(gòu)型相加，就得到新的構(gòu)型了。 舉例說明： #HF/3-21G opt freq test 0 2 H 0.0 0.0 0.0 H 0.0 0.0 1.0 H 0.0 0.0 -1.0 上面的計算，是優(yōu)化一個 H-H-H的直線構(gòu)型然后算頻率。我們知道，3個H的穩(wěn)定構(gòu)型，應(yīng)該是一個H2分子，加上一個H原子。但是由于我們輸入的時候，中間H原子兩邊的H的鍵長相等，因此，在這個對稱性的限制下，結(jié)果給出的“穩(wěn)定構(gòu)型”是： H 0.000000 0.000000 0.000000 H 0.000000 0.000000 0.934091 H 0.000000 0.000000 -0.934091 頻率分析，有一個虛頻： 1 2 3 SGU PIU PIU Frequencies -- -2291.1909 1121.1486 1121.1486 Red. masses -- 1.0078 1.0078 1.0078 Frc consts -- 3.1172 0.7464 0.7464 IR Inten -- 42.9706 8.8167 8.8167 Raman Activ -- 0.0000 0.0000 0.0000 Depolar (P) -- 0.0000 0.0000 0.0000 Depolar (U) -- 0.0000 0.0000 0.0000 Atom AN X Y Z X Y Z X Y Z 1 1 0.00 0.00 0.82 0.82 0.03 0.00 -0.03 0.82 0.00 2 1 0.00 0.00 -0.41 -0.41 -0.02 0.00 0.02 -0.41 0.00 3 1 0.00 0.00 -0.41 -0.41 -0.02 0.00 0.02 -0.41 0.00 第一個是虛頻�？梢钥闯觯湔駝幽Ｊ绞牵褐虚g的H的z坐標(biāo)變大，兩邊的兩個H的z坐標(biāo)變小，（然后是中間的H的z坐標(biāo)變小，兩邊的H的z坐標(biāo)變大，完成一次振動），這個如果還不清楚，畫一下圖就知道了。 現(xiàn)在，我們就把原來的坐標(biāo)，加上振動模式對應(yīng)的坐標(biāo)： H 0.000000 0.000000 0.000000 H 0.000000 0.000000 0.934091 ＋ H 0.000000 0.000000 -0.934091 0.00 0.00 0.82 0.00 0.00 -0.41 0.00 0.00 -0.41 直接加的結(jié)果，是 H 0.000000 0.000000 0.820000 H 0.000000 0.000000 0.524091 H 0.000000 0.000000 -1.344091 畫圖知道，顯然對原來的構(gòu)型，變化太多了。對這個例子，如果取一個系數(shù)0.1乘以振動坐標(biāo)，再求和，那么結(jié)果就是： H 0.000000 0.000000 0.082000 H 0.000000 0.000000 0.893091 H 0.000000 0.000000 -0.975091 顯然就合理多了。用這個坐標(biāo)去重新優(yōu)化，就可以得到真正的沒有虛頻的穩(wěn)定結(jié)構(gòu)了。 這個系數(shù)0.1，只是個經(jīng)驗的數(shù)值，對于不同體系，可以自己設(shè)定。設(shè)置過小，會得到同樣的虛頻，設(shè)置過大，可能會得到別的構(gòu)型（當(dāng)然也可能碰巧得到更好的構(gòu)型）。 頻率分析只能在勢能面的穩(wěn)定點進(jìn)行，這樣，頻率分析就必須在已經(jīng)優(yōu)化好的結(jié)構(gòu)上進(jìn)行。頻率分析的另外一個用處是判斷穩(wěn)定點的本質(zhì)。穩(wěn)定點表述的是在勢能面上力為零的點，它即可能是極小值，也可能是鞍點。極小值在勢能面的各個方向都是極小的。而鞍點則是在某些方向上是極小的，但在某一個方向上是極大的，因為鞍點是連接兩個極小值的點。有關(guān)鞍點的信息: 1.負(fù)的頻率2.頻率相應(yīng)簡正振動的模式 當(dāng)一個結(jié)構(gòu)產(chǎn)生負(fù)的振動頻率時，可以表明在該振動方向可能存在著能量更低的結(jié)構(gòu)。判斷所得鞍點是不是需要的鞍點的方法，就是察看它的簡正振動模式，分析是不是可以導(dǎo)向所需要的產(chǎn)物或反應(yīng)物。進(jìn)一步的，更好的辦法是通過IRC計算來判斷反應(yīng)物，產(chǎn)物與得到的鞍點是否有關(guān)系。 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 高斯/MOPAC | 量化計算有意義的討論 | 量化 | Gaussian |
| computational chemistry （gaussian） | 自己感覺好的東西 | Gaussian | 論文圖表專輯 |
| 精華 | gaussian計算中常見問題 | Photonic Crystal |
» 本帖已獲得的紅花（最新10朵）
zi_cathy» 猜你喜歡
細(xì)胞培養(yǎng)，這22個細(xì)節(jié)一定要注意�。ㄒ唬�
已經(jīng)有0人回復(fù)
-
CSC訪學(xué)博后項目獲批，外方學(xué)校暫停合作該怎么辦？
已經(jīng)有51人回復(fù)
-
物理化學(xué)論文潤色/翻譯怎么收費?
已經(jīng)有94人回復(fù)
-
請問四氫呋喃溶解的聚合物用甲醇沉淀時，如何使沉淀過程加速？
已經(jīng)有2人回復(fù)
-
Chemical Bonding at Surfaces and Interfaces,最經(jīng)典的一本表面上化學(xué)相互作用教材
已經(jīng)有0人回復(fù)
-
分子動力學(xué)模擬合作需求
已經(jīng)有2人回復(fù)
-
誠邀有志之士加入江蘇大學(xué)環(huán)境學(xué)院合成生物學(xué)課題組！(長期接收學(xué)碩專碩、碩博連讀)
已經(jīng)有16人回復(fù)
-
一志愿南大理論與計算化學(xué)方向（070300）342分，求調(diào)劑
已經(jīng)有1人回復(fù)
-
歡迎化學(xué)、化工、物理、材料、或計算機(jī)背景的考研咨詢
已經(jīng)有0人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 筆記本運行Gaussian,裝linux系統(tǒng)能加快計算速度么
已經(jīng)有16人回復(fù)
- 使用chem3D調(diào)用gaussian
已經(jīng)有3人回復(fù)
- Gaussian中原子凈電荷的指定
已經(jīng)有3人回復(fù)
- 象棋入門六大名著-棋壇大家嘔心瀝血之作【原創(chuàng)】
已經(jīng)有28人回復(fù)
- gaussion09manual及gaussian入門（非09）
已經(jīng)有95人回復(fù)
- 【轉(zhuǎn)載】Gaussian入門技巧四（關(guān)于對稱性）
已經(jīng)有7人回復(fù)
- 【轉(zhuǎn)載】Gaussian入門技巧二（開殼層和閉殼層）
已經(jīng)有12人回復(fù)
- 【轉(zhuǎn)載】Gaussian入門技巧一（自旋多重度）
已經(jīng)有8人回復(fù)
honey寶貝
金蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1432.5
- 紅花: 4
- 帖子: 42
- 在線: 15.6小時
- 蟲號: 1163884
- 注冊: 2010-12-06
- 專業(yè): 理論和計算化學(xué)
3樓2011-04-06 09:12:01
zmcommon
木蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1685.3
- 紅花: 1
- 帖子: 145
- 在線: 53.6小時
- 蟲號: 990007
- 注冊: 2010-04-06
- 性別: MM
- 專業(yè): 理論和計算化學(xué)
正好用到，謝謝�。� 2樓2011-04-02 13:39:58
chujiajie
木蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1960.6
- 帖子: 107
- 在線: 66.7小時
- 蟲號: 485571
- 注冊: 2007-12-27
- 專業(yè): 理論和計算化學(xué)
5樓2011-06-28 17:35:54
keaiweiwei
銅蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 173.9
- 散金: 40
- 帖子: 102
- 在線: 36.4小時
- 蟲號: 1092208
- 注冊: 2010-09-07
- 性別: MM
- 專業(yè): 理論和計算化學(xué)
★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
|
再請教一下樓主，消除虛頻我之前看的說是把虛頻對應(yīng)的震動坐標(biāo)的1/2疊加到原始坐標(biāo)再優(yōu)化，在這里我怎么找虛頻對應(yīng)的震動坐標(biāo)？就是你舉的例子中的這一部分對應(yīng)的坐標(biāo)嗎？ Atom AN X Y Z X Y Z X Y Z 1 1 0.00 0.00 0.82 0.82 0.03 0.00 -0.03 0.82 0.00 2 1 0.00 0.00 -0.41 -0.41 -0.02 0.00 0.02 -0.41 0.00 3 1 0.00 0.00 -0.41 -0.41 -0.02 0.00 0.02 -0.41 0.00 新手，對這些都不是太清楚，謝謝啊。 |
8樓2012-07-31 10:18:22
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 281求調(diào)劑 +5 | 亞克西good 2026-03-26 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 生物學(xué) 296 求調(diào)劑 +5 | 朵朵- 2026-03-26 | 7/350 |
|
|
[考研] 085600，專業(yè)課化工原理，320分求調(diào)劑 +5 | 大饞小子 2026-03-29 | 5/250 |
|
|
|
[考研] 293求調(diào)劑 +3 | 末未mm 2026-03-30 | 5/250 |
|
|
|
[考研] 329求調(diào)劑，一志愿西北工業(yè)大學(xué)，材料工程（085601） +5 | 小小機(jī)靈蟲 2026-03-29 | 11/550 |
|
|
|
[基金申請] 面上5B能上會嗎？ +7 | redcom 2026-03-29 | 7/350 |
|
|
|
[考研] 生物技術(shù)與工程 +7 | 1294608413 2026-03-25 | 8/400 |
|
|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +25 | RichLi_ 2026-03-25 | 25/1250 |
|
|
|
[考研] 求化學(xué)調(diào)劑 +11 | wulanna 2026-03-28 | 11/550 |
|
|
|
[考研] 一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑 +6 | niko- 2026-03-27 | 6/300 |
|
|
|
[考研] 總分293求調(diào)劑 +8 | 加一一九 2026-03-25 | 11/550 |
|
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
|
[考研] 張芳銘-中國農(nóng)業(yè)大學(xué)-環(huán)境工程專碩-298 +4 | 手機(jī)用戶 2026-03-26 | 4/200 |
|
|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
|
[考研] 085600，材料與化工321分調(diào)劑 +4 | 大饞小子 2026-03-27 | 6/300 |
|
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
|
[考研] 各位老師您好：本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|