Is the secondary component of GW190814 the lightest black hole or the heaviest neutron star ever discovered in a double compact-object system [R. Abbott et al., ApJ Lett., 896, L44 (2020)]? This is the central question animating this letter. Covariant density functional theory provides a unique framework to investigate both the properties of finite nuclei and neutron stars, while enforcing causality at all densities. By tuning existing energy density functionals we were able to: (a) account for a 2.6 Msun neutron star, (b) satisfy the original constraint on the tidal deformability of a 1.4 Msun neutron star, and (c) reproduce ground-state properties of finite nuclei. Yet, for the class of models explored in this work, we find that the stiffening of the equation of state required to support super-massive neutron stars is inconsistent with either constraints obtained from energetic heavy-ion collisions or from the low deformability of medium-mass stars. Thus, we speculate that the maximum neutron star mass can not be significantly higher than the existing observational limit and that the 2.6 Msun compact object is likely to be the lightest black hole ever discovered.

Uncertainty quantification has emerged as a rapidly growing field in nuclear science. Theoretical predictions of physical observables often involve extrapolations to regions that are poorly constrained by laboratory experiments and astrophysical observations. Without properly quantified theoretical errors, such model predictions are of limited value. Also, one often deals with theoretical constructs that involve fundamental quantities that are not accessible to experiment or observation. Particularly relevant in this context is the pressure of pure neutron matter. In this contribution we develop an analytic framework to answer the question of "How can new data reduce uncertainties of current theoretical models?" [P.-G. Reinhard and W. Nazarewicz, Phys. Rev. C81, 051303(R) (2010)]. Simple and insightful expressions are obtained to quantify the impact of one or two new observables on theoretical uncertainties in two critical quantities: the slope of the symmetry energy at saturation density and the pressure of pure neutron matter at twice nuclear matter saturation density.

How does subatomic matter organize itself? Neutron stars are cosmic laboratories uniquely poised to answer this fundamental question that lies at the heart of nuclear science. Newly commissioned rare isotope facilities, telescopes operating across the entire electromagnetic spectrum, and ever more sensitive gravitational wave detectors will probe the properties of neutron-rich matter with unprecedented precision over an enormous range of densities. Yet, a coordinated effort between observation, experiment, and theoretical research is of paramount importance for realizing the full potential of these investments. Theoretical nuclear physics provides valuable insights into the properties of neutron-rich matter in regimes that are not presently accessible to experiment or observation. In particular, nuclear density functional theory is likely the only tractable framework that can bridge the entire nuclear landscape by connecting finite nuclei to neutron stars. This compelling connection is the main scope of the present review.

Despite a length-scale difference of 18 orders of magnitude, the internal structure of neutron stars and the spatial distribution of neutrons in atomic nuclei are profoundly connected.

The complex nature of the nuclear forces generates a broad range and diversity of observational phenomena. Heavy nuclei, though orders of magnitude less massive than neutron stars, are governed by the same underlying physics, which is enshrined in the nuclear equation of state. Heavy nuclei are expected to develop a neutron-rich skin where many neutrons collect near the surface. Such a skin thickness is strongly sensitive to the poorly-known density dependence of the symmetry energy near saturation density. An accurate and model-independent determination of the neutron-skin thickness of heavy nuclei would provide a significant first constraint on the density dependence of the nuclear symmetry energy. The determination of the neutron-skin thickness of heavy nuclei has far reaching consequences in many areas of physics as diverse as heavy-ion collisions, polarized electron and proton scattering off nuclei, precision tests of the standard model using atomic parity violation, and nuclear astrophysics. While a systematic and concerted experimental effort has been made to measure the neutron-skin thickness of heavy nuclei, a precise and model-independent determination remains elusive. How to move forward at a time when many new facilities are being commissioned and how to strengthen the synergy with other areas of physics are primary goals of this review.

The first detection of a binary neutron star merger has opened the brand new era of multimessenger astronomy. This historic detection has been instrumental in providing constraints on the tidal polarizability of neutron stars. In turn, the tidal polarizability has been used to impose limits on stellar radii and ultimately on the equation of state (EOS). The tidal polarizability is also sensitive to the second tidal Love number k2. It is the main purpose of this work to perform a detailed study of k2 which, for a given compactness parameter, encodes the entire sensitivity of the tidal polarizability to the EOS. In particular, we examine the role that the crustal component of the EOS plays in the determination of k2. A set of realistic models of the equation of state that yield an accurate description of the properties of finite nuclei and support neutron stars of two solar masses is used. Given that the tidal polarizability scales as the fifth power of the compactness parameter, a universal relation exists among the tidal polarizability and the compactness parameter that is highly insensitive to the underlying EOS. Thus, besides an extraction of the tidal polarizabilities, a measurement of the individual stellar masses is also required to impact the mass-radius relation. However, we observe a strong sensitivity of k2 to the EOS, particularly to the contribution from the inner crust. Although by itself the tidal polarizability can not contribute to the determination of the mass-radius relation, future detections of binary neutron star mergers are poised to provide significant constraints on both the tidal polarizabilities and masses of the individual stars, and thus ultimately on the mass-radius relation. Yet, subleading corrections to the tidal polarizability are encoded in the second Love number k2 which displays a large sensitivity to the entire (crust-plus-core) EOS.

