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.

We investigate the nuclear collective excitations of Wigner-Seitz cells containing nuclear clusters immersed in a gas of neutrons. This baryonic non-uniform system is specific to the structure of inner crust matter of neutron stars. The collective excitations are studied in the framework of a spherical Hartree-Fock-Bogoliubov + Quasiparticle Random Phase Approximation, formulated in coordinate representation. The calculations are done for two representative Wigner-Seitz cells with baryonic density equal to 0.02 fm$^{-3}$ and 0.08 fm$^{-3}$. It is shown that the excitations with low multipolarities are concentrated almost entirely in one strongly collective mode which exhausts a very large fraction of the energy-weighted sum rule. Since these collective modes are located at very low energies compared to the giant resonances in standard nuclei, they may affect significantly the specific heat of baryonic inner crust matter of neutron stars.

We analyse the temperature dependence of pairing correlations in the inner crust matter of neutron stars. The study is done in a finite-temperature HFB approach and by using a zero range pairing force adjusted to the pairing properties of infinite neutron matter. Within the same approach we investigate how the specific heat of the inner crust depends on temperature, matter inhomogeneity, and the assumption used for the pairing force. It is shown that in a physical relevant range of densities the pairing properties of inner crust matter depend significantly on temperature. The finite-temperature HFB calculations show also that the specific heat is rather sensitive to the presence of nuclear clusters inside the inner crust. However, the most dramatic change of the specific heat is determined by the scenario used for the neutron matter superfluidity.