• ### Models of gravitational lens candidates from Space Warps CFHTLS(1711.07297)

Nov. 21, 2017 astro-ph.GA
We report modelling follow-up of recently-discovered gravitational-lens candidates in the Canada France Hawaii Telescope Legacy Survey. Lens modelling was done by a small group of specially-interested volunteers from the SpaceWarps citizen-science community who originally found the candidate lenses. Models are categorised according to seven diagnostics indicating (a) the image morphology and how clear or indistinct it is, (b) whether the mass map and synthetic lensed image appear to be plausible, and (c) how the lens-model mass compares with the stellar mass and the abundance-matched halo mass. The lensing masses range from ~10^11 Msun to >10^13 Msun. Preliminary estimates of the stellar masses show a smaller spread in stellar mass (except for two lenses): a factor of a few below or above ~10^11 Msun. Therefore, we expect the stellar-to-total mass fraction to decline sharply as lensing mass increases. The most massive system with a convincing model is J1434+522 (SW05). The two low-mass outliers are J0206-095 (SW19) and J2217+015 (SW42); if these two are indeed lenses, they probe an interesting regime of very low star-formation efficiency. Some improvements to the modelling software (SpaghettiLens), and discussion of strategies regarding scaling to future surveys with more and frequent discoveries, are included.
• ### Large Synoptic Survey Telescope Galaxies Science Roadmap(1708.01617)

The Large Synoptic Survey Telescope (LSST) will enable revolutionary studies of galaxies, dark matter, and black holes over cosmic time. The LSST Galaxies Science Collaboration has identified a host of preparatory research tasks required to leverage fully the LSST dataset for extragalactic science beyond the study of dark energy. This Galaxies Science Roadmap provides a brief introduction to critical extragalactic science to be conducted ahead of LSST operations, and a detailed list of preparatory science tasks including the motivation, activities, and deliverables associated with each. The Galaxies Science Roadmap will serve as a guiding document for researchers interested in conducting extragalactic science in anticipation of the forthcoming LSST era.
• ### Maximizing Science in the Era of LSST: A Community-Based Study of Needed US Capabilities(1610.01661)

Oct. 5, 2016 astro-ph.IM
The Large Synoptic Survey Telescope (LSST) will be a discovery machine for the astronomy and physics communities, revealing astrophysical phenomena from the Solar System to the outer reaches of the observable Universe. While many discoveries will be made using LSST data alone, taking full scientific advantage of LSST will require ground-based optical-infrared (OIR) supporting capabilities, e.g., observing time on telescopes, instrumentation, computing resources, and other infrastructure. This community-based study identifies, from a science-driven perspective, capabilities that are needed to maximize LSST science. Expanding on the initial steps taken in the 2015 OIR System Report, the study takes a detailed, quantitative look at the capabilities needed to accomplish six representative LSST-enabled science programs that connect closely with scientific priorities from the 2010 decadal surveys. The study prioritizes the resources needed to accomplish the science programs and highlights ways that existing, planned, and future resources could be positioned to accomplish the science goals.
• ### Gravitational lens modelling in a citizen science context(1502.00008)

Jan. 30, 2015 astro-ph.CO, astro-ph.IM
We develop a method to enable collaborative modelling of gravitational lenses and lens candidates, that could be used by non-professional lens enthusiasts. It uses an existing free-form modelling program (glass), but enables the input to this code to be provided in a novel way, via a user-generated diagram that is essentially a sketch of an arrival-time surface. We report on an implementation of this method, SpaghettiLens, which has been tested in a modelling challenge using 29 simulated lenses drawn from a larger set created for the Space Warps citizen science strong lens search. We find that volunteers from this online community asserted the image parities and time ordering consistently in some lenses, but made errors in other lenses depending on the image morphology. While errors in image parity and time ordering lead to large errors in the mass distribution, the enclosed mass was found to be more robust: the model-derived Einstein radii found by the volunteers were consistent with those produced by one of the professional team, suggesting that given the appropriate tools, gravitational lens modelling is a data analysis activity that can be crowd-sourced to good effect. Ideas for improvement are discussed, these include (a) overcoming the tendency of the models to be shallower than the correct answer in test cases, leading to systematic overestimation of the Einstein radius by 10 per cent at present, and (b) detailed modelling of arcs.
• ### Strong Lens Time Delay Challenge: II. Results of TDC1(1409.1254)

