计算化学公社

 找回密码 Forget password
 注册 Register
Views: 5394|回复 Reply: 3

[算法与编程] MP2闭壳层一阶导数公式

[复制链接 Copy URL]

268

帖子

3

威望

4329

eV
积分
4657

Level 6 (一方通行)

发表于 Post on 2019-2-20 16:48:19 | 显示全部楼层 Show all |阅读模式 Reading model
本帖最后由 Warm_Cloud 于 2019-2-27 14:39 编辑

半年前写了MP2的导数的代码,由于文献上的公式有些乱,当时踩了不少坑,现在整理了一下MP2闭壳层一阶导数公式,有兴趣的朋友可以照着写一下,其中用到的电子积分的一阶导数需要自己提供,如果里面有什么错误敬请大家指正。开壳层的和冻芯的情况,以后会继续更新。
MP2闭壳层一阶导数公式.pdf (88.45 KB, 下载次数 Times of downloads: 71)

评分 Rate

参与人数
Participants 7
eV +38 收起 理由
Reason
linqiaosong + 4 大师大法好
zyniso + 5 GJ!
bnulk + 5 大师大法好!
978142355 + 5 大师大法好!
sobereva + 10
让你变成回忆 + 4 好物!
lerel + 5 好物!

查看全部评分 View all ratings

459

帖子

1

威望

3999

eV
积分
4478

Level 6 (一方通行)

发表于 Post on 2019-2-28 09:05:39 | 显示全部楼层 Show all
暖大师的CPHF那里可否再详细一些,你实际编程的时候应该用了z-vector吧?方程组是用的什么方法解的?

268

帖子

3

威望

4329

eV
积分
4657

Level 6 (一方通行)

 楼主 Author| 发表于 Post on 2019-2-28 10:02:30 | 显示全部楼层 Show all
niobium 发表于 2019-2-28 09:05
暖大师的CPHF那里可否再详细一些,你实际编程的时候应该用了z-vector吧?方程组是用的什么方法解的?

对,z-vector。迭代求解的,给一个Pia的初值,比如Lia/(ei-ea) ,然后迭代求解。

268

帖子

3

威望

4329

eV
积分
4657

Level 6 (一方通行)

 楼主 Author| 发表于 Post on 2019-2-28 10:06:14 | 显示全部楼层 Show all
niobium 发表于 2019-2-28 09:05
暖大师的CPHF那里可否再详细一些,你实际编程的时候应该用了z-vector吧?方程组是用的什么方法解的?

我这个文档里面没有推导,完全是用来编程的。等有时间了应该会补充一下推导。

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2023-2-2 23:19 , Processed in 0.190736 second(s), 25 queries .

快速回复 返回顶部 返回列表 Return to list