 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 1880 | 回復(fù): 2 | |||
1466036890銅蟲 (正式寫手)
|
[求助]
求助利用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ì)算，希望蟲友們給予幫助。 |
回復(fù)此樓» 猜你喜歡
cif求助下載
已經(jīng)有1人回復(fù)
-
高效微納米儲(chǔ)能材料與產(chǎn)業(yè)化推進(jìn)實(shí)驗(yàn)室 蘇州 招收碩士研究生（含調(diào)劑）
已經(jīng)有5人回復(fù)
-
物理化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有105人回復(fù)
-
基金正文30頁(yè)指的是報(bào)告正文還是整個(gè)申請(qǐng)書
已經(jīng)有9人回復(fù)
-
高效微納米儲(chǔ)能材料與產(chǎn)業(yè)化推進(jìn)實(shí)驗(yàn)室 蘇州 招收碩士研究生（含調(diào)劑）
已經(jīng)有0人回復(fù)
-
高效微納米儲(chǔ)能材料與產(chǎn)業(yè)化推進(jìn)實(shí)驗(yàn)室 蘇州 招收碩士研究生（含調(diào)劑）
已經(jīng)有2人回復(fù)
-
化學(xué)化工學(xué)院 招收化工、化學(xué)、材料等相關(guān)方向研究生（學(xué)碩、專碩都有調(diào)劑名額）
已經(jīng)有1人回復(fù)
-
高效微納米儲(chǔ)能材料與產(chǎn)業(yè)化推進(jìn)實(shí)驗(yàn)室 蘇州 招收碩士研究生（含調(diào)劑）
已經(jīng)有0人回復(fù)
-
溫州大學(xué)新能源材料微觀溯源課題組招收2026級(jí)博士研究生（材料、化學(xué)方向）
已經(jīng)有6人回復(fù)
paramecium86
版主 (著名寫手)
- QC強(qiáng)帖: 5
- 應(yīng)助: 587 (博士)
- 貴賓: 0.01
- 金幣: 26380.1
- 紅花: 349
- 帖子: 2025
- 在線: 1442小時(shí)
- 蟲號(hào): 2466443
- 注冊(cè): 2013-05-16
- 性別: GG
- 專業(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)盤超級(jí)會(huì)員V9的分享 可以看看這個(gè)例子 之前算過(guò)的一個(gè)示例。 |
2樓2023-05-10 21:48:47
1234zou
金蟲 (小有名氣)
- 應(yīng)助: 7 (幼兒園)
- 金幣: 1941.5
- 紅花: 5
- 帖子: 92
- 在線: 31.6小時(shí)
- 蟲號(hào): 2914121
- 注冊(cè): 2014-01-04
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
|
還有另一種途徑是不學(xué)習(xí)Molpro各種語(yǔ)法和功能，直接采用MOKIT程序黑箱、自動(dòng)式的多參考計(jì)算，這里以兩根O-H鍵均被拉長(zhǎng)的水分子為例，寫一個(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格式，想要什么理論方法直接寫（哪怕高斯不支持該方法也行，因?yàn)椴⒉皇莾H調(diào)用高斯算的），符合書寫直覺(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ù)等大量信息，直觀清楚 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直接打開（不管你用什么量化程序算，都會(huì)給出fch文件，十分通用） CASSCF初始軌道 h2o_uhf_uno_asrot2gvb4_s.fch CASSCF收斂軌道，也是自然軌道 h2o_uhf_gvb4_CASSCF_NO.fch MOKIT程序免費(fè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ā)表 | |
|---|---|---|---|---|
|
[考研] 材料工程,326分，求調(diào)劑 +5 | KRSLSR 2026-03-10 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-11 | 4/200 |
|
|
[考博] 售SCI一區(qū)T0P文章，我:8.O.5.5.1.O.5.4,科目齊全,可+急 +3 | 7x4iz20zqa 2026-03-11 | 8/400 |
|
|
|
[考研] 一志愿中科大材料與化工。353分能調(diào)劑到哪？ +5 | MayUxw1 2026-03-08 | 7/350 |
|
|
|
[考研] 311求調(diào)劑 +6 | zchqwer 2026-03-10 | 6/300 |
|
|
|
[考研] 材料工程專碩，一志愿中國(guó)礦業(yè)大學(xué)，總分314，求調(diào)劑 +4 | 無(wú)懈可擊的巨人 2026-03-10 | 4/200 |
|
|
|
[考研] 265求調(diào)劑 +8 | 小木蟲085600 2026-03-09 | 11/550 |
|
|
|
[考研] 2026考研求調(diào)劑-材料類-本科211一志愿985-初試301分 +10 | 蟲友233 2026-03-07 | 10/500 |
|
|
|
[考研] 085602化工求調(diào)劑 +7 | 董boxing 2026-03-10 | 7/350 |
|
|
|
[考研] 328，0703考生求調(diào)劑，一志愿為東北師范大學(xué) +3 | 觀素律 2026-03-09 | 4/200 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 三dd. 2026-03-10 | 3/150 |
|
|
|
[考研] 08工科 +5 | li李樂(lè)成 2026-03-06 | 5/250 |
|
|
|
[考研] 一志愿天大化工（085600）調(diào)劑總分338 +5 | 蔡大美女 2026-03-09 | 5/250 |
|
|
|
[考研] 考研一志愿長(zhǎng)安大學(xué)材料與化工309分請(qǐng)求調(diào)劑 +6 | dtdxzxx 2026-03-06 | 8/400 |
|
|
|
[考研] 家人們 調(diào)劑不迷路 看這里 +8 | likeihood 2026-03-09 | 13/650 |
|
|
|
[考研] 0832食品科學(xué)與工程293調(diào)劑 +3 | 東東不東 2026-03-07 | 3/150 |
|
|
|
[考研] 求調(diào)劑，數(shù)一英一274分 +4 | 小菲會(huì)努力 2026-03-08 | 4/200 |
|
|
|
[考研] 296求調(diào)劑 +4 | Xinyu Wu311 2026-03-09 | 4/200 |
|
|
|
[考博] 26年博士申請(qǐng) +4 | 科研狗111 2026-03-07 | 4/200 |
|
|
|
[考研] 303求調(diào)劑 +8 | forgman95 2026-03-05 | 10/500 |
|