Showing 1–50 of 129 results for author: Keeton, C R

Comparisons Between Resolved Star Formation Rate and Gas Tracers in the Strongly Lensed Galaxy SDSS J0901+1814 at Cosmic Noon

Authors: Qingxiang Chen, Chelsea E. Sharon, Hiddo S. Algera, Andrew J. Baker, Charles R. Keeton, Dieter Lutz, Daizhong Liu, Anthony J. Young, Amit Tagore, Jesus Rivera, Erin K. Hicks, Sahar S. Allam, Douglas L. Tucker

Abstract: We report new radio observations of SDSS J090122.37+181432.3, a strongly lensed star-forming galaxy at $z=2.26$. We image 1.4 GHz (L-band) and 3 GHz (S-band) continuum using the VLA and 1.2 mm (band 6) continuum with ALMA, in addition to the CO(7-6) and CI(${\rm ^3P_2\rightarrow ^3\!P_1}$) lines, all at $\lesssim1.^{\prime\prime}7$ resolution. Based on the VLA integrated flux densities, we decompo… ▽ More We report new radio observations of SDSS J090122.37+181432.3, a strongly lensed star-forming galaxy at $z=2.26$. We image 1.4 GHz (L-band) and 3 GHz (S-band) continuum using the VLA and 1.2 mm (band 6) continuum with ALMA, in addition to the CO(7-6) and CI(${\rm ^3P_2\rightarrow ^3\!P_1}$) lines, all at $\lesssim1.^{\prime\prime}7$ resolution. Based on the VLA integrated flux densities, we decompose the radio spectrum into its free-free (FF) and non-thermal components. The infrared-radio correlation (IRRC) parameter $q_{\rm TIR}=2.65_{-0.31}^{+0.24}$ is consistent with expectations for star forming galaxies. We obtain radio continuum-derived SFRs that are free of dust extinction, finding $\rm {620}_{-220}^{+280}\,M_\odot\,yr^{-1}$, $\rm {230}_{-160}^{+570}\,M_\odot\,yr^{-1}$, and $\rm {280}_{-120}^{+460}\,M_\odot\,yr^{-1}$ from the FF emission, non-thermal emission, and when accounting for both emission processes, respectively, in agreement with previous results. We estimate the gas mass from the CI(${\rm ^3P_2\rightarrow ^3\!P_1}$) line as $M_{\rm gas}=(1.2\pm0.2)\times10^{11}\,M_\odot$, which is consistent with prior CO(1-0)-derived gas masses. Using our new IR and radio continuum data to map the SFR, we assess the dependence of the Schmidt-Kennicutt relation on choices of SFR and gas tracer for $\sim{\rm kpc}$ scales. The different SFR tracers yield different slopes, with the IR being the steepest, potentially due to highly obscured star formation in J0901. The radio continuum maps have the lowest slopes and overall fidelity for mapping the SFR, despite producing consistent total SFRs. We also find that the Schmidt-Kennicutt relation slope is flattest when using CO(7-6) or CI(${\rm ^3P_2\rightarrow ^3\!P_1}$) to trace gas mass, suggesting that those transitions are not suitable for tracing the bulk molecular gas in galaxies like J0901. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 28 pages, 13 figures, 9 tables

Statistical Bias in the Hubble Constant and Mass Power Law Slope for Mock Strong Lenses

Authors: Dilys Ruan, Charles R. Keeton

