The combination of multipartitioning of the Hamiltonian with canonical Van Vleck perturbation theory leads to a Hermitian variant of quasidegenerate N-electron valence perturbation theory.

Many recent developments in the area of multistate multireference perturbation theories focused on methods that use a state-averaged 0th order Hamiltonian. We recently found that the dynamic correlation dressed complete active space method fails in describing ligand field and charge transfer states in a balanced way precisely because it uses a state-averaged 0th order Hamiltonian [L. Lang and F. Neese, J. Chem. Phys. 150, 104104 (2019)]. The multipartitioning idea allows the use of state-specific 0th order Hamiltonians in a multistate framework and could therefore alleviate the mentioned problem. However, the effective Hamiltonian is non-Hermitian in the traditional formulation of multipartitioning, which can lead to unphysical behavior, especially for nearly degenerate states. In order to achieve a more balanced treatment of states with different physical character and at the same time have a Hermitian effective Hamiltonian, we combine in this work multipartitioning with canonical Van Vleck perturbation theory. At the 2nd order, the result is a Hermitian variant of multipartitioning quasidegenerate N-electron valence state perturbation theory. The effect of model space noninvariance of the method is discussed and the benefit of a Hermitian formulation is highlighted with numerical examples. The method is shown to give good results for the calculation of electronic transitions of the [CuCl4]2 -complex and for the calculation of electron paramagnetic resonance parameters, which are two examples where the balance between ligand field and charge transfer configurations is of utmost importance.