24小时热门版块排行榜    

查看: 4910  |  回复: 8

Yolanda7

捐助贵宾 (正式写手)


[交流] MS中各基组的含义

大家好。我在学习MS。先用dmol3计算了H2O分子的算例:DNP基组。我对于基组的设置不太明白,读不懂.outmol中的关于基组的部分
“Hydrogen     nbas=  1, z=   1, nrfn=  3, rcut=  6.99, e_ref= -0.043138 Ha  rcore=  0.00  zval=   1.00   1.00
n主量子数  L 轨道量子数 occ:占据程度   e:此轨道的本证能量
  n=1  L=0  occ= 1.00 e=      -0.232994Ha        -6.3401eV
  n=1  L=0  occ= 0.00 e=      -0.844913Ha       -22.9913eV
  n=2  L=1  occ= 0.00 e=      -1.999991Ha       -54.4226eV
Oxygen       nbas=  2, z=   8, nrfn=  7, rcut=  6.99, e_ref= -0.057357 Ha  rcore=  0.00  zval=   8.00   8.00
  n=1  L=0  occ= 2.00 e=     -18.907598Ha      -514.5021eV
  n=2  L=0  occ= 2.00 e=      -0.876373Ha       -23.8473eV
  n=2  L=1  occ= 4.00 e=      -0.330740Ha        -8.9999eV
  n=2  L=0  occ= 0.00 e=      -2.144799Ha       -58.3630eV
  n=2  L=1  occ= 0.00 e=      -1.593898Ha       -43.3722eV
  n=3  L=2  occ= 0.00 e=      -2.722210Ha       -74.0751eV
  n=3  L=2  occ= 0.00 e=      -1.387900Ha       -37.7667eV  eliminated

我想问下,nbas(?) e_ref(the calculated atomic reference energy)具体是什么意思,如何得到的?nrfn好像是给定原子的基函数的数目。本例中O的第七个函数怎么是eliminated?帮助文件中好像说: “The atomic program generates all atomic numeric basis sets, typically with several polarization functions. These functions can be used in a DMol3 calculation by using the extended or all keyword. You can select specific orbitals to be used in the calculation by using the user keyword.” 所以,程序一开始就产生所有的原子基组,包括极化函数;然后,使用时,根据DNP DN DND等类型的指定,确定哪些基组是有用的 或 eliminated?

Minimal、DNP、DN、DND基组具体如何理解?
我查了下资料:
MIN (Minimal)  极小基组,对每一个占据分子轨道用自由原子的一个原子轨道展开;
DN (Double Numeric)  双数字基组,对每一个占据分子轨道近似地用自由原子的两个原子轨道展开。我看到MS帮助文件中对O元素采用DN基组后,
  “Oxygen       nbas=  2, z=   8, nrfn=  7, rcut=  6.24, e_ref= -0.054 Ha
                                  rcore=  0.00  zval=   8.00   8.00
   n=1  L=0  occ= 2.00 e=     -18.758046Ha      -510.4326eV
   n=2  L=0  occ= 2.00 e=      -0.871142Ha       -23.7050eV
   n=2  L=1  occ= 4.00 e=      -0.338180Ha        -9.2023eV
   n=2  L=0  occ= 0.00 e=      -2.130079Ha       -57.9624eV
   n=2  L=1  occ= 0.00 e=      -1.593478Ha       -43.3608eV
   n=3  L=2  occ= 0.00 e=      -2.722222Ha       -74.0755eV  eliminated
   n=3  L=2  occ= 0.00 e=      -1.388886Ha       -37.7935eV  eliminated”
那么,自由原子的两个原子轨道好像是第二组价电子原子轨道出现两次吗。一个轨道怎么会出现两次,且有不同的本证能量?

DNP (Double Numeric with Polarization)  增加了极化函数的双数值基组(极化函数:函数中存在高于自由原子中的最高占据轨道角动量的角动量);这里极化函数,我不太明白。假如某原子最高占据轨道为1s,那么极化函数为2p,还是1p?某原子最高占据轨道为2s,那么极化函数为2p,还是3p?某原子最高占据轨道为2p,那么极化函数为2d,还是3d?某原子最高占据轨道为3d,那么极化函数为多少?最高占据轨道4f,极化函数是多少?

