| 5 | 1/1 | 返回列表 |
| 查看: 1710 | 回復(fù): 10 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
fangyongxinxi新蟲(chóng) (正式寫(xiě)手)
|
[交流]
【求助】Q2: 如何自寫(xiě)程序?qū)崿F(xiàn)GTO 1s 型重疊矩陣 已有2人參與
|
||
|
在《量子化學(xué)--基本原理和從頭算法》,240頁(yè),講到: 分子自洽場(chǎng)計(jì)算 氨分子 其中,N1s 型的STO-3G表達(dá)式為: 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 問(wèn)題是:如何自寫(xiě)程序---積分程序,來(lái)實(shí)現(xiàn)重疊積分, 即如何實(shí)現(xiàn): 自己嘗試了Python語(yǔ)言中的integrate 和numpy.trapz函數(shù) 得到的結(jié)果都不正確,郁悶中。 如果用乘積定理來(lái)實(shí)現(xiàn)的話,倒是可以,但指數(shù)函數(shù)的無(wú)窮邊界如何積分,不清楚。 大家是怎么實(shí)現(xiàn)計(jì)算這類積分的? |
木蟲(chóng) (著名寫(xiě)手)
未來(lái)國(guó)家凍涼

新蟲(chóng) (正式寫(xiě)手)
至尊木蟲(chóng) (著名寫(xiě)手)
本人已永久離開(kāi)小木蟲(chóng)
|
這要用高斯-厄米積分,高斯積分在量化計(jì)算中用處很大。 重疊積分在《量子化學(xué)中的計(jì)算方法》的4.1.4節(jié)推導(dǎo)很簡(jiǎn)明清楚,很容易就能寫(xiě)出代碼來(lái)。但注意4.1-28式是錯(cuò)的,√ζ應(yīng)改成1/√ζ GTO之間的重疊積分會(huì)寫(xiě)了,寫(xiě)GTO收縮的基函數(shù)之間的重疊積分就十分容易,套上循環(huán)做累加就可以了。 [ Last edited by sobereva on 2010-9-1 at 13:28 ] |
銅蟲(chóng) (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 327求調(diào)劑 +4 | 汲亦昊 2026-03-29 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 312,生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
[考研] 求化學(xué)調(diào)劑 +4 | wulanna 2026-03-28 | 4/200 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +9 | 奶油草莓. 2026-03-22 | 10/500 |
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
[考研] 266分,求材料冶金能源化工等調(diào)劑 +7 | 哇呼哼呼哼 2026-03-27 | 9/450 |
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
[考研] 一志愿陜師大生物學(xué)071000,298分,求調(diào)劑 +5 | SYA! 2026-03-23 | 5/250 |
|
|
[考研] 324求調(diào)劑 +8 | hanamiko 2026-03-26 | 10/500 |
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
[考研] 341求調(diào)劑 +7 | 青檸檬1 2026-03-26 | 7/350 |
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
[考研] 考研一志愿蘇州大學(xué)初始315(英一)求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|
|
[考博] 26申博自薦 +3 | whh869393 2026-03-24 | 3/150 |
|