Double-Core Hole (DCH) states of small molecules are assessed with the restricted<br>active space self-consistent field (RASSCF) and multi-state restricted active space perturbation<br>theory of second order (MS-RASPT2) approximations. To ensure an unbiased<br>description of the relaxation and correlation effects on the DCH states, the neutral<br>ground state and DCH wave functions are optimized separately, whereas the spectral<br>intensities are computed with a biorthonormalized set of molecular orbitals within the<br>state-interaction (SI) approximation. Accurate shake-up satellites binding energies and<br>intensities of double-core-ionized states (K<sup>-2</sup>) are obtained for H<sub>2</sub>O, N<sub>2</sub>, CO and C<sub>2</sub>H<sub>2n</sub><br>(n=1–3). The results are analyzed in details and show excellent agreement with recent<br>experimental data.