24小时热门版块排行榜    

查看: 337  |  回复: 2
当前主题已经存档。

wangyujia

木虫 (正式写手)

[交流] 【求助】DMOL算出来的这些东东是原子轨道吗?

头一次用DMOL计算分子轨道,算了一下水分子,倒是把它的五个轨道给算出来了,但是在结果文件里还有这样一些东东,我觉得是原子轨道,但是又不太像,请高手指教。
Hydrogen     nbas=  1, z=   1, nrfn=  3, rcut=  6.05, e_ref= -0.037940 Ha
   n=1  L=0  occ= 1.00 e=      -0.215807Ha        -5.8724eV
   n=1  L=0  occ= 0.00 e=      -0.844563Ha       -22.9817eV
   n=2  L=1  occ= 0.00 e=      -1.999921Ha       -54.4206eV  eliminated
Oxygen       nbas=  2, z=   8, nrfn=  7, rcut=  6.05, e_ref= -0.058400 Ha
   n=1  L=0  occ= 2.00 e=     -18.743105Ha      -510.0261eV
   n=2  L=0  occ= 2.00 e=      -0.860853Ha       -23.4250eV
   n=2  L=1  occ= 4.00 e=      -0.327086Ha        -8.9005eV
   n=2  L=0  occ= 0.00 e=      -2.129958Ha       -57.9591eV
   n=2  L=1  occ= 0.00 e=      -1.593333Ha       -43.3568eV
   n=3  L=2  occ= 0.00 e=      -2.722097Ha       -74.0721eV
   n=3  L=2  occ= 0.00 e=      -1.384180Ha       -37.6655eV  eliminated

就拿H来说,n=1及L=0时不就一个轨道吗?怎么给出了两行啊?这两行的能量还相差这么大??
再请高手顺便解释一下nrfn,rcut,e_ref各是什么意思?前两个的单位各是什么?谢谢!!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

study163

木虫 (正式写手)

天生我材

★ ★ ★ ★ ★
wangyujia(金币+4,VIP+0):谢谢啊!还望继续解决我的问题啊!! 12-5 09:21
aylayl08(金币+1,VIP+0):谢谢参与 12-5 09:47
你选 的基组是极化的吧,所以你的氢考虑到了P轨道(functions with angular momentum one higher than that of highest occupied orbital in free atom. )
dnp :  Use DN basis with polarization functions, i.e., functions with angular momentum one higher than that of highest occupied orbital in free atom.
dnd: dnd Double-numerical + d-DNP basis except that no p functions are used on hydrogen.
我的计算我做主
2楼2009-12-04 10:21:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyujia

木虫 (正式写手)

引用回帖:
Originally posted by study163 at 2009-12-4 10:21:
你选 的基组是极化的吧,所以你的氢考虑到了P轨道(functions with angular momentum one higher than that of highest occupied orbital in free atom. )
dnp :  Use DN basis with polarization functions, i ...

谢谢!!
在你的提示下,我把说明书中的basis set部分看了一遍,我选的是DND,所以p轨道那一行的后面还有eliminated的字样。
说明书中说DN是在MIN的基础上加了a second set of valence AOs,所以H有1s和1s‘,DND是在DN的基础上加了d-function。加上的这两个东东我都不知道是什么东东,大概知道的人不会太多,毕竟没几个人做计算方法的。我也不去管它算了!!

另外再问一个问题,我算了一下H原子,DMOL给出的能量有总能,还有一个binding energy,总能大约就是H原子的基态能量,-13.6eV,那个binding energy是什么啊?怎么那么小啊?才0.05eV!!和以前在XPS里提到了那个binding energy不太一样啊!

[ Last edited by wangyujia on 2009-12-5 at 09:20 ]
3楼2009-12-04 11:22:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wangyujia 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 12/600 2026-08-01 20:23 by Kddn1pB6pAzX
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 9zw1gfyi7x 2026-07-31 10/500 2026-08-01 20:11 by Kddn1pB6pAzX
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 9zw1gfyi7x 2026-07-31 9/450 2026-08-01 20:10 by Kddn1pB6pAzX
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 jzg3xk1po8 2026-07-31 13/650 2026-08-01 20:02 by Kddn1pB6pAzX
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 jzg3xk1po8 2026-07-31 10/500 2026-08-01 19:57 by Kddn1pB6pAzX
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 jzg3xk1po8 2026-07-31 10/500 2026-08-01 19:53 by Kddn1pB6pAzX
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 9zw1gfyi7x 2026-08-01 7/350 2026-08-01 18:42 by Kddn1pB6pAzX
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 9zw1gfyi7x 2026-08-01 6/300 2026-08-01 18:37 by Kddn1pB6pAzX
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 9zw1gfyi7x 2026-08-01 6/300 2026-08-01 18:37 by Kddn1pB6pAzX
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 9zw1gfyi7x 2026-08-01 8/400 2026-08-01 18:31 by Kddn1pB6pAzX
[考博] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 9zw1gfyi7x 2026-08-01 8/400 2026-08-01 18:31 by Kddn1pB6pAzX
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 9zw1gfyi7x 2026-07-31 11/550 2026-08-01 18:22 by Kddn1pB6pAzX
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 9zw1gfyi7x 2026-07-31 10/500 2026-08-01 18:20 by Kddn1pB6pAzX
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 9zw1gfyi7x 2026-07-31 11/550 2026-08-01 18:17 by Kddn1pB6pAzX
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jzg3xk1po8 2026-07-31 7/350 2026-08-01 17:42 by Kddn1pB6pAzX
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 8/400 2026-07-31 18:13 by beefly
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[基金申请] 今年的WR进展到哪一步了? +5 wsgjhwz 2026-07-25 9/450 2026-07-30 12:59 by 可淡不可忘
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[分析] 水蒸气蒸馏法乳化现象 +3 月光娇娇 2026-07-26 3/150 2026-07-26 23:25 by caiyun
信息提示
请填处理意见