 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 2344 | 回復(fù): 5 | ||||||
hairan木蟲(chóng) (著名寫(xiě)手)
瞌睡木蟲(chóng)
|
[交流]
Gaussian和MolPro計(jì)算，對(duì)位于分子不同區(qū)域的同一種原子使用不同基組的定義方式 已有4人參與
|
|
1、計(jì)算的需求： 我們?cè)谟?jì)算中經(jīng)常會(huì)遇到一種情況，就是分子比較大但又不至于太大而必須一部分使用類(lèi)似QM/MM的方式處理，此時(shí)我們希望對(duì)整個(gè)分子用同一種計(jì)算方法（比如HF或MP2或CCSD(T)等），但想對(duì)分子的不同區(qū)域使用不同的基組，這與平常我們根據(jù)元素符號(hào)來(lái)定義基組有些不同，我們需要使用原子默認(rèn)編號(hào)或自定義編號(hào)來(lái)定義基組。 我在網(wǎng)上搜索了一番，用Gaussian如何處理這個(gè)問(wèn)題還有些答案，但不完全，用Molpro如何處理則完全沒(méi)有答案，現(xiàn)將如何使用這兩個(gè)軟件按照原子編號(hào)來(lái)定義基組分享給大家。 2、Gaussian軟件的設(shè)置方式： 2.1 關(guān)鍵詞：gen 2.2 輸入文件示例：甲醛分子，C和O用一種基組，兩個(gè)H分別用另外兩個(gè)不同的基組 ------------------------------------------------------------------------------------------ # MP2/gen MP2 calculation for formaldehyde 0 1 C 0.0000000000 0.0000000000 -0.5904227440 O 0.0000000000 0.0000000000 0.5915772560 H 0.0000000000 0.9327210617 -1.1773089520 H 0.0000000000 -0.9427210617 -1.2773089520 1 2 0 aug-cc-pVTZ **** 3 0 cc-pVTZ **** 4 0 aug-cc-pVQZ **** ------------------------------------------------------------------------------------------ 2.3 說(shuō)明 2.3.1 在輸入文件#打頭的那一行要有一個(gè)gen關(guān)鍵詞告訴程序去輸入文件結(jié)尾讀取基組信息，這與按照元素符號(hào)定義基組是一樣的。 2.3.2 一個(gè)需要知道的知識(shí)是Gaussian軟件默認(rèn)分子坐標(biāo)所在的行數(shù)當(dāng)做分子的編號(hào)，比如第三行的原子的編號(hào)為3，在GaussView中稱(chēng)作Lebel，大多數(shù)Gaussian使用者不知道的是原子編號(hào)是可以用來(lái)定義基組的。 2.3.3 下面說(shuō)到關(guān)鍵部分，即輸入文件末尾的基組信息定義。 a. 最常用的情況：對(duì)不同的元素用不同的基組，此時(shí)具體基組上面一行出現(xiàn)的是元素符號(hào)，那一行以數(shù)字0結(jié)尾。 b. 需要不同區(qū)域的同一種原子用不同的基組時(shí)：可以將元素符號(hào)替換為原子編號(hào)。如果多個(gè)原子都用同一種基組可以寫(xiě)在一行中，如果原子編號(hào)是連續(xù)的，可以用“-”來(lái)定義一個(gè)范圍，比如“7-16”表示編號(hào)從7到16的所有原子。 c. 注意事項(xiàng)：原子編號(hào)那一行仍然以數(shù)字“0”結(jié)尾；贗勢(shì)的定義中關(guān)于原子的這一行與基組的定義方式一致。 3、MolPro軟件的設(shè)置方式 與Gaussian不同的是，MolPro軟件對(duì)這個(gè)功能的支持還不是很完善，我也是與軟件作者H.-J. Werner教授通信之后才確定我的處理方式是可用的。 3.1 坐標(biāo)的預(yù)處理 需要專(zhuān)門(mén)定義基組的原子需要自己手動(dòng)定義編號(hào)，例如： C 0.0000000000 0.0000000000 -0.5904227440 O 0.0000000000 0.0000000000 0.5915772560 H1 0.0000000000 0.9327210617 -1.1773089520 H2 0.0000000000 -0.9427210617 -1.2773089520 請(qǐng)注意代表氫的兩個(gè)字母H后面分別加了一個(gè)數(shù)字，這個(gè)就是我們手動(dòng)定義的編號(hào)，這個(gè)編號(hào)將在基組的定義中用到。 3.2 輸入文件示例：甲醛分子，C和O用一種基組，兩個(gè)H分別用另外兩個(gè)不同的基組。 ------------------------------------------------------------------------------------------ *** Geomtyp=xyz geometry 4 MP2 calculation for formaldehyde C 0.0000000000 0.0000000000 -0.5904227440 O 0.0000000000 0.0000000000 0.5915772560 H1 0.0000000000 0.9327210617 -1.1773089520 H2 0.0000000000 -0.9427210617 -1.2773089520 } gprint,basis Basis Default=AVTZ H1=VTZ H2=AVQZ end HF MP2 ------------------------------------------------------------------------------------------ 3.3 說(shuō)明 3.3.1 用MolPro處理時(shí)，可以將對(duì)應(yīng)了原子數(shù)最多的一種基組設(shè)置為默認(rèn)，就是上面的Default=AVTZ; 3.3.2 需要專(zhuān)門(mén)定義基組的原子依次寫(xiě)下來(lái)，比如下面的“H1=VTZ”，每個(gè)一行，也可以寫(xiě)在一行中，用逗號(hào)隔開(kāi)； 3.3.3 注意事項(xiàng)： a. MolPro用這種方式處理時(shí)有個(gè)限制：基組定義部分好像不能超過(guò)50個(gè)字符，因此需要盡量用最省字符的方式； b. MolPro支持一些基組的簡(jiǎn)寫(xiě)，可以節(jié)省字符的數(shù)量，比如相關(guān)一致基組aug-cc-pVnZ基組可以簡(jiǎn)寫(xiě)為AVnZ，cc-pVnZ可以簡(jiǎn)寫(xiě)為VnZ （n代表zeta數(shù)目，可以為D, T, Q, 5, 6......）。 c. Werner教授說(shuō)如此定義基組在計(jì)算結(jié)果方面是沒(méi)問(wèn)題的，但如果你要在輸入出文件中輸出基組信息，可能會(huì)輸出的不對(duì)，比如當(dāng)前的輸入文件對(duì)應(yīng)的輸出文件中，只輸出了H1對(duì)應(yīng)的VTZ信息，沒(méi)有輸出H2的AVQZ的信息。MolPro的作者將去解決這個(gè)問(wèn)題。 |
回復(fù)此樓» 收錄本帖的淘帖專(zhuān)輯推薦
 經(jīng)驗(yàn) | 計(jì)算化學(xué) | 自己感覺(jué)好的東西 |
