 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 10 | 1/1 | 返回列表 |
| 查看: 1589 | 回復: 9 | ||
[求助]
用MOLPRO中的基組外推發(fā)現(xiàn)結果有躍變？ 已有2人參與
|
|
計算了一體系的勢能曲線，根據(jù)avqz，5z，6z的結果用molpro內(nèi)置的外推方程外推，參考能量用E(n)=ECBS+A*exp(-C*n)，相關能用(n+ 0.5)^{-X} with X= 4.0。 發(fā)現(xiàn)qz，5z，6z各自的曲線都是光滑的，cbs結果除了一個小區(qū)域內(nèi)的外推結果相對整個曲線的變化趨勢有明顯偏離外，其它地方均光滑而連續(xù)。如圖， Graph1.jpg Graph2.jpg 有沒有人也遇到過這種問題，應該怎樣解決？ |
回復此樓» 猜你喜歡
CSC訪學博后項目獲批，外方學校暫停合作該怎么辦？
已經(jīng)有51人回復
-
[電子教材]天津大學物理化學（第七版）上、下冊
已經(jīng)有3人回復
-
物理化學論文潤色/翻譯怎么收費?
已經(jīng)有73人回復
-
半透性容器進行高濕或者低濕研究
已經(jīng)有0人回復
-
陰離子交換膜電解二氧化碳還原
已經(jīng)有0人回復
-
分子動力學模擬合作需求
已經(jīng)有2人回復
-
九江學院2026年最新高層次人才招聘公告
已經(jīng)有0人回復
-
深圳信息職業(yè)技術學院2026年教師招聘公告（最新）
已經(jīng)有0人回復
-
研究生招生
已經(jīng)有0人回復
-
山西大學資源與環(huán)境工程研究所化工分離/環(huán)境功能材料小組招收2026年博士生
已經(jīng)有2人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 用MOLPRO構造原子自然軌道（ANO）基函數(shù)
已經(jīng)有3人回復
- 高精度量子化學計算軟件Molpro 入門及理論精講課程
已經(jīng)有3人回復
- Molpro 講習班 2012年第一期
已經(jīng)有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不光滑，本質(zhì)上還是因為reference不光滑。 不過，通常即使reference不光滑，加上correlation之后整體也是光滑的。所以，參考能和相關能分別做外推，可能不是一個good idea. 比如使用RHF或者UHF計算閉殼層體系的UCCSD能量，鍵稍微長一點，RHF和UHF是不一樣的，但最后總的UCCSD能量都是一樣的。 一般的，reference到5Z也就收斂了，QZ也基本收斂，所以可以對總的能量，用外推correlation能量的方法外推兩個最大基組的能量。 Science 342:1499，他們新的FH2勢能面，就是外推AV5Z和AV6Z的能量，因為HF能量在AV5Z已經(jīng)收斂，所以可以跟相關能放在一起做外推。 [ 發(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看看能量差多少。 至于基組外推，分開來推和合起來推我都試過，對于小分子來說，從絕對外推值來看，二者相差不大。 印象中分開推的人多一點。就我個人看，從本質(zhì)上說，外推公式都是經(jīng)驗公式，公式的形式必須要擬合能量曲線的趨勢，所以只要擬合的好，怎么推都可以。但問題是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
相關版塊跳轉(zhuǎn)
我要訂閱樓主 kent1022 的主題更新
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿蘇州大學材料求調(diào)劑，總分315（英一） +3 | sbdksD 2026-03-19 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 296求調(diào)劑 +3 | www_q 2026-03-18 | 6/300 |
|
|
|
[考研] 266求調(diào)劑 +5 | 陽陽哇塞 2026-03-14 | 10/500 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-18 | 3/150 |
|
|
|
[考研] 085600材料與化工求調(diào)劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
|
[考研] 一志愿天大材料與化工（085600）總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
|
[考研] 0703化學 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
|
[考研] 0703化學調(diào)劑，求各位老師收留 +10 | 秋有木北 2026-03-14 | 10/500 |
|
|
|
[考研] 085600材料與化工 +5 | 安全上岸！ 2026-03-16 | 5/250 |
|
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
|
[考研] 070300化學319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
|
[考研] 0703化學求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
|
[基金申請] 被我言中：新模板不強調(diào)格式了，假專家開始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[考研] 一志愿蘇州大學材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學123 2026-03-15 | 8/400 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|