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.

The masses and radii of neutron stars are discussed in general relativity and scalar-tensor theory of gravity and the differences are compared with the current uncertainties stemming from the nuclear equation of state in the relativistic mean-field framework. It is shown that astrophysical and gravitational waves observations of radii of neutron stars with masses $M \lesssim 1.4 M_{\odot}$ constrain only the nuclear equation of state, and in particular the density dependence of the nuclear symmetry energy. Future observations of massive neutron stars may constrain the coupling parameters of the scalar-tensor theory provided that a general consensus on the dense nuclear matter equation of state is reached.

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 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.

We present a new mechanism for deep crustal heating in accreting neutron stars. Charged pions ($\pi^+$) are produced in nuclear collisions on the neutron star surface during active accretion and upon decay they provide a flux of neutrinos into the neutron star crust. For massive and/or compact neutron stars, neutrinos deposit $\approx 1\textrm{--} 2 \, \mathrm{MeV}$ of heat per accreted nucleon into the inner crust. The strength of neutrino heating is comparable to the previously known sources of deep crustal heating, such as from pycnonuclear fusion reactions, and is relevant for studies of cooling neutron stars. We model the thermal evolution of a transient neutron star in a low-mass X-ray binary, and in the particular case of the neutron star MXB~1659-29 we show that additional deep crustal heating requires a higher thermal conductivity for the neutron star inner crust. A better knowledge of pion production cross sections near threshold would improve the accuracy of our predictions.

The limiting rotational frequency of neutron stars may be determined by the strength of their crusts. As a star spins up from accretion, centrifugal forces will cause the crust to fail. If the crust breaks unevenly, a rotating mass quadrupole moment will radiate gravitational waves (GW). This radiation can prevent further spin up and may be a promising source for continuous GW searches. We calculate that for a breaking strain (strength) of neutron star crust that is consistent with molecular dynamics simulations, the crust may fail at rotational frequencies in agreement with observations.

We show that the neutron star in the transient system MXB~1659-29 has a core neutrino luminosity that substantially exceeds that of the modified Urca reactions (i.e., $n+n\to n+p+e^{-}+\bar{\nu}_{e}$ and inverse) and is consistent with the direct Urca reactions ($n\to p+e^{-}+\bar{\nu}_{e}$ and inverse) occurring in a small fraction of the core. Observations of the thermal relaxation of the neutron star crust following 2.5 years of accretion allow us to measure the energy deposited into the core during accretion, which is then reradiated as neutrinos, and infer the core temperature. For a nucleonic core, this requires that the nucleons are unpaired and that the proton fraction exceed a critical value to allow the direct Urca reaction to proceed. The neutron star in MXB~1659-29 is the first with a firmly detected thermal component in its X-ray spectrum that needs a fast neutrino cooling process. Measurements of the temperature variation of the neutron star core during quiescence would place an upper limit on the core specific heat and serve as a check on the fraction of the neutron star core in which nucleons are unpaired.

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.

Complex and exotic nuclear geometries are expected to appear naturally in dense nuclear matter found in the crust of neutron stars and supernovae environment collectively referred to as nuclear pasta. The pasta geometries depend on the average baryon density, proton fraction and temperature and are critically important in the determination of many transport properties of matter in supernovae and the crust of neutron stars. Using a set of self-consistent microscopic nuclear energy density functionals we present the first results of large scale quantum simulations of pasta phases at baryon densities $0.03 \leq \rho \leq 0.10$ fm$^{-3}$, proton fractions $0.05 \leq Y_p \leq 0.40$, and zero temperature. The full quantum simulations, in particular, allow us to thoroughly investigate the role and impact of the nuclear symmetry energy on pasta configurations. We use the Sky3D code that solves the Skyrme Hartree-Fock equations on a three-dimensional Cartesian grid. For the nuclear interaction we use the state of the art UNEDF1 parametrization, which was introduced to study largely deformed nuclei, hence is suitable for studies of the nuclear pasta. Density dependence of the nuclear symmetry energy is simulated by tuning two purely isovector observables that are insensitive to the current available experimental data. We find that a minimum total number of nucleons $A=2000$ is necessary to prevent the results from containing spurious shell effects and to minimize finite size effects. We find that a variety of nuclear pasta geometries are present in the neutron star crust and the result strongly depends on the nuclear symmetry energy. The impact of the nuclear symmetry energy is less pronounced as the proton fractions increase. Quantum nuclear pasta calculations at $T=0$ MeV are shown to get easily trapped in meta-stable states, and possible remedies to avoid meta-stable solutions are discussed.

