ghcacj

[交流] 【求助】用Orca計算一個片簇的單點能，不收斂。

用Orca計算一個片簇的單點能，怎么算都不收斂，之前也算過其他的片簇，一般20次迭代就可以收斂了，這次計算的這個算了125次都不收斂，dE發(fā)生振蕩。
用的是ubuntu 10.04+orca 2.8
計算輸入文件如下：

! ECP{LANLDZ=Zn}
! UKS B3LYP/G 6-31G TightSCF Grid4 SP DFT
! PAL4
* xyz 0 1
C    14.7569    6.7818    3.1858
H    14.913    6.0981    3.7958
N    15.4572    7.9572    3.1056
C    14.7569    10.2092    13.8052
H    14.913    10.8929    13.1952
N    15.4572    9.0338    13.8854
C    13.8052    2.2341    6.7818
H    13.1952    2.0780    6.0981
N    13.8854    1.5338    7.9572
C    13.8052    14.7569    10.2092
H    13.1952    14.913    10.8929
N    13.8854    15.4572    9.0338
C    15.2773    11.6813    6.2614
H    14.5936    12.2913    6.4175
N    16.4527    11.6011    6.9617
C    15.2773    5.3097    10.7296
H    14.5936    4.6997    10.5735
N    16.4527    5.3899    10.0293
C    15.2773    6.2614    11.6813
H    14.5936    6.4175    12.2913
N    16.4527    6.9617    11.6011
H    15.6351    9.1972    9.7579
H    14.9045    10.0961    10.8131
H    15.2256    10.6636    9.3875
C    15.2773    10.7296    5.3097
H    14.5936    10.5735    4.6997
N    16.4527    10.0293    5.3899
H    15.6351    7.7938    7.2331
H    14.9045    6.8949    6.1779
H    15.2256    6.3274    7.6035
C    14.7569    3.1858    6.7818
H    14.913    3.7958    6.0981
H    16.2893    1.2624    9.8514
H    15.3904    2.3176    10.582
H    14.8229    0.8920    10.2609
N    15.4572    3.1056    7.9572
C    14.7569    13.8052    10.2092
H    14.913    13.1952    10.8929
H    16.2893    15.7286    7.1396
H    15.3904    14.6734    6.4090
H    14.8229    16.099    6.7301
N    15.4572    13.8854    9.0338
C    13.8052    6.7818    2.2341
H    13.1952    6.0981    2.0780
H    15.7286    9.8514    0.7017
H    14.6734    10.582    1.6006
H    16.099    10.2609    2.1681
N    13.8854    7.9572    1.5338
C    13.8052    10.2092    14.7569
H    13.1952    10.8929    14.913
H    15.7286    7.1396    16.2893
H    14.6734    6.4090    15.3904
H    16.099    6.7301    14.8229
N    13.8854    9.0338    15.4572
C    14.902    8.6297    2.0890
C    15.3955    9.9482    1.5955
C    14.902    8.3613    14.902
C    15.3955    7.0428    15.3955
C    14.902    2.0890    8.6297
C    15.3955    1.5955    9.9482
C    14.902    14.902    8.3613
C    15.3955    15.3955    7.0428
C    16.8568    10.5845    10.5845
C    15.5383    10.091    10.091
C    16.8568    6.4065    6.4065
C    15.5383    6.9000    6.9000
C    19.2251    10.2092    3.1858
H    19.069    10.8929    3.7958
N    18.5248    9.0338    3.1056
C    19.2251    6.7818    13.8052
H    19.069    6.0981    13.1952
N    18.5248    7.9572    13.8854
C    20.1768    14.7569    6.7818
H    20.7868    14.913    6.0981
N    20.0966    15.4572    7.9572
C    20.1768    2.2341    10.2092
H    20.7868    2.0780    10.8929
N    20.0966    1.5338    9.0338
C    18.7047    11.6813    10.7296
H    19.3884    12.2913    10.5735
N    17.5293    11.6011    10.0293
C    18.7047    5.3097    6.2614
H    19.3884    4.6997    6.4175
N    17.5293    5.3899    6.9617
H    18.3469    7.7938    9.7579
H    19.0775    6.8949    10.8131
H    18.7564    6.3274    9.3875
C    18.7047    10.7296    11.6813
H    19.3884    10.5735    12.2913
N    17.5293    10.0293    11.6011
H    18.3469    9.1972    7.2331
H    19.0775    10.0961    6.1779
H    18.7564    10.6636    7.6035
C    18.7047    6.2614    5.3097
H    19.3884    6.4175    4.6997
N    17.5293    6.9617    5.3899
C    19.2251    3.1858    10.2092
H    19.069    3.7958    10.8929
H    17.6927    1.2624    7.1396
H    18.5916    2.3176    6.4090
H    19.1591    0.8920    6.7301
N    18.5248    3.1056    9.0338
C    19.2251    13.8052    6.7818
H    19.069    13.1952    6.0981
H    17.6927    15.7286    9.8514
H    18.5916    14.6734    10.582
H    19.1591    16.099    10.2609
N    18.5248    13.8854    7.9572
C    20.1768    6.7818    14.7569
H    20.7868    6.0981    14.913
H    18.2534    9.8514    16.2893
H    19.3086    10.582    15.3904
H    17.883    10.2609    14.8229
N    20.0966    7.9572    15.4572
C    20.1768    10.2092    2.2341
H    20.7868    10.8929    2.0780
H    18.2534    7.1396    0.7017
H    19.3086    6.4090    1.6006
H    17.883    6.7301    2.1681
N    20.0966    9.0338    1.5338
C    19.08       8.3613    2.0890
C    18.5865    7.0428    1.5955
C    19.08       8.6297    14.902
C    18.5865    9.9482    15.3955
C    19.08    14.902    8.6297
C    18.5865    15.3955    9.9482
C    19.08       2.0890    8.3613
C    18.5865    1.5955    7.0428
C    17.1252    10.5845    6.4065
C    18.4437    10.091    6.9000
C    17.1252    6.4065    10.5845
C    18.4437    6.9000    10.091
Zn    16.991    12.7432    8.4955
Zn    16.991    4.2477    8.4955
Zn    16.991    8.4955    12.7432
Zn    16.991    8.4955    4.2477
H    13.3084    16.2196    8.7409
H    20.6736    16.2196    8.2501
H    13.3084    8.7409    16.2196
H    20.6736    8.7409    0.7714
H    13.3084    8.2501    0.7714
H    20.6736    0.7714    8.7409
H    20.6736    8.2501    16.2196
H    13.3084    0.7714    8.2501
*
片簇結(jié)構(gòu)如下：

