 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2455 | 回復(fù): 4 | ||
daniell251木蟲 (著名寫手)
|
[求助]
采用PBE方法計(jì)算 如何輸入？
|
|
一篇文獻(xiàn)中說明了他們所采用的方法： ...Calculations used the exchange and correlation functionals of Perdew, Bueke, and Ernzerhof, and the TZVP basis set of Godbelt, Andzelm, and co-workers for carbon, hydrogen, oxygen and phosphorus. Calculations on gold employ the Stuttart 97 effective core potential and basis...... 我采用了如下的輸入 結(jié)果報(bào)錯(cuò) 說有# opt=modredundant freq pbe/genecp 這一行存在語法錯(cuò)誤 哪位達(dá)人能夠幫助我解答嗎 還有 我的整個(gè)輸入是不是都存在問題 謝謝啦 # opt=modredundant freq pbe/genecp 文件名 1 1 結(jié)構(gòu)坐標(biāo) Au 0 sddall **** C H P O 0 tzvp **** Au sddall ------------- 謝謝各位了 --------------- 問題貌似解決了 應(yīng)該是在輸入PBE 方法時(shí) 不能只輸入一次 pbe 而是要采用 pbepbe . 可為什么會(huì)是這樣呢 還有 我的其他輸入有沒有問題 謝謝大家了 [ Last edited by daniell251 on 2012-9-4 at 02:36 ] |
回復(fù)此樓» 猜你喜歡
PbS量子點(diǎn)如何確定濃度
已經(jīng)有11人回復(fù)
-
細(xì)胞培養(yǎng)，這22個(gè)細(xì)節(jié)一定要注意�。ㄒ唬�
已經(jīng)有0人回復(fù)
-
物理化學(xué)論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有288人回復(fù)
-
CSC訪學(xué)博后項(xiàng)目獲批，外方學(xué)校暫停合作該怎么辦？
已經(jīng)有51人回復(fù)
-
大灣區(qū)大學(xué)劉天輝課題組招收2026級(jí)博士生（與中山大學(xué)物理學(xué)院聯(lián)合培養(yǎng)）
已經(jīng)有22人回復(fù)
-
PbS量子點(diǎn)紫外吸收
已經(jīng)有0人回復(fù)
-
比利時(shí)魯汶大學(xué)與國家留學(xué)基金委共同資助博士研究生CSC-KU Leuven PhD Scholarship
已經(jīng)有0人回復(fù)
-
倫敦大學(xué)學(xué)院Benjamin Moss博士招收人工智能、原位光譜、催化相關(guān)26年博士生
已經(jīng)有11人回復(fù)
-
分子動(dòng)力學(xué)模擬合作需求
已經(jīng)有2人回復(fù)
-
中國科學(xué)院成都山地災(zāi)害與環(huán)境研究所2026年特別研究助理（博士后）招聘啟事
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 如何采用castep計(jì)算E-V曲線，請專家指示
已經(jīng)有20人回復(fù)
- =.= Zr_sv的PBE-PAW贗勢在計(jì)算中出現(xiàn)了個(gè)跳躍點(diǎn)……
已經(jīng)有21人回復(fù)
- matlab反求未知輸入?yún)?shù)
已經(jīng)有4人回復(fù)
- 求助 Origin 自定義函數(shù)輸入問題
已經(jīng)有4人回復(fù)
- 如何計(jì)算水凝膠的藥物緩釋測試方法和緩釋量計(jì)算方法
已經(jīng)有7人回復(fù)
- 如何產(chǎn)生Sr的自旋-軌道耦合的PBE贗勢
已經(jīng)有12人回復(fù)
- 如何在origin中輸入一個(gè)函數(shù)（比較復(fù)雜的函數(shù)），然后得到曲線
已經(jīng)有13人回復(fù)
- 【求助】PBE0 計(jì)算Si、Ge能帶結(jié)構(gòu)
已經(jīng)有9人回復(fù)
youji
木蟲 (正式寫手)
- 應(yīng)助: 7 (幼兒園)
- 金幣: 3588.2
- 散金: 3481
- 紅花: 4
- 帖子: 485
- 在線: 630.9小時(shí)
- 蟲號(hào): 408052
- 注冊: 2007-06-20
- 專業(yè): 能源化工
2樓2012-09-05 14:22:13
daniell251
木蟲 (著名寫手)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 2398.9
- 紅花: 7
- 帖子: 1215
- 在線: 298.1小時(shí)
- 蟲號(hào): 681951
- 注冊: 2008-12-26
- 性別: GG
- 專業(yè): 金屬有機(jī)化學(xué)
3樓2012-09-06 00:35:48
mearlt
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 7.3
- 帖子: 16
- 在線: 19.7小時(shí)
- 蟲號(hào): 2285575
- 注冊: 2013-02-17
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
4樓2015-01-23 21:43:15
小范范1989
木蟲 (著名寫手)
- 應(yīng)助: 123 (高中生)
- 金幣: 3821.3
- 散金: 968
- 紅花: 243
- 帖子: 1979
- 在線: 452.8小時(shí)
- 蟲號(hào): 2013403
- 注冊: 2012-09-19
- 性別: GG
- 專業(yè): 原子和分子物理
|
贗勢基組 簡單來說，贗勢就是不計(jì)算內(nèi)層電子，而是把內(nèi)層電子的貢獻(xiàn)用一個(gè)勢來描述，放在哈密頓里面。適用于重元素。 http://blog.163.com/ares_119/blog/static/1157851132009513378305/ 贗勢基組，實(shí)際上包括 贗勢和基組兩個(gè)部分，內(nèi)層電子采用贗勢，即effective core potential (ECP)，外層價(jià)電子采用一般的基組。 比如： LanL2DZ: D95V on first row, Los Alamos ECP plus DZ on Na-Bi. 就是對第一行原子是 D95V (這個(gè)是非贗勢基組)，對Na-Bi是使用一個(gè)叫做Los Alamos的有效核勢加上一個(gè)DZ基組。所以Lanl2dz就是對前面的原子全電子基組，對后面的原子是贗勢基組。 （再次說明，量化里面，C，O那一行，算周期表的第一行） 使用贗勢的3個(gè)原因： 1。沒有相應(yīng)的全電子基組。 2。減少計(jì)算量。 3。贗勢可以包含重金屬相對論效應(yīng)的修正。 在高斯中，lanl2dz基組，在手冊中可以查到其定義為： LanL2DZ: D95V on first row , Los Alamos ECP plus DZ on Na-Bi，也就是說，對于C，O等元素來說（量化中，H大概算第0周期，C，O才是第一周期），Lanl2dz實(shí)際上還是全電子基組，而對于Na以后才是對內(nèi)層電子用Los Alamos ECP贗勢，外層電子用DZ 基組。 使用贗勢的輸入文件： 1.所有原子使用lanl2dz -------------------------- #HF/lanl2dz opt lanl2dz for all atoms 0 2 O 0.0 0.0 0.0 C 0.0 0.0 1.2 Cu 0.0 0.0 3.2 -------------------------- 2.所有原子使用lanl2dz的另一種輸入方法。 -------------------------- #HF/genecp opt lanl2dz for all atoms 0 2 O 0.0 0.0 0.0 C 0.0 0.0 1.2 Cu 0.0 0.0 3.2 C O 0 lanl2dz **** Cu 0 lanl2dz //定義價(jià)電子的基組, C O 0 是碳，氧，零，其中零用作終止符號(hào)。 **** Cu 0 lanl2dz //定義內(nèi)層電子的贗勢 -------------------------- 3.混和基組，即有的使用全電子，有的使用Lanl2dz。格式同2。 #HF/genecp opt lanl2dz for Cu, 6-31G(d) for C and O 0 2 O 0.0 0.0 0.0 C 0.0 0.0 1.2 Cu 0.0 0.0 3.2 C O 0 6-31G(d) //另一種全電子基組 **** Cu 0 lanl2dz //定義價(jià)電子的基組 **** Cu 0 lanl2dz //定義內(nèi)層電子的贗勢 Gaussian和MolPro計(jì)算，對位于分子不同區(qū)域的同一種原子使用不同基組的定義方式 http://www.gaoyang168.com/bbs/viewthread.php?tid=8278533&pid=2#pid2 1、計(jì)算的需求： 我們在計(jì)算中經(jīng)常會(huì)遇到一種情況，就是分子比較大但又不至于太大而必須一部分使用類似QM/MM的方式處理，此時(shí)我們希望對整個(gè)分子用同一種計(jì)算方法（比如HF或MP2或CCSD(T)等），但想對分子的不同區(qū)域使用不同的基組，這與平常我們根據(jù)元素符號(hào)來定義基組有些不同，我們需要使用原子默認(rèn)編號(hào)或自定義編號(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. 最常用的情況：對不同的元素用不同的基組，此時(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)于原子的這一行與基組的定義方式一致。 |
5樓2015-01-24 11:16:39
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 daniell251 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 322求調(diào)劑 +4 | X123- 2026-03-14 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑 +11 | 妮妮ninicgb 2026-03-15 | 15/750 |
|
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
|
[考研] 初試 317 +6 | 半拉月丙 2026-03-20 | 6/300 |
|
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +3 | 北道巷 2026-03-18 | 3/150 |
|
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
|
[考研] 299求調(diào)劑 +4 | 某某某某位 2026-03-21 | 4/200 |
|
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
|
[考研] 316求調(diào)劑 +6 | 梁茜雯 2026-03-19 | 6/300 |
|
|
|
[考研] 初始318分求調(diào)劑（有工作經(jīng)驗(yàn)） +3 | 1911236844 2026-03-17 | 3/150 |
|
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
|
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +9 | 大可digkids 2026-03-16 | 9/450 |
|
|
|
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|