查看原文
其他

离线安装PySCF-1.7.6

jxzou 量子化学 2023-12-29

由于PySCF程序更新较为频繁,不少小伙伴反映按照本公众号之前发过的离线安装教程《离线安装PySCF程序(1.5及更高版本)》经常会出现编译无法通过等问题,因此再次更新一下安装教程。作为一款量子化学软件,PySCF现在十分强大,支持的功能非常多,此处随意列举几个:


  • 密度拟合加速的HF和CASSCF功能

  • 二分量、四分量相对论Hartree-Fock

  • 与其他程序有方便的接口进行DMRG、FCIQMC、SHCI等计算

  • 多参考微扰方法SC-NEVPT2

  • 各类EOM方法IP/EA/EE-EOM-CCSD

  • CCSD(T)解析梯度

  • DFT二阶解析导数、TDDFT解析梯度

  • ADC激发态方法

  • G0W0近似方法

  • 支持分子和周期性体系的NAO、IAO、IBO轨道及PM局域化等方法

 

在可预见的将来很可能会成为使用人数仅次于Gaussian和ORCA的量子化学软件。当然,PySCF还是免费开源的,有不少开发者基于它开发程序或插件。关于PySCF更多特点介绍请见

https://github.com/pyscf/pyscf/blob/master/FEATURES

 

注意PySCF联网在线安装只需pip install pyscf一行命令即可。本文介绍的是离线安装步骤,适合不允许联网或很难联网的内部节点。读者在开始编译前需确认自己机子上有gcc和g++编译器,有MKL数学库,以及cmake软件。运行如下命令可查看自己机子上是否存在

which gccwhich g++which cmakecmake --version (查看版本号)echo $MKLROOT

 

笔者撰文时用的gcc版本为4.8.5(更高版本当然也可以),cmake版本为3.19(不能低于3.5),MKL数学库用的是Intel Parallel Studio XE 2019 update 5里的。读者机子上若缺少任一前提条件请自行安装。

 

1. 下载

先到GitHub网站下载PySCF压缩包

https://github.com/pyscf/pyscf/releases

此处我们以1.7.6版本为例,下载压缩包pyscf-1.7.6.tar.gz。

 

解压,进入lib目录

tar -zxf pyscf-1.7.6.tar.gzcd pyscf-1.7.6/pyscf/lib

 

打开此目录下的CMakeLists.txt文件,搜索URL或GIT_REPOSITORY可以看到所需三个库的网址和版本号。笔者将其摘选出,展示如下

GIT_REPOSITORY https://github.com/sunqm/libcint.gitGIT_TAG v4.0.7
URL https://gitlab.com/libxc/libxc/-/archive/4.3.4/libxc-4.3.4.tar.gz
GIT_REPOSITORY https://github.com/fishjojo/xcfun.gitGIT_TAG cmake-3.5

 

显然我们到对应网址去下载对应版本的压缩包即可。有GIT_TAG标签的表示需要到Tag界面点击下载对应的版本。具体做法是先点击Releases,然后点击Tags




就能找到相应的版本了。以后面对更新的PySCF版本,只要按照这种方式寻找所需库,就不会在安装过程中碰到缺库或无法兼容的报错。例如,笔者按照这些网址下载下来的压缩包分别是

  • libcint-4.0.7.tar.gz

  • libxc-4.3.4.tar.gz

  • xcfun-cmake-3.5.tar.gz

其实这三个库的安装步骤在CMakeLists.txt文件里也有写。为了更清晰地展示安装过程,以下也逐一贴出。

 

2.编译libcint

到存放压缩包的目录下,依次执行

tar -zxf libcint-4.0.7.tar.gzcd libcint-4.0.7mkdir build && cd build
cmake -DWITH_F12=1 -DWITH_RANGE_COULOMB=1 -DWITH_COULOMB_ERF=1 \-DMIN_EXPCUTOFF=20 -DKEEP_GOING=1 \-DCMAKE_INSTALL_PREFIX:PATH=/home/$USER/software/cint_and_xc \-DCMAKE_INSTALL_LIBDIR:PATH=lib ..
make && make install

 

