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
*
片簇結構如下：

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

回復此樓

» 猜你喜歡

求助立方晶體In(OH)3的CIF卡片已經有2人回復
[電子教材]天津大學物理化學（第七版）上、下冊已經有4人回復
物理化學論文潤色/翻譯怎么收費? 已經有72人回復
陰離子交換膜電解二氧化碳還原已經有0人回復
九江學院2026年最新高層次人才招聘公告已經有0人回復
深圳信息職業(yè)技術學院2026年教師招聘公告（最新）已經有0人回復
研究生招生已經有0人回復
山西大學資源與環(huán)境工程研究所化工分離/環(huán)境功能材料小組招收2026年博士生已經有3人回復
一志愿南師大0703化學 275求調劑已經有4人回復
A 區(qū)過線就可申！重慶雙非公辦高校 0856調劑，接受跨專業(yè)！納米材料 / 潤滑熱門方向已經有0人回復

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

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

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

查看全部散金貼

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

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

recoli

金蟲 (正式寫手)

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

ghcacj(金幣+12): 謝謝 2011-03-04 10:13:27

就你的結構來說，我認為電荷是4，閉殼層，所以輸入文件寫*xyz 4 1。

把C換成甲基是可以的，但是在N上添加H會影響整體的電荷。

贊一下

回復此樓

7樓2011-03-04 00:33:15

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

查看全部 8 個回答

ghcacj

榮譽版主 (著名寫手)

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

贊一下

回復此樓

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

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

recoli

金蟲 (正式寫手)

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

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

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

贊一下

回復此樓

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

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

ghcacj

榮譽版主 (著名寫手)

引用回帖:

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

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

贊一下

回復此樓

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

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

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

[考研] 一志愿北化085600材料專碩275\|有文章專利｜求調劑 +12	Micky11223 2026-03-25	13/650	2026-03-31 09:36 by longlotian
[考研] 0703化學 +17	妮妮ninicgb 2026-03-27	17/850	2026-03-31 09:36 by 有機小民工
[考研] 22408 359分調劑 +4	Qshers 2026-03-27	8/400	2026-03-31 08:53 by Qshers
[考研] 085602化工求調劑（331分） +8	111@127 2026-03-30	8/400	2026-03-30 21:23 by 研究僧導導
[考研] 277跪求調劑 +8	1915668 2026-03-27	12/600	2026-03-30 21:01 by dophin1985
[考研] 327求調劑 +5	小卡不卡. 2026-03-29	5/250	2026-03-30 19:30 by Wang200018
[考研] 085600，專業(yè)課化工原理，320分求調劑 +5	大饞小子 2026-03-29	5/250	2026-03-30 19:02 by 源_2020
[考研] 291求調劑 +8	HanBeiNingZC 2026-03-24	8/400	2026-03-30 17:01 by 無際的草原
[考研] 求調劑 +7	青春裁為三截 2026-03-29	7/350	2026-03-30 13:14 by laoshidan
[考研] 0856求調劑 +8	楒桉 2026-03-28	8/400	2026-03-30 10:00 by wzy-lxz
[考研] 一志愿南航 335分 \| 0856 \| GPA 4.07 \| 有科研經歷 +8	cccchenso 2026-03-29	8/400	2026-03-29 23:53 by 我是小康
[考研] 總分293求調劑 +8	加一一九 2026-03-25	11/550	2026-03-29 19:53 by 無際的草原
[考研] 340求調劑 +6	Amber00 2026-03-26	6/300	2026-03-29 12:06 by 無際的草原
[考研] 調劑考研 +3	王杰一 2026-03-29	3/150	2026-03-29 08:09 by fmesaito
[考研] 356求調劑 +3	gysy?s?a 2026-03-28	3/150	2026-03-29 00:33 by 544594351
[考研] 316求調劑 +7	江辭666 2026-03-26	7/350	2026-03-28 21:28 by sanrepian
[考研] 308求調劑 +7	墨墨漠 2026-03-27	7/350	2026-03-28 07:43 by 熱情沙漠
[考研] 314求調劑 +3	溪云珂 2026-03-26	3/150	2026-03-27 11:20 by sanrepian
[論文投稿] Journal of Mechanical Science and Technology +3	Russ_ss 2026-03-25	5/250	2026-03-27 10:49 by 陸小果畫大餅
[考研] 【2026考研調劑】制藥工程 284分求相關專業(yè)調劑名額 +4	袁奐奐 2026-03-25	8/400	2026-03-25 14:32 by lbsjt