Exact Solution of the Isovector Proton Neutron Pairing Hamiltonian
The complete exact solution of the T=1 neutron-proton pairing Hamiltonian is presented in the context of the SO(5) Richardson-Gaudin model with non-degenerate single-particle levels and including isospin-symmetry breaking terms. The power of the method is illustrated with a numerical calculation for Ge for a model space which is out of reach of modern shell-model codes.