Source author record

Philip J. Marshall

Philip J. Marshall appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

39works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

39 published item(s)

preprint2022arXiv

GREAT3 results I: systematic errors in shear estimation and the impact of real galaxy morphology

We present first results from the third GRavitational lEnsing Accuracy Testing (GREAT3) challenge, the third in a sequence of challenges for testing methods of inferring weak gravitational lensing shear distortions from simulated galaxy images. GREAT3 was divided into experiments to test three specific questions, and included simulated space- and ground-based data with constant or cosmologically-varying shear fields. The simplest (control) experiment included parametric galaxies with a realistic distribution of signal-to-noise, size, and ellipticity, and a complex point spread function (PSF). The other experiments tested the additional impact of realistic galaxy morphology, multiple exposure imaging, and the uncertainty about a spatially-varying PSF; the last two questions will be explored in Paper II. The 24 participating teams competed to estimate lensing shears to within systematic error tolerances for upcoming Stage-IV dark energy surveys, making 1525 submissions overall. GREAT3 saw considerable variety and innovation in the types of methods applied. Several teams now meet or exceed the targets in many of the tests conducted (to within the statistical errors). We conclude that the presence of realistic galaxy morphology in simulations changes shear calibration biases by $\sim 1$ per cent for a wide range of methods. Other effects such as truncation biases due to finite galaxy postage stamps, and the impact of galaxy type as measured by the Sérsic index, are quantified for the first time. Our results generalize previous studies regarding sensitivities to galaxy size and signal-to-noise, and to PSF properties such as seeing and defocus. Almost all methods' results support the simple model in which additive shear biases depend linearly on PSF ellipticity.

preprint2020arXiv

Time Delay Lens Modeling Challenge: I. Experimental Design

Strong gravitational lenses with measured time delay are a powerful tool to measure cosmological parameters, especially the Hubble constant ($H_0$). Recent studies show that by combining just three multiply-imaged AGN systems, one can determine $H_0$ to 2.4% precision. Furthermore, the number of time-delay lens systems is growing rapidly, enabling the determination of $H_0$ to 1% precision in the near future. However, as the precision increases it is important to ensure that systematic errors and biases remain subdominant. For this purpose, challenges with simulated datasets are a key component in this process. Following the experience of the past challenge on time delay, where it was shown that time delays can indeed be measured precisely and accurately at the sub-percent level, we now present the "Time Delay Lens Modeling Challenge" (TDLMC). The goal of this challenge is to assess the present capabilities of lens modeling codes and assumptions and test the level of accuracy of inferred cosmological parameters given realistic mock datasets. We invite scientists to model a set of simulated HST observations of 50 mock lens systems. The systems are organized in rungs, with the complexity and realism increasing going up the ladder. The goal of the challenge is to infer $H_0$ for each rung, given the HST images, the time delay, and stellar velocity dispersion of the deflector for a fixed background cosmology. The TDLMC challenge starts with the mock data release on 2018 January 8th. The deadline for blind submission is different for each rung. The deadline for Rung 0-1 is 2018 September 8; the deadline for Rung 2 is 2019 April 8 and the one for Rung 3 is 2019 September 8. This first paper gives an overview of the challenge including the data design, and a set of metrics to quantify the modeling performance and challenge details.

preprint2019arXiv

H0LiCOW XII. Lens mass model of WFI2033-4723 and blind measurement of its time-delay distance and $H_0$

We present the lens mass model of the quadruply-imaged gravitationally lensed quasar WFI2033-4723, and perform a blind cosmographical analysis based on this system. Our analysis combines (1) time-delay measurements from 14 years of data obtained by the COSmological MOnitoring of GRAvItational Lenses (COSMOGRAIL) collaboration, (2) high-resolution $\textit{Hubble Space Telescope}$ imaging, (3) a measurement of the velocity dispersion of the lens galaxy based on ESO-MUSE data, and (4) multi-band, wide-field imaging and spectroscopy characterizing the lens environment. We account for all known sources of systematics, including the influence of nearby perturbers and complex line-of-sight structure, as well as the parametrization of the light and mass profiles of the lensing galaxy. After unblinding, we determine the effective time-delay distance to be $4784_{-248}^{+399}~\mathrm{Mpc}$, an average precision of $6.6\%$. This translates to a Hubble constant $H_{0} = 71.6_{-4.9}^{+3.8}~\mathrm{km~s^{-1}~Mpc^{-1}}$, assuming a flat $Λ$CDM cosmology with a uniform prior on $Ω_\mathrm{m}$ in the range [0.05, 0.5]. This work is part of the $H_0$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) collaboration, and the full time-delay cosmography results from a total of six strongly lensed systems are presented in a companion paper (H0LiCOW XIII).

preprint2019arXiv

H0LiCOW XIII. A 2.4% measurement of $H_{0}$ from lensed quasars: $5.3σ$ tension between early and late-Universe probes

We present a measurement of the Hubble constant ($H_{0}$) and other cosmological parameters from a joint analysis of six gravitationally lensed quasars with measured time delays. All lenses except the first are analyzed blindly with respect to the cosmological parameters. In a flat $Λ$CDM cosmology, we find $H_{0} = 73.3_{-1.8}^{+1.7}$, a 2.4% precision measurement, in agreement with local measurements of $H_{0}$ from type Ia supernovae calibrated by the distance ladder, but in $3.1σ$ tension with $Planck$ observations of the cosmic microwave background (CMB). This method is completely independent of both the supernovae and CMB analyses. A combination of time-delay cosmography and the distance ladder results is in $5.3σ$ tension with $Planck$ CMB determinations of $H_{0}$ in flat $Λ$CDM. We compute Bayes factors to verify that all lenses give statistically consistent results, showing that we are not underestimating our uncertainties and are able to control our systematics. We explore extensions to flat $Λ$CDM using constraints from time-delay cosmography alone, as well as combinations with other cosmological probes, including CMB observations from $Planck$, baryon acoustic oscillations, and type Ia supernovae. Time-delay cosmography improves the precision of the other probes, demonstrating the strong complementarity. Allowing for spatial curvature does not resolve the tension with $Planck$. Using the distance constraints from time-delay cosmography to anchor the type Ia supernova distance scale, we reduce the sensitivity of our $H_0$ inference to cosmological model assumptions. For six different cosmological models, our combined inference on $H_{0}$ ranges from $\sim73$-$78~\mathrm{km~s^{-1}~Mpc^{-1}}$, which is consistent with the local distance ladder constraints.

preprint2016arXiv

Detection of lensing substructure using ALMA observations of the dusty galaxy SDP.81

