应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » 其他 » 关于HF方法和DFT方法中轨道能量的差别
121/2返回列表12下一页
查看: 3406  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

zhangguangping

木虫 (著名写手)

[交流] 关于HF方法和DFT方法中轨道能量的差别 已有5人参与

弱弱的问一下量化高手们：
在HF计算和DFT计算中，对于同一个体系的轨道能量的值差别很大，但是看到两者最后的单电子方程差别很小，形式上DFT就是多一个关联能，但是这个能量不至于差别这么大。
比如：分别用HF和DFT去计算一下H原子，会发现，HF得到的第一个占据轨道的能量就是-0.5 Hatree,而DFT只有-0.3 Hatree左右。但是两者最后得到的体系的能量几乎相同。现在就是说两者对于总能量的定义也不同。

所以想请教大家：他们的这个轨道能量的差别原因是什么？总能量的定义有什么区别？
回复此楼

» 收录本帖的淘帖专辑推荐

学习备用 VASP 量子化学扫盲贴

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
弘德明志博学笃行
1楼 2012-04-23 00:43:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

zhangguangping

木虫 (著名写手)
zhangguangping: 回帖置顶 2012-04-23 08:48:27
再具体说明一下问题：
我用高斯分别在HF水平和DFT水平计算了H原子的轨道和总能量。结果如下：
输入文件：
# hf/6-31g (or B3LYP/6-31g)

Title Card Required

0 2
H                 -0.59433960   -0.05660377    0.00000000


计算结果：
HF结果:

The electronic state is 2-A1G.
Alpha  occ. eigenvalues --   -0.49823
Alpha virt. eigenvalues --    0.91619
  Beta virt. eigenvalues --    0.09502   1.10807
...
1|1|UNPC-abc|SP|UHF|6-31G|H1(2)|ADMINISTRATOR|23-Apr-2012|0
||# hf/6-31g||Title Card Required||0,2|H,0,-0.5943396,-0.05660377,0.||
Version=IA32W-G09RevA.02|State=2-A1G|HF=-0.4982329|S2=0.75|S2-1=0.|S2A
=0.75|RMSD=0.000e+000|Dipole=0.,0.,0.|Quadrupole=0.,0.,0.,0.,0.,0.|PG=
OH [O(H1)]||@


DFT结果:

The electronic state is 2-A1G.
Alpha  occ. eigenvalues --   -0.31617
Alpha virt. eigenvalues --    0.71090
  Beta virt. eigenvalues --    0.06729   1.07729
...

1|1|UNPC-abc|SP|UB3LYP|6-31G|H1(2)|ADMINISTRATOR|23-Apr-201
2|0||# B3LYP/6-31g||Title Card Required||0,2|H,0,-0.5943396,-0.0566037
7,0.||Version=IA32W-G09RevA.02|State=2-A1G|HF=-0.5002728|S2=0.75|S2-1=
0.|S2A=0.75|RMSD=3.909e-012|Dipole=0.,0.,0.|Quadrupole=0.,0.,0.,0.,0.,
0.|PG=OH [O(H1)]||@

问题：DFT和HF算出来的轨道能量差别很大，而文件最后给出的总能量确几乎相同。用DFT比HF多考虑了一个关联能好像不能解释两者计算出来的轨道能量的差别吧？

@yumu0411 @lishijunzong
赞一下 回复此楼
弘德明志博学笃行
4楼2012-04-23 08:43:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

yumu0411

金虫 (著名写手)

打酱油
★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
gmy1990: 金币+1 2012-04-26 10:10:21
我感觉最主要的原因就是DFT多了关联函数了吧！这样DFT考虑的因素就比HF多！就会导致分子轨道能量又差别！
你说的总能是什么总能！内能？因为每个分子总能中多有HF总能
赞一下 回复此楼
学习学习再学习
2楼2012-04-23 07:20:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lishijunzong

铁杆木虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
一般而言，对于多电子体系，HF的能量会比DFT高，这是因为没有电子与电子考虑斥力的作用
赞一下 回复此楼
3楼2012-04-23 07:31:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lishijunzong

铁杆木虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
氢原子，只有一个电子，怎样来相关啊
赞一下 回复此楼
5楼2012-04-23 12:26:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangmt

至尊木虫 (著名写手)

我叫MT
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
顶楼上的。。。
H原子只有一个电子，哪来的相关？跟谁相关啊。。。。。
好歹你也换个体系来试验啊。。
赞一下 回复此楼
一群自以为正义凛然的年轻人将一切不能以科学解释的事情定性为封建迷信并大刀阔斧地进行消灭，其实这是修养不足学识浅薄的一种体现，也是可恶的偏执和愚蠢的自以
6楼2012-04-23 12:36:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yumu0411

金虫 (著名写手)

打酱油
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
空轨道和占有轨道只是相对的！难道一个氢的一个电子只在1s轨道！电子也有可能在空轨道存在吧？所以相关函数是有影响的！
赞一下 回复此楼
学习学习再学习
7楼2012-04-23 14:04:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
现在的问题是：不管对于一个H原子有没有相关能。为什么轨道的能量计算出来不同呢？(但是轨道的空间分布非常相似)
并且最后体系的能量是相同的。那么DFT中总能量的计算公式是什么？
我想知道导致这个差异的原因是什么。从单电子方程上看，两者的单电子方程是非常相似的。应该算出来的轨道能量也是相近的。
@beefly @yumu0411 @lishijunzong @zhangmt
赞一下 回复此楼
弘德明志博学笃行
8楼2012-04-23 15:14:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独唱团

金虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
dft算出来的轨道能量，没有实际意义，没法跟hf比的。
你用koopmans定理算算电离能，hf还将就，dft结果直接没法用。
赞一下 回复此楼
9楼2012-04-24 03:24:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
引用回帖:
9楼: Originally posted by 独唱团 at 2012-04-24 03:24:15:
dft算出来的轨道能量，没有实际意义，没法跟hf比的。
你用koopmans定理算算电离能，hf还将就，dft结果直接没法用。

谢谢。我也明白DFT算出来的轨道没有物理意义，轨道能量也不对应物理意义。我的意思是说，形式上很相近的两个方程，为什么解确差别这么大呢？主要问题是：波函数的空间分布差别不大，但是能量本征值差别很大。
赞一下 回复此楼
弘德明志博学笃行
10楼2012-04-24 09:20:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhangguangping 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026国自然函评费到账 +19 羊腰板 2026-08-21 22/1100 2026-08-25 10:38 by popular289
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 今天基金会出结果吗？20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 只有每年这种时候来逛逛小木虫 +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +7 scalable 2026-08-24 7/350 2026-08-24 23:01 by anata1209
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[基金申请] filecode，4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤？
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 科研孤儿太难了 +18 我4大白菜 2026-08-20 19/950 2026-08-24 09:47 by 凯恩广盛催化分�
[基金申请] 什么时候开奖？ +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 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 紫杉醇
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
关闭
请填处理意见
关闭