计算化学公社

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

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

[复制链接 Copy URL]

311

帖子

3

威望

6507

eV
积分
6878

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: 99)

评分 Rate

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

查看全部评分 View all ratings

欢迎使用量子化学软件Amesp

538

帖子

1

威望

5917

eV
积分
6475

Level 6 (一方通行)

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

311

帖子

3

威望

6507

eV
积分
6878

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) ,然后迭代求解。
欢迎使用量子化学软件Amesp

311

帖子

3

威望

6507

eV
积分
6878

Level 6 (一方通行)

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

我这个文档里面没有推导,完全是用来编程的。等有时间了应该会补充一下推导。
欢迎使用量子化学软件Amesp

本版积分规则 Credits rule

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

GMT+8, 2026-8-8 03:24 , Processed in 0.212923 second(s), 25 queries , Gzip On.

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