The shear modulus of the neutron star crust and nonradial oscillations of neutron stars

Shear moduli are calculated for bcc crystalline and rapidly quenched Coulomb solids produced by the Monte Carlo simulation method. The shear moduli are calculated up to the transition temperature and include the effects of thermal fluctuations. An effective shear modulus appropriate to an approximate 'isotropic' body is introduced. It is found that the values of the 'average shear modulus' for the quenched solids remain about the same as those for the corresponding bcc crystals, although the individual shear moduli of the former, disordered solids deviate considerably from the cubic symmetry of the latter. These results are applied to analyses of neutron star oscillations. It is found that the periods of the two interfacial modes are increased by about 10 percent compared to previous results, and that s-mode periods are increased by about 30 percent. The periods of the f and p modes are hardly affected at all. The surface g-mode periods are not greatly affected, while the t-mode periods are increased by 20-25 percent.