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 nondegenerate single-particle levels and including isospin symmetry-breaking terms. The power of the method is illustrated with a numerical calculation for 64 Ge for a pf g 9=2 model space which is out of reach of modern shell-model codes. DOI: 10.1103/PhysRevLett.96.072503 PACS numbers: 21.60.Fw, 03.65.ÿw, 27.50.+e, 74.20.Rp Exactly solvable models (ESM) provide important insights into the structure of many-body quantum systems. The two main advantages of ESMs are: (1) they can describe in an analytical or exact numerical way a wide variety of elementary phenomena. (2) They can be and have been used as a testing ground for various many-body approaches.A particular class of ESMs, extensively used in nuclear physics, are the dynamical-symmetry models. In this case the Hamiltonian can be expressed in terms of Casimir operators of a chain of nested algebras. An example often used to introduce nuclear superconductivity [see, e.g., Ref.[1] ] is the rank-1 (Lie) algebra SU(2). Examples of dynamical-symmetry models associated with a rank-2 algebra are Elliott's SU(3) model of nuclear deformation [2] and the SO(5) model of T 1 isovector pairing between neutrons and protons [3] which has found many applications in nuclei [see, e.g., Ref. [4] ].The concept of quantum integrability, closely linked with exact solvability, goes beyond the limits of the dynamical-symmetry approach. A quantum system is integrable if there exist as many commuting Hermitian operators (integrals of motion) as quantum degrees of freedom [5]. The set of Casimir operators of a chain of nested algebras satisfies this condition.Dynamical-symmetry models are usually defined for degenerate single-particle levels. Lifting this degeneracy breaks the dynamical symmetry but may still preserve integrability. The pairing model with nondegenerate single-particle levels, of which an exact solution was found by Richardson in the 1960s [6], represents an example of an ESM with such characteristics. Recently, more general exactly solvable pairing models, both for fermions and for bosons, called Richardson-Gaudin (RG) models, have been proposed [7,8].The RG pairing models are based on rank-1 algebras: SU(2) for fermions and SU(1,1) for bosons. In this Letter we carry out the first step in extending the RG models to higher-rank algebras by considering a RG model based on the rank-2 algebra SO(5). The model Hamiltonian describes a two-component system consisting of neutrons and protons interacting through an isovector (T 1) pairing force and distributed over nondegenerate orbits. This neutron-proton (np) pairing Hamiltonian with nondegenerate orbits has been studied by Richardson [9] who proposed an exact solution. However, it was shown subsequently that Richardson's solution is incorrect for more than two nucleon pairs [10] by explicitly solving the case of three-nucleon pairs. Independently, Links et al. derived an ex...