| 5 | 1/1 | 返回列表 |
| 查看: 1709 | 回復(fù): 10 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
fangyongxinxi新蟲 (正式寫手)
|
[交流]
【求助】Q2: 如何自寫程序?qū)崿F(xiàn)GTO 1s 型重疊矩陣 已有2人參與
|
||
|
在《量子化學(xué)--基本原理和從頭算法》,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ù) 得到的結(jié)果都不正確,郁悶中。 如果用乘積定理來實現(xiàn)的話,倒是可以,但指數(shù)函數(shù)的無窮邊界如何積分,不清楚。 大家是怎么實現(xiàn)計算這類積分的? |
金蟲 (正式寫手)
新蟲 (正式寫手)
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
|
這要用高斯-厄米積分,高斯積分在量化計算中用處很大。 重疊積分在《量子化學(xué)中的計算方法》的4.1.4節(jié)推導(dǎo)很簡明清楚,很容易就能寫出代碼來。但注意4.1-28式是錯的,√ζ應(yīng)改成1/√ζ GTO之間的重疊積分會寫了,寫GTO收縮的基函數(shù)之間的重疊積分就十分容易,套上循環(huán)做累加就可以了。 [ Last edited by sobereva on 2010-9-1 at 13:28 ] |
銅蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311(085601)求調(diào)劑 +4 | liziyeyeye 2026-03-28 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 312,生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
[考研] 085600,材料與化工321分求調(diào)劑 +9 | 大饞小子 2026-03-28 | 9/450 |
|
|
[考研]
|
y7czhao 2026-03-26 | 10/500 |
|
|
[考研] 347求調(diào)劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
[考研] 調(diào)劑 +3 | 好好讀書。 2026-03-28 | 3/150 |
|
|
[考研] 藥學(xué)105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +9 | vv迷 2026-03-22 | 9/450 |
|
|
[考研] 復(fù)試調(diào)劑,一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
[考研] 材料求調(diào)劑 +5 | .m.. 2026-03-25 | 5/250 |
|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +4 | 柳溪邊 2026-03-26 | 4/200 |
|
|
[考研] 281求調(diào)劑 +3 | 亞克西good 2026-03-26 | 5/250 |
|
|
[考研] 中國科學(xué)院深圳先進技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-26 | 6/300 |
|
|
[考研] 總分293求調(diào)劑 +6 | 加一一九 2026-03-25 | 8/400 |
|
|
[考研] 考研一志愿蘇州大學(xué)初始315(英一)求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
[考研] 一志愿南航材料專317分求調(diào)劑 +5 | 炸呀炸呀炸薯條 2026-03-23 | 5/250 |
|