| 查看: 2754 | 回復: 10 | ||||
[求助]
gaussian計算頻率時輸出的hessian是那一部分啊
|
||||
|
我用水做一例子計算hessian矩陣。以下是輸出文件中的一部分,但是不太清楚哪部分是。是最后的 Hessian after L703:后面的嗎。是的話,水的話其hessian矩陣是9*9的嗎,那么它上三角的內容是什么呢。還有下面出現(xiàn)的D是什么意思啊。是我們平常說的指數(shù)e嗎? Hessian from J: 1 2 3 4 5 1 0.000000D+00 2 0.000000D+00 0.000000D+00 3 0.000000D+00 0.000000D+00 0.000000D+00 4 0.000000D+00 0.000000D+00 0.000000D+00 -0.441362D+01 5 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00 -0.948350D+00 6 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00 7 0.000000D+00 0.000000D+00 0.000000D+00 0.220681D+01 0.000000D+00 8 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00 0.474175D+00 9 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00 0.121877D+01 10 0.000000D+00 0.000000D+00 0.000000D+00 0.220681D+01 0.000000D+00 11 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00 0.474175D+00 12 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00 -0.121877D+01 6 7 8 9 10 6 -0.273959D+01 7 0.000000D+00 -0.214794D+01 8 0.121674D+01 0.000000D+00 -0.602963D+00 9 0.136980D+01 0.000000D+00 -0.121776D+01 -0.129566D+01 10 0.000000D+00 -0.588633D-01 0.000000D+00 0.000000D+00 -0.214794D+01 11 -0.121674D+01 0.000000D+00 0.128788D+00 -0.101712D-02 0.000000D+00 12 0.136980D+01 0.000000D+00 0.101712D-02 -0.741375D-01 0.000000D+00 11 12 11 -0.602963D+00 12 0.121776D+01 -0.129566D+01 Force constants from FoF: 1 2 3 4 5 1 0.000000D+00 2 0.000000D+00 0.000000D+00 3 0.000000D+00 0.000000D+00 0.000000D+00 4 0.000000D+00 0.000000D+00 0.000000D+00 -0.441362D+01 5 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00 -0.948350D+00 6 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00 7 0.000000D+00 0.000000D+00 0.000000D+00 0.220681D+01 0.000000D+00 8 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00 0.474175D+00 9 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00 0.121877D+01 10 0.000000D+00 0.000000D+00 0.000000D+00 0.220681D+01 0.000000D+00 11 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00 0.474175D+00 12 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00 -0.121877D+01 6 7 8 9 10 6 -0.273959D+01 7 0.000000D+00 -0.214794D+01 8 0.121674D+01 0.000000D+00 -0.602963D+00 9 0.136980D+01 0.000000D+00 -0.121776D+01 -0.129566D+01 10 0.000000D+00 -0.588633D-01 0.000000D+00 0.000000D+00 -0.214794D+01 11 -0.121674D+01 0.000000D+00 0.128788D+00 -0.101712D-02 0.000000D+00 12 0.136980D+01 0.000000D+00 0.101712D-02 -0.741375D-01 0.000000D+00 11 12 11 -0.602963D+00 12 0.121776D+01 -0.129566D+01 Second derivatives from L703 IAtSt=0: 1 2 3 4 5 1 0.000000D+00 2 0.000000D+00 0.000000D+00 3 0.000000D+00 0.000000D+00 0.000000D+00 4 0.000000D+00 0.000000D+00 0.000000D+00 -0.441362D+01 5 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00 -0.948350D+00 6 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00 7 0.000000D+00 0.000000D+00 0.000000D+00 0.220681D+01 0.000000D+00 8 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00 0.474175D+00 9 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00 0.121877D+01 10 0.000000D+00 0.000000D+00 0.000000D+00 0.220681D+01 0.000000D+00 11 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00 0.474175D+00 12 0.000000D+00 0.000000D+00 0.000000D+00 0.000000D+00 -0.121877D+01 6 7 8 9 10 6 -0.273959D+01 7 0.000000D+00 -0.214794D+01 8 0.121674D+01 0.000000D+00 -0.602963D+00 9 0.136980D+01 0.000000D+00 -0.121776D+01 -0.129566D+01 10 0.000000D+00 -0.588633D-01 0.000000D+00 0.000000D+00 -0.214794D+01 11 -0.121674D+01 0.000000D+00 0.128788D+00 -0.101712D-02 0.000000D+00 12 0.136980D+01 0.000000D+00 0.101712D-02 -0.741375D-01 0.000000D+00 11 12 11 -0.602963D+00 12 0.121776D+01 -0.129566D+01 Polarizability into L703: 1 2 3 1 0.783222D+00 2 0.000000D+00 0.680812D+01 3 0.000000D+00 0.000000D+00 0.362564D+01 Polarizability after L703: 1 2 3 1 0.783222D+00 2 0.000000D+00 0.680812D+01 3 0.000000D+00 0.000000D+00 0.362564D+01 Dipole Derivatives into L703: 1 2 3 4 5 1 -0.871467D+00 0.000000D+00 0.000000D+00 0.435734D+00 0.000000D+00 2 0.000000D+00 -0.387147D+00 0.000000D+00 0.000000D+00 0.193574D+00 3 0.000000D+00 0.000000D+00 -0.329661D+00 0.000000D+00 0.103814D+00 6 7 8 9 1 0.000000D+00 0.435734D+00 0.000000D+00 0.000000D+00 2 0.176832D+00 0.000000D+00 0.193574D+00 -0.176832D+00 3 0.164831D+00 0.000000D+00 -0.103814D+00 0.164831D+00 Dipole Derivatives after L703: 1 2 3 4 5 1 -0.871467D+00 0.000000D+00 0.000000D+00 0.435734D+00 0.000000D+00 2 0.000000D+00 -0.387147D+00 0.000000D+00 0.000000D+00 0.193574D+00 3 0.000000D+00 0.000000D+00 -0.329661D+00 0.000000D+00 0.103814D+00 6 7 8 9 1 0.000000D+00 0.435734D+00 0.000000D+00 0.000000D+00 2 0.176832D+00 0.000000D+00 0.193574D+00 -0.176832D+00 3 0.164831D+00 0.000000D+00 -0.103814D+00 0.164831D+00 Hessian into L703: 1 2 3 4 5 1 0.441373D+01 2 0.000000D+00 0.166319D+01 3 0.000000D+00 0.000000D+00 0.317440D+01 4 -0.220687D+01 0.000000D+00 0.000000D+00 0.214806D+01 5 0.000000D+00 -0.831595D+00 -0.103436D+01 0.000000D+00 0.997733D+00 6 0.000000D+00 -0.957818D+00 -0.158720D+01 0.000000D+00 0.996090D+00 7 -0.220687D+01 0.000000D+00 0.000000D+00 0.588069D-01 0.000000D+00 8 0.000000D+00 -0.831595D+00 0.103436D+01 0.000000D+00 -0.166139D+00 9 0.000000D+00 0.957818D+00 -0.158720D+01 0.000000D+00 0.382723D-01 6 7 8 9 6 0.149973D+01 7 0.000000D+00 0.214806D+01 8 -0.382723D-01 0.000000D+00 0.997733D+00 9 0.874759D-01 0.000000D+00 -0.996090D+00 0.149973D+01 Hessian after L703: 1 2 3 4 5 1 0.117370D-03 2 0.000000D+00 0.714839D+00 3 0.000000D+00 0.000000D+00 0.434810D+00 4 -0.586851D-04 0.000000D+00 0.000000D+00 0.115086D-03 5 0.000000D+00 -0.357419D+00 0.182376D+00 0.000000D+00 0.394771D+00 6 0.000000D+00 0.260955D+00 -0.217405D+00 0.000000D+00 -0.221665D+00 7 -0.586851D-04 0.000000D+00 0.000000D+00 -0.564006D-04 0.000000D+00 8 0.000000D+00 -0.357419D+00 -0.182376D+00 0.000000D+00 -0.373513D-01 9 0.000000D+00 -0.260955D+00 -0.217405D+00 0.000000D+00 0.392894D-01 6 7 8 9 6 0.204067D+00 7 0.000000D+00 0.115086D-03 8 -0.392894D-01 0.000000D+00 0.394771D+00 9 0.133384D-01 0.000000D+00 0.221665D+00 0.204067D+00 |
自己感覺好的東西 |
專家顧問 (職業(yè)作家)
地溝油冶煉專家
![]() |
專家經驗: +458 |
|
M^1/2 * hessian * M^1/2就是質量加權的hessian。M是原子量矩陣。Matlab語言的寫法是 sqrtm(M) * hessian * sqrtm(M)。 有的對角化程序不需要提前計算質量加權的hessian。Matlab命令: [p,q]=eig(hessian,M); q=diag(q); p是質量加權的簡正模式(不是Gaussian程序打印的那種),q是原子單位的特征值。 5140.48715246*sqrt(q)是波數(shù)單位的振動頻率。 |