We study the abundance of substructure in the matter density near galaxies using ALMA Science Verification observations of the strong lensing system SDP.81. We present a method to measure the abundance of subhalos around galaxies using interferometric observations of gravitational lenses. Using simulated ALMA observations, we explore the effects of various systematics, including antenna phase errors and source priors, and show how such errors may be measured or marginalized. We apply our formalism to ALMA observations of SDP.81. We find evidence for the presence of a $M=10^{8.96\pm 0.12} M_{\odot}$ subhalo near one of the images, with a significance of $6.9σ$ in a joint fit to data from bands 6 and 7; the effect of the subhalo is also detected in both bands individually. We also derive constraints on the abundance of dark matter subhalos down to $M\sim 2\times 10^7 M_{\odot}$, pushing down to the mass regime of the smallest detected satellites in the Local Group, where there are significant discrepancies between the observed population of luminous galaxies and predicted dark matter subhalos. We find hints of additional substructure, warranting further study using the full SDP.81 dataset (including, for example, the spectroscopic imaging of the lensed carbon monoxide emission). We compare the results of this search to the predictions of $Λ$CDM halos, and find that given current uncertainties in the host halo properties of SDP.81, our measurements of substructure are consistent with theoretical expectations. Observations of larger samples of gravitational lenses with ALMA should be able to improve the constraints on the abundance of galactic substructure.

preprint2015arXiv

CHITAH: Strong-gravitational-lens hunter in imaging surveys

Strong gravitationally lensed quasars provide powerful means to study galaxy evolution and cosmology. Current and upcoming imaging surveys will contain thousands of new lensed quasars, augmenting the existing sample by at least two orders of magnitude. To find such lens systems, we built a robot, CHITAH, that hunts for lensed quasars by modeling the configuration of the multiple quasar images. Specifically, given an image of an object that might be a lensed quasar, CHITAH first disentangles the light from the supposed lens galaxy and the light from the multiple quasar images based on color information. A simple rule is designed to categorize the given object as a potential four-image (quad) or two-image (double) lensed quasar system. The configuration of the identified quasar images is subsequently modeled to classify whether the object is a lensed quasar system. We test the performance of CHITAH using simulated lens systems based on the Canada-France-Hawaii Telescope Legacy Survey. For bright quads with large image separations (with Einstein radius $r_{\rm ein}>1.1"$) simulated using Gaussian point-spread functions, a high true-positive rate (TPR) of $\sim$90% and a low false-positive rate of $\sim$$3\%$ show that this is a promising approach to search for new lens systems. We obtain high TPR for lens systems with $r_{\rm ein}\gtrsim0.5"$, so the performance of CHITAH is set by the seeing. We further feed a known gravitational lens system, COSMOS 5921$+$0638, to CHITAH, and demonstrate that CHITAH is able to classify this real gravitational lens system successfully. Our newly built CHITAH is omnivorous and can hunt in any ground-based imaging surveys.

preprint2015arXiv

Correcting the z~8 Galaxy Luminosity Function for Gravitational Lensing Magnification Bias

We present a Bayesian framework to account for the magnification bias from both strong and weak gravitational lensing in estimates of high-redshift galaxy luminosity functions. We illustrate our method by estimating the $z\sim8$ UV luminosity function using a sample of 97 Y-band dropouts (Lyman break galaxies) found in the Brightest of Reionizing Galaxies (BoRG) survey and from the literature. We find the luminosity function is well described by a Schechter function with characteristic magnitude of $M^\star = -19.85^{+0.30}_{-0.35}$, faint-end slope of $α= -1.72^{+0.30}_{-0.29}$, and number density of $\log_{10} Ψ^\star [\textrm{Mpc}^{-3}] = -3.00^{+0.23}_{-0.31}$. These parameters are consistent within the uncertainties with those inferred from the same sample without accounting for the magnification bias, demonstrating that the effect is small for current surveys at $z\sim8$, and cannot account for the apparent overdensity of bright galaxies compared to a Schechter function found recently by Bowler et al. (2014a,b) and Finkelstein et al. (2014). We estimate that the probability of finding a strongly lensed $z\sim8$ source in our sample is in the range $\sim 3-15 \%$ depending on limiting magnitude. We identify one strongly-lensed candidate and three cases of intermediate lensing in BoRG (estimated magnification $μ>1.4$) in addition to the previously known candidate group-scale strong lens. Using a range of theoretical luminosity functions we conclude that magnification bias will dominate wide field surveys -- such as those planned for the Euclid and WFIRST missions -- especially at $z>10$. Magnification bias will need to be accounted for in order to derive accurate estimates of high-redshift luminosity functions in these surveys and to distinguish between galaxy formation models.

preprint2015arXiv

Discovery of two gravitationally lensed quasars in the Dark Energy Survey

We present spectroscopic confirmation of two new lensed quasars via data obtained at the 6.5m Magellan/Baade Telescope. The lens candidates have been selected from the Dark Energy Survey (DES) and WISE based on their multi-band photometry and extended morphology in DES images. Images of DES J0115-5244 show two blue point sources at either side of a red galaxy. Our long-slit data confirm that both point sources are images of the same quasar at $z_{s}=1.64.$ The Einstein Radius estimated from the DES images is $0.51$". DES J2200+0110 is in the area of overlap between DES and the Sloan Digital Sky Survey (SDSS). Two blue components are visible in the DES and SDSS images. The SDSS fiber spectrum shows a quasar component at $z_{s}=2.38$ and absorption compatible with Mg II and Fe II at $z_{l}=0.799$, which we tentatively associate with the foreground lens galaxy. The long-slit Magellan spectra show that the blue components are resolved images of the same quasar. The Einstein Radius is $0.68$" corresponding to an enclosed mass of $1.6\times10^{11}\,M_{\odot}.$ Three other candidates were observed and rejected, two being low-redshift pairs of starburst galaxies, and one being a quasar behind a blue star. These first confirmation results provide an important empirical validation of the data-mining and model-based selection that is being applied to the entire DES dataset.

preprint2015arXiv

Precision cosmology with time delay lenses: high resolution imaging requirements

Lens time delays are a powerful probe of cosmology, provided that the gravitational potential of the main deflector can be modeled with sufficient precision. Recent work has shown that this can be achieved by detailed modeling of the host galaxies of lensed quasars, which appear as "Einstein Rings" in high resolution images. We carry out a systematic exploration of the high resolution imaging required to exploit the thousands of lensed quasars that will be discovered by current and upcoming surveys with the next decade. Specifically, we simulate realistic lens systems as imaged by the Hubble Space Telescope (HST), James Webb Space Telescope (JWST), and ground based adaptive optics images taken with Keck or the Thirty Meter Telescope (TMT). We compare the performance of these pointed observations with that of images taken by the Euclid (VIS), Wide-Field Infrared Survey Telescope (WFIRST) and Large Synoptic Survey Telescope (LSST) surveys. We use as our metric the precision with which the slope $γ'$ of the total mass density profile $ρ_{tot}\propto r^{-γ'}$ for the main deflector can be measured. Ideally, we require that the statistical error on $γ'$ be less than 0.02, such that it is subdominant to other sources of random and systematic uncertainties. We find that survey data will likely have sufficient depth and resolution to meet the target only for the brighter gravitational lens systems, comparable to those discovered by the SDSS survey. For fainter systems, that will be discovered by current and future surveys, targeted follow-up will be required. However, the exposure time required with upcoming facilitites such as JWST, the Keck Next Generation Adaptive Optics System, and TMT, will only be of order a few minutes per system, thus making the follow-up of hundreds of systems a practical and efficient cosmological probe.

