 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 2343 | 回復(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» 猜你喜歡
意大利博洛尼亞大學(xué)招聘研究助理2名（2D材料/鋰離子與鈉離子電池方向）
已經(jīng)有1人回復(fù)
-
關(guān)于碳量子點(diǎn)抗菌的研究
已經(jīng)有15人回復(fù)
-
物理化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有165人回復(fù)
-
哈爾濱工業(yè)大學(xué)（深圳）綠色化學(xué)團(tuán)隊(duì)招收2026年秋季入學(xué)博士生
已經(jīng)有11人回復(fù)
-
南大化學(xué)調(diào)劑
已經(jīng)有11人回復(fù)
-
澳大利亞昆士蘭科技大學(xué)PhD博士獎(jiǎng)學(xué)金
已經(jīng)有0人回復(fù)
-
0703調(diào)劑
已經(jīng)有2人回復(fù)
-
1
已經(jīng)有20人回復(fù)
-
0703 調(diào)劑
已經(jīng)有5人回復(fù)
-
西南某211招27屆博士
已經(jīng)有23人回復(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
- 在線(xiàn): 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
- 在線(xiàn): 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
- 在線(xiàn): 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
- 在線(xiàn): 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
- 在線(xiàn): 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ā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +3 | 暗涌afhb 2026-03-16 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +4 | 梁梁校校 2026-03-19 | 4/200 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +3 | sbdksD 2026-03-19 | 3/150 |
|
|
|
[考研] 317求調(diào)劑 +3 | 申子申申 2026-03-19 | 6/300 |
|
|
|
[考研] 287求調(diào)劑 +3 | 晨昏線(xiàn)與星海 2026-03-19 | 4/200 |
|
|
|
[考研] 0817調(diào)劑 +3 | 沒(méi)有答案_ 2026-03-14 | 3/150 |
|
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過(guò)A區(qū)線(xiàn) +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
|
[考研] 085601專(zhuān)碩，總分342求調(diào)劑，地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
|
[考研] 297求調(diào)劑 +8 | 戲精丹丹丹 2026-03-17 | 8/400 |
|
|
|
[考博] 環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
|
[基金申請(qǐng)] 被我言中：新模板不強(qiáng)調(diào)格式了，假專(zhuān)家開(kāi)始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[考研] 材料專(zhuān)碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽(yáng)花天天開(kāi)心 2026-03-16 | 4/200 |
|
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
|
[考研] 機(jī)械專(zhuān)碩325，尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy！ 2026-03-13 | 6/300 |
|