exabyss916

新蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 177
散金: 65
帖子: 78
在線: 57小時(shí)
蟲號(hào): 507882
注冊: 2008-02-20
專業(yè): 理論和計(jì)算化學(xué)

[交流] 【求助】關(guān)于UGBS和ANO-RCC用在DKH計(jì)算的精度已有3人參與

我用兩種方法分別計(jì)算了RbCs系統(tǒng)，由于計(jì)算量太大只計(jì)算了單點(diǎn)能量，但既然方法都一樣只是Basis set不一樣即使單點(diǎn)能量應(yīng)該也可以直接比較。

兩種方法能量相差幾百個(gè)Hatree（計(jì)算的能量為10000個(gè)Hatree左右）。
請問到底哪種方法是合理的呢？

回復(fù)此樓

» 猜你喜歡

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

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

使用Multiwfn繪制原子軌道圖形、研究原子殼層結(jié)構(gòu)及相對論效應(yīng)的影響已經(jīng)有17人回復(fù)

1樓 2010-05-17 12:08:31

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

exabyss916

新蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 177
散金: 65
帖子: 78
在線: 57小時(shí)
蟲號(hào): 507882
注冊: 2008-02-20
專業(yè): 理論和計(jì)算化學(xué)

更多計(jì)算細(xì)節(jié)：

基組1：
# MP2/gen Integral=DKH
這里gen產(chǎn)生的是ano-rcc基組，結(jié)果如下：
1\1\GINC-NODE04\SP\RMP2-FC\Gen\Cs1Rb1\name\08-May-2010\0\\# MP2/gen
Integral=DKH\\\0,1\Rb\Cs,1,4.3\\Version=EM64L-G09RevA.01\State=1-SG\HF
=-10488.3784426\MP2=-10488.3888615\RMSD=7.597e-09\PG=C*V [C*(Rb1Cs1)]\

基組2：
# MP2/UGBS guess(indo) Integral=DKH
結(jié)果如下：
1\1\GINC-NODE05\SP\RMP2-FC\UGBS\Cs1Rb1\name\10-May-2010\0\\# MP2/UGB
S guess(indo) Integral=DKH\\\0,1\Rb\Cs,1,4.3\\Version=EM64L-G09RevA.01
\State=1-SG\HF=-10759.4475783\MP2=-10759.4590548\RMSD=2.655e-09\PG=C*V [C*(Rb1Cs1)]\\@

贊一下

回復(fù)此樓

2樓2010-05-18 15:04:33

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

天空空

木蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 179.5
散金: 6
紅花: 4
帖子: 585
在線: 160.8小時(shí)
蟲號(hào): 602252
注冊: 2008-09-14
性別: MM
專業(yè): 數(shù)學(xué)物理

★ ★ ★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
aylayl08(金幣+3):謝謝指點(diǎn) 2010-05-18 16:29:46

一般而言，Primitive (uncontracted) basis set 是最精確的。 UGBS就是Primitive (uncontracted) basis set。但是你最好check一下，我的建議是你不要只比較能量，看看兩者對其它性質(zhì)的描述有區(qū)別沒？

贊一下(2人)

回復(fù)此樓

好好學(xué)習(xí)，天天向上！嘿嘿

3樓2010-05-18 16:01:41

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

beefly

專家顧問 (職業(yè)作家)

地溝油冶煉專家

專家經(jīng)驗(yàn): +458
QC強(qiáng)帖: 12
應(yīng)助: 88 (初中生)
金幣: 13000.7
散金: 27307
紅花: 300
帖子: 4672
在線: 926.9小時(shí)
蟲號(hào): 408372
注冊: 2007-06-21
性別: GG
專業(yè): 理論和計(jì)算化學(xué)
管轄: 計(jì)算模擬

★ ★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
aylayl08(金幣+2):謝謝提示 2010-05-19 15:10:27

不可能差這么大，即使和最小基組比較，最多只能差（1 hartree＊重原子數(shù)）。檢查一下HF軌道占據(jù)是否一致。還可能發(fā)生了基組依賴

贊一下(2人)

回復(fù)此樓

beefly《西太平洋大學(xué)現(xiàn)代英漢詞典》[bi:fli]牛肉一般地