preprint2015arXiv

Probing the inner kpc of massive galaxies with strong gravitational lensing

We examine the prospects of detecting demagnified images of gravitational lenses in observations of strongly lensed mm-wave molecular emission lines with ALMA. We model the lensing galaxies as a superposition of a dark matter component, a stellar component, and a central supermassive black hole and assess the detectability of the central images for a range of relevant parameters (e.g., stellar core, black hole mass, and source size). We find that over a large range of plausible parameters, future deep observations of lensed molecular lines with ALMA should enable detection of the central images at $\gtrsim 3σ$ significance. We use a Fisher analysis to examine the constraints that could be placed on these parameters in various scenarios and find that for large stellar cores, both the core size and the mass of the central SMBHs can be accurately measured. We also study the prospects for detecting binary SMBHs with such observations and find that only under rare conditions and with very long integrations ($\sim$40-hr) the masses of both SMBHs may be measured using the distortions of central images.

preprint2015arXiv

Space Warps II. New Gravitational Lens Candidates from the CFHTLS Discovered through Citizen Science

We report the discovery of 29 promising (and 59 total) new lens candidates from the CFHT Legacy Survey (CFHTLS) based on about 11 million classifications performed by citizen scientists as part of the first Space Warps lens search. The goal of the blind lens search was to identify lens candidates missed by robots (the RingFinder on galaxy scales and ArcFinder on group/cluster scales) which had been previously used to mine the CFHTLS for lenses. We compare some properties of the samples detected by these algorithms to the Space Warps sample and find them to be broadly similar. The image separation distribution calculated from the Space Warps sample shows that previous constraints on the average density profile of lens galaxies are robust. SpaceWarps recovers about 65% of known lenses, while the new candidates show a richer variety compared to those found by the two robots. This detection rate could be increased to 80% by only using classifications performed by expert volunteers (albeit at the cost of a lower purity), indicating that the training and performance calibration of the citizen scientists is very important for the success of Space Warps. In this work we present the SIMCT pipeline, used for generating in situ a sample of realistic simulated lensed images. This training sample, along with the false positives identified during the search, has a legacy value for testing future lens finding algorithms. We make the pipeline and the training set publicly available.

preprint2015arXiv

Space Warps: I. Crowd-sourcing the Discovery of Gravitational Lenses

We describe Space Warps, a novel gravitational lens discovery service that yields samples of high purity and completeness through crowd-sourced visual inspection. Carefully produced colour composite images are displayed to volunteers via a web- based classification interface, which records their estimates of the positions of candidate lensed features. Images of simulated lenses, as well as real images which lack lenses, are inserted into the image stream at random intervals; this training set is used to give the volunteers instantaneous feedback on their performance, as well as to calibrate a model of the system that provides dynamical updates to the probability that a classified image contains a lens. Low probability systems are retired from the site periodically, concentrating the sample towards a set of lens candidates. Having divided 160 square degrees of Canada-France-Hawaii Telescope Legacy Survey (CFHTLS) imaging into some 430,000 overlapping 82 by 82 arcsecond tiles and displaying them on the site, we were joined by around 37,000 volunteers who contributed 11 million image classifications over the course of 8 months. This Stage 1 search reduced the sample to 3381 images containing candidates; these were then refined in Stage 2 to yield a sample that we expect to be over 90% complete and 30% pure, based on our analysis of the volunteers performance on training images. We comment on the scalability of the SpaceWarps system to the wide field survey era, based on our projection that searches of 10$^5$ images could be performed by a crowd of 10$^5$ volunteers in 6 days.

preprint2015arXiv

The Sloan Lens ACS Survey. XII. Extending Strong Lensing to Lower Masses

We present observational results from a new Hubble Space Telescope (HST) Snapshot program to extend the methods of the Sloan Lens ACS (SLACS) Survey to lower lens-galaxy masses. We discover 40 new galaxy-scale strong lenses, which we supplement with 58 previously discovered SLACS lenses. In addition, we determine the posterior PDFs of the Einstein radius for 33 galaxies (18 new and 15 from legacy SLACS data) based on single lensed images. We find a less-than-unity slope of $0.64\pm0.06$ for the $\log_{10} σ_*$-$\log_{10} σ_{\rm SIE}$ relation, which corresponds to a 6-$σ$ evidence that the total mass-density profile of early-type galaxies varies systematically in the sense of being shallower at higher lens-galaxy velocity dispersions. The trend is only significant when single-image systems are considered, highlighting the importance of including both "lenses" and "non-lenses" for an unbiased treatment of the lens population when extending to lower mass ranges. By scaling simple stellar population models to the HST I-band data, we identify a strong trend of increasing dark-matter fraction at higher velocity dispersions, which can be alternatively interpreted as a trend in the stellar initial mass function (IMF) normalization. Consistent with previous findings and the suggestion of a non-universal IMF, we find that a Salpeter IMF is ruled out for galaxies with velocity dispersion less than $180$ km/s. Considered together, our mass-profile and dark-matter-fraction trends with increasing galaxy mass could both be explained by an increasing relative contribution on kiloparsec scales from a dark-matter halo with a spatial profile more extended than that of the stellar component.

preprint2014arXiv

Data Mining for Gravitationally Lensed Quasars

Gravitationally lensed (GL) quasars are brighter than their unlensed counterparts and produce images with distinctive morphological signatures. Past searches and target selection algorithms, in particular the Sloan Quasar Lens Search (SQLS), have relied on basic morphological criteria, which were applied to samples of bright, spectroscopically confirmed quasars. The SQLS techniques are not sufficient for searching into new surveys (e.g. DES, PS1, LSST), because spectroscopic information is not readily available and the large data volume requires higher purity in target/candidate selection. We carry out a systematic exploration of machine learning techniques and demonstrate that a two step strategy can be highly effective. In the first step we use catalog-level information ($griz$+WISE magnitudes, second moments) to preselect targets, using artificial neural networks. The accepted targets are then inspected with pixel-by-pixel pattern recognition algorithms (Gradient-Boosted Trees), to form a final set of candidates. The results from this procedure can be used to further refine the simpler SQLS algorithms, with a twofold (or threefold) gain in purity and the same (or $80\%$) completeness at target-selection stage, or a purity of $70\%$ and a completeness of $60\%$ after the candidate-selection step. Simpler photometric searches in $griz$+WISE based on colour cuts would provide samples with $7\%$ purity or less. Our technique is extremely fast, as a list of candidates can be obtained from a stage III experiment (e.g. DES catalog/database) in {a few} CPU hours. The techniques are easily extendable to Stage IV experiments like LSST with the addition of time domain information.

preprint2014arXiv

Hierarchical probabilistic inference of cosmic shear