Dec. 12, 2014 astro-ph.CO, astro-ph.IM
We present the results of the first strong lens time delay challenge. The motivation, experimental design, and entry level challenge are described in a companion paper. This paper presents the main challenge, TDC1, which consisted of analyzing thousands of simulated light curves blindly. The observational properties of the light curves cover the range in quality obtained for current targeted efforts (e.g.,~COSMOGRAIL) and expected from future synoptic surveys (e.g.,~LSST), and include simulated systematic errors. \nteamsA\ teams participated in TDC1, submitting results from \nmethods\ different method variants. After a describing each method, we compute and analyze basic statistics measuring accuracy (or bias) $A$, goodness of fit $\chi^2$, precision $P$, and success rate $f$. For some methods we identify outliers as an important issue. Other methods show that outliers can be controlled via visual inspection or conservative quality control. Several methods are competitive, i.e., give $|A|<0.03$, $P<0.03$, and $\chi^2<1.5$, with some of the methods already reaching sub-percent accuracy. The fraction of light curves yielding a time delay measurement is typically in the range $f =$20--40\%. It depends strongly on the quality of the data: COSMOGRAIL-quality cadence and light curve lengths yield significantly higher $f$ than does sparser sampling. Taking the results of TDC1 at face value, we estimate that LSST should provide around 400 robust time-delay measurements, each with $P<0.03$ and $|A|<0.01$, comparable to current lens modeling uncertainties. In terms of observing strategies, we find that $A$ and $f$ depend mostly on season length, while P depends mostly on cadence and campaign duration.
• ### Evidence for the accelerated expansion of the Universe from weak lensing tomography with COSMOS(0911.0053)

March 24, 2010 astro-ph.CO
We present a tomographic cosmological weak lensing analysis of the HST COSMOS Survey. Applying our lensing-optimized data reduction, principal component interpolation for the ACS PSF, and improved modelling of charge-transfer inefficiency, we measure a lensing signal which is consistent with pure gravitational modes and no significant shape systematics. We carefully estimate the statistical uncertainty from simulated COSMOS-like fields obtained from ray-tracing through the Millennium Simulation. We test our pipeline on simulated space-based data, recalibrate non-linear power spectrum corrections using the ray-tracing, employ photometric redshifts to reduce potential contamination by intrinsic galaxy alignments, and marginalize over systematic uncertainties. We find that the lensing signal scales with redshift as expected from General Relativity for a concordance LCDM cosmology, including the full cross-correlations between different redshift bins. For a flat LCDM cosmology, we measure sigma_8(Omega_m/0.3)^0.51=0.75+-0.08 from lensing, in perfect agreement with WMAP-5, yielding joint constraints Omega_m=0.266+0.025-0.023, sigma_8=0.802+0.028-0.029 (all 68% conf.). Dropping the assumption of flatness and using HST Key Project and BBN priors only, we find a negative deceleration parameter q_0 at 94.3% conf. from the tomographic lensing analysis, providing independent evidence for the accelerated expansion of the Universe. For a flat wCDM cosmology and prior w in [-2,0], we obtain w<-0.41 (90% conf.). Our dark energy constraints are still relatively weak solely due to the limited area of COSMOS. However, they provide an important demonstration for the usefulness of tomographic weak lensing measurements from space. (abridged)
• A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.
• ### An Atlas of Predicted Exotic Gravitational Lenses(0904.1454)

