We briefly examine the properties of dense plasmas characteristic of the atmospheres of neutron stars and of the interior of massive white dwarfs. These astrophysical bodies are natural laboratories to study respectively the problem of pressure ionization of hydrogen in a strong magnetic field and the crystallization of the quantum one-component-plasma at finite temperature.

Equilibrium configurations of cold neutron stars near the minimum mass are studied, using the recent equation of state SLy, which describes in a unified, physically consistent manner, both the solid crust and the liquid core of neutron stars. Results are compared with those obtained using an older FPS equation of state of cold catalyzed matter. The value of M_min\simeq 0.09M_sun depends very weakly on the equation of state of cold catalyzed matter: it is 0.094 M_sun for the SLy model, and 0.088 M_sun for the FPS one. Central density at M_min is significantly lower than the normal nuclear density: for the SLy equation of state we get central density 1.7 10^{14} g/cm^3, to be compared with 2.3 10^{14} g/cm^3 obtained for the FPS one. Even at M_min, neutron stars have a small liquid core of radius of about 4 km, containing some 2-3% of the stellar mass. Neutron stars with 0.09 M_sun <M<0.17 M_sun are bound with respect to dispersed configuration of the hydrogen gas, but are unbound with respect to dispersed Fe^56. The effect of uniform rotation on the minimum-mass configuration of cold neutron stars is studied. Rotation increases the value of M_min; at rotation period of 10 ms the minimum mass of neutron stars increases to 0.13 M_sun, and corresponds to the mass-shedding (Keplerian) configuration. In the case of the shortest observed rotation period of radio pulsars 1.56 ms, minimum mass of uniformly rotating cold neutron stars corresponds to the mass-shedding limit, and is found at 0.61 M_sun for the SLy EOS and 0.54 M_sun for the FPS EOS.

An equation of state (EOS) of neutron star matter, describing both the neutron star crust and the liquid core, is calculated. It is based on the effective nuclear interaction SLy of the Skyrme type, which is particularly suitable for the application to the calculation of the properties of very neutron rich matter (Chabanat et al. 1997, 1998). The structure of the crust, and its EOS, is calculated in the T=0 approximation, and under the assumption of the ground state composition. The crust-core transition is a very weakly first-order phase transition, with relative density jump of about one percent. The EOS of the liquid core is calculated assuming (minimal) n-p-e-mu composition. Parameters of static neutron stars are calculated and compared with existing observational data on neutron stars. The minimum and maximum masses of static neutron stars are 0.094 M_sun and 2.05 M_sun, respectively. Effects of rotation on the minimum and the maximum mass of neutron stars are briefly discussed.

The boundary between the solid crust, assumed to be in complete thermodynamic equilibrium (cold catalyzed matter), and the liquid core of a neutron star is studied using Skyrme SLy effective N-N interactions. An approximate value of the density at the inner edge of the crust is obtained from the threshold for the instability of homogeneous $npe$ matter with respect to small periodic density perturbations. Calculations of the bottom layer of the crust are performed within the Compressible Liquid Drop Model. Spherical nuclei are energetically preferred over exotic ones (cylinders, plates, tubes, bubbles), down to the inner edge of the crust, found at the density 0.08 fm^{-3}.

We address the question concerning the maximum density, at which nuclei (and more generally -- nuclear structures) can exist in neutron star interiors. An absolute upper bound to such a density is obtained using the bulk approximation, in which surface and Coulomb effects are neglected. A very good approximation to the limiting density is given by the threshold for the instability of a uniform n-p-e plasma with respect to density modulations; this threshold is calculated using the Extended Thomas-Fermi approximation for the Skyrme energy functionals. For recent SLy Skyrme forces, which are particularly suitable for the description of very neutron rich nucleon systems, the maximum density at which nuclear structures can exist is 0.08 nucleons/fm^3 (i.e., half of the normal nuclear density); at this density protons constitute only 4\% of nucleons.