24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 10000  |  回复: 65
本帖产生 1 个 QC强帖 ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

coolrainbow

木虫 (著名写手)

未来国家冻凉

[交流] 【coolrainbow个人文集】从手算理解量子化学:H2体系的量子化学计算(pdf+程序) 已有64人参与

本文以最小基组的氢气分子为例,以简单的理论分析和编程计算的方法,手动计算了H2分子的HF,FCI,MP2,DFT和优化,频率分析,TDDFT等,以帮助大家深入的理解量子化学的基本概念,特别是通过对Hartree-Fock波函数的远程行为与稳定性做了定量分析,详细的阐述了H2解离行为的物理意义。对TDDFT/TDHF也做了手动的计算,可以使大家对Gaussian中的TDDFT有更进一步的认识。本人编了一个小程序,可以复现本文中的全部手动计算结果,这个程序的源码我也提供了,你可以根据你的需要修改代码,以便让它输出你所需要的东西。同时,我在文章中也提供了Gaussian的输入文件,如果你跟着我的文章计算的话,会发现你算出的结果和Gaussian中是一样的。

这个工作是我春节在家期间的无聊之作,由于看林志玲姐姐受到了刺激而决定动笔写的,从写pdf到调试程序大概花了五六天的时间。我自认为理论分析写的还是不错的,不过程序代码写的很水,是专门针对H2的体系写的,并且TDDFT的代码我没有写进去。

如果大家觉得我写的东西还算有点价值的话,我就非常开心了,不白费我打了那么多的字。如果哪位能指出其中的错误,那我就更开心了。如果大家觉得有用的话,我将来还可以继续写下去,把本文中没有写的内容,如耦合簇计算,BSSE校正原理,经典轨迹动力学,相对论计算的加进去,或者某些大家觉得不太好的地方将来重写。

鸣谢

1本文的DFT计算的格点积分的技术参考了sobereva的文章:
http://hi.baidu.com/sobereva/blo ... 16fd56ccbf1ab8.html
感谢sobereva对知识无私分享的精神。
2 DFT计算的某些技术使用了CCL上的源代码:
http://www.ccl.net/cca/software/ ... bedev-Laikov-Grids/
对作者表示感谢。

版权

1 网络上的东西大家可以自由传播,不过做人要厚道,看在我打了那么多字的份上。大家在注明原作者(也就是我:coolrainbow)的条件下可以自由传播本pdf。
2 本文附带的程序代码中,除LBWT.f以外,全部代码可以自由传播与修改。不过说真的,这些代码写的真的不怎么样!
3 本文附带的程序代码中,LBWT.f的版权归原作者所有,具体版权可以参考:
http://www.ccl.net/cca/software/ ... bedev-Laikov-Grids/


附件中,我上传了完整的pdf和程序文件及源代码(2000多行我的代码+6000行LBWT.f),下面是目录和部分内容的截图


目录
§ 1 体系构建
§ 2 分子积分计算
§ 3 单组态RHF 计算
§ 4 验证… …
§ 5 RHF 的物理解释
   § 5.1 电子态
   § 5.2 分子轨道
   § 5.3 单、双电子密度分布
   § 5.4 波函数的进一步分析
   § 5.5 分子轨道基下能量表达式
   § 5.6 RHF 波函数的解离极限
§ 6 组态相互作用计算
§ 7 还不相信的话… …
§ 8 FCI 波函数的物理解释
   § 8.1 解离极限
   § 8.2 成分分析
   § 8.3 相关效应
   § 8.4 H2 的布居分析
§ 9 H2 的绝热态曲线;不相交原理
§ 10 精确能量的微扰展开:MP2 能量
§ 12 耦合簇方法:一句话
§ 13 势能面:优化和频率分析
    § 13.1 RHF 优化和频率分析
    § 13.2 一般的优化和频率分析