We show that observations of the core temperature of transiently-accreting neutron stars combined with observations of an accretion outburst give a lower limit to the neutron star core heat capacity. For the neutron stars in the low mass X-ray binaries KS 1731-260, MXB 1659-29, and XTE J1701-462, we show that the lower limit is a factor of a few below the core heat capacity expected if neutrons and protons in the core are paired, so that electrons provide the dominant contribution to the heat capacity. This limit rules out a core dominated by a quark color-flavor-locked (CFL) phase, which would have a much lower heat capacity. Future observations of or limits on cooling during quiescence will further constrain the core heat capacity.

Neutron star and supernova matter at densities just below the nuclear matter saturation density is expected to form a lattice of exotic shapes. These so-called nuclear pasta phases are caused by Coulomb frustration. Their elastic and transport properties are believed to play an important role for thermal and magnetic field evolution, rotation and oscillation of neutron stars. Furthermore, they can impact neutrino opacities in core-collapse supernovae. In this work, we present proof-of-principle 3D Skyrme Hartree-Fock (SHF) simulations of nuclear pasta with the Multi-resolution ADaptive Numerical Environment for Scientific Simulations (MADNESS). We perform benchmark studies of $^{16} \mathrm{O}$, $^{208} \mathrm{Pb}$ and $^{238} \mathrm{U}$ nuclear ground states and calculate binding energies via 3D SHF simulations. Results are compared with experimentally measured binding energies as well as with theoretically predicted values from an established SHF code. The nuclear pasta simulation is initialized in the so-called waffle geometry as obtained by the Indiana University Molecular Dynamics (IUMD) code. The size of the unit cell is 24\:fm with an average density of about $\rho = 0.05 \:\mathrm{fm}^{-3}$, proton fraction of $Y_p = 0.3$ and temperature of $T=0\:$MeV. Our calculations reproduce the binding energies and shapes of light and heavy nuclei with different geometries. For the pasta simulation, we find that the final geometry is very similar to the initial waffle state. In the present pasta calculations spin-orbit forces are not included but will be added in the future. Within the MADNESS framework, we can successfully perform calculations of inhomogeneous nuclear matter. By using pasta configurations from IUMD it is possible to explore different geometries and test the impact of self-consistent calculations on the latter.

A thorough understanding of properties of neutron stars requires both a reliable knowledge of the equation of state (EOS) of super-dense nuclear matter and the strong-field gravity theories simultaneously. To provide information that may help break this EOS-gravity degeneracy, we investigate effects of nuclear symmetry energy on the gravitational binding energy of neutron stars within GR and the scalar-tensor subset of alternative gravity models. We focus on effects of the slope $L$ of nuclear symmetry energy at saturation density and the high-density behavior of nuclear symmetry energy. We find that the variation of either the density slope $L$ or the high-density behavior of nuclear symmetry energy leads to large changes in the binding energy of neutron stars. The difference in predictions using the GR and the scalar-tensor theory appears only for massive neutron stars, and even then is significantly smaller than the difference resulting from variations in the symmetry energy.

At present, neutron star radii from both observations and model predictions remain very uncertain. Whereas different models can predict a wide range of neutron star radii, it is not possible for most models to predict radii that are smaller than about 10 km, thus if such small radii are established in the future they will be very difficult to reconcile with model estimates. By invoking a new term in the equation of state that enhances the energy density, but leaves the pressure unchanged we simulate the current uncertainty in the neutron star radii. This new term can be possibly due to the exchange of the weakly interacting light U-boson with appropriate in-medium parameters, which does not compromise the success of the conventional nuclear models. The validity of this new scheme will be tested eventually by more precise measurements of neutron star radii.

