 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (5631)
- >導(dǎo)師招生 (1085)
- >考研 (964)
- >蟲(chóng)友互識(shí) (688)
- >文獻(xiàn)求助 (351)
- >考博 (192)
- >休閑灌水 (178)
- >碩博家園 (97)
- >基金申請(qǐng) (69)
- >論文投稿 (60)
- >公派出國(guó) (59)
- >論文道賀祈福 (53)
- >教師之家 (49)
- >招聘信息布告欄 (47)
- >博后之家 (41)
- >SciFinder/Reaxys (38)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 1879 | 回復(fù): 2 | ||
1466036890銅蟲(chóng) (正式寫(xiě)手)
|
[求助]
求助利用molpro做CASPT2計(jì)算的輸入文件 已有1人參與
|
| 求采用molpro軟件的CASPT2輸入文件，想采用CASPT2計(jì)算無(wú)能壘反應(yīng)最低能量路徑上各點(diǎn)的單點(diǎn)能，之前只是采用molpro做過(guò)顯相關(guān)CC計(jì)算，從來(lái)沒(méi)有做過(guò)CASPT2計(jì)算，希望蟲(chóng)友們給予幫助。 |
回復(fù)此樓» 猜你喜歡
量子點(diǎn)電致發(fā)光器件（QLED）：藍(lán)光量子點(diǎn)（QDs）和納米氧化鋅電子傳輸材料
已經(jīng)有0人回復(fù)
-
需要一個(gè)圓偏振光誘導(dǎo)化學(xué)反應(yīng)，請(qǐng)問(wèn)使用圓偏振光作為光源，應(yīng)該需要哪些設(shè)備搭建
已經(jīng)有0人回復(fù)
-
物理化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有176人回復(fù)
-
[電子教材]天津大學(xué)物理化學(xué)（第七版）上、下冊(cè)
已經(jīng)有2人回復(fù)
-
陰離子交換膜電解二氧化碳還原
已經(jīng)有0人回復(fù)
-
九江學(xué)院2026年最新高層次人才招聘公告
已經(jīng)有0人回復(fù)
-
深圳信息職業(yè)技術(shù)學(xué)院2026年教師招聘公告（最新）
已經(jīng)有0人回復(fù)
-
求ENDNOTE破解軟件安裝包
已經(jīng)有4人回復(fù)
-
研究生招生
已經(jīng)有0人回復(fù)
-
河海大學(xué) （211.雙一流高校）鈣鈦礦光電實(shí)驗(yàn)室丁勇團(tuán)隊(duì)歡迎有志青年！
已經(jīng)有3人回復(fù)
-
黑龍江省原子能研究院化學(xué)課題組招收2026級(jí)化學(xué)專(zhuān)業(yè)或化學(xué)相關(guān)專(zhuān)業(yè)碩士生
已經(jīng)有0人回復(fù)
paramecium86
版主 (著名寫(xiě)手)
- QC強(qiáng)帖: 5
- 應(yīng)助: 587 (博士)
- 貴賓: 0.01
- 金幣: 26380.1
- 紅花: 349
- 帖子: 2025
- 在線(xiàn): 1442小時(shí)
- 蟲(chóng)號(hào): 2466443
- 注冊(cè): 2013-05-16
- 性別: GG
- 專(zhuān)業(yè): 理論和計(jì)算化學(xué)
- 管轄: 量子化學(xué)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
1466036890: 金幣+20, ★★★很有幫助, 已收到，非常感謝 2023-05-11 10:50:11
感謝參與，應(yīng)助指數(shù) +1
1466036890: 金幣+20, ★★★很有幫助, 已收到，非常感謝 2023-05-11 10:50:11
|
鏈接：https://pan.baidu.com/s/1yf4zrbVU8ZPnDrkKg-o-mQ?pwd=hgfk 提取碼：hgfk --來(lái)自百度網(wǎng)盤(pán)超級(jí)會(huì)員V9的分享 可以看看這個(gè)例子 之前算過(guò)的一個(gè)示例。 |
2樓2023-05-10 21:48:47
1234zou
金蟲(chóng) (小有名氣)
- 應(yīng)助: 7 (幼兒園)
- 金幣: 1941.5
- 紅花: 5
- 帖子: 92
- 在線(xiàn): 31.6小時(shí)
- 蟲(chóng)號(hào): 2914121
- 注冊(cè): 2014-01-04
- 性別: GG
- 專(zhuān)業(yè): 理論和計(jì)算化學(xué)
|
還有另一種途徑是不學(xué)習(xí)Molpro各種語(yǔ)法和功能，直接采用MOKIT程序黑箱、自動(dòng)式的多參考計(jì)算，這里以?xún)筛鵒-H鍵均被拉長(zhǎng)的水分子為例，寫(xiě)一個(gè)高斯gjf文件 %mem=4GB %nprocshared=2 #p CASPT2/cc-pVDZ mokit{CASPT2_prog=Molpro} 0 1 O -0.23497692 0.90193619 -0.068688 H 1.26502308 0.90193619 -0.068688 H -0.73568721 2.31589843 -0.068688 輸入文件使用最常用的gjf格式，想要什么理論方法直接寫(xiě)（哪怕高斯不支持該方法也行，因?yàn)椴⒉皇莾H調(diào)用高斯算的），符合書(shū)寫(xiě)直覺(jué)。中間CASPT2_prog=Molpro表示調(diào)用Molpro做CASPT2計(jì)算，一樣符合直覺(jué)。用automr提交，命令如下 automr h2o.gjf >h2o.out 2>&1 僅20秒結(jié)束，不用肉眼看軌道、挑選軌道、做軌道局域化、大小基組投影、擔(dān)心收斂性、調(diào)試關(guān)鍵詞等各種曲折操作，直接算出來(lái)就是CASSCF(4,4)和CASPT2(4,4)，程序會(huì)自動(dòng)確定活性軌道是兩根O-H的成鍵、反鍵軌道。h2o.out輸出文件里有各個(gè)方法的電子能量、雙自由基特征、未成對(duì)電子數(shù)等大量信息，直觀(guān)清楚 E(CASCI) = -75.90802109 a.u. E(CASSCF) = -75.90823063 a.u. ----------------------- Radical index ----------------------- biradical character (2c^2) y0= 0.165 tetraradical character(2c^2) y1= 0.139 Yamaguchi's unpaired electrons (sum_n n(2-n) ): 1.121 Head-Gordon's unpaired electrons(sum_n min(n,(2-n))): 0.607 Head-Gordon's unpaired electrons(sum_n (n(2-n))^2 ): 0.315 ------------------------------------------------------------- Leave subroutine do_cas at Sat Aug 12 14:46:05 2023 Enter subroutine do_mrpt2... CASPT2 based on CASSCF orbitals. Frozen_core = F, CASPT2 using program molpro $molpro -t 1 -n 2 -m 250m h2o_uhf_gvb4_CASSCF_CASPT2.com E(ref) = -75.90823031 a.u. IP-EA shift = 0.25 (default) E(corr) = -0.13847648 a.u. E(CASPT2) = -76.04670678 a.u. Leave subroutine do_mrpt2 at Sat Aug 12 14:46:16 2023 另有fch文件可用GaussView和Multiwfn直接打開(kāi)（不管你用什么量化程序算，都會(huì)給出fch文件，十分通用） CASSCF初始軌道 h2o_uhf_uno_asrot2gvb4_s.fch CASSCF收斂軌道，也是自然軌道 h2o_uhf_gvb4_CASSCF_NO.fch MOKIT程序免費(fèi)、開(kāi)源，見(jiàn)https://gitlab.com/jxzou/mokit 中文簡(jiǎn)介見(jiàn)https://gitlab.com/jxzou/mokit/-/blob/master/README_zh.md |
3樓2023-08-12 14:54:16
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 1466036890 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[文學(xué)芳草園] 伙伴們，祝我生日快樂(lè)吧 +3 | myrtle 2026-03-10 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-10 | 9/450 |
|
|
|
[考研] 274環(huán)境工程求調(diào)劑 +8 | 扶柳盈江 2026-03-05 | 9/450 |
|
|
|
[考研]
一志愿天津大學(xué)材料與化工275求調(diào)劑
 10+5
