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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

brafid

鐵桿木蟲 (正式寫手)

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

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

這是已經(jīng)優(yōu)化好的結(jié)構(gòu)，我加上個電子算單點能，基組稍微大了一些，算了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

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

用的基函數(shù)也多了點
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.

計算結(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個循環(huán)了，難道單純增加循環(huán)數(shù)到1000么？

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

回復此樓

» 猜你喜歡

中國科學院成都山地災害與環(huán)境研究所2026年特別研究助理（博士后）招聘啟事已經(jīng)有0人回復
全獎博士英國利物浦大學 × 臺灣清華大學聯(lián)合培養(yǎng) 已經(jīng)有0人回復
物理化學論文潤色/翻譯怎么收費? 已經(jīng)有177人回復
尋博導已經(jīng)有13人回復
基金正文30頁指的是報告正文還是整個申請書已經(jīng)有9人回復
哈爾濱工業(yè)大學（深圳）綠色化學團隊招收2026年秋季入學博士生已經(jīng)有11人回復
澳大利亞昆士蘭科技大學PhD博士獎學金已經(jīng)有0人回復

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

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

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

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

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

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

回帖支持 ( 顯示支持度最高的前 50 名 )

yypcnbe

木蟲 (小有名氣)

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

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
brafid: 金幣+38, ★★★★★最佳答案, 感謝幫助，基本都能算下來了，剩余金幣奉送 2014-05-11 16:43:57

引用回帖:

7樓: Originally posted by brafid at 2014-04-15 17:05:40
那要是算絕熱電離能呢？我也正在算絕熱的電離能作比較，也是好幾天了沒跑下來呢...

你這個計算如果只要一個單點能，就可以不用重新算了，都已經(jīng)收斂到10-7了，一般設置SCF收斂為6就可以了，除非對于特別小的分子的精確計算。

贊一下(1人)

回復此樓

8樓2014-04-15 20:47:23

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

普通回帖

youyno

銀蟲 (正式寫手)

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

【答案】應助回帖

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

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

贊一下

回復此樓

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

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

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

槍下游魂

木蟲 (著名寫手)

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

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

贊一下

回復此樓

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

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

brafid

鐵桿木蟲 (正式寫手)

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

引用回帖:

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

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

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

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

贊一下

回復此樓

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

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

brafid

鐵桿木蟲 (正式寫手)

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

引用回帖:

3樓: Originally posted by 槍下游魂 at 2014-04-11 19:29:11
LS說的在理，而且加了一個電子可能結(jié)構(gòu)也會變化，最后的結(jié)果可能還是要先跑優(yōu)化

跑優(yōu)化的是絕熱親和能與電離能，不跑優(yōu)化的是垂直電離能，看的兩篇做抗氧化性的文獻上用的垂直，不知道哪個用的多？

贊一下

回復此樓

5樓2014-04-11 23:31:34

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

yypcnbe

木蟲 (小有名氣)

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

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
brafid: 金幣+10, ★★★很有幫助, 先10個金幣致謝，如果有效剩余金幣附贈 2014-04-15 17:04:47

你SCF設置為scfcon=5或者6就行了，默認的是8，對于垂直電離能，收斂到5已經(jīng)足夠了。

贊一下

回復此樓

6樓2014-04-15 09:58:00

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

brafid

鐵桿木蟲 (正式寫手)

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

引用回帖:

6樓: Originally posted by yypcnbe at 2014-04-15 09:58:00
你SCF設置為scfcon=5或者6就行了，默認的是8，對于垂直電離能，收斂到5已經(jīng)足夠了。

那要是算絕熱電離能呢？我也正在算絕熱的電離能作比較，也是好幾天了沒跑下來呢

贊一下

回復此樓

7樓2014-04-15 17:05:40

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

brafid

鐵桿木蟲 (正式寫手)

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

引用回帖:

8樓: Originally posted by yypcnbe at 2014-04-15 20:47:23
你這個計算如果只要一個單點能，就可以不用重新算了，都已經(jīng)收斂到10-7了，一般設置SCF收斂為6就可以了，除非對于特別小的分子的精確計算。...

能否再咨詢下，我如果是用稍微低一點的基組方法算這個陰離子體系的絕熱電離勢，比如6-311+G來優(yōu)化陰離子結(jié)構(gòu)，然后用大機組6-311++G（3df，2p）來算單點能，是不是結(jié)果一樣，而且會快些？

贊一下

回復此樓

9樓2014-05-11 16:48:36

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

yypcnbe

木蟲 (小有名氣)

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

引用回帖:

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

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

贊一下

回復此樓

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

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

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

返回列表

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

[考研] 311求調(diào)劑 +4	冬十三 2026-03-18	4/200	2026-03-18 21:47 by 盡舜堯1
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5	share_joy 2026-03-16	5/250	2026-03-18 14:48 by haxia
[考研] 材料與化工一志愿南昌大學327求調(diào)劑推薦 +8	Ncdx123456 2026-03-13	9/450	2026-03-18 14:40 by haxia
[考研] 331求調(diào)劑（0703有機化學 +7	ZY-05 2026-03-13	8/400	2026-03-18 14:13 by 007_lilei
[考研] 303求調(diào)劑 +4	睿08 2026-03-17	6/300	2026-03-18 11:01 by Iveryant
[考研] 材料，紡織，生物（0856、0710），化學招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考研] 268求調(diào)劑 +6	簡單點0 2026-03-17	6/300	2026-03-18 09:04 by 無際的草原
[考研] 278求調(diào)劑 +5	煙火先于春 2026-03-17	5/250	2026-03-18 08:43 by 星空星月
[考研] 296求調(diào)劑 +5	大口吃飯身體健 2026-03-13	5/250	2026-03-17 21:05 by 不惑可樂
[考研] 梁成偉老師課題組歡迎你的加入 +8	一鴨鴨喲 2026-03-14	10/500	2026-03-17 15:07 by 一鴨鴨喲
[考研] 290求調(diào)劑 +6	孔志浩 2026-03-12	11/550	2026-03-17 14:41 by 周舟舟77
[考研] 211本，11408一志愿中科院277分，曾在中科院自動化所實習 +6	Losir 2026-03-12	7/350	2026-03-17 12:09 by danranxie
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 26考研一志愿中國石油大學(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 289求調(diào)劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper
[考研] 304求調(diào)劑 +6	Mochaaaa 2026-03-12	7/350	2026-03-13 22:18 by 星空星月
[考研] 求調(diào)劑 +5	一定有學上- 2026-03-12	5/250	2026-03-13 18:31 by ms629
[考研] 材料301分求調(diào)劑 +5	Liyouyumairs 2026-03-12	5/250	2026-03-13 14:42 by JourneyLucky
[考研] 289求調(diào)劑 +3	李政瑩 2026-03-12	3/150	2026-03-13 11:02 by 求調(diào)劑zz

24小時熱門版塊排行榜

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

» 猜你喜歡

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

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

【答案】應助回帖

【答案】應助回帖

【答案】應助回帖

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

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