We present a Bayesian framework to account for the magnification bias from both strong and weak gravitational lensing in estimates of high-redshift galaxy luminosity functions. We illustrate our method by estimating the $zsim8$ UV luminosity function
using a sample of 97 Y-band dropouts (Lyman break galaxies) found in the Brightest of Reionizing Galaxies (BoRG) survey and from the literature. We find the luminosity function is well described by a Schechter function with characteristic magnitude of $M^star = -19.85^{+0.30}_{-0.35}$, faint-end slope of $alpha = -1.72^{+0.30}_{-0.29}$, and number density of $log_{10} Psi^star [textrm{Mpc}^{-3}] = -3.00^{+0.23}_{-0.31}$. These parameters are consistent within the uncertainties with those inferred from the same sample without accounting for the magnification bias, demonstrating that the effect is small for current surveys at $zsim8$, and cannot account for the apparent overdensity of bright galaxies compared to a Schechter function found recently by Bowler et al. (2014a,b) and Finkelstein et al. (2014). We estimate that the probability of finding a strongly lensed $zsim8$ source in our sample is in the range $sim 3-15 %$ depending on limiting magnitude. We identify one strongly-lensed candidate and three cases of intermediate lensing in BoRG (estimated magnification $mu>1.4$) in addition to the previously known candidate group-scale strong lens. Using a range of theoretical luminosity functions we conclude that magnification bias will dominate wide field surveys -- such as those planned for the Euclid and WFIRST missions -- especially at $z>10$. Magnification bias will need to be accounted for in order to derive accurate estimates of high-redshift luminosity functions in these surveys and to distinguish between galaxy formation models.
Gravitational lensing - the deflection of light rays by gravitating matter - has become a major tool in the armoury of the modern cosmologist. Proposed nearly a hundred years ago as a key feature of Einsteins theory of General Relativity, we trace th
e historical development since its verification at a solar eclipse in 1919. Einstein was apparently cautious about its practical utility and the subject lay dormant observationally for nearly 60 years. Nonetheless there has been rapid progress over the past twenty years. The technique allows astronomers to chart the distribution of dark matter on large and small scales thereby testing predictions of the standard cosmological model which assumes dark matter comprises a massive weakly-interacting particle. By measuring distances and tracing the growth of dark matter structure over cosmic time, gravitational lensing also holds great promise in determining whether the dark energy, postulated to explain the accelerated cosmic expansion, is a vacuum energy density or a failure of General Relativity on large scales. We illustrate the wide range of applications which harness the power of gravitational lensing, from searches for the earliest galaxies magnified by massive clusters to those for extrasolar planets which temporarily brighten a background star. We summarise the future prospects with dedicated ground and space-based facilities designed to exploit this remarkable physical phenomenon.
In 1964, Refsdal hypothesized that a supernova whose light traversed multiple paths around a strong gravitational lens could be used to measure the rate of cosmic expansion. We report the discovery of such a system. In Hubble Space Telescope imaging,
we have found four images of a single supernova forming an Einstein cross configuration around a redshift z=0.54 elliptical galaxy in the MACS J1149.6+2223 cluster. The clusters gravitational potential also creates multiple images of the z=1.49 spiral supernova host galaxy, and a future appearance of the supernova elsewhere in the cluster field is expected. The magnifications and staggered arrivals of the supernova images probe the cosmic expansion rate, as well as the distribution of matter in the galaxy and cluster lenses.
We investigate the cosmic evolution of the internal structure of massive early-type galaxies over half of the age of the Universe. We perform a joint lensing and stellar dynamics analysis of a sample of 81 strong lenses from the SL2S and SLACS survey
s and combine the results with a hierarchical Bayesian inference method to measure the distribution of dark matter mass and stellar IMF across the population of massive early-type galaxies. Lensing selection effects are taken into account. We find that the dark matter mass projected within the inner 5 kpc increases for increasing redshift, decreases for increasing stellar mass density, but is roughly constant along the evolutionary tracks of early-type galaxies. The average dark matter slope is consistent with that of an NFW profile, but is not well constrained. The stellar IMF normalization is close to a Salpeter IMF at $log{M_*} = 11.5$ and scales strongly with increasing stellar mass. No dependence of the IMF on redshift or stellar mass density is detected. The anti-correlation between dark matter mass and stellar mass density supports the idea of mergers being more frequent in more massive dark matter halos.
We create a baseline of the black hole (BH) mass (MBH) - stellar-velocity dispersion (sigma) relation for active galaxies, using a sample of 66 local (0.02<z<0.09) Seyfert-1 galaxies, selected from the Sloan Digital Sky Survey (SDSS). Analysis of SDS
S images yields AGN luminosities free of host-galaxy contamination and morphological classification. 51/66 galaxies have spiral morphology. 28 bulges have Sersic index n<2 and are considered candidate pseudo bulges, with eight being definite pseudo bulges based on multiple classification criteria met. Only 4/66 galaxies show sign of interaction/merging. High signal-to-noise ratio Keck spectra provide the width of the broad Hbeta emission line free of FeII emission and stellar absorption. AGN luminosity and Hbeta line widths are used to estimate MBH. The Keck-based spatially-resolved kinematics is used to determine stellar-velocity dispersion within the spheroid effective radius. We find that sigma can vary on average by up to 40% across definitions commonly used in the literature, emphasizing the importance of using self-consistent definitions in comparisons and evolutionary studies. The MBH-sigma relation for our Seyfert-1 galaxies has the same intercept and scatter as that of reverberation-mapped AGNs as well as quiescent galaxies, consistent with the hypothesis that our single epoch MBH estimator and sample selection do not introduce significant biases. Barred galaxies, merging galaxies, and those hosting pseudo bulges do not represent outliers in the MBH-sigma relation. This is in contrast with previous work, although no firm conclusion can be drawn due to the small sample size and limited resolution of the SDSS images.
We present an investigation about the shape of the initial mass function (IMF) of early-type galaxies (ETGs), based on a joint lensing and dynamical analysis, and on stellar population synthesis models, for a sample of 55 lens ETGs identified by the
Sloan Lens ACS (SLACS) Survey. We construct axisymmetric dynamical models based on the Jeans equations which allow for orbital anisotropy and include a dark matter halo. The models reproduce in detail the observed textit{HST} photometry and are constrained by the total projected mass within the Einstein radius and the stellar velocity dispersion ($sigma$) within the SDSS fibers. Comparing the dynamically-derived stellar mass-to-light ratios $(M_*/L)_{rm dyn}$, obtained for an assumed halo slope $rho_{rm h}propto r^{-1}$, to the stellar population ones $(M_*/L)_{rm pop}$, derived from full-spectrum fitting and assuming a Salpeter IMF, we infer the mass normalization of the IMF. Our results confirm the previous analysis by the SLACS team that the mass normalization of the IMF of high $sigma$ galaxies is consistent on average with a Salpeter slope. Our study allows for a fully consistent study of the trend between IMF and $sigma$ for both the SLACS and ATLAS samples, which explore quite different $sigma$ ranges. The two samples are highly complementary, the first being essentially $sigma$ selected, and the latter volume-limited and nearly mass selected. We find that the two samples merge smoothly into a single trend of the form $logalpha =(0.38pm0.04)timeslog(sigma_{rm e}/200,mathrm{km~s}^{-1})+(-0.06pm0.01)$, where $alpha=(M_*/L)_{rm dyn}/(M_*/L)_{rm pop}$ and $sigma_{rm e}$ is the luminosity averaged $sigma$ within one effective radius $R_{rm e}$. This is consistent with a systematic variation of the IMF normalization from Kroupa to Salpeter in the interval $sigma_{rm e}approx90-270,mathrm{km~s}^{-1}$.
We present an improved and expanded simply parameterized phenomenological model of the broad line region (BLR) in active galactic nuclei (AGN) for modeling reverberation mapping data. By modeling reverberation mapping data directly, we can constrain
the geometry and dynamics of the BLR and measure the black hole mass without relying on the normalization factor needed in the traditional analysis. For realistic simulated reverberation mapping datasets of high-quality, we can recover the black hole mass to $0.05-0.25$ dex uncertainty and distinguish between dynamics dominated by elliptical orbits and inflowing gas. While direct modeling of the integrated emission line light curve allows for measurement of the mean time lag, other details of the geometry of the BLR are better constrained by the full spectroscopic dataset of emission line profiles. We use this improved model of the BLR to explore possible sources of uncertainty in measurements of the time lag using cross-correlation function (CCF) analysis and in measurements of the black hole mass using the virial product. Sampling the range of geometries and dynamics in our model of the BLR suggests that the theoretical uncertainty in black hole masses measured using the virial product is on the order of 0.25 dex. These results support the use of the CCF to measure time lags and the virial product to measure black hole masses when direct modeling techniques cannot be applied, provided the uncertainties associated with the interpretation of the results are taken into account.
We present dynamical modeling of the broad line region (BLR) for a sample of five Seyfert 1 galaxies using reverberation mapping data taken by the Lick AGN Monitoring Project in 2008. By modeling the AGN continuum light curve and H$beta$ line profile
s directly we are able to constrain the geometry and kinematics of the BLR and make a measurement of the black hole mass that does not depend upon the virial factor, $f$, needed in traditional reverberation mapping analysis. We find that the geometry of the BLR is generally a thick disk viewed close to face-on. While the H$beta$ emission is found to come preferentially from the far side of the BLR, the mean size of the BLR is consistent with the lags measured with cross-correlation analysis. The BLR kinematics are found to be consistent with either inflowing motions or elliptical orbits, often with some combination of the two. We measure black hole masses of $log_{10}(M_{rm,BH}/M_odot)=6.62^{+0.10}_{-0.13}$ for Arp 151, $7.42^{+0.26}_{-0.27}$ for Mrk 1310, $7.51^{+0.23}_{-0.14}$ for NGC 5548, $6.42^{+0.24}_{-0.18}$ for NGC 6814, and $6.99^{+0.32}_{-0.25}$ for SBS 1116+583A. The $f$ factors measured individually for each AGN are found to correlate with inclination angle, although not with $M_{rm,BH}$, $L_{5100}$, or FWHM/$sigma$ of the emission line profile.
We present high-quality Keck/LRIS longslit spectroscopy of a pilot sample of 25 local active galaxies selected from the SDSS (0.02<z<0.1; MBH>10^7 M_sun) to study the relations between black hole mass (MBH) and host-galaxy properties. We determine st
ellar kinematics of the host galaxy, deriving stellar-velocity dispersion profiles and rotation curves from three spectral regions (including CaH&K, MgIb triplet, and CaII triplet). In addition, we perform surface photometry on SDSS images, using a newly developed code for joint multi-band analysis. BH masses are estimated from the width of the Hbeta emission line and the host-galaxy free 5100A AGN luminosity. Combining results from spectroscopy and imaging allows us to study four MBH scaling relations: MBH-sigma, MBH-L(sph), MBH-M(sph,*), MBH-M(sph,dyn). We find the following results. First, stellar-velocity dispersions determined from aperture spectra (e.g. SDSS fiber spectra or unresolved data from distant galaxies) can be biased, depending on aperture size, AGN contamination, and host-galaxy morphology. However, such a bias cannot explain the offset seen in the MBH-sigma relation at higher redshifts. Second, while the CaT region is the cleanest to determine stellar-velocity dispersions, both the MgIb region, corrected for FeII emission, and the CaHK region, although often swamped by the AGN powerlaw continuum and emission lines, can give results accurate to within a few percent. Third, the MBH scaling relations of our pilot sample agree in slope and scatter with those of other local active and inactive galaxies. In the next papers of the series we will quantify the scaling relations, exploiting the full sample of ~100 objects.
(Abridged) We complete the census of nuclear X-ray activity in 100 early type Virgo galaxies observed by the Chandra X-ray Telescope as part of the AMUSE-Virgo survey, down to a (3sigma) limiting luminosity of 3.7E+38 erg/s over 0.5-7 keV. The stella
r mass distribution of the targeted sample, which is mostly composed of formally `inactive galaxies, peaks below 1E+10 M_Sun, a regime where the very existence of nuclear super-massive black holes (SMBHs) is debated. Out of 100 objects, 32 show a nuclear X-ray source, including 6 hybrid nuclei which also host a massive nuclear cluster as visible from archival HST images. After carefully accounting for contamination from nuclear low-mass X-ray binaries based on the shape and normalization of their X-ray luminosity function, we conclude that between 24-34% of the galaxies in our sample host a X-ray active SMBH (at the 95% C.L.). This sets a firm lower limit to the black hole occupation fraction in nearby bulges within a cluster environment. At face value, the active fraction -down to our luminosity limit- is found to increase with host stellar mass. However, taking into account selection effects, we find that the average Eddington-scaled X-ray luminosity scales with black hole mass as M_BH^(-0.62^{+0.13}_{-0.12}), with an intrinsic scatter of 0.46^({+0.08}_{-0.06}) dex. This finding can be interpreted as observational evidence for `down-sizing of black hole accretion in local early types, that is, low mass black holes shine relatively closer to their Eddington limit than higher mass objects. As a consequence, the fraction of active galaxies, defined as those above a fixed X-ray Eddington ratio, decreases with increasing black hole mass.
