The crust of a neutron star is thought to be comprised of a lattice of nuclei immersed in a sea of free electrons and neutrons. As the neutrons are superfluid their angular momentum is carried by an array of quantized vortices. These vortices can pin to the nuclear lattice and prevent the neutron superfluid from spinning down, allowing it to store angular momentum which can then be released catastrophically, giving rise to a pulsar glitch. A crucial ingredient for this model is the maximum pinning force that the lattice can exert on the vortices, as this allows us to estimate the angular momentum that can be exchanged during a glitch. In this paper we perform, for the first time, a detailed and quantitative calculation of the pinning force \emph{per unit length} acting on a vortex immersed in the crust and resulting from the mesoscopic vortex-lattice interaction. We consider realistic vortex tensions, allow for displacement of the nuclei and average over all possible orientation of the crystal with respect to the vortex. We find that, as expected, the mesoscopic pinning force becomes weaker for longer vortices and is generally much smaller than previous estimates, based on vortices aligned with the crystal. Nevertheless the forces we obtain still have maximum values of order $f_{\rm{pin}}\approx 10^{15}$ dyn/cm, which would still allow for enough angular momentum to be stored in the crust to explain large Vela glitches, if part of the star is decoupled during the event.

The cell structure of clusters in the inner crust of a cold \beta-equilibrium neutron star is studied within a Thomas Fermi approach and compared with other approaches which include shell effects. Relativistic nuclear models are considered. We conclude that the symmetry energy slope L may have quite dramatic effects on the cell structure if it is very large or small. Rod-like and slab-like pasta clusters have been obtained in all models except one with a large slope L.

We analyse how the structure of the inner curst is influenced by the pairing correlations. The inner-crust matter, formed by nuclear clusters immersed in a superfluid neutron gas and ultra-relativistic electrons, is treated in the Wigner-Seitz approximation. The properties of the Wigner-Seitz cells, i.e., their neutron to proton ratio and their radius at a given baryonic density, are obtained from the energy minimization at beta equilibrium. To obtain the binding energy of baryonic matter we perform Skyrme-HFB calculations with zero-range density-dependent pairing forces of various intensities. We find that the Wigner-Seitz cells have much smaller numbers of protons compared to previous calculations. For the dense cells the binding energy of the configurations with small proton numbers do not converge to a well-defined minimum value which precludes the determination of their structure. We show that for these cells there is a significant underestimation of the binding energy due to the boundary conditions at the border of the cells imposed through the Wigner-Seitz approximation.

We discuss the thermalisation process of the neutron stars crust described by solving the heat transport equation with a microscopic input for the specific heat of baryonic matter. The heat equation is solved with initial conditions specific to a rapid cooling of the core. To calculate the specific heat of inner crust baryonic matter, i.e., nuclear clusters and unbound neutrons, we use the quasiparticle spectrum provided by the Hartree-Fock-Bogoliubov approach at finite temperature. In this framework we analyse the dependence of the crust thermalisation on pairing properties and on cluster structure of inner crust matter. It is shown that the pairing correlations reduce the crust thermalisation time by a very large fraction. The calculations show also that the nuclear clusters have a non-negligible influence on the time evolution of the surface temperature of the neutron star.