Nov. 11, 2009 astro-ph.CO
Wide-field optical imaging surveys will contain tens of thousands of new strong gravitational lenses. Some of these will have new and unusual image configurations, and so will enable new applications: for example, systems with high image multiplicity will allow more detailed study of galaxy and group mass distributions, while high magnification is needed to super-resolve the faintest objects in the high redshift universe. Inspired by a set of six unusual lens systems [including five selected from the Sloan Lens ACS (SLACS) and Strong Lensing Legacy (SL2S) surveys, plus the cluster Abell 1703], we consider several types of multi-component, physically-motivated lens potentials, and use the ray-tracing code "glamroc" to predict exotic image configurations. We also investigate the effects of galaxy source profile and size, and use realistic sources to predict observable magnifications and estimate very approximate relative cross-sections. We find that lens galaxies with misaligned disks and bulges produce swallowtail and butterfly catastrophes, observable as "broken" Einstein rings. Binary or merging galaxies show elliptic umbilic catastrophes, leading to an unusual Y-shaped configuration of 4 merging images. While not the maximum magnification configuration possible, it offers the possibility of mapping the local small-scale mass distribution. We estimate the approximate abundance of each of these exotic galaxy-scale lenses to be ~1 per all-sky survey. In higher mass systems, a wide range of caustic structures are expected, as already seen in many cluster lens systems. We interpret the central ring and its counter-image in Abell 1703 as a "hyperbolic umbilic" configuration, with total magnification ~100 (depending on source size). The abundance of such configurations is also estimated to be ~1 per all-sky survey.
• ### Direct Observation of Cosmic Strings via their Strong Gravitational Lensing Effect: II. Results from the HST/ACS Image Archive(0908.0602)

Aug. 5, 2009 astro-ph.CO
We have searched 4.5 square degrees of archival HST/ACS images for cosmic strings, identifying close pairs of similar, faint galaxies and selecting groups whose alignment is consistent with gravitational lensing by a long, straight string. We find no evidence for cosmic strings in five large-area HST treasury surveys (covering a total of 2.22 square degrees), or in any of 346 multi-filter guest observer images (1.18 square degrees). Assuming that simulations ccurately predict the number of cosmic strings in the universe, this non-detection allows us to place upper limits on the unitless Universal cosmic string tension of G mu/c^2 < 2.3 x 10^-6, and cosmic string density of Omega_s < 2.1 x 10^-5 at the 95% confidence level (marginalising over the other parameter in each case). We find four dubious cosmic string candidates in 318 single filter guest observer images (1.08 square degrees), which we are unable to conclusively eliminate with existing data. The confirmation of any one of these candidates as cosmic strings would imply G mu/c^2 ~ 10^-6 and Omega_s ~ 10^-5. However, we estimate that there is at least a 92% chance that these string candidates are random alignments of galaxies. If we assume that these candidates are indeed false detections, our final limits on G mu/c^2 and Omega_s fall to 6.5 x 10^-7 and 7.3 x 10^-6. Due to the extensive sky coverage of the HST/ACS image archive, the above limits are universal. They are quite sensitive to the number of fields being searched, and could be further reduced by more than a factor of two using forthcoming HST data.
• ### Strong gravitational lensing probes of the particle nature of dark matter(0902.3219)

Feb. 18, 2009 astro-ph.CO
There is a vast menagerie of plausible candidates for the constituents of dark matter, both within and beyond extensions of the Standard Model of particle physics. Each of these candidates may have scattering (and other) cross section properties that are consistent with the dark matter abundance, BBN, and the most scales in the matter power spectrum; but which may have vastly different behavior at sub-galactic "cutoff" scales, below which dark matter density fluctuations are smoothed out. The only way to quantitatively measure the power spectrum behavior at sub-galactic scales at distances beyond the local universe, and indeed over cosmic time, is through probes available in multiply imaged strong gravitational lenses. Gravitational potential perturbations by dark matter substructure encode information in the observed relative magnifications, positions, and time delays in a strong lens. Each of these is sensitive to a different moment of the substructure mass function and to different effective mass ranges of the substructure. The time delay perturbations, in particular, are proving to be largely immune to the degeneracies and systematic uncertainties that have impacted exploitation of strong lenses for such studies. There is great potential for a coordinated theoretical and observational effort to enable a sophisticated exploitation of strong gravitational lenses as direct probes of dark matter properties. This opportunity motivates this white paper, and drives the need for: a) strong support of the theoretical work necessary to understand all astrophysical consequences for different dark matter candidates; and b) tailored observational campaigns, and even a fully dedicated mission, to obtain the requisite data.
• ### Analytic models of plausible gravitational lens potentials(0705.0682)