On August 17, 2017 the LIGO-Virgo collaboration detected for the first time gravitational waves from the binary merger of two neutron stars (GW170817). Unlike the merger of two black holes, the associated electromagnetic radiation was also detected by a host of telescopes operating over a wide range of frequencies---opening a brand new era of multimessenger astronomy. This historical detection is providing fundamental new insights into the astrophysical site for the r-process and on the nature of dense matter. In this contribution we examine the impact of GW170817 on the equation of state of neutron rich matter, particularly on the density dependence of the symmetry energy. Limits on the tidal polarizability extracted from GW170817 seem to suggest that the symmetry energy is soft, thereby excluding models that predict overly large stellar radii.

Neutron stars are unique cosmic laboratories for the exploration of matter under extreme conditions of density and neutron-proton asymmetry. Due to their enormous dynamic range, neutron stars display a myriad of exotic states of matter that are impossible to recreate under normal laboratory conditions. In these three lectures I will discuss how the strong synergy that has developed between nuclear physics and astrophysics will uncover some of the deepest secrets behind these fascinating objects. In particular, I will highlight the enormous impact that the very first detection of gravitational waves from the binary neutron-star merger GW170817 is having in constraining the composition, structure, and dynamics of neutron stars.

The historical first detection of a binary neutron star merger by the LIGO-Virgo collaboration [B. P. Abbott et al. Phys. Rev. Lett. 119, 161101 (2017)] is providing fundamental new insights into the astrophysical site for the $r$-process and on the nature of dense matter. A set of realistic models of the equation of state (EOS) that yield an accurate description of the properties of finite nuclei, support neutron stars of two solar masses, and provide a Lorentz covariant extrapolation to dense matter are used to confront its predictions against tidal polarizabilities extracted from the gravitational-wave data. Given the sensitivity of the gravitational-wave signal to the underlying EOS, limits on the tidal polarizability inferred from the observation translate into constraints on the neutron-star radius. Based on these constraints, models that predict a stiff symmetry energy, and thus large stellar radii, can be ruled out. Indeed, we deduce an upper limit on the radius of a $1.4\,M_{\odot}$ neutron star of $R_{\star}^{1.4}\!<\!13.76\,{\rm km}$. Given the sensitivity of the neutron-skin thickness of ${}^{208}$Pb to the symmetry energy, albeit at a lower density, we infer a corresponding upper limit of about $R_{\rm skin}^{208}\!\lesssim\!0.25\,{\rm fm}$. However, if the upcoming PREX-II experiment measures a significantly thicker skin, this may be evidence of a softening of the symmetry energy at high densities---likely indicative of a phase transition in the interior of neutron stars.

It has been recently suggested that differences in the charge radii of mirror nuclei are proportional to the neutron-skin thickness of neutron-rich nuclei and to the slope of the symmetry energy $L$ [B.A. Brown Phys. Rev. Lett. 119, 122502 (2017)]. The determination of the neutron skin has important implications for nuclear physics and astrophysics. Although the use of electroweak probes provides a largely model-independent determination of the neutron skin, the experimental challenges are enormous. Thus, the possibility that differences in the charge radii of mirror nuclei may be used as a surrogate for the neutron skin is a welcome alternative. To test the validity of this assumption we perform calculations based on a set of relativistic energy density functionals that span a wide region of values of $L$. Our results confirm that the difference in charge radii of various neutron-deficient nickel isotopes and their corresponding mirror nuclei is indeed strongly correlated to both the neutron-skin thickness and $L$. Moreover, given that various neutron-star properties are also sensitive to $L$, a data-to-data relation emerges between the difference in charge radii of mirror nuclei and the radius of low-mass neutron stars.

Exotic nuclei, particularly those near the driplines, are at the core of one of the fundamental questions driving nuclear structure and astrophysics today: what are the limits of nuclear binding? Exotic nuclei play a critical role in both informing theoretical models as well as in our understanding of the origin of the heavy elements. Our purpose is to refine existing mass models through the training of an artificial neural network that will mitigate the large model discrepancies far away from stability. The basic paradigm of our two-pronged approach is an existing mass model that captures as much as possible of the underlying physics followed by the implementation of a Bayesian Neural Network (BNN) refinement to account for the missing physics. Bayesian inference is employed to determine the parameters of the neural network so that model predictions may be accompanied by theoretical uncertainties. Despite the undeniable quality of the mass models adopted in this work, we observe a significant improvement (of about 40%) after the BNN refinement is implemented. Indeed, in the specific case of the Duflo-Zuker mass formula, we find that the rms deviation relative to experiment is reduced from rms =0.503MeV to rms=0.286 MeV. These newly refined mass tables are used to map the neutron drip lines (or rather "drip bands") and to study a few critical r-process nuclei. The BNN approach is highly successful in refining the predictions of existing mass models. In particular, the large discrepancy displayed by the original "bare" models in regions where experimental data is unavailable is considerably quenched after the BNN refinement. This lends credence to our approach and has motivated us to publish refined mass tables that we trust will be helpful for future astrophysical applications.