[ Last edited by ghcacj on 2011-3-2 at 16:05 ]

回復(fù)此樓

» 猜你喜歡

西南某211招27屆博士已經(jīng)有33人回復(fù)
大連工業(yè)杰青/長江團隊-生物質(zhì)材料-儲能電池方向招收2026級博士生已經(jīng)有0人回復(fù)
物理化學(xué)論文潤色/翻譯怎么收費? 已經(jīng)有204人回復(fù)
深圳大學(xué)高等研究院趙偉課題組招收化學(xué)調(diào)劑生已經(jīng)有3人回復(fù)
青島大學(xué)固態(tài)電池研究團隊招收材料、化學(xué)、物理等方向博士研究生已經(jīng)有32人回復(fù)
化學(xué)工程085602 305分求調(diào)劑已經(jīng)有17人回復(fù)
26博士申請已經(jīng)有0人回復(fù)
武漢紡織大學(xué)全國重點實驗室楊應(yīng)奎教授團隊擬招17-19位碩博生！已經(jīng)有0人回復(fù)
材料與化工專碩接受跨專業(yè)調(diào)劑（重慶市屬高校）已經(jīng)有0人回復(fù)
武漢紡織大學(xué)全國重點實驗室楊應(yīng)奎教授團隊擬招17-19位碩博生！已經(jīng)有0人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

如何處理自旋計算產(chǎn)生的能帶已經(jīng)有3人回復(fù)
學(xué)位學(xué)歷證書和成績單的公證都是2010年12月20日的，現(xiàn)在簽證能用嗎？已經(jīng)有6人回復(fù)
請教：orca 中用CCSD(T)計算單點能出錯--提示內(nèi)存不足該怎么辦？已經(jīng)有10人回復(fù)
大家說說管理集群時的奇聞怪事吧已經(jīng)有22人回復(fù)
難道每次計算前都咬進行收斂性測試嗎？已經(jīng)有14人回復(fù)
收斂精度已經(jīng)有7人回復(fù)
計算時候總是不能收斂如何解決？已經(jīng)有5人回復(fù)
MOLCAS中采用CASPT2計算CO2二價陽離子1a2態(tài)不收斂怎么辦已經(jīng)有4人回復(fù)
vasp計算：優(yōu)化能收斂，但同樣的設(shè)置靜態(tài)計算不收斂已經(jīng)有3人回復(fù)
【求助】Orca 2.8的一個計算贗勢的例子報錯已經(jīng)有12人回復(fù)
【ghcacj個人文集】單機環(huán)境下安裝openmpi使用Orca并行運算的方法介紹已經(jīng)有12人回復(fù)
【求助】請教使用orca計算出的這個自旋密度圖是怎么畫的？已經(jīng)有5人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2011-03-02 16:01:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