Dec. 11, 2008 astro-ph
Gravitational lenses on galaxy scales are plausibly modelled as having ellipsoidal symmetry and a universal dark matter density profile, with a Sersic profile to describe the distribution of baryonic matter. Predicting all lensing effects requires knowledge of the total lens potential: in this work we give analytic forms for that of the above hybrid model. Emphasising that complex lens potentials can be constructed from simpler components in linear combination, we provide a recipe for attaining elliptical symmetry in either projected mass or lens potential. We also provide analytic formulae for the lens potentials of Sersic profiles for integer and half-integer index. We then present formulae describing the gravitational lensing effects due to smoothly-truncated universal density profiles in cold dark matter model. For our isolated haloes the density profile falls off as radius to the minus fifth or seventh power beyond the tidal radius, functional forms that allow all orders of lens potential derivatives to be calculated analytically, while ensuring a non-divergent total mass. We show how the observables predicted by this profile differ from that of the original infinite-mass NFW profile. Expressions for the gravitational flexion are highlighted. We show how decreasing the tidal radius allows stripped haloes to be modelled, providing a framework for a fuller investigation of dark matter substructure in galaxies and clusters. Finally we remark on the need for finite mass halo profiles when doing cosmological ray-tracing simulations, and the need for readily-calculable higher order derivatives of the lens potential when studying catastrophes in strong lenses.
• ### Dark Matter and Baryons in the Most X-ray Luminous and Merging Galaxy Cluster RX J1347.5-1145(0711.4850)

March 27, 2008 astro-ph
The galaxy cluster RX J1347-1145 is one of the most X-ray luminous and most massive clusters known. Its extreme mass makes it a prime target for studying issues addressing cluster formation and cosmology. In this paper we present new high-resolution HST/ACS and Chandra X-ray data. The high resolution and sensitivity of ACS enabled us to detect and quantify several new multiply imaged sources, we now use a total of eight for the strong lensing analysis. Combining this information with shape measurements of weak lensing sources in the central regions of the cluster, we derive a high-resolution, absolutely-calibrated mass map. This map provides the best available quantification of the total mass of the central part of the cluster to date. We compare the reconstructed mass with that inferred from the new Chandra X-ray data, and conclude that both mass estimates agree extremely well in the observed region, namely within 400 / h_70 kpc of the cluster center. In addition we study the major baryonic components (gas and stars) and hence derive the dark matter distribution in the center of the cluster. We find that the dark matter and baryons are both centered on the BCG within the uncertainties (alignment is better than <10 kpc). We measure the corresponding 1-D profiles and find that dark matter distribution is consistent with both NFW and cored profiles, indicating that a more extended radial analysis is needed to pinpoint the concentration parameter, and hence the inner slope of the dark matter profile.
• ### AMUSE-Virgo I. Super-massive black holes in low-mass spheroids(0711.2073)

March 4, 2008 astro-ph
We present the first results from the AGN Multiwavelength Survey of Early-type galaxies in the Virgo cluster (AMUSE-Virgo). This large program targets 100 early-type galaxies with the Advanced CCD Imaging Spectrometer on board the Chandra X-ray Observatory and the Multi-band Imaging Photometer on board the Spitzer Space Telescope, with the aim of providing an unbiased census of low-level super-massive black hole (SMBH) activity in the local universe. Here we report on the Chandra observations of the first 16 targets, and combine them with results from archival data of another, typically more massive, 16 targets. Point-like X-ray emission from a position coincident with the optical nucleus is detected in 50% of the galaxies (down to our completeness limit of ~4E+38 erg/sec). Two of the X-ray nuclei are hosted by galaxies (VCC1178=N4464 and VCC1297=N4486B) with absolute B magnitudes fainter than -18, where nuclear star clusters are known to become increasingly common. After carefully accounting for possible contamination from low mass X-ray binaries, we argue that the detected nuclear X-ray sources are most likely powered by low-level accretion on to a SMBH, with a <11% chance contamination in VCC1178, where a star cluster is barely resolvable in archival Hubble Space Telescope images. Based on black hole mass estimates from the global properties of the host galaxies, all the detected nuclei are highly sub-Eddington, with luminosities in the range -8.4<log(L_0.3-10keV/L_Edd)<-5.9. The incidence of nuclear X-ray activity increases with the stellar mass M_star of the host galaxy: only between 3-44% of the galaxies with M_star<1E+10 M_Sun harbor an X-ray active SMBH. The fraction rises to between 49-87% in galaxies with stellar mass above 1E+10 M_Sun (at the 95% confidence level).
• ### A Bayesian approach to strong lensing modelling of galaxy clusters(0706.0048)

