
Q-Chem是一款功能齐全的从头算量子化学程序包,最早由诺贝尔化学奖获得者John A. Pople主导开发。全面支持从DFT/HF到各种高级的post-HF相关能计算方法。并提供广泛的解决方案,适用于不同的研究目标,如单分子磁体的自旋轨道耦合效应、高通量计算、酶的量子力学/分子力学(QM/MM)研究等。被广泛用于处理工业界、学术界和国家实验室的各类理论模拟研究,在理论化学、药物设计、材料科学、生物化学以及相关领域的教学和研究中发挥了极大作用。可以包括研究分子结构、化学反应、分子振动、电子光谱、NMR谱和溶剂化效应等。
Q-Chem自1999年第一个版本1.2发布至今,已经陪伴广大科研用户26年了。自1999年以来,Q-Chem引领量子化学计算的潮流,不断突破科技边界,26年的专注研发,每年至少三次的版本迭代,Q-Chem始终站在量子化学软件的最前沿。 从1.2版本到最新的6.2,Q-Chem的每一次版本升级都是对科学的一次深刻致敬。
Q-Chem7.0新功能
Q-Chem 7.0版本在周期性体系、生物大分子、激发态理论、电子相关方法及量子核效应等方向均有实质性扩展。
QC-PBC:周期性边界条件下的全电子计算
Q-Chem 7.0 引入QC-PBC模块,支持以Gaussian 型基组开展周期性体系的全电子计算。目前实现的方法包
基态与激发态的 DFT 及 TDDFT
MP2、LT-MP2 与 MP3,支持Γ点与 k点取样
CCSD、CCSD(T) 及 CCSDT,限于Γ点
几何结构优化、溶剂化模型、解析Hessian及声子计算
Python 接口,支持外部工作流调用
该模块的引入使 Q-Chem 从分子气相/溶液体系扩展至晶体、表面及材料体系,且保持了全电子精度。
M-Chem:生物大分子体系模拟
M-Chem 模块面向生物大分子体系,主要技术优势如下:
基于 MPI/OpenMP 混合并行的分子动力学(MD)引擎
支持固定电荷(Amber)、可极化(AMOEBA)及ReaxFF力场
Nose-Hoover 热浴与压浴、共轭梯度自洽场、扩展 Lagrangian 多体力求解方案
周期性边界条件与 Particle Mesh Ewald(PME)长程静电处理
模块化前端,用于标准力场及蛋白-水体系的参数分配与体系构建
QM/MM 与 ReaxFF 的耦合框架
兼容外部分析代码的轨迹输出格式

密度泛函理论的方法学进展
COACH 泛函
COACH为新型范围分离杂化(RSH)meta-GGA 泛函。与现有 RSH meta-GGA(如ωB97M-V)相比,其在多种体系上的精度与可迁移性均有提升(Jiashu Liang, Martin Head-Gordon)。
线性标度局域相关
针对 MP2 及双杂化 DFT,Q-Chem 7.0 实现了新的线性标度局域相关算法。该算法在计算标度上可与 DLPNO 竞争,同时采用单一数值阈值控制精度,误差控制更为直接(Zhenling Wang, Yao Shen, Martin Head-Gordon)。

