Free Form Lensing Implications for the Collision of Dark Matter and Gas in the Frontier Fields Cluster MACSJ0416.1-2403

We present a free form mass reconstruction of the massive lensing cluster MACSJ0416.1-2403 using the latest Hubble Frontier Fields data. Our model independent method finds that the extended lensing pattern is generated by two elongated, closely projected clusters of similar mass. Our lens model identifies new lensed images with which we improve the accuracy of the dark matter distribution. We find that the bimodal mass distribution is nearly coincident with the bimodal X-ray emission, but with the two dark matter peaks lying closer together than the centroids of the X-ray emisison. We show this can be achieved if the collision has occurred close to the plane and such that the cores are deflected around each other. The projected mass profiles of both clusters are well constrained because of the many interior lensed images, leading to surprisingly flat mass profiles of both components in the region 15-100 kpc. We discuss the extent to which this may be generated by tidal forces in our dynamical model which are large during an encounter of this type as the cores "graze" each other. The relative velocity between the two cores is estimated to be about 1200 km/s and mostly along the line of sight so that our model is consistent with the relative redshift difference between the two cD galaxies (dz = 0.04).Comment: 22 pages, 18 figures, 2 table

Dimensional renormalization: ladders to rainbows

Renormalization factors are most easily extracted by going to the massless limit of the quantum field theory and retaining only a single momentum scale. We derive factors and renormalized Green functions to all orders in perturbation theory for rainbow graphs and vertex (or scattering diagrams) at zero momentum transfer, in the context of dimensional renormalization, and we prove that the correct anomalous dimensions for those processes emerge in the limit D -> 4.Comment: RevTeX, no figure

Master integrals for massive two-loop Bhabha scattering in QED

We present a set of scalar master integrals (MIs) needed for a complete treatment of massive two-loop corrections to Bhabha scattering in QED, including integrals with arbitrary fermionic loops. The status of analytical solutions for the MIs is reviewed and examples of some methods to solve MIs analytically are worked out in more detail. Analytical results for the pole terms in epsilon of so far unknown box MIs with five internal lines are given.Comment: 23 pages, 5 tables, 12 figures, references added, appendix B enlarge

The Surprisingly Steep Mass Profile of Abell 1689, from a Lensing Analysis of Subaru Images

Subaru observations of A1689 (z=0.183) are used to derive an accurate, model-independent mass profile for the entire cluster, r<2 Mpc/h, by combining magnification bias and distortion measurements. The projected mass profile steepens quickly with increasing radius, falling away to zero at r~1.0 Mpc/h, well short of the anticipated virial radius. Our profile accurately matches onto the inner profile, r<200 kpc/h, derived from deep HST/ACS images. The combined ACS and Subaru information is well fitted by an NFW profile with virial mass, (1.93 \pm 0.20)10^15 M_sun, and surprisingly high concentration, c_vir=13.7^{+1.4}_{-1.1}, significantly larger than theoretically expected (c_vir~4), corresponding to a relatively steep overall profile. A slightly better fit is achieved with a steep power-law model that has its 2D logarithmic slope -3 and core radius theta_c~1.7' (r_c~210 kpc/h), whereas an isothermal profile is strongly rejected. These results are based on a reliable sample of background galaxies selected to be redder than the cluster E/S0 sequence. By including the faint blue galaxy population a much smaller distortion signal is found, demonstrating that blue cluster members significantly dilute the true signal for r~400 kpc/h. This contamination is likely to affect most weak lensing results to date.Comment: 5 pages, 3 figures, to appear in ApJ

Lectures on multiloop calculations

I discuss methods of calculation of propagator diagrams (massless, those of Heavy Quark Effective Theory, and massive on-shell diagrams) up to 3 loops. Integration-by-parts recurrence relations are used to reduce them to linear combinations of basis integrals. Non-trivial basis integrals have to be calculated by some other method, e.g., using Gegenbauer polynomial technique. Many of them are expressed via hypergeometric functions; in the massless and HQET cases, their indices tend to integers at $\epsilon\to0$. I discuss the algorithm of their expansion in $\epsilon$, in terms of multiple $\zeta$ values. These lectures were given at Calc-03 school, Dubna, 14--20 June 2003.Comment: 52 pages, 49 figures. Lectures at Calc-03 school, Dubna, 14--20 June 2003. v2: 2 references added, minor typos corrected. v3: methodical improvements, typo in eq. (3.19) corrected, 2 references adde

Power Spectrum of Velocity Fluctuations in the Universe