recoli

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 2547.2
帖子: 553
在線: 92.9小時
蟲號: 253003

ghcacj(金幣+3): 謝謝 2011-03-03 13:59:08

據(jù)我所知orca默認是開啟RI近似的。但是6-31G這種基組沒有輔助基組，所以應(yīng)該關(guān)閉RI近似
%method
RI off
end

贊一下

回復(fù)此樓

3樓2011-03-02 17:29:36

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 8 個回答

ghcacj

榮譽版主 (著名寫手)

我想請教下各位，單點能計算不能收斂，有哪些常用的解決辦法？比如Opt？還是調(diào)整片簇的切割方法？還是增大基組？

贊一下

回復(fù)此樓

2樓2011-03-02 16:39:11

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ghcacj

榮譽版主 (著名寫手)

引用回帖:

Originally posted by recoli at 2011-03-02 17:29:36:
據(jù)我所知orca默認是開啟RI近似的。但是6-31G這種基組沒有輔助基組，所以應(yīng)該關(guān)閉RI近似
%method
RI off
end

按照你的設(shè)置，基組也改為6-31G*，但是還是不行，有沒有更好的建議？

贊一下

回復(fù)此樓

4樓2011-03-03 13:59:41

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

recoli

金蟲 (正式寫手)

應(yīng)助: 1 (幼兒園)
金幣: 2547.2
帖子: 553
在線: 92.9小時
蟲號: 253003

ghcacj(金幣+5): 謝謝 2011-03-03 20:14:38

我看了一下你的結(jié)構(gòu)，Zn應(yīng)該都是Zn2+離子吧？中間作為橋的四個配體每個少一個H是-1價，那么總體電荷應(yīng)該是+4？樓主認為呢？

電荷和多重度是很重要的，一定要設(shè)置正確，不然即使收斂了得到的也不是正確結(jié)果。

贊一下

回復(fù)此樓

5樓2011-03-03 17:50:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +4	零八# 2026-03-28	4/200	2026-03-28 15:01 by 17865157980
[考研] 求調(diào)劑 +6	爭取九點睡 2026-03-28	6/300	2026-03-28 14:14 by Ainin_
[考研] 311求調(diào)劑 +9	lin0039 2026-03-26	9/450	2026-03-28 13:05 by 唐沐兒
[考研] 085602 307分求調(diào)劑 +7	不知道叫什么！ 2026-03-26	7/350	2026-03-28 09:57 by 神馬都不懂
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +9	丹青奶蓋 2026-03-26	10/500	2026-03-28 07:45 by barnett0632
[考研] 295求調(diào)劑 +5	1428151015 2026-03-27	6/300	2026-03-28 04:04 by fmesaito
[考研] 0856材料化工調(diào)劑總分330 +10	zhubinhao 2026-03-27	10/500	2026-03-28 03:34 by fmesaito
[考研] 287求調(diào)劑 +10	land xuxu 2026-03-26	10/500	2026-03-27 15:33 by 帕爾馬拉特
[考研] 348求調(diào)劑 +4	小懶蟲不懶了 2026-03-27	5/250	2026-03-27 12:47 by 果果媽咪
[考研] 材料求調(diào)劑 +5	.m.. 2026-03-25	5/250	2026-03-27 11:08 by 不吃魚的貓
[考研] 求調(diào)劑323材料與化工 +7	1124361 2026-03-24	7/350	2026-03-27 10:22 by wangjy2002
[考研] 325求調(diào)劑 +5	李嘉圖·S·路 2026-03-23	5/250	2026-03-27 00:42 by wxiongid
[考研] 一志愿天津大學(xué)339材料與化工求調(diào)劑 +3	江往賣魚 2026-03-26	3/150	2026-03-26 09:42 by 王小欠i
[考研] 求調(diào)劑 +3	李李不服輸 2026-03-25	3/150	2026-03-25 13:03 by cmz0325
[考研] 調(diào)劑 +4	13853210211 2026-03-24	4/200	2026-03-24 19:44 by ms629
[考研] 305分求調(diào)劑（食品工程） +5	Sxy112 2026-03-21	7/350	2026-03-24 12:27 by 544594351
[考研] 341求調(diào)劑(一志愿湖南大學(xué)070300) +5	番茄頭--- 2026-03-22	6/300	2026-03-23 23:45 by Txy@872106
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長期招生公告 +3	李子目 2026-03-22	3/150	2026-03-22 21:01 by 怎么釋懷
[考研] 315分，誠求調(diào)劑，材料與化工085600 +3	13756423260 2026-03-22	3/150	2026-03-22 20:11 by edmund7
[考研] 275求調(diào)劑 +6	shansx 2026-03-22	8/400	2026-03-22 15:27 by barlinike