版塊導(dǎo)航: 正在加載中...

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

24小時(shí)熱門版塊排行榜

返回列表

exabyss916

新蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 177
散金: 65
帖子: 78
在線: 57小時(shí)
蟲號(hào): 507882
注冊(cè): 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左右）。
請(qǐng)問到底哪種方法是合理的呢？

回復(fù)此樓

» 猜你喜歡

求調(diào)劑，本科吉大，一志愿華工339 已經(jīng)有2人回復(fù)
意大利博洛尼亞大學(xué)招聘研究助理2名（2D材料/鋰離子與鈉離子電池方向）已經(jīng)有1人回復(fù)
物理化學(xué)論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有253人回復(fù)
關(guān)于碳量子點(diǎn)抗菌的研究已經(jīng)有15人回復(fù)
南大化學(xué)調(diào)劑已經(jīng)有11人回復(fù)
0703調(diào)劑已經(jīng)有2人回復(fù)
1 已經(jīng)有20人回復(fù)
0703 調(diào)劑已經(jīng)有5人回復(fù)
西南某211招27屆博士已經(jīng)有29人回復(fù)
大連工業(yè)杰青/長江團(tuán)隊(duì)-生物質(zhì)材料-儲(chǔ)能電池方向招收2026級(jí)博士生已經(jīng)有0人回復(fù)

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

使用Multiwfn繪制原子軌道圖形、研究原子殼層結(jié)構(gòu)及相對(duì)論效應(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
注冊(cè): 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
注冊(cè): 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一下，我的建議是你不要只比較能量，看看兩者對(duì)其它性質(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 (初中生)
金幣: 13001.7
散金: 27307
紅花: 300
帖子: 4672
在線: 926.9小時(shí)
蟲號(hào): 408372
注冊(cè): 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
注冊(cè): 2008-02-20
專業(yè): 理論和計(jì)算化學(xué)

引用回帖:

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

確實(shí)不太一樣，軌道能級(jí)也不一樣：
基組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 (初中生)
金幣: 13001.7
散金: 27307
紅花: 300
帖子: 4672
在線: 926.9小時(shí)
蟲號(hào): 408372
注冊(cè): 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
注冊(cè): 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的回帖

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

返回列表

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

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

[考研] 316求調(diào)劑 +5	梁茜雯 2026-03-19	5/250	2026-03-20 21:26 by 你好你好666
[考研] 265求調(diào)劑 +12	梁梁校校 2026-03-19	13/650	2026-03-20 21:01 by 無際的草原
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +3	困困困困坤坤 2026-03-20	3/150	2026-03-20 20:38 by 學(xué)員8dgXkO
[考研] 289求調(diào)劑 +6	懷瑾握瑜l 2026-03-20	6/300	2026-03-20 20:30 by 學(xué)員8dgXkO
[考研] 0703化學(xué)調(diào)劑，六級(jí)已過，有科研經(jīng)歷 +13	曦熙兮 2026-03-15	13/650	2026-03-20 19:35 by Dream007008
[考研] 271材料工程求調(diào)劑 +7	.6lL 2026-03-18	7/350	2026-03-20 09:10 by xingguangj
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4	xbxudjdn 2026-03-18	4/200	2026-03-20 09:07 by 不168
[論文投稿] 申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒變，給編輯又寫郵件了，沒回復(fù) 10+3	wangf9518 2026-03-17	4/200	2026-03-19 23:55 by babero
[考研] 一志愿福大288有機(jī)化學(xué)，求調(diào)劑 +3	小木蟲200408204 2026-03-18	3/150	2026-03-19 13:31 by houyaoxu
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 一志愿武理材料305分求調(diào)劑 +5	想上岸的鯉魚 2026-03-18	6/300	2026-03-18 17:53 by 無際的草原
[考研] 297求調(diào)劑 +8	戲精丹丹丹 2026-03-17	8/400	2026-03-18 14:30 by laoshidan
[考研] 材料專碩306英一數(shù)二 +10	z1z2z3879 2026-03-16	13/650	2026-03-18 14:20 by 007_lilei
[考研] 302求調(diào)劑 +10	呼呼呼。。。。 2026-03-17	10/500	2026-03-18 12:45 by Linda Hu
[考研] 303求調(diào)劑 +4	睿08 2026-03-17	6/300	2026-03-18 11:01 by Iveryant
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 304求調(diào)劑 +3	曼殊2266 2026-03-14	3/150	2026-03-16 16:39 by houyaoxu
[考研] 288求調(diào)劑 +4	奇點(diǎn)0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky
[考研] 289求調(diào)劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper

24小時(shí)熱門版塊排行榜

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

» 猜你喜歡

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

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

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