| 7 | 1/1 | 返回列表 |
| 查看: 1222 | 回復(fù): 6 | ||||
exabyss916新蟲 (小有名氣)
|
[交流]
【求助】關(guān)于UGBS和ANO-RCC用在DKH計算的精度 已有3人參與
|
|
我用兩種方法分別計算了RbCs系統(tǒng),由于計算量太大只計算了單點能量,但既然方法都一樣只是Basis set不一樣即使單點能量應(yīng)該也可以直接比較。 兩種方法能量相差幾百個Hatree(計算的能量為10000個Hatree左右)。 請問到底哪種方法是合理的呢? |
新蟲 (小有名氣)
|
更多計算細(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)]\\@ |

專家顧問 (職業(yè)作家)
地溝油冶煉專家
![]() |
專家經(jīng)驗: +458 |

新蟲 (小有名氣)
|
確實不太一樣,軌道能級也不一樣: 基組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 |
專家顧問 (職業(yè)作家)
地溝油冶煉專家
![]() |
專家經(jīng)驗: +458 |
|
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é)果是錯的。 Gaussian使用一般收縮基組時效率非常低,例如Rb作ANO收縮和不收縮,UHF計算時間之比大約20:1,而molpro、molcas等程序為1:3。ANO-RCC改用非收縮形式,也許能得到和UGBS接近的結(jié)果。 還有個可能存在的問題。Cs ANO-RCC的s高斯函數(shù)指數(shù)遠(yuǎn)遠(yuǎn)大于UGBS,Gaussian在做積分時可能有問題。 |

新蟲 (小有名氣)
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工272求調(diào)劑 +9 | 阿斯蒂芬2004 2026-03-28 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 調(diào)劑求院校招收 +4 | 鶴鯨鴿 2026-03-28 | 4/200 |
|
|
[考研] 311求調(diào)劑 +9 | lin0039 2026-03-26 | 9/450 |
|
|
[考研] 數(shù)一英一271專碩(085401)求調(diào)劑,可跨 +4 | 前行必有光 2026-03-28 | 5/250 |
|
|
[考研] 277跪求調(diào)劑 +5 | 1915668 2026-03-27 | 9/450 |
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +9 | 多久上課 2026-03-27 | 9/450 |
|
|
[考研] 求調(diào)劑 +4 | 零八# 2026-03-27 | 4/200 |
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +8 | 北道巷 2026-03-24 | 8/400 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[碩博家園] 招收生物學(xué)/細(xì)胞生物學(xué)調(diào)劑 +3 | IceGuo 2026-03-26 | 4/200 |
|
|
[考研] 317求調(diào)劑 +7 | 蛋黃咸肉粽 2026-03-26 | 7/350 |
|
|
[考研] 325求調(diào)劑 +3 | Aoyijiang 2026-03-23 | 3/150 |
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +4 | 李李不服輸 2026-03-26 | 4/200 |
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|
|
[考研] 282求調(diào)劑 +3 | wcq131415 2026-03-24 | 3/150 |
|
|
[考研] 318求調(diào)劑 +5 | plum李子 2026-03-21 | 8/400 |
|
|
[考研] 340求調(diào)劑 +5 | 話梅糖111 2026-03-24 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑,求導(dǎo)師收 +7 | 天天好運來上岸?/a> 2026-03-24 | 7/350 |
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|