專家顧問 (職業(yè)作家)
地溝油冶煉專家
![]() |
專家經驗: +458 |

木蟲 (著名寫手)
木蟲 (著名寫手)
專家顧問 (職業(yè)作家)
地溝油冶煉專家
![]() |
專家經驗: +458 |

|
感謝您的回復!麻煩再問一下您公式中的(hessian,M)是社么意思,質量加權的hessian矩陣嗎,是Hessian after L703后面的那個矩陣嗎。那么Hessian after L703后面的那個矩陣跟Force constants in Cartesian coordinates后面的那個矩陣有什么不同。因為本身計算頻率的hessian矩陣就是力常數(shù)矩陣啊。還有,gaussian輸出文件里的力常數(shù)單位是mdyn*A-1*amu-1嗎,轉換成波數(shù)的單位是利用公式sqrt(q)/2*pi*c吧,我怎么得不到5140.48715246這個值呢。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調劑一志愿武漢理工大學材料工程(085601) +3 | WW.' 2026-03-23 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 081700 調劑 267分 +5 | 迷人的哈哈 2026-03-23 | 5/250 |
|
|
[考研] 306求0703調劑一志愿華中師范 +6 | 紙魚ly 2026-03-21 | 6/300 |
|
|
[考研] 085600材料與化工 +10 | 安全上岸! 2026-03-16 | 10/500 |
|
|
[考研] 一志愿西安交通大學材料工程專業(yè) 282分求調劑 +11 | 楓橋ZL 2026-03-18 | 13/650 |
|
|
[考研] 275求調劑 +6 | shansx 2026-03-22 | 8/400 |
|
|
[考研] 生物學一志愿985,分數(shù)349求調劑 +4 | zxts12 2026-03-21 | 7/350 |
|
|
[考研] 材料學碩301分求調劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
[考研]
|
Grand777 2026-03-21 | 3/150 |
|
|
[考研] 0703化學調劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
[考研] 266求調劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
[考研] 085700資源與環(huán)境308求調劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
[考研] 304求調劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
[考研] 材料專業(yè)求調劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
[考研] 274求調劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
[考研] 一志愿 西北大學 ,070300化學學碩,總分287,雙非一本,求調劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
[考研]
|
然11 2026-03-19 | 4/200 |
|
|
[考研] 求調劑一志愿南京航空航天大學289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] 288求調劑,一志愿華南理工大學071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|