 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 1705 | 回復(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)計算這類積分的？ |
回復(fù)此樓» 猜你喜歡
大連海事大學(xué)船舶潔凈能源研究中心2026年博士研究生招生啟事
已經(jīng)有10人回復(fù)
-
2026年湘潭大學(xué)接收化學(xué)工程與技術(shù)相關(guān)方向?qū)W術(shù)型碩士研究生調(diào)劑生
已經(jīng)有1人回復(fù)
-
物理化學(xué)論文潤色/翻譯怎么收費?
已經(jīng)有102人回復(fù)
-
2026湖北師范大學(xué)化學(xué)化工學(xué)院招收碩士研究生調(diào)劑
已經(jīng)有0人回復(fù)
-
2026湖北師范大學(xué)化學(xué)化工學(xué)院招收研究生調(diào)劑
已經(jīng)有0人回復(fù)
-
26年材料/化學(xué)相關(guān)專業(yè)博士研究生招生 溫州大學(xué)化學(xué)與材料工程學(xué)院(26年秋季入學(xué))
已經(jīng)有0人回復(fù)
-
0703調(diào)劑
已經(jīng)有2人回復(fù)
-
2026年博士研究生招生
已經(jīng)有1人回復(fù)
-
304求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求大蝦指導(dǎo) 一個關(guān)于atoi函數(shù)的問題
已經(jīng)有3人回復(fù)
- 利用概率密度函數(shù)生成隨機變量過程中出現(xiàn)超幾何分布函數(shù)問題
已經(jīng)有4人回復(fù)
- 程序員必讀書單【轉(zhuǎn)載】
已經(jīng)有22人回復(fù)
- 博士點基金課題申請人和主要參加人員論文或論著情況表如何填寫
已經(jīng)有4人回復(fù)
- 雅思寫作流程圖【轉(zhuǎn)載】
已經(jīng)有55人回復(fù)
- 請高手指點一下怎么確定一個符號矩陣中的等值元素
已經(jīng)有6人回復(fù)
- c++編程基本思想
已經(jīng)有176人回復(fù)
- 縮短程序運行時間的方法
已經(jīng)有7人回復(fù)
- 【求助】存款證明上寫多長時間？(具體解答已經(jīng)編輯更新，望對大家有幫助)
已經(jīng)有42人回復(fù)
- Common Technical Writing Issues, 一個關(guān)于寫英文文章細節(jié)的很好的PPT
已經(jīng)有25人回復(fù)
fangyongxinxi
新蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 41
- 散金: 69
- 紅花: 3
- 帖子: 876
- 在線: 138.5小時
- 蟲號: 437418
- 注冊: 2007-08-26
- 性別: MM
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
2樓2010-09-01 09:03:07
sobereva
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
- QC強帖: 17
- 應(yīng)助: 226 (大學(xué)生)
- 貴賓: 5.4
- 金幣: 2180.5
- 紅花: 217
- 帖子: 1241
- 在線: 366.8小時
- 蟲號: 281271
- 注冊: 2006-09-23
- 專業(yè): 宇宙學(xué)
★ ★ ★
erylingjet(金幣+3):非常感謝！ 2010-09-01 14:14:36
fangyongxinxi(金幣+1): 2010-09-01 22:19:29
erylingjet(金幣+3):非常感謝！ 2010-09-01 14:14:36
fangyongxinxi(金幣+1): 2010-09-01 22:19:29
|
這要用高斯-厄米積分，高斯積分在量化計算中用處很大。 重疊積分在《量子化學(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 ] |
3樓2010-09-01 13:24:09
fkingdeath
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1.2
- 帖子: 22
- 在線: 16.6小時
- 蟲號: 1078366
- 注冊: 2010-08-20
- 性別: GG
- 專業(yè): 固體與表面物理化學(xué)
4樓2010-09-01 14:08:58
fangyongxinxi
新蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 41
- 散金: 69
- 紅花: 3
- 帖子: 876
- 在線: 138.5小時
- 蟲號: 437418
- 注冊: 2007-08-26
- 性別: MM
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
5樓2010-09-01 22:20:03
songyanyan
新蟲 (正式寫手)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1264.4
- 散金: 24
- 帖子: 337
- 在線: 166.7小時
- 蟲號: 929458
- 注冊: 2009-12-16
- 性別: MM
- 專業(yè): 催化化學(xué)
6樓2010-09-04 09:31:04
coolrainbow
木蟲 (著名寫手)
未來國家凍涼
- QC強帖: 11
- 應(yīng)助: 1 (幼兒園)
- 貴賓: 0.1
- 金幣: 3979.5
- 散金: 10
- 紅花: 50
- 帖子: 1137
- 在線: 272.3小時
- 蟲號: 239560
- 注冊: 2006-04-07
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
7樓2010-09-04 11:32:26
sandf
木蟲 (正式寫手)
- 應(yīng)助: 15 (小學(xué)生)
- 金幣: 1631.4
- 散金: 121
- 紅花: 2
- 帖子: 604
- 在線: 136.6小時
- 蟲號: 305516
- 注冊: 2006-12-10
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
8樓2010-09-05 08:40:30
maomao1210
金蟲 (正式寫手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 1432.8
- 散金: 242
- 紅花: 16
- 沙發(fā): 1
- 帖子: 991
- 在線: 442.1小時
- 蟲號: 253215
- 注冊: 2006-05-20
- 性別: MM
- 專業(yè): 考古理論
9樓2010-09-05 10:12:34
fangyongxinxi
新蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 41
- 散金: 69
- 紅花: 3
- 帖子: 876
- 在線: 138.5小時
- 蟲號: 437418
- 注冊: 2007-08-26
- 性別: MM
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
10樓2010-09-23 13:09:09
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 fangyongxinxi 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 311（085601）求調(diào)劑 +3 | liziyeyeye 2026-03-28 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
[考研] 0703一志愿9，初試成績：338，四六級已過，有科研經(jīng)歷，求調(diào)劑！ +4 | Zuhui0306 2026-03-25 | 4/200 |
|
|
|
[考研] 085602 307分 求調(diào)劑 +7 | 不知道叫什么！ 2026-03-26 | 7/350 |
|
|
|
[考研] 材料求調(diào)劑一志愿哈工大324 +7 | 閆旭東 2026-03-28 | 9/450 |
|
|
|
[考研] 331環(huán)境科學(xué)與工程求調(diào)劑 +3 | 熠然好運氣 2026-03-27 | 3/150 |
|
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
|
[考研] 材料與化工085600，總分304，本科有兩篇sci參與，求調(diào)劑 +10 | 幸運的醬醬 2026-03-22 | 12/600 |
|
|
|
[考研] 材料292調(diào)劑 +12 | 橘頌思美人 2026-03-23 | 12/600 |
|
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
|
[考研] 085600，材料與化工321分，求調(diào)劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學(xué)習(xí) 2026-03-23 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +6 | 林之夕 2026-03-24 | 6/300 |
|
|
|
[考研] 325求調(diào)劑 +5 | 李嘉圖·S·路 2026-03-23 | 5/250 |
|
|
|
[考研] 081700 調(diào)劑 267分 +11 | 迷人的哈哈 2026-03-23 | 11/550 |
|
|
|
[考研] 材料與化工328分調(diào)劑 +6 | 。，。，。，。i 2026-03-23 | 6/300 |
|
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
|
[考研] 300分，材料，求調(diào)劑，英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|