No Arabic abstract
The functional form of the nuclear symmetry energy in the whole range of densities relevant for the neutron stars is still unknown. Discrepancies concern both the low as well as the high density behaviour of this function. By use of Bezier curves three different families of the symmetry energy shapes, relevant for different density range were introduced. Their consequences for the crustal properties of neutron stars are presented.
The form of the nuclear symmetry energy $E_s$ around saturation point density leads to a different crust-core transition point in the neutron star and affect the crust properties. We show that the knowledge about $E_s$ close to the saturation point is not sufficient, because the very low density behaviour is relevant. We also claim that crust properties are strongly influenced by the very high density behavior of $E_s$, so in order to conclude about the form of low density part of the symmetry energy one must isolate properly the high density part.
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 breaking stress (the maximum of the stress-strain curve) of neutron star crust is important for neutron star physics including pulsar glitches, emission of gravitational waves from static mountains, and flares from star quakes. We perform many molecular dynamic simulations of the breaking stress at different coupling parameters (inverse temperatures) and strain rates. We describe our results with the Zhurkov model of strength. We apply this model to estimate the breaking stress for timescales ~1 s - 1 year, which are most important for applications, but much longer than can be directly simulated. At these timescales the breaking stress depends strongly on the temperature. For coupling parameter <200, matter breaks at very small stress, if it is applied for a few years. This viscoelastic creep can limit the lifetime of mountains on neutron stars. We also suggest an alternative model of timescale-independent breaking stress, which can be used to estimate an upper limit on the breaking stress.
The strength of neutron star crust is crucial for modelling magnetar flares, pulsar glitches and gravitational wave emission. We aim to shed some light on this problem by analysing uniaxial stretch deformation (elongation and contraction) of perfect body-centered cubic Coulomb crystals, paying special attention to the inherent anisotropy of this process. Our analysis is based on the semi-analytical approach of Baiko and Kozhberov (2017), which, for any uniform deformation, allows one to calculate, in fully non-linear regime, critical deformation parameters beyond which the lattice loses its dynamic stability. We determine critical strain, pressure anisotropy and deformation energy for any stretch direction with respect to the crystallographic axes. These quantities are shown to be strongly anisotropic: they vary by a factor of almost 10 depending on the orientation of the deformation axis. For polycrystalline crust, we argue that the maximum strain for the stretch deformation sustainable elastically is 0.04. It is lower than the breaking strain of 0.1 obtained in molecular dynamic simulations of a shear deformation by Horowitz and Kadau (2009). The maximum pressure anisotropy of polycrystalline matter is estimated to be in the range from 0.005 to 0.014 $nZ^2e^2/a$, where $n$ is the ion number density, $Ze$ is the ion charge, and $a$ is the ion-sphere radius. We discuss possible mechanisms of plastic motion and formation of large crystallites in neutron star crust as well as analyse energy release associated with breaking of such crystallites in the context of magnetic field evolution and magnetar flaring activity.
We demonstrate that observations of glitches in the Vela pulsar can be used to investigate the strength of the crust-core coupling in a neutron star, and suggest that recovery from the glitch is dominated by torque exerted by the re-coupling of superfluid components of the core that were decoupled from the crust during the glitch. Assuming that the recoupling is mediated by mutual friction between the superfluid neutrons and the charged components of the core, we use the observed magnitudes and timescales of the shortest timescale components of the recoveries from two recent glitches in the Vela pulsar to infer the fraction of the core that is coupled to the crust during the glitch, and hence spun up by the glitch event. Within the framework of a two-fluid hydrodynamic model of glitches, we analyze whether crustal neutrons alone are sufficient to drive the glitch activity observed in the Vela pulsar. We use two sets of neutron star equations of state (EOSs), both of which span crust and core consistently and cover a range of the slope of the symmetry energy at saturation density $30 < L <120$ MeV. One set produces maximum masses $approx$2.0$M_{odot}$, the second $approx$2.6$M_{odot}$. We also include the effects of entrainment of crustal neutrons by the superfluid lattice. We find that for medium to stiff EOSs, observations imply $>70%$ of the moment of inertia of the core is coupled to the crust during the glitch, though for softer EOSs $Lapprox 30$MeV as little as $5%$ could be coupled. No EOS is able to reproduce the observed glitch activity with crust neutrons alone, but extending the region where superfluid vortices are strongly pinned into the core by densities as little as 0.016fm$^{-3}$ above the crust-core transition density restores agreement with the observed glitch activity.