计算化学公社

 找回密码 Forget password
 注册 Register
Views: 7589|回复 Reply: 3
打印 Print 上一主题 Last thread 下一主题 Next thread

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

[复制链接 Copy URL]

308

帖子

3

威望

6251

eV
积分
6619

Level 6 (一方通行)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 Warm_Cloud 于 2019-2-27 14:39 编辑

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




评分 Rate

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

查看全部评分 View all ratings

欢迎使用量子化学软件Amesp

529

帖子

1

威望

5495

eV
积分
6044

Level 6 (一方通行)

2#
发表于 Post on 2019-2-28 09:05:39 | 只看该作者 Only view this author
暖大师的CPHF那里可否再详细一些,你实际编程的时候应该用了z-vector吧?方程组是用的什么方法解的?

308

帖子

3

威望

6251

eV
积分
6619

Level 6 (一方通行)

3#
 楼主 Author| 发表于 Post on 2019-2-28 10:02:30 | 只看该作者 Only view this author
niobium 发表于 2019-2-28 09:05
暖大师的CPHF那里可否再详细一些,你实际编程的时候应该用了z-vector吧?方程组是用的什么方法解的?

对,z-vector。迭代求解的,给一个Pia的初值,比如Lia/(ei-ea) ,然后迭代求解。
欢迎使用量子化学软件Amesp

308

帖子

3

威望

6251

eV
积分
6619

Level 6 (一方通行)

4#
 楼主 Author| 发表于 Post on 2019-2-28 10:06:14 | 只看该作者 Only view this author
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, 2025-8-14 16:38 , Processed in 1.313722 second(s), 26 queries , Gzip On.

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