We present a method of calculating the strong-field gravitational lensing caused by many analytic and numerical spacetimes. We use this procedure to calculate the distortion caused by isolated black holes and by numerically evolved black hole binaries. We produce both demonstrative images illustrating details of the spatial distortion and realistic images of collections of stars taking both lensing amplification and redshift into account. On large scales the lensing from inspiraling binaries resembles that of single black holes, but on small scales the resulting images show complex and in some cases self-similar structure across different angular scales.

We develop the formalism for computing the magnetic field within an axisymmetric neutron star with a strong Type II superconductor core surrounded by a normal conductor. The formalism takes full account of the constraints imposed by hydrostatic equilibrium with a barotropic equation of state. We specialize to purely poloidal magnetic fields and develop the "most dipolar case" for which we find that the surface field strength is $\simeq H_b\epsilon_b/3\simeq 3\times 10^{12}$ G, where $H_b$ is the magnetic field at the outer edge of the core and $\epsilon_b R$ is the thickness of the normal shell. This value only depends on the equation of state of nuclear matter. We also find that the quadrupolar distortion of the star is $\sim 10^{-9}$.

For highly compact neutron stars, constructing numerical initial data for black hole-neutron star binary evolutions is very difficult. We describe improvements to an earlier method that enable it to handle these more challenging cases. We examine the case of a 6:1 mass ratio system in inspiral close to merger, where the star is governed by a polytropic $\Gamma=2$, an SLy, or an LS220 equation of state. In particular, we are able to obtain a solution with a realistic LS220 equation of state for a star with compactness 0.26 and mass 1.98 $M_\odot$, which is representative of the highest reliably determined neutron star masses. For the SLy equation of state, we can obtain solutions with a comparable compactness of 0.25, while for a family of polytropic equations of state, we obtain solutions with compactness up to 0.21, the largest compactness that is stable in this family. These compactness values are significantly higher than any previously published results. We find that improvements in adapting the computational domain to the neutron star surface and in accounting for the center of mass drift of the system are the key ingredients allowing us to obtain these solutions.