| 24小時(shí)熱門(mén)版塊排行榜 |
| 查看: 1673 | 回復(fù): 10 | |||||
[交流]
【求助】gaussian中簡(jiǎn)正坐標(biāo)的輸出 已有2人參與
|
|
在一篇文章《Vibrational Analysis in Gaussian》中說(shuō)到“The purpose of this document is to describe how Gaussian calculates the reduced mass, frequencies, force constants, and normal coordinates which are printed out at the end of a frequency calculation.”從這句話看normal coordinates應(yīng)該在計(jì)算頻率的時(shí)候打印出來(lái)了啊,可是不知道哪地方給的是normal coordinates?以H2o為例,計(jì)算頻率時(shí)的部分輸出為: (Enter /home/gau03/g03/l716.exe) Dipole =-3.44478646D-17-3.96506834D-17-6.72452773D-01 Polarizability= 4.00613449D-02 3.18281559D-16 5.50808430D+00 1.58802971D-16 1.97181875D-16 2.56587227D+00 HyperPolar =-3.59200065D-17-6.27316867D-18-2.96845704D-16 -1.10949459D-14 7.62420726D-02-3.00230546D-16 -7.01237208D+00-1.56024452D-16 2.91713493D-15 -4.11420102D+00 Full mass-weighted force constant matrix: Low frequencies --- -6.0514 -1.0843 -0.0019 -0.0012 -0.0009 4.5788 Low frequencies --- 2170.0372 4140.0019 4391.0737 Diagonal vibrational polarizability: 0.0000000 0.0418305 0.1108974 Harmonic frequencies (cm**-1), IR intensities (KM/Mole), Raman scattering activities (A**4/AMU), depolarization ratios for plane and unpolarized incident light, reduced masses (AMU), force constants (mDyne/A), and normal coordinates:(這里也顯示normal coordinates,但是不知道怎么看,我怎么找不到啊? 1 2 3 A1 A1 B2 Frequencies -- 2170.0372 4140.0019 4391.0737 Red. masses -- 1.0785 1.0491 1.0774 Frc consts -- 2.9923 10.5941 12.2392 IR Inten -- 7.2379 44.2862 29.9710 Atom AN X Y Z X Y Z X Y Z 1 1 0.00 0.45 0.54 0.00 -0.57 0.42 0.00 0.54 -0.45 2 1 0.00 -0.45 0.54 0.00 0.57 0.42 0.00 0.54 0.45 3 8 0.00 0.00 -0.07 0.00 0.00 -0.05 0.00 -0.07 0.00 謝謝指點(diǎn)。 |
量化 | Gaussian |
至尊木蟲(chóng) (著名寫(xiě)手)
本人已永久離開(kāi)小木蟲(chóng)
|
Atom AN X Y Z X Y Z X Y Z 1 1 0.00 0.45 0.54 0.00 -0.57 0.42 0.00 0.54 -0.45 2 1 0.00 -0.45 0.54 0.00 0.57 0.42 0.00 0.54 0.45 3 8 0.00 0.00 -0.07 0.00 0.00 -0.05 0.00 -0.07 0.00 這個(gè)就是。例如2號(hào)原子的X,Y,Z坐標(biāo)構(gòu)成第1種簡(jiǎn)正坐標(biāo)的系數(shù)為0.00 -0.45 0.54。 你把這里每個(gè)振動(dòng)模式的3*3矩陣排列成9*1就容易看出簡(jiǎn)正坐標(biāo)與原先3N坐標(biāo)的關(guān)系了。如第一個(gè)簡(jiǎn)正坐標(biāo): 0.00 0.45 0.54 0.00 -0.45 0.54 0.00 0.00 -0.07 [ Last edited by sobereva on 2010-8-1 at 16:18 ] |
至尊木蟲(chóng) (著名寫(xiě)手)
本人已永久離開(kāi)小木蟲(chóng)
| 簡(jiǎn)正振動(dòng)坐標(biāo)和原始3N個(gè)笛卡爾坐標(biāo)是一個(gè)線性變換的關(guān)系,前者可寫(xiě)為后者的線性組合。每個(gè)振動(dòng)模式對(duì)應(yīng)一簡(jiǎn)正坐標(biāo),由于簡(jiǎn)正坐標(biāo)在每個(gè)原子笛卡爾坐標(biāo)上都有分量,所以振動(dòng)自然會(huì)體現(xiàn)在原子笛卡爾坐標(biāo)的變化上。 |
至尊木蟲(chóng) (著名寫(xiě)手)
本人已永久離開(kāi)小木蟲(chóng)
至尊木蟲(chóng) (著名寫(xiě)手)
本人已永久離開(kāi)小木蟲(chóng)
|
若不考慮質(zhì)權(quán)、笛卡爾-內(nèi)坐標(biāo)轉(zhuǎn)換這些細(xì)節(jié)問(wèn)題,L就是Hessian矩陣的本征向量矩陣,也就是Vibrational Analysis in Gaussian的公式7。 關(guān)于變換關(guān)系的細(xì)節(jié)問(wèn)題在Vibrational Analysis in Gaussian里面都詳細(xì)說(shuō)明。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703 化學(xué) 求調(diào)劑,一志愿山東大學(xué) 342 分 +5 | Shern—- 2026-03-28 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 一志愿北京工業(yè)大學(xué),324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
[考研] 343求調(diào)劑 +6 | 愛(ài)羈絆 2026-03-29 | 6/300 |
|
|
[考研] 332求調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
[考研] 本科雙非材料,跨考一志愿華電085801電氣,283求調(diào)劑,任何專業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
[考研] 322求調(diào)劑 +7 | 宋明欣 2026-03-27 | 7/350 |
|
|
[考研] 283求調(diào)劑 +3 | A child 2026-03-28 | 3/150 |
|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +12 | 路癡小琪 2026-03-27 | 12/600 |
|
|
[考研] 312,生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分,英語(yǔ)已過(guò)四六級(jí),六級(jí)521,求調(diào)劑 +4 | 美味可樂(lè)雞翅 2026-03-26 | 4/200 |
|
|
[考研] 化學(xué)308分求調(diào)劑 +8 | 你好明天你好 2026-03-23 | 9/450 |
|
|
[考研] 一志愿吉大071010,316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
[考研] 325求調(diào)劑 +5 | 李嘉圖·S·路 2026-03-23 | 5/250 |
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑!。! +4 | Acaciad 2026-03-25 | 4/200 |
|
|
[考研]
|
黃粱一夢(mèng)千年 2026-03-24 | 3/150 |
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|