|
穿只靴子 2026-03-07 | 20/1000 |
|
|
|
[考研] 接受26屆調(diào)劑生 +22 | 豬豬豬毛 2026-03-06 | 23/1150 |
|
|
|
[考研] 288求調(diào)劑 +13 | 王曉陽(yáng)- 2026-03-09 | 18/900 |
|
|
|
[考研] 0857 資源與環(huán)境 285分 +5 | 未名考生 2026-03-09 | 5/250 |
|
|
|
[考研] 279求調(diào)劑 +3 | 莫xiao 2026-03-10 | 4/200 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-08 | 7/350 |
|
|
|
[考研] 一志愿：武漢理工，材料工程，英二數(shù)二 總分314 +3 | 2202020125 2026-03-10 | 4/200 |
|
|
|
[考研] 307求調(diào)劑 +3 | 辛仁豆腐 2026-03-08 | 5/250 |
|
|
|
[考研] 296求調(diào)劑 +4 | Xinyu Wu311 2026-03-09 | 4/200 |
|
|
|
[教師之家] 交大前校長(zhǎng)王樹(shù)國(guó)：現(xiàn)在最先進(jìn)的科技并不在大學(xué)實(shí)驗(yàn)室，而是在企業(yè)研究院 +4 | zju2000 2026-03-08 | 6/300 |
|
|
|
[考研] 求調(diào)劑推薦 +4 | 微辣不吃 2026-03-06 | 4/200 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-06 | 4/200 |
|
|
|
[考研] 268求調(diào)劑 +4 | 劉合華 2026-03-05 | 4/200 |
|
|
|
[考研] 材料調(diào)劑 +5 | ws 上岸鴨 2026-03-05 | 5/250 |
|
|
|
[考研] 材料調(diào)劑 +4 | L9370 2026-03-05 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +4 | 曼殊2266 2026-03-05 | 4/200 |
|
|
|
[考研] 一志愿西交化工專(zhuān)碩288專(zhuān)業(yè)課93求調(diào)劑 +6 | 好運(yùn)好運(yùn)接接 2026-03-04 | 6/300 |
|