Dec. 26, 2007 astro-ph
In this paper, we describe a procedure for modelling strong lensing galaxy clusters with parametric methods, and to rank models quantitatively using the Bayesian evidence. We use a publicly available Markov chain Monte-Carlo (MCMC) sampler ('Bayesys'), allowing us to avoid local minima in the likelihood functions. To illustrate the power of the MCMC technique, we simulate three clusters of galaxies, each composed of a cluster-scale halo and a set of perturbing galaxy-scale subhalos. We ray-trace three light beams through each model to produce a catalogue of multiple images, and then use the MCMC sampler to recover the model parameters in the three different lensing configurations. We find that, for typical Hubble Space Telescope (HST)-quality imaging data, the total mass in the Einstein radius is recovered with ~1-5% error according to the considered lensing configuration. However, we find that the mass of the galaxies is strongly degenerated with the cluster mass when no multiple images appear in the cluster centre. The mass of the galaxies is generally recovered with a 20% error, largely due to the poorly constrained cut-off radius. Finally, we describe how to rank models quantitatively using the Bayesian evidence. We confirm the ability of strong lensing to constrain the mass profile in the central region of galaxy clusters in this way. Ultimately, such a method applied to strong lensing clusters with a very large number of multiple images may provide unique geometrical constraints on cosmology. The implementation of the MCMC sampler used in this paper has been done within the framework of the Lenstool software package, which is publicly available.
• ### Direct Observation of Cosmic Strings via their Strong Gravitational Lensing Effect: I. Predictions for High Resolution Imaging Surveys(0710.5544)