» 本帖已獲得的紅花（最新10朵）
小范范1989» 猜你喜歡
高效微納米儲(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é)論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有115人回復(fù)
-
化學(xué)化工學(xué)院 招收化工、化學(xué)、材料等相關(guān)方向研究生（學(xué)碩、專(zhuān)碩都有調(diào)劑名額）
已經(jīng)有1人回復(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ù)
-
高效微納米儲(chǔ)能材料與產(chǎn)業(yè)化推進(jìn)實(shí)驗(yàn)室 蘇州 招收碩士研究生（含調(diào)劑）
已經(jīng)有0人回復(fù)
-
溫州大學(xué)新能源材料微觀溯源課題組招收2026級(jí)博士研究生（材料、化學(xué)方向）
已經(jīng)有6人回復(fù)
-
材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流
已經(jīng)有5人回復(fù)
-
26申博
已經(jīng)有2人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- Multiwfn支持的弱相互作用的分析方法概覽
已經(jīng)有10人回復(fù)
- 大家有沒(méi)有遇到過(guò)用同樣方法算同一個(gè)東西，兩次結(jié)果不一樣的？
已經(jīng)有12人回復(fù)
- Molpro和Gaussian算cas，我該信任那一個(gè)呢？
已經(jīng)有7人回復(fù)
- 【討論】molpro, molcas和gaussian的優(yōu)缺點(diǎn)
已經(jīng)有3人回復(fù)
小范范1989
木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 123 (高中生)
- 金幣: 3821.3
- 散金: 968
- 紅花: 243
- 帖子: 1979
- 在線: 452.8小時(shí)
- 蟲(chóng)號(hào): 2013403
- 注冊(cè): 2012-09-19
- 性別: GG
- 專(zhuān)業(yè): 原子和分子物理
2樓2014-12-11 08:12:36
761267
金蟲(chóng) (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 429.7
- 帖子: 83
- 在線: 61.3小時(shí)
- 蟲(chóng)號(hào): 632122
- 注冊(cè): 2008-10-21
- 專(zhuān)業(yè): 結(jié)構(gòu)化學(xué)
3樓2014-12-16 08:33:04
kong890109
金蟲(chóng) (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 976.4
- 紅花: 10
- 帖子: 139
- 在線: 133.6小時(shí)
- 蟲(chóng)號(hào): 2382781
- 注冊(cè): 2013-03-28
- 性別: GG
- 專(zhuān)業(yè): 化學(xué)動(dòng)力學(xué)
4樓2014-12-16 08:49:21
小范范1989
木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 123 (高中生)
- 金幣: 3821.3
- 散金: 968
- 紅花: 243
- 帖子: 1979
- 在線: 452.8小時(shí)
- 蟲(chóng)號(hào): 2013403
- 注冊(cè): 2012-09-19
- 性別: GG
- 專(zhuān)業(yè): 原子和分子物理
5樓2014-12-26 16:23:29
henrylv
木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 6 (幼兒園)
- 金幣: 5078.6
- 散金: 2469
- 紅花: 8
- 帖子: 1008
- 在線: 1219.5小時(shí)
- 蟲(chóng)號(hào): 537362
- 注冊(cè): 2008-04-01
- 專(zhuān)業(yè): 海洋遙感
6樓2017-01-02 10:38:55
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 hairan 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[論文投稿]
申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變，給編輯又寫(xiě)郵件了，沒(méi)回復(fù)
 10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +4 | 脫穎而出 2026-03-16 | 12/600 |
|
|
|
[考研]
 271材料工程求調(diào)劑
+6
|
.6lL 2026-03-18 | 6/300 |
|
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
|
[考研] 295求調(diào)劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
|
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +7 | 手機(jī)用戶(hù) 2026-03-17 | 7/350 |
|
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
|
[考研] 293求調(diào)劑 +11 | zjl的號(hào) 2026-03-16 | 16/800 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專(zhuān)業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
|
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂(lè) 2026-03-14 | 5/250 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|