24小时热门版块排行榜    

查看: 11507  |  回复: 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的回帖

wang0912302

铁虫 (正式写手)

赞一个啊,辛苦了啊,,,非常感谢啊!


小木虫(金币+0.5):给个红包,谢谢回帖交流
赞一个啊,辛苦了啊,,,非常感谢啊!
41楼2011-04-16 08:01:50
已阅   回复此楼   关注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 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +46 医学老男孩 2026-08-13 103/5150 2026-08-18 15:33 by dcqxinyang
[基金申请] 明天放榜? +3 Shxjjxjkx 2026-08-18 3/150 2026-08-18 15:32 by a441914426
[教师之家] 两姐妹,一个辛辛苦苦30岁工作当大学教师,另一个23岁就工作当开车教练 +3 瞬息宇宙 2026-08-11 5/250 2026-08-18 14:55 by LNP@mRNA
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 8/400 2026-08-18 11:47 by gltch
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[基金申请] 时间戳变了,能看出什么问题? +5 基诺咪客 2026-08-17 6/300 2026-08-18 11:06 by hunter无悔
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
信息提示
请填处理意见