 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 10 | 1/1 | 返回列表 |
| 查看: 1588 | 回復: 9 | ||
[求助]
用MOLPRO中的基組外推發(fā)現(xiàn)結果有躍變？ 已有2人參與
|
|
計算了一體系的勢能曲線，根據(jù)avqz，5z，6z的結果用molpro內置的外推方程外推，參考能量用E(n)=ECBS+A*exp(-C*n)，相關能用(n+ 0.5)^{-X} with X= 4.0。 發(fā)現(xiàn)qz，5z，6z各自的曲線都是光滑的，cbs結果除了一個小區(qū)域內的外推結果相對整個曲線的變化趨勢有明顯偏離外，其它地方均光滑而連續(xù)。如圖， Graph1.jpg Graph2.jpg 有沒有人也遇到過這種問題，應該怎樣解決？ |
回復此樓» 猜你喜歡
求催化作用基礎第三版 （甄開吉) 課件
已經有3人回復
-
26年材料/化學相關專業(yè)博士研究生招生 溫州大學化學與材料工程學院(26年秋季入學)
已經有0人回復
-
物理化學論文潤色/翻譯怎么收費?
已經有63人回復
-
哈爾濱工業(yè)大學（深圳）綠色化學團隊招收2026年秋季入學博士生
已經有11人回復
-
澳大利亞昆士蘭科技大學PhD博士獎學金
已經有0人回復
-
求助電催化劑制備，怎么在碳布表面生長PdCo合金
已經有0人回復
-
求助電催化劑制備，怎么在碳布表面生長PdCo合金
已經有0人回復
-
304求調劑
已經有3人回復
-
上海華東理工大學機械與動力工程學院先進材料與制造研究組招聘公告
已經有0人回復
-
304求調劑
已經有3人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 用MOLPRO構造原子自然軌道（ANO）基函數(shù)
已經有3人回復
- 高精度量子化學計算軟件Molpro 入門及理論精講課程
已經有3人回復
- Molpro 講習班 2012年第一期
已經有3人回復
luqing6879
木蟲 (著名寫手)
- 應助: 8 (幼兒園)
- 金幣: 3567.5
- 散金: 249
- 紅花: 3
- 帖子: 1907
- 在線: 808.1小時
- 蟲號: 614833
- 注冊: 2008-09-28
- 專業(yè): 理論和計算化學
2樓2014-05-10 00:51:43
3樓2014-05-10 22:06:55
孫1986
金蟲 (正式寫手)
- 應助: 1 (幼兒園)
- 金幣: 1046.9
- 散金: 1087
- 紅花: 19
- 帖子: 455
- 在線: 243.1小時
- 蟲號: 1405621
- 注冊: 2011-09-17
- 專業(yè): 原子和分子物理
4樓2014-05-11 12:48:32
【答案】應助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
kent1022: 金幣+20 2014-05-12 16:14:37
感謝參與，應助指數(shù) +1
kent1022: 金幣+20 2014-05-12 16:14:37
|
correlation不光滑，本質上還是因為reference不光滑。 不過，通常即使reference不光滑，加上correlation之后整體也是光滑的。所以，參考能和相關能分別做外推，可能不是一個good idea. 比如使用RHF或者UHF計算閉殼層體系的UCCSD能量，鍵稍微長一點，RHF和UHF是不一樣的，但最后總的UCCSD能量都是一樣的。 一般的，reference到5Z也就收斂了，QZ也基本收斂，所以可以對總的能量，用外推correlation能量的方法外推兩個最大基組的能量。 Science 342:1499，他們新的FH2勢能面，就是外推AV5Z和AV6Z的能量，因為HF能量在AV5Z已經收斂，所以可以跟相關能放在一起做外推。 [ 發(fā)自手機版 http://www.gaoyang168.com/3g ] |
5樓2014-05-11 13:55:08
luqing6879
木蟲 (著名寫手)
- 應助: 8 (幼兒園)
- 金幣: 3567.5
- 散金: 249
- 紅花: 3
- 帖子: 1907
- 在線: 808.1小時
- 蟲號: 614833
- 注冊: 2008-09-28
- 專業(yè): 理論和計算化學
|
加入相關能后是可以使能量曲線看上去更平滑。但是如virtuallzx之前在樓主mrci里的那個帖子所說，看上去光滑的片段在放大后可能還是有間斷。所以如4 5樓所判斷的，問題還是出現(xiàn)在hf。 首先看看那些點的t1是多少吧。如果大的話可能說明在該點的附近有臨近的激發(fā)態(tài)，做個2state casscf看看能量差多少。 至于基組外推，分開來推和合起來推我都試過，對于小分子來說，從絕對外推值來看，二者相差不大。 印象中分開推的人多一點。就我個人看，從本質上說，外推公式都是經驗公式，公式的形式必須要擬合能量曲線的趨勢，所以只要擬合的好，怎么推都可以。但問題是hf能和相關能的趨勢往往是不盡相同的。所以用不同的公式分別擬合二者，才更符合科學精神。合起來推的結果，由于相關比hf能小幾個數(shù)量級，所以相關能的信息完全被hf淹沒了。這也是為什么用兩種方法推差別不大的原因。 |
6樓2014-05-12 00:33:27
7樓2014-05-12 14:20:41
luqing6879
木蟲 (著名寫手)
- 應助: 8 (幼兒園)
- 金幣: 3567.5
- 散金: 249
- 紅花: 3
- 帖子: 1907
- 在線: 808.1小時
- 蟲號: 614833
- 注冊: 2008-09-28
- 專業(yè): 理論和計算化學
8樓2014-05-12 14:59:39
9樓2014-05-12 20:24:12
10樓2014-05-12 22:52:23
相關版塊跳轉
我要訂閱樓主 kent1022 的主題更新
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 化學求調劑 +3 | 臨澤境llllll 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 328求調劑，英語六級551，有科研經歷 +4 | 生物工程調劑 2026-03-16 | 12/600 |
|
|
|
[考研] 085410人工智能專碩317求調劑（0854都可以） +3 | xbxudjdn 2026-03-18 | 3/150 |
|
|
[考研] 311求調劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
|
[考研] 295求調劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
|
[考研]
 0854，計算機類招收調劑
+3
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 293求調劑 +11 | zjl的號 2026-03-16 | 16/800 |
|
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[基金申請] 被我言中：新模板不強調格式了，假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[考研] 一志愿蘇州大學材料工程（085601）專碩有科研經歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
|
[考研] 一志愿南京大學，080500材料科學與工程，調劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
|
[考研] 274求調劑 +5 | 時間點 2026-03-13 | 5/250 |
|
|
|
[考研] 機械專碩325，尋找調劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
|
[基金申請]
今年的國基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 一志愿211 0703方向310分求調劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
|
[考研] 327求調劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研] 本科南京大學一志愿川大藥學327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|