The nature of the interaction between superfluid vortices and the neutron star crust, conjectured by Anderson and Itoh in 1975 to be at the heart vortex creep and the cause of glitches, has been a long-standing question in astrophysics. Using a qualitatively new approach, we follow the dynamics as superfluid vortices move in response to the presence of "nuclei" (nuclear defects in the crust). The resulting motion is perpendicular to the force, similar to the motion of a spinning top when pushed. We show that nuclei repel vortices in the neutron star crust, and characterize the force as a function of the vortex-nucleus separation.

We investigate the dynamics of a quantized vortex and a nuclear impurity immersed in a neutron superfluid within a fully microscopic time-dependent three-dimensional approach. The magnitude and even the sign of the force between the quantized vortex and the nuclear impurity have been a matter of debate for over four decades. We determine that the vortex and the impurity repel at neutron densities, 0.014 fm$^{-3}$ and 0.031 fm$^{-3}$, which are relevant to the neutron star crust and the origin of glitches, while previous calculations have concluded that the force changes its sign between these two densities and predicted contradictory signs. The magnitude of the force increases with the density of neutron superfluid, while the magnitude of the pairing gap decreases in this density range.

It is shown that in an inhomogeneous neutron matter the pairing field bounds neutrons around the Fermi level leading to the formation of localized Andreev states. In the case of the inner crust of neutron stars the localization length has been determined as a function of the nuclear density.

The inner crust of a neutron star consists of nuclei immersed in a superfluid neutron liquid. As these nuclei move through the fermionic medium they bring it into motion as well. As a result their mass is strongly renormalized and the spectrum of the ion lattice vibrations is significantly affected. Consequently we predict that the specific heat and the lattice thermal energy of the Coulomb crystal at these densities are noticeably modified.

In the inner crust of a neutron star, due to the high density and pressure, nuclei which are still present, are immersed in a neutron superfluid. One then expects that the dynamical properties of nuclei are significantly affected. In order to estimate the magnitude of the effect associated with the presence of a superfluid medium, we formulate the hydrodynamical approach to the nuclear dynamics in the inner crust of neutron stars. We calculate the renormalized nuclear mass and the strength of the medium-induced interaction between nuclei. We argue that these effects noticeably modify the properties of the Coulomb crystal in the inner crust.

The bottom part of the neutron star crust is investigated using the Skyrme-Hartree-Fock approach with the Coulomb interaction treated beyond the Wigner-Seitz approximation. A variety of nuclear phases is found to coexist in this region. Their stability and relative energies are governed by the Coulomb, surface and shell energies. We have also found that a substantial contribution is coming from the spin-orbit interaction.

We estimate the quantum corrections to the ground state energy in neutron matter (which could be termed as well either shell correction energy or Casimir energy) at subnuclear densities, where various types of inhomogeneities (bubbles, rods, plates) are energetically favorable. We show that the magnitude of these energy corrections are comparable to the energy differences between various types of inhomogeneous phases. We discuss the dependence of these corrections on a number of physical parameters (density, filling factor, temperature, lattice distortions).

We investigate the inner crust structure of neutron stars using the Skyrme-Hartree-Fock approach with the Coulomb interaction treated beyond the Wigner-Seitz approximation. Our results suggest that the shell effects associated with unbound neutrons play an important role and, in particular, lead to complicated phase transition pattern between various nuclear phases (as a function of the density). Namely, we show that the relative energies of different phases are rapidly oscillating functions of the neutron density. In the semiclassical approach this behavior is explained as an interference effect due to periodic orbits of similar lengths. We discuss also the dependence of the shell effects on pairing correlations.

In the present study we investigate the static properties of nuclei in the inner crust of neutron stars. Using the Hartree-Fock method in coordinate space, together with the semiclassical approximation, we examine the patterns of phase transitions.

The inner crust of neutron stars consists of nuclei of various shapes immersed in a neutron gas and stabilized by the Coulomb interaction in the form of a crystal lattice. The scattering of neutrons on nuclear inhomegeneities leads to the quantum correction to the total energy of the system. This correction resembles the Casimir energy and turns out to have a large influence on the structure of the crust.

Homogeneous neutron matter at subnuclear densities becomes unstable towards the formation of inhomogeneities. Depending on the average value of the neutron density one can observe the appearance of either bubbles, rods, tubes or plates embeded in a neutron gas. We estimate the quantum corrections to the ground state energy (which could be termed either shell correction or Casimir energy) of such phases of neutron matter. The calculations are performed by evaluating the contribution of the shortest periodic orbits in the Gutzwiller trace formula for the density of states. The magnitude of the quantum corrections to the ground state energy of neutron matter are of the same order as the energy differences between various phases.