4樓2010-05-18 21:33:09

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

exabyss916

新蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 177
散金: 65
帖子: 78
在線: 57小時(shí)
蟲號(hào): 507882
注冊: 2008-02-20
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

Originally posted by beefly at 2010-05-18 21:33:09:
不可能差這么大，即使和最小基組比較，最多只能差（1 hartree＊重原子數(shù)）。檢查一下HF軌道占據(jù)是否一致。還可能發(fā)生了基組依賴

確實(shí)不太一樣，軌道能級也不一樣：
基組1：
   Occupied  (SG) (SG) (SG) (SG) (PI) (PI) (SG) (SG) (SG) (PI)
               (PI) (SG) (PI) (PI) (SG) (PI) (PI) (DLTA) (DLTA)
               (SG) (SG) (SG) (PI) (PI) (DLTA) (DLTA) (PI) (PI)
               (SG) (SG) (PI) (PI) (DLTA) (DLTA) (PI) (PI) (SG)
               (SG) (SG) (PI) (PI) (SG) (SG) (PI) (PI) (SG)
Alpha  occ. eigenvalues -- **********-589.04971-212.16515-190.33495-190.33491
Alpha  occ. eigenvalues -- -190.33491 -89.74429 -45.77127 -38.89039 -38.87668
Alpha  occ. eigenvalues --  -38.87668 -38.54419 -38.54412 -38.54412 -27.83521
Alpha  occ. eigenvalues --  -27.83519 -27.83519 -27.83516 -27.83516 -15.37358
Alpha  occ. eigenvalues -- -9.34989  -7.61111  -7.61066  -7.61066  -7.60985
Alpha  occ. eigenvalues -- -7.60985  -6.93297  -6.93297  -6.93286  -5.70744
Alpha  occ. eigenvalues -- -5.70351  -5.70351  -3.27885  -3.27885  -3.27865
Alpha  occ. eigenvalues -- -3.27865  -3.27861  -2.11817  -1.33481  -0.70938
Alpha  occ. eigenvalues -- -0.70938  -0.70614  -0.49512  -0.49075  -0.49075
Alpha  occ. eigenvalues -- -0.15553

基組2：
   Occupied  (SG) (SG) (SG) (SG) (PI) (PI) (SG) (SG) (PI) (PI)
               (SG) (SG) (PI) (PI) (SG) (PI) (PI) (DLTA) (DLTA)
               (SG) (SG) (PI) (PI) (SG) (SG) (PI) (PI) (SG) (PI)
               (PI) (DLTA) (DLTA) (SG) (PI) (PI) (DLTA) (DLTA)
               (SG) (SG) (SG) (PI) (PI) (SG) (PI) (PI) (SG)
Alpha  occ. eigenvalues -- **********-561.64098-212.13438-190.30259-190.30255
Alpha  occ. eigenvalues -- -190.30255 -77.44271 -68.46967 -68.46962 -68.46962
Alpha  occ. eigenvalues --  -45.73965 -38.51214 -38.51201 -38.51201 -27.80330
Alpha  occ. eigenvalues --  -27.80327 -27.80327 -27.80318 -27.80318 -12.54071
Alpha  occ. eigenvalues -- -9.56621  -9.56605  -9.56605  -9.31841  -6.90121
Alpha  occ. eigenvalues -- -6.90103  -6.90103  -4.62803  -4.62798  -4.62798
Alpha  occ. eigenvalues -- -4.62785  -4.62785  -3.24712  -3.24705  -3.24705
Alpha  occ. eigenvalues -- -3.24690  -3.24690  -1.55366  -1.30469  -0.79716
Alpha  occ. eigenvalues -- -0.79562  -0.79562  -0.67816  -0.67732  -0.67732
Alpha  occ. eigenvalues -- -0.11975

贊一下

回復(fù)此樓

5樓2010-05-19 08:44:17

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

beefly

專家顧問 (職業(yè)作家)

地溝油冶煉專家

專家經(jīng)驗(yàn): +458
QC強(qiáng)帖: 12
應(yīng)助: 88 (初中生)
金幣: 13000.7
散金: 27307
紅花: 300
帖子: 4672
在線: 926.9小時(shí)
蟲號(hào): 408372
注冊: 2007-06-21
性別: GG
專業(yè): 理論和計(jì)算化學(xué)
管轄: 計(jì)算模擬