§ 14 RHF 的稳定性
§ 15 基态DFT 计算:Xα 方法
§ 16 激发态计算
    § 16.1 TDHF,TDDFT 计算
    § 16.2 验证
    § 16.3 直接计算激发
§ 17 程序









[ Last edited by 375642546 on 2011-4-19 at 15:37 ]
回复此楼
技术博客:http://hi.baidu.com/coolrainbow/blog
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qianyang

铁杆木虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
,LZ很强大。这个要顶。感谢!
39楼2011-04-15 21:44:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 66 个回答

qzhaosdu

金虫 (著名写手)

强烈支持
Anewday,anewhour,anewminute,anewpeople.
3楼2011-02-27 18:07:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hustyh0801

铁杆木虫 (著名写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
感谢楼主的无私奉献!!!
4楼2011-02-27 18:56:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ahehyh

木虫 (正式写手)

过客


小木虫(金币+0.5):给个红包,谢谢回帖交流
不错,可以学习下,呵呵!
我是天空的一片云
5楼2011-02-27 19:53:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 调剂 +4 GK72 2026-03-30 4/200 2026-03-30 20:32 by dick_runner
[考研] 一志愿郑州大学,080500学硕,总分317分求调剂 +10 举个栗子oi 2026-03-24 11/550 2026-03-30 19:31 by michael2011
[考研] 085601材料工程找调剂 +17 oatmealR 2026-03-29 18/900 2026-03-30 19:21 by Wang200018
[考研] 化学0703 调剂 306分 一志愿211 +8 26要上岸 2026-03-28 8/400 2026-03-30 16:39 by herarysara
[考研] 考研调剂 +5 Sanmu-124 2026-03-26 5/250 2026-03-30 13:36 by chemdavid
[考研] 求调剂,一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +6 @taotao 2026-03-26 7/350 2026-03-30 10:43 by 我是小康
[考研] 一志愿211,335分,0856,求调剂院校和导师 +7 倾____萧 2026-03-27 8/400 2026-03-30 09:37 by longlotian
[考研] 356求调剂 +3 gysy?s?a 2026-03-28 3/150 2026-03-29 00:33 by 544594351
[考研] 315求调剂 +4 akie... 2026-03-28 5/250 2026-03-28 21:05 by zhq0425
[考研] 本科新能源科学与工程,一志愿华理能动285求调剂 +3 AZMK 2026-03-27 5/250 2026-03-28 16:19 by xxxsssccc
[考研] 药学105500求调剂 +3 Ssun。。 2026-03-28 3/150 2026-03-28 11:24 by lxf170613
[考研] 085405 考的11408求各位老师带走 +3 Qiu学ing 2026-03-28 3/150 2026-03-28 09:19 by 乐呵呵的追梦人
[考研] 材料求调剂一志愿哈工大324 +7 闫旭东 2026-03-28 9/450 2026-03-28 08:51 by Xu de nuo
[考研] 274求调剂 +17 顾九笙要谦虚 2026-03-24 23/1150 2026-03-27 15:16 by caszguilin
[考研] 314求调剂 +3 溪云珂 2026-03-26 3/150 2026-03-27 11:20 by sanrepian
[考研] 一志愿吉大071010,316分求调剂 +3 xgbiknn 2026-03-27 3/150 2026-03-27 10:36 by guoweigw
[考研] 总分322求生物学/生化与分子/生物信息学相关调剂 +5 星沉uu 2026-03-26 6/300 2026-03-26 19:02 by macy2011
[考研] 一志愿天津大学339材料与化工求调剂 +3 江往卖鱼 2026-03-26 3/150 2026-03-26 09:42 by 王小欠i
[考研] 【2026考研调剂】制药工程 284分 求相关专业调剂名额 +4 袁奂奂 2026-03-25 8/400 2026-03-25 14:32 by lbsjt
[考研] 材料专硕331求调剂 +4 鲜当牛 2026-03-24 4/200 2026-03-24 15:58 by JourneyLucky
信息提示
请填处理意见