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

计算化学公社

 找回密码
 现在注册!
查看: 1404|回复: 8

[Multiwfn资源与经验] 在Multiwfn中分析比CCSD更高级别波函数的方法

[复制链接]

1万

帖子

25

威望

1万

eV
积分
34206

管理员

公社社长

发表于 2017-11-28 00:13:26 | 显示全部楼层 |阅读模式
在Multiwfn中分析比CCSD更高级别波函数的方法

文/Sobereva @北京科音  2017-Nov-28


1 前言

Gaussian程序能产生的最高级别的波函数是CCSD,虽然这已经非常精确了,比如在Science, 355, 49 (2017)中作者以CCSD密度为金标准,考察了不同泛函对电子密度的重现性,但是有时候由于特殊原因,就是希望在更高级别下做波函数分析。例如在J. Chem. Theory Comput., 13, 4753 (2017)中,作者对CCSD密度还不知足,于是以CCSDTQ密度为参考考察了不同理论方法对电子密度及衍生量的计算精度;再比如,在J. Chem. Theory Comput., 13, 2705 (2017)中,作者提出了基于自然轨道占据数将动态和静态相关图形化展现的方法(此方法在Multiwfn中支持,看3.4.1及之后版本手册4.A.6节的例子),文中用到了FCI级别自然轨道,以使得对文中所考虑的电子结构复杂的情况,静态和动态相关都能完美展现。

能够产生比CCSD更高级别波函数的程序较少。PSI4可以产生最高到CCSD(T)级别波函数,MRCC可以产生无穷高阶耦合簇和无穷高阶CI的波函数(包括FCI。但是产生不了微扰近似版本的耦合簇的波函数)。Multiwfn程序从3.5(dev)版开始,能够结合这两个程序在这些比CCSD更高的理论级别做波函数分析,下面就说一下具体方法。本文使用的PSI4为1.1版,MRCC为Sept 25, 2017版。

