New Features in Q-Chem 7
Q-Chem 7 is coming in July 2026! Q-Chem 7 includes all of the performance and usability features you've come to expect from Q-Chem (like our parallel RI algorithms for fast DFT performance, and our Robust SCF algorithm for easy convergence), along with a plethora of new features.
Get a sneak peek at some of the anticipated new features below!
QC-PBC
QC-PBC is a new module for handling solid-state systems and materials by all-electron calculations with periodic boundary conditions and Gaussian basis sets.
- DFT and TDDFT
- MP2, LT-MP2, and MP3 for gamma and k-point calculations
- CCSD, CCSD(T), and CCSDT for gamma-point calculations
- Python interfacing
- Geometry optimization
- Solvation
- Analytic frequency and phonon calculations
M-Chem
M-Chem is a new module for handling large biomolecular systems.
- High-performance MPI/OpenMP hybrid implementation for molecular dynamics simulations
- AMOEBA and MBUCB force fields
- Nose-Hoover thermostat and barostat
- MD in vacuum and single-point calculations with AMOEBA
Density Functional Theory
- Extended tight binding DFT (xtb) energy and gradient (Rebecca Tomann, Martin Head-Gordon) (Feature)
- COACH functional: A new range-separated hybrid (RSH) meta-GGA that is more accurate and transferable than the best existing RSH meta-GGAs, such as ωB97M-V, for a wide variety of systems. (Jiashu Liang, Martin Head-Gordon) (https://arxiv.org/abs/2603.23466)
- Faster MP2 and double-hybrid DFT: New algorithm provides competitive scaling with DLPNO for MP2 and double-hybrid DFT (was this merged in [48085/qchem]?), with improved error control for higher accuracy. (Zhenling Wang, Yao Shen, Martin Head-Gordon) (https://pubs.acs.org/doi/10.1021/acs.jctc.5c01358)
- New complex-variable DFT functionals for handling electronically metastable states, including LDA (Slater X, VWN5 C, VWN1RPA C), GGA (B88 X, PBE X, PBE C, LYP C) and hybrid (all hybrid functionals composed of the above LDAs and GGAs and possibly exact HF exchange, including PBE0, B3LYP and BH&HLYP). (Charlotte Titeca, Yifan Jiang, Thomas-C. Jagau) (https://doi.org/10.1021/acs.jpclett.5c04034)
- ECD with TDDFT (Xunkun Huang, WanZhen Liang)
Correlated Methods
- MRSF-TDDFT with properties: MRSF-TDDFT is an improved version of spin-flip TDDFT that enables effectively spin-pure treatments of doubly excited states, bond-breaking, conical intersections, and some other cases of strongly correlated systems. Q-Chem 7 presents effective MRSF-TDDFT implementation including calculation of state and transition properties, such as state and transition dipole moments, oscillator strengths, spin–orbit couplings, and density-matrix based analyses of the MRSF states and transitions. (Arnab Chakraborty, Zheng Pei, Yihan Shao, Anna I. Krylov) (https://doi.org/10.26434/chemrxiv.15002538/v2)
- New excited-state analysis features based on Earth Mover's Distance (Zhe Wang, Martin Head-Gordon) (https://doi.org/10.1021/acs.jctc.3c00894)
- Charge-displacement metrics of TDDFT in libwfa. These metrics are based on are rigorously invariant with respect to orbital rotations, unlike earlier metrics such as Tozer’s Lambda, and metrics widely used in the Gaussian program (John Herbert) (https://doi.org/10.1021/acs.jctc.4c01085)
- Large speedups in RI-CC2 ground state code in libgmbpt (Hrishikesh Ram and Martin Head-Gordon)
- CC2 and RI-CC2 Dyson orbitals for EA, IP, EE-EA, and EE-IP (Mauro Gascón Navas, Thomas C. Jagau, Robin E. Moorby, Simen Camps, Tianyi Gao)
- THC-sRI Improvements: Includes RI-CC2 and sRI-CC2 oscillator strengths; THC-sRI-CC2 ground and excited state energies and properties; and THC-sRI-CCSD ground and excited state energies (Chongxiao Zhao, Ruihao Bi, Qi Ou, Joohno Lee, Chenyang Li, Wenjie Dou)
- Performance improvements for RI and CD CCSD and EOM-CCSD
- New tools for calculating Auger decay rates using CC/EOM-CC
- NMR chemical shifts improvements (Xiao Liu, Martin Head-Gordon)
- New MPI/OpenMP parallel finite first differences algorithm using RI gives very large speedups over Q-Chem’s legacy NMR code
- Support for all modern density functionals through hybrids
Molecular Dynamics, Non-adiabatic Dynamics, Embedding, and Solvation
- New embedding interface for external packages like CP2K (Ar Fonlon, Elena Kolodzeiski, Dustin R. Broderick, Kay Carter-Fenk, Christopher J. Mundy, Christopher J. Stein, John M. Herbert)
- An order of magnitude total speed up for large QM/EFP jobs via parallelization of one-electron integrals (Lyudmila Slipchenko)
- Pairwise harmonic confiner for optimization in internal coordinates (Chance Brandt, John Herbert)
Fragment and Energy Decomposition Analysis
- Correlated wavefunction EDA: EDA-II is available for MP2, BW-s2, κ-MP2 using the linear-scaling codebase (Zhenling Wang, Hengyuan Shen, Martin Head-Gordon) (https://doi.org/10.1021/acs.jctc.4c01301)
- EDA OVOCV analysis and OODFT OVOCV analysis: OVOCV can be used to replace or complement NOCV analysis in EDA; it provides clearer identification of donor and acceptor orbitals than NOCV methods. (Hengyuan Shen, Martin Head-Gordon) (https://doi.org/10.1021/acs.jpca.4c02364)
- Energy Decomposition Analysis for Restricted Open-shell Kohn-Sham Theory (ROKS-EDA) (Haobo Ling, Hengyuan Shen, Martin Head-Gordon) (https://doi.org/10.26434/chemrxiv-2025-5t48c)
- Two new MBD dispersion models for XSAPT+MBD:
- MBDrev, based on machine learning applied to the SAPT10k data set (Corentin Villot and Ka Un Lao) (https://doi.org/10.1063/5.0204064)
- oC8raim, with improved performance for ions and extended to heavy elements (Keegan Paice, John Herbert)
Incorporation of Quantum Nuclear Effects (NEO Suite)
- NEO-AIMD methods for NEO-HF and NEO-DFT, including:
- NEO Born-Oppenheimer MD (NEO-BOMD) (Joseph Dickinson, Sharon Hammes-Schiffer)
- Extended Lagrangian MD (NEO-ELMD) (Joseph Dickinson, Sharon Hammes-Schiffer)
- Constrained NEO MD (CNEO-MD) (Joseph Dickinson, Sharon Hammes-Schiffer)
- NEO-PCM and NEO-QM/MM (Joseph Dickinson, Sharon Hammes-Schiffer)