Besides their intrinsic nuclear-structure value, nuclear mass models are essential for astrophysical applications, such as r-process nucleosynthesis and neutron-star structure. To overcome the intrinsic limitations of existing "state-of-the-art" mass models, we propose a refinement based on a Bayesian Neural Network (BNN) formalism. A novel BNN approach is implemented with the goal of optimizing mass residuals between theory and experiment. A significant improvement (of about 40%) in the mass predictions of existing models is obtained after BNN refinement. Moreover, these improved results are now accompanied by proper statistical errors. Finally, by constructing a "world average" of these predictions, a mass model is obtained that is used to predict the composition of the outer crust of a neutron star. The power of the Bayesian neural network method has been successfully demonstrated by a systematic improvement in the accuracy of the predictions of nuclear masses. Extension to other nuclear observables is a natural next step that is currently under investigation.

Recent progress in the determination of both masses and radii of neutron stars are starting to place stringent constraints on the dense matter equation of state. In particular, new theoretical developments together with improved statistical tools seem to favor stellar radii that are significantly smaller than those predicted by models using purely nucleonic equations of state. Given that the underlying equation of state must also account for the observation of $2M_{\odot}$ neutron stars, theoretical approaches to the study of the dense matter equation of state are facing serious challenges. In response to this challenge, we compute in a model-independent way the underlying equation of state associated with an assumed mass-radius template similar to the "common radius" assumption used in recent studies. Once such a mass-radius template is adopted, the equation of state follows directly from the implementation of Lindblom's algorithm; assumptions on the nature or composition of the dense stellar core are not required. By analyzing mass-radius profiles with a maximum mass consistent with observation and common radii in the 8 to 11 km range, a lower limit on the stellar radius of a $1.4M_{\odot}$ neutron star of $R_{\rm NS}\!\gtrsim\!10.7$ km is obtained in order to prevent the equation of state from violating causality.

The main goal of the present contribution is a pedagogical introduction to the fascinating world of neutron stars by relying on relativistic density functional theory. Density functional theory provides a powerful--and perhaps unique--framework for the calculation of both the properties of finite nuclei and neutron stars. Given the enormous densities that may be reached in the core of neutron stars, it is essential that such theoretical framework incorporates from the outset the basic principles of Lorentz covariance and special relativity. After a brief historical perspective, we present the necessary details required to compute the equation of state of dense, neutron-rich matter. As the equation of state is all that is needed to compute the structure of neutron stars, we discuss how nuclear physics--particularly certain kind of laboratory experiments--can provide significant constrains on the behavior of neutron-rich matter.

Theoretical predictions of physical observables often involve extrapolations to regions that are poorly constrained by laboratory experiments and astrophysical observations. Without properly quantified theoretical errors, such model predictions are of very limited utility. In this contribution we use maximum-likelihood estimation to compute theoretical errors and assess correlations between physical observables. We illustrate the power and elegance of these methods using examples of both pedagogical and realistic interest. In particular, we implement a gaussian approximation to the likelihood function to develop a new relativistic effective interaction constrained by ground-state properties of finite nuclei, their monopole response, and masses of neutron stars.

Pulsar glitches-the sudden spin-up in the rotational frequency of a neutron star-suggest the existence of an angular-momentum reservoir confined to the inner crust of the neutron star. Large and regular glitches observed in the Vela pulsar have originally constrained the fraction of the stellar moment of inertia that must reside in the solid crust to about 1.4%. However, crustal entrainment-which until very recently has been ignored-suggests that in order to account for the Vela glitches, the fraction of the moment of inertia residing in the crust must increase to about 7%. This indicates that the required angular momentum reservoir may exceed that which is available in the crust. We explore the possibility that uncertainties in the equation of state provide enough flexibility for the construction of models that predict a large crustal thickness and consequently a large crustal moment of inertia. Given that analytic results suggest that the crustal moment of inertia is sensitive to the transition pressure at the crust-core interface, we tune the parameters of the model to maximize the transition pressure, while still providing an excellent description of nuclear observables. In this manner we are able to obtain fractional moments of inertia as large as 7% for neutron stars with masses below 1.6 solar masses. In particular, we find that if the neutron-skin thickness of 208Pb falls within the (0.20-0.26) fm range, large enough transition pressures can be generated to explain the large Vela glitches without invoking an additional angular-momentum reservoir beyond that confined to the solid crust. Our results suggest that the crust may be enough.

