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.

We present the simplest nuclear energy density functional (NEDF) to date, determined by only 4 significant phenomenological parameters, yet capable of fitting measured nuclear masses with better accuracy than the Bethe-Weizs\"acker mass formula, while also describing density structures (charge radii, neutron skins etc.) and time-dependent phenomena (induced fission, giant resonances, low energy nuclear collisions, etc.). The 4 significant parameters are necessary to describe bulk nuclear properties (binding energies and charge radii); an additional 2 to 3 parameters have little influence on the bulk nuclear properties, but allow independent control of the density dependence of the symmetry energy and isovector excitations, in particular the Thomas-Reiche-Kuhn sum rule. This Hohenberg-Kohn-style of density functional theory successfully realizes Weizs\"acker's ideas and provides a computationally tractable model for a variety of static nuclear properties and dynamics, from finite nuclei to neutron stars, where it will also provide a new insight into the physics of the r-process, nucleosynthesis, and neutron star crust structure. This new NEDF clearly separates the bulk geometric properties - volume, surface, symmetry, and Coulomb energies which amount to 8MeV per nucleon or up to 2000MeV per nucleus for heavy nuclei - from finer details related to shell effects, pairing, isospin breaking, etc. which contribute at most a few MeV for the entire nucleus. Thus it provides a systematic framework for organizing various contributions to the NEDF. Measured and calculated physical observables - symmetry and saturation properties, the neutron matter equation of state, and the frequency of giant dipole resonances - lead directly to new terms not considered in current NEDF parameterizations.

We present variational Monte Carlo calculations of the neutron matter equation of state using chiral nuclear forces. The ground-state wavefunction of neutron matter, containing non-perturbative many-body correlations, is obtained from auxiliary-field quantum Monte Carlo simulations of up to about 340 neutrons interacting on a 10^3 discretized lattice. The evolution Hamiltonian is chosen to be attractive and spin-independent in order to avoid the fermion sign problem and is constructed to best reproduce broad features of the chiral nuclear force. This is facilitated by choosing a lattice spacing of 1.5 fm, corresponding to a momentum-space cutoff of Lambda = 414 MeV/c, a resolution scale at which strongly repulsive features of nuclear two-body forces are suppressed. Differences between the evolution potential and the full chiral nuclear interaction (Entem and Machleidt Lambda = 414 MeV) are then treated perturbatively. Our results for the equation of state are compared to previous quantum Monte Carlo simulations which employed chiral two-body forces at next-to-next-to-leading order (N2LO). In addition we include the effects of three-body forces at N2LO, which provide important repulsion at densities higher than 0.02 fm^-3, as well as two-body forces at N3LO.

We present an efficient and general method to compute vortex-pinning interactions - which arise in neutron stars, superconductors, and trapped cold atoms - at arbitrary separations using real-time dynamics. This method overcomes uncertainties associated with matter redistribution by the vortex position and the related choice of ensemble that plague the typical approach of comparing energy differences between stationary pinned and unpinned configurations: uncertainties that prevent agreement in the literature on the sign and magnitude of the vortex-nucleus interaction in the crust of neutron stars. We demonstrate and validate the method with Gross-Pitaevskii-like equations for the unitary Fermi gas, and demonstrate how the technique of adiabatic state preparation with time-dependent simulation can be used to calculate vortex-pinning interactions in fermionic systems such as the vortex-nucleus interaction in the crust of neutron stars.

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.

We study in a fully selfconsistent approach the structure of a vortex in low density superfluid neutron matter. We determine that the matter density profile of a vortex shows a significant depletion in the region of the core, a feature never reported for a vortex state in a Fermi superfluid.

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).

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.

Using ensembles of two, three and four spheres immersed in a fermionic background we evaluate the (integrated) density of states and the Casimir energy. We thus infer that for sufficiently smooth objects, whose various geometric characteristic lengths are larger then the Fermi wave length one can use the simplest semiclassical approximation (the contribution due shortest periodic orbits only) to evaluate the Casimir energy. We also show that the Casimir energy for several objects can be represented fairly accurately as a sum of pairwise Casimir interactions between pairs of objects.

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.