DND Double Numerical plus d-functions。除了对氢原子不使用 p轨道函数外,与 DNP基组相同?;又或者,DNP基组是在DND基组上在所有氢原子上附加了极化p函数。DND是对所有非氢原子加了d函数,那么,DNP中所有非氢原子的极化函数都是(主量子数+1)的d轨道吗?


[ Last edited by Yolanda7 on 2012-4-28 at 09:41 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Yolanda7

捐助贵宾 (正式写手)


★ ★ ★
fzx2008: 金币+3, 谢谢交流分享 2012-04-28 10:00:20
看了下扬州大学MS售后培训,对基组稍微理解:
• Min:一个占位原子轨道只使用一个数值轨道基组。速度更快,精度最低。
• DN:双数值轨道基组。在Min基组的基础上加入了第二组价电子原子轨道。   
• DND:双数值轨道基组+d轨道极化函数。与DN基组类似,但是为所有非氢原子加入了d轨道函数极化。为DMol3模块默认基组,在可接受的计算时内,确保精度在可接受的范围内。也是DMol3计算的最低可信精度基组。
• DNP:双数值轨道基组+p轨道极化函数。与DND基组类似,但是为所有氢原子加入了p轨道函数极化。 精度更高,计算时间更长。对于氢键计算更为重要。
• TNP:三数值轨道基组+轨道极化函数。 与DNP基组类似,为所有原子加入极化函数。目前只能够使用与H到Cl的元素,且不包括He和Ne。 精度更高,计算时间更长。
2楼2012-04-28 09:42:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liushilei1111

木虫 (著名写手)


不要鄙视,我是送祝福
3楼2012-04-28 09:58:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Yolanda7

捐助贵宾 (正式写手)


引用回帖:
3楼: Originally posted by liushilei1111 at 2012-04-28 09:58:16:
不要鄙视,我是送祝福

谢谢啦~
4楼2012-04-28 10:56:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Yolanda7

捐助贵宾 (正式写手)


但是,我后面的问题还是不很理解。有咩有高手啊?
5楼2012-04-29 08:29:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lq6865387

木虫 (著名写手)


6楼2012-04-29 23:00:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Yolanda7

捐助贵宾 (正式写手)


引用回帖:
6楼: Originally posted by lq6865387 at 2012-04-29 23:00:20:
顶起

谢谢
7楼2012-04-30 11:20:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lq6865387

木虫 (著名写手)


引用回帖:
7楼: Originally posted by Yolanda7 at 2012-04-30 11:20:33:
谢谢

呵呵 ~~~
8楼2012-05-01 15:26:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgygdx

新虫 (小有名气)


求普及物理意义
9楼2015-04-24 10:42:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Yolanda7 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +19 木水思豆 2026-08-25 22/1100 2026-08-26 06:24 by 玲珑133
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 Qj98d49dHFb9 2026-08-25 5/250 2026-08-26 04:57 by cNXvBfCpiZOM
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 7K1CJE38xLG4 2026-08-25 4/200 2026-08-26 02:04 by cNXvBfCpiZOM
[基金申请] 今天务委会开完了,明天出结果吗 +18 angus9576 2026-08-25 22/1100 2026-08-26 00:41 by merchancy
[基金申请] 有没有大神帮我看看基金代码 28+4 1234567wang 2026-08-24 10/500 2026-08-25 19:15 by lfy8008
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 19/950 2026-08-25 17:57 by chick875
[基金申请] 2026国自然函评费到账 +20 羊腰板 2026-08-21 23/1150 2026-08-25 16:34 by zsna
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 今天基金会出结果吗?20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 范进中举一文的中心思想 +7 炎黄贵胄 2026-08-22 8/400 2026-08-25 08:48 by ZJTJZ
[基金申请] 只有每年这种时候来逛逛小木虫 +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
信息提示
请填处理意见