Abstract: Strong gravitational lensing offers constraints on the Hubble constant that are independent of other methods. However, those constraints are subject to uncertainties in lens models. Previous studies suggest that using an elliptical power law + external shear (EPL+XS) for the lensing galaxy can yield results that are precise but inaccurate. We examine such models by generating and fitting mock lens… ▽ More Strong gravitational lensing offers constraints on the Hubble constant that are independent of other methods. However, those constraints are subject to uncertainties in lens models. Previous studies suggest that using an elliptical power law + external shear (EPL+XS) for the lensing galaxy can yield results that are precise but inaccurate. We examine such models by generating and fitting mock lenses which produces multiple images of a background quasar-like point source. Despite using the same model for input and output, we find statistical bias in the Hubble constant on the order of 3% to 5%, depending on whether the elliptical lenses have noise or not. The phase space distribution has a `flared' shape that causes the mass power law slope to be underestimated and the Hubble constant to be overestimated. The bias varies with image configuration, which we quantify through annulus length between images with the first and second time delays ($Δr_{1,2}$). The statistical bias is worse for configurations that have narrow annuli (e.g., symmetric cross configurations). Assuming a source at redshift 2.0 and an EPL+XS lens at redshift 0.3, we find that the bias can be reduced, but not eliminated, if we limit the sample to systems with annulus lengths $Δr_{1,2} \gtrsim 0.3$ arcsec. As lens samples grow, it may be helpful to prioritize this range of image configurations for follow-up observation and analysis. △ Less

Submitted 1 July, 2024; v1 submitted 28 September, 2023; originally announced September 2023.

A New Method for the Reconstruction of Strongly Lensed Galaxies with Resolved Kinematics

Authors: Anthony J. Young, Charles R. Keeton, Andrew J. Baker

Abstract: Integral field spectroscopy of high-redshift galaxies has become a powerful tool for understanding their dynamics and evolutionary states. However, in the case of gravitationally lensed systems, it has proved difficult to model both lensing and intrinsic kinematics in a way that takes full advantage of the information available in the spectral domain. In this paper, we introduce a new method for p… ▽ More Integral field spectroscopy of high-redshift galaxies has become a powerful tool for understanding their dynamics and evolutionary states. However, in the case of gravitationally lensed systems, it has proved difficult to model both lensing and intrinsic kinematics in a way that takes full advantage of the information available in the spectral domain. In this paper, we introduce a new method for pixel-based source reconstruction that alters standard regularization schemes for two-dimensional data in a way that leverages kinematic information in a physically motivated but flexible fashion, and that is better suited to the three-dimensional nature of integral field data. To evaluate the performance of this method, we compare its results to those of a more traditional two-dimensional non-parametric approach using mock ALMA observations of a typical high-redshift dusty star-forming galaxy. We find that 3D regularization applied to an entire data cube reconstructs a source's intensity and velocity structure more accurately than 2D regularization applied to separate velocity channels. Cubes reconstructed with 3D regularization also have more uniform noise and resolution properties and are less sensitive to the signal-to-noise ratio of individual velocity channels than the results of 2D regularization. Our new approach to modeling integral field observations of lensed systems can be implemented without making restrictive a priori assumptions about intrinsic kinematics, and opens the door to new observing strategies that prioritize spectral resolution over spatial resolution (e.g., for multi-configuration arrays like ALMA). △ Less

Submitted 11 April, 2022; v1 submitted 28 February, 2022; originally announced February 2022.

Comments: Updated to match published version; 22 pages, 11 figures

Journal ref: ApJ, 2022, 929, 6

Deep XMM-Newton Observations of an X-ray Weak, Broad Absorption Line Quasar at $z=6.5$

Authors: Jinyi Yang, Xiaohui Fan, Feige Wang, Giorgio Lanzuisi, Riccardo Nanni, Massimo Cappi, George Chartas, Mauro Dadina, Roberto Decarli, Xiangyu Jin, Charles R. Keeton, Bram P. Venemans, Fabian Walter, Ran Wang, Xue-Bing Wu, Minghao Yue, Ann Zabludoff

Abstract: We report X-ray observations of the most distant known gravitationally lensed quasar, J0439+1634 at $z=6.52$, which is also a broad absorption line (BAL) quasar, using the XMM-Newton Observatory. With a 130 ks exposure, the quasar is significantly detected as a point source at the optical position with a total of 358$^{+19}_{-19}$ net counts using the EPIC instrument. By fitting a power-law plus G… ▽ More We report X-ray observations of the most distant known gravitationally lensed quasar, J0439+1634 at $z=6.52$, which is also a broad absorption line (BAL) quasar, using the XMM-Newton Observatory. With a 130 ks exposure, the quasar is significantly detected as a point source at the optical position with a total of 358$^{+19}_{-19}$ net counts using the EPIC instrument. By fitting a power-law plus Galactic absorption model to the observed spectra, we obtain a spectral slope of $Γ=1.45^{+0.10}_{-0.09}$. The derived optical-to-X-ray spectral slope $α_{\rm{ox}}$ is $-2.07^{+0.01}_{-0.01}$, suggesting that the X-ray emission of J0439+1634 is weaker by a factor of 18 than the expectation based on its 2500 Angstrom luminosity and the average $α_{\rm{ox}}$ vs. luminosity relationship. This is the first time that an X-ray weak BAL quasar at $z>6$ has been observed spectroscopically. Its X-ray weakness is consistent with the properties of BAL quasars at lower redshift. By fitting a model including an intrinsic absorption component, we obtain intrinsic column densities of $N_{\rm{H}}=2.8^{+0.7}_{-0.6}\times10^{23}\,\rm{cm}^{-2}$ and $N_{\rm{H}}= 4.3^{+1.8}_{-1.5}\times10^{23}\,\rm{cm}^{-2}$, assuming a fixed $Γ$ of 1.9 and a free $Γ$, respectively. The intrinsic rest-frame 2--10 keV luminosity is derived as $(9.4-15.1)\times10^{43}\,\rm{erg\,s}^{-1}$, after correcting for lensing magnification ($μ=51.3$). The absorbed power-law model fitting indicates that J0439+1634 is the highest redshift obscured quasar with a direct measurement of the absorbing column density. The intrinsic high column density absorption can reduce the X-ray luminosity by a factor of $3-7$, which also indicates that this quasar could be a candidate of intrinsically X-ray weak quasar. △ Less

Submitted 20 December, 2021; originally announced December 2021.

Comments: 8 pages, 3 figures, 1 table; Accepted for publication in ApJL

ALMA Observations of the Sub-kpc Structure of the Host Galaxy of a z= 6.5 Lensed Quasar: A Rotationally-Supported Hyper-Starburst System at the Epoch of Reionization

Authors: Minghao Yue, Jinyi Yang, Xiaohui Fan, Feige Wang, Justin Spilker, Iskren Y. Georgiev, Charles R. Keeton, Katrina C. Litke, Daniel P. Marrone, Fabian Walter, Ran Wang, Xue-Bing Wu, Bram P. Venemans, Ann Zabludoff

Abstract: We report ALMA observations of the dust continuum and {\cii} emission of the host galaxy of J0439+1634, a gravitationally lensed quasar at $z=6.5$. Gravitational lensing boosts the source-plane resolution to $\sim0\farcs15$ $(\sim0.8\text{ kpc})$. The lensing model derived from the ALMA data is consistent with the fiducial model in \citet{fan19} based on {\it HST} imaging. The host galaxy of J0439… ▽ More We report ALMA observations of the dust continuum and {\cii} emission of the host galaxy of J0439+1634, a gravitationally lensed quasar at $z=6.5$. Gravitational lensing boosts the source-plane resolution to $\sim0\farcs15$ $(\sim0.8\text{ kpc})$. The lensing model derived from the ALMA data is consistent with the fiducial model in \citet{fan19} based on {\it HST} imaging. The host galaxy of J0439+1634 can be well-fitted by a Sérsic profile consistent with an exponential disk, both in the far-infrared (FIR) continuum and the {\cii} emission. The overall magnification is $4.53\pm0.05$ for the continuum and $3.44\pm0.05$ for the {\cii} line. The host galaxy of J0439+1634 is a compact ultra-luminous infrared galaxy, with a total star formation rate (SFR) of $1.56\times10^{3}M_\odot/\text{year}$ after correcting for lensing and an effective radius of $0.74$ kpc. The resolved regions in J0439+1634 follow the ``{\cii} deficit," where the {\cii}-to-FIR ratio decreases with FIR surface brightness. The reconstructed velocity field of J0439+1634 appears to be rotation-like. The maximum line-of-sight rotation velocity of 130 km/s at a radius of 2 kpc. However, our data cannot be fit by an axisymmetric thin rotating disk, and the inclination of the rotation axis, $i$, remains unconstrained. We estimate the dynamical mass of the host galaxy to be $7.9\sin^{-2}(i)\times10^{9}M_\odot$. J0439+1634 is likely to have a high gas-mass fraction and an oversized SMBH compared to local relations. The SFR of J0439+1634 reaches the maximum possible values, and the SFR surface density is close to the highest value seen in any star-forming galaxy currently known in the universe. △ Less

Submitted 15 June, 2021; v1 submitted 14 June, 2021; originally announced June 2021.

Comments: 14 pages, 7 figures. Accepted by ApJ

Systematic vs. Statistical Uncertainties in Masses and Magnifications of the Hubble Frontier Fields

Authors: Catie A Raney, Charles R Keeton, Sean Brennan, Hsin Fan

Abstract: The Hubble Frontier Fields data, along with multiple data sets obtained by other telescopes, have provided some of the most extensive constraints on cluster lenses to date. Multiple lens modeling teams analyzed the fields and made public a number of deliverables. By comparing these results, we can then undertake a unique and vital test of the state of cluster lens modeling. Specifically, we see ho… ▽ More The Hubble Frontier Fields data, along with multiple data sets obtained by other telescopes, have provided some of the most extensive constraints on cluster lenses to date. Multiple lens modeling teams analyzed the fields and made public a number of deliverables. By comparing these results, we can then undertake a unique and vital test of the state of cluster lens modeling. Specifically, we see how well the different teams can reproduce similar magnifications and mass profiles. We find that the circularly averaged mass profiles of the fields are remarkably constrained (scatter <5%) at distances of 1 arcmin from the cluster core, yet magnifications can vary significantly. Averaged across the six fields, we find a bias of -6% (-17%) and a scatter of ~40% (~65%) at a modest magnification of 3 (10). Statistical errors reported by individual teams are often significantly smaller than the differences among all the teams, indicating the importance of continued systematics studies in cluster lensing. △ Less

Submitted 13 April, 2020; originally announced April 2020.

Comments: 25 pages, 21 figures. Accepted for publication in MNRAS

Exploring Effects on Magnifications due to Line-of-Sight Galaxies in the Hubble Frontier Fields

Authors: Catie A. Raney, Charles R. Keeton, Sean Brennan

Abstract: Cluster lensing has become an important tool in the search for high redshift galaxies through its ability to magnify sources. In order to determine the intrinsic properties of these galaxies, lensing mass models must be constructed to determine the magnification of the images. These models are traditionally two-dimensional, focusing on the mass within the cluster and either ignoring or approximati… ▽ More Cluster lensing has become an important tool in the search for high redshift galaxies through its ability to magnify sources. In order to determine the intrinsic properties of these galaxies, lensing mass models must be constructed to determine the magnification of the images. These models are traditionally two-dimensional, focusing on the mass within the cluster and either ignoring or approximating any contribution from line-of-sight galaxies. In this paper, we present the first full set of three-dimensional mass models of the six Hubble Frontier Fields and use them to test for systematic biases in magnifications due to using the traditional 2D approach. We find that omitting foreground or background galaxies causes image position offsets between 0.1-0.4", a non-negligible fraction of the typical 0.3-0.7" residuals of current state-of-the-art models. We also find that median image magnifications can shift by up to 6%, though it is dependent on the field. This can be alleviated in some cases by approximating the mass in the lensing plane, but a 5% magnification bias still exists in other cases; image position offsets are also improved, but are still present at 0.10". △ Less

Submitted 5 November, 2019; originally announced November 2019.

Comments: 26 pages, 9 figures, 14 tables. Accepted for publication in MNRAS

Far-Infrared Properties of the Bright, Gravitationally Lensed Quasar J0439+1634 at z=6.5

Authors: Jinyi Yang, Bram Venemans, Feige Wang, Xiaohui Fan, Mladen Novak, Roberto Decarli, Fabian Walter, Minghao Yue, Emmanuel Momjian, Charles R. Keeton, Ran Wang, Ann Zabludoff, Xue-Bing Wu, Fuyan Bian

Abstract: We present IRAM/NOEMA, JCMT/SCUBA-2 and VLA observations of the most distant known gravitationally lensed quasar J0439+1634 at z = 6.5. We detect strong dust emission, [CII] 158 $μ$m, [CI] 369 $μ$m, [OI] 146 $μ$m, CO(6-5), CO(7-6), CO(9-8), CO(10-9), H$_{\rm 2}$O $3_{\rm 1,2}-2_{\rm 2,1}$, and H$_{\rm 2}$O $3_{\rm 2,1}-3_{\rm 1,2}$ lines as well as a weak radio continuum. The strong [CII] line yie… ▽ More We present IRAM/NOEMA, JCMT/SCUBA-2 and VLA observations of the most distant known gravitationally lensed quasar J0439+1634 at z = 6.5. We detect strong dust emission, [CII] 158 $μ$m, [CI] 369 $μ$m, [OI] 146 $μ$m, CO(6-5), CO(7-6), CO(9-8), CO(10-9), H$_{\rm 2}$O $3_{\rm 1,2}-2_{\rm 2,1}$, and H$_{\rm 2}$O $3_{\rm 2,1}-3_{\rm 1,2}$ lines as well as a weak radio continuum. The strong [CII] line yields a systemic redshift of the host galaxy to be z=6.5188$\pm$0.0001. The magnification makes J0439+1634 the far-infrared (FIR) brightest quasar at z > 6 known, with the brightest [CII] line yet detected at this redshift. The FIR luminosity is (3.4$\pm$0.2)$\times$10$^{13}$ $μ^{-1}$ $L_{\odot}$, where $μ$ $\sim$ 2.6 - 6.6 is the magnification of the host galaxy, estimated based on the lensing configuration from HST imaging. We estimate the dust mass to be (2.2$\pm$0.1)$\times$10$^{9}$ $μ^{-1}$ $M_{\odot}$. The CO Spectral Line Energy Distribution using four CO lines is best fit by a two-component model of the molecular gas excitation. The estimates of molecular gas mass derived from CO lines and atomic carbon mass are consistent, in the range of 3.9 - 8.9 $\times$10$^{10} μ^{-1}$ $M_{\odot}$. The [CII]/[CI], [CII]/CO, and [OI]/[CII] line luminosity ratios suggest a photodissociation region model with more than one component. The ratio of H$_{\rm 2}$O $3_{\rm 2,1}-3_{\rm 1,2}$ line luminosity to $L_{\rm TIR}$ is consistent with values in local and high redshift ultra-/hyper-luminous infrared galaxies. The VLA observations reveal an unresolved radio continuum source, and indicate that J0439+1634 is a radio quiet quasar with R = 0.05 - 0.17. △ Less

Submitted 30 June, 2019; originally announced July 2019.

Comments: 13 pages, 6 figures. Accepted by ApJ

Magnification Cross Sections for the Elliptic Umbilic Caustic Surface

Authors: Amir B. Aazami, Charles R. Keeton, Arlie O. Petters

Abstract: We show that the asymptotic scaling of the magnification volume cross section corresponding to an elliptic umbilic caustic surface is $μ^{-2.5}$ in the two-image region and $μ^{-2}$ in the four-image region, where $μ$ is the total unsigned magnification. We show that the asymptotic scaling of the magnification volume cross section corresponding to an elliptic umbilic caustic surface is $μ^{-2.5}$ in the two-image region and $μ^{-2}$ in the four-image region, where $μ$ is the total unsigned magnification. △ Less

Submitted 26 July, 2019; v1 submitted 28 May, 2019; originally announced May 2019.

Comments: v2: 6 pages, 6 figures, references added, expanded discussion in Section II

Journal ref: Universe 2019, 5(7), 161

Resolved Molecular Gas and Star Formation Properties of the Strongly Lensed z=2.26 Galaxy SDSS J0901+1814

Authors: Chelsea E. Sharon, Amitpal S. Tagore, Andrew J. Baker, Jesus Rivera, Charles R. Keeton, Dieter Lutz, Reinhard Genzel, David J. Wilner, Erin K. S. Hicks, Sahar S. Allam, Douglas L. Tucker

Abstract: We present ~1" resolution (~2 kpc in the source plane) observations of the CO(1-0), CO(3-2), Halpha, and [N II] lines in the strongly-lensed z=2.26 star-forming galaxy SDSS J0901+1814. We use these observations to constrain the lensing potential of a foreground group of galaxies, and our source-plane reconstructions indicate that SDSS J0901+1814 is a nearly face-on (i~30 degrees) massive disk with… ▽ More We present ~1" resolution (~2 kpc in the source plane) observations of the CO(1-0), CO(3-2), Halpha, and [N II] lines in the strongly-lensed z=2.26 star-forming galaxy SDSS J0901+1814. We use these observations to constrain the lensing potential of a foreground group of galaxies, and our source-plane reconstructions indicate that SDSS J0901+1814 is a nearly face-on (i~30 degrees) massive disk with r_{1/2}>~4 kpc for its molecular gas. Using our new magnification factors (mu_tot~30), we find that SDSS J0901+1814 has a star formation rate (SFR) of 268^{+63}_{-61} M_sun/yr, M_gas=(1.6^{+0.3}_{-0.2})x10^11x(alpha_CO/4.6) M_sun, and M_star=(9.5^{+3.8}_{-2.8})x10^10 M_sun, which places it on the star-forming galaxy "main sequence." We use our matched high-angular resolution gas and SFR tracers (CO and Halpha, respectively) to perform a spatially resolved (pixel-by-pixel) analysis of SDSS J0901+1814 in terms of the Schmidt-Kennicutt relation. After correcting for the large fraction of obscured star formation (SFR_Halpha/SFR_TIR=0.054^{+0.015}_{-0.014}), we find SDSS J0901+1814 is offset from "normal" star-forming galaxies to higher star formation efficiencies independent of assumptions for the CO-to-H_2 conversion factor. Our mean best-fit index for the Schmidt-Kennicutt relation for SDSS J0901+1814, evaluated with different CO lines and smoothing levels, is n=1.54+/-0.13; however, the index may be affected by gravitational lensing, and we find n=1.24+/-0.02 when analyzing the source-plane reconstructions. While the Schmidt-Kennicutt index largely appears unaffected by which of the two CO transitions we use to trace the molecular gas, the source-plane reconstructions and dynamical modeling suggest that the CO(1-0) emission is more spatially extended than the CO(3-2) emission. △ Less

Submitted 23 May, 2019; originally announced May 2019.

Comments: ApJ accepted; 40 pages, 24 figures, 8 tables

Dark Matter Science in the Era of LSST

Authors: Keith Bechtol, Alex Drlica-Wagner, Kevork N. Abazajian, Muntazir Abidi, Susmita Adhikari, Yacine Ali-Haïmoud, James Annis, Behzad Ansarinejad, Robert Armstrong, Jacobo Asorey, Carlo Baccigalupi, Arka Banerjee, Nilanjan Banik, Charles Bennett, Florian Beutler, Simeon Bird, Simon Birrer, Rahul Biswas, Andrea Biviano, Jonathan Blazek, Kimberly K. Boddy, Ana Bonaca, Julian Borrill, Sownak Bose, Jo Bovy , et al. (155 additional authors not shown)

Abstract: Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We… ▽ More Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 11 pages, 2 figures, Science Whitepaper for Astro 2020, more information at https://lsstdarkmatter.github.io

Probing the Fundamental Nature of Dark Matter with the Large Synoptic Survey Telescope

Authors: Alex Drlica-Wagner, Yao-Yuan Mao, Susmita Adhikari, Robert Armstrong, Arka Banerjee, Nilanjan Banik, Keith Bechtol, Simeon Bird, Kimberly K. Boddy, Ana Bonaca, Jo Bovy, Matthew R. Buckley, Esra Bulbul, Chihway Chang, George Chapline, Johann Cohen-Tanugi, Alessandro Cuoco, Francis-Yan Cyr-Racine, William A. Dawson, Ana Díaz Rivero, Cora Dvorkin, Denis Erkal, Christopher D. Fassnacht, Juan García-Bellido, Maurizio Giannotti , et al. (75 additional authors not shown)

Abstract: Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We d… ▽ More Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ . △ Less

Submitted 24 April, 2019; v1 submitted 4 February, 2019; originally announced February 2019.

Comments: 96 pages, 22 figures, 1 table

Report number: FERMILAB-PUB-19-048-A-AE

The Discovery of a Gravitationally Lensed Quasar at z = 6.51

Authors: Xiaohui Fan, Feige Wang, Jinyi Yang, Charles R. Keeton, Minghao Yue, Ann Zabludoff, Fuyan Bian, Marco Bonaglia, Iskren Y. Georgiev, Joseph F. Hennawi, Jiangtao Li, Ian D. McGreer, Rohan Naidu, Fabio Pacucci, Sebastian Rabien, David Thompson, Bram Venemans, Fabian Walter, Xue-Bing Wu, Ran Wang

Abstract: Strong gravitational lensing provides a powerful probe of the physical properties of quasars and their host galaxies. A high fraction of the most luminous high-redshift quasars was predicted to be lensed due to magnification bias. However, no multiple imaged quasar was found at z>5 in previous surveys. We report the discovery of J043947.08+163415.7, a strongly lensed quasar at z=6.51, the first su… ▽ More Strong gravitational lensing provides a powerful probe of the physical properties of quasars and their host galaxies. A high fraction of the most luminous high-redshift quasars was predicted to be lensed due to magnification bias. However, no multiple imaged quasar was found at z>5 in previous surveys. We report the discovery of J043947.08+163415.7, a strongly lensed quasar at z=6.51, the first such object detected at the epoch of reionization, and the brightest quasar yet known at z>5. High-resolution HST imaging reveals a multiple imaged system with a maximum image separation theta ~ 0.2", best explained by a model of three quasar images lensed by a low luminosity galaxy at z~0.7, with a magnification factor of ~50. The existence of this source suggests that a significant population of strongly lensed, high redshift quasars could have been missed by previous surveys, as standard color selection techniques would fail when the quasar color is contaminated by the lensing galaxy. △ Less

Submitted 28 October, 2018; originally announced October 2018.

Comments: 8 pages, 4 figures, submitted to ApJL

Quantifying the power spectrum of small-scale structure in semi-analytic galaxies

Authors: Sean Brennan, Andrew J. Benson, Francis-Yan Cyr-Racine, Charles R. Keeton, Leonidas A. Moustakas, Anthony R. Pullen

Abstract: In the cold dark matter (CDM) picture of structure formation, galaxy mass distributions are predicted to have a considerable amount of structure on small scales. Strong gravitational lensing has proven to be a useful tool for studying this small-scale structure. Much of the attention has been given to detecting individual dark matter subhalos through lens modeling, but recent work has suggested th… ▽ More In the cold dark matter (CDM) picture of structure formation, galaxy mass distributions are predicted to have a considerable amount of structure on small scales. Strong gravitational lensing has proven to be a useful tool for studying this small-scale structure. Much of the attention has been given to detecting individual dark matter subhalos through lens modeling, but recent work has suggested that the full population of subhalos could be probed using a power spectrum analysis. In this paper we quantify the power spectrum of small-scale structure in simulated galaxies, with the goal of understanding theoretical predictions and setting the stage for using measurements of the power spectrum to test dark matter models. We use a sample of simulated galaxies generated from the \texttt{Galacticus} semi-analytic model to determine the power spectrum distribution first in the CDM paradigm and then in a warm dark matter scenario. We find that a measurement of the slope and amplitude of the power spectrum on galaxy strong lensing scales ($k\sim 1$ kpc$^{-1}$) could be used to distinguish between CDM and alternate dark matter models, especially if the most massive subhalos can be directly detected via gravitational imaging. △ Less

Submitted 29 August, 2019; v1 submitted 10 August, 2018; originally announced August 2018.

Comments: 8 pages, 10 figures

The Atacama Cosmology Telescope: CO(J = 3 - 2) mapping and lens modeling of an ACT-selected dusty star-forming galaxy

Authors: J. Rivera, A. J. Baker, P. A. Gallardo, M. Gralla, A. I. Harris, K. M. Huffenberger, J. P. Hughes, C. R. Keeton, C. H. Lopez-Caraballo, T. A. Marriage, B. Partridge, J. Sievers, A. S. Tagore, F. Walter, A. Weiss, E. J. Wollack

Abstract: We report Northern Extended Millimeter Array (NOEMA) CO($J = 3 - 2$) observations of the dusty star-forming galaxy ACT-S\,J020941+001557 at $z = 2.5528$, which was detected as an unresolved source in the Atacama Cosmology Telescope (ACT) equatorial survey. Our spatially resolved spectral line data support the derivation of a gravitational lens model from 37 independent velocity channel maps using… ▽ More We report Northern Extended Millimeter Array (NOEMA) CO($J = 3 - 2$) observations of the dusty star-forming galaxy ACT-S\,J020941+001557 at $z = 2.5528$, which was detected as an unresolved source in the Atacama Cosmology Telescope (ACT) equatorial survey. Our spatially resolved spectral line data support the derivation of a gravitational lens model from 37 independent velocity channel maps using a pixel-based algorithm, from which we infer a velocity-dependent magnification factor $μ\approx 7-22$ with a luminosity-weighted mean $\left<μ\right>\approx 13$. The resulting source-plane reconstruction is consistent with a rotating disk, although other scenarios cannot be ruled out by our data. After correction for lensing, we derive a line luminosity $L^{\prime}_{\rm CO(3-2)}= (5.53\pm 0.69) \times 10^{10}\,{\rm \,K\,km\,s^{-1}\,pc^{2}}$, a cold gas mass $M_{\rm gas}= (3.86 \pm 0.33) \times 10^{10}\,M_{\odot}$, a dynamical mass $M_{\rm dyn}\,{\rm sin}^2\,i = 3.9^{+1.8}_{-1.5} \times 10^{10}\,M_{\odot}$, and a gas mass fraction $f_{\rm gas}\,{\rm csc}^2\,i = 1.0^{+0.8}_{-0.4}$. The line brightness temperature ratio of $r_{3,1}\approx 1.6$ relative to a Green Bank Telescope CO($J=1-0$) detection may be elevated by a combination of external heating of molecular clouds, differential lensing, and/or pointing errors. △ Less

Submitted 3 June, 2019; v1 submitted 23 July, 2018; originally announced July 2018.

Comments: 8 pages, 5 figures, accepted to ApJ

Beyond subhalos: Probing the collective effect of the Universe's small-scale structure with gravitational lensing

Authors: Francis-Yan Cyr-Racine, Charles R. Keeton, Leonidas A. Moustakas

Abstract: Gravitational lensing has emerged as a powerful probe of the matter distribution on subgalactic scales, which itself may contain important clues about the fundamental origins and properties of dark matter. Broadly speaking, two different approaches have been taken in the literature to map the small-scale structure of the Universe using strong lensing, with one focused on measuring the position and… ▽ More Gravitational lensing has emerged as a powerful probe of the matter distribution on subgalactic scales, which itself may contain important clues about the fundamental origins and properties of dark matter. Broadly speaking, two different approaches have been taken in the literature to map the small-scale structure of the Universe using strong lensing, with one focused on measuring the position and mass of a small number of discrete massive subhalos appearing close in projection to lensed images, and the other focused on detecting the collective effect of all the small-scale structure between the lensed source and the observer. In this paper, we follow the latter approach and perform a detailed study of the sensitivity of galaxy-scale gravitational lenses to the ensemble properties of small-scale structure. As in some previous studies, we adopt the language of the substructure power spectrum to characterize the statistical properties of the small-scale density field. We present a comprehensive theory that treats lenses with extended sources as well as those with time-dependent compact sources (such as quasars) in a unified framework for the first time. Our approach uses mode functions to provide both computational advantages and insights about couplings between the lens and source. The goal of this paper is to develop the theory and gain the intuition necessary to understand how the sensitivity to the substructure power spectrum depends on the source and lens properties, with the eventual aim of identifying the most promising targets for such studies. △ Less

Submitted 27 August, 2019; v1 submitted 20 June, 2018; originally announced June 2018.

Comments: 24 pages + appendices, 12 figures. v2: matches published version

Journal ref: Phys. Rev. D 100, 023013 (2019)

Probing the nature of dark matter by forward modeling flux ratios in strong gravitational lenses

Authors: Daniel Gilman, Simon Birrer, Tommaso Treu, Charles R. Keeton, Anna Nierenberg

Abstract: The free streaming length of dark matter particles determines the abundance of structure on sub-galactic scales. We present a statistical technique, amendable to any parameterization of subhalo density profile and mass function, to probe dark matter on these scales with quadrupole image lenses. We consider a warm dark matter particle with a mass function characterized by a normalization and free s… ▽ More The free streaming length of dark matter particles determines the abundance of structure on sub-galactic scales. We present a statistical technique, amendable to any parameterization of subhalo density profile and mass function, to probe dark matter on these scales with quadrupole image lenses. We consider a warm dark matter particle with a mass function characterized by a normalization and free streaming scale $m_{\rm{hm}}$. We forecast bounds on dark matter warmth for 120-180 lenses, attainable with future surveys, at typical lens (source) redshifts of 0.5 (1.5) in early-type galaxies with velocity dispersions of 220-270 km/sec. We demonstrate that limits on $m_{\rm{hm}}$ deteriorate rapidly with increasing uncertainty in image fluxes, underscoring the importance of precise measurements and accurate lens models. For our forecasts, we assume the deflectors in the lens sample do not exhibit complex morphologies, so we neglect systematic errors in their modeling. Omitting the additional signal from line of sight halos, our constraints underestimate the true power of the method. Assuming cold dark matter, for a low normalization, corresponding the destruction of all subhalos within the host scale radius, we forecast $2σ$ bounds on $m_{\rm{hm}}$ (thermal relic mass) of $10^{7.5} \ (5.0)$, $10^{8} \ (3.6)$, and $10^{8.5} \ (2.7) \ M_{\odot} \left(\rm{keV}\right)$ for flux errors of $2\%$, $4\%$, and $8\%$. With a higher normalization, these constraints improve to $10^{7.2} \ (6.6)$, $10^{7.5} \ (5.3) $, and $10^{7.8} \ (4.3) \ M_{\odot} \left(\rm{keV}\right)$ with 120 systems. We are also able to measure the normalization of the mass function, which has implications for baryonic feedback models and tidal stripping. △ Less

Submitted 16 August, 2018; v1 submitted 13 December, 2017; originally announced December 2017.

Comments: accepted to MNRAS; added a new simulation to illustrate how constraints on free streaming length scale with mass function normalization; 17 pages, 12 figures, 1 table, 5 appendices

A Spectroscopic Survey of the Fields of 28 Strong Gravitational Lenses: The Group Catalog

Authors: Michelle L. Wilson, Ann I. Zabludoff, S. Mark Ammons, Ivelina G. Momcheva, Kurtis A. Williams, Charles R. Keeton

Abstract: With a large, unique spectroscopic survey in the fields of 28 galaxy-scale strong gravitational lenses, we identify groups of galaxies in the 26 adequately-sampled fields. Using a group finding algorithm, we find 210 groups with at least five member galaxies; the median number of members is eight. Our sample spans redshifts of 0.04 $\le z_{grp} \le$ 0.76 with a median of 0.31, including 174 groups… ▽ More With a large, unique spectroscopic survey in the fields of 28 galaxy-scale strong gravitational lenses, we identify groups of galaxies in the 26 adequately-sampled fields. Using a group finding algorithm, we find 210 groups with at least five member galaxies; the median number of members is eight. Our sample spans redshifts of 0.04 $\le z_{grp} \le$ 0.76 with a median of 0.31, including 174 groups with $0.1 < z_{grp} < 0.6$. Groups have radial velocity dispersions of 60 $\le σ_{grp} \le$ 1200 km s$^{-1}$ with a median of 350 km s$^{-1}$. We also discover a supergroup in field B0712+472 at $z =$ 0.29 consisting of three main groups. We recover groups similar to $\sim$ 85% of those previously reported in these fields within our redshift range of sensitivity and find 187 new groups with at least five members. The properties of our group catalog, specifically 1) the distribution of $σ_{grp}$, 2) the fraction of all sample galaxies that are group members, and 3) the fraction of groups with significant substructure, are consistent with those for other catalogs. The distribution of group virial masses agrees well with theoretical expectations. Of the lens galaxies, 12 of 26 (46%) (B1422+231, B1600+434, B2114+022, FBQS J0951+2635, HE0435-1223, HST J14113+5211, MG0751+2716, MGJ1654+1346, PG 1115+080, Q ER 0047-2808, RXJ1131-1231, and WFI J2033-4723) are members of groups with at least five galaxies, and one more (B0712+472) belongs to an additional, visually identified group candidate. There are groups not associated with the lens that still are likely to affect the lens model; in six of 25 (24%) fields (excluding the supergroup), there is at least one massive ($σ_{grp} \ge$ 500 km s$^{-1}$) group or group candidate projected within 2$^{\prime}$ of the lens. △ Less

Submitted 26 October, 2017; originally announced October 2017.

Comments: 87 pages, 8 figures, a version of this was published in ApJ

Journal ref: 2016, ApJ, 833, 194

A Spectroscopic Survey of the Fields of 28 Strong Gravitational Lenses: Implications for $H_0$

Authors: Michelle L. Wilson, Ann I. Zabludoff, Charles R. Keeton, Kenneth C. Wong, Kurtis A. Williams, K. Decker French, Ivelina G. Momcheva

Abstract: Strong gravitational lensing provides an independent measurement of the Hubble parameter ($H_0$). One remaining systematic is a bias from the additional mass due to a galaxy group at the lens redshift or along the sightline. We quantify this bias for more than 20 strong lenses that have well-sampled sightline mass distributions, focusing on the convergence $κ$ and shear $γ$. In 23% of these fields… ▽ More Strong gravitational lensing provides an independent measurement of the Hubble parameter ($H_0$). One remaining systematic is a bias from the additional mass due to a galaxy group at the lens redshift or along the sightline. We quantify this bias for more than 20 strong lenses that have well-sampled sightline mass distributions, focusing on the convergence $κ$ and shear $γ$. In 23% of these fields, a lens group contributes a $\ge$1% convergence bias; in 57%, there is a similarly significant line-of-sight group. For the nine time delay lens systems, $H_0$ is overestimated by 11$^{+3}_{-2}$% on average when groups are ignored. In 67% of fields with total $κ\ge$ 0.01, line-of-sight groups contribute $\gtrsim 2\times$ more convergence than do lens groups, indicating that the lens group is not the only important mass. Lens environment affects the ratio of four (quad) to two (double) image systems; all seven quads have lens groups while only three of 10 doubles do, and the highest convergences due to lens groups are in quads. We calibrate the $γ$-$κ$ relation: $\log(κ_{\rm{tot}}) = (1.94 \pm 0.34) \log(γ_{\rm{tot}}) + (1.31 \pm 0.49)$ with a rms scatter of 0.34 dex. Shear, which, unlike convergence, can be measured directly from lensed images, can be a poor predictor of $κ$; for 19% of our fields, $κ$ is $\gtrsim 2γ$. Thus, accurate cosmology using strong gravitational lenses requires precise measurement and correction for all significant structures in each lens field. △ Less

Submitted 26 October, 2017; originally announced October 2017.

Comments: 34 pages, 11 figures, accepted for publication in ApJ

Discovery of the first quadruple gravitationally lensed quasar candidate with Pan-STARRS

Authors: C. T. Berghea, G. J. Nelson, C. E. Rusu, C. R. Keeton, R. P. Dudik

Abstract: We report the serendipitous discovery of the first gravitationally lensed quasar candidate from Pan-STARRS. The grizy images reveal four point-like images with magnitudes between 14.9 mag and 18.1 mag. The colors of the point sources are similar, and they are more consistent with quasars than with stars or galaxies. The lensing galaxy is detected in the izy bands, with an inferred photometric reds… ▽ More We report the serendipitous discovery of the first gravitationally lensed quasar candidate from Pan-STARRS. The grizy images reveal four point-like images with magnitudes between 14.9 mag and 18.1 mag. The colors of the point sources are similar, and they are more consistent with quasars than with stars or galaxies. The lensing galaxy is detected in the izy bands, with an inferred photometric redshift of ~0.6, lower than that of the point sources. We successfully model the system with a singular isothermal ellipsoid with shear, using the relative positions of the five objects as constraints. While the brightness ranking of the point sources is consistent with that of the model, we find discrepancies between the model-predicted and observed fluxes, likely due to microlensing by stars and millilensing due to the dark matter substructure. In order to fully confirm the gravitational lens nature of this system and add it to the small but growing number of the powerful probes of cosmology and astrophysics represented by quadruply lensed quasars, we require further spectroscopy and high-resolution imaging. △ Less

Submitted 11 July, 2017; v1 submitted 23 May, 2017; originally announced May 2017.

Comments: 29 pages, 7 figures, accepted in ApJ

Joint Strong and Weak Lensing Analysis of the Massive Cluster Field J0850+3604

Authors: Kenneth C. Wong, Catie Raney, Charles R. Keeton, Keiichi Umetsu, Ann I. Zabludoff, S. Mark Ammons, K. Decker French

Abstract: We present a combined strong and weak lensing analysis of the J085007.6+360428 (J0850) field, which was selected by its high projected concentration of luminous red galaxies and contains the massive cluster Zwicky 1953. Using Subaru/Suprime-Cam $BVR_{c}I_{c}i^{\prime}z^{\prime}$ imaging and MMT/Hectospec spectroscopy, we first perform a weak lensing shear analysis to constrain the mass distributio… ▽ More We present a combined strong and weak lensing analysis of the J085007.6+360428 (J0850) field, which was selected by its high projected concentration of luminous red galaxies and contains the massive cluster Zwicky 1953. Using Subaru/Suprime-Cam $BVR_{c}I_{c}i^{\prime}z^{\prime}$ imaging and MMT/Hectospec spectroscopy, we first perform a weak lensing shear analysis to constrain the mass distribution in this field, including the cluster at $z = 0.3774$ and a smaller foreground halo at $z = 0.2713$. We then add a strong lensing constraint from a multiply-imaged galaxy in the imaging data with a photometric redshift of $z \approx 5.03$. Unlike previous cluster-scale lens analyses, our technique accounts for the full three-dimensional mass structure in the beam, including galaxies along the line of sight. In contrast with past cluster analyses that use only lensed image positions as constraints, we use the full surface brightness distribution of the images. This method predicts that the source galaxy crosses a lensing caustic such that one image is a highly-magnified "fold arc", which could be used to probe the source galaxy's structure at ultra-high spatial resolution ($< 30$ pc). We calculate the mass of the primary cluster to be $\mathrm{M_{vir}} = 2.93_{-0.65}^{+0.71} \times 10^{15}~\mathrm{M_{\odot}}$ with a concentration of $\mathrm{c_{vir}} = 3.46_{-0.59}^{+0.70}$, consistent with the mass-concentration relation of massive clusters at a similar redshift. The large mass of this cluster makes J0850 an excellent field for leveraging lensing magnification to search for high-redshift galaxies, competitive with and complementary to that of well-studied clusters such as the HST Frontier Fields. △ Less

Submitted 25 June, 2017; v1 submitted 10 May, 2017; originally announced May 2017.

Comments: Accepted for publication in The Astrophysical Journal; 14 pages, 13 figures, 3 tables

Probing dark matter substructure in the gravitational lens HE0435-1223 with the WFC3 grism

Authors: A. M. Nierenberg, T. Treu, G. Brammer, A. H. G. Peter, C. D. Fassnacht, C. R. Keeton, C. S. Kochanek, K. B. Schmidt, D. Sluse, S. A. Wright

Abstract: Strong gravitational lensing provides a powerful test of Cold Dark Matter (CDM) as it enables the detection and mass measurement of low mass haloes even if they do not contain baryons. Compact lensed sources such as Active Galactic Nuclei (AGN) are particularly sensitive to perturbing subhalos, but their use as a test of CDM has been limited by the small number of systems which have significant ra… ▽ More Strong gravitational lensing provides a powerful test of Cold Dark Matter (CDM) as it enables the detection and mass measurement of low mass haloes even if they do not contain baryons. Compact lensed sources such as Active Galactic Nuclei (AGN) are particularly sensitive to perturbing subhalos, but their use as a test of CDM has been limited by the small number of systems which have significant radio emission which is extended enough avoid significant lensing by stars in the plane of the lens galaxy, and red enough to be minimally affected by differential dust extinction. Narrow-line emission is a promising alternative as it is also extended and, unlike radio, detectable in virtually all optically selected AGN lenses. We present first results from a WFC3 grism narrow-line survey of lensed quasars, for the quadruply lensed AGN HE0435-1223. Using a forward modelling pipeline which enables us to robustly account for spatial blending, we measure the [OIII] 5007 Å~ flux ratios of the four images. We find that the [OIII] fluxes and positions are well fit by a simple smooth mass model for the main lens. Our data rule out a $M_{600}>10^{8} (10^{7.2}) M_\odot$ NFW perturber projected within $\sim$1\farcs0 (0\farcs1) arcseconds of each of the lensed images, where $M_{600}$ is the perturber mass within its central 600 pc. The non-detection is broadly consistent with the expectations of $Λ$CDM for a single system. The sensitivity achieved demonstrates that powerful limits on the nature of dark matter can be obtained with the analysis of $\sim20$ narrow-line lenses. △ Less

Submitted 13 June, 2017; v1 submitted 18 January, 2017; originally announced January 2017.

Comments: Accepted for publication in MNRAS, 15 pages, 8 figures

Strong lensing signatures of luminous structure and substructure in early-type galaxies

Authors: Daniel Gilman, Adriano Agnello, Tommaso Treu, Charles R. Keeton, Anna M. Nierenberg

Abstract: The arrival times, positions, and fluxes of multiple images in strong lens systems can be used to infer the presence of dark subhalos in the deflector, and thus test predictions of cold dark matter models. However, gravitational lensing does not distinguish between perturbations to a smooth gravitational potential arising from baryonic and non-baryonic mass. In this work, we quantify the extent to… ▽ More The arrival times, positions, and fluxes of multiple images in strong lens systems can be used to infer the presence of dark subhalos in the deflector, and thus test predictions of cold dark matter models. However, gravitational lensing does not distinguish between perturbations to a smooth gravitational potential arising from baryonic and non-baryonic mass. In this work, we quantify the extent to which the stellar mass distribution of a deflector can reproduce flux ratio and astrometric anomalies typically associated with the presence of a dark matter subhalo. Using Hubble Space Telescope images of nearby galaxies, we simulate strong lens systems with real distributions of stellar mass as they would be observed at redshift $z_d=0.5$. We add a dark matter halo and external shear to account for the smooth dark matter field, omitting dark substructure, and use a Monte Carlo procedure to characterize the distributions of image positions, time delays, and flux ratios for a compact background source of diameter 5 pc. By convolving high-resolution images of real galaxies with a Gaussian PSF, we simulate the most detailed smooth potential one could construct given high quality data, and find scatter in flux ratios of $\approx 10\%$, which we interpret as a typical deviation from a smooth potential caused by large and small scale structure in the lensing galaxy. We demonstrate that the flux ratio anomalies arising from galaxy-scale baryonic structure can be minimized by selecting the most massive and round deflectors, and by simultaneously modeling flux ratio and astrometric data. △ Less

Submitted 26 October, 2016; originally announced October 2016.

Comments: 26 pages, 22 figures, submitted to MNRAS

Quantifying Environmental and Line-of-Sight Effects in Models of Strong Gravitational Lens Systems

Authors: Curtis McCully, Charles R. Keeton, Kenneth C. Wong, Ann I. Zabludoff

Abstract: Matter near a gravitational lens galaxy or projected along the line of sight (LOS) can affect strong lensing observables by more than contemporary measurement errors. We simulate lens fields with realistic three-dimensional mass configurations (self-consistently including voids), and then fit mock lensing observables with increasingly complex lens models to quantify biases and uncertainties associ… ▽ More Matter near a gravitational lens galaxy or projected along the line of sight (LOS) can affect strong lensing observables by more than contemporary measurement errors. We simulate lens fields with realistic three-dimensional mass configurations (self-consistently including voids), and then fit mock lensing observables with increasingly complex lens models to quantify biases and uncertainties associated with different ways of treating the lens environment (ENV) and LOS. We identify the combination of mass, projected offset, and redshift that determines the importance of a perturbing galaxy for lensing. Foreground structures have a stronger effect on the lens potential than background structures, due to nonlinear effects in the foreground and downweighting in the background. There is dramatic variation in the net strength of ENV/LOS effects across different lens fields; modeling fields individually yields stronger priors for $H_0$ than ray tracing through $N$-body simulations. Models that ignore mass outside the lens yield poor fits and biased results. Adding external shear can account for tidal stretching from galaxies at redshifts $z \ge z_{\rm lens}$, but it requires corrections for external convergence and cannot reproduce nonlinear effects from foreground galaxies. Using the tidal approximation is reasonable for most perturbers as long as nonlinear redshift effects are included. Even then, the scatter in $H_0$ is limited by the lens profile degeneracy. Asymmetric image configurations produced by highly elliptical lens galaxies are less sensitive to the lens profile degeneracy, so they offer appealing targets for precision lensing analyses in future surveys like LSST and Euclid. △ Less

Submitted 10 April, 2017; v1 submitted 20 January, 2016; originally announced January 2016.

Comments: 26 pages, 19 figures

Journal ref: ApJ, 836, 141 (2017)

A Dark Census: Statistically Detecting the Satellite Populations of Distant Galaxies

Authors: Francis-Yan Cyr-Racine, Leonidas A. Moustakas, Charles R. Keeton, Kris Sigurdson, Daniel A. Gilman

Abstract: In the standard structure formation scenario based on the cold dark matter paradigm, galactic halos are predicted to contain a large population of dark matter subhalos. While the most massive members of the subhalo population can appear as luminous satellites and be detected in optical surveys, establishing the existence of the low mass and mostly dark subhalos has proven to be a daunting task. Ga… ▽ More In the standard structure formation scenario based on the cold dark matter paradigm, galactic halos are predicted to contain a large population of dark matter subhalos. While the most massive members of the subhalo population can appear as luminous satellites and be detected in optical surveys, establishing the existence of the low mass and mostly dark subhalos has proven to be a daunting task. Galaxy-scale strong gravitational lenses have been successfully used to study mass substructures lying close to lensed images of bright background sources. However, in typical galaxy-scale lenses, the strong lensing region only covers a small projected area of the lens's dark matter halo, implying that the vast majority of subhalos cannot be directly detected in lensing observations. In this paper, we point out that this large population of dark satellites can collectively affect gravitational lensing observables, hence possibly allowing their statistical detection. Focusing on the region of the galactic halo outside the strong lensing area, we compute from first principles the statistical properties of perturbations to the gravitational time delay and position of lensed images in the presence of a mass substructure population. We find that in the standard cosmological scenario, the statistics of these lensing observables are well approximated by Gaussian distributions. The formalism developed as part of this calculation is very general and can be applied to any halo geometry and choice of subhalo mass function. Our results significantly reduce the computational cost of including a large substructure population in lens models and enable the use of Bayesian inference techniques to detect and characterize the distributed satellite population of distant lens galaxies. △ Less

Submitted 29 July, 2016; v1 submitted 4 June, 2015; originally announced June 2015.

Comments: 21 pages + appendices, 7 figures. v2: Some derivations streamlined, extended appendices. Matches version published in PRD

Journal ref: Phys. Rev. D 94, 043505 (2016)

A Spectroscopic Survey of the Fields of 28 Strong Gravitational Lenses: The Redshift Catalog

Authors: Ivelina Momcheva, Kurtis A. Williams, Richard J. Cool, Charles R. Keeton, Ann I. Zabludoff

Abstract: We present the spectroscopic redshift catalog from a wide-field survey of the fields of 28 galaxy-mass strong gravitational lenses. We discuss the acquisition and reduction of the survey data, collected over 40 nights of 6.5m MMT and Magellan time, employing four different multi-object spectrographs. We determine that no biases are introduced by combining datasets obtained with different instrumen… ▽ More We present the spectroscopic redshift catalog from a wide-field survey of the fields of 28 galaxy-mass strong gravitational lenses. We discuss the acquisition and reduction of the survey data, collected over 40 nights of 6.5m MMT and Magellan time, employing four different multi-object spectrographs. We determine that no biases are introduced by combining datasets obtained with different instrument/spectrograph combinations. Special care is taken to determine redshift uncertainties using repeat observations. The redshift catalog consists of 9768 new and unique galaxy redshifts. 82.4% of the catalog redshifts are between z=0.1 and z=0.7, and the catalog median redshift is z=0.36. The data from this survey will be used to study the lens environments and line-of-sight structures to gain a better understanding of the effects of large scale structure on lens statistics and lens-derived parameters. △ Less

Submitted 6 March, 2015; originally announced March 2015.

Comments: Accepted for publication in ApJS. All catalogs are available for download at this URL: http://www.astro.yale.edu/iva/#downloads

Characterizing the Best Cosmic Telescopes with the Millennium Simulations

Authors: K. Decker French, Kenneth C. Wong, Ann I. Zabludoff, S. Mark Ammons, Charles R. Keeton, Raul E. Angulo

Abstract: Certain configurations of massive structures projected along the line of sight maximize the number of detections of gravitationally lensed $z\sim10$ galaxies. We characterize such lines of sight with the étendue $σ_μ$, the area in the source plane magnified over some threshold $μ$. We use the Millennium I and Millennium XXL cosmological simulations to determine the frequency of high $σ_μ$ beams on… ▽ More Certain configurations of massive structures projected along the line of sight maximize the number of detections of gravitationally lensed $z\sim10$ galaxies. We characterize such lines of sight with the étendue $σ_μ$, the area in the source plane magnified over some threshold $μ$. We use the Millennium I and Millennium XXL cosmological simulations to determine the frequency of high $σ_μ$ beams on the sky, their properties, and efficient selection criteria. We define the best beams as having $σ_{μ>3} >2000$ arcsec$^2$, for a $z\sim10$ source plane, and predict $477 \pm 21$ such beams on the sky. The total mass in the beam and $σ_{μ>3}$ are strongly correlated. After controlling for total mass, we find a significant residual correlation between $σ_{μ>3}$ and the number of cluster-scale halos ($>10^{14} M_\odot h^{-1}$) in the beam. Beams with $σ_{μ>3} >2000$ arcsec$^2$, which should be best at lensing $z\sim10$ galaxies, are ten times more likely to contain multiple cluster-scale halos than a single cluster-scale halo. Beams containing an Abell 1689-like massive cluster halo often have additional structures along the line of sight, including at least one additional cluster-scale ($M_{200}>10^{14}M_\odot h^{-1}$) halo 28% of the time. Selecting beams with multiple, massive structures will lead to enhanced detection of the most distant and intrinsically faint galaxies. △ Less

Submitted 25 February, 2014; originally announced February 2014.

Comments: 14 pages, 13 figures, Accepted for publication in ApJ

A New Hybrid Framework to Efficiently Model Lines of Sight to Gravitational Lenses

Authors: Curtis McCully, Charles R. Keeton, Kenneth C. Wong, Ann I. Zabludoff

Abstract: In strong gravitational lens systems, the light bending is usually dominated by one main galaxy, but may be affected by other mass along the line of sight (LOS). Shear and convergence can be used to approximate the contributions from less significant perturbers (e.g. those that are projected far from the lens or have a small mass), but higher order effects need to be included for objects that are… ▽ More In strong gravitational lens systems, the light bending is usually dominated by one main galaxy, but may be affected by other mass along the line of sight (LOS). Shear and convergence can be used to approximate the contributions from less significant perturbers (e.g. those that are projected far from the lens or have a small mass), but higher order effects need to be included for objects that are closer or more massive. We develop a framework for multiplane lensing that can handle an arbitrary combination of tidal planes treated with shear and convergence and planes treated exactly (i.e., including higher order terms). This framework addresses all of the traditional lensing observables including image positions, fluxes, and time delays to facilitate lens modelling that includes the non-linear effects due to mass along the LOS. It balances accuracy (accounting for higher-order terms when necessary) with efficiency (compressing all other LOS effects into a set of matrices that can be calculated up front and cached for lens modelling). We identify a generalized multiplane mass sheet degeneracy, in which the effective shear and convergence are sums over the lensing planes with specific, redshift-dependent weighting factors. △ Less

Submitted 5 September, 2014; v1 submitted 31 December, 2013; originally announced January 2014.

Comments: 13 pages, 2 figures

Journal ref: MNRAS, 443, 3631 (2014)

Mapping Compound Cosmic Telescopes Containing Multiple, Projected Cluster-Scale Halos

Authors: S. Mark Ammons, Kenneth C. Wong, Ann I. Zabludoff, Charles R. Keeton

Abstract: Lines of sight with multiple, projected, cluster-scale halos have high total masses and complex lens plane interactions that can boost the area of magnification, or etendue, making detection of faint background sources more likely than elsewhere. To identify these new "compound" cosmic telescopes, we have found lines-of-sight with the highest integrated mass densities, as traced by the projected c… ▽ More Lines of sight with multiple, projected, cluster-scale halos have high total masses and complex lens plane interactions that can boost the area of magnification, or etendue, making detection of faint background sources more likely than elsewhere. To identify these new "compound" cosmic telescopes, we have found lines-of-sight with the highest integrated mass densities, as traced by the projected concentrations of Luminous Red Galaxies (LRGs). We use 1151 MMT Hectospec spectra to derive preliminary magnification maps for two such lines of sight with total mass exceeding ~ 3 x 10$^{15}$ Msun -- J0850+3604 (0850) and J1306+4632 (1306). We identify 2-3 group- and cluster-scale halos in each beam over 0.1 < z < 0.7, all of which are well-traced by LRGs. In Subaru Suprime-Cam imaging of beam 0850, we discover serendipitously a candidate multiply-imaged V-dropout source at z = 5.03, whose location is consistent with the critical curves for a source plane of $z_s$ = 5.03 predicted by our mass model. Incorporating the position of the candidate multiply-imaged galaxy as a constraint on the critical curve location in 0850 narrows the 68% confidence band on lens plane area with mu > 10 for a source plane of $z_s$ = 10 to [1.8, 4.2] square arcminutes, comparable to that of MACS 0717+3745 and El Gordo, two of the most powerful known single cluster lenses. The 68% confidence intervals on the lens plane area with mu > 10 for 1306 are [2.3, 6.7] square arcminutes. The significant lensing power of our beams makes them powerful probes of reionization and galaxy formation in the early Universe. △ Less

Submitted 11 November, 2013; originally announced November 2013.

Comments: Accepted to ApJ, 23 pages, 15 figures. PDF version with full-scale figures available at http://www.u.arizona.edu/~ammons81/ms.pdf until publication

Report number: LLNL release number: LLNL-JRNL-609412

Dark energy with gravitational lens time delays

Authors: T. Treu, P. J. Marshall, F. -Y. Cyr-Racine, C. D. Fassnacht, C. R. Keeton, E. V. Linder, L. A. Moustakas, M. Bradac, E. Buckley-Geer, T. Collett, F. Courbin, G. Dobler, D. A. Finley, J. Hjorth, C. S. Kochanek, E. Komatsu, L. V. E. Koopmans, G. Meylan, P. Natarajan, M. Oguri, S. H. Suyu, M. Tewes, K. C. Wong, A. I. Zabludoff, D. Zaritsky , et al. (13 additional authors not shown)

Abstract: Strong lensing gravitational time delays are a powerful and cost effective probe of dark energy. Recent studies have shown that a single lens can provide a distance measurement with 6-7 % accuracy (including random and systematic uncertainties), provided sufficient data are available to determine the time delay and reconstruct the gravitational potential of the deflector. Gravitational-time delays… ▽ More Strong lensing gravitational time delays are a powerful and cost effective probe of dark energy. Recent studies have shown that a single lens can provide a distance measurement with 6-7 % accuracy (including random and systematic uncertainties), provided sufficient data are available to determine the time delay and reconstruct the gravitational potential of the deflector. Gravitational-time delays are a low redshift (z~0-2) probe and thus allow one to break degeneracies in the interpretation of data from higher-redshift probes like the cosmic microwave background in terms of the dark energy equation of state. Current studies are limited by the size of the sample of known lensed quasars, but this situation is about to change. Even in this decade, wide field imaging surveys are likely to discover thousands of lensed quasars, enabling the targeted study of ~100 of these systems and resulting in substantial gains in the dark energy figure of merit. In the next decade, a further order of magnitude improvement will be possible with the 10000 systems expected to be detected and measured with LSST and Euclid. To fully exploit these gains, we identify three priorities. First, support for the development of software required for the analysis of the data. Second, in this decade, small robotic telescopes (1-4m in diameter) dedicated to monitoring of lensed quasars will transform the field by delivering accurate time delays for ~100 systems. Third, in the 2020's, LSST will deliver 1000's of time delays; the bottleneck will instead be the aquisition and analysis of high resolution imaging follow-up. Thus, the top priority for the next decade is to support fast high resolution imaging capabilities, such as those enabled by the James Webb Space Telescope and next generation adaptive optics systems on large ground based telescopes. △ Less

Submitted 5 June, 2013; originally announced June 2013.

Comments: White paper submitted to SNOWMASS2013

A New Approach to Identifying the Most Powerful Gravitational Lensing Telescopes

Authors: Kenneth C. Wong, Ann I. Zabludoff, S. Mark Ammons, Charles R. Keeton, David W. Hogg, Anthony H. Gonzalez

Abstract: The best gravitational lenses for detecting distant galaxies are those with the largest mass concentrations and the most advantageous configurations of that mass along the line of sight. Our new method for finding such gravitational telescopes uses optical data to identify projected concentrations of luminous red galaxies (LRGs). LRGs are biased tracers of the underlying mass distribution, so line… ▽ More The best gravitational lenses for detecting distant galaxies are those with the largest mass concentrations and the most advantageous configurations of that mass along the line of sight. Our new method for finding such gravitational telescopes uses optical data to identify projected concentrations of luminous red galaxies (LRGs). LRGs are biased tracers of the underlying mass distribution, so lines of sight with the highest total luminosity in LRGs are likely to contain the largest total mass. We apply this selection technique to the Sloan Digital Sky Survey and identify the 200 fields with the highest total LRG luminosities projected within a 3.5' radius over the redshift range 0.1 < z < 0.7. The redshift and angular distributions of LRGs in these fields trace the concentrations of non-LRG galaxies. These fields are diverse; 22.5% contain one known galaxy cluster and 56.0% contain multiple known clusters previously identified in the literature. Thus, our results confirm that these LRGs trace massive structures and that our selection technique identifies fields with large total masses. These fields contain 2-3 times higher total LRG luminosities than most known strong-lensing clusters and will be among the best gravitational lensing fields for the purpose of detecting the highest redshift galaxies. △ Less

Submitted 19 April, 2013; v1 submitted 10 April, 2013; originally announced April 2013.

Comments: Accepted for publication in The Astrophysical Journal; emulateapj format; 35 pages, 8 figures, 3 tables

The Discovery of the Most Distant Known Type Ia Supernova at Redshift 1.914

Authors: David O. Jones, Steven A. Rodney, Adam G. Riess, Bahram Mobasher, Tomas Dahlen, Curtis McCully, Teddy F. Frederiksen, Stefano Casertano, Jens Hjorth, Charles R. Keeton, Anton Koekemoer, Louis-Gregory Strolger, Tommy G. Wiklind, Peter Challis, Or Graur, Brian Hayden, Brandon Patel, Benjamin J. Weiner, Alexei V. Filippenko, Peter Garnavich, Saurabh W. Jha, Robert P. Kirshner, Henry C. Ferguson, Norman A. Grogin, Dale Kocevski

Abstract: We present the discovery of a Type Ia supernova (SN) at redshift $z = 1.914$ from the CANDELS multi-cycle treasury program on the \textit{Hubble Space Telescope (HST)}. This SN was discovered in the infrared using the Wide-Field Camera 3, and it is the highest-redshift Type Ia SN yet observed. We classify this object as a SN\,Ia by comparing its light curve and spectrum with those of a large sampl… ▽ More We present the discovery of a Type Ia supernova (SN) at redshift $z = 1.914$ from the CANDELS multi-cycle treasury program on the \textit{Hubble Space Telescope (HST)}. This SN was discovered in the infrared using the Wide-Field Camera 3, and it is the highest-redshift Type Ia SN yet observed. We classify this object as a SN\,Ia by comparing its light curve and spectrum with those of a large sample of Type Ia and core-collapse supernovae (SNe). Its apparent magnitude is consistent with that expected from the $Λ$CDM concordance cosmology. We discuss the use of spectral evidence for classification of $z > 1.5$ SNe\,Ia using {\it HST} grism simulations, finding that spectral data alone can frequently rule out SNe\,II, but distinguishing between SNe\,Ia and SNe\,Ib/c can require prohibitively long exposures. In such cases, a quantitative analysis of the light curve may be necessary for classification. Our photometric and spectroscopic classification methods can aid the determination of SN rates and cosmological parameters from the full high-redshift CANDELS SN sample. △ Less

Submitted 15 April, 2013; v1 submitted 2 April, 2013; originally announced April 2013.

Comments: 10 pages, 5 figures; accepted for publication in ApJ; minor changes

Optimal Mass Configurations for Lensing High-Redshift Galaxies

Authors: Kenneth C. Wong, S. Mark Ammons, Charles R. Keeton, Ann I. Zabludoff

Abstract: We investigate the gravitational lensing properties of lines of sight containing multiple cluster-scale halos, motivated by their ability to lens very high-redshift (z ~ 10) sources into detectability. We control for the total mass along the line of sight, isolating the effects of distributing the mass among multiple halos and of varying the physical properties of the halos. Our results show that… ▽ More We investigate the gravitational lensing properties of lines of sight containing multiple cluster-scale halos, motivated by their ability to lens very high-redshift (z ~ 10) sources into detectability. We control for the total mass along the line of sight, isolating the effects of distributing the mass among multiple halos and of varying the physical properties of the halos. Our results show that multiple-halo lines of sight can increase the magnified source-plane region compared to the single cluster lenses typically targeted for lensing studies, and thus are generally better fields for detecting very high-redshift sources. The configurations that result in optimal lensing cross sections benefit from interactions between the lens potentials of the halos when they overlap somewhat on the sky, creating regions of high magnification in the source plane not present when the halos are considered individually. The effect of these interactions on the lensing cross section can even be comparable to changing the total mass of the lens from 10^15 M_sun to 3x10^15 M_sun. The gain in lensing cross section increases as the mass is split into more halos, provided that the lens potentials are projected close enough to interact with each other. A nonzero projected halo angular separation, equal halo mass ratio, and high projected halo concentration are the best mass configurations, whereas projected halo ellipticity, halo triaxiality, and the relative orientations of the halos are less important. Such high mass, multiple-halo lines of sight exist in the SDSS. △ Less

Submitted 4 March, 2013; v1 submitted 12 March, 2012; originally announced March 2012.

Comments: Accepted for publication in ApJ; emulateapj format; 24 pages, 13 figures, 1 table; plots updated to reflect erratum

Substructure in the lens HE 0435-1223

Authors: Ross Fadely, Charles R. Keeton

Abstract: We investigate the properties of dark matter substructure in the gravitational lens HE 0435-1223 (z_l=0.455) via its effects on the positions and flux ratios of the quadruply-imaged background quasar (z_s=1.689). We start with a smooth mass model, add individual, truncated isothermal clumps near the lensed images, and use the Bayesian evidence to compare the quality of different models. Compared w… ▽ More We investigate the properties of dark matter substructure in the gravitational lens HE 0435-1223 (z_l=0.455) via its effects on the positions and flux ratios of the quadruply-imaged background quasar (z_s=1.689). We start with a smooth mass model, add individual, truncated isothermal clumps near the lensed images, and use the Bayesian evidence to compare the quality of different models. Compared with smooth models, models with at least one clump near image A are strongly favored. The mass of this clump within its Einstein radius is log(Mein/Msun) = 7.65 +0.87/-0.84. The Bayesian evidence provides weaker support for a second clump near image B, with log(Mein/Msun) = 6.55 +1.01/-1.51. We also examine models with a full population of substructure, and find the mass fraction in substructure at the Einstein radius to be f_sub > 0.00077, assuming the total clump masses follow a mass function dN/dM proportional to M^(-1.9) over the range M = 10^7-10^10 Msun. Few-clump and population models produce similar Bayesian evidence values, so neither type of model is objectively favored. △ Less

Submitted 2 September, 2011; originally announced September 2011.

Comments: 18 pages, 12 figures, 4 tables. MNRAS accepted

Lensing by Kerr Black Holes. II: Analytical Study of Quasi-Equatorial Lensing Observables

Authors: Amir B. Aazami, Charles R. Keeton, A. O. Petters

Abstract: In this second paper, we develop an analytical theory of quasi-equatorial lensing by Kerr black holes. In this setting we solve perturbatively our general lens equation with displacement given in Paper I, going beyond weak-deflection Kerr lensing to third order in our expansion parameter epsilon, which is the ratio of the angular gravitational radius to the angular Einstein radius. We obtain new f… ▽ More In this second paper, we develop an analytical theory of quasi-equatorial lensing by Kerr black holes. In this setting we solve perturbatively our general lens equation with displacement given in Paper I, going beyond weak-deflection Kerr lensing to third order in our expansion parameter epsilon, which is the ratio of the angular gravitational radius to the angular Einstein radius. We obtain new formulas and results for the bending angle, image positions, image magnifications, total unsigned magnification, and centroid, all to third order in epsilon and including the displacement. New results on the time delay between images are also given to second order in epsilon, again including displacement. For all lensing observables we show that the displacement begins to appear only at second order in epsilon. When there is no spin, we obtain new results on the lensing observables for Schwarzschild lensing with displacement. △ Less

Submitted 12 October, 2011; v1 submitted 21 February, 2011; originally announced February 2011.

Comments: 23 pages; final published version

Journal ref: J. Math. Phys. 52, 102501 (2011)

Lensing by Kerr Black Holes. I: General Lens Equation and Magnification Formula

Authors: Amir B. Aazami, Charles R. Keeton, A. O. Petters

Abstract: We develop a unified, analytic framework for gravitational lensing by Kerr black holes. In this first paper we present a new, general lens equation and magnification formula governing lensing by a compact object. Our lens equation assumes that the source and observer are in the asymptotically flat region and does not require a small angle approximation. Furthermore, it takes into account the displ… ▽ More We develop a unified, analytic framework for gravitational lensing by Kerr black holes. In this first paper we present a new, general lens equation and magnification formula governing lensing by a compact object. Our lens equation assumes that the source and observer are in the asymptotically flat region and does not require a small angle approximation. Furthermore, it takes into account the displacement that occurs when the light ray's tangent lines at the source and observer do not meet on the lens plane. We then explore our lens equation in the case when the compact object is a Kerr black hole. Specifically, we give an explicit expression for the displacement when the observer is in the equatorial plane of the Kerr black hole as well as for the case of spherical symmetry. △ Less

Submitted 12 October, 2011; v1 submitted 21 February, 2011; originally announced February 2011.

Comments: 11 pages; final published version

Journal ref: J. Math. Phys. 52, 092502 (2011)

On statistical uncertainty in nested sampling

Authors: Charles R. Keeton

Abstract: Nested sampling has emerged as a valuable tool for Bayesian analysis, in particular for determining the Bayesian evidence. The method is based on a specific type of random sampling of the likelihood function and prior volume of the parameter space. I study the statistical uncertainty in the evidence computed with nested sampling. I examine the uncertainty estimator from Skilling (2004, 2006) and i… ▽ More Nested sampling has emerged as a valuable tool for Bayesian analysis, in particular for determining the Bayesian evidence. The method is based on a specific type of random sampling of the likelihood function and prior volume of the parameter space. I study the statistical uncertainty in the evidence computed with nested sampling. I examine the uncertainty estimator from Skilling (2004, 2006) and introduce a new estimator based on a detailed analysis of the statistical properties of nested sampling. Both perform well in test cases and make it possible to obtain the statistical uncertainty in the evidence with no additional computational cost. △ Less

Submitted 4 February, 2011; originally announced February 2011.

Comments: 9 pages, 4 figures; accepted to MNRAS

Near-Infrared K and L' Flux Ratios in Six Lensed Quasars

Authors: Ross Fadely, Charles R. Keeton

Abstract: We examine the wavelength dependence of flux ratios for six gravitationally lensed quasars using K and L' images obtained at the Gemini North 8m telescope. We select lenses with source redshifts z_s < 2.8 so that K-band images probe rest-frame optical emission from accretion disks, while L'-band images probe rest-frame near-infrared flux emitted (in part) from the more extended surrounding torus.… ▽ More We examine the wavelength dependence of flux ratios for six gravitationally lensed quasars using K and L' images obtained at the Gemini North 8m telescope. We select lenses with source redshifts z_s < 2.8 so that K-band images probe rest-frame optical emission from accretion disks, while L'-band images probe rest-frame near-infrared flux emitted (in part) from the more extended surrounding torus. Since the observations correspond to different source sizes, the K and L' flux ratios are sensitive to structure on different scales and may be useful for studying small-structure in the lens galaxies. Four of the six lenses show differences between K and L' flux ratios. In HE 0435$-1223, SDSS 0246-0825, and HE 2149-2745 the differences may be attributable to known microlensing and/or intrinsic variability. In SDSS 0806+2006 the wavelength dependence is not easily attributed to known variations, and may indicate the presence of substructure. By contrast, in Q0142-100 and SBS 0909+523 the K and L' flux ratios are consistent within the uncertainties. We discuss the utility of the current data for studying chromatic effects related to microlensing, dust extinction, and dark matter substructure. △ Less

Submitted 10 January, 2011; originally announced January 2011.

Comments: 19 pages, 4 Figures, 4 Tables, AJ accepted

The Effect of Environment on Shear in Strong Gravitational Lenses

Authors: Kenneth C. Wong, Charles R. Keeton, Kurtis A. Williams, Ivelina G. Momcheva, Ann I. Zabludoff

Abstract: Using new photometric and spectroscopic data in the fields of nine strong gravitational lenses that lie in galaxy groups, we analyze the effects of both the local group environment and line-of-sight galaxies on the lens potential. We use Monte Carlo simulations to derive the shear directly from measurements of the complex lens environment, providing the first detailed independent check of the shea… ▽ More Using new photometric and spectroscopic data in the fields of nine strong gravitational lenses that lie in galaxy groups, we analyze the effects of both the local group environment and line-of-sight galaxies on the lens potential. We use Monte Carlo simulations to derive the shear directly from measurements of the complex lens environment, providing the first detailed independent check of the shear obtained from lens modeling. We account for possible tidal stripping of the group galaxies by varying the fraction of total mass apportioned between the group dark matter halo and individual group galaxies. The environment produces an average shear of gamma = 0.08 (ranging from 0.02 to 0.17), significant enough to affect quantities derived from lens observables. However, the direction and magnitude of the shears do not match those obtained from lens modeling in three of the six 4-image systems in our sample (B1422, RXJ1131, and WFI2033). The source of this disagreement is not clear, implying that the assumptions inherent in both the environment and lens model approaches must be reconsidered. If only the local group environment of the lens is included, the average shear is gamma = 0.05 (ranging from 0.01 to 0.14), indicating that line-of-sight contributions to the lens potential are not negligible. We isolate the effects of various theoretical and observational uncertainties on our results. Of those uncertainties, the scatter in the Faber-Jackson relation and error in the group centroid position dominate. Future surveys of lens environments should prioritize spectroscopic sampling of both the local lens environment and objects along the line of sight, particularly those bright (I < 21.5) galaxies projected within 5' of the lens. △ Less

Submitted 10 November, 2010; originally announced November 2010.

Comments: Accepted for publication in The Astrophysical Journal; 28 pages, 9 figures, 5 tables

LSST Science Book, Version 2.0

Authors: LSST Science Collaboration, Paul A. Abell, Julius Allison, Scott F. Anderson, John R. Andrew, J. Roger P. Angel, Lee Armus, David Arnett, S. J. Asztalos, Tim S. Axelrod, Stephen Bailey, D. R. Ballantyne, Justin R. Bankert, Wayne A. Barkhouse, Jeffrey D. Barr, L. Felipe Barrientos, Aaron J. Barth, James G. Bartlett, Andrew C. Becker, Jacek Becla, Timothy C. Beers, Joseph P. Bernstein, Rahul Biswas, Michael R. Blanton, Joshua S. Bloom , et al. (223 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 1 December, 2009; originally announced December 2009.

Comments: 596 pages. Also available at full resolution at http://www.lsst.org/lsst/scibook

Identifying Anomalies in Gravitational Lens Time Delays

Authors: Arthur B. Congdon, Charles R. Keeton, C. Erik Nordgren

Abstract: We examine the ability of gravitational lens time delays to reveal complex structure in lens potentials. In Congdon, Keeton & Nordgren (2008), we predicted how the time delay between the bright pair of images in a "fold" lens scales with the image separation, for smooth lens potentials. Here we show that the proportionality constant increases with the quadrupole moment of the lens potential, and… ▽ More We examine the ability of gravitational lens time delays to reveal complex structure in lens potentials. In Congdon, Keeton & Nordgren (2008), we predicted how the time delay between the bright pair of images in a "fold" lens scales with the image separation, for smooth lens potentials. Here we show that the proportionality constant increases with the quadrupole moment of the lens potential, and depends only weakly on the position of the source along the caustic. We use Monte Carlo simulations to determine the range of time delays that can be produced by realistic smooth lens models consisting of isothermal ellipsoid galaxies with tidal shear. We can then identify outliers as "time delay anomalies". We find evidence for anomalies in close image pairs in the cusp lenses RX J1131$-$1231 and B1422+231. The anomalies in RX J1131$-$1231 provide strong evidence for substructure in the lens potential, while at this point the apparent anomalies in B1422+231 mainly indicate that the time delay measurements need to be improved. We also find evidence for time delay anomalies in larger-separation image pairs in the fold lenses, B1608+656 and WFI 2033$-$4723, and the cusp lens RX J0911+0551. We suggest that these anomalies are caused by some combination of substructure and a complex lens environment. Finally, to assist future monitoring campaigns we use our smooth models with shear to predict the time delays for all known four-image lenses. △ Less

Submitted 21 December, 2009; v1 submitted 12 November, 2009; originally announced November 2009.

Comments: 24 pages, 10 figures, 5 tables, accepted by ApJ, minor edits to match text of published version

Journal ref: Astrophys.J.709:552-571,2010

Improved Constraints on the Gravitational Lens Q0957+561. II. Strong Lensing

Authors: Ross Fadely, Charles R. Keeton, Reiko Nakajima, Gary M. Bernstein

Abstract: We present a detailed strong lensing analysis of an HST/ACS legacy dataset for the first gravitational lens, Q0957+561. With deep imaging we identify 24 new strongly lensed features, which we use to constrain mass models. We model the stellar component of the lens galaxy using the observed luminosity distribution, and the dark matter halo using several different density profiles. We draw on the… ▽ More We present a detailed strong lensing analysis of an HST/ACS legacy dataset for the first gravitational lens, Q0957+561. With deep imaging we identify 24 new strongly lensed features, which we use to constrain mass models. We model the stellar component of the lens galaxy using the observed luminosity distribution, and the dark matter halo using several different density profiles. We draw on the weak lensing analysis by Nakajima et al. (2009) to constrain the mass sheet and environmental terms in the lens potential. Adopting the well-measured time delay, we find H_0 = 85 (+14/-13) km/s/Mpc (68% CL) using lensing constraints alone. The principal uncertainties in H_0 are tied to the stellar mass-to-light ratio (a variant of the radial profile degeneracy in lens models). Adding constraints from stellar population synthesis models, we obtain H_0 = 79.3 (+6.7/-8.5) km/s/Mpc (68% CL). We infer that the lens galaxy has a rising rotation curve and a dark matter distribution with an inner core. Intriguingly, we find the quasar flux ratios predicted by our models to be inconsistent with existing radio measurements, suggesting the presence of substructure in the lens. △ Less

Submitted 8 January, 2010; v1 submitted 9 September, 2009; originally announced September 2009.

Comments: 23 pages, 19 figures, 7 tables; ApJ accepted. Some additional comments and minor changes from v1

Journal ref: Astrophys.J.711:246-267,2010

Gravitational lensing with stochastic substructure: Effects of the clump mass function and spatial distribution

Authors: Charles R. Keeton

Abstract: Mass clumps in gravitational lens galaxies can perturb lensed images in characteristic ways. Strong lens flux ratios have been used to constrain the amount of dark matter substructure in lens galaxies, and various other observables have been considered as additional probes of substructure. We study the general theory of lensing with stochastic substructure in order to understand how lensing obse… ▽ More Mass clumps in gravitational lens galaxies can perturb lensed images in characteristic ways. Strong lens flux ratios have been used to constrain the amount of dark matter substructure in lens galaxies, and various other observables have been considered as additional probes of substructure. We study the general theory of lensing with stochastic substructure in order to understand how lensing observables depend on the mass function and spatial distribution of clumps. We find that magnification perturbations are mainly sensitive to the total mass in substructure projected near the lensed images; when the source is small, flux ratios are not very sensitive to the shape of the clump mass function. Position perturbations are mainly sensitive to a characteristic clump mass scale, namely m_eff = <m^2>/<m>, with some mild dependence on other mass moments when the spatial distribution is not uniform. They have contributions from both "local" and "global" populations of clumps (i.e., those projected near the images, and those farther away). Time delay perturbations are sensitive to the same characteristic mass, m_eff, and mainly driven by the global population of clumps. While there is significant scatter in all lensing quantities, there are some non-trivial correlations that may contain further information about the clump population. Our results indicate that a joint analysis of multiple lens observables will offer qualitatively new constraints on the mass function and spatial distribution of dark matter substructure in distant galaxies. △ Less

Submitted 20 August, 2009; originally announced August 2009.

Comments: 18 pages, 7 figures, submitted to MNRAS

Gravitational Lensing of Anisotropic Sources

Authors: Rosalba Perna, Charles R. Keeton

Abstract: In strong gravitational lensing, the multiple images we see correspond to light rays that leave the source in slightly different directions. If the source emission is anisotropic, the images may differ from conventional lensing predictions (which assume isotropy). To identify scales on which source anisotropy may be important, we study the angle delta between the light rays emerging from the sou… ▽ More In strong gravitational lensing, the multiple images we see correspond to light rays that leave the source in slightly different directions. If the source emission is anisotropic, the images may differ from conventional lensing predictions (which assume isotropy). To identify scales on which source anisotropy may be important, we study the angle delta between the light rays emerging from the source, for different lensing configurations. If the lens has a power law profile M propto R^gamma, the angle delta initially increases with lens redshift and then either diverges (for a steep profile gamma<1), remains constant (for an isothermal profile gamma=1), or vanishes (for a shallow profile gamma>1) as zl approaches zs. The scaling with lens mass is roughly delta propto M^(1/(2-gamma)). The results for an NFW profile are qualitatively similar to those for a shallow power law, with delta peaking at about half the redshift of the source (not half the distance). In practice, beaming could modify the statistics of beamed sources lensed by massive clusters: for an opening angle theta_jet, there is a probability as high as P ~ 0.02-0.07 (theta_jet/{0.5} deg)^-1 that one of the lensed images may be missed (for 2 < zs < 6). Differential absorption within Active Galactic Nuclei could modify the flux ratios of AGNs lensed by clusters; a sample of AGNs lensed by clusters could provide further constraints on the sizes of absorbing regions. Source anisotropy is not likely to be a significant effect in galaxy-scale strong lensing. △ Less

Submitted 14 May, 2009; v1 submitted 24 April, 2009; originally announced April 2009.

Comments: minor changes, accepted to MNRAS

Planetesimal Disk Microlensing

Authors: Kevin Heng, Charles R. Keeton

Abstract: Motivated by debris disk studies, we investigate the gravitational microlensing of background starlight by a planetesimal disk around a foreground star. We use dynamical survival models to construct a plausible example of a planetesimal disk and study its microlensing properties using established ideas of microlensing by small bodies. When a solar-type source star passes behind a planetesimal di… ▽ More Motivated by debris disk studies, we investigate the gravitational microlensing of background starlight by a planetesimal disk around a foreground star. We use dynamical survival models to construct a plausible example of a planetesimal disk and study its microlensing properties using established ideas of microlensing by small bodies. When a solar-type source star passes behind a planetesimal disk, the microlensing light curve may exhibit short-term, low-amplitude residuals caused by planetesimals several orders of magnitude below Earth mass. The minimum planetesimal mass probed depends on the photometric sensitivity and the size of the source star, and is lower when the planetesimal lens is located closer to us. Planetesimal lenses may be found more nearby than stellar lenses because the steepness of the planetesimal mass distribution changes how the microlensing signal depends on the lens/source distance ratio. Microlensing searches for planetesimals require essentially continuous monitoring programs that are already feasible and can potentially set constraints on models of debris disks, the supposed extrasolar analogues of Kuiper belts. △ Less

Submitted 19 October, 2009; v1 submitted 30 March, 2009; originally announced March 2009.

Comments: 15 pages, 8 figures (emulateapj). Accepted by ApJ

Journal ref: Astrophys.J.707:621-631,2009

Improved Constraints on the Gravitational Lens Q0957+561. I. Weak Lensing

Authors: R. Nakajima, G. M. Bernstein, R. Fadely, C. R. Keeton, T. Schrabback

Abstract: Attempts to constrain the Hubble constant using the strong gravitational lens system Q0957+561 are limited by systematic uncertainties in the mass model, since the time delay is known very precisely. One important systematic effect is the mass sheet degeneracy, which arises because strong lens modeling cannot constrain the presence or absence of a uniform mass sheet $κ$, which rescales $H_0$ by… ▽ More Attempts to constrain the Hubble constant using the strong gravitational lens system Q0957+561 are limited by systematic uncertainties in the mass model, since the time delay is known very precisely. One important systematic effect is the mass sheet degeneracy, which arises because strong lens modeling cannot constrain the presence or absence of a uniform mass sheet $κ$, which rescales $H_0$ by the factor $(1-κ)$. In this paper we present new constraints on the mass sheet derived from a weak-lensing analysis of the Hubble Space Telescope imaging of a 6 arcmin square region surrounding the lensed quasar. The average mass sheet within a circular aperture (the strong lens model region) is constrained by integrating the tangential weak gravitational shear over the surrounding area. We find the average convergence within a $30"$ radius around the lens galaxy to be $κ(<30") = 0.166\pm0.056$ ($1σ$ confidence level), normalized to the quasar redshift. This includes contributions from both the lens galaxy and the surrounding cluster. We also constrain a few other low-order terms in the lens potential by applying a multipole aperture mass formalism to the gravitational shear in an annulus around the strong lensing region. Implications for strong lens models and the Hubble constant are discussed in an accompanying paper. △ Less

Submitted 8 April, 2009; v1 submitted 24 March, 2009; originally announced March 2009.

Comments: 13 pages, 8 figures, accepted for publication in ApJ

Journal ref: Astrophys.J.697:1793-1804,2009

Strong gravitational lensing probes of the particle nature of dark matter

Authors: Leonidas A. Moustakas, Kevork Abazajian, Andrew Benson, Adam S. Bolton, James S. Bullock, Jacqueline Chen, Edward Cheng, Dan Coe, Arthur B. Congdon, Neal Dalal, Juerg Diemand, Benjamin M. Dobke, Greg Dobler, Olivier Dore, Aaron Dutton, Richard Ellis, Chris D. Fassnacht, Henry Ferguson, Douglas Finkbeiner, Raphael Gavazzi, Fredrick William High, Tesla Jeltema, Eric Jullo, Manoj Kaplinghat, Charles R. Keeton , et al. (22 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 18 February, 2009; originally announced February 2009.

Comments: Science white paper submitted to the Astro2010 Decadal Cosmology & Fundamental Physics Science Frontier Panel

Dark Matter Structures in the Universe: Prospects for Optical Astronomy in the Next Decade

Authors: P. J. Marshall, M. Auger, J. G. Bartlett, M. Bradac, A. Cooray, N. Dalal, G. Dobler, C. D. Fassnacht, B. Jain, C. R. Keeton, R. Mandelbaum, L. A. Moustakas, M. A. Strauss, J. A. Tyson, D. Wittman, S. A. Wright

Abstract: The Cold Dark Matter theory of gravitationally-driven hierarchical structure formation has earned its status as a paradigm by explaining the distribution of matter over large spans of cosmic distance and time. However, its central tenet, that most of the matter in the universe is dark and exotic, is still unproven; the dark matter hypothesis is sufficiently audacious as to continue to warrant a… ▽ More The Cold Dark Matter theory of gravitationally-driven hierarchical structure formation has earned its status as a paradigm by explaining the distribution of matter over large spans of cosmic distance and time. However, its central tenet, that most of the matter in the universe is dark and exotic, is still unproven; the dark matter hypothesis is sufficiently audacious as to continue to warrant a diverse battery of tests. While local searches for dark matter particles or their annihilation signals could prove the existence of the substance itself, studies of cosmological dark matter in situ are vital to fully understand its role in structure formation and evolution. We argue that gravitational lensing provides the cleanest and farthest-reaching probe of dark matter in the universe, which can be combined with other observational techniques to answer the most challenging and exciting questions that will drive the subject in the next decade: What is the distribution of mass on sub-galactic scales? How do galaxy disks form and bulges grow in dark matter halos? How accurate are CDM predictions of halo structure? Can we distinguish between a need for a new substance (dark matter) and a need for new physics (departures from General Relativity)? What is the dark matter made of anyway? We propose that the central tool in this program should be a wide-field optical imaging survey, whose true value is realized with support in the form of high-resolution, cadenced optical/infra-red imaging, and massive-throughput optical spectroscopy. △ Less

Submitted 17 February, 2009; originally announced February 2009.

Comments: White paper submitted to the 2010 Astronomy & Astrophysics Decadal Survey

Galaxy density profiles and shapes -- II. selection biases in strong lensing surveys

Authors: Rachel Mandelbaum, Glenn van de Ven, Charles R. Keeton

Abstract: [Abridged] Many current and future astronomical surveys will rely on samples of strong gravitational lens systems to draw conclusions about galaxy mass distributions. We use a new strong lensing pipeline (presented in Paper I of this series) to explore selection biases that may cause the population of strong lensing systems to differ from the general galaxy population. Our focus is on point-sour… ▽ More [Abridged] Many current and future astronomical surveys will rely on samples of strong gravitational lens systems to draw conclusions about galaxy mass distributions. We use a new strong lensing pipeline (presented in Paper I of this series) to explore selection biases that may cause the population of strong lensing systems to differ from the general galaxy population. Our focus is on point-source lensing by early-type galaxies with two mass components (stellar and dark matter) that have a variety of density profiles and shapes motivated by observational and theoretical studies of galaxy properties. We seek not only to quantify but also to understand the physics behind selection biases related to: galaxy mass, orientation and shape; dark matter profile parameters such as inner slope and concentration; and adiabatic contraction. We study how all of these properties affect the lensing Einstein radius, total cross-section, quad/double ratio, and image separation distribution. We find significant (factors of several) selection biases with mass; orientation, for a given galaxy shape at fixed mass; cusped dark matter profile inner slope and concentration; concentration of the stellar and dark matter deprojected Sersic models. Interestingly, the intrinsic shape of a galaxy does not strongly influence its lensing cross-section when we average over viewing angles. Our results are an important first step towards understanding how strong lens systems relate to the general galaxy population. △ Less

Submitted 18 June, 2009; v1 submitted 19 August, 2008; originally announced August 2008.

Comments: 26 pages, 15 figures; paper I at arXiv:0808.2493; accepted for publication in MNRAS (minor revisions); PDF file with full resolution figures at http://www.sns.ias.edu/~rmandelb/paper2.pdf

Galaxy density profiles and shapes -- I. simulation pipeline for lensing by realistic galaxy models

Authors: Glenn van de Ven, Rachel Mandelbaum, Charles R. Keeton

Abstract: Studies of strong gravitational lensing in current and upcoming wide and deep photometric surveys, and of stellar kinematics from (integral-field) spectroscopy at increasing redshifts, promise to provide valuable constraints on galaxy density profiles and shapes. However, both methods are affected by various selection and modelling biases, whch we aim to investigate in a consistent way. In this… ▽ More Studies of strong gravitational lensing in current and upcoming wide and deep photometric surveys, and of stellar kinematics from (integral-field) spectroscopy at increasing redshifts, promise to provide valuable constraints on galaxy density profiles and shapes. However, both methods are affected by various selection and modelling biases, whch we aim to investigate in a consistent way. In this first paper in a series we develop a flexible but efficient pipeline to simulate lensing by realistic galaxy models. These galaxy models have separate stellar and dark matter components, each with a range of density profiles and shapes representative of early-type, central galaxies without significant contributions from other nearby galaxies. We use Fourier methods to calculate the lensing properties of galaxies with arbitrary surface density distributions, and Monte Carlo methods to compute lensing statistics such as point-source lensing cross-sections. Incorporating a variety of magnification bias modes lets us examine different survey limitations in image resolution and flux. We rigorously test the numerical methods for systematic errors and sensitivity to basic assumptions. We also determine the minimum number of viewing angles that must be sampled in order to recover accurate orientation-averaged lensing quantities. We find that for a range of non-isothermal stellar and dark matter density profiles typical of elliptical galaxies, the combined density profile and corresponding lensing properties are surprisingly close to isothermal around the Einstein radius. The converse implication is that constraints from strong lensing and/or stellar kinematics, which are indeed consistent with isothermal models near the Einstein radius, cannot trivially be extrapolated to smaller and larger radii. △ Less

Submitted 18 June, 2009; v1 submitted 19 August, 2008; originally announced August 2008.

Comments: 31 pages, 15 figures; paper II at arXiv:0808.2497; accepted for publication in MNRAS; PDF file with full resolution figures at http://www.sns.ias.edu/~glenn/paper1.pdf