The critical densities and impact of forming \D resonances in neutron stars are investigated within an extended nonlinear relativistic mean-field (RMF) model. The critical densities for the formation of four different charge states of \D are found to depend differently on the separate kinetic and potential parts of nuclear symmetry energy, the first example of a microphysical property of neutron stars to do so. Moreover, they are sensitive to the in-medium Delta mass $m_{\Delta}$ and the completely unknown $\Delta$-$\rho$ coupling strength $g_{\rho\Delta}$. In the universal baryon-meson coupling scheme where the respective $\Delta$-meson and nucleon-meson coupling constants are assumed to be the same, the critical density for the first $\Delta^-(1232)$ to appear is found to be \rc=$(2.08\pm0.02)\rho_0$ using RMF model parameters consistent with current constraints on all seven macroscopic parameters usually used to characterize the equation of state (EoS) of isospin-asymmetric nuclear matter (ANM) at saturation density $\rho_0$. Moreover, the composition and the mass-radius relation of neutron stars are found to depend significantly on the values of the $g_{\rho\Delta}$ and $m_{\Delta}$.

A number of observed phenomena associated with individual neutron star systems or neutron star populations find explanations in models in which the neutron star crust plays an important role. We review recent work examining the sensitivity to the slope of the symmetry energy $L$ of such models, and constraints extracted on $L$ from confronting them with observations. We focus on six sets of observations and proposed explanations: (i) The cooling rate of the neutron star in Cassiopeia A, confronting cooling models which include enhanced cooling in the nuclear pasta regions of the inner crust, (ii) the upper limit of the observed periods of young X-ray pulsars, confronting models of magnetic field decay in the crust caused by the high resistivity of the nuclear pasta layer, (iii) glitches from the Vela pulsar, confronting the paradigm that they arise due to a sudden re-coupling of the crustal neutron superfluid to the crustal lattice after a period during which they were decoupled due to vortex pinning, (iv) The frequencies of quasi-periodic oscillations in the X-ray tail of light curves from giant flares from soft gamma-ray repeaters, confronting models of torsional crust oscillations, (v) the upper limit on the frequency to which millisecond pulsars can be spun-up due to accretion from a binary companion, confronting models of the r-mode instability arising above a threshold frequency determined in part by the viscous dissipation timescale at the crust-core boundary, and (vi) the observations of precursor electromagnetic flares a few seconds before short gamma-ray bursts, confronting a model of crust shattering caused by resonant excitation of a crustal oscillation mode by the tidal gravitational field of a companion neutron star just before merger.

We perform a systematic assessment of models for the equation of state (EOS) of dense matter in the context of recent neutron star mass and radius measurements to obtain a broad picture of the structure of neutron stars. We demonstrate that currently available neutron star mass and radius measurements provide strong constraints on moments of inertia, tidal deformabilities, and crust thicknesses. A measurement of the moment of inertia of PSR J0737-3039A with 10% error, without any other information from observations, will constrain the EOS over a range of densities to within 50%$-$60%. We find tidal deformabilities between 0.6 and $6\times 10^{36}$ g cm$^{2}$ s$^{2}$ (to 95% confidence) for $M=1.4~\mathrm{M}_{\odot}$, and any measurement which constrains this range will provide an important constraint on dense matter. The crustal fraction of the moment of inertia can be as large as 10% for $M=1.4~\mathrm{M}_{\odot}$ permitting crusts to have a large enough moment of inertia reservoir to explain glitches in the Vela pulsar even with a large amount of superfluid entrainment. Finally, due to the uncertainty in the equation of state, there is at least a 40% variation in the thickness of the crust for a fixed mass and radius, which implies that future simulations of the cooling of a neutron star crust which has been heated by accretion will need to take this variation into account.

There is a degeneracy between the equation of state (EOS) of superdense neutron-rich nuclear matter and the strong-field gravity in understanding properties of neutron stars. While the EOS is still poorly known, there are also longstanding ambiguities in choosing General Relativity or alternative gravity theories in the not-so-well tested strong-field regime. Besides possible appearance of hyperons and new phases, the most uncertain part of the nucleonic EOS is currently the density dependence of nuclear symmetry energy. To provide information that may help break the EOS-gravity degeneracy, we investigate effects of symmetry energy within its uncertain range determined by terrestrial nuclear laboratory experiments on the gravitational binding energy and spacetime curvature of neutron stars within GR and the scalar-tensor (ST) theory of gravity. In particular, we focus on effects of the following parameters characterizing the EOS of neutron-rich nucleonic matter: (1) the incompressibility $K_0$ of symmetric nuclear matter, (2) the slope $L$ of symmetry energy at saturation density and (3) the high-density behavior of symmetry energy. We find that the variation of either the density slope $L$ or the high-density behavior of symmetry energy leads to large changes in both the binding energy and curvature of neutron stars while effects of varying the more constrained $K_0$ are negligibly small. The difference in predictions using the GR and the ST theory appears only for massive neutron stars, and is significantly smaller than the differences resulting from variations in the symmetry energy. We conclude that within the ST theory of gravity, the EOS-gravity degeneracy has been broken by the recent relativistic pulsar measurements, and that measurements of neutron star properties sensitive to the compactness constrain mainly the density dependence of the symmetry energy.

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.

According to the Hugenholtz-Van Hove theorem, the nuclear symmetry energy $S(\rho)$ and its slope $L(\rho)$ at arbitrary densities can be decomposed in terms of the density and momentum dependence of the single-nucleon potentials in isospin-asymmetric nuclear matter which are potentially accessible to experiment. We quantify the correlations between several well-known isovector observables and $L(\rho)$ to locate the density range in which each isovector observable is most sensitive to the density dependence of the $S(\rho)$. We then study the correlation coefficients between those isovector observables and all the components of the $L(\rho)$. The neutron skin thickness of $^{208}$Pb is found to be strongly correlated with the $L(\rho)$ at a subsaturation density of $\rho = 0.59 \rho_0$ through the density dependence of the first-order symmetry potential. Neutron star radii are found to be strongly correlated with the $L(\rho)$ over a wide range of supra-saturation densities mainly through both the density and momentum dependence of the first-order symmetry potential. Finally, we find that although the crust-core transition pressure has a complex correlation with the $L(\rho)$, it is strongly correlated with the momentum derivative of the first-order symmetry potential, and the density dependence of the second-order symmetry potential.

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.

Gravitational wave (GW) astronomy opens up an entirely new window on the Universe to probe the equations of state (EOS) of neutron-rich matter. With the advent of next generation GW detectors, measuring the gravitational radiation from coalescing binary neutron star systems, mountains on rotating neutron stars, and stellar oscillation modes may become possible in the near future. Using a set of model EOSs satisfying the latest constraints from terrestrial nuclear experiments, state of the art nuclear many-body calculations of the pure neutron matter EOS, and astrophysical observations consistently, we study various GW signatures of the high-density behavior of the nuclear symmetry energy, which is considered among the most uncertain properties of dense neutron-rich nucleonic matter. In particular, we find the tidal polarizability of neutron stars, potentially measurable in binary systems just prior to merger, is more sensitive to the high density component of the nuclear symmetry energy than the symmetry energy at nuclear saturation density. We also find that the upper limit on the GW strain amplitude from elliptically deformed stars is very sensitive to the density dependence of the symmetry energy. This suggests that future developments in modeling of the neutron star crust, and direct gravitational wave signals from accreting binaries will provide a wealth of information on the EOS of neutron-rich matter. We also review the sensitivity of the $r$-mode instability window to the density dependence of the symmetry energy. Whereas models with larger values of the density slope of the symmetry energy at saturation seem to be disfavored by the current observational data, within a simple $r$-mode model, we point out that a subsequent softer behavior of the symmetry energy at high densities (hinted at by recent observational interpretations) could rule them in.

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.

Significant progress has been made in recent years in constraining nuclear symmetry energy at and below the saturation density of nuclear matter using data from both terrestrial nuclear experiments and astrophysical observations. However, many interesting questions remain to be studied especially at supra-saturation densities. In this lecture note, after a brief summary of the currently available constraints on nuclear symmetry energy near the saturation density we first discuss the relationship between the symmetry energy and the isopin and momentum dependence of the single-nucleon potential in isospin-asymmetric nuclear medium. We then discuss several open issues regarding effects of the tensor force induced neutron-proton short-range correlation (SRC) on nuclear symmetry energy. Finally, as an example of the impacts of nuclear symmetry energy on properties of neutron stars and gravitational waves, we illustrate effects of the high-density symmetry energy on the tidal polarizability of neutron stars in coalescing binaries.

Using a set of model equations of state satisfying the latest constraints from both terrestrial nuclear experiments and astrophysical observations as well as state-of-the-art nuclear many-body calculations of the pure neutron matter equation of state, the tidal polarizability of canonical neutron stars in coalescing binaries is found to be a very sensitive probe of the high-density behavior of nuclear symmetry energy which is among the most uncertain properties of dense neutron-rich nucleonic matter. Moreover, it changes less than $\pm 10%$ by varying various properties of symmetric nuclear matter and symmetry energy around the saturation density within their respective ranges of remaining uncertainty.

In this review, we will discuss the results of our recent work to study the general optimization of the pure isovector parameters of the popular relativistic mean-field (RMF) and Skyrme-Hartree-Fock (SHF) nuclear energy-density functionals (EDFs), using constraints on the pure neutron matter (PNM) equation of state (EoS) from recent {\sl ab initio} calculations. By using RMF and SHF parameterizations that give equivalent predictions for ground-state properties of doubly magic nuclei and properties of symmetric nuclear matter (SNM) and PNM, we found that such optimization leads to broadly consistent symmetry energy $J$ and its slope parameter $L$ at saturation density within a tight range of $\sigma(J) < 2$ MeV and $\sigma(L) < 6$ MeV. We demonstrate that a clear model dependence shows up (a) in the curvature parameter of the symmetry energy $K_{\rm sym}$, (b) the symmetry energy at supra-saturation densities, and (c) the radius of neutron stars.

We study the nuclear symmetry energy S(rho) and related quantities of nuclear physics and nuclear astrophysics predicted generically by relativistic mean-field (RMF) and Skyrme-Hartree-Fock (SHF) models. We establish a simple prescription for preparing equivalent RMF and SHF parametrizations starting from a minimal set of empirical constraints on symmetric nuclear matter, nuclear binding energy and charge radii, enforcing equivalence of their Lorenz effective masses, and then using the pure neutron matter (PNM) equation of state (EoS) obtained from ab-initio calculations to optimize the pure isovector parameters in the RMF and SHF models. We find the resulting RMF and SHF parametrizations give broadly consistent predictions of the symmetry energy J and its slope parameter L at saturation density within a tight range of <~2 MeV and <~6 MeV respectively, but that clear model dependence shows up in the predictions of higher-order symmetry energy parameters, leading to important differences in (a) the slope of the correlation between J and L from the confidence ellipse, (b) the isospin-dependent part of the incompressibility of nuclear matter K_tau, (c) the symmetry energy at supra-saturation densities, and (d) the predicted neutron star radii. The model dependence can lead to about 1-2 km difference in predictions of the neutron star radius given identical predicted values of J, L and symmetric nuclear matter (SNM) saturation properties. Allowing the full freedom in the effective masses in both models leads to constraints of 30<~J<~31.5 MeV, 35<~L<~60 MeV, -330<~K_tau<~-216 MeV for the RMF model as a whole and 30<~J<~33 MeV, 28<~L<~65 MeV, -420<~K_tau<~-325 MeV for the SHF model as a whole. Notably, given PNM constraints, these results place RMF and SHF models as a whole at odds with some constraints on K_tau inferred from giant monopole resonance and neutron skin experimental results.

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.

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.

We study the dipolar magnetic field configuration in dependence on brane tension and present solutions of Maxwell equations in the internal and external background spacetime of a magnetized spherical star in a Randall-Sundrum II type braneworld. The star is modelled as sphere consisting of perfect highly magnetized fluid with infinite conductivity and frozen-in dipolar magnetic field. With respect to solutions for magnetic fields found in the Schwarzschild spacetime brane tension introduces enhancing corrections both to the interior and the exterior magnetic field. These corrections could be relevant for the magnetic fields of magnetized compact objects as pulsars and magnetars and may provide the observational evidence for the brane tension through the modification of formula for magneto-dipolar emission which gives amplification of electromagnetic energy loss up to few orders depending on the value of the brane tension.

The general principles needed to compute the effect of a stationary gravitational field on the quasistationary electromagnetic phenomena in normal conductors and superconductors are formulated from general relativistic point of view. Generalization of the skin effect, that is the general relativistic modification of the penetration depth (of the time-dependent magnetic field in the conductor) due to its relativistic coupling to the gravitational field is obtained. The effect of the gravitational field on the penetration and coherence depths in superconductors is also studied. As an illustration of the foregoing general results, we discuss their application to superconducting systems in the outer core of neutron stars. The relevance of these effects to electrodynamics of magnetized neutron stars has been shown.