请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 580|回复: 3

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

[复制链接]

212

帖子

3

威望

2099

eV
积分
2371

Level 5 (御坂)

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

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

评分

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

查看全部评分

186

帖子

1

威望

1858

eV
积分
2064

Level 5 (御坂)

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

212

帖子

3

威望

2099

eV
积分
2371

Level 5 (御坂)

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

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

212

帖子

3

威望

2099

eV
积分
2371

Level 5 (御坂)

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

我这个文档里面没有推导,完全是用来编程的。等有时间了应该会补充一下推导。
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2019-8-21 09:46 , Processed in 0.159700 second(s), 27 queries .

快速回复 返回顶部 返回列表