复变量 DFT
新增复变量 DFT 泛函,用于处理电子亚稳态(如共振态、散射态)。已实现的泛函涵盖:
LDA:Slater 交换、VWN5/VWN1RPA 关联
GGA:B88、PBE 交换,PBE、LYP 关联
杂化泛函:由上述 LDA/GGA 组分与精确 HF 交换构建,包括 PBE0、B3LYP、BH&HLYP 等(Charlotte Titeca, Yifan Jiang, Thomas-C. Jagau)。
扩展紧束缚(xtb)
新增 xtb 方法的能量与梯度计算,适用于大体系快速预优化或初筛(Rebecca Tomann, Martin Head-Gordon)。
激发态与电子相关方法
MRSF-TDDFTMRSF-TDDFT 是 Spin-Flip TDDFT 的改进形式,可对双重激发态、键解离、锥形交叉及强关联体系实现有效自旋纯处理。Q-Chem 7.0 实现了该方法的完整性质计算框架,包括:
态与跃迁偶极矩
振子强度
自旋-轨道耦合
基于密度矩阵的态与跃迁分析(Arnab Chakraborty, Zheng Pei, Yihan Shao, Anna I. Krylov)。
激发态分析
基于 Earth Mover's Distance 的激发态分析工具(Zhe Wang, Martin Head-Gordon)
libwfa 中的电荷位移度量,严格满足轨道旋转不变性,区别于传统 Tozer's Lambda 及 Gaussian 程序中广泛使用的早期度量(John Herbert)
耦合簇方法加速与扩展
RI-CC2 基态代码在 libgmbpt 中实现显著加速(Hrishikesh Ram, Martin Head-Gordon)
CC2 与 RI-CC2 的 Dyson 轨道计算,支持 EA、IP、EE-EA、EE-IP 过程(Mauro Gascón Navas, Thomas C. Jagau, Robin E. Moorby, Simen Camps, Tianyi Gao)
THC-sRI 系列改进:RI-CC2 与 sRI-CC2 振子强度;THC-sRI-CC2 基态与激发态能量及性质;THC-sRI-CCSD 基态与激发态能量(Chongxiao Zhao, Ruihao Bi, Qi Ou, Joohno Lee, Chenyang Li, Wenjie Dou)
RI 与 CD 框架下 CCSD 及 EOM-CCSD 的性能优化
基于 CC/EOM-CC 的俄歇衰变速率计算工具
NMR 化学位移
全新 MPI/OpenMP 并行有限差分算法,结合 RI 技术,相对旧版 NMR 代码实现大幅加速
支持所有现代杂化密度泛函(Xiao Liu, Martin Head-Gordon)
能量分解与片段分析
相关波函数
EDAEDA-II 现已支持 MP2、BW-s2 及 κ-MP2,基于线性标度代码库实现(Zhenling Wang, Hengyuan Shen, Martin Head-Gordon)。
OVOCV 分析
EDA 与 OODFT 中引入 OVOCV(Orbital-Valence Occupation Change based on Valence)分析,可作为 NOCV 的替代或补充,在给体-受体轨道识别上具有更高分辨率(Hengyuan Shen, Martin Head-Gordon)。
ROKS-EDA
实现限制性开壳层 Kohn-Sham(ROKS)理论框架下的能量分解分析(Haobo Ling, Hengyuan Shen, Martin Head-Gordon)。
力分解分析(FDA)
新增 FDA 功能,与 EDA 互补。
色散模型
MBDrev:基于 SAPT10k 数据集的机器学习优化模型(Corentin Villot, Ka Un Lao)
oC8raim:针对离子体系优化并扩展至重元素(Keegan Paice, John Herbert)
动力学、嵌入与溶剂化
多尺度嵌入
新增外部程序接口(如 CP2K),支持密度嵌入与 QM/MM 耦合(Ar Fonlon, Elena Kolodzeiski, Dustin R. Broderick, Kay Carter-Fenk, Christopher J. Mundy, Christopher J. Stein, John M. Herbert)
QM/EFP 框架中实现单电子积分的并行化,大型体系总计算速度提升约一个数量级(Lyudmila Slipchenko)
几何优化
内坐标优化中引入成对谐波限制器(pairwise harmonic confiner),改善优化稳定性(Chance Brandt, John Herbert)。
量子核效应:NEO 套件
Q-Chem 7.0 系统性地纳入核量子效应(Nuclear Quantum Effects),在 NEO(Nuclear-Electronic Orbital)框架下实现:
NEO-HF 与 NEO-DFT 的 Born-Oppenheimer 分子动力学(NEO-BOMD)
扩展 Lagrangian 分子动力学(NEO-ELMD)
约束 NEO 分子动力学(CNEO-MD)
NEO-PCM 与 NEO-QM/MM
上述方法由 NEO-AIMD 统一驱动,适用于质子转移、氢键动力学及同位素效应等核量子效应敏感过程(Joseph Dickinson, Sharon Hammes-Schiffer)。
光谱学工具
电子圆二色谱(ECD):支持 TDDFT(Xunkun Huang, WanZhen Liang)与 EOM-CCSD
振动圆二色谱(VCD):限制性与非限制性形式
CVS-XCIS:用于 X 射线光谱建模
俄歇衰变:新增 CC/EOM-CC 计算路径
ΔSCF:新驱动程序,简化输入并增强分析与可视化工具
Q-Chem 7.0 的方法学扩展覆盖了从周期性固体到生物大分子、从基态到激发态、从电子结构到核量子效应的多个维度。QC-PBC 与 M-Chem 的引入填补了Q-Chem在材料与生物体系模拟上的空白;线性标度局域相关、MRSF-TDDFT、THC-sRI-CC 及 NEO 套件的加入则进一步强化了其在高精度电子结构与动力学模拟上的能力。
Q-Chem6.4新功能
Q-Chem 6.4 震撼来袭!立即升级,体验显著提升的计算性能与易用性,解锁全新的化学与光谱研究工具.依托这些创新功能,您可以部署大规模计算工作流,更快获得更精确的结果,并将研究范围拓展至以往无法触及的体系与科学问题.
光谱模拟新突破
核心-价电子分离(CVS)方案支持 ROCIS、XCIS 与 QROCIS
精确计算开壳体系的核心能级光谱
EOM-IP/EE-CCSD 求解器新增内价层投影选项
采用类 CVS 方法,精准计算由内价层电子主导的高能激发态与电离态
全新 ΔSCF 驱动程序
更快捷、更简便地执行 ΔSCF 计算,界面简洁流畅,并配备全新的结果分析工具.
大体系计算性能飞跃
迭代式 CC-in-DFT 嵌入方法用于性质计算
CC/EOM-CC 能量与梯度的两步 Cholesky 分解
显著加速 CD 计算,梯度计算效果更佳
随机 RI-CC2 解析梯度与导数耦合
DFT 方法重大更新
全新 DFT 泛函体系:
B97MV-D3s(BJ), B97MV-D4, wB97MV-D4, wB97MV-D3s(BJ), wB97XV-D3s(BJ), wB97XV-D4, Pr2SCAN69-D4,Pr2SCAN50-D4, wPr2SCAN50-D4,revDOD-PBEP86-D4, revDSD-PBEP86-D4, xDSD75-PBEP86-D4,wDSD72-PBEP86-D4,B2NC-PLYP, mPW2-PLYP,mPW2NC-PLYP,SOS0-PBE0-2
TAO-DFT 支持 B97 类泛函
PBEh-3c 与 HF-3c 半数值解析频率计算
精度与解析 Hessian 相当,计算效率更高。HF-3c 尤其适用于超高效几何优化
(在最小基组 HF 成本下,对非共价相互作用具有良好的精度).支持频率计算!
DFT-D4 全面升级
多参考态自旋翻转 DFT(MR-SF-DFT)
提升精度,解决原始 SF-TD-DFT 的自旋污染问题
NEO 方法扩展
NEO-CC2 方法用于激发态计算 NEO-CC 有限差分梯度几何优化
NEO 方法支持冻结核近似 NEO-CCSD(T) 与 NEO-CCSDTeep,epp 方法
Q-Chem6.3新功能
速度再突破
密度拟合 DFT 核心重写:中型到大型(>2000 基函数)能量与受力平均提速 20–30 %,16 核测试 Cys-11 肽链 B3LYP/def2-TZVP 单点快 23 %。
MPI 级并行覆盖有限差分频率、Numerical Hessian、Many-Body Expansion,节点数几乎线性扩展,百原子频率 overnight。
SCF 永不“卡壳”
Robust SCF:内置多阶段收敛流水线,自动识别 DIIS 漂移→切换 EDIIS/GDM/RCA/Damping,夜间作业零人工重启,专为高通量材料筛选与大数据训练设计。
光谱工具箱升级
Auger 衰减速率:全自动 Coulomb 波 + Feshbach-Fano 框架,无需外部脚本。
开壳层 Auger:Projection EOMIP-CCSD 支持冻结核心与双重态参考。
动态极化量:EOM-EE/SF-CCSD 虚频极化率,后续 CAS–London 轨道可直接读入做 Verdet 常数、CARS 强度。
开壳层 & 共振体系专属
自旋-轨道耦合:EOM-DEA、EOM-DIP 波函数均可后处理,无需重新计算。
梯度:EOM-DIP/DEA-CCSD 解析梯度,势能面扫描、束缚态-共振态混合优化稳定。
复数 EOM-CCSD:基于 RI 分解,衰变态能量、宽度一次给出,支持 200 电子级体系。
多体方法“全垒打”
EOM-CCSDT 正式版:EE/SF/IP/EA/DIP/DEA 六通道同精度,<0.05 eV 误差挑战光电子谱。
CC2-BW:Brillouin-Wigner 尺寸一致修正,克服标准 CC2 对三激发低估,激发能误差 ↓ 30 %。
化学解析看得见
Broken Bond Orbitals (BBO):一键生成,能量分解(EDA)自动区分泡利排斥、静电、轨道互作。
CT Metrics:libwfa 模块提供轨道不变电荷转移距离、激子尺寸,避免传统 Δq·Δr 伪影。
ALMO(MSDFT) Diabatic 耦合:ΔSCF 激发态与基态间电子耦合,用于单分子电导、激子裂变速率。
极端条件与反应操控
力学化学:新增 PV 能量项,外压 0–10 GPa 连续可调,几何优化与过渡态搜索自动包含体积功。
柔性约束:扭角、平底势、键长、角度、质心距、二面角、r12⊥r34 等 10 余种内置限制,支持 PES 扫描与过渡态限制。
真实环境不再“一刀切”
HetPCM:同一 QM 集群可分区赋予介电常数(表面 ε=78、疏水核 ε=4),酶催化模型更逼真。
SMD 半数值频率:解析梯度 + 数值微分,成本≈纯解析,100 原子有机溶剂频率 1–2 h 完成。
GMBE-DM:密度矩阵碎片化 SCF,可搭配 def2-TZVPP 等高阶基组,片段间自洽到 1 μEh,线性寡聚物 10× 提速。
量子核效应(NEO)套件
NEO-PCM 解析 Hessian:溶液相量子质子频率、同位素效应、质子隧道分裂全可算。
泛函“零改”支持:ωB97X-D、ωB97X-D3、B3LYP-D3(BJ)…色散校正已内置,直接做 NEO-SCF 能量与梯度。
多质子同步优化:GDM+DIIS 混合算法支持任意数目量子质子,初始核猜测改进,收敛步数减半。
Q-Chem6.2新功能
& Q-Chem 4.1新功能
& Q-Chem4.0在取得巨大成功的3.2版的优异功能的基础上,增加了:

