當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖

brafid

應(yīng)助: 8 (幼兒園)
金幣: 6344
散金: 10
帖子: 440
在線: 245.5小時(shí)
蟲號: 122753
注冊: 2005-12-04
性別: GG
專業(yè): 先進(jìn)晶體材料

[求助] 已經(jīng)結(jié)構(gòu)優(yōu)化好，加了一個(gè)電子，算單點(diǎn)能，為啥就算不下來了？已有2人參與

這是已經(jīng)優(yōu)化好的結(jié)構(gòu)，我加上個(gè)電子算單點(diǎn)能，基組稍微大了一些，算了10天，居然告訴我不收斂，請教高手幫忙看看怎么能算下來
---------------------------------------------
# ub3lyp/6-311++g(3df,2p) sp scf=maxcycle=600
---------------------------------------------
1/38=1/1;
2/12=2,17=6,18=5,40=1/2;
3/5=4,6=6,7=1216,11=2,16=1,25=1,30=1,74=-5,116=2/1,2,3;
4//1;
5/5=2,7=600,38=5/2;
6/7=2,8=2,9=2,10=2,28=1/1;
99/5=1,9=1/99;
-------------------
Title Card Required
-------------------
Symbolic Z-matrix:
Charge = -1 Multiplicity = 2
O                   -1.42336 6.52227  -0.44806
H                   -2.27078 6.94039  -0.24277

當(dāng)然可能分子稍大了點(diǎn)
Stoichiometry C29H28O10(1-,2)
Framework group  C1[X(C29H28O10)]
Deg. of freedom 195

用的基函數(shù)也多了點(diǎn)
Two-electron integral symmetry is turned on.
  1801 basis functions,  2598 primitive gaussians,  2035 cartesian basis functions
142 alpha electrons    141 beta electrons
   nuclear repulsion energy    4152.3093305335 Hartrees.

計(jì)算結(jié)果如下：
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB3LYP) =  -1874.82223793    A.U. after  601 cycles
         Convg  = 0.8130D-07          -V/T =  2.0041
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 0.7611 S= 0.5055
<L.S>= 0.000000000000E+00
Annihilation of the first spin contaminant:
S**2 before annihilation    0.7611, after    0.7501
DiagDN has N=  1801 LTot=  5437 but NE2=    3 cannot use DSYEVD.
Convergence failure -- run terminated.
Error termination via Lnk1e in /share/apps/gauss/g09/l502.exe at Fri Apr  4 11:20:03 2014.
Job cpu time: 61 days 14 hours 54 minutes 33.5 seconds.

哪位高手給提題建議，看看怎么能算下來？已經(jīng)是600個(gè)循環(huán)了，難道單純增加循環(huán)數(shù)到1000么？

謝謝各位的不吝賜教�。�！

回復(fù)此樓

» 猜你喜歡

085404（計(jì)算機(jī)技術(shù)）293分希望調(diào)劑到材料相關(guān)專業(yè) 已經(jīng)有0人回復(fù)
汕頭大學(xué)化學(xué)招調(diào)劑，相關(guān)內(nèi)容如下。已經(jīng)有4人回復(fù)
物理化學(xué)論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有110人回復(fù)
大灣區(qū)大學(xué)（籌）劉天輝課題組招收博士生（與中山大學(xué)物理學(xué)院聯(lián)合培養(yǎng)）已經(jīng)有16人回復(fù)
PbS量子點(diǎn)如何確定濃度已經(jīng)有11人回復(fù)
大灣區(qū)大學(xué)劉天輝課題組招收2026級博士生（與中山大學(xué)物理學(xué)院聯(lián)合培養(yǎng)）已經(jīng)有22人回復(fù)
PbS量子點(diǎn)紫外吸收已經(jīng)有0人回復(fù)
比利時(shí)魯汶大學(xué)與國家留學(xué)基金委共同資助博士研究生CSC-KU Leuven PhD Scholarship 已經(jīng)有0人回復(fù)
半透性容器進(jìn)行高濕或者低濕研究已經(jīng)有0人回復(fù)
倫敦大學(xué)學(xué)院Benjamin Moss博士招收人工智能、原位光譜、催化相關(guān)26年博士生已經(jīng)有11人回復(fù)
中國科學(xué)院成都山地災(zāi)害與環(huán)境研究所2026年特別研究助理（博士后）招聘啟事已經(jīng)有0人回復(fù)

