WEIGHING EL GORDO WITH A PRECISION SCALE: HUBBLE SPACE TELESCOPE WEAK-LENSING ANALYSIS OF THE MERGING GALAXY CLUSTER ACT-CL J0102-4915 AT z=0.87

We present a Hubble Space Telescope weak-lensing study of the merging galaxy cluster El Gordo (ACT-CL J0102-4915) at z = 0.87 discovered by the Atacama Cosmology Telescope (ACT) collaboration as the strongest Sunyaev-Zel'dovich decrement in its similar to 1000 deg(2) survey. Our weak-lensing analysis confirms that ACT-CL J0102-4915 is indeed an extreme system consisting of two massive (greater than or similar to 10(15) M-circle dot each) subclusters with a projected separation of similar to 0.7 h(70)(-1) Mpc. This binary mass structure revealed by our lensing study is consistent with the cluster galaxy distribution and the dynamical study carried out with 89 spectroscopic members. We estimate the mass of ACT-CL J0102-4915 by simultaneously fitting two axisymmetric Navarro-Frenk-White (NFW) profiles allowing their centers to vary. We use only a single parameter for the NFW mass profile by enforcing the mass-concentration relation from numerical simulations. Our Markov-Chain-Monte-Carlo analysis shows that the masses of the northwestern (NW) and the southeastern (SE) components are M-200c = (1.38 +/- 0.22) Chi 10(15) h(70)(-1) M-circle dot and (0.78 +/- 0.20) Chi 10(15) h(70)(-1) M-circle dot, respectively, where the quoted errors include only 1 sigma statistical uncertainties determined by the finite number of source galaxies. These mass estimates are subject to additional uncertainties (20%-30%) due to the possible presence of triaxiality, correlated/uncorrelated large scale structure, and departure of the cluster profile from the NFW model. The lensing-based velocity dispersions are 1133(-61)(+58) km s(-1) and 1064(-66)(+62) km s(-1) for the NW and SE components, respectively, which are consistent with their spectroscopic measurements (1290 +/- 134 km s(-1) and 1089 +/- 200 km s(-1), respectively). The centroids of both components are tightly constrained (similar to 4 '') and close to the optical luminosity centers. The X-ray and mass peaks are spatially offset by similar to 8 '' (similar to 62 h(70)(-1) kpc), which is significant at the similar to 2 sigma confidence level. The mass peak, however, does not lead the gas peak in the direction expected if we are viewing the cluster soon after first core passage during a high speed merger. Under the assumption that the merger is happening in the plane of the sky, extrapolation of the two NFW halos to a radius r(200a) = 2.4 h(70)(1) Mpc yields a combined mass of M-200a = (3.13 +/- 0.56) Chi 10(15) h(70)(1) M-circle dot. This extrapolated total mass is consistent with our two-component-based dynamical analysis and previous X-ray measurements, projecting ACT-CL J0102-4915 to be the most massive cluster at z > 0.6 known to date.