这里的库存放路径/home/$USER/software/cint_and_xc是笔者的个人偏好,读者可以按照自己的喜好或需求更改。

 

3. 编译libxc

到存放压缩包的目录下,依次执行

tar -zxf libxc-4.3.4.tar.gzcd libxc-4.3.4/build
cmake -DCMAKE_BUILD_TYPE=RELEASE -DBUILD_SHARED_LIBS=1 \-DCMAKE_INSTALL_PREFIX:PATH=/home/$USER/software/cint_and_xc \-DCMAKE_INSTALL_LIBDIR:PATH=lib ..
make -j8make install

 

其中-j8表示8核并行编译,若读者机子上没这么多核,请减小核数或直接删去这一参数。

 

4. 编译xcfun

到存放压缩包的目录下,依次执行

tar -zxf xcfun-cmake-3.5.tar.gzcd xcfun-cmake-3.5mkdir build && cd build
cmake -DCMAKE_BUILD_TYPE=RELEASE -DBUILD_SHARED_LIBS=1 \-DXCFUN_MAX_ORDER=3 \-DCMAKE_INSTALL_PREFIX:PATH=/home/$USER/software/cint_and_xc \-DCMAKE_INSTALL_LIBDIR:PATH=lib ..
make -j8make install


上述四步完成后,可以在/home/$USER/software/cint_and_xc下发现有bin、include、share和lib四个文件夹。接着将lib路径添加进环境变量,即打开~/.bashrc文件写入

export LD_LIBRARY_PATH=/home/$USER/software/cint_and_xc/lib:$LD_LIBRARY_PATH

然后source使之生效。

 

5.安装PySCF

我们回到本文一开始提到的pyscf-1.7.6/pyscf/lib目录,再次打开CMakeLists.txt文件,找到第一个# set(BLAS_LIBRARIES"-L/这一行,在此处删除注释符号“# ”,更改mkl库路径为当前系统下的mkl路径,例如笔者机子上的是/opt/intel/mkl/lib/intel64。注意别把前头的-L删了。下面几行“or”部分不用动。接着执行

mkdir build && cd build
cmake -DBUILD_LIBCINT=0 -DBUILD_LIBXC=0 -DBUILD_XCFUN=0 \-DCMAKE_INSTALL_PREFIX:PATH=/home/$USER/software/cint_and_xc ..
make


安装完成后将PySCF的路径添加进~/.bashrc文件:

export PYTHONPATH=/home/$USER/software/pyscf-1.7.6:$PYTHONPATH


完成安装。安装包和解压出的文件夹都可以删除,只留下cint_and_xc和pyscf-1.7.6文件夹即可。最后同样要记得执行source ~/.bashrc,或者退出重登。

 

6.测试例子

随便测试个CCSD(T)计算。新建一个文件a.py,写入

from pyscf import gto, scf, ccmol = gto.M(atom='H 0 0 0; F 0 0 1.1',basis='ccpvdz')mf = scf.RHF(mol)mf.kernel()
mycc = cc.CCSD(mf)mycc.frozen = 1mycc.kernel()
mycc.ccsd_t()


保存。运行

python a.py

输出

converged SCF energy = -99.9873974403487E(CCSD) = -100.2018830963878 E_corr = -0.2144856560390928CCSD(T) correction = -0.00239622928398965


读者可以自行用高斯做个计算对比。注意高斯在电子相关计算中默认冻结芯轨道,而PySCF默认不冻结。对于氟化氢这个例子需要冻结的轨道只有1个,即F原子的1s轨道,因此这里显式地设定了mycc.frozen = 1让其与高斯一致。

继续滑动看下一个

您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存