The Facility for Rare Isotope Beams (FRIB) will be a world-leading laboratory for the study of nuclear structure, reactions and astrophysics. Experiments with intense beams of rare isotopes produced at FRIB will guide us toward a comprehensive description of nuclei, elucidate the origin of the elements in the cosmos, help provide an understanding of matter in neutron stars, and establish the scientific foundation for innovative applications of nuclear science to society. FRIB will be essential for gaining access to key regions of the nuclear chart, where the measured nuclear properties will challenge established concepts, and highlight shortcomings and needed modifications to current theory. Conversely, nuclear theory will play a critical role in providing the intellectual framework for the science at FRIB, and will provide invaluable guidance to FRIB's experimental programs. This article overviews the broad scope of the FRIB theory effort, which reaches beyond the traditional fields of nuclear structure and reactions, and nuclear astrophysics, to explore exciting interdisciplinary boundaries with other areas. \keywords{Nuclear Structure and Reactions. Nuclear Astrophysics. Fundamental Interactions. High Performance Computing. Rare Isotopes. Radioactive Beams.

We discuss several pairing-related phenomena in nuclear systems, ranging from superfluidity in neutron stars to the gradual breaking of pairs in finite nuclei. We focus on the links between many-body pairing as it evolves from the underlying nucleon-nucleon interaction and the eventual experimental and theoretical manifestations of superfluidity in infinite nuclear matter and of pairing in finite nuclei. We analyse the nature of pair correlations in nuclei and their potential impact on nuclear structure experiments. We also describe recent experimental evidence that points to a relation between pairing and phase transitions (or transformations) in finite nuclear systems. Finally, we discuss recent investigations of ground-state properties of random two-body interactions where pairing plays little role although the interactions yield interesting nuclear properties such as 0+ ground states in even-even nuclei.

Using the most recent realistic effective interactions for nuclear matter with a smooth extrapolation to high densities including causality, we constrain the equation of state and calculate maximum masses of rotating neutron stars. First and second order phase transitions to, e.g., quark matter at high densities are included. If neutron star masses of $\sim 2.3M_\odot$ from quasi-periodic oscillations in low mass X-ray binaries are confirmed, a soft equation of state as well as strong phase transitions can be excluded in neutron star cores.

Recent equations of state for dense nuclear matter are discussed with possible phase transitions arising in neutron stars such as pion, kaon and hyperon kondensation, superfluidity and quark matter. Specifically, we treat the nuclear to quark matter phase transition, the possible mixed phase and its structure. A number of numerical calculations of rotating neutron stars with and without phase transitions are given and compared to observed masses, radii, temperatures and glitches.

As rotating neutron stars slow down, the pressure and the density in the core region increase due to the decreasing centrifugal forces and phase transitions may occur in the center. We extract the analytic behavior near the critical angular velocity $\Omega_0$, where the phase transitions occur in the center of a neutron star, and calculate the moment of inertia, angular velocity, rate of slow down, braking index, etc. For a first order phase transition these quantities have a characteristic behavior, e.g., the braking index diverges as $\sim(\Omega_0-\Omega)^{-1/2}$. Observational consequences for first, second and other phase transitions are discussed.