How XMCQDPT Works¶
The Extended Multi-Configuration Quasi-Degenerate Perturbation Theory (XMCQDPT) method builds upon a converged State-Averaged CASSCF calculation by introducing dynamic electron correlation through second-order perturbation theory.
Unlike CASSCF, XMCQDPT does not optimize molecular orbitals or modify the active space. Instead, it uses the optimized CASSCF wavefunctions as reference states and computes perturbative corrections arising from excitations outside the active space.
The overall workflow is illustrated below.
Hartree–Fock
│
▼
State-Averaged CASSCF
│
▼
CAS-CI Reference States
│
▼
Construct Zeroth-Order Hamiltonian
│
▼
Rotate to Intermediate Basis
│
▼
Evaluate Second-Order Perturbation Corrections
│
▼
Construct Effective Hamiltonian
│
▼
Diagonalize Effective Hamiltonian
│
▼
Final XMCQDPT Energies
Step 1 – Read the CASSCF Reference¶
The XMCQDPT calculation begins by reading the optimized molecular orbitals and CI wavefunctions from the previous CASSCF calculation using
The reference contains
- optimized molecular orbitals,
- optimized active-space occupations,
- configuration interaction (CI) coefficients,
- state-averaged electronic wavefunctions.
Unlike RHF or CASSCF, no additional SCF or orbital optimization is performed.
Step 2 – Construct the CAS-CI Reference States¶
Using the optimized orbitals, GAMESS performs one final CI diagonalization within the active space.
Optimized Orbitals
│
▼
Generate Configuration State Functions
│
▼
Diagonalize the CI Hamiltonian
│
▼
CAS-CI Reference States
These are the states printed under
At this stage, the wavefunction already contains static electron correlation, but dynamic correlation is still absent.
Step 3 – Construct the Zeroth-Order Hamiltonian¶
The next step is the construction of the zeroth-order effective Hamiltonian, denoted as
Unlike single-reference perturbation methods such as MP2, XMCQDPT treats multiple electronic states simultaneously.
Instead of working with independent states, GAMESS constructs an effective Hamiltonian of the form
The diagonal elements correspond to the reference-state energies, while the off-diagonal elements describe coupling between nearby electronic states.
This matrix is printed in the output as
Step 4 – Rotate to the Intermediate Basis¶
Electronic states that are close in energy often interact strongly.
Direct application of perturbation theory may therefore become numerically unstable because of small energy denominators.
To overcome this problem, XMCQDPT first diagonalizes the zeroth-order Hamiltonian,
and rotates the original CAS-CI states into an intermediate basis.
The transformation matrix is printed under
The intermediate states are
- orthogonal,
- weakly coupled,
- better suited for perturbation theory.
This intermediate-state formalism is the defining feature of XMCQDPT and greatly improves the treatment of nearly degenerate electronic states.
Step 5 – Evaluate Dynamic Correlation¶
Once the intermediate states have been constructed, GAMESS evaluates the interaction between the reference wavefunctions and the external configuration space.
These interactions arise from electronic excitations outside the active space.
Rather than explicitly including millions of excited configurations in the CI expansion, XMCQDPT estimates their contributions using second-order perturbation theory.
This step accounts for dynamic electron correlation, which is absent in the CASSCF description.
Step 6 – Construct the Second-Order Effective Hamiltonian¶
After all perturbative corrections have been evaluated, GAMESS constructs the second-order effective Hamiltonian
This matrix contains
- second-order energy corrections,
- corrected interstate couplings,
- correlation effects arising from the external configuration space.
It is printed in the output as
Step 7 – Obtain the Final XMCQDPT States¶
The second-order effective Hamiltonian is finally diagonalized to obtain the corrected electronic states.
The resulting energies are reported under
while the corrected wavefunctions are printed in
These are the final electronic energies that should be used for
- excitation energies,
- potential energy surfaces,
- photochemical studies,
- conical intersection analysis,
- spectroscopy calculations.
Static vs Dynamic Correlation¶
It is useful to distinguish the roles of CASSCF and XMCQDPT.
| Method | Captures |
|---|---|
| RHF | Mean-field description |
| CASSCF | Static (strong) electron correlation |
| XMCQDPT | Dynamic electron correlation |
CASSCF correctly describes situations where multiple electronic configurations contribute significantly to the wavefunction, while XMCQDPT accounts for the remaining correlation arising from instantaneous electron-electron interactions.
Together, these methods provide a balanced and highly accurate description of ground and excited electronic states.
Summary¶
The XMCQDPT procedure can be summarized as follows:
- Read the optimized CASSCF wavefunctions.
- Construct the CAS-CI reference states.
- Build the zeroth-order effective Hamiltonian.
- Rotate the reference states into an intermediate basis.
- Evaluate second-order perturbative corrections from external excitations.
- Construct the second-order effective Hamiltonian.
- Diagonalize the effective Hamiltonian to obtain the final XMCQDPT energies and wavefunctions.
Because it combines the multiconfigurational description of CASSCF with dynamic correlation from perturbation theory, XMCQDPT has become one of the most reliable methods for studying excited-state potential energy surfaces, photochemical reactions, and systems with strong electronic near-degeneracy.