The symmetry energy describes how the energy of nuclear matter rises as one goes away from equal numbers of neutrons and protons. This is very important to describe neutron rich matter in astrophysics. This article reviews our knowledge of the symmetry energy from theoretical calculations, nuclear structure measurements, heavy ion collisions, and astronomical observations. We then present a roadmap to make progress in areas of relevance to the symmetry energy that promotes collaboration between the astrophysics and the nuclear physics communities.

We explore the unique and fascinating structure of neutron stars. Although neutron stars are of interest in many areas of Physics, our aim is to provide an intellectual bridge between Nuclear Physics and Astrophysics. We argue against the naive perception of a neutron star as a uniform assembly of neutrons packed to enormous densities. Rather, by focusing on the many exotic phases that are speculated to exist in a neutron star, we show how the reality is different and far more interesting.

The Lead Radius Experiment (PREX) has provided the first model-independent evidence in favor of a neutron-rich skin in ${}^{208}$Pb. Although the error bars are large, the reported large central value of 0.33\,fm is particularly intriguing. To test whether such a thick neutron-skin in ${}^{208}$Pb is already incompatible with laboratory experiments or astrophysical observations, we employ relativistic models with neutron-skin thickness in ${}^{208}$Pb ranging from 0.16 to 0.33 fm to compute ground state properties of finite nuclei, their collective monopole and dipole response, and mass-{\sl vs}-radius relations for neutron stars. No compelling reason was found to rule out models with large neutron skins in ${}^{208}$Pb from the set of observables considered in this work.

The neutron-skin thickness of heavy nuclei provides a fundamental link to the equation of state of neutron-rich matter, and hence to the properties of neutron stars. The Lead Radius Experiment ("PREX") at Jefferson Laboratory has recently provided the first model-independence evidence on the existence of a neutron-rich skin in 208Pb. In this contribution we examine how the increased accuracy in the determination of neutron skins expected from the commissioning of intense polarized electron beams may impact the physics of neutron stars.

Laboratory experiments sensitive to the density dependence of the symmetry energy may place stringent constraints on the equation of state of neutron-rich matter and, thus, on the structure, dynamics, and composition of neutron stars. Understanding the equation of state of neutron-rich matter is a central goal of nuclear physics that cuts across a variety of disciplines. In this contribution I focus on how laboratory experiments on neutron skins and on both Pygmy and Giant resonances can help us elucidate the structure of neutron stars.

Background: The inner crust of a neutron star is believed to consist of Coulomb-frustrated complex structures known as "nuclear pasta" that display interesting and unique low-energy dynamics. Purpose: To elucidate the structure and composition of the neutron-star crust as a function of temperature, density, and proton fraction. Methods: A new lattice-gas model, the "Charged-Ising Model" (CIM), is introduced to simulate the behavior of neutron-star matter. Preliminary Monte Carlo simulations on 30^3 lattices are performed for a variety of temperatures, densities, and proton fractions. Results: Results are obtained for the heat capacity, pair-correlation function, and static structure factor for a variety of conditions appropriate to the inner stellar crust. Conclusions: Although relatively simple, the CIM captures the essence of Coulomb frustration that is required to simulate the subtle dynamics of the inner stellar crust. Moreover, the computationally demanding long-range Coulomb interactions have been pre-computed at the appropriate lattice sites prior to the start of the simulation resulting in enormous computational gains. This work demonstrates the feasibility of future CIM simulations involving a large number of particles as a function of density, temperature, and proton fraction.

Understanding the equation of state (EOS) of neutron-rich matter is a central goal of nuclear physics that cuts across a variety of disciplines. Indeed, the limits of nuclear existence, the collision of energetic heavy ions, the structure of neutron stars, and the dynamics of core-collapse supernova all depend critically on the nuclear-matter EOS. In this contribution I focus on the EOS of cold baryonic matter with special emphasis on its impact on the structure, dynamics, and composition of neutron stars. In particular, I discuss how laboratory experiments on neutron skins as well as on Pygmy and Giant resonances can help us elucidate the structure of these fascinating objects.

Background: The neutron skin of a heavy nucleus as well as many neutron-star properties are highly sensitive to the poorly constrained density dependence of the symmetry energy. Purpose: To provide for the first time meaningful theoretical errors and to assess the degree of correlation between the neutron-skin thickness of ${}^{208}$Pb and several neutron-star properties. Methods: A proper covariance analysis based on the predictions of an accurately-calibrated relativistic functional {\sl "FSUGold"} is used to quantify theoretical errors and correlation coefficients. Results: We find correlation coefficients of nearly one (or minus one) between the neutron-skin thickness of ${}^{208}$Pb and a host of observables of relevance to the structure, dynamics, and composition of neutron stars. Conclusions: We suggest that a follow-up PREX measurement, ideally with a 0.5% accuracy, could significantly constrain the equation of state of neutron-star matter.