★ ★ ★ ★ ★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
aylayl08(金幣+5):謝謝提示 2010-05-19 15:10:37

UHF/UGBS能量
Rb: -2978.86122825
Cs: -7780.59443740
總和：-10759.45566565，和你算的RHF能量差不多。

ANO-RCC主頁上的CASPT2能量：
Rb: -2979.02
Cs: -7780.95

如果忽略相關(guān)能，應(yīng)該和UGBS更接近。所以你的ANO-RCC結(jié)果是錯(cuò)的。

Gaussian使用一般收縮基組時(shí)效率非常低，例如Rb作ANO收縮和不收縮，UHF計(jì)算時(shí)間之比大約20:1，而molpro、molcas等程序?yàn)?:3。ANO-RCC改用非收縮形式，也許能得到和UGBS接近的結(jié)果。

還有個(gè)可能存在的問題。Cs ANO-RCC的s高斯函數(shù)指數(shù)遠(yuǎn)遠(yuǎn)大于UGBS，Gaussian在做積分時(shí)可能有問題。

贊一下(2人)

回復(fù)此樓

beefly《西太平洋大學(xué)現(xiàn)代英漢詞典》[bi:fli]牛肉一般地

6樓2010-05-19 13:42:44

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

exabyss916

新蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 177
散金: 65
帖子: 78
在線: 57小時(shí)
蟲號(hào): 507882
注冊: 2008-02-20
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

Originally posted by beefly at 2010-05-19 13:42:44:
UHF/UGBS能量
Rb: -2978.86122825
Cs: -7780.59443740
總和：-10759.45566565，和你算的RHF能量差不多。

ANO-RCC主頁上的CASPT2能量：
Rb: -2979.02
Cs: -7780.95

如果忽略相關(guān)能，應(yīng)該和UGBS更接近。 ...

非常感謝詳細(xì)的回復(fù)和相關(guān)數(shù)據(jù)，我也問了gaussian技術(shù)支持，他們也說用應(yīng)該用非收縮基組計(jì)算，那我以后就都用UGBS系列計(jì)算。

贊一下

回復(fù)此樓

7樓2010-05-19 14:30:51

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

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

[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +21	rare12345 2026-03-18	21/1050	2026-03-20 14:31 by 無懈可擊111
[考研] 一志愿中南化學(xué)337求調(diào)劑 +4	niko- 2026-03-19	5/250	2026-03-20 13:49 by 促天成
[考研] 296求調(diào)劑 +4	www_q 2026-03-18	8/400	2026-03-20 12:12 by 學(xué)員8dgXkO
[考博] 招收博士1-2人 +3	QGZDSYS 2026-03-18	3/150	2026-03-20 11:58 by 呱呱呱呱叫
[考研] 一志愿西南交通專碩材料355 本科雙非求調(diào)劑 +4	西南交通專材355 2026-03-19	4/200	2026-03-20 11:39 by 花開富貴幸福人?/a>
[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3	晨昏線與星海 2026-03-19	3/150	2026-03-19 13:36 by houyaoxu
[考研] 0703化學(xué) 305求調(diào)劑 +4	FY_yy 2026-03-14	4/200	2026-03-19 05:54 by anny19840123
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +4	守候夕陽CF 2026-03-18	4/200	2026-03-18 22:16 by li123456789.
[考研] 08工科 320總分求調(diào)劑 +5	梨花珞晚風(fēng) 2026-03-17	5/250	2026-03-18 14:49 by haxia
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國級證書40+數(shù)一英一287 +8	小李0854 2026-03-16	8/400	2026-03-18 14:35 by 搏擊518
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 材料與化工專碩調(diào)劑 +5	heming3743 2026-03-16	5/250	2026-03-17 14:03 by 勇敢太監(jiān)王公公
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 中科大材料與化工319求調(diào)劑 +3	孟鑫材料 2026-03-14	3/150	2026-03-14 20:10 by ms629
[考研] 297一志愿上交085600求調(diào)劑 +5	指尖八千里 2026-03-14	5/250	2026-03-14 17:26 by a不易