In the late inspiral phase, gravitational waves from binary neutron star mergers carry the imprint of the equation of state due to the tidally deformed structure of the components. If the stars contain solid crusts, then their shear modulus can affect the deformability of the star and, thereby, modify the emitted signal. Here, we investigate the effect of realistic equations of state (EOSs) of the crustal matter, with a realistic model for the shear modulus of the stellar crust in a fully general relativistic framework. This allows us to systematically study the deviations that are expected from fluid models. In particular, we use unified EOSs, both relativistic and non-relativistic, in our calculations. We find that realistic EOSs of crusts cause a small correction, of $\sim 1\%$, in the second Love number. This correction will likely be subdominant to the statistical error expected in LIGO-Virgo observations at their respective advanced design sensitivities, but rival that error in third generation detectors. For completeness, we also study the effect of crustal shear on the magnetic-type Love number and find it to be much smaller.

We use gravitational wave (GW) and electromagnetic (EM) observations of GW170817 to constrain the extent of pressure anisotropy in it. While it is quite likely that the pressure inside a neutron star is mostly isotropic, certain physical processes or characteristics, such as phase transitions in nuclear matter or the presence of strong magnetic fields, can introduce pressure anisotropy. In this work, we show that anisotropic pressure in neutron stars can reduce their tidal deformability substantially. For the anisotropy-pressure model of Bowers and Liang and a couple of relativistic EOSs -- DDH$\delta$ and GM1 -- we demonstrate that this reduction in spherical neutron stars with masses in the range of 1 to 2 $M_\odot$ can be 23% to 46%. This suggests that certain EOSs that are ruled out by GW170817 observations, under assumptions of pressure isotropy, can become viable if the stars had a significant enough anisotropic pressure component, but do not violate causality. We also show how the inference of the star radius can be used to rule out certain EOSs (such as GM1), even for high anisotropic pressure, because their radii are larger than what the observations find.