Point estimators for the shearing of galaxy images induced by gravitational lensing involve a complex inverse problem in the presence of noise, pixelization, and model uncertainties. We present a probabilistic forward modeling approach to gravitational lensing inference that has the potential to mitigate the biased inferences in most common point estimators and is practical for upcoming lensing surveys. The first part of our statistical framework requires specification of a likelihood function for the pixel data in an imaging survey given parameterized models for the galaxies in the images. We derive the lensing shear posterior by marginalizing over all intrinsic galaxy properties that contribute to the pixel data (i.e., not limited to galaxy ellipticities) and learn the distributions for the intrinsic galaxy properties via hierarchical inference with a suitably flexible conditional probabilitiy distribution specification. We use importance sampling to separate the modeling of small imaging areas from the global shear inference, thereby rendering our algorithm computationally tractable for large surveys. With simple numerical examples we demonstrate the improvements in accuracy from our importance sampling approach, as well as the significance of the conditional distribution specification for the intrinsic galaxy properties when the data are generated from an unknown number of distinct galaxy populations with different morphological characteristics.

preprint2014arXiv

Ideas for Citizen Science in Astronomy

We review the relatively new, internet-enabled, and rapidly-evolving field of citizen science, focusing on research projects in stellar, extragalactic and solar system astronomy that have benefited from the participation of members of the public, often in large numbers. We find these volunteers making contributions to astronomy in a variety of ways: making and analyzing new observations, visually classifying features in images and light curves, exploring models constrained by astronomical datasets, and initiating new scientific enquiries. The most productive citizen astronomy projects involve close collaboration between the professionals and amateurs involved, and occupy scientific niches not easily filled by great observatories or machine learning methods: citizen astronomers are most strongly motivated by being of service to science. In the coming years we expect participation and productivity in citizen astronomy to increase, as survey datasets get larger and citizen science platforms become more efficient. Opportunities include engaging the public in ever more advanced analyses, and facilitating citizen-led enquiry by designing professional user interfaces and analysis tools with citizens in mind.

preprint2014arXiv

RingFinder: automated detection of galaxy-scale gravitational lenses in ground-based multi-filter imaging data

We present RingFinder, a tool for finding galaxy-scale strong gravitational lenses in multiband imaging data. By construction, the method is sensitive to configurations involving a massive foreground early-type galaxy and a faint, background, blue source. RingFinder detects the presence of blue residuals embedded in an otherwise smooth red light distribution by difference imaging in two bands. The method is automated for efficient application to current and future surveys, having originally been designed for the 150-deg2 Canada France Hawaii Telescope Legacy Survey (CFHTLS). We describe each of the steps of RingFinder. We then carry out extensive simulations to assess completeness and purity. For sources with magnification mu>4, RingFinder reaches 42% (resp. 25%) completeness and 29% (resp. 86%) purity before (resp. after) visual inspection. The completeness of RingFinder is substantially improved in the particular range of Einstein radii 0.8 < REin < 2. and lensed images brighter than g = 22.5, where it can be as high as 70%. RingFinder does not introduce any significant bias in the source or deflector population. We conclude by presenting the final catalog of RingFinder CFHTLS galaxy-scale strong lens candidates. Additional information obtained with Hubble Space Telescope and Keck Adaptive Optics high resolution imaging, and with Keck and Very Large Telescope spectroscopy, is used to assess the validity of our classification, and measure the redshift of the foreground and the background objects. From an initial sample of 640,000 early type galaxies, RingFinder returns 2500 candidates, which we further reduce by visual inspection to 330 candidates. We confirm 33 new gravitational lenses from the main sample of candidates, plus an additional 16 systems taken from earlier versions of RingFinder. First applications are presented in the SL2S galaxy-scale Lens Sample paper series.

preprint2014arXiv

The SL2S Galaxy-scale Lens Sample. V. Dark Matter Halos and Stellar IMF of Massive Early-type Galaxies out to Redshift 0.8

We investigate the cosmic evolution of the internal structure of massive early-type galaxies over half of the age of the Universe. We perform a joint lensing and stellar dynamics analysis of a sample of 81 strong lenses from the SL2S and SLACS surveys and combine the results with a hierarchical Bayesian inference method to measure the distribution of dark matter mass and stellar IMF across the population of massive early-type galaxies. Lensing selection effects are taken into account. We find that the dark matter mass projected within the inner 5 kpc increases for increasing redshift, decreases for increasing stellar mass density, but is roughly constant along the evolutionary tracks of early-type galaxies. The average dark matter slope is consistent with that of an NFW profile, but is not well constrained. The stellar IMF normalization is close to a Salpeter IMF at $\log{M_*} = 11.5$ and scales strongly with increasing stellar mass. No dependence of the IMF on redshift or stellar mass density is detected. The anti-correlation between dark matter mass and stellar mass density supports the idea of mergers being more frequent in more massive dark matter halos.

preprint2013arXiv

Dark Matter Astrometry: Accuracy of sub-halo positions for the measurement of self-interaction cross sections

Direct evidence for the existence of dark matter and measurements of its interaction cross-section have been provided by the physical offset between dark matter and intra- cluster gas in merging systems like the Bullet Cluster. Although a smaller signal, this effect is more abundant in minor mergers where infalling substructure dark matter and gas are segregated. In such low-mass systems the gravitational lensing signal comes primarily from weak lensing. A fundamental step in determining such an offset in sub- structure is the ability to accurately measure the positions of dark matter sub-peaks. Using simulated Hubble Space Telescope observations, we make a first assessment of the precision and accuracy with which we can measure infalling groups using weak gravitational lensing. We demonstrate that using an existing and well-used mass re- construction algorithm can measure the positions of 1.5x1013 M substructures that have parent halos ten times more massive with a bias of less than 0.3 . In this regime, our analysis suggests the precision is sufficient to detect (at 3 σ statistical significance) the expected mean offset between dark matter and baryonic gas in infalling groups from a sample of 50 massive clusters.

preprint2013arXiv

Improving the precision of time-delay cosmography with observations of galaxies along the line of sight

In order to use strong gravitational lens time delays to measure precise and accurate cosmological parameters the effects of mass along the line of sight must be taken into account. We present a method to achieve this by constraining the probability distribution function of the effective line of sight convergence k_ext. The method is based on matching the observed overdensity in the weighted number of galaxies to that found in mock catalogs with k_ext obtained by ray-tracing through structure formation simulations. We explore weighting schemes based on projected distance, mass, luminosity, and redshift. This additional information reduces the uncertainty of k_ext from sigma_k $0.06 to ~0.04 for very overdense lines of sight like that of the system B1608+656. For more common lines of sight, sigma_k is reduced to ~<0.03, corresponding to an uncertainty of ~<3% on distance. This uncertainty has comparable effects on cosmological parameters to that arising from the mass model of the deflector and its immediate environment. Photometric redshifts based on g, r, i and K photometries are sufficient to constrain k_ext almost as well as with spectroscopic redshifts. As an illustration, we apply our method to the system B1608+656. Our most reliable k_ext estimator gives sigma_k=0.047 down from 0.065 using only galaxy counts. Although deeper multi-band observations of the field of B1608+656 are necessary to obtain a more precise estimate, we conclude that griK photometry, in addition to spectroscopy to characterize the immediate environment, is an effective way to increase the precision of time-delay cosmography.

