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ù)此樓

» 猜你喜歡

200年來人們一直無法破解的難題：氣液臨界點狀態(tài)方程已經(jīng)有7人回復(fù)
細(xì)胞培養(yǎng)，這22個細(xì)節(jié)一定要注意�。ㄒ唬� 已經(jīng)有0人回復(fù)
物理化學(xué)論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有76人回復(fù)
關(guān)于引力延遲下雙星系統(tǒng)角動量不守恒的疑惑---求高手解答！已經(jīng)有10人回復(fù)
已知面上掛了，小學(xué)校實驗方向無平臺，感覺太難了已經(jīng)有26人回復(fù)
CSC訪學(xué)博后項目獲批，外方學(xué)校暫停合作該怎么辦？已經(jīng)有51人回復(fù)
求助各位大佬，球球了，在這科薩頭已經(jīng)有0人回復(fù)
分子動力學(xué)模擬合作需求已經(jīng)有2人回復(fù)
招生材料與化工、環(huán)境科學(xué)調(diào)劑研究生已經(jīng)有4人回復(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ìn)行收斂性測試嗎？已經(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個人文集】單機(jī)環(huán)境下安裝openmpi使用Orca并行運(yùn)算的方法介紹已經(jīng)有12人回復(fù)
【求助】請教使用orca計算出的這個自旋密度圖是怎么畫的？已經(jīng)有5人回復(fù)

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

查看全部散金貼

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

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

ghcacj

榮譽(yù)版主 (著名寫手)

引用回帖:

Originally posted by recoli at 2011-03-02 17:29:36:
據(jù)我所知orca默認(rèn)是開啟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的回帖

查看全部 8 個回答

ghcacj

榮譽(yù)版主 (著名寫手)

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

贊一下

回復(fù)此樓

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

已閱回復(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默認(rèn)是開啟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的回帖

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？樓主認(rèn)為呢？

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

贊一下

回復(fù)此樓

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

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

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

[考研] 一志愿南開大學(xué)0710生物學(xué)359求調(diào)劑 +3	兔兔兔111223314 2026-03-29	3/150	2026-03-29 23:14 by Evan_Liu
[考研] 290求調(diào)劑 +3	dfffsar 2026-03-29	3/150	2026-03-29 22:38 by 毛毛毛阿莫2
[考研] 307求調(diào)劑 +9	超級伊昂大王 2026-03-24	10/500	2026-03-29 20:26 by 永字號
[考研] 調(diào)劑310 +12	溫柔的晚安 2026-03-25	13/650	2026-03-29 20:01 by 無際的草原
[考研] 343求調(diào)劑085601 +3	要努力學(xué)習(xí)x 2026-03-29	3/150	2026-03-29 18:35 by wxiongid
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +4	小張zxy張 2026-03-26	8/400	2026-03-28 19:25 by lbsjt
[考研] 291求調(diào)劑 +6	HanBeiNingZC 2026-03-24	6/300	2026-03-28 07:55 by baoball
[考研] 考研調(diào)劑 +4	Sanmu-124 2026-03-26	4/200	2026-03-27 17:49 by kiokin
[考研] 材料292調(diào)劑 +12	橘頌思美人 2026-03-23	12/600	2026-03-27 15:44 by caszguilin
[考研] 一志愿鄭大085600，310分求調(diào)劑 +5	李瀟可 2026-03-26	5/250	2026-03-27 11:14 by 不吃魚的貓
[考研] 材料求調(diào)劑 +5	.m.. 2026-03-25	5/250	2026-03-27 11:08 by 不吃魚的貓
[考研] 292求調(diào)劑 +4	求求了收下我吧?/a> 2026-03-26	4/200	2026-03-27 10:37 by zhshch
[考研] 材料調(diào)劑 5+4	想要一壺桃花水 2026-03-25	10/500	2026-03-26 19:56 by 不吃魚的貓
[考研] 289求調(diào)劑 +17	碩星赴 2026-03-23	17/850	2026-03-26 16:18 by 不吃魚的貓
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4	gdlf9999 2026-03-24	4/200	2026-03-25 23:01 by boxking200
[考研] 考研一志愿蘇州大學(xué)初始315（英一）求調(diào)劑 +3	sbdksD 2026-03-24	4/200	2026-03-25 18:16 by xcjcqu
[考研] 0854人工智能方向招收調(diào)劑 +4	章小魚567 2026-03-24	4/200	2026-03-25 13:29 by 2177681040
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3	開開心心沒煩惱 2026-03-23	4/200	2026-03-24 00:06 by 開開心心沒煩惱
[考研] 336化工調(diào)劑 +4	王大坦1 2026-03-23	5/250	2026-03-23 18:32 by allen-yin
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4	allen-yin 2026-03-23	6/300	2026-03-23 15:04 by 汪��？！