MPPT — Multipartitioning perturbation theory#

MPPT / Electronic-structure calculations

Multipartitioning
perturbation theory

Quasi-degenerate multireference perturbation theory for incomplete model spaces, with spin–orbit coupling and GRECP / GRPP relativistic pseudopotentials.

A custom PySCF interface connected to MPPT’s
Fortran and C++ computational core.

Diagrammatic MRPT

⟨Ψμ0|=∑I=1NPCIμ*⟨I|
|Ψν0⟩=∑J=1NPCJν|J⟩
Closed second-order contribution: i, j label holes; a, b label particles. This term is diagonal in the determinant basis. Diagram and indices ↗
01 / MAIN FEATURES

Methods and capabilities

Iterative growth of a selected determinant space A small seed branches through three expansions, each in a different colour. The retained tree occupies only part of a much larger complete active space. Seed+1+2+3
Select the most important determinants from enormous CAS spaces.

Determinant selection

Incomplete
model spaces

CIPSIXX builds selected or complete determinant spaces. Selection reduces the cost of both the variational calculation and the subsequent perturbative treatment.

Selection algorithm
Assembly and rediagonalization of the second-order effective Hamiltonian First diagonalize the model-space Hamiltonian to construct auxiliary states and denominators. Add the symmetric second-order correction W to PHP, then diagonalize the corrected effective Hamiltonian. +=PHPW(2)Heff[2] Diagonalize→Perturb→Diagonalize

The MPPT formulation

Multipartitioning
perturbation theory

  • Multiple partitionings. Each model determinant defines its own zeroth-order Hamiltonian.
  • Density-based denominators. Auxiliary-state density and Fock-like matrices give separate electron-removal and attachment energies: an IP/EA-like construction.
  • Wave-function relaxation. Diagonalize → perturb → diagonalize updates the state mixing, as well as the energies.
Method and equations
Scalar states coupled into spin–orbit-split energy levels HSO Scalar statesSpin–orbit levels
Relativistic pseudopotentials and spin–orbit coupling in the effective Hamiltonian.

GRECP / GRPP

Relativistic
spin–orbit calculations

GRECP, also called GRPP, accounts for scalar-relativistic effects and spin–orbit coupling. Suitable Gatchina pseudopotential parameter sets also incorporate Breit and QED contributions, including self-energy and vacuum polarization.

Relativistic formulation

Orbital generation / PySCF interface

Orbitals for compact
model spaces

Computed virtual orbitals can be spatially delocalized: the local atomic character of interest is spread across many MOs. AO projections identify linear combinations that concentrate this character in a few orbitals, giving a more compact representation of the corresponding excitation channel.

Orbital generation and projections

Delocalized virtual MOs

Illustrative local target character spread over many spatially delocalized virtual orbitals |χtarget⟩≈∑rwr|φrvir⟩

Character spread over many MOs

Rotated virtual MOs

Illustrative target character concentrated in a few rotated virtual orbitals; the full orbital span is preserved |χtarget⟩≈∑k∈Tw~k|φ~kvir⟩

Character concentrated in a few MOs

χtarget is the chosen local target function; T contains the dominant projected orbitals. The full virtual space is preserved.
02 / USER INTERFACE

PySCF preparation.
MPPT calculations.

A custom YAML interface wraps PySCF for SCF and CASSCF, with AVAS initial guesses and reusable orbital checkpoints.

MPPT’s own Fortran and C++ programs perform determinant selection, correlation and spin–orbit calculations. Run the stages separately through a shared HDF5 file or enable the complete workflow in one input.

ce.yaml / input and basis from the quick start

# Prepare the reference; stop before MPPT
pyscf2mppt ce.yaml \
  --output-dir runs/ce-reference \
  --stop-after dump

# Read saved orbital checkpoints
pyscf2mppt --inspect-checkpoint \
  runs/ce-reference/checkpoints
03 / YOUR FIRST CALCULATION

From SCF to
spin–orbit energy levels

The Ce³⁺ example takes you from reference orbitals to correlated, relativistic spin–orbit-split levels. Run preparation and downstream stages separately, then use one input for the complete calculation.

FULL RUN / WITH ce.yaml FROM THE GUIDE
pyscf2mppt ce.yaml --output-dir runs/ce-full

First result: runs/ce-full/summary.md

04 / THEORETICAL FOUNDATIONS

Diagrammatic formulation
of MPPT

The diagrammatic formulation connects general model spaces, determinant-specific partitionings and a second-order effective Hamiltonian. The literature guide follows the method from its foundations to relativistic applications.

Equations and literature
2026

A Quasirelativistic Perturbation Theory As a Tool for Modeling f–f and f–d Transitions in Lanthanide Ions

M. M. Seregin, I. O. Glebov, N. S. Mosyagin,
A. V. Oleynichenko & A. V. Zaitsevskii

THESIS

Mikhail Seregin · diploma thesis

Manuscript link to be added.

Placeholder
1995

Multi-partitioning quasidegenerate perturbation theory

A. Zaitsevskii & J.-P. Malrieu

1999

Diagrammatic formulation of the second-order many-body multipartitioning perturbation theory

A. Zaitsevskii & R. Cimiraglia

MPPT / Looking ahead

Roadmap

Current capabilities and next directions

  • Faster diagrams, less memory

    Selected DF CVCV factor streaming is available. Further diagram acceleration remains in development.

  • Magnetic properties

    Native moments, applicable principal g and field/thermal response are available. Model and controls ↗

  • Energy-guided selection

    New determinant-selection algorithms with feedback from MPPT energies.

  • Energy gradients