Multiple image gravitational lens systems, and especially quads are invaluable in determining the amount and distribution of mass in galaxies. This is usually done by mass modeling using parametric or free-form methods. An alternative way of extracti
ng information about lens mass distribution is to use lensing degeneracies and invariants. Where applicable, they allow one to make conclusions about whole classes of lenses without model fitting. Here, we use approximate, but observationally useful invariants formed by the three relative polar angles of quad images around the lens center to show that many smooth elliptical+shear lenses can reproduce the same set of quad image angles within observational error. This result allows us to show in a model-free way what the general class of smooth elliptical+shear lenses looks like in the three dimensional (3D) space of image relative angles, and that this distribution does not match that of the observed quads. We conclude that, even though smooth elliptical+shear lenses can reproduce individual quads, they cannot reproduce the quad population. What is likely needed is substructure, with clump masses larger than those responsible for flux ratio anomalies in quads, or luminous or dark nearby perturber galaxies.
It has been shown in previous work that DARKexp, which is a theoretically derived, maximum entropy, one shape parameter model for isotropic collisionless systems, provides very good fits to simulated and observed dark-matter halos. Specifically, it f
its the energy distribution, N(E), and the density profiles, including the central cusp. Here, we extend DARKexp N(E) to include the distribution in angular momentum, L^2, for spherically symmetric systems. First, we argue, based on theoretical, semi-analytical, and simulation results, that while dark-matter halos are relaxed in energy, they are not nearly as relaxed in angular momentum, which precludes using maximum entropy to uniquely derive N(E,L^2). Instead, we require that when integrating N(E,L^2) over squared angular momenta one retrieves the DARKexp N(E). Starting with a general expression for N(E,L^2) we show how the distribution of particles in L^2 is related to the shape of the velocity distribution function, VDF, and velocity anisotropy profile, beta(r). We then demonstrate that astrophysically realistic halos, as judged by the VDF shape and beta(r), must have linear or convex distributions in L^2, for each separate energy bin. The distribution in energy of the most bound particles must be nearly flat, and become more tilted in favor of radial orbits for less bound particles. These results are consistent with numerical simulations and represent an important step towards deriving the full distribution function for spherically symmetric dark-matter halos.
