| 5 | 1/1 | 返回列表 |
| 查看: 2770 | 回復: 10 | ||||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||
[求助]
gaussian計算頻率時輸出的hessian是那一部分啊
|
||||
|
我用水做一例子計算hessian矩陣。以下是輸出文件中的一部分,但是不太清楚哪部分是。是最后的 Hessian after L703:后面的嗎。是的話,水的話其hessian矩陣是9*9的嗎,那么它上三角的內(nèi)容是什么呢。還有下面出現(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è)作家)
地溝油冶煉專家
![]() |
專家經(jīng)驗: +458 |
|
M^1/2 * hessian * M^1/2就是質(zhì)量加權(quán)的hessian。M是原子量矩陣。Matlab語言的寫法是 sqrtm(M) * hessian * sqrtm(M)。 有的對角化程序不需要提前計算質(zhì)量加權(quán)的hessian。Matlab命令: [p,q]=eig(hessian,M); q=diag(q); p是質(zhì)量加權(quán)的簡正模式(不是Gaussian程序打印的那種),q是原子單位的特征值。 5140.48715246*sqrt(q)是波數(shù)單位的振動頻率。 |

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

木蟲 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 298求B區(qū)調(diào)劑 +4 | zzz,,r 2026-04-02 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 292分,材料與化工,申請調(diào)劑 +19 | 程晴之 2026-04-01 | 23/1150 |
|
|
[考研] 材料調(diào)劑 +3 | 一樣YWY 2026-04-02 | 3/150 |
|
|
[考研] 366求調(diào)劑一志愿東北大學 +6 | 運氣來得若有似?/a> 2026-04-02 | 6/300 |
|
|
[考研] 07生物學求調(diào)劑 一志愿同濟大學359分 +3 | LAMC. 2026-03-30 | 3/150 |
|
|
[考研] 085600,材料與化工321分調(diào)劑 +8 | 大饞小子 2026-03-27 | 10/500 |
|
|
[考研] 266求調(diào)劑 +4 | 學員97LZgn 2026-04-02 | 4/200 |
|
|
[考研] 生物學327,求調(diào)劑 +5 | 書上的梅子 2026-04-01 | 6/300 |
|
|
[考研] 298求調(diào)劑 +4 | 什么是胖頭魚 2026-03-30 | 6/300 |
|
|
[考研] 070300化學279求調(diào)劑 +15 | 哈哈哈^_^ 2026-03-31 | 17/850 |
|
|
[考研] 0855機械初試280求調(diào)劑 +3 | kazenotori 2026-03-31 | 3/150 |
|
|
[考研] 考研材料工程351分調(diào)劑 +5 | 整個好的 2026-03-31 | 5/250 |
|
|
[考研] 一志愿:西北大學,英一數(shù)一408-284分求調(diào)劑 +7 | 12.27 2026-03-27 | 7/350 |
|
|
[考研] 0856求調(diào)劑 +9 | 楒桉 2026-03-28 | 9/450 |
|
|
[考研] 254材料與化工求調(diào)劑 +3 | 翰冬林楠 2026-03-30 | 4/200 |
|
|
[考研] 343求調(diào)劑 +8 | 愛羈絆 2026-03-28 | 8/400 |
|
|
[考研] 313求調(diào)劑 +6 | 賣個關子吧 2026-03-31 | 6/300 |
|
|
[考研] 求調(diào)劑 +10 | 家佳佳佳佳佳 2026-03-29 | 10/500 |
|
|
[考研] 312,生物學求調(diào)劑 +3 | 小譯同學abc 2026-03-28 | 3/150 |
|
|
[考研] 一志愿吉大071010,316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|