» 本主題相關(guān)商家推薦: (我也要在這里推廣)

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

CASTEP幾何結(jié)構(gòu)優(yōu)化的問題雖然收斂了感覺不是很對請大家?guī)兔φ乙幌略?/a> 已經(jīng)有11人回復(fù)

vasp結(jié)構(gòu)優(yōu)化出現(xiàn)問題：算完一個(gè)離子步后就停止已經(jīng)有12人回復(fù)
大俠，我用高斯對同一個(gè)結(jié)構(gòu)進(jìn)行優(yōu)化，用高精度算沒有虛頻，為什么半經(jīng)驗(yàn)AM1算有虛呀已經(jīng)有5人回復(fù)
castep計(jì)算單點(diǎn)能的幾何優(yōu)化問題已經(jīng)有12人回復(fù)
【求助】VASP中，“結(jié)構(gòu)優(yōu)化”和“電子自洽計(jì)算”的關(guān)系？已經(jīng)有8人回復(fù)
【求助】PM3優(yōu)化的基礎(chǔ)上采用6-31G*算單點(diǎn)能得到的結(jié)果與PM3下算得的大小次序不一致已經(jīng)有15人回復(fù)

1樓 2014-04-11 16:20:16

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

yypcnbe

木蟲 (小有名氣)

應(yīng)助: 36 (小學(xué)生)
金幣: 3500.2
紅花: 5
帖子: 198
在線: 2185.2小時(shí)
蟲號: 1245634
注冊: 2011-03-26
性別: GG
專業(yè): 催化化學(xué)

引用回帖:

9樓: Originally posted by brafid at 2014-05-11 16:48:36
能否再咨詢下，我如果是用稍微低一點(diǎn)的基組方法算這個(gè)陰離子體系的絕熱電離勢，比如6-311+G來優(yōu)化陰離子結(jié)構(gòu)，然后用大機(jī)組6-311++G（3df，2p）來算單點(diǎn)能，是不是結(jié)果一樣，而且會快些？...

你對于前面的中性體系用的什么基組去優(yōu)化的？如果你用的是6-311++G（3df，2p)基組優(yōu)化的中性體系的話，那么最好不要用你說的方法去做，就算你用的還是6-311+G基組優(yōu)化中性體系，你也沒必要用后面的那個(gè)大基組優(yōu)化陰離子了。還有如果你用的是6-311++G（3df，2p)基組優(yōu)化的中性體系的話，為了節(jié)省時(shí)間，提高效率，我建議你可以先用6-311+G*基組去優(yōu)化陰離子，然后再用那個(gè)大基組進(jìn)一步去優(yōu)化，基本上很快會收斂。

贊一下

回復(fù)此樓

10樓2014-05-11 17:37:23

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

查看全部 12 個(gè)回答

youyno

銀蟲 (正式寫手)

應(yīng)助: 114 (高中生)
金幣: 157.5
散金: 3774
紅花: 57
帖子: 666
在線: 217.1小時(shí)
蟲號: 1817680
注冊: 2012-05-15
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
gmy1990: 金幣+2 2014-04-11 20:51:21
brafid: 金幣+2, ★有幫助, 提出合理化建議 2014-04-11 22:18:02

為什么要加一個(gè)電子呢？加了一個(gè)電子之后自旋多重度也發(fā)生改變，體系的模式都不一樣，以之前的結(jié)構(gòu)算單點(diǎn)能就算成功了，又有什么作用呢？實(shí)在不知道lz想干什么，我們又怎么能夠幫到你呢？

贊一下

回復(fù)此樓

平生多磨礪，男兒自橫行！

2樓2014-04-11 17:10:31

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

槍下游魂

木蟲 (著名寫手)

