| 查看: 1708 | 回復: 10 | |||
fangyongxinxi新蟲 (正式寫手)
|
[交流]
【求助】Q2: 如何自寫程序實現(xiàn)GTO 1s 型重疊矩陣 已有2人參與
|
|
在《量子化學--基本原理和從頭算法》,240頁,講到: 分子自洽場計算 氨分子 其中,N1s 型的STO-3G表達式為: y = C11*(exp(-A11*x*x))+C12*(exp(-A12*x*x))+C13*(exp(-A13*x*x)) 其中, C11 = 0.1543*22.5377; A11=99.9997 C12 = 0.5353*6.284; A12=18.2151 C13 = 0.4446*2.3579; A13=4.9297 問題是:如何自寫程序---積分程序,來實現(xiàn)重疊積分, 即如何實現(xiàn): 自己嘗試了Python語言中的integrate 和numpy.trapz函數(shù) 得到的結果都不正確,郁悶中。 如果用乘積定理來實現(xiàn)的話,倒是可以,但指數(shù)函數(shù)的無窮邊界如何積分,不清楚。 大家是怎么實現(xiàn)計算這類積分的? |
新蟲 (正式寫手)
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
|
這要用高斯-厄米積分,高斯積分在量化計算中用處很大。 重疊積分在《量子化學中的計算方法》的4.1.4節(jié)推導很簡明清楚,很容易就能寫出代碼來。但注意4.1-28式是錯的,√ζ應改成1/√ζ GTO之間的重疊積分會寫了,寫GTO收縮的基函數(shù)之間的重疊積分就十分容易,套上循環(huán)做累加就可以了。 [ Last edited by sobereva on 2010-9-1 at 13:28 ] |
銅蟲 (初入文壇)
新蟲 (正式寫手)
新蟲 (正式寫手)
木蟲 (著名寫手)
未來國家凍涼

木蟲 (正式寫手)

金蟲 (正式寫手)
新蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 289求調劑 +12 | 新時代材料 2026-03-27 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 316求調劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
[考研] 343求調劑 +5 | 愛羈絆 2026-03-28 | 5/250 |
|
|
[考研] 0856材料化工調劑 總分330 +11 | zhubinhao 2026-03-27 | 11/550 |
|
|
[考研] 一志愿華理,數(shù)一英一285求A區(qū)調劑 +8 | AZMK 2026-03-25 | 12/600 |
|
|
[考研] 0703化學求調劑 +9 | 奶油草莓. 2026-03-22 | 10/500 |
|
|
[考研]
|
孅華 2026-03-22 | 17/850 |
|
|
[考研] 340求調劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
[考研] 331環(huán)境科學與工程求調劑 +3 | 熠然好運氣 2026-03-27 | 3/150 |
|
|
[考研] 291求調劑 +7 | 孅華 2026-03-22 | 7/350 |
|
|
[考研] 272求調劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 一志愿南師大0703化學 275求調劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
[考研] 化學308分求調劑 +8 | 你好明天你好 2026-03-23 | 9/450 |
|
|
[考研] 0703化學338求調劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 349求調劑 +4 | 李木子啊哈哈 2026-03-25 | 4/200 |
|
|
[考研] 329求調劑 +5 | 1() 2026-03-22 | 5/250 |
|
|
[考研] 化學調劑一志愿上海交通大學336分-本科上海211 +4 | 小魚愛有機 2026-03-25 | 4/200 |
|
|
[考研] 344求調劑 +3 | desto 2026-03-24 | 3/150 |
|
|
[考研] 315分,誠求調劑,材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|
|
[考研] 275求調劑 +6 | shansx 2026-03-22 | 8/400 |
|