The recent high-resolution measurement of the electric dipole (E1) polarizability (alphad) in 208Pb [Phys. Rev. Lett. 107, 062502 (2011)] provides a unique constraint on the neutron-skin thickness of this nucleus. The neutron-skin thickness (rskin) of 208Pb is a quantity of critical importance for our understanding of a variety of nuclear and astrophysical phenomena. To assess the model dependence of the correlation between alphad and rskin, we carry out systematic calculations for 208Pb, 132Sn, and 48Ca based on the nuclear density functional theory (DFT) using both non-relativistic and relativistic energy density functionals (EDFs). Our analysis indicates that whereas individual models exhibit a linear dependence between alphad and rskin, this correlation is not universal when one combines predictions from a host of different models. By averaging over these model predictions, we provide estimates with associated systematic errors for rskin and alphad for the nuclei under consideration. We conclude that precise measurements of rskin in both 48Ca and 208Pb---combined with the recent measurement of alphad---should significantly constrain the isovector sector of the nuclear energy density functional.

We review and extend with nonrelativistic nuclear mean field calculations a previous study of the impact of the nuclear symmetry energy on the structure and composition of the outer crust of nonaccreting neutron stars. First, we develop a simple "toy model" to understand the most relevant quantities determining the structure and composition of the outer crust: the nuclear symmetry energy and the pressure of the electron gas. While the latter is a well determined quantity, the former ---specially its density dependence--- still lacks an accurate characterization. We thus focus on the influence of the nuclear symmetry energy on the crustal composition. For that, we employ different nuclear models that are accurate in the description of terrestrial nuclei. We show that those models with stiffer symmetry energies ---namely, those that generate thicker neutron skins in heavy nuclei and have smaller symmetry energies at subnormal nuclear densities--- generate more exotic isotopes in the stellar crust than their softer counterparts.

Monte Carlo simulations of neutron-rich matter of relevance to the inner neutron-star crust are performed for a system of A=5,000 nucleons. To determine the proton fraction in the inner crust, numerical simulations are carried out for a variety of densities and proton fractions. We conclude---as others have before us using different techniques---that the proton fraction in the inner stellar crust is very small. Given that the purported "nuclear pasta" phase in stellar crusts develops as a consequence of the long-range Coulomb interaction among protons, we question whether pasta formation is possible in such proton-poor environments. To answer this question, we search for physical observables sensitive to the transition between spherical nuclei and exotic pasta structures. Of particular relevance is the static structure factor S(k)---an observable sensitive to density fluctuations. However, no dramatic behavior was observed in S(k). We regard the identification of physical observables sensitive to the existence---or lack-thereof---of a pasta phase in proton-poor environments as an open problem of critical importance.

Nuclear effective interactions are useful tools in astrophysical applications especially if one can guide the extrapolations to the extremes regions of isospin and density that are required to simulate dense, neutron-rich systems. Isospin extrapolations may be constrained in the laboratory by measuring the neutron skin thickness of a heavy nucleus, such as 208Pb. Similarly, future observations of massive neutron stars will constrain the extrapolations to the high-density domain. In this contribution we introduce a new relativistic effective interaction that is simultaneously constrained by the properties of finite nuclei, their collective excitations, and neutron-star properties. By adjusting two of the empirical parameters of the theory, one can efficiently tune the neutron skin thickness of 208Pb and the maximum neutron star mass. We illustrate this procedure in response to the recent interpretation of X-ray observations by Steiner, Lattimer, and Brown that suggests that the FSUGold effective interaction predicts neutron star radii that are too large and a maximum stellar mass that is too small. The new effective interaction is fitted to a neutron skin thickness in 208Pb of only R_n - R_p = 0.16 fm and a moderately large maximum neutron star mass of 1.94 Msun.

The sensitivity of the stellar moment of inertia to the neutron-star matter equation of state is examined using accurately-calibrated relativistic mean-field models. We probe this sensitivity by tuning both the density dependence of the symmetry energy and the high density component of the equation of state, properties that are at present poorly constrained by existing laboratory data. Particularly attractive is the study of the fraction of the moment of inertia contained in the solid crust. Analytic treatments of the crustal moment of inertia reveal a high sensitivity to the transition pressure at the core-crust interface. This may suggest the existence of a strong correlation between the density dependence of the symmetry energy and the crustal moment of inertia. However, no correlation was found. We conclude that constraining the density dependence of the symmetry energy - through, for example, the measurement of the neutron skin thickness in 208Pb - will place no significant bound on either the transition pressure or the crustal moment of inertia.

Motivated by a recent astrophysical measurement of the pressure of cold matter above nuclear-matter saturation density, we compute the equation of state of neutron star matter using accurately calibrated relativistic models. The uniform stellar core is assumed to consist of nucleons and leptons in beta equilibrium; no exotic degrees of freedom are included. We found the predictions of these models to be in fairly good agreement with the measured equation of state. Yet the Mass-vs-Radius relations predicted by these same models display radii that are consistently larger than the observations.

Understanding the equation of state (EOS) of cold nuclear matter, namely, the relation between the pressure and energy density, is a central goal of nuclear physics that cuts across a variety of disciplines. Indeed, the limits of nuclear existence, the collision of heavy ions, the structure of neutron stars, and the dynamics of core-collapse supernova, all depend critically on the equation of state of hadronic matter. In this contribution I will concentrate on the special role that nuclear physics plays in constraining the EOS of cold baryonic matter and its impact on the properties of neutron stars.

