低标度量子化学程序LSQC2.4新版本发布
发布时间:2019-11-10 15:02:59

低标度量子化学程序(Low Scaling Quantum Chemistry, LSQC)是用于大体系线性标度或低标度电子结构计算的程序包,由南京大学黎书华教授、李伟副教授课题组开发。LSQC 1.0版于2006420日发布,当前版本为2019101日发布的2.4版。本文简单介绍LSQC程序包的功能和特点,欢迎有大分子体系计算需求的老师和同学下载使用。

 

LSQC程序目前主要包含黎书华教授课题组开发的普适的基于能量分块方法Generalized Energy-Based Fragmentation, GEBF)和分子中的簇Cluster in Molecule, CIM)局域相关方法,在程序中分别对应GEBF模块和CIM模块。这两种方法均为线性标度方法,计算时间随体系增大仅线性增加,且可以高度并行化。

 

GEBF模块简介

LSQC中的GEBF模块是GEBF方法的一种高效的实现。

1-191110151240221.png 

GEBF的流程简述如下

(1) 对目标体系进行分块。

(2) 组合为静电嵌入的初始子体系(primitive subsystem),并用氢原子封端。

(3) 根据容斥原理推出衍生子体系(derived subsystem),保证目标体系的每个原子不重不漏。

(4) 用传统量子化学程序计算子体系的能量或能量梯度(其他性质)。

(5) 组合这些能量或能量梯度(其他性质)得到目标体系的能量或能量梯度(其他性质)。

目标体系的能量表达为


1-191110152239405.png 

(1)

其中1-191110151459410.png表示第m个子体系(包含背景点电荷)的能量,Cm表示第m个子体系的线性组合系数,1-191110151431335.png表示原子A上的电荷,M表示子体系的总个数。

 

目标体系的能量梯度为


1-191110151514155.png


 

(2)

其中A表示给定子体系中的原子, a  b 表示点电荷中心,1-191110152251960.png表示第m个子体系作用在电荷中心a的电场, 1-191110152335B5.png表示ab上电荷的库仑相互作用


1-191110152342D0.png 

(3)

为了省时,程序实际采用简化公式


1-191110152356128.png 

(4)

对原子上的梯度进行计算,并用之进行结构优化、振动频率计算等。目标体系的性质为


1-19111015240U24.png 

(5)

其中Ωm是第m个子体系的物理性质。

 

当子体系使用X方法(如B3LYP)时,记为GEBF-X方法(如GEBF-B3LYP)。子体系的计算可由现有的量子化学软件完成,当前版本只支持Gaussian程序(G09G16),在后续版本中会支持更多的程序。

 

使用GEBF方法,目前可以实现对数百乃至数千原子的体系进行全量子计算包括分子团簇、多肽链、蛋白质和折叠体等多种大体系的计算。

 

LSQC 2.4版支持的方法如下:

半经验方法(AM1PM3PM6,取决于高斯支持的半经验方法

Hartree-Fock

DFT(包括B3LYPM06-2X,取决于高斯支持的泛函,支持Grimme色散校正

电子相关方法(MP2MP3MP4CCSDCCSD(T)

 

支持以上方法水平下,如下性质的计算

单点能

结构优化

频率计算

红外强度

Raman强度

零点能

Gibbs自由能

偶极矩

静态极化率

化率

核磁共振

 

应用实例:

1-191110152424954.png 

GEBF-M06-2X/6-31G(d)水平结构优化

582原子,6126基函数

计算时间:15小时 (16)

 

 

CIM模块简介

CIM模块是CIM方法的一种高效的实现。CIM局域相关方法适合于对大体系进行post-HF计算。在CIM框架下,体系的电子相关能被表达为所有占据轨道贡献的加和:


1-191110152553564.png 

(6)

MP2CCSD方法为例,第i个占据轨道的贡献为


1-191110152602634.png 

(7)

其中1-191110152615222.png是双电子积分, MP2方法中1-19111015262c08.pngCCSD方法中1-191110152646419.png上式对正则分子轨道和局域分子轨道都是成立的。

 

CIM方法中,首先对整个分子进行HF计算,并将所得占据轨道局域化。每个局域占据分子轨道分别充当一次中心分子轨道,以其为中心,加上与之空间相邻的占据轨道和空轨道,共同构成一个轨道簇,记为簇P。这样,对于给定的局域占据分子轨道,其相关能的贡献可以通过解簇PMP2CC方程而获得:


1-191110152J2296.png

(1.8)

 

轨道相关能在X方法(如MP2)下求得,则记为CIM-X方法(如CIM-MP2方法)。

 

CIM方法与GEBF方法类似,都可进行大体系的后HF计算。但由于CIM方法基于轨道对体系进行分片,相比GEBF方法更加黑箱化。部分体系,如金属氧化物、石墨烯分子片等难于进行基于原子的分片操作的体系,可以使用CIM方法计算。目前HF计算依赖免费开源程序PySCF

 

LSQC 2.4版支持的方法如下:

CIM-MP2

CIM-RI-MP2

目前仅支持以上两种方法的单点计算。同时当前版本也支持中等大小体系的传统MP2RI-MP2计算。

 

应用实例:

1-191110152QD04.png 

CIM-RI-MP2/def2-TZVP 单点

884原子,19660基函数

计算时间:36小时 (24)

 

获取方式:

LSQC为免费程序,可以在官方网站:http://itcc.nju.edu.cn/lsqc  上下载使用。官方还提供了详细的安装和使用文档,欢迎大家使用。