Dec. 21, 2007 astro-ph
We use current theoretical estimates for the density of long cosmic strings to predict the number of strong gravitational lensing events in astronomical imaging surveys as a function of angular resolution and survey area. We show that angular resolution is the most important factor, and that interesting limits on the dimensionless string tension Gmu/c^2 can be obtained by existing and planned surveys. At the resolution of the Hubble Space Telescope (0.14"), it is sufficient to survey of order a few square degrees -- well within reach of the current HST archive -- to probe the regime Gmu/c^2 ~ 10^{-7}. If lensing by cosmic strings is not detected, such a survey would improve the limit on the string tension by a factor of two over that available from the cosmic microwave background. Future high resolution imaging surveys, covering a few hundred square degrees or more, either from space in the optical or from large-format radio telescopes on the ground, would be able to further lower this limit to Gmu/c^2 < 10^{-8}.
• ### Physical component analysis of galaxy cluster weak gravitational lensing data(astro-ph/0511287)

Oct. 30, 2007 astro-ph
We present a novel approach for reconstructing the projected mass distribution of clusters of galaxies from sparse and noisy weak gravitational lensing shear data. The reconstructions are regularised using knowledge gained from numerical simulations of clusters: trial mass distributions are constructed from N physically-motivated components, each of which has the universal density profile and characteristic geometry observed in simulated clusters. The parameters of these components are assumed to be distributed \emph{a priori} in the same way as they are in the simulated clusters. Sampling mass distributions from the components' parameters' posterior probability density function allows estimates of the mass distribution to be generated, with error bars. The appropriate number of components is inferred from the data itself via the Bayesian evidence, and is typically found to be small, reflecting the quality of the simulated data used in this work. Ensemble average mass maps are found to be robust to the details of the noise realisation, and succeed in recovering the input mass distribution (from a realistic simulated cluster) over a wide range of scales. We comment on the residuals of the reconstruction and their implications, and discuss the extension of the method to include strong lensing information.
• ### Bayesian evidence as a tool for comparing datasets(astro-ph/0412535)

Oct. 30, 2007 astro-ph
We introduce a new conservative test for quantifying the consistency of two or more datasets. The test is based on the Bayesian answer to the question, How much more probable is it that all my data were generated from the same model system than if each dataset were generated from an independent set of model parameters?''. We make explicit the connection between evidence ratios and the differences in peak chi-squared values, the latter of which are more widely used and more cheaply calculated. Calculating evidence ratios for three cosmological datasets (recent CMB data (WMAP, ACBAR, CBI, VSA), SDSS and the most recent SNe Type 1A data) we find that concordance is favoured and the tightening of constraints on cosmological parameters is indeed justified.
• ### Strong and weak lensing united III: Measuring the mass distribution of the merging galaxy cluster 1E0657-56(astro-ph/0608408)

Aug. 18, 2006 astro-ph, hep-ph, gr-qc
The galaxy cluster 1E0657-56 (z = 0.296) is remarkably well-suited for addressing outstanding issues in both galaxy evolution and fundamental physics. We present a reconstruction of the mass distribution from both strong and weak gravitational lensing data. Multi-color, high-resolution HST ACS images allow detection of many more arc candidates than were previously known, especially around the subcluster. Using the known redshift of one of the multiply imaged systems, we determine the remaining source redshifts using the predictive power of the strong lens model. Combining this information with shape measurements of "weakly" lensed sources, we derive a high-resolution, absolutely-calibrated mass map, using no assumptions regarding the physical properties of the underlying cluster potential. This map provides the best available quantification of the total mass of the central part of the cluster. We also confirm the result from Clowe et al. (2004,2006a).
• ### The SNAP Strong Lens Survey(astro-ph/0501328)

Jan. 17, 2005 astro-ph
Basic considerations of lens detection and identification indicate that a wide field survey of the types planned for weak lensing and Type Ia SNe with SNAP are close to optimal for the optical detection of strong lenses. Such a piggy-back'' survey might be expected even pessimistically to provide a catalogue of a few thousand new strong lenses, with the numbers dominated by systems of faint blue galaxies lensed by foreground ellipticals. After sketching out our strategy for detecting and measuring these galaxy lenses using the SNAP images, we discuss some of the scientific applications of such a large sample of gravitational lenses: in particular we comment on the partition of information between lens structure, the source population properties and cosmology. Understanding this partitioning is key to assessing strong lens cosmography's value as a cosmological probe.
• ### Very Small Array observations of the Sunyaev-Zel'dovich effect in nearby galaxy clusters(astro-ph/0405582)

Dec. 2, 2004 astro-ph
We present VSA observations (~34GHz) on scales ~20 arcmin towards a complete, X-ray-flux-limited sample of seven clusters at redshift z<0.1. Four have significant SZ detections in the presence of CMB primordial anisotropy. We use a bayesian MCMC method for inference from the VSA data, with X-ray priors on cluster positions and temperatures, and radio priors on sources. We make assumptions of beta-model gas distributions and of hydrostatic equilibrium, to evaluate probability densities for the gas mass and total mass out to r_200. Our combined estimate of the gas fraction is 0.08^{+0.06}_{-0.04}h^{-1} The random errors are poor (note that the errors are higher than would have been obtained with the usual chi-squared method) but the control of bias is good. We have described the MCMC analysis method specifically in terms of SZ but hope the description will be of more general use. We find that the effects of primordial CMB contamination tend to be similar in the estimates of both the gas mass and total mass over our narrow range of angular scales, so that there is little effect of primordials on the gas fraction determination. Using our total mass estimates we find a normalisation of the mass-temperature relation based on the profiles from the VSA cluster pressure maps that is in good agreement with recent M-T determinations from X-ray cluster measurements.