The composition and equation of state of the outer crust of non-accreting neutron stars is computed using accurate nuclear mass tables. The main goal of the present study is to understand the impact of the symmetry energy on the structure of the outer crust. First, a simple "toy model" is developed to illustrate the competition between the electronic density and the symmetry energy. Then, realistic mass tables are used to show that models with a stiff symmetry energy - those that generate large neutron skins for heavy nuclei - predict a sequence of nuclei that are more neutron-rich than their softer counterparts. This result may be phrased in the form of a correlation: the larger the neutron skin of 208Pb, the more exotic the composition of the outer crust.

Relativistic mean-field models of nuclear structure have been enormously successful at reproducing ground-state properties of finite nuclei throughout the periodic table using a handful of accurately calibrated parameters. In this contribution we use powerful theoretical, experimental, and observational constraints -- not employed in the calibration procedure -- to validate two such models: NL3 and FSUGold. It is observed that FSUGold is consistent with all these constraints, except perhaps for a high density equation of state that appears mildly softer than required by astronomical observations. It is argued that incorporating such constrains goes a long way in removing much of the ambiguity left over from the standard calibrating procedure.

The neutron radius of a heavy nucleus is a fundamental nuclear-structure observable that remains elusive. Progress in this arena has been limited by the exclusive use of hadronic probes that are hindered by large and controversial uncertainties in the reaction mechanism. The Parity Radius Experiment at the Jefferson Laboratory offers an attractive electro-weak alternative to the hadronic program and promises to measure the neutron radius of 208Pb accurately and model independently via parity-violating electron scattering. In this contribution we examine the far-reaching implications that such a determination will have in areas as diverse as nuclear structure, atomic parity violation, and astrophysics.

The nucleus of 208Pb -- a system that is 18 order of magnitudes smaller and 55 orders of magnitude lighter than a neutron star -- may be used as a miniature surrogate to establish important correlations between its neutron skin and several neutron-star properties. Indeed, a nearly model-independent correlation develops between the neutron skin of 208Pb and the liquid-to-solid transition density in a neutron star. Further, we illustrate how a measurement of the neutron skin in 208Pb may be used to place important constraints on the cooling mechanism operating in neutron stars and may help elucidate the existence of quarks stars.

The nuclear pasta -- a novel state of matter having nucleons arranged in a variety of complex shapes -- is expected to be found in the crust of neutron stars and in core-collapse supernovae at subnuclear densities of about $10^{14}$ g/cm$^3$. Due to frustration, a phenomenon that emerges from the competition between short-range nuclear attraction and long-range Coulomb repulsion, the nuclear pasta displays a preponderance of unique low-energy excitations. These excitations could have a strong impact on many transport properties, such as neutrino propagation through stellar environments. The excitation spectrum of the nuclear pasta is computed via a molecular-dynamics simulation involving up to 100,000 nucleons. The dynamic response of the pasta displays a classical plasma oscillation in the 1-2 MeV region. In addition, substantial strength is found at low energies. Yet this low-energy strength is missing from a simple ion model containing a single-representative heavy nucleus. The low-energy strength observed in the dynamic response of the pasta is likely to be a density wave involving the internal degrees of freedom of the clusters.

Based upon previous discussions on the structure of compact stars geared towards undergraduate physics students, a real experiment involving two upper-level undergraduate physics students, a beginning physics graduate, and two advanced graduate students was conducted. A recent addition to the physics curriculum at Florida State University, "The Physics of Stars", sparked quite a few students' interests in the subject matter involving stellar structure. This, coupled with "Stars and Statistical Physics" by Balian and Blaizot and "Neutron Stars for Undergraduates" by Silbar and Reddy, is the cornerstone of this small research group who tackled solving the structure equations for compact objects in the Summer of 2004. Through the use of a simple finite-difference algorithm coupled to Microsoft Excel and Maple, solutions to the equations for stellar structure are presented in the Newtonian regime appropriate to the physics of white dwarf stars.

An accurately calibrated relativistic parametrization is introduced to compute the ground state properties of finite nuclei, their linear response, and the structure of neutron stars. While similar in spirit to the successful NL3 parameter set, it produces an equation of state that is considerably softer -- both for symmetric nuclear matter and for the symmetry energy. This softening appears to be required for an accurate description of several collective modes having different neutron-to-proton ratios. Among the predictions of this model are a symmetric nuclear-matter incompressibility of K=230 MeV and a neutron skin thickness in 208Pb of Rn-Rp=0.21 fm. Further, the impact of such a softening on the properties of neutron stars is as follows: the model predicts a limiting neutron star mass of Mmax=1.72 Msun, a radius of R=12.66 km for a ``canonical'' M=1.4 Msun neutron star, and no (nucleon) direct Urca cooling in neutrons stars with masses below M=1.3 Msun.

