| 5 | 1/1 | 返回列表 |
| 查看: 1706 | 回復(fù): 10 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
fangyongxinxi新蟲 (正式寫手)
|
[交流]
【求助】Q2: 如何自寫程序?qū)崿F(xiàn)GTO 1s 型重疊矩陣 已有2人參與
|
||
|
在《量子化學(xué)--基本原理和從頭算法》,240頁,講到: 分子自洽場計算 氨分子 其中,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 問題是:如何自寫程序---積分程序,來實現(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ā)表 | |
|---|---|---|---|---|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 300求調(diào)劑,材料科學(xué)英一數(shù)二 +6 | leaflight 2026-03-24 | 6/300 |
|
|
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +3 | nhj_ 2026-03-25 | 3/150 |
|
|
[考研] 291求調(diào)劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
[考研]
|
孅華 2026-03-22 | 17/850 |
|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +9 | 丹青奶蓋 2026-03-26 | 10/500 |
|
|
[考研] 330一志愿中國海洋大學(xué) 化學(xué)工程 085602 有讀博意愿 求調(diào)劑 +3 | wywy.. 2026-03-27 | 4/200 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 283求調(diào)劑(080500) +4 | A child 2026-03-27 | 4/200 |
|
|
[考研]
材料調(diào)劑
5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
[考研] 281求調(diào)劑 +3 | 亞克西good 2026-03-26 | 5/250 |
|
|
[考研] 340求調(diào)劑 +3 | Amber00 2026-03-26 | 3/150 |
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +4 | 李李不服輸 2026-03-26 | 4/200 |
|
|
[考研] 086000生物與醫(yī)藥292求調(diào)劑 +6 | 小小陳小小 2026-03-22 | 9/450 |
|
|
[考研] 26考研-291分-廈門大學(xué)(085601)-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +3 | min3 2026-03-24 | 4/200 |
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
[考研] 一志愿南航材料專317分求調(diào)劑 +5 | 炸呀炸呀炸薯條 2026-03-23 | 5/250 |
|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
[考研] 306求調(diào)劑 +5 | 來好運來來來 2026-03-22 | 5/250 |
|