We investigate the power spectrum of velocity fluctuations in the universe, $V^2(k)$, starting from four different measures of velocity: (1) the power spectrum of velocity fluctuations from peculiar velocities of galaxies; (2) the rms peculiar velocity of galaxy clusters; (3) the power spectrum of velocity fluctuations from the power spectrum of density fluctuations in the galaxy distribution; (4) and the bulk velocity from peculiar velocities of galaxies. We show that measures (1) and (2) are not consistent with each other and either the power spectrum from peculiar velocities of galaxies is overestimated or the rms cluster peculiar velocity is underestimated. The amplitude of velocity fluctuations derived from the galaxy distribution (measure 3) depends on the parameter $\beta$. We estimate the parameter $\beta$ on the basis of measures (2) and (4). The power spectrum of velocity fluctuations from the galaxy distribution in the Stromlo-APM redshift survey is consistent with the observed rms cluster velocity and with the observed large-scale bulk flow when the parameter $\beta$ is in the range 0.4-0.5. In this case the value of the function $V(k)$ at wavelength $\lambda=120h^{-1}$Mpc is $\sim 350$ km s$^{-1}$ and the rms amplitude of the bulk flow at the radius $r=60h^{-1}$ Mpc is $\sim 340$ km s$^{-1}$. The velocity dispersion of galaxy systems originates mostly from the large-scale velocity fluctuations with wavelengths $\lambda >100h^{-1}$ Mpc.Comment: Astrophysical Journal, Vol. 493, in press: 23 pages, uses AAS Latex, and 14 separate postscript figure

BIMA and Keck Imaging of the Radio Ring PKS 1830-211

We discuss BIMA (Berkeley Illinois Maryland Association) data and present new high quality optical and near-IR Keck images of the bright radio ring PKS 1830-211. Applying a powerful new deconvolution algorithm we have been able to identify both images of the radio source. In addition we recover an extended source in the optical, consistent with the expected location of the lensing galaxy. The source counterparts are very red, I-K=7, suggesting strong Galactic absorption with additional absorption by the lensing galaxy at z=0.885, and consistent with the detection of high redshift molecules in the lens.Comment: To be published in the ASP Conference Proceedings, 'Highly Redshifted Radio Lines', Greenbank, W

Free-form lens model and mass estimation of the high redshift galaxy cluster ACT-CL J0102-4915, "El Gordo"

We examine the massive colliding cluster El Gordo, one of the most massive clusters at high redshift. We use a free-form lensing reconstruction method that avoids making assumptions about the mass distribution. We use data from the RELICS program and identify new multiply lensed system candidates. The new set of constraints and free-form method provides a new independent mass estimate of this intriguing colliding cluster. Our results are found to be consistent with earlier parametric models, indirectly confirming the assumptions made in earlier work. By fitting a double gNFW profile to the lens model, and extrapolating to the virial radius, we infer a total mass for the cluster of $M_{200c}=(1.08^{+0.65}_{-0.12})\times10^{15}$M$_{\odot}$. We estimate the uncertainty in the mass due to errors in the photometric redshifts, and discuss the uncertainty in the inferred virial mass due to the extrapolation from the lens model. We also find in our lens map a mass overdensity corresponding to the large cometary tail of hot gas, reinforcing its interpretation as a large tidal feature predicted by hydrodynamical simulations that mimic El Gordo. Finally, we discuss the observed relation between the plasma and the mass map, finding that the peak in the projected mass map may be associated with a large concentration of colder gas, exhibiting possible star formation. El Gordo is one of the first clusters that will be observed with JWST, which is expected to unveil new high redshift lensed galaxies around this interesting cluster, and provide a more accurate estimation of its mass.Comment: 19 pages, 10 figures. Updated figure

Analytical and numerical methods for massive two-loop self-energy diagrams

Motivated by the precision results in the electroweak theory studies of two-loopFeynman diagrams are performed. Specifically this paper gives a contribution to the knowledge of massive two-loop self-energy diagrams in arbitrary and especially four dimensions.This is done in three respects firstly results in terms of generalized, multivariable hypergeometric functions are presented giving explicit series for small and large momenta. Secondly the imaginary parts of these integrals are expressed as complete elliptic integrals.Finally one-dimensional integral representations with elementary functions are derived.They are very well suited for the numerical evaluations.Comment: 24 page

Cloning Hubble Deep Fields I: A Model-Independent Measurement of Galaxy Evolution

We present a model-independent method of quantifying galaxy evolution in high-resolution images, which we apply to the Hubble Deep Field (HDF). Our procedure is to k-correct all pixels belonging to the images of a complete set of bright galaxies and then to replicate each galaxy image to higher redshift by the product of its space density, 1/V_{max}, and the cosmological volume. The set of bright galaxies is itself selected from the HDF, because presently the HDF provides the highest quality UV images of a redshift-complete sample of galaxies (31 galaxies with I<21.9, \bar{z}=0.5, and for which V/V_{max} is spread fairly). These galaxies are bright enough to permit accurate pixel-by-pixel k-corrections into the restframe UV (\sim 2000 A). We match the shot noise, spatial sampling and PSF smoothing of the HDF data, resulting in entirely empirical and parameter-free ``no-evolution'' deep fields of galaxies for direct comparison with the HDF. In addition, the overcounting rate and the level of incompleteness can be accurately quantified by this procedure. We obtain the following results. Faint HDF galaxies (I>24) are much smaller, more numerous, and less regular than our ``no-evolution'' extrapolation, for any interesting geometry. A higher proportion of HDF galaxies ``dropout'' in both U and B, indicating that some galaxies were brighter at higher redshifts than our ``cloned'' z\sim0.5 population.Comment: 51 pages, 23 figures, replacement includes figures not previously include