Nonuniform neutron-rich matter present in both core-collapse supernovae and neutron-star crusts is described in terms of a semiclassical model that reproduces nuclear-matter properties and includes long-range Coulomb interactions. The neutron-neutron correlation function and the corresponding static structure factor are calculated from molecular dynamics simulations involving 40,000 to 100,000 nucleons. The static structure factor describes coherent neutrino scattering which is expected to dominate the neutrino opacity. At low momentum transfers the static structure factor is found to be small because of ion screening. In contrast, at intermediate momentum transfers the static structure factor displays a large peak due to coherent scattering from all the neutrons in a cluster. This peak moves to higher momentum transfers and decreases in amplitude as the density increases. A large static structure factor at zero momentum transfer, indicative of large density fluctuations during a first-order phase transition, may increase the neutrino opacity. However, no evidence of such an increase has been found. Therefore, it is unlikely that the system undergoes a simple first-order phase transition. It is found that corrections to the commonly used single heavy nucleus approximation first appear at a density of the order of $10^{13}$ g/cm$^3$ and increase rapidly with increasing density. Thus, neutrino opacities are overestimated in the single heavy nucleus approximation relative to the complete molecular dynamics simulations.

Neutron-rich matter at subnuclear densities may involve complex structures displaying a variety of shapes, such as spherical, slablike, and/or rodlike shapes. These phases of the {\it nuclear pasta} are expected to exist in the crust of neutron stars and in core-collapse supernovae. The dynamics of core-collapse supernovae is very sensitive to the interactions between neutrinos and nucleons/nuclei. Indeed, neutrino excitation of the low-energy modes of the pasta may allow for a significant energy transfer to the nuclear medium, thereby reviving the stalled supernovae shock. The linear response of the nuclear pasta to neutrinos is modeled via a simple semi-classical simulation. The transport mean-free path for $\mu$ and $\tau$ neutrinos (and antineutrinos) is expressed in terms of the static structure factor of the pasta, which is evaluated using Metropolis Monte Carlo simulations.

Accurately calibrated (or ``best fit'') relativistic mean-field models are used to compute the distribution of isoscalar monopole strength in 90Zr and 208Pb, and the isovector dipole strength in 208Pb using a continuum random-phase-approximation approach. It is shown that the distribution of isoscalar monopole strength in 208Pb--but not in 90Zr--is sensitive to the density dependence of the symmetry energy. This sensitivity hinders the extraction of the compression modulus of symmetric nuclear matter from the isoscalar giant monopole resonance (ISGMR) in 208Pb. Thus, one relies on 90Zr, a nucleus with both a small neutron-proton asymmetry and a well developed ISGMR peak, to constrain the compression modulus of symmetric nuclear matter to the range K=(248 +/- 6) MeV. In turn, the sensitivity of the ISGMR in 208Pb to the density dependence of the symmetry energy is used to constrain its neutron skin to the range Rn-Rp<=0.22 fm. The impact of this result on the enhanced cooling of neutron stars is briefly addressed.

Neutron-star radii provide useful information on the equation of state of neutron rich matter. Particularly interesting is the density dependence of the equation of state (EOS). For example, the softening of the EOS at high density, where the pressure rises slower than anticipated, could signal a transition to an exotic phase. However, extracting the density dependence of the EOS requires measuring the radii of neutron stars for a broad range of masses. A ``normal'' 1.4 solar mass neutron star has a central density of a few times nuclear-matter saturation density. In contrast, low mass (of the order of 0.5 solar masses) neutron stars have central densities near nuclear-matter saturation density so its radius provides information on the EOS at low density. Unfortunately, low-mass stars are rare because they may be hard to form. Instead, a precision measurement of nuclear radii on atomic nuclei may contain similar information. Indeed, we find a strong correlation between the neutron radius of 208Pb and the radius of a 0.5 solar-mass neutron star. Thus, the radius of such a neutron star can be inferred from a measurement of the the neutron radius of 208Pb. Comparing this value to the measured radius of a 1.4 solar-mass neutron star should provide the strongest constraint to date on the density dependence of the equation of state.

Recent observations by the Chandra observatory suggest that some neutron stars may cool rapidly, perhaps by the direct URCA process which requires a high proton fraction. The proton fraction is determined by the nuclear symmetry energy whose density dependence may be constrained by measuring the neutron radius of a heavy nucleus, such as 208Pb. Such a measurement is necessary for a reliable extrapolation of the proton fraction to the higher densities present in the neutron star. A large neutron radius in 208Pb implies a stiff symmetry energy that grows rapidly with density, thereby favoring a high proton fraction and allowing direct URCA cooling. Predictions for the neutron radius in 208Pb are correlated to the proton fraction in dense matter by using a variety of relativistic effective field-theory models. Models that predict a neutron (Rn) minus proton (Rp) root-mean-square radius in 208Pb to be Rn-Rp<0.20 fm have proton fractions too small to allow the direct URCA cooling of 1.4 solar-mass neutron stars. Conversely, if Rn-Rp>0.25 fm, the direct URCA process is allowed (by all models) to cool down a 1.4 solar-mass neutron star. The Parity Radius Experiment at Jefferson Laboratory aims to measure the neutron radius in 208Pb accurately and model independently via parity-violating electron scattering. Such a measurement would greatly enhance our ability to either confirm or dismiss the direct URCA cooling of neutron stars.

