Investigation of Multiple-Bond Dissociation Using Brillouin–Wigner Perturbation with Improved Virtual Orbitals

Repartitioned Brillouin-Wigner perturbation theory with a size-consistent second-order correlation energy

Second-order Møller-Plesset perturbation theory (MP2) often breaks down catastrophically in small-gap systems, leaving much to be desired in its performance for myriad chemical applications such as noncovalent interactions, thermochemistry, and dative bonding in transition metal complexes. This divergence problem has reignited interest in Brillouin-Wigner perturbation theory (BWPT), which is regular at all orders but lacks size consistency and extensivity, severely limiting its application to chemistry. In this work, we propose an alternative partitioning of the Hamiltonian that leads to a regular BWPT perturbation series that, through the second order, is size-extensive, size-consistent (provided its Hartree-Fock reference is also), and orbital invariant. Our second-order size-consistent Brillouin-Wigner (BW-s2) approach can describe the exact dissociation limit of H2 in a minimal basis set, regardless of the spin polarization of the reference orbitals. More broadly, we find that BW-s2 offers improvements relative to MP2 for covalent bond breaking, noncovalent interaction energies, and metal/organic reaction energies, although rivaling coupled-cluster with single and double substitutions for thermochemical properties.

Block Effective Hamiltonian Theory and Its Application

Block effective Hamiltonian theory (BEHT) is presented in this work. Configuration interaction functions are divided into P, Q, and R spaces. Effective Hamiltonian is constructed with the partitioning technique within the P space. The eigenvalue problem of the effective Hamiltonian is then solved iteratively. It is demonstrated that the ground-state energies of N₂, HF, and F₂ calculated with BEHT converge to the multireference configuration interaction energies from below and the iteration number becomes smaller as BEHT energy becomes closer to the exact energy. The accuracy of BEHT is better than that of the second-order multireference perturbation theory, although the matrix elements in both methods are the same. The ionization potentials of the singlet state of HF, the doublet state of F, and the triplet state of NH and the potential energy curves of CH₄, C₂, and N₂ are calculated with BEHT and compared with experiments and results of CASSCF, CCSD, and CCSD(T) and the results of the full configuration interaction if available. The iteration numbers are all less than 10 in this study. These calculations show the good performances of BEHT in comparison with other theoretical approximation methods.

Single-Root Multireference Brillouin-Wigner Perturbative Approach to Excitation Energies

The state-specific Brillouin-Wigner multireference perturbation theory [which employs Jeziorski-Monkhorst parametrization of the wave function] using improved virtual orbitals, denoted as IVO-BWMRPT, is applied to calculate excitation energies (EEs) for methylene, ethylene, trimethylenemethane, and benzyne systems exhibiting various degrees of diradical character. In IVO-BWMRPT, all of the parameters appearing in the wave function ansatz are optimized for a specific electronic state. For these systems, the IVO-BWMRPT method provides EEs that are in close agreement with the benchmark results and experiments, where available, indicating that the method does not introduce imbalance in the target-specific treatment of closed- and open-shell states involved. The good performance of the present methodology is primarily related to structural compactness of the formalism. Overall, present findings are encouraging for both further development of the approach and chemical applications on the energy differences of strongly correlated systems.