24小时热门版块排行榜

返回列表

beefly

专家顾问 (职业作家)

地沟油冶炼专家

专家经验: +459
QC强帖: 12
应助: 88 (初中生)
金币: 13377.5
散金: 27307
红花: 302
帖子: 4685
在线: 929小时
虫号: 408372
注册: 2007-06-21
性别: GG
专业: 理论和计算化学
管辖: 计算模拟

[交流] 让MOLPRO和CFOUR共享MRCC模块

MOLPRO和CFOUR都具有MRCC的接口，可以进行一些高级别的从头计算。但是，由于两个接口的功能不尽相同，以及MOLPRO和CFOUR两个程序各自的优缺点，不同的MRCC计算往往要用不同的接口，才能实现最佳效果。目前，MOLPRO/MRCC无法进行以下计算：UHF参考态的CI和CC，CI和CC的解析一阶导数（梯度和特性），CI和CC的解析二阶导数（Hessian和特性），这些计算需要用CFOUR/MRCC完成。但是另一方面，当存在弥散函数或重原子的情况下，CFOUR的Hartree-Fock经常不收敛，或得到奇怪的HF总能量和HOMO-LUMO能隙，需要改用MOLPRO。

MRCC的源代码需要单独申请，然后编译。此外，从MOLPRO2009开始，在某些官方发布的MOLPRO二进制串行程序中，也包含了编译好的MRCC。例如，Molpro 2010.1 x86_64 linux串行版安装之后，在Molpro根目录/bin之下，除了有molpro的程序和脚本外，还可以看到四个MRCC的程序：xmrcc，mrcc，dmrcc，goldstone。因此如果购买了MOLPRO的版权，就无须再申请和编译MRCC了，只要安装合适的串行版MOLPRO即可。

在CFOUR下面调用MRCC（无论是自己申请的，还是MOLPRO自带的），需要设定MRCC程序的路径。下面的运行脚本设定了CFOUR路径和来自MOLPRO的MRCC路径。此外，每个CFOUR计算还需要创建一个单独的临时文件目录，以避免不同的CFOUR计算互相干扰。

CODE:

#!/bin/bash

# CFour

export PATH=/data/share/software/CFOUR/CFOUR_2.00_beta/bin:$PATH

# MRCC

export PATH=/data/share/software/MOLPRO-2010S/bin:$PATH

rm -rf /tmp/C41

mkdir /tmp/C41/

cp ./GENBAS /tmp/C41/GENBAS

cp ./ECPDATA /tmp/C41/

cp ./01.inp /tmp/C41/ZMAT

cd /tmp/C41/

xcfour > /home/user1/CFOUR/01.out

输入文件01.inp定义了CCSDTQ/3-21G级别的几何优化，使用解析梯度。其中CC_PROGRAM=MRCC表示调用外部的MRCC模块。

CODE:

TEST

H

F 1 R1*

R1 =   0.9594259177

*ACES2(CALC=CCSDTQ,BASIS=3-21G

REFERENCE=RHF

CC_PROGRAM=MRCC

MULTIPLICITY=1

FROZEN_CORE=ON

MEM_UNIT=GB,MEMORY=6)

最新版本的DIRAC，ORCA和PSI4也可以调用MRCC，但是二进制MOLPRO自带的MRCC可能版本比较老，未必支持。也许要下载并安装最新的MRCC。

回复此楼