 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2104 | 回復(fù): 11 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
xiao_qft新蟲 (初入文壇)
|
[求助]
為何三線態(tài)氧分子(Triplet-O2)勢能面掃描PES曲線存在斷點(diǎn)？ 已有1人參與
|
||
|
輸入文件為（三線態(tài)）： %mem=1920MB %nprocshared=8 #T B3LYP/aug-cc-pVQZ Scan Test O2 PES Scan 0 3 O O 1 R R 0.8 84 0.05 單線態(tài)僅3改為1（其它不變）： %mem=1920MB %nprocshared=8 #T B3LYP/aug-cc-pVQZ Scan Test O2 PES Scan 0 1 O O 1 R R 0.8 84 0.05 結(jié)果附件圖片。 哪兒出的問題？謝謝 O2-PES.png |
回復(fù)此樓» 猜你喜歡
意大利博洛尼亞大學(xué)招聘研究助理2名（2D材料/鋰離子與鈉離子電池方向）
已經(jīng)有1人回復(fù)
-
關(guān)于碳量子點(diǎn)抗菌的研究
已經(jīng)有15人回復(fù)
-
物理化學(xué)論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有188人回復(fù)
-
南大化學(xué)調(diào)劑
已經(jīng)有11人回復(fù)
-
0703調(diào)劑
已經(jīng)有4人回復(fù)
-
1
已經(jīng)有21人回復(fù)
-
0703 調(diào)劑
已經(jīng)有5人回復(fù)
-
山西大學(xué)資源與環(huán)境工程研究所化工分離/環(huán)境功能材料小組招收2026年博士生
已經(jīng)有2人回復(fù)
-
西南某211招27屆博士
已經(jīng)有30人回復(fù)
-
大連工業(yè)杰青/長江團(tuán)隊(duì)-生物質(zhì)材料-儲(chǔ)能電池方向招收2026級(jí)博士生
已經(jīng)有0人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
virtualzx
木蟲 (著名寫手)
- QC強(qiáng)帖: 2
- 應(yīng)助: 263 (大學(xué)生)
- 金幣: 7161.3
- 紅花: 54
- 帖子: 1605
- 在線: 317.6小時(shí)
- 蟲號(hào): 2069080
- 注冊: 2012-10-18
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
【答案】應(yīng)助回帖
|
這個(gè)結(jié)果你之前計(jì)算得到的比好得多了，至少結(jié)果是定性正確的。不論單線態(tài)或三線態(tài)，RODFT或是CASSCF是可以收斂到定性正確的結(jié)果的，但是都需要有足夠好的初始猜測。 進(jìn)行掃描之前最好要先進(jìn)行單點(diǎn)計(jì)算驗(yàn)證結(jié)果是正確的，也就是進(jìn)行一個(gè)單點(diǎn)計(jì)算，然后查看軌道，檢查是否得到正確的電子組態(tài)；收斂不代表結(jié)果正確。檢查軌道正確后，直接讀入這組軌道作為第一步的初始猜測，然后從這里開始掃描。開殼層單線態(tài)往往不易得到正確結(jié)果，這時(shí)用三線態(tài)收斂得到的軌道作為初始猜測一般會(huì)得到比較好的結(jié)果。 O2解離至少需要8電子6軌道，但是最好12電子8軌道。就如牛老師說的，如果你要進(jìn)行很準(zhǔn)確的計(jì)算，最好使用更專業(yè)一些的程序進(jìn)行多組態(tài)電子相關(guān)計(jì)算。 [ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ] |
7樓2015-08-03 22:49:51
xiao_qft
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 381
- 帖子: 22
- 在線: 7.2小時(shí)
- 蟲號(hào): 3927172
- 注冊: 2015-06-16
- 專業(yè): 納米醫(yī)學(xué)
|
從log文件中看到斷點(diǎn)前后有如下信息： 斷點(diǎn)前： <Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 1.0000 <S**2>= 2.0153 S= 1.0051 <L.S>= 0.000000000000E+00 Annihilation of the first spin contaminant: S**2 before annihilation 2.0153, after 2.0001 斷點(diǎn)后： <Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 1.0000 <S**2>= 3.0138 S= 1.3066 <L.S>= 0.000000000000E+00 Annihilation of the first spin contaminant: S**2 before annihilation 3.0138, after 2.0279 該如何處理？ |
2樓2015-07-31 14:44:09
virtualzx
木蟲 (著名寫手)
- QC強(qiáng)帖: 2
- 應(yīng)助: 263 (大學(xué)生)
- 金幣: 7161.3
- 紅花: 54
- 帖子: 1605
- 在線: 317.6小時(shí)
- 蟲號(hào): 2069080
- 注冊: 2012-10-18
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
xiao_qft: 金幣+10, ★★★★★最佳答案 2015-08-04 16:53:50
感謝參與，應(yīng)助指數(shù) +1
xiao_qft: 金幣+10, ★★★★★最佳答案 2015-08-04 16:53:50
|
這里你的兩個(gè)曲線都不對(duì)，O2解離能應(yīng)該是0.26hartree左右，這倆都太高。三線態(tài)和單線態(tài)在解離處應(yīng)該是簡并的，你這里能量不同。 單線態(tài)是閉殼層計(jì)算，因此沒法正確描述解離后的開殼層結(jié)構(gòu)（解離后有四個(gè)未成對(duì)電子） 用unrestrict 方法進(jìn)行解離計(jì)算經(jīng)常會(huì)得到由于自旋污染造成的間斷點(diǎn)，即使是基態(tài)也會(huì)有斷點(diǎn)，稱為 Coulson–Fischer點(diǎn)(Phys. Rev. B 13, 4274)，這個(gè)點(diǎn)之后出現(xiàn)對(duì)稱破缺解。因此最好用Restricted Open-shell來算。當(dāng)然要真正正確的結(jié)果就最好進(jìn)行多組態(tài)計(jì)算。 另外鍵解離要從遠(yuǎn)到近掃描，否則即使沒有自旋耦合，斷點(diǎn)也很常見，即使沒有斷點(diǎn)，也往往得到錯(cuò)誤的能量曲線。掃描時(shí)是用前一個(gè)點(diǎn)的軌道作為下一個(gè)點(diǎn)的初始猜測的，解離時(shí)化學(xué)鍵重點(diǎn)電子從成對(duì)到不成對(duì)，程序往往無法及時(shí)發(fā)現(xiàn)。所以從遠(yuǎn)到近掃描。首先在遠(yuǎn)處進(jìn)行單點(diǎn)計(jì)算，查看軌道，正確之后從這里開始向短處進(jìn)行掃描。 [ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ] |
3樓2015-08-01 01:35:01
xiao_qft
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 381
- 帖子: 22
- 在線: 7.2小時(shí)
- 蟲號(hào): 3927172
- 注冊: 2015-06-16
- 專業(yè): 納米醫(yī)學(xué)
|
本人量化新手，您的回復(fù)是否表明對(duì)于O2分子，簡單的Gaussian Scan關(guān)鍵字是無法得到正確的PES的？ 按您的提示，采用由遠(yuǎn)及近掃描、以及改為ROB3LYP關(guān)鍵字 對(duì)于單線態(tài)O2： 1、B3LYP&正向掃描 2、B3LYP&反向掃描 3、ROB3LYP&反向掃描 三個(gè)方法獲得的PES完全重合 對(duì)于三線態(tài)O2： 1、B3LYP&正向掃描 2、B3LYP&反向掃描 兩個(gè)方法獲得的PES僅僅是間斷點(diǎn)移位，如下圖 3、ROB3LYP&反向掃描 無法獲得結(jié)果，見下log文件 Convergence failure -- run terminated. Error termination via Lnk1e in /home/gate/GAUSSIAN/g09/l502.exe “O2解離能應(yīng)該是0.26hartree左右” 至少對(duì)于基態(tài)Triplet-O2，Gaussian 輸入文件該如何編寫才能得到正確的結(jié)果？ 謝謝！  O2-Tripelet_Reverse_Scan.png |
» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : O2T_PES_03.log
2015-08-03 14:14:46, 17.29 K
4樓2015-08-03 14:18:31
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +5 | 困困困困坤坤 2026-03-20 | 6/300 |
|
|
|
[考研] 287求調(diào)劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
[考研] 考研調(diào)劑 +4 | 來好運(yùn)來來來 2026-03-21 | 4/200 |
|
|
|
[考研] 生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑 +4 | zxts12 2026-03-21 | 7/350 |
|
|
|
[基金申請(qǐng)] 山東省面上項(xiàng)目限額評(píng)審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
|
[考研] 材料求調(diào)劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|
|
[考研] 一志愿南大，0703化學(xué)，分?jǐn)?shù)336，求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|
|
|
[考研] 工科0856求調(diào)劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
|
[考研] 0856調(diào)劑，是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|