对Multiwfn不了解者请参考《Multiwfn入门tips》(http://sobereva.com/167)、《Multiwfn波函数分析程序的意义、功能与用途》(http://sobereva.com/184)。在阅读本文之前应当先阅读此文《详谈Multiwfn支持的输入文件类型、产生方法以及相互转换》(http://sobereva.com/379)了解一些基本信息。


2 结合PSI4在CCSD(T)级别下做波函数分析

PSI4网址为http://www.psicode.org。Linux下安装方法:去官网下载比如psi4conda-1.1-py36-Linux-x86_64.sh,然后运行之,按照提示操作即可,会自动下载所有要用的小程序和库。最后按照提示在用户目录下的.bashrc里添加环境变量,如
export PATH=/sob/psi4/bin:$PATH
export PSI_SCRATCH=/sob/psi4scr

下面示例输入文件用于产生氟化氢的CCSD(T)级别的密度矩阵,并在当前目录下输出HF_CCSDpT.fch文件。PSI4默认是不冻核的。
molecule HF {
H        0.0        0.0       -0.831975
F        0.0        0.0        0.092442
}

set basis cc-pVTZ
grad, wfn = gradient('CCSD(T)', return_wfn=True)
fchk_writer = psi4.FCHKWriter(wfn)
fchk_writer.write('HF_CCSDpT.fchk')


将上面输入文件的内容存到比如test.inp里,然后运行psi4 test.inp test.out,即开始计算并产生输出文件test.out。

得到的HF_CCSDpT.fch里面记录的轨道是参考态HF级别的,而记录的密度矩阵是CCSD(T)的。Multiwfn是基于轨道计算的,从输入文件里载入的也是轨道而非密度矩阵。因此,若以正常方式令Multiwfn载入这个.fch,程序只会载入HF轨道,因此之后分析的也都是HF级别的情况。为了能让Multiwfn分析CCSD(T)密度,需要做以下步骤:
(1)照常将HF_CCSDpT.fch载入Multiwfn
(2)进入主功能1000(隐藏选项)的子功能98。此时Multiwfn会从.fch中寻找Total SCF Density字段并读取密度矩阵(虽然名字是SCF,但PSI4是将CCSD(T)的密度矩阵写入了这里)
(3)Multiwfn立刻输出了将CCSD(T)密度矩阵对角化得到的自然轨道占据数,并且此时内存里的密度矩阵、基函数的系数矩阵也都是CCSD(T)自然轨道的了。然而,内存中的GTF信息仍然是HF级别的,因此之后做与实空间函数有关的分析也依然是HF级别的。如果想分析CCSD(T)级别下的与实空间函数有关的问题,需要做下面一步
(4)进入主功能100的子功能2,选择导出.molden文件,比如叫t.molden。然后重启Multiwfn,载入t.molden,此时内存里所有波函数信息就都对应于CCSD(T)了,之后任何分析也都是CCSD(T)级别的了

如果.fch文件是开壳层体系,在步骤(2)时还会从.fch文件中读取自旋密度矩阵,并得到Alpha和Beta密度矩阵。程序会问你是否想产生UNO轨道,如果选n,则会分别产生Alpha和Beta自然轨道;如果选y,那么会基于总密度矩阵得到空间自然轨道(即不区分自旋)。


3 结合MRCC在高阶耦合簇/CI级别下做波函数分析

MRCC网址为https://www.mrcc.hu,下载预编译的包,解压到本机,将目录添加到PATH环境变量后即可直接使用。输入文件名必须叫MINP。运行诸如dmrcc |tee out.txt,就会启动MRCC,读取当前目录下的MINP,并将信息同时输出到屏幕和out.txt文件中。

下面的输入文件内容用于在CCSDT/cc-pVTZ下计算氟化氢单点,并且产生相应级别的密度矩阵(若不写dens=1则只计算能量)。
basis=cc-pvtz
calc=CCSDT
mem=2500MB
dens=1

geom=xyz
2

H        0.0        0.0       -0.831975
F        0.0        0.0        0.092442

计算完毕后,会在当前目录下产生MOLDEN文件,这是Molden输入文件,里面记录的是参考态HF轨道。还输出了CCDENSITIES文件,可以用文本编辑器打开,第一部分记录了二阶约化密度矩阵(2RDM),有四个标号;之后记录的是一阶约化密度矩阵(1RDM),有两个标号(另外两个标号都是0)。这些密度矩阵以MO为基,标号是从第一个考虑电子相关的MO开始计的。我们也得像上一节那样,把1RDM转化成自然轨道存到.molden文件里,才能被Multiwfn利用并在当前级别下做波函数分析。

将MOLDEN改名为比如mrcc.molden,就可以被Multiwfn读取。载入它之后,进入主功能1000的子功能97,输入CCDENSITIES文件的路径,然后输入有多少个轨道被冻结。MRCC默认是冻核的,当前体系有两个内核电子(输出文件中可看到Number of core electrons:     2),当前体系又是闭壳层的,每个占据轨道有俩电子,所以冻结的是1个轨道,因此输入1。之后Multiwfn马上会输出自然轨道占据数,会看到那个冻核的轨道占据数是精确的2.0,其它都是非整数。之后,如果打算分析与实空间函数有关的问题,也是先按照上一节说的方式导出.molden文件,然后再载入之,之后的各种分析就都是CCSDT级别的了。

MRCC做各种CI,以及FCI,也都是把密度矩阵输出到CCDENSITIES里。基于CI密度进行分析和上述完全相同。MRCC里做FCI很简单,把calc=后面写FCI即可。下面是在FCI/aug-cc-pVDZ下计算拉长的LiH的输入文件,不做冻核近似
basis=aug-cc-pvdz
calc=fci
mem=2500MB
dens=1
core=0

geom=xyz
2

H        0.0        0.0      0.0
Li       0.0        0.0      3.0

评分

参与人数 7eV +35 收起 理由
ABetaCarw + 5 赞!
ZCSco + 5 GJ!
muxijiao + 5 好物!
冰释之川 + 5 好物!
978142355 + 5 好物!
让你变成回忆 + 5 GJ!
captain + 5 好物!

查看全部评分

北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

26

帖子

0

威望

152

eV
积分
178

Level 3 能力者

发表于 2018-1-13 17:19:49 | 显示全部楼层
sob老师,multiwfn可以分析基于casscf密度矩阵的波函数吗?

1万

帖子

25

威望

1万

eV
积分
34206

管理员

公社社长

 楼主| 发表于 2018-1-14 21:32:33 | 显示全部楼层
snow_wl 发表于 2018-1-13 17:19
sob老师,multiwfn可以分析基于casscf密度矩阵的波函数吗?

可以
只要密度矩阵信息以自然轨道形式存在于输入文件即可
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

6

帖子

0

威望

23

eV
积分
29

Level 2 能力者

发表于 2018-2-9 11:01:43 | 显示全部楼层
您好,能发一下MRCC的安装包吗?我去官网上没下载到

1万

帖子

25

威望

1万

eV
积分
34206

管理员

公社社长

 楼主| 发表于 2018-2-9 13:32:04 | 显示全部楼层
lidanhui2 发表于 2018-2-9 11:01
您好,能发一下MRCC的安装包吗?我去官网上没下载到

去量化版的资源分享版
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

6

帖子

0

威望

23

eV
积分
29

Level 2 能力者

发表于 2018-2-9 15:51:30 | 显示全部楼层
老师,我的阅读权限达不到,您能发我邮箱里吗?940469319@qq.com十分感谢

1

帖子

0

威望

7

eV
积分
8

Level 1 能力者

发表于 2018-9-17 20:18:59 | 显示全部楼层
老师,您好,用psi4 计算CCSD(T)电子密度时,出现的错误如下dpd_block_matrix: No memory left. 默认内存是256MB,节点内存够。但还是出错,您知道这是什么原因吗?

1万

帖子

25

威望

1万

eV
积分
34206

管理员

公社社长

 楼主| 发表于 2018-9-18 20:27:40 | 显示全部楼层
lidanhui2 发表于 2018-2-9 15:51
老师,我的阅读权限达不到,您能发我邮箱里吗?十分感谢

直接去官网申请,不要钱
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

1万

帖子

25

威望

1万

eV
积分
34206

管理员

公社社长

 楼主| 发表于 2018-9-18 20:28:16 | 显示全部楼层
Danhuili 发表于 2018-9-17 20:18
老师,您好,用psi4 计算CCSD(T)电子密度时,出现的错误如下dpd_block_matrix: No memory left. 默认内存是 ...

我不确定。PSI4有官方论坛,你可以把输入输出文件都发到那里看开发者怎么说
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

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

GMT+8, 2018-9-26 13:47 , Processed in 0.243911 second(s), 25 queries .

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