and

We present a new method for constraining the Milky Way halo gravitational potential by simul-taneously fitting multiple tidal streams. This method requires full three-dimensional positions and velocities for all stars in the streams, but does not require identification of any specific stream, nor de-termination of stream membership for any star. We exploit the principle that the action distribution of stream stars is most clustered—that is, most informative—when the potential used to calculate the actions is closest to the true potential. We measure the amount of clustering with the Kullback-Leibler Divergence (KLD) or relative entropy, a statistical measure of information which also provides uncertainties for our parameter estimates. We show, for toy Gaia-like data in a spherical isochrone potential, that maximizing the KLD of the action distribution relative to a smoother distribution recovers the true values of the potential parameters. The precision depends on the observational errors and the number and type of streams in the sample; we find that with the phase-space structure and observational uncertainties expected in the Gaia red-giant-star data set, we measure the enclosed mass at the average radius of the sample stars accurate to 3 % and precise to a factor two. Recovery of the scale radius is also precise to roughly a factor two, and is biased 50 % high by the small galac