版塊導航: 正在加載中...

登錄注冊

應《網(wǎng)絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

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 ]

回復此樓

» 猜你喜歡

cif求助下載已經(jīng)有1人回復
高效微納米儲能材料與產(chǎn)業(yè)化推進實驗室蘇州招收碩士研究生（含調(diào)劑）已經(jīng)有5人回復
物理化學論文潤色/翻譯怎么收費? 已經(jīng)有52人回復
高效微納米儲能材料與產(chǎn)業(yè)化推進實驗室蘇州招收碩士研究生（含調(diào)劑）已經(jīng)有0人回復
高效微納米儲能材料與產(chǎn)業(yè)化推進實驗室蘇州招收碩士研究生（含調(diào)劑）已經(jīng)有2人回復
化學化工學院招收化工、化學、材料等相關(guān)方向研究生（學碩、專碩都有調(diào)劑名額）已經(jīng)有1人回復
高效微納米儲能材料與產(chǎn)業(yè)化推進實驗室蘇州招收碩士研究生（含調(diào)劑）已經(jīng)有0人回復
溫州大學新能源材料微觀溯源課題組招收2026級博士研究生（材料、化學方向）已經(jīng)有6人回復
材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流已經(jīng)有5人回復
26申博已經(jīng)有2人回復

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

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

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

查看全部散金貼

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

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

ghcacj

榮譽版主 (著名寫手)

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

贊一下

回復此樓

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

已閱回復此樓關(guān)注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這種基組沒有輔助基組，所以應該關(guān)閉RI近似
%method
RI off
end

贊一下

回復此樓

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

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

ghcacj

榮譽版主 (著名寫手)

引用回帖:

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

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

贊一下

回復此樓

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

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

recoli

金蟲 (正式寫手)

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

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

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

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

贊一下

回復此樓

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

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

ghcacj

榮譽版主 (著名寫手)

引用回帖:

Originally posted by recoli at 2011-03-03 17:50:42:
我看了一下你的結(jié)構(gòu)，Zn應該都是Zn2+離子吧？中間作為橋的四個配體每個少一個H是-1價，那么總體電荷應該是+4？樓主認為呢？

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

麻煩幫我看看呢，量子化學方面我是入門選手。
我的結(jié)構(gòu)是一個Zn作為節(jié)點，連接四個甲基咪唑，形成一種類似于沸石中Si-O西面體的結(jié)果，為了計算電荷，按照文獻的方法，截取出這樣的一個片段進行量化計算，因為是截取的片段，因此終端需要掛載一定的原子來補充，之前碰到的結(jié)果都是末端原本連接一個C原子，所以填上一個甲基，現(xiàn)在末端應該連接Zn，末端采用H作為補充不知道是否合適，但是看到文獻里有這么做，沒想到怎么計算都不收斂，按理說總體電荷應該為0，多重度為1，可是就是不行，正在嘗試是否是開殼層導致的，有的人建議將同樣的結(jié)構(gòu)先用閉殼層計算，然后拿這個結(jié)果作為intial guess進行開殼層計算，但是結(jié)果還沒出來。祈禱了。
謝謝你的關(guān)注。

[ Last edited by ghcacj on 2011-3-3 at 20:21 ]

贊一下

回復此樓

6樓2011-03-03 20:14:25

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

recoli

金蟲 (正式寫手)

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

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

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

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

贊一下

回復此樓

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

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

ghcacj

榮譽版主 (著名寫手)

引用回帖:

Originally posted by recoli at 2011-03-04 00:33:15:
就你的結(jié)構(gòu)來說，我認為電荷是4，閉殼層，所以輸入文件寫*xyz 4 1。

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

不得不說，你的判斷非常正確，認真分析了我要計算的結(jié)構(gòu)之后，修改計算文件，真的很快就收斂了，13次循環(huán)，非常感謝。

贊一下

回復此樓

8樓2011-03-04 10:13:17

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 ghcacj 的主題更新

返回列表

普通表情龍兔虎貓高級回復 (可上傳附件)

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

[考博] 環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名 +3	QGZDSYS 2026-03-13	3/150	2026-03-18 00:13 by PY-Hzb
[考研] 301求調(diào)劑 +6	yy要上岸呀 2026-03-17	6/300	2026-03-17 23:58 by 星空星月
[考研] 296求調(diào)劑 +5	大口吃飯身體健 2026-03-13	5/250	2026-03-17 21:05 by 不惑可樂
[考研] 299求調(diào)劑 +4	△小透明* 2026-03-17	4/200	2026-03-17 20:09 by peike
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 0854可跨調(diào)劑，一作一項核心論文五項專利，省、國級證書40+數(shù)一英一287 +3	小李0854 2026-03-16	3/150	2026-03-17 13:40 by 熱情沙漠
[考研] 211本，11408一志愿中科院277分，曾在中科院自動化所實習 +6	Losir 2026-03-12	7/350	2026-03-17 12:09 by danranxie
[考研] 271求調(diào)劑 +12	生如夏花… 2026-03-11	14/700	2026-03-17 10:56 by lovewei0727
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 286求調(diào)劑 +3	lemonzzn 2026-03-16	5/250	2026-03-16 20:43 by lemonzzn
[考研] 一志愿985，本科211，0817化學工程與技術(shù)319求調(diào)劑 +5	Liwangman 2026-03-15	5/250	2026-03-16 17:10 by 我的船我的海
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 085600調(diào)劑 +5	漾漾123sun 2026-03-12	6/300	2026-03-16 15:58 by 漾漾123sun
[考研] 277材料科學與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運氣yunqi
[考研] 0856求調(diào)劑 +3	劉夢微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流 +4	自由的_飛翔 2026-03-13	5/250	2026-03-14 19:39 by hmn_wj
[考研] 308 085701 四六級已過求調(diào)劑 +7	溫喬喬喬喬 2026-03-12	14/700	2026-03-14 10:49 by JourneyLucky
[考研] 0703，333分求調(diào)劑一志愿鄭州大學-物理化學 +3	李魔女斗篷 2026-03-11	3/150	2026-03-13 22:24 by JourneyLucky
[考研] 301求調(diào)劑 +6	Liyouyumairs 2026-03-11	6/300	2026-03-13 20:11 by JourneyLucky
[考研] 290求調(diào)劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike

24小時熱門版塊排行榜

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

» 猜你喜歡

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

» 搶金幣啦！回帖就可以得到: 查看全部散金貼

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

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

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

查看全部散金貼