 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2392 | 回復(fù): 5 | ||||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||||
hairan木蟲 (著名寫手)
瞌睡木蟲
|
[交流]
Gaussian和MolPro計(jì)算，對(duì)位于分子不同區(qū)域的同一種原子使用不同基組的定義方式 已有4人參與
|
|||||
|
1、計(jì)算的需求： 我們在計(jì)算中經(jīng)常會(huì)遇到一種情況，就是分子比較大但又不至于太大而必須一部分使用類似QM/MM的方式處理，此時(shí)我們希望對(duì)整個(gè)分子用同一種計(jì)算方法（比如HF或MP2或CCSD(T)等），但想對(duì)分子的不同區(qū)域使用不同的基組，這與平常我們根據(jù)元素符號(hào)來定義基組有些不同，我們需要使用原子默認(rèn)編號(hào)或自定義編號(hào)來定義基組。 我在網(wǎng)上搜索了一番，用Gaussian如何處理這個(gè)問題還有些答案，但不完全，用Molpro如何處理則完全沒有答案，現(xiàn)將如何使用這兩個(gè)軟件按照原子編號(hào)來定義基組分享給大家。 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 說明 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中稱作Lebel，大多數(shù)Gaussian使用者不知道的是原子編號(hào)是可以用來定義基組的。 2.3.3 下面說到關(guān)鍵部分，即輸入文件末尾的基組信息定義。 a. 最常用的情況：對(duì)不同的元素用不同的基組，此時(shí)具體基組上面一行出現(xiàn)的是元素符號(hào)，那一行以數(shù)字0結(jié)尾。 b. 需要不同區(qū)域的同一種原子用不同的基組時(shí)：可以將元素符號(hào)替換為原子編號(hào)。如果多個(gè)原子都用同一種基組可以寫在一行中，如果原子編號(hào)是連續(xù)的，可以用“-”來定義一個(gè)范圍，比如“7-16”表示編號(hào)從7到16的所有原子。 c. 注意事項(xiàng)：原子編號(hào)那一行仍然以數(shù)字“0”結(jié)尾；贗勢的定義中關(guān)于原子的這一行與基組的定義方式一致。 3、MolPro軟件的設(shè)置方式 與Gaussian不同的是，MolPro軟件對(duì)這個(gè)功能的支持還不是很完善，我也是與軟件作者H.-J. Werner教授通信之后才確定我的處理方式是可用的。 3.1 坐標(biāo)的預(yù)處理 需要專門定義基組的原子需要自己手動(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 請注意代表氫的兩個(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 說明 3.3.1 用MolPro處理時(shí)，可以將對(duì)應(yīng)了原子數(shù)最多的一種基組設(shè)置為默認(rèn)，就是上面的Default=AVTZ; 3.3.2 需要專門定義基組的原子依次寫下來，比如下面的“H1=VTZ”，每個(gè)一行，也可以寫在一行中，用逗號(hào)隔開； 3.3.3 注意事項(xiàng)： a. MolPro用這種方式處理時(shí)有個(gè)限制：基組定義部分好像不能超過50個(gè)字符，因此需要盡量用最省字符的方式； b. MolPro支持一些基組的簡寫，可以節(jié)省字符的數(shù)量，比如相關(guān)一致基組aug-cc-pVnZ基組可以簡寫為AVnZ，cc-pVnZ可以簡寫為VnZ （n代表zeta數(shù)目，可以為D, T, Q, 5, 6......）。 c. Werner教授說如此定義基組在計(jì)算結(jié)果方面是沒問題的，但如果你要在輸入出文件中輸出基組信息，可能會(huì)輸出的不對(duì)，比如當(dāng)前的輸入文件對(duì)應(yīng)的輸出文件中，只輸出了H1對(duì)應(yīng)的VTZ信息，沒有輸出H2的AVQZ的信息。MolPro的作者將去解決這個(gè)問題。 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 經(jīng)驗(yàn) | 計(jì)算化學(xué) | 自己感覺好的東西 |
» 本帖已獲得的紅花（最新10朵）
小范范1989» 猜你喜歡
請問四氫呋喃溶解的聚合物用甲醇沉淀時(shí)，如何使沉淀過程加速？
已經(jīng)有2人回復(fù)
-
七嗪類物質(zhì)合成求助
已經(jīng)有0人回復(fù)
-
物理化學(xué)論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有54人回復(fù)
-
Chemical Bonding at Surfaces and Interfaces,最經(jīng)典的一本表面上化學(xué)相互作用教材
已經(jīng)有0人回復(fù)
-
分子動(dòng)力學(xué)模擬合作需求
已經(jīng)有3人回復(fù)
-
誠邀有志之士加入江蘇大學(xué)環(huán)境學(xué)院合成生物學(xué)課題組！(長期接收學(xué)碩專碩、碩博連讀)
已經(jīng)有16人回復(fù)
-
一志愿南大理論與計(jì)算化學(xué)方向（070300）342分，求調(diào)劑
已經(jīng)有1人回復(fù)
-
歡迎化學(xué)、化工、物理、材料、或計(jì)算機(jī)背景的考研咨詢
已經(jīng)有0人回復(fù)
-
調(diào)劑招生
已經(jīng)有5人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- Multiwfn支持的弱相互作用的分析方法概覽
已經(jīng)有10人回復(fù)
- 大家有沒有遇到過用同樣方法算同一個(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
木蟲 (著名寫手)
- 應(yīng)助: 123 (高中生)
- 金幣: 3821.3
- 散金: 968
- 紅花: 243
- 帖子: 1979
- 在線: 452.8小時(shí)
- 蟲號(hào): 2013403
- 注冊: 2012-09-19
- 性別: GG
- 專業(yè): 原子和分子物理
5樓2014-12-26 16:23:29
761267
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 444.2
- 帖子: 83
- 在線: 61.3小時(shí)
- 蟲號(hào): 632122
- 注冊: 2008-10-21
- 專業(yè): 結(jié)構(gòu)化學(xué)
3樓2014-12-16 08:33:04
henrylv
木蟲 (著名寫手)
- 應(yīng)助: 6 (幼兒園)
- 金幣: 5078.6
- 散金: 2469
- 紅花: 8
- 帖子: 1008
- 在線: 1219.5小時(shí)
- 蟲號(hào): 537362
- 注冊: 2008-04-01
- 專業(yè): 海洋遙感
6樓2017-01-02 10:38:55
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 330分求調(diào)劑 +13 | qzenlc 2026-03-29 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 285求調(diào)劑 +5 | AZMK 2026-04-02 | 7/350 |
|
|
|
[考研] 312 化工或制藥調(diào)劑 +4 | 小小墨123 2026-04-02 | 5/250 |
|
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3 | 崔wj 2026-04-02 | 3/150 |
|
|
|
[考研] 279求調(diào)劑 +4 | 傅文秋 2026-04-02 | 4/200 |
|
|
|
[考研] 290求調(diào)劑 +3 | 1314捧花 2026-04-02 | 3/150 |
|
|
|
[考研] 261求B區(qū)調(diào)劑 +5 | 明仔· 2026-04-01 | 7/350 |
|
|
|
[考研] 0710生物學(xué)，325求調(diào)劑 +3 | mkkkkkl 2026-04-01 | 3/150 |
|
|
|
[考研] 化學(xué)工程專碩324分，一志愿中國礦業(yè)大學(xué)求調(diào)劑 +7 | 耿耿1314 2026-04-01 | 7/350 |
|
|
|
[考研] 085600，321分求調(diào)劑 +13 | 大饞小子 2026-03-31 | 13/650 |
|
|
|
[考研] 318求調(diào)劑 +8 | 七憶77 2026-04-01 | 8/400 |
|
|
|
[考研] 0855機(jī)械初試280求調(diào)劑 +3 | kazenotori 2026-03-31 | 3/150 |
|
|
|
[考研] 材料科學(xué)與工程求調(diào)劑 +13 | 深V宿舍吧 2026-03-29 | 13/650 |
|
|
|
[基金申請] 面上5B能上會(huì)嗎？ +8 | redcom 2026-03-29 | 8/400 |
|
|
|
[考研] 313求調(diào)劑 +6 | 賣個(gè)關(guān)子吧 2026-03-31 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +10 | 家佳佳佳佳佳 2026-03-29 | 10/500 |
|
|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
|
[考研] 316求調(diào)劑 +7 | 江辭666 2026-03-26 | 7/350 |
|
|
|
[考研] 265求調(diào)劑 +8 | 小木蟲085600 2026-03-27 | 8/400 |
|