 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1707 | 回復: 10 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
fangyongxinxi新蟲 (正式寫手)
|
[交流]
【求助】Q2: 如何自寫程序?qū)崿F(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ù) 得到的結(jié)果都不正確，郁悶中。 如果用乘積定理來實現(xiàn)的話，倒是可以，但指數(shù)函數(shù)的無窮邊界如何積分，不清楚。 大家是怎么實現(xiàn)計算這類積分的？ |
回復此樓» 猜你喜歡
細胞培養(yǎng)，這22個細節(jié)一定要注意�。ㄒ唬�
已經(jīng)有0人回復
-
CSC訪學博后項目獲批，外方學校暫停合作該怎么辦？
已經(jīng)有51人回復
-
物理化學論文潤色/翻譯怎么收費?
已經(jīng)有83人回復
-
大灣區(qū)大學劉天輝課題組招收2026級博士生（與中山大學物理學院聯(lián)合培養(yǎng)）
已經(jīng)有22人回復
-
PbS量子點紫外吸收
已經(jīng)有0人回復
-
比利時魯汶大學與國家留學基金委共同資助博士研究生CSC-KU Leuven PhD Scholarship
已經(jīng)有0人回復
-
七嗪類物質(zhì)合成求助
已經(jīng)有0人回復
-
倫敦大學學院Benjamin Moss博士招收人工智能、原位光譜、催化相關(guān)26年博士生
已經(jīng)有11人回復
-
分子動力學模擬合作需求
已經(jīng)有2人回復
-
中國科學院成都山地災害與環(huán)境研究所2026年特別研究助理（博士后）招聘啟事
已經(jīng)有0人回復
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求大蝦指導 一個關(guān)于atoi函數(shù)的問題
已經(jīng)有3人回復
- 利用概率密度函數(shù)生成隨機變量過程中出現(xiàn)超幾何分布函數(shù)問題
已經(jīng)有4人回復
- 程序員必讀書單【轉(zhuǎn)載】
已經(jīng)有22人回復
- 博士點基金課題申請人和主要參加人員論文或論著情況表如何填寫
已經(jīng)有4人回復
- 雅思寫作流程圖【轉(zhuǎn)載】
已經(jīng)有55人回復
- 請高手指點一下怎么確定一個符號矩陣中的等值元素
已經(jīng)有6人回復
- c++編程基本思想
已經(jīng)有176人回復
- 縮短程序運行時間的方法
已經(jīng)有7人回復
- 【求助】存款證明上寫多長時間？(具體解答已經(jīng)編輯更新，望對大家有幫助)
已經(jīng)有42人回復
- Common Technical Writing Issues, 一個關(guān)于寫英文文章細節(jié)的很好的PPT
已經(jīng)有25人回復
fkingdeath
銅蟲 (初入文壇)
- 應助: 0 (幼兒園)
- 金幣: 1.2
- 帖子: 22
- 在線: 16.6小時
- 蟲號: 1078366
- 注冊: 2010-08-20
- 性別: GG
- 專業(yè): 固體與表面物理化學
4樓2010-09-01 14:08:58
fangyongxinxi
新蟲 (正式寫手)
- 應助: 0 (幼兒園)
- 金幣: 41
- 散金: 69
- 紅花: 3
- 帖子: 876
- 在線: 138.5小時
- 蟲號: 437418
- 注冊: 2007-08-26
- 性別: MM
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學和
2樓2010-09-01 09:03:07
sobereva
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
- QC強帖: 17
- 應助: 226 (大學生)
- 貴賓: 5.4
- 金幣: 2180.5
- 紅花: 217
- 帖子: 1241
- 在線: 366.8小時
- 蟲號: 281271
- 注冊: 2006-09-23
- 專業(yè): 宇宙學
★ ★ ★
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
|
這要用高斯-厄米積分，高斯積分在量化計算中用處很大。 重疊積分在《量子化學中的計算方法》的4.1.4節(jié)推導很簡明清楚，很容易就能寫出代碼來。但注意4.1-28式是錯的，√ζ應改成1/√ζ GTO之間的重疊積分會寫了，寫GTO收縮的基函數(shù)之間的重疊積分就十分容易，套上循環(huán)做累加就可以了。 [ Last edited by sobereva on 2010-9-1 at 13:28 ] |
3樓2010-09-01 13:24:09
fangyongxinxi
新蟲 (正式寫手)
- 應助: 0 (幼兒園)
- 金幣: 41
- 散金: 69
- 紅花: 3
- 帖子: 876
- 在線: 138.5小時
- 蟲號: 437418
- 注冊: 2007-08-26
- 性別: MM
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學和
5樓2010-09-01 22:20:03
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 343求調(diào)劑 +5 | 愛羈絆 2026-03-28 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 321求調(diào)劑 +6 | 璞玉~~ 2026-03-25 | 7/350 |
|
|
[考研] 數(shù)一英一271專碩（085401）求調(diào)劑，可跨 +4 | 前行必有光 2026-03-28 | 5/250 |
|
|
|
[考研] 277跪求調(diào)劑 +5 | 1915668 2026-03-27 | 9/450 |
|
|
|
[考研] 一志愿南師大0703化學 275求調(diào)劑 +4 | Ripcord上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 287求調(diào)劑 +10 | land xuxu 2026-03-26 | 10/500 |
|
|
|
[考研] 一志愿鄭大085600，310分求調(diào)劑 +5 | 李瀟可 2026-03-26 | 5/250 |
|
|
|
[碩博家園] 招收生物學/細胞生物學調(diào)劑 +3 | IceGuo 2026-03-26 | 4/200 |
|
|
|
[考研]
材料調(diào)劑
 5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
|
[考研] 340求調(diào)劑 +3 | Amber00 2026-03-26 | 3/150 |
|
|
|
[考研]
 085602 289分求調(diào)劑
+8
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
|
[考研] 086000生物與醫(yī)藥292求調(diào)劑 +6 | 小小陳小小 2026-03-22 | 9/450 |
|
|
|
[考研] 081700 調(diào)劑 267分 +11 | 迷人的哈哈 2026-03-23 | 11/550 |
|
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑 +6 | jiajunX 2026-03-22 | 6/300 |
|
|
|
[考研] 296求調(diào)劑 +4 | 汪��？！ 2026-03-25 | 7/350 |
|
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|
|
|
[考研] 0854電子信息求調(diào)劑 324 +4 | Promise-jyl 2026-03-23 | 4/200 |
|
|
|
[考研] 一志愿河北工業(yè)大學0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 12/600 |
|
|
|
[考研] 315分，誠求調(diào)劑，材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|