preprint2013arXiv

Reconstructing the Lensing Mass in the Universe from Photometric Catalogue Data

High precision cosmological distance measurements towards individual objects such as time delay gravitational lenses or type Ia supernovae are affected by weak lensing perturbations by galaxies and groups along the line of sight. In time delay gravitational lenses, "external convergence," kappa, can dominate the uncertainty in the inferred distances and hence cosmological parameters. In this paper we attempt to reconstruct kappa, due to line of sight structure, using a simple halo model. We use mock catalogues from the Millennium Simulation, and calibrate and compare our reconstructed P(kappa) to ray-traced kappa "truth" values; taking into account realistic observational uncertainties. We find that the reconstruction of kappa provides an improvement in precision of ~50% over galaxy number counts. We find that the lowest-kappa lines of sight have the best constrained P(kappa). In anticipation of large future samples of lenses, we find that selecting the third of the systems with the highest precision kappa estimates gives a sample of unbiased time delay distance measurements with just ~1% uncertainty due to line of sight external convergence effects. Photometric data are sufficient to pre-select the best-constrained lines of sight, and can be done before investment in light-curve monitoring. Conversely, we show that selecting lines of sight with high external shear could, with the reconstruction model presented, induce biases of up to 1% in time delay distance. We find that a major potential source of systematic error is uncertainty in the high mass end of the stellar mass-halo mass relation; this could introduce ~2% biases on the time-delay distance if completely ignored. We suggest areas for the improvement of this general analysis framework (including more sophisticated treatment of high mass structures) that should allow yet more accurate cosmological inferences to be made.

preprint2013arXiv

The SL2S Galaxy-scale Lens Sample. III. Lens Models, Surface Photometry and Stellar Masses for the final sample

We present Hubble Space Telescope (HST) imaging data and CFHT Near IR ground-based images for the final sample of 56 candidate galaxy-scale lenses uncovered in the CFHT Legacy Survey as part of the Strong Lensing in the Legacy Survey (SL2S) project. The new images are used to perform lens modeling, measure surface photometry, and estimate stellar masses of the deflector early-type galaxies. Lens modeling is performed on the HST images (or CFHT when HST is not available) by fitting the spatially extended light distribution of the lensed features assuming a singular isothermal ellipsoid mass profile and by reconstructing the intrinsic source light distribution on a pixelized grid. Based on the analysis of systematic uncertainties and comparison with inference based on different methods we estimate that our Einstein Radii are accurate to \sim3%. HST imaging provides a much higher success rate in confirming gravitational lenses and measuring their Einstein radii than CFHT imaging does. Lens modeling with ground-based images however, when successful, yields Einstein radius measurements that are competitive with spaced-based images. Information from the lens models is used together with spectroscopic information from the companion paper IV to classify the systems, resulting in a final sample of 39 confirmed (grade-A) lenses and 17 promising candidates. The redshifts of the main deflector span a range 0.3<zd< 0.8, providing an excellent sample for the study of the cosmic evolution of the mass distribution of early-type galaxies over the second half of the history of the Universe.

preprint2013arXiv

The SL2S Galaxy-scale Lens Sample. IV. The dependence of the total mass density profile of early-type galaxies on redshift, stellar mass, and size

We present optical and near infrared spectroscopy obtained at Keck, VLT, and Gemini for a sample of 36 secure strong gravitational lens systems and 17 candidates identified as part of the SL2S survey. The deflectors are massive early-type galaxies in the redshift range z_d=0.2-0.8, while the lensed sources are at z_s=1-3.5. We combine this data with photometric and lensing measurements presented in the companion paper III and with lenses from the SLACS and LSD surveys to investigate the cosmic evolution of the internal structure of massive early-type galaxies over half the age of the universe. We study the dependence of the slope of the total mass density profile γ' (ρ(r)\propto r^{-γ'}) on stellar mass, size, and redshift. We find that two parameters are sufficent to determine γ' with less than 6% residual scatter. At fixed redshift, γ' depends solely on the surface stellar mass density \partial γ'/ \partial Σ_*=0.38\pm 0.07, i.e. galaxies with denser stars also have steeper slopes. At fixed M_* and R_{eff}, γ' depends on redshift, in the sense that galaxies at a lower redshift have steeper slopes (\partial γ' / \partial z = -0.31\pm 0.10). However, the mean redshift evolution of γ' for an individual galaxy is consistent with zero dγ'/dz=-0.10\pm0.12. This result is obtained by combining our measured dependencies of γ' on z,M_*,R_{eff} with the evolution of the R_{eff}-M_* taken from the literature, and is broadly consistent with current models of the formation and evolution of massive early-type galaxies. Detailed quantitative comparisons of our results with theory will provide qualitatively new information on the detailed physical processes at work.

preprint2013arXiv

The SWELLS Survey. VI. hierarchical inference of the initial mass functions of bulges and discs

The long-standing assumption that the stellar initial mass function (IMF) is universal has recently been challenged by a number of observations. Several studies have shown that a "heavy" IMF (e.g., with a Salpeter-like abundance of low mass stars and thus normalisation) is preferred for massive early-type galaxies, while this IMF is inconsistent with the properties of less massive, later-type galaxies. These discoveries motivate the hypothesis that the IMF may vary (possibly very slightly) across galaxies and across components of individual galaxies (e.g. bulges vs discs). In this paper we use a sample of 19 late-type strong gravitational lenses from the SWELLS survey to investigate the IMFs of the bulges and discs in late-type galaxies. We perform a joint analysis of the galaxies' total masses (constrained by strong gravitational lensing) and stellar masses (constrained by optical and near-infrared colours in the context of a stellar population synthesis [SPS] model, up to an IMF normalisation parameter). Using minimal assumptions apart from the physical constraint that the total stellar mass within any aperture must be less than the total mass within the aperture, we find that the bulges of the galaxies cannot have IMFs heavier (i.e. implying high mass per unit luminosity) than Salpeter, while the disc IMFs are not well constrained by this data set. We also discuss the necessity for hierarchical modelling when combining incomplete information about multiple astronomical objects. This modelling approach allows us to place upper limits on the size of any departures from universality. More data, including spatially resolved kinematics (as in paper V) and stellar population diagnostics over a range of bulge and disc masses, are needed to robustly quantify how the IMF varies within galaxies.

preprint2012arXiv

Constraining the dark energy equation of state with double source plane strong lenses

We investigate the possibility of constraining the dark energy equation of state by measuring the ratio of Einstein radii in a strong gravitational lens system with two source planes. This quantity is independent of the Hubble parameter and directly measures the growth of angular diameter distances as a function of redshift. We investigate the prospects for a single double source plane system and for a forecast population of systems discovered by re-observing a population of single source lenses already known from a photometrically selected catalogue such as CASSOWARY or from a spectroscopically selected catalogue such as SLACS. We find that constraints comparable to current data-sets (15% uncertainty on the dark equation of state at 68%CL) are possible with a handful of double source plane systems. We also find that the method's degeneracy between Omega_M and w is almost orthogonal to that of CMB and BAO measurements, making this method highly complimentary to current probes.

preprint2012arXiv

Evidence for dark matter contraction and a Salpeter IMF in a massive early-type galaxy

Stars and dark matter account for most of the mass of early-type galaxies, but uncertainties in the stellar population and the dark matter profile make it challenging to distinguish between the two components. Nevertheless, precise observations of stellar and dark matter are extremely valuable for testing the many models of structure formation and evolution. We present a measurement of the stellar mass and inner slope of the dark matter halo of a massive early-type galaxy at $z=0.222$. The galaxy is the foreground deflector of the double Einstein ring gravitational lens system SDSSJ0946+1006, also known as the "Jackpot". By combining the tools of lensing and dynamics, we first constrain the mean slope of the total mass density profile ($ρ_{\rm{tot}}\propto r^{-γ'}$) within the radius of the outer ring to be $γ' = 1.98\pm0.02\pm0.01$. Then we obtain a bulge-halo decomposition, assuming a power-law form for the dark matter halo. Our analysis yields $γ_{\rm{DM}} = 1.7\pm0.2$ for the inner slope of the dark matter profile, in agreement with theoretical findings on the distribution of dark matter in ellipticals, and a stellar mass from lensing and dynamics $M_*^{\rm{LD}} = 5.5_{-1.3}^{+0.4}\times10^{11}M_\Sun$. By comparing this measurement with stellar masses inferred from stellar population synthesis fitting we find that a Salpeter IMF provides a good description of the stellar population of the lens while the probability of the IMF being heavier than Chabrier is 95%. Our data suggest that growth by accretion of small systems from a compact red nugget is a plausible formation scenario for this object.

preprint2012arXiv

Nomads of the Galaxy

We estimate that there may be up to ~10^5 compact objects in the mass range 10^{-8} -10^{-2} solar mass per main sequence star that are unbound to a host star in the Galaxy. We refer to these objects as nomads; in the literature a subset of these are sometimes called free-floating or rogue planets. Our estimate for the number of Galactic nomads is consistent with a smooth extrapolation of the mass function of unbound objects above the Jupiter-mass scale, the stellar mass density limit, and the metallicity of the interstellar medium. We analyze the prospects for detecting nomads via Galactic microlensing. The Wide-Field Infrared Survey Telescope (WFIRST) will measure the number of nomads per main sequence star greater than the mass of Jupiter to ~ 13%, and the corresponding number greater than the mass of Mars to ~25%. All-sky surveys such as GAIA and LSST can identify nomads greater than about the mass of Jupiter. We suggest a dedicated drift scanning telescope that covers approximately 100 square degrees in the Southern hemisphere could identify nomads as small as 10^{-8} solar mass via microlensing of bright stars with characteristic lightcurve timescales of a few seconds.

preprint2012arXiv

The SL2S Galaxy-scale Gravitational Lens Sample. I. The alignment of mass and light in massive early-type galaxies at z=0.2-0.9

We study the relative alignment of mass and light in a sample of 16 massive early-type galaxies at z=0.2-0.9 that act as strong gravitational lenses. The sample was identified from deep multi-band images obtained as part of the Canada France Hawaii Telescope Legacy Survey as part of the Strong Lensing Legacy Survey (SL2S). Higher resolution follow-up imaging is available for a subset of 10 systems. We construct gravitational lens models and infer total enclosed mass, elongation, and position angle of the mass distribution. By comparison with the observed distribution of light we infer that there is a substantial amount of external shear with mean value $<γ_{\rm ext}>= 0.12 \pm 0.05$, arising most likely from the environment of the SL2S lenses. In a companion paper, we combine these measurements with follow-up Keck spectroscopy to study the evolution of the stellar and dark matter content of early-type galaxies as a function of cosmic time.

preprint2012arXiv

The SWELLS survey. III. Disfavouring "heavy" initial mass functions for spiral lens galaxies

We present gravitational lens models for 20 strong gravitational lens systems observed as part of the Sloan WFC Edge-on Late-type Lens Survey (SWELLS) project. Fifteen of the lenses are taken from paper I while five are newly discovered systems. The systems are galaxy-galaxy lenses where the foreground deflector has an inclined disc, with a wide range of morphological types, from late-type spiral to lenticular. For each system, we compare the total mass inside the critical curve inferred from gravitational lens modelling to the stellar mass inferred from stellar population synthesis (SPS) models, computing the stellar mass fraction f* = M(SPS)/M(lens). We find that, for the lower mass SWELLS systems, adoption of a Salpeter stellar initial mass function (IMF) leads to estimates of f* that exceed 1. This is unphysical, and provides strong evidence against the Salpeter IMF being valid for these systems. Taking the lower mass end of the SWELLS sample sigma(SIE) < 230 km/s, we find that the IMF is lighter (in terms of stellar mass-to-light ratio) than Salpeter with 98% probability, and consistent with the Chabrier IMF and IMFs between the two. This result is consistent with previous studies of spiral galaxies based on independent techniques. In combination with recent studies of massive early-type galaxies that have favoured a heavier Salpeter-like IMF, this result strengthens the evidence against a universal stellar IMF.

preprint2012arXiv

The SWELLS survey. IV. Precision measurements of the stellar and dark matter distributions in a spiral lens galaxy

We construct a fully self-consistent mass model for the lens galaxy J2141 at z=0.14, and use it to improve on previous studies by modelling its gravitational lensing effect, gas rotation curve and stellar kinematics simultaneously. We adopt a very flexible axisymmetric mass model constituted by a generalized NFW dark matter halo and a stellar mass distribution obtained by deprojecting the MGE fit to the high-resolution K'-band LGSAO imaging data of the galaxy, with the (spatially constant) M/L ratio as a free parameter. We model the stellar kinematics by solving the anisotropic Jeans equations. We find that the inner logarithmic slope of the dark halo is weakly constrained (gamma = 0.82^{+0.65}_{-0.54}), and consistent with an unmodified NFW profile. We infer the galaxy to have (i) a dark matter fraction within 2.2 disk radii of 0.28^{+0.15}_{-0.10}, independent of the galaxy stellar population, implying a maximal disk for J2141; (ii) an apparently uncontracted dark matter halo, with concentration c_{-2} = 7.7_{-2.5}^{+4.2} and virial velocity v_{vir} = 242_{-39}^{+44} km/s, consistent with LCDM predictions; (iii) a slightly oblate halo (q_h = 0.75^{+0.27}_{-0.16}), consistent with predictions from baryon-affected models. Comparing the stellar mass inferred from the combined analysis (log_{10} Mstar/Msun = 11.12_{-0.09}^{+0.05}) with that inferred from SPS modelling of the galaxies colours, and accounting for a cold gas fraction of 20+/-10%, we determine a preference for a Chabrier IMF over Salpeter IMF by a Bayes factor of 5.7 (substantial evidence). We infer a value beta_{z} = 1 - sigma^2_{z}/sigma^2_{R} = 0.43_{-0.11}^{+0.08} for the orbital anisotropy parameter in the meridional plane, in agreement with most studies of local disk galaxies, and ruling out at 99% CL that the dynamics of this system can be described by a two-integral distribution function. [Abridged]

preprint2011arXiv

The Mass Distribution and Assembly of the Milky Way from the Properties of the Magellanic Clouds

We present a new measurement of the mass of the Milky Way (MW) based on observed properties of its largest satellite galaxies, the Magellanic Clouds (MCs), and an assumed prior of a ΛCDM universe. The large, high-resolution Bolshoi cosmological simulation of this universe provides a means to statistically sample the dynamical properties of bright satellite galaxies in a large population of dark matter halos. The observed properties of the MCs, including their circular velocity, distance from the center of the MW, and velocity within the MW halo, are used to evaluate the likelihood that a given halo would have each or all of these properties; the posterior PDF for any property of the MW system can thus be constructed. This method provides a constraint on the MW virial mass, 1.2 +0.7 -0.3(stat.) +0.3 -0.4 (sys.) x 10^12 M\odot (68% confidence), which is consistent with recent determinations that involve very different assumptions. In addition, we calculate the posterior PDF for the density profile of the MW and its satellite accretion history. Although typical satellites of 10^12 M\odot halos are accreted over a wide range of epochs over the last 10 Gyr, we find a \sim72% probability that the Magellanic Clouds were accreted within the last Gyr, and a 50% probability that they were accreted together.

preprint2011arXiv

The Sloan Lens ACS Survey. XI. Beyond Hubble resolution: size, luminosity and stellar mass of compact lensed galaxies at intermediate redshift

We exploit the strong lensing effect to explore the properties of intrinsically faint and compact galaxies at intermediate redshift, at the highest possible resolution at optical wavelengths. Our sample consists of 46 strongly-lensed emission line galaxies discovered by the Sloan Lens ACS (SLACS) Survey. The galaxies have been imaged at high resolution with HST in three bands (V_HST, I_814 and H_160), allowing us to infer their size, luminosity, and stellar mass using stellar population synthesis models. Lens modeling is performed using a new fast and robust code, klens, which we test extensively on real and synthetic non-lensed galaxies, and also on simulated galaxies multiply-imaged by SLACS- like galaxy-scale lenses. Our tests show that our measurements of galaxy size, flux, and Sersic index are robust and accurate, even for objects intrinsically smaller than the HST point spread function. The median magnification is 8.8, with a long tail that extends to magnifications above 40. Modeling the SLACS sources reveals a population of galaxies with colors and Sersic indices (median n ~ 1) consistent with the objects detected in the field with HST in the GEMS survey, but that are (typically) ~ 2 magnitudes fainter and ~ 5 times smaller in apparent size. The closest analog are ultracompact emission line galaxies identified by HST grism surveys. The lowest mass galaxies in our sample are comparable to the brightest Milky Way satellites in stellar mass (10^7 solar masses) and have well-determined half light radii of 0."05 (~0.3 kpc).

preprint2010arXiv

AMUSE-Virgo II. Down-sizing in black hole accretion

(Abridged) We complete the census of nuclear X-ray activity in 100 early type Virgo galaxies observed by the Chandra X-ray Telescope as part of the AMUSE-Virgo survey, down to a (3sigma) limiting luminosity of 3.7E+38 erg/s over 0.5-7 keV. The stellar mass distribution of the targeted sample, which is mostly composed of formally `inactive' galaxies, peaks below 1E+10 M_Sun, a regime where the very existence of nuclear super-massive black holes (SMBHs) is debated. Out of 100 objects, 32 show a nuclear X-ray source, including 6 hybrid nuclei which also host a massive nuclear cluster as visible from archival HST images. After carefully accounting for contamination from nuclear low-mass X-ray binaries based on the shape and normalization of their X-ray luminosity function, we conclude that between 24-34% of the galaxies in our sample host a X-ray active SMBH (at the 95% C.L.). This sets a firm lower limit to the black hole occupation fraction in nearby bulges within a cluster environment. At face value, the active fraction -down to our luminosity limit- is found to increase with host stellar mass. However, taking into account selection effects, we find that the average Eddington-scaled X-ray luminosity scales with black hole mass as M_BH^(-0.62^{+0.13}_{-0.12}), with an intrinsic scatter of 0.46^({+0.08}_{-0.06}) dex. This finding can be interpreted as observational evidence for `down-sizing' of black hole accretion in local early types, that is, low mass black holes shine relatively closer to their Eddington limit than higher mass objects. As a consequence, the fraction of active galaxies, defined as those above a fixed X-ray Eddington ratio, decreases with increasing black hole mass.

preprint2010arXiv

Cosmic Evolution of Virial and Stellar Mass in Massive Early-Type Galaxies

We measure the average mass properties of a sample of 41 strong gravitational lenses at moderate redshift (z ~ 0.4 - 0.9), and present the lens redshift for 6 of these galaxies for the first time. Using the techniques of strong and weak gravitational lensing on archival data obtained from the Hubble Space Telescope, we determine that the average mass overdensity profile of the lenses can be fit with a power-law profile (Delta_Sigma prop. to R^{-0.86 +/- 0.16}) that is within 1-sigma of an isothermal profile (Delta_Sigma prop. to R^{-1}) with velocity dispersion sigma_v = 260 +/- 20 km/s. Additionally, we use a two-component de Vaucouleurs+NFW model to disentangle the total mass profile into separate luminous and dark matter components, and determine the relative fraction of each component. We measure the average rest frame V-band stellar mass-to-light ratio (Upsilon_V = 4.0 +/- 0.6 h M_sol/L_sol) and virial mass-to-light ratio (tau_V = 300 +/- 90 h M_sol/L_sol) for our sample, resulting in a virial-to-stellar mass ratio of M_vir/M_* = 75 +/- 25. Finally, we compare our results to a previous study using low redshift lenses, to understand how galaxy mass profiles evolve over time. We investigate the evolution of M_vir/M_*(z) = alpha(1+z)^{beta}, and find best fit parameters of alpha = 51 +/- 36 and beta = 0.9 +/- 1.8, constraining the growth of virial to stellar mass ratio over the last ~7 Gigayears. We note that, by using a sample of strong lenses, we are able to constrain the growth of M_vir/M_*(z) without making any assumptions about the IMF of the stellar population.

preprint2010arXiv

Gravitationally lensed quasars and supernovae in future wide-field optical imaging surveys

Cadenced optical imaging surveys in the next decade will be capable of detecting time-varying galaxy-scale strong gravitational lenses in large numbers, increasing the size of the statistically well-defined samples of multiply-imaged quasars by two orders of magnitude, and discovering the first strongly-lensed supernovae. We carry out a detailed calculation of the likely yields of several planned surveys, using realistic distributions for the lens and source properties and taking magnification bias and image configuration detectability into account. We find that upcoming wide-field synoptic surveys should detect several thousand lensed quasars. In particular, the LSST should find 8000 lensed quasars, 3000 of which will have well-measured time delays, and also ~130 lensed supernovae, which is compared with ~15 lensed supernovae predicted to be found by the JDEM. We predict the quad fraction to be ~15% for the lensed quasars and ~30% for the lensed supernovae. Generating a mock catalogue of around 1500 well-observed double-image lenses, we compute the available precision on the Hubble constant and the dark energy equation parameters for the time delay distance experiment (assuming priors from Planck): the predicted marginalised 68% confidence intervals are σ(w_0)=0.15, σ(w_a)=0.41, and σ(h)=0.017. While this is encouraging in the sense that these uncertainties are only 50% larger than those predicted for a space-based type-Ia supernova sample, we show how the dark energy figure of merit degrades with decreasing knowledge of the the lens mass distribution. (Abridged)

preprint2010arXiv

Selecting Quasars by their Intrinsic Variability

We present a new and simple technique for selecting extensive, complete and pure quasar samples, based on their intrinsic variability. We parametrize the single-band variability by a power-law model for the light-curve structure function, with amplitude A and power-law index gamma. We show that quasars can be efficiently separated from other non-variable and variable sources by the location of the individual sources in the A-gamma plane. We use ~60 epochs of imaging data, taken over ~5 years, from the SDSS stripe 82 (S82) survey, where extensive spectroscopy provides a reference sample of quasars, to demonstrate the power of variability as a quasar classifier in multi-epoch surveys. For UV-excess selected objects, variability performs just as well as the standard SDSS color selection, identifying quasars with a completeness of 90% and a purity of 95%. In the redshift range 2.5<z<3, where color selection is known to be problematic, variability can select quasars with a completeness of 90% and a purity of 96%. This is a factor of 5-10 times more pure than existing color-selection of quasars in this redshift range. Selecting objects from a broad griz color box without u-band information, variability selection in S82 can afford completeness and purity of 92%, despite a factor of 30 more contaminants than quasars in the color-selected feeder sample. This confirms that the fraction of quasars hidden in the 'stellar locus' of color-space is small. To test variability selection in the context of Pan-STARRS 1 (PS1) we created mock PS1 data by down-sampling the S82 data to just 6 epochs over 3 years. Even with this much sparser time sampling, variability is an encouragingly efficient classifier. For instance, a 92% pure and 44% complete quasar candidate sample is attainable from the above $griz$-selected catalog.

preprint2010arXiv

The SL2S Galaxy-scale Lens Sample. II. Cosmic evolution of dark and luminous mass in early-type galaxies

We present a joint gravitational lensing and stellar-dynamical analysis of 11 early-type galaxies (median deflector redshift $\zd=0.5$) from Strong Lenses in the Legacy Survey (SL2S). Using newly measured redshifts and stellar velocity dispersions from Keck spectroscopy with lens models from Paper I, we derive the total mass density slope inside the Einstein radius for each of the 11 lenses. The average total density slope is found to be $<γ'>= 2.16^{+0.09}_{-0.09}$ ($ρ_{\rm tot}\propto r^{-γ'}$), with an intrinsic scatter of $0.25^{+0.10}_{-0.07}$. We also determine the dark matter fraction for each lens within half the effective radius, and find the average projected dark matter mass fraction to be $0.42^{+0.08}_{-0.08}$ with a scatter of $0.20^{+0.09}_{-0.07}$ for a Salpeter IMF. By combining the SL2S results with those from the Sloan Lens ACS Survey (median $\zd=0.2$) and the Lenses Structure and Dynamics survey (median $\zd=0.8$), we investigate cosmic evolution of $γ'$ and find a mild trend $\partial<γ'>/\partial\zd = -0.25^{+0.10}_{-0.12}$. This suggests that the total density profile of massive galaxies has become slightly steeper over cosmic time. If this result is confirmed by larger samples, it would indicate that dissipative processes played some role in the growth of massive galaxies since $z\sim1$.

preprint2008arXiv

Automated detection of galaxy-scale gravitational lenses in high resolution imaging data

Lens modeling is the key to successful and meaningful automated strong galaxy-scale gravitational lens detection. We have implemented a lens-modeling "robot" that treats every bright red galaxy (BRG) in a large imaging survey as a potential gravitational lens system. Using a simple model optimized for "typical" galaxy-scale lenses, we generate four assessments of model quality that are used in an automated classification. The robot infers the lens classification parameter H that a human would have assigned; the inference is performed using a probability distribution generated from a human-classified training set, including realistic simulated lenses and known false positives drawn from the HST/EGS survey. We compute the expected purity, completeness and rejection rate, and find that these can be optimized for a particular application by changing the prior probability distribution for H, equivalent to defining the robot's "character." Adopting a realistic prior based on the known abundance of lenses, we find that a lens sample may be generated that is ~100% pure, but only ~20% complete. This shortfall is due primarily to the over-simplicity of the lens model. With a more optimistic robot, ~90% completeness can be achieved while rejecting ~90% of the candidate objects. The remaining candidates must be classified by human inspectors. We are able to classify lens candidates by eye at a rate of a few seconds per system, suggesting that a future 1000 square degree imaging survey containing 10^7 BRGs, and some 10^4 lenses, could be successfully, and reproducibly, searched in a modest amount of time. [Abridged]

preprint2008arXiv

Enhanced lensing rate by clustering of massive galaxies: newly discovered systems in the SLACS fields

[Abridged] We exploit the clustering of massive galaxies to perform a high efficiency imaging search for gravitational lenses. Our dataset comprises 44 fields imaged by the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS), each of which is centered on a lens discovered by the Strong Lens ACS Survey (SLACS). We compare four different search methods: 1) automated detection with the HST Archive Galaxy-scale Gravitational Lens Survey (HAGGLeS) robot, 2) examining cutout images of bright galaxies (BGs) after subtraction of a smooth galaxy light distribution, 3) examining the unsubtracted BG cutouts, and 4) performing a full-frame visual inspection of the ACS images. We compute purity and completeness and consider investigator time for the four algorithms, using the main SLACS lenses as a testbed. The first and second algorithms perform the best. We present the four new lens systems discovered during this comprehensive search, as well as one other likely candidate. For each new lens we use the fundamental plane to estimate the lens velocity dispersion and predict, from the resulting lens geometry, the redshifts of the lensed sources. Two of these new systems are found in galaxy clusters, which include the SLACS lenses in the two respective fields. Overall we find that the enhanced lens abundance (30^{+24}_{-8} lenses/degree^2) is higher than expected for random fields (12^{+4}_{-2} lenses/degree^2 for the COSMOS survey). Additionally, we find that the gravitational lenses we detect are qualitatively different from those in the parent SLACS sample: this imaging survey is largely probing higher-redshift, and lower-mass, early-type galaxies.