& Q-Chem 的特色

& 基态自洽场方法
1. Hartree-Fock 方法
2. 密度泛函理论
3. 线性标度方法
4. AOINTS 包用于双电子积分
5. SCF 改进
& 最大重叠方法

& 基组
& QM/MM
& 基于波函的电子关联处理
1. MP 微扰理论
2. 局域MP2 方法
3. RI-MP2
4. 耦合簇方法
5. 优化轨道的耦合簇方法
& 激发态方法
1. 支持的计算类型
2. CIS 方法
3. TDDFT
4. 基于耦合簇的激发态方法
自旋反转激发态方法:改善了双、三自由基体系的处理;结合单行列式波函处理键断裂问题;可用于OD和CCSD 理论级别
5. 分解分析

& 特性分析
1. 自动结构优化和过渡态优化
2. 振动光谱
3. NMR 屏蔽张量
4. 自然键轨道NBO 分析
5. Stewart 原子
6. 动量密度
7. Intracules
8. 分子中的原子
9. 溶剂模型
10.基于Dirac-Fock 理论的相对论能量校正
11.对角绝热校正
& Q-Chem 创新点

& Q-Chem效用示例
- Linux机群的快速DFT计算
对于有103个原子的紫杉醇分子,在有8CPU的Xeon机群上进行BLYP/6-31(df,pd)等级(1925个基函数)的能量和梯度计算,使用Q-Chem仅需28分钟。
- MP2能达到HF的计算速度
高效的RI算法使中等的大基组的MP2计算更为有效。左图是丙氨酸四肽相对构象能(有27个稳定结构)的比较的示例。
- 线性标度NMR计算
右图和表格给出计算的分子的NMR位移相对于环绕的水分子的数目的收敛性。这是迄今进行的最大体系的NMR计算,在HF/6-31G**等级共有3000个基函数。
- 新发现
Q-Chem独有的EOM-(EE, IP, EA, SF)-CCSD方法能有效地描述双或三自由基的多个组态。使用该方法发现了首个具有开壳层双重态基态的有机分子的实例,DMX三自由基的三个不成对电子呈反铁磁性式偶合。