應(yīng)助: 718 (博后)
金幣: 8336.9
散金: 175
紅花: 80
帖子: 1498
在線: 542.6小時(shí)
蟲號: 2620258
注冊: 2013-08-28
專業(yè): 理論和計(jì)算化學(xué)

LS說的在理，而且加了一個(gè)電子可能結(jié)構(gòu)也會變化，最后的結(jié)果可能還是要先跑優(yōu)化

贊一下

回復(fù)此樓

3樓2014-04-11 19:29:11

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

brafid

鐵桿木蟲 (正式寫手)

應(yīng)助: 8 (幼兒園)
金幣: 6344
散金: 10
帖子: 440
在線: 245.5小時(shí)
蟲號: 122753
注冊: 2005-12-04
性別: GG
專業(yè): 先進(jìn)晶體材料

引用回帖:

2樓: Originally posted by youyno at 2014-04-11 17:10:31
為什么要加一個(gè)電子呢？加了一個(gè)電子之后自旋多重度也發(fā)生改變，體系的模式都不一樣，以之前的結(jié)構(gòu)算單點(diǎn)能就算成功了，又有什么作用呢？實(shí)在不知道lz想干什么，我們又怎么能夠幫到你呢？

我是想算電子親和能和電離能的，加一個(gè)電子算親和能，減少一個(gè)電子，算電離能，結(jié)果都差不多，600圈之后都不轉(zhuǎn)了。

優(yōu)化好的結(jié)構(gòu)，直接加電子和減少電子，算的都是垂直電離勢和親和勢，而如果是再優(yōu)化了結(jié)構(gòu)的算的都是絕熱電離勢，我看文獻(xiàn)上這么說的，而且一些最新的文獻(xiàn)報(bào)道的這類分子的抗氧化性的計(jì)算，也都是算的垂直電離勢和親和勢。

我就是想得到一個(gè)垂直電離勢，其實(shí)就是算出失去一個(gè)電子的單點(diǎn)能和得到一個(gè)電子的單點(diǎn)能，不知道能否給出點(diǎn)合理的建議？

贊一下

回復(fù)此樓

4樓2014-04-11 23:19:35

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

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

[考研] 321求調(diào)劑 +6	何潤采123 2026-03-18	6/300	2026-03-19 10:55 by daliangogo
[考研] 一志愿天大材料與化工（085600）總分338 +5	蔡大美女 2026-03-13	5/250	2026-03-19 10:44 by 是小劉呀～
[考研] 材料工程專碩調(diào)劑 +5	204818@lcx 2026-03-17	6/300	2026-03-18 22:55 by 204818@lcx
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +3	xbxudjdn 2026-03-18	3/150	2026-03-18 22:14 by zhq0425
[考研] 材料專業(yè)求調(diào)劑 +5	hanamiko 2026-03-18	5/250	2026-03-18 20:19 by 星空星月
[考研] 266求調(diào)劑 +5	陽陽哇塞 2026-03-14	9/450	2026-03-18 15:05 by stone_128
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7	ZY-05 2026-03-13	8/400	2026-03-18 14:13 by 007_lilei
[考研] 0703化學(xué)調(diào)劑 +4	pupcoco 2026-03-17	7/350	2026-03-18 12:14 by djl2006
[考研] 0854，計(jì)算機(jī)類招收調(diào)劑 +3	胡辣湯放糖 2026-03-15	6/300	2026-03-18 12:09 by 上岸上岸……..
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考博] 26申博 +4	八6八68 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 278求調(diào)劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 294求調(diào)劑 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 中科大材料專碩319求調(diào)劑 +3	孟鑫材料 2026-03-13	3/150	2026-03-14 18:10 by houyaoxu
[基金申請] 現(xiàn)在如何回避去年的某一個(gè)專家，不知道名字 +3	zk200107 2026-03-12	6/300	2026-03-14 17:13 by zk200107
[考研] 復(fù)試調(diào)劑 +4	z1z2z3879 2026-03-14	5/250	2026-03-14 16:30 by JourneyLucky
[考研] 學(xué)碩285求調(diào)劑 +13	Wisjxn 2026-03-12	46/2300	2026-03-14 10:33 by JourneyLucky
[考研] 290求調(diào)劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike