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.

The modeling of many neutron star observables incorporates the microphysics of both the stellar crust and core, which is tied intimately to the properties of the nuclear matter equation of state (EoS). We explore the predictions of such models over the range of experimentally constrained nuclear matter parameters, focusing on the slope of the symmetry energy at nuclear saturation density $L$. We use a consistent model of the composition and EoS of neutron star crust and core matter to model the binding energy of pulsar B of the double pulsar system J0737-3039, the frequencies of torsional oscillations of the neutron star crust and the instability region for r-modes in the neutron star core damped by electron-electron viscosity at the crust-core interface. By confronting these models with observations, we illustrate the potential of astrophysical observables to offer constraints on poorly known nuclear matter parameters complementary to terrestrial experiments, and demonstrate that our models consistently predict $L<70$ MeV.

We present a systematic survey the range of predictions of the neutron star inner crust composition, crust-core transition densities and pressures, and density range of the nuclear `pasta' phases at the bottom of the crust provided by the compressible liquid drop model in the light of current experimental and theoretical constraints on model parameters. Using a Skyrme-like model for nuclear matter, we construct baseline sequences of crust models by consistently varying the density dependence of the bulk symmetry energy at nuclear saturation density, $L$, under two conditions: (i) that the magnitude of the symmetry energy at saturation density $J$ is held constant, and (ii) $J$ correlates with $L$ under the constraint that the pure neutron matter (PNM) EoS satisfies the results of ab-initio calculations at low densities. Such baseline crust models facilitate consistent exploration of the $L$ dependence of crustal properties. The remaining surface energy and symmetric nuclear matter parameters are systematically varied around the baseline, and different functional forms of the PNM EoS at sub-saturation densities implemented, to estimate theoretical `error bars' for the baseline predictions. Inner crust composition and transition densities are shown to be most sensitive to the surface energy at very low proton fractions and to the behavior of the sub-saturation PNM EoS. Recent calculations of the energies of neutron drops suggest that the low-proton-fraction surface energy might be higher than predicted in Skyrme-like models, which our study suggests may result in a greatly reduced volume of pasta in the crust than conventionally predicted.

The structure and composition of the inner crust of neutron stars, as well as global stellar properties such as radius and moment of inertia, have been shown to correlate with parameters characterizing the symmetry energy of nuclear matter such as its magnitude J and density dependence L at saturation density. It is thus mutually beneficial to nuclear physicists and astrophysicists to examine the combined effects of such correlations on potential neutron star observables in the light of recent experimental and theoretical constraints on J, L, and relationships between them. We review some basic correlations between these nuclear and astrophysical observables, and illustrate the impact of recent progress in constraining the J-L parameter space on the composition of the inner crust, crust-core transition density and pressure, and extent of the hypothesized pasta region. We use a simple compressible liquid drop model in conjunction with a simple model of nuclear matter which allows for independent, smooth, variation of the J and L. We extend the model into the core using the same nuclear matter model to explore the effects on global crust and core properties, and on potential observables such as crust oscillation frequencies and mechanically supported crust deformation. Throughout we illustrate the importance of the relationship between J and L implicit in a particular model of nuclear matter to the predictions of neutron star properties.

The effects of the existence of exotic nuclear shapes at the bottom of the neutron star inner crust - nuclear `pasta' - on observational phenomena are estimated by comparing the limiting cases that those phases have a vanishing shear modulus and that they have the shear modulus of a crystalline solid . We estimate the effect on torsional crustal vibrations and on the maximum quadrupole ellipticity sustainable by the crust. The crust composition and transition densities are calculated consistently with the global properties, using a liquid drop model with a bulk nuclear equation of state (EoS) which allows a systematic variation of the nuclear symmetry energy. The symmetry energy J and its density dependence L at nuclear saturation density are the dominant nuclear inputs which determine the thickness of the crust, the range of densities at which pasta might appear, as well as global properties such as the radius and moment of inertia. We show the importance of calculating the global neutron star properties on the same footing as the crust EoS, and demonstrate that in the range of experimentally acceptable values of L, the pasta phase can alter the crust frequencies by up to a factor of three, exceeding the effects of superfluidity on the crust modes, and decrease the maximum quadrupole ellipticity sustainable by the crust by up to an order of magnitude. The signature of the pasta phases and the density dependence of the symmetry energy on the potential observables highlights the possibility of constraining the EoS of dense, neutron-rich matter and the properties of the pasta phases using astrophysical observations.

Significant progress has been made in recent years in constraining the density dependence of nuclear symmetry energy using terrestrial nuclear laboratory data. Around and below the nuclear matter saturation density, the experimental constraints start to merge in a relatively narrow region. At supra-saturation densities, there are, however, still large uncertainties. After summarizing the latest experimental constraints on the density dependence of nuclear symmetry energy, we highlight a few recent studies examining imprints of nuclear symmetry energy on the binding energy, energy release during hadron-quark phase transitions as well as the $w$-mode frequency and damping time of gravitational wave emission of neutron stars.