24小时热门版块排行榜    

查看: 11505  |  回复: 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 ]
回复此楼

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

第一性原理 第一性原理和电化学 计算MS等 计算化学
量子化学和计算化学 mulwiwfn 理论计算 Gaussian
微观分子理论和模拟

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

技术博客:http://hi.baidu.com/coolrainbow/blog
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cationly

木虫 (正式写手)

很赞!
2楼2011-02-27 17:49:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

recoli

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
膜拜一下,码了这么多字和代码分享给大家,实在是我辈楷模
P.S. 志玲姐姐真给力
6楼2011-02-27 20:14:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongma2008

木虫 (著名写手)

论文?苦!主。

膜拜一下。。。。
7楼2011-02-27 23:16:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuanan027

木虫 (著名写手)

强悍!
Timeandtidewaitfornoman.Comeon,baby!
8楼2011-02-28 00:22:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xulisonghai

铜虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
感谢!膜拜!但不知这么笨的自己能不能学会呢?
9楼2011-02-28 08:55:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hbinjm

木虫 (著名写手)

小柔的爸


小木虫(金币+0.5):给个红包,谢谢回帖交流
很强大
虽然现在转行了 但还是要顶
10楼2011-02-28 09:14:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 coolrainbow 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 重要消息,中午系统在维护 +4 yuleib84 2026-08-18 5/250 2026-08-18 14:25 by 学员tGYWa5
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +45 医学老男孩 2026-08-13 101/5050 2026-08-18 13:40 by 医学老男孩
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 KZyU1THERRVo 2026-08-18 4/200 2026-08-18 13:28 by vHWVV9D5Deb5
[基金申请] 今天维护系统维护 祝所有人 高中 +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
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 5/250 2026-08-18 11:25 by 無關想念
[基金申请] 时间戳变了,能看出什么问题? +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
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +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
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
信息提示
请填处理意见