Differences in the density dependence of the symmetry energy predicted by nonrelativistic and relativistic models are suggested, at least in part, as the culprit for the discrepancy in the values of the compression modulus of symmetric nuclear matter extracted from the energy of the giant monopole resonance in 208Pb. ``Best-fit'' relativistic models, with stiffer symmetry energies than Skyrme interactions, consistently predict higher compression moduli than nonrelativistic approaches. Relativistic models with compression moduli in the physically acceptable range of K=200-300 MeV are used to compute the distribution of isoscalar monopole strength in 208Pb. When the symmetry energy is artificially softened in one of these models, in an attempt to simulate the symmetry energy of Skyrme interactions, a lower value for the compression modulus is indeed obtained. It is concluded that the proposed measurement of the neutron skin in 208Pb, aimed at constraining the density dependence of the symmetry energy and recently correlated to the structure of neutron stars, will also become instrumental in the determination of the compression modulus of nuclear matter.

Parity violating electron scattering can measure the neutron density of a heavy nucleus accurately and model independently. This is because the weak charge of the neutron is much larger then that of the proton. The Parity Radius Experiment (PREX) at Jefferson Laboratory aims to measure the root mean square neutron radius of $^{208}$Pb with an absolute accuracy of 1% ($\pm 0.05$ Fm). This is more accurate then past measurements with hadronic probes, which all suffer from controversial strong interaction uncertainties. PREX should clearly resolve the neutron-rich skin. Furthermore, this benchmark value for $^{208}$Pb will provide a calibration for hadronic probes, such as proton scattering, which can then be used to measure neutron densities of many exotic nuclei. The PREX result will also have many implications for neutron stars. The neutron radius of Pb depends on the pressure of neutron-rich matter: the greater the pressure, the larger the radius as neutrons are pushed out against surface tension. The same pressure supports a neutron star against gravity. The Pb radius is sensitive to the equation of state at normal densities while the radius of a 1.4 solar mass neutron star also depends on the equation of state at higher densities. Measurements of the radii of a number of isolated neutron stars such as Geminga and RX J185635-3754 should soon improve significantly. By comparing the equation of state information from the radii of both Pb and neutron stars one can search for a softening of the high density equation of state from a phase transition to an exotic state. Possibilities include kaon condensates, strange quark matter or color superconductors.

A new relation between the neutron skin of a heavy nucleus and the radius of a neutron star is proposed: the larger the neutron skin of the nucleus the larger the radius of the star. Relativistic models that reproduce a variety of ground-state observables can not determine uniquely the neutron skin of a heavy nucleus. Thus, a large range of neutron skins is generated by supplementing the models with nonlinear couplings between isoscalar and isovector mesons. We illustrate how the correlation between the neutron skin and the radius of the star can be used to place important constraints on the equation of state and how it may help elucidate the existence of a phase transition in the interior of the neutron star.

We study relationships between the neutron-rich skin of a heavy nucleus and the properties of neutron-star crusts. Relativistic effective field theories with a thicker neutron skin in $^{208}$Pb have a larger electron fraction and a lower liquid-to-solid transition density for neutron-rich matter. These properties are determined by the density dependence of the symmetry energy which we vary by adding nonlinear couplings between isoscalar and isovector mesons. An accurate measurement of the neutron radius in $^{208}$Pb---via parity violating electron scattering---may have important implications for the structure of neutron stars.

We calculate the longitudinal response of 208Pb using a relativistic random-phase approximation to three different parameterizations of the Walecka model with scalar self-interactions. From a nonspectral calculation of the response-that automatically includes the mixing between positive- and negative-energy states-we extract the distribution of strength for the isoscalar monopole, dipole, and high-energy octupole resonances. We employ a consistent formalism that uses the same interaction in the calculation of the ground state as in the calculation of the response. As a result, the conservation of the vector current is strictly maintained throughout the calculation. Further, at small momentum transfers the spurious dipole strength-associated with the uniform translation of the center-of-mass-gets shifted to zero excitation energy and is cleanly separated from the sole remaining physical fragment located at an excitation energy of about 24 MeV; no additional dipole strength is observed. The best description of the collective modes is obtained using a ``soft'' parameterization having a compression modulus of K=224 MeV.

Core collapse supernovae are dominated by weakly interacting neutrinos. This provides a unique opportunity for macroscopic parity violation. We speculate that parity violation in a strong magnetic field can lead to an asymmetry in the explosion and a recoil of the newly formed neutron star. We estimate the asymmetry from neutrino-polarized-neutron elastic scattering, polarized electron capture and neutrino-nucleus elastic scattering in a (partially) polarized electron gas.