It has been clear for some time now that super-critical surface magnetic fields, exceeding 4 x 10^13 G, exist on a subset of neutron stars. These magnetars may harbor interior fields many orders of magnitude larger, potentially reaching equipartition values. However, the impact of these strong fields on stellar structure has been largely ignored, potentially complicating attempts to infer the high density nuclear equation of state. Here we assess the effect of these strong magnetic fields on the mass-radius relationship of neutron stars. We employ an effective field theory model for the nuclear equation of state that includes the impact of hyperons, anomalous magnetic moments, and the physics of the crust. We consider two magnetic field geometries, bounding the likely magnitude of the impact of magnetic fields: a statistically isotropic, tangled field and a force-free configuration. In both cases even equipartition fields have at most a 30% impact on the maximum mass. However, the direction of the effect of the magnetic field depends on the geometry employed - force-free fields leading to reductions in the maximum neutron star mass and radius while tangled fields increase both - challenging the common intuition in the literature on the impact of magnetic fields.

The gravitational aether theory is a modification of general relativity that decouples vacuum energy from gravity, and thus can potentially address the cosmological constant problem. The classical theory is distinguishable from general relativity only in the presence of relativistic pressure (or vorticity). Since the interior of neutron stars has high pressure and as their mass and radius can be measured observationally, they are the perfect laboratory for testing the validity of the aether theory. In this paper, we solve the equations of stellar structure for the gravitational aether theory and find the predicted mass-radius relation of non-rotating neutron stars using two different realistic proposals for the equation of state of nuclear matter. We find that the maximum neutron star mass predicted by the aether theory is 12% - 16% less than the maximum mass predicted by general relativity assuming these two equations of state. We also show that the effect of aether is similar to modifying the equation of state in general relativity. The effective pressure of the neutron star given by the aether theory at a fiducial density differs from the values given by the two nuclear equations of state to an extent that can be constrained using future gravitational wave observations of neutron stars in compact systems. This is a promising way to test the aether theory if further progress is made in constraining the equation of state of nuclear matter in densities above the nuclear saturation density.