Gravitationally lensed quasars are useful for studying astrophysics and cosmology, and enlarging the sample size of lensed quasars is important for multiple studies. In this work, we develop a lens search algorithm for four-image (quad) lensed quasars based on their time variability. In the development of the lens search algorithm, we constructed a pipeline simulating multi-epoch images of lensed quasars in cadenced surveys, accounting for quasar variabilities, quasar hosts, lens galaxies, and the PSF variation. Applying the simulation pipeline to the Hyper Suprime-Cam (HSC) transient survey, we generated HSC-like difference images of the mock lensed quasars from Oguri & Marshall's lens catalog. We further developed a lens search algorithm that picks out variable objects as lensed quasar candidates based on their spatial extent in the difference images. We tested our lens search algorithm with the mock lensed quasars and variable objects from the HSC transient survey. Using difference images from multiple epochs, our lens search algorithm achieves a high true-positive rate (TPR) of 90.1% and a low false-positive rate (FPR) of 2.3% for the bright quads with wide separation. With a preselection of the number of blobs in the difference image, we obtain a TPR of 97.6% and a FPR of 2.6% for the bright quads with wide separation. Even when difference images are only available in one single epoch, our lens search algorithm can still detect the bright quads with wide separation at high TPR of 97.6% and low FPR of 2.4% in the optimal seeing scenario, and at TPR of $\sim94%$ and FPR of $\sim5%$ in typical scenarios. Therefore, our lens search algorithm is promising and is applicable to ongoing and upcoming cadenced surveys, particularly the HSC transient survey and the Rubin Observatory Legacy Survey of Space and Time, for finding new lensed quasar systems. [abridged]

Time-delay cosmography of lensed quasars has achieved 2.4% precision on the measurement of the Hubble constant, $H_0$. As part of an ongoing effort to uncover and control systematic uncertainties, we investigate three potential sources: 1- stellar kinematics, 2- line-of-sight effects, and 3- the deflector mass model. To meet this goal in a quantitative way, we reproduced the H0LiCOW/SHARP/STRIDES (hereafter TDCOSMO) procedures on a set of real and simulated data, and we find the following. First, stellar kinematics cannot be a dominant source of error or bias since we find that a systematic change of 10% of measured velocity dispersion leads to only a 0.7% shift on $H_0$ from the seven lenses analyzed by TDCOSMO. Second, we find no bias to arise from incorrect estimation of the line-of-sight effects. Third, we show that elliptical composite (stars + dark matter halo), power-law, and cored power-law mass profiles have the flexibility to yield a broad range in $H_0$ values. However, the TDCOSMO procedures that model the data with both composite and power-law mass profiles are informative. If the models agree, as we observe in real systems owing to the "bulge-halo" conspiracy, $H_0$ is recovered precisely and accurately by both models. If the two models disagree, as in the case of some pathological models illustrated here, the TDCOSMO procedure either discriminates between them through the goodness of fit, or it accounts for the discrepancy in the final error bars provided by the analysis. This conclusion is consistent with a reanalysis of six of the TDCOSMO (real) lenses: the composite model yields $H_0=74.2^{+1.6}_{-1.6}$ $km.s^{-1}.Mpc^{-1}$, while the power-law model yields $74.0^{+1.7}_{-1.8}$ $km.s^{-1}.Mpc^{-1}$. In conclusion, we find no evidence of bias or errors larger than the current statistical uncertainties reported by TDCOSMO.

We present spectroscopic confirmation of three new two-image gravitationally lensed quasars, compiled from existing strong lens and X-ray catalogs. Images of HSC J091843.27-022007.5 show a red galaxy with two blue point sources at either side, separated by 2.26 arcsec. This system has a source and a lens redshifts $z_s = 0.804$ and $z_l = 0.459$, respectively, as obtained by our follow-up spectroscopic data. CXCO J100201.50+020330.0 shows two point sources separated by 0.85 arcsec on either side of an early-type galaxy. The follow-up spectroscopic data confirm the fainter quasar has the same redshift with the brighter quasar from the SDSS fiber spectrum at $z_s = 2.016$. The deflecting foreground galaxy is a typical early-type galaxy at a redshift of $z_l = 0.439$. SDSS J135944.21+012809.8 has two point sources with quasar spectra at the same redshift $z_s = 1.096$, separated by 1.05 arcsec, and fits to the HSC images confirm the presence of a galaxy between these. These discoveries demonstrate the power of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP)'s deep imaging and wide sky coverage. Combined with existing X-ray source catalogues and follow-up spectroscopy, the HSC-SSP provides us unique opportunities to find multiple-image quasars lensed by a foreground galaxy.

We study the dynamical properties of massive quiescent galaxies at $1.4 < z < 2.1$ using deep Hubble Space Telescope WFC3/F160W imaging and a combination of literature stellar velocity dispersion measurements and new near-infrared spectra obtained using KMOS on the ESO VLT. We use these data to show that the typical dynamical-to-stellar mass ratio has increased by $\sim$0.2 dex from $z = 2$ to the present day, and investigate this evolution in the context of possible changes in the stellar initial mass function (IMF) and/or fraction of dark matter contained within the galaxy effective radius, $f_\mathrm{DM}$. Comparing our high-redshift sample to their likely descendants at low-redshift, we find that $f_\mathrm{DM}$ has increased by a factor of more than 4 since $z \approx 1.8$, from $f_\mathrm{DM}$ = $6.6\pm1.0$% to $\sim$24%. The observed increase appears robust to changes in the methods used to estimate dynamical masses or match progenitors and descendants. We quantify possible variation of the stellar IMF through the offset parameter $\alpha$, defined as the ratio of dynamical mass in stars to the stellar mass estimated using a Chabrier IMF. We demonstrate that the correlation between stellar velocity dispersion and $\alpha$ reported among quiescent galaxies at low-redshift is already in place at $z = 2$, and argue that subsequent evolution through (mostly minor) merging should act to preserve this relation while contributing significantly to galaxies overall growth in size and stellar mass.

We present six new time-delay measurements obtained from $R_c$-band monitoring data acquired at the MPIA 2.2m telescope at La Silla observatory between October 2016 and February 2020. The lensed quasars HE 0047-1756, WG 0214-2105, DES 0407-5006, 2M 1134-2103, PSJ 1606-2333 and DES 2325-5229 were observed almost daily at high signal-to-noise ratio to obtain high-quality light curves where we can record fast and small-amplitude variations of the quasars. We measure time delays between all pairs of multiple images with only one or two seasons of monitoring at the exception of the time delays relative to image D of PSJ 1606-2333. The most precise estimate is obtained for the delay between image A and image B of DES 0407-5006, with $\tau_{AB} = -128.4^{+3.5}_{-3.8}$ days (2.8% precision) including systematics due to extrinsic variability in the light curves. For HE 0047-1756, we combine our high-cadence data with measurements from decade-long light-curves from previous COSMOGRAIL campaigns, and reach a precision of 0.9 day on the final measurement. The present work demonstrates the feasibility of measuring time delays in lensed quasars in only one or two seasons, provided high signal-to-noise ratio data are obtained at a cadence close to daily.

We study the stellar mass functions (SMFs) of star-forming and quiescent galaxies in 11 galaxy clusters at 1.0<z<1.4, drawn from the Gemini Observations of Galaxies in Rich Early Environments (GOGREEN) survey. Based on more than 500 hours of Gemini/GMOS spectroscopy, and deep multi-band photometry taken with a range of observatories, we probe the SMFs down to a stellar mass limit of 10^9.7 Msun (10^9.5 Msun for star-forming galaxies). At this early epoch, the fraction of quiescent galaxies is already highly elevated in the clusters compared to the field at the same redshift. The quenched fraction excess (QFE) represents the fraction of galaxies that would be star-forming in the field, but are quenched due to their environment. The QFE is strongly mass dependent, and increases from ~30% at Mstar=10^9.7 Msun, to ~80% at Mstar=10^11.0 Msun. Nonetheless, the shapes of the SMFs of the two individual galaxy types, star-forming and quiescent galaxies, are identical between the clusters and the field - to high statistical precision. Yet, along with the different quiescent fractions is the total galaxy SMF environmentally dependent, with a relative deficit of low-mass galaxies in the clusters. These results are in stark contrast with findings in the local Universe, and thus require a substantially different quenching mode to operate at early times. We discuss these results in the light of several popular quenching models.

We report the largest sample of candidate strong gravitational lenses belonging to the Survey of Gravitationally-lensed Objects in HSC Imaging for group-to-cluster scale (SuGOHI-c) systems. These candidates are compiled from the S18A data release of the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) Survey. We visually inspect $\sim39,500$ galaxy clusters, selected from several catalogs, overlapping with the Wide, Deep, and UltraDeep fields, spanning the cluster redshift range $0.05<z_{cl}<1.38$. We discover 641 candidate lens systems, of which 536 are new. From the full sample, 47 are almost certainly bonafide lenses, 181 of them are highly probable lenses and 413 are possible lens systems. Additionally, we present 131 lens candidates at galaxy-scale serendipitously discovered during the inspection. We obtained spectroscopic follow-up of 10 candidates using the X-shooter. With this follow-up, we confirm 8 systems as strong gravitational lenses. Of the remaining two, one of the sources is too faint to detect any emission, and the other has a tentative redshift close to the lens redshift, but additional arcs in this system are yet to be observed spectroscopically. Since the HSC-SSP is an ongoing survey, we expect to find $\sim600$ definite or probable lenses using this procedure and even more if combined with other lens finding methods.

We present a blind time-delay cosmographic analysis for the lens system DES J0408$-$5354. This system is extraordinary for the presence of two sets of multiple images at different redshifts, which provide the opportunity to obtain more information at the cost of increased modelling complexity with respect to previously analyzed systems. We perform detailed modelling of the mass distribution for this lens system using three band Hubble Space Telescope imaging. We combine the measured time delays, line-of-sight central velocity dispersion of the deflector, and statistically constrained external convergence with our lens models to estimate two cosmological distances. We measure the "effective" time-delay distance corresponding to the redshifts of the deflector and the lensed quasar $D_{\Delta t}^{\rm eff}=3382^{+146}_{-115}$ Mpc and the angular diameter distance to the deflector $D_{\rm d}=1711^{+376}_{-280}$ Mpc, with covariance between the two distances. From these constraints on the cosmological distances, we infer the Hubble constant $H_0 = 74.2^{+2.7}_{-3.0}$ km s$^{-1}$ Mpc$^{-1}$ assuming a flat $\Lambda$CDM cosmology and a uniform prior for $\Omega_{\rm m}$ as $\Omega_{\rm m} \sim \mathcal{U}(0.05, 0.5)$. This measurement gives the most precise constraint on $H_0$ to date from a single lens. Our measurement is consistent with that obtained from the previous sample of six lenses analyzed by the $H_0$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) collaboration. It is also consistent with measurements of $H_0$ based on the local distance ladder, reinforcing the tension with the inference from early Universe probes, for example, with 2.2$\sigma$ discrepancy from the cosmic microwave background measurement.

Strong lenses are extremely useful probes of the distribution of matter on galaxy and cluster scales at cosmological distances, but are rare and difficult to find. The number of currently known lenses is on the order of 1,000. We wish to use crowdsourcing to carry out a lens search targeting massive galaxies selected from over 442 square degrees of photometric data from the Hyper Suprime-Cam (HSC) survey. We selected a sample of $\sim300,000$ galaxies with photometric redshifts in the range $0.2 < z_{phot} < 1.2$ and photometrically inferred stellar masses $\log{M_*} > 11.2$. We crowdsourced lens finding on this sample of galaxies on the Zooniverse platform, as part of the Space Warps project. The sample was complemented by a large set of simulated lenses and visually selected non-lenses, for training purposes. Nearly 6,000 citizen volunteers participated in the experiment. In parallel, we used YattaLens, an automated lens finding algorithm, to look for lenses in the same sample of galaxies. Based on a statistical analysis of classification data from the volunteers, we selected a sample of the most promising $\sim1,500$ candidates which we then visually inspected: half of them turned out to be possible (grade C) lenses or better. Including lenses found by YattaLens or serendipitously noticed in the discussion section of the Space Warps website, we were able to find 14 definite lenses, 129 probable lenses and 581 possible lenses. YattaLens found half the number of lenses discovered via crowdsourcing. Crowdsourcing is able to produce samples of lens candidates with high completeness and purity, compared to currently available automated algorithms. A hybrid approach, in which the visual inspection of samples of lens candidates pre-selected by discovery algorithms and/or coupled to machine learning is crowdsourced, will be a viable option for lens finding in the 2020s.

We present results on the environmental dependence of the star-forming galaxy main sequence in 11 galaxy cluster fields at $1.0 < z < 1.5$ from the Gemini Observations of Galaxies in Rich Early Environments Survey (GOGREEN) survey. We use a homogeneously selected sample of field and cluster galaxies whose membership is derived from dynamical analysis. Using [OII]-derived star formation rates (SFRs), we find that cluster galaxies have suppressed SFRs at fixed stellar mass in comparison to their field counterparts by a factor of 1.4 $\pm$ 0.1 ($\sim3.3\sigma$) across the stellar mass range: $9.0 < \log(M_{*} /M_{\odot}) < 11.2$. We also find that this modest suppression in the cluster galaxy star-forming main sequence is mass and redshift dependent: the difference between cluster and field increases towards lower stellar masses and lower redshift. When comparing the distribution of cluster and field galaxy SFRs to the star-forming main sequence, we find an overall shift towards lower SFRs in the cluster population, and note the absence of a tail of high SFR galaxies as seen in the field. Given this observed suppression in the cluster galaxy star-forming main sequence, we explore the implications for several scenarios such as formation time differences between cluster and field galaxies, and environmentally-induced star formation quenching and associated timescales.

Minor mergers have been proposed as the driving mechanism for the size growth of quiescent galaxies with decreasing redshift. The process whereby large star-forming galaxies quench and join the quiescent population at the large size end has also been suggested as an explanation for this size growth. Given the clear association of quenching with clusters, we explore this mechanism by studying the structural properties of 23 spectroscopically identified recently quenched (or "poststarburst" (PSB)) cluster galaxies at $z\sim1$. Despite clear PSB spectral signatures implying rapid and violent quenching, 87\% of these galaxies have symmetric, undisturbed morphologies in the stellar continuum. Remarkably, they follow a mass-size relation lying midway between the star-forming and quiescent field relations, with sizes $0.1$ dex smaller than $z\sim1$ star-forming galaxies at log$(M_{*}/M_{\odot})=10.5$. This implies a rapid change in the light profile without directly effecting the stellar distribution, suggesting changes in the mass-to-light ratio gradients across the galaxy are responsible. We develop fading toy models to explore how star-forming galaxies move across the mass-size plane as their stellar populations fade to match those of the PSBs. "Outside-in" fading has the potential to reproduce the contraction in size and increase in bulge-dominance observed between star-forming and PSB cluster galaxies. Since cluster PSBs lie on the large size end of the quiescent mass-size relation, and our previous work shows cluster galaxies are smaller than field galaxies, the sizes of quiescent galaxies must grow both from the quenching of star-forming galaxies and dry minor mergers.

We report the serendipitous discovery of HSC J0904$-$0102, a quadruply-lensed Lyman break galaxy (LBG) in the Survey of Gravitationally-lensed Objects in Hyper Suprime-Cam Imaging (SuGOHI). Owing to its point-like appearance, the source was thought to be a lensed active galactic nucleus. We obtained follow-up spectroscopic data with the Gemini Multi-Object Spectrographs on the Gemini South Telescope, which confirmed this to be a lens system. The deflecting foreground galaxy is a typical early-type galaxy at a high redshift of $z_{\ell} = 0.957$ with stellar velocity dispersion $\sigma_v=259\pm56$ km~s$^{-1}$. The lensed source is identified as an LBG at $z_{\rm s} = 3.403$, based on the sharp drop bluewards of Ly$\alpha$ and other absorption features. A simple lens mass model for the system, assuming a singular isothermal ellipsoid, yields an Einstein radius of $\theta_{\rm Ein} = 1. 23^{\prime\prime}$ and a total mass within the Einstein radius of $M_{\rm Ein} = (5.55\pm 0.24) \times 10^{11}M_{\odot}$ corresponding to a velocity dispersion of $\sigma_{\rm SIE}= 283\pm 3$ km~s$^{-1}$, which is in good agreement with the value derived spectroscopically. The most isolated lensed LBG image has a magnification of $\sim 6.5$. In comparison with other lensed LBGs and typical $z\sim4$ LBG populations, HSC J0904$-$0102 is unusually compact, an outlier at $>2\sigma$ confidence. Together with a previously discovered SuGOHI lens, HSC J1152$+$0047, that is similarly compact, we believe that the HSC Survey is extending LBG studies down to smaller galaxy sizes.

Extremely light bosonic wave dark matter ($\psi$DM) is an emerging dark matter candidate contesting the conventional cold dark matter paradigm and a model subject to intense scrutiny of late. This work for the first time reports testable salient features pertinent to gravitational lenses of $\psi$DM halos. $\psi$DM halos are distinctly filled with large-amplitude, small-scale density fluctuations with $\delta\rho/\rho_{\rm halo}\sim 1$ in form of density granules. This halo characteristically yields fuzzy gravitational lenses that can yield ubiquitous flux ratio anomalies of few tens percent levels often found in lensed quasar multiple images. These fuzzy lenses may also produce rare hexad images and octad images, formed when the sources are located in some well-defined zones of highly distorted caustics. We have found new critical features appearing in the highly de-magnified lens center when the halo has sufficiently high surface density near a very compact massive core.

We present half-light sizes measured from H${\alpha}$ emission tracing star-formation in 281 star-forming galaxies from the KMOS3D survey at 0.7 < z < 2.7. Sizes are derived by fitting 2D exponential disk models, with bootstrap errors averaging 20%. H${\alpha}$ sizes are a median (mean) of 1.19 (1.26) times larger than those of the stellar continuum, which due to radial dust gradients places an upper limit on the growth in stellar size via star formation, with just 43% intrinsic scatter. At fixed continuum size the H${\alpha}$ size shows no residual trend with stellar mass, star formation rate, redshift or morphology. The only significant residual trend is with the excess obscuration of H${\alpha}$ by dust, at fixed continuum obscuration. The scatter in continuum size at fixed stellar mass is likely driven by the scatter in halo spin parameters. The stability of the ratio of H${\alpha}$ size to continuum size demonstrates a high degree of stability in halo spin and in the transfer of angular momentum to the disk over a wide range of physical conditions and cosmic time. This may require local regulation by feedback processes. The implication of our results, as we demonstrate using a toy model, is that our upper limit on star-formation driven growth is sufficient only to evolve star-forming galaxies approximately along the observed size-mass relation, consistent with the size growth of galaxies at constant cumulative co-moving number density. To explain the observed evolution of the size-mass relation of star-forming disk galaxies other processes, such as the preferential quenching of compact galaxies or galaxy mergers, may be required.

We present the HOLISMOKES project on strong gravitational lensing of supernovae (SNe) as a probe of SN physics and cosmology. We investigate the effects of microlensing on early-phase SN Ia spectra, and find that within 10 rest-frame days after SN explosion, distortions of SN Ia spectra due to microlensing are typically negligible (<1% distortion within the 68% credible region, and ~10% distortion within the 95% credible region). This shows great prospects of using lensed SNe Ia to obtain intrinsic early-phase SN spectra for deciphering SN Ia progenitors. As a demonstration of the usefulness of lensed SNe Ia for cosmology, we simulate a sample of mock lensed SN Ia systems that are expected to have accurate and precise time-delay measurements in the era of the Legacy Survey of Space and Time (LSST). Adopting realistic yet conservative uncertainties on their time-delay distances and lens angular diameter distances (of 6.6% and 5%, respectively), we find that a sample of 20 lensed SNe Ia would allow a constraint on the Hubble constant ($H_0$) with 1.3% uncertainty in the flat $\Lambda$CDM cosmology. We find a similar constraint on $H_0$ in an open $\Lambda$CDM cosmology, while the constraint degrades to 3% in a flat wCDM cosmology. We anticipate lensed SNe to be an independent and powerful probe of SN physics and cosmology in the upcoming LSST era.

We present the results of 15 years of monitoring of lensed quasars conducted by the COSMOGRAIL program at the Leonhard Euler 1.2m Swiss Telescope. The decade-long light curves of 23 lensed QSOs are presented for the first time. We complement our data set with other monitoring data available in the literature to measure the time delays in 18 systems, among which 9 reach a relative precision better than 15% for at least one time delay. To achieve this, we develop an automated version of the curve-shifting toolbox PyCS to ensure robust estimation of the time delay in the presence of microlensing while accounting for the errors due to imperfect representation of microlensing. We also reanalyze the previously published time delays of RX J1131$-$1231 and HE 0435$-$1223, adding respectively 6 and 2 new seasons of monitoring and confirming the previous time-delay measurements. When the time delay measurement is possible, we correct the light curves of the lensed images from their time delay and present the difference curves to highlight the microlensing signal contained in the data. This is to date the largest sample of decade-long lens monitoring data, useful to measure $H_0$, to measure the size of quasar accretion disks with microlensing, and to study quasar variability.

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

A double source plane (DSP) system is a precious probe for the density profile of distant galaxies and cosmological parameters. However, these measurements could be affected by the surrounding environment of the lens galaxy. Thus, it is important to evaluate the cluster-scale mass for detailed mass modeling. We observed the {\it Eye of Horus}, a DSP system discovered by the Subaru HSC--SSP, with XMM--Newton. We detected two X-ray extended emissions, originating from two clusters, one centered at the {\it Eye of Horus}, and the other located $\sim100$ arcsec northeast to the {\it Eye of Horus}. We determined the dynamical mass assuming hydrostatic equilibrium, and evaluated their contributions to the lens mass interior of the Einstein radius. The contribution of the former cluster is $1.1^{+1.2}_{-0.5}\times10^{12}~M_{\odot}$, which is $21-76\%$ of the total mass within the Einstein radius. The discrepancy is likely due to the complex gravitational structure along the line of sight. On the other hand, the contribution of the latter cluster is only $\sim2\%$ on the {\it Eye of Horus}. Therefore, the influence associated with this cluster can be ignored.

Strong gravitationally lensed quasars provide powerful means to study galaxy evolution and cosmology. We use Chitah to hunt for new lens systems in the Hyper Suprime$-$Cam Subaru Strategic Program (HSC SSP) S16A. We present 46 lens candidates, of which 3 are previously known. Including 2 additional lenses found by YattaLens, we obtain X-shooter spectra of 6 promising candidates for lens confirmation and redshift measurements. We report new spectroscopic redshift measurements for both the lens and source galaxies in 4 lens systems. We apply the lens modeling software Glee to model our 6 X-shooter lenses uniformly. Through our analysis of the HSC images, we find that HSCJ022622$-$042522, HSCJ115252$+$004733, and HSCJ141136$-$010216 have point-like lensed images, and that the lens light distribution is well aligned with mass distribution within 6 deg. Thanks to the X-shooter spectra, we estimate fluxes on the Baldwin- Phillips-Terlevich (BPT) diagram, and find that HSCJ022622$-$042522 has a probable quasar source, based on the upper limit of the Nii flux intensity. We also measure the FWHM of Ly$\alpha$ emission of HSCJ141136$-$010216 to be $\sim$254 km/s, showing that it is a probable Lyman-$\alpha$ emitter.

With the advent of high-cadence and multi-band photometric monitoring facilities, continuum reverberation mapping is becoming of increasing importance to measure the physical size of quasar accretion disks. The method is based on the measurement of the time it takes for a signal to propagate from the center to the outer parts of the central engine, assuming the continuum light curve at a given wavelength has a time shift of the order of a few days with respect to light curves obtained at shorter wavelengths. We show that with high-quality light curves, this assumption is not valid anymore and that light curves at different wavelengths are not only shifted in time but also distorted: in the context of the lamp-post model and thin-disk geometry, the multi-band light curves are in fact convolved by a transfer function whose size increase with wavelength. We illustrate the effect with simulated light curves in the LSST ugrizy bands and examine the impact on the delay measurements when using three different methods, namely JAVELIN, CREAM, and PyCS. We find that current accretion disk sizes estimated from JAVELIN and PyCS are underestimated by $\sim30\%$ and that unbiased measurement are only obtained with methods that properly take the skewed transfer functions into account, as the CREAM code does. With the LSST-like light curves, we expect to achieve measurement errors below $5\%$ with typical 2-day photometric cadence.

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

We present spectra of the most massive quiescent galaxy yet discovered at $z>3$, spectroscopically confirmed via the detection of Balmer absorption features in the $H-$ and $K-$bands of Keck/MOSFIRE. The spectra confirm a galaxy with no significant ongoing star formation, consistent with the lack of rest-frame UV flux and overall photometric spectral energy distribution. With a stellar mass of $3.1^{+0.1}_{-0.2} \times 10^{11} ~\rm{M}_\odot$ at $z = 3.493$, this galaxy is nearly three times more massive than the highest redshift spectroscopically confirmed absorption-line identified galaxy known. The star-formation history of this quiescent galaxy implies that it formed $>1000 ~\rm{M}_\odot$/yr for almost 0.5 Gyr beginning at $z\sim7.2$, strongly suggestive that it is the descendant of massive dusty star-forming galaxies at $5<z<7$ recently observed with ALMA. While galaxies with similarly extreme stellar masses are reproduced in some simulations at early times, such a lack of ongoing star formation is not seen there. This suggests the need for a more rapid quenching process than is currently prescribed, challenging our current understanding of how ultra-massive galaxies form and evolve in the early Universe.

The upcoming Large Synoptic Survey Telescope (LSST) will detect many strongly lensed Type Ia supernovae (LSNe Ia) for time-delay cosmography. This will provide an independent and direct way for measuring the Hubble constant $H_0$, which is necessary to address the current $4.4 \sigma$ tension in $H_0$ between the local distance ladder and the early Universe measurements. We present a detailed analysis of different observing strategies for the LSST, and quantify their impact on time-delay measurement between multiple images of LSNe Ia. For this, we produced microlensed mock-LSST light curves for which we estimated the time delay between different images. We find that using only LSST data for time-delay cosmography is not ideal. Instead, we advocate using LSST as a discovery machine for LSNe Ia, enabling time delay measurements from follow-up observations from other instruments in order to increase the number of systems by a factor of 2 to 16 depending on the observing strategy. Furthermore, we find that LSST observing strategies, which provide a good sampling frequency (the mean inter-night gap is around two days) and high cumulative season length (ten seasons with a season length of around 170 days per season), are favored. Rolling cadences subdivide the survey and focus on different parts in different years; these observing strategies trade the number of seasons for better sampling frequency. In our investigation, this leads to half the number of systems in comparison to the best observing strategy. Therefore rolling cadences are disfavored because the gain from the increased sampling frequency cannot compensate for the shortened cumulative season length. We anticipate that the sample of lensed SNe Ia from our preferred LSST cadence strategies with rapid follow-up observations would yield an independent percent-level constraint on $H_0$.

Faraday rotation measure at radio wavelengths is commonly used to diagnose large-scale magnetic fields. It is argued that the length-scales on which magnetic fields vary in large-scale diffuse astrophysical media can be inferred from correlations in the observed RM. RM is a variable which can be derived from the polarised radiative transfer equations in restrictive conditions. This paper assesses the usage of RMF (rotation measure fluctuation) analyses for magnetic field diagnostics in the framework of polarised radiative transfer. We use models of various magnetic field configurations and electron density distributions to show how density fluctuations could affect the correlation length of the magnetic fields inferred from the conventional RMF analyses. We caution against interpretations of RMF analyses when a characteristic density is ill-defined, e.g. in cases of log-normal distributed and fractal-like density structures. As the spatial correlations are generally not the same in the line-of-sight longitudinal direction and the sky plane direction, one also needs to clarify the context of RMF when inferring from observational data. In complex situations, a covariant polarised radiative transfer calculation is essential to capture all aspects of radiative and transport processes, which would otherwise ambiguate the interpretations of magnetism in galaxy clusters and larger-scale cosmological structures.

We present the completed KMOS$^\mathrm{3D}$ survey $-$ an integral field spectroscopic survey of 739, $\log(M_{\star}/M_{\odot})>9$, galaxies at $0.6<z<2.7$ using the K-band Multi Object Spectrograph (KMOS) at the Very Large Telescope (VLT). KMOS$^\mathrm{3D}$ provides a population-wide census of kinematics, star formation, outflows, and nebular gas conditions both on and off the star-forming galaxy main sequence through the spatially resolved and integrated properties of H$\alpha$, [N II], and [S II] emission lines. We detect H$\alpha$ emission for 91% of galaxies on the main sequence of star-formation and 79% overall. The depth of the survey has allowed us to detect galaxies with star-formation rates below 1 M$_{\odot}$/ yr$^{-1}$, as well as to resolve 81% of detected galaxies with $\geq3$ resolution elements along the kinematic major axis. The detection fraction of H$\alpha$ is a strong function of both color and offset from the main sequence, with the detected and non-detected samples exhibiting different SED shapes. Comparison of H$\alpha$ and UV+IR star formation rates (SFRs) reveal that dust attenuation corrections may be underestimated by 0.5 dex at the highest masses ($\log(M_{\star}/M_{\odot})>10.5$). We confirm our first year results of a high rotation dominated fraction (monotonic velocity gradient and $v_\mathrm{rot}$/$\sigma_0 > \sqrt{3.36}$) of 77% for the full KMOS$^\mathrm{3D}$ H$\alpha$sample. The rotation-dominated fraction is a function of both stellar mass and redshift with the strongest evolution measured over the redshift range of the survey for galaxies with $\log(M_{\star}/M_{\odot})<10.5$. With this paper we include a final data release of all 739 observed objects.

Context: The determination of the stellar initial mass function (IMF) of massive galaxies is one of the open problems in cosmology. Strong gravitational lensing is one of the few methods that allow us to constrain the IMF outside of the Local Group. Aims: The goal of this study is to statistically constrain the distribution in the IMF mismatch parameter, defined as the ratio between the true stellar mass of a galaxy and that inferred assuming a reference IMF, of massive galaxies from the BOSS CMASS sample. Methods: We take 23 strong lenses drawn from the CMASS sample, measure their Einstein radii and stellar masses using multi-band photometry from the Hyper Suprime-Cam survey, then fit a model distribution for the IMF mismatch parameter and dark matter halo mass to the whole sample. We use a prior on halo mass from weak lensing measurements and account for strong lensing selection effects in our model. Results: Assuming an NFW density profile for the dark matter distribution, we infer a value $\mu_{\mathrm{IMF}} = -0.04\pm0.11$ for the average base-10 logarithm of the IMF mismatch parameter, defined with respect to a Chabrier IMF. A Salpeter IMF is in tension with our measurements. Conclusions: Our results are consistent with a scenario in which the region of massive galaxies where the IMF normalization is significantly heavier than that of the Milky Way is much smaller than the scales $5\sim10$~kpc probed by the Einstein radius of the lenses in our sample, as recent spatially resolved studies of the IMF in massive galaxies suggest. The MCMC chains describing the posterior probability distribution of the model are available online, together with the code used to obtain them.

We present results on the rest-frame $H$-band luminosity functions (LF) of red sequence galaxies in seven clusters at 1.0 < z < 1.3 from the Gemini Observations of Galaxies in Rich Early Environments Survey (GOGREEN). Using deep GMOS-z' and IRAC $3.6 \mu$m imaging, we identify red sequence galaxies and measure their LFs down to $M_{H} \sim M_{H}^{*} + (2.0 - 3.0)$. By stacking the entire sample, we derive a shallow faint end slope of $ \alpha \sim -0.35^{+0.15}_{-0.15} $ and $ M_{H}^{*} \sim -23.52^{+0.15}_{-0.17} $, suggesting that there is a deficit of faint red sequence galaxies in clusters at high redshift. By comparing the stacked red sequence LF of our sample with a sample of clusters at z~0.6, we find an evolution in the faint end of the red sequence over the ~2.6 Gyr between the two samples, with the mean faint end red sequence luminosity growing by more than a factor of two. The faint-to-luminous ratio of our sample ($0.78^{+0.19}_{-0.15}$) is consistent with the trend of decreasing ratio with increasing redshift as proposed in previous studies. A comparison with the field shows that the faint-to-luminous ratios in clusters are consistent with the field at z~1.15 and exhibit a stronger redshift dependence. Our results support the picture that the build up of the faint red sequence galaxies occurs gradually over time and suggest that faint cluster galaxies, similar to bright cluster galaxies, experience the quenching effect induced by environment already at z~1.15.

We present the measurement of the Hubble Constant, $H_0$, with three strong gravitational lens systems. We describe a blind analysis of both PG1115+080 and HE0435-1223 as well as an extension of our previous analysis of RXJ1131-1231. For each lens, we combine new adaptive optics (AO) imaging from the Keck Telescope, obtained as part of the SHARP AO effort, with Hubble Space Telescope (HST) imaging, velocity dispersion measurements, and a description of the line-of-sight mass distribution to build an accurate and precise lens mass model. This mass model is then combined with the COSMOGRAIL measured time delays in these systems to determine $H_{0}$. We do both an AO-only and an AO+HST analysis of the systems and find that AO and HST results are consistent. After unblinding, the AO-only analysis gives $H_{0}=82.8^{+9.4}_{-8.3}~\rm km\,s^{-1}\,Mpc^{-1}$ for PG1115+080, $H_{0}=70.1^{+5.3}_{-4.5}~\rm km\,s^{-1}\,Mpc^{-1}$ for HE0435-1223, and $H_{0}=77.0^{+4.0}_{-4.6}~\rm km\,s^{-1}\,Mpc^{-1}$ for RXJ1131-1231. The joint AO-only result for the three lenses is $H_{0}=75.6^{+3.2}_{-3.3}~\rm km\,s^{-1}\,Mpc^{-1}$. The joint result of the AO+HST analysis for the three lenses is $H_{0}=76.8^{+2.6}_{-2.6}~\rm km\,s^{-1}\,Mpc^{-1}$. All of the above results assume a flat $\Lambda$ cold dark matter cosmology with a uniform prior on $\Omega_{\textrm{m}}$ in [0.05, 0.5] and $H_{0}$ in [0, 150] $\rm km\,s^{-1}\,Mpc^{-1}$. This work is a collaboration of the SHARP and H0LiCOW teams, and shows that AO data can be used as the high-resolution imaging component in lens-based measurements of $H_0$. The full time-delay cosmography results from a total of six strongly lensed systems are presented in a companion paper.

We present the $0.6<z<2.6$ evolution of the ionized gas velocity dispersion in 175 star-forming disk galaxies based on data from the full KMOS$^{\rm 3D}$ integral field spectroscopic survey. In a forward-modelling Bayesian framework including instrumental effects and beam-smearing, we fit simultaneously the observed galaxy velocity and velocity dispersion along the kinematic major axis to derive the intrinsic velocity dispersion $\sigma_0$. We find a reduction of the average intrinsic velocity dispersion of disk galaxies as a function of cosmic time, from $\sigma_0\sim45$ km s$^{-1}$ at $z\sim2.3$ to $\sigma_0\sim30$ km s$^{-1}$ at $z\sim0.9$. There is substantial intrinsic scatter ($\sigma_{\sigma_0, {\rm int}}\approx10$ km s$^{-1}$) around the best-fit $\sigma_0-z$-relation beyond what can be accounted for from the typical measurement uncertainties ($\delta\sigma_0\approx12$ km s$^{-1}$), independent of other identifiable galaxy parameters. This potentially suggests a dynamic mechanism such as minor mergers or variation in accretion being responsible for the scatter. Putting our data into the broader literature context, we find that ionized and atomic+molecular velocity dispersions evolve similarly with redshift, with the ionized gas dispersion being $\sim10-15$ km s$^{-1}$ higher on average. We investigate the physical driver of the on average elevated velocity dispersions at higher redshift, and find that our galaxies are at most marginally Toomre-stable, suggesting that their turbulent velocities are powered by gravitational instabilities, while stellar feedback as a driver alone is insufficient. This picture is supported through comparison with a state-of-the-art analytical model of galaxy evolution.

We present a blind time-delay strong lensing (TDSL) cosmographic analysis of the doubly imaged quasar SDSS 1206+4332. We combine the relative time delay between the quasar images, Hubble Space Telescope imaging, the Keck stellar velocity dispersion of the lensing galaxy, and wide-field photometric and spectroscopic data of the field to constrain two angular diameter distance relations. The combined analysis is performed by forward modelling the individual data sets through a Bayesian hierarchical framework, and it is kept blind until the very end to prevent experimenter bias. After unblinding, the inferred distances imply a Hubble constant $H_0 = 68.8^{+5.4}_{-5.1}$ kms$^{-1}$Mpc$^{-1}$, assuming a flat Lambda cold dark matter cosmology with uniform prior on $\Omega_{\rm m}$ in [0.05, 0.5]. The precision of our cosmographic measurement with the doubly imaged quasar SDSS 1206+4332 is comparable with those of quadruply imaged quasars and opens the path to perform on selected doubles the same analysis as anticipated for quads. Our analysis is based on a completely independent lensing code than our previous three H0LiCOW systems and the new measurement is fully consistent with those. We provide the analysis scripts paired with the publicly available software to facilitate independent analysis. The consistency between blind measurements with independent codes provides an important sanity check on lens modelling systematics. By combining the likelihoods of the four systems under the same prior, we obtain $H_0 = 72.5^{+2.1}_{-2.3}$kms$^{-1}$Mpc$^{-1}$. This measurement is independent of the distance ladder and other cosmological probes.

We present new measurements of the time delays of WFI2033-4723. The data sets used in this work include 14 years of data taken at the 1.2m Leonhard Euler Swiss telescope, 13 years of data from the SMARTS 1.3m telescope at Las Campanas Observatory and a single year of high-cadence and high-precision monitoring at the MPIA 2.2m telescope. The time delays measured from these different data sets, all taken in the R-band, are in good agreement with each other and with previous measurements from the literature. Combining all the time-delay estimates from our data sets results in Dt_AB = 36.2-0.8+0.7 days (2.1% precision), Dt_AC = -23.3-1.4+1.2 days (5.6%) and Dt_BC = -59.4-1.3+1.3 days (2.2%). In addition, the close image pair A1-A2 of the lensed quasars can be resolved in the MPIA 2.2m data. We measure a time delay consistent with zero in this pair of images. We also explore the prior distributions of microlensing time-delay potentially affecting the cosmological time-delay measurements of WFI2033-4723. There is however no strong indication in our measurements that microlensing time delay is neither present nor absent. This work is part of a H0LiCOW series focusing on measuring the Hubble constant from WFI2033-4723.

It has recently been proposed that gravitationally lensed type-Ia supernovae can provide microlensing-free time-delay measurements provided that the measurement is taken during the achromatic expansion phase of the explosion and that color light curves are used rather than single-band light curves. If verified, this would provide both precise and accurate time-delay measurements, making lensed type-Ia supernovae a new golden standard for time-delay cosmography. However, the 3D geometry of the expanding shell can introduce an additional bias that has not yet been fully explored. In this work, we present and discuss the impact of this effect on time-delay cosmography with lensed supernovae and find that on average it leads to a bias of a few tenths of a day for individual lensed systems. This is negligible in view of the cosmological time delays predicted for typical lensed type-Ia supernovae but not for the specific case of the recently discovered type-Ia supernova iPTF16geu, whose time delays are expected to be smaller than a day.

Strong-gravitational lens systems with quadruply-imaged quasars (quads) are unique probes to address several fundamental problems in cosmology and astrophysics. Although they are intrinsically very rare, ongoing and planned wide-field deep-sky surveys are set to discover thousands of such systems in the next decade. It is thus paramount to devise a general framework to model strong-lens systems to cope with this large influx without being limited by expert investigator time. We propose such a general modelling framework (implemented with the publicly available software Lenstronomy) and apply it to uniformly model three-band Hubble Space Telescope Wide Field Camera 3 images of 13 quads. This is the largest uniformly modelled sample of quads to date and paves the way for a variety of studies. To illustrate the scientific content of the sample, we investigate the alignment between the mass and light distribution in the deflectors. The position angles of these distributions are well-aligned, except when there is strong external shear. However, we find no correlation between the ellipticity of the light and mass distributions. We also show that the observed flux-ratios between the images depart significantly from the predictions of simple smooth models. The departures are strongest in the bluest band, consistent with microlensing being the dominant cause in addition to millilensing. Future papers will exploit this rich dataset in combination with ground based spectroscopy and time delays to determine quantities such as the Hubble constant, the free streaming length of dark matter, and the normalization of the initial stellar mass function.

Polarization of radiation is a powerful tool to study cosmic magnetism and analysis of polarization can be used as a diagnostic tool for large-scale structures. In this paper, we present a solid theoretical foundation for using polarized light to investigate large-scale magnetic field structures: the cosmological polarized radiative transfer (CPRT) formulation. The CPRT formulation is fully covariant. It accounts for cosmological and relativistic effects in a self-consistent manner and explicitly treats Faraday rotation, as well as Faraday conversion, emission, and absorption processes. The formulation is derived from the first principles of conservation of phase-space volume and photon number. Without loss of generality, we consider a flat Friedmann-Robertson-Walker (FRW) space-time metric and construct the corresponding polarized radiative transfer equations. We propose an all-sky CPRT calculation algorithm, based on a ray-tracing method, which allows cosmological simulation results to be incorporated and, thereby, model templates of polarization maps to be constructed. Such maps will be crucial in our interpretation of polarized data, such as those to be collected by the Square Kilometer Array (SKA). We describe several tests which are used for verifying the code and demonstrate applications in the study of the polarization signatures in different distributions of electron number density and magnetic fields. We present a pencil-beam CPRT calculation and an all-sky calculation, using a simulated galaxy cluster or a model magnetized universe obtained from GCMHD+ simulations as the respective input structures. The implications on large-scale magnetic field studies are discussed; remarks on the standard methods using rotation measure are highlighted.

Minor mergers are thought to be responsible for the size growth of quiescent field galaxies with decreasing redshift. We test this hypothesis using the cluster environment as a laboratory. Satellite galaxies in clusters move at high velocities, making mergers between them rare. The stellar mass-size relation in ten clusters and in the field is measured and compared at $z~\mathtt{\sim}~1$. Our cluster sample contains 344 spectroscopically-confirmed cluster members with Gemini/GMOS and 182 confirmed with HST WFC3 G141 grism spectroscopy. On average, quiescent and star-forming cluster galaxies are smaller than their field counterparts by ($0.08\pm0.04$) dex and ($0.07\pm0.01$) dex respectively. These size offsets are consistent with the average sizes of quiescent and star-forming field galaxies between $1.2\leqslant z\leqslant1.5$, implying the cluster environment has inhibited size growth between this period and $z~\mathtt{\sim}~1$. The negligible differences measured between the $z~\mathtt{\sim}~0$ field and cluster quiescent mass-size relations in other works imply that the average size of quiescent cluster galaxies must rise with decreasing redshift. Using a toy model, we show that the disappearance of the compact cluster galaxies might be explained if, on average, $\mathtt{\sim}40\%$ of them merge with their brightest cluster galaxies (BCGs) and $\mathtt{\sim}60\%$ are tidally destroyed into the intra-cluster light (ICL) between $0\leqslant z\leqslant1$. This is in agreement with the observed stellar mass growth of BCGs between $0\leqslant z\leqslant1$ and the observed ICL stellar mass fraction at $z~\mathtt{\sim}~0$. Our results support minor mergers as the cause for the size growth in quiescent field galaxies, with cluster-specific processes responsible for the similarity between the field and cluster quiescent mass-size relations at low redshift.

We investigate the local and line-of-sight overdensities of strong gravitational lens galaxies using wide-area multiband imaging from the Hyper Suprime-Cam Subaru Strategic Program. We present 41 new definite or probable lens candidates discovered in Data Release 2 of the survey. Using a combined sample of 87 galaxy-scale lenses out to a lens redshift of $z_{\mathrm{L}} \sim 0.8$, we compare galaxy number counts in lines of sight toward known and newly-discovered lenses in the survey to those of a control sample consisting of random lines of sight. We also compare the local overdensity of lens galaxies to a sample of "twin" galaxies that have a similar redshift and velocity dispersion to test whether lenses lie in different environments from similar non-lens galaxies. We find that lens fields contain higher number counts of galaxies compared to the control fields, but this effect arises from the local environment of the lens. Once galaxies in the lens plane are removed, the lens lines of sight are consistent with the control sample. The local environments of the lenses are overdense compared to the control sample, and are slightly overdense compared to those of the twin sample, although the significance is marginal. There is no significant evidence of the evolution of the local overdensity of lens environments with redshift.

The primary goals of the STRong lensing Insights into the Dark Energy Survey (STRIDES) collaboration are to measure the dark energy equation of state parameter and the free streaming length of dark matter. To this aim, STRIDES is discovering strongly lensed quasars in the imaging data of the Dark Energy Survey and following them up to measure time delays, high resolution imaging, and spectroscopy sufficient to construct accurate lens models. In this paper, we first present forecasts for STRIDES. Then, we describe the STRIDES classification scheme, and give an overview of the Fall 2016 follow-up campaign. We continue by detailing the results of two selection methods, the Outlier Selection Technique and a morphological algorithm, and presenting lens models of a system, which could possibly be a lensed quasar in an unusual configuration. We conclude with the summary statistics of the Fall 2016 campaign. Including searches presented in companion papers (Anguita et al.; Ostrovski et al.), STRIDES followed up 117 targets identifying 7 new strongly lensed systems, and 7 nearly identical quasars (NIQs), which could be confirmed as lenses by the detection of the lens galaxy. 76 candidates were rejected and 27 remain otherwise inconclusive, for a success rate in the range 6-35\%. This rate is comparable to that of previous searches like SQLS even though the parent dataset of STRIDES is purely photometric and our selection of candidates cannot rely on spectroscopic information.

Time-delay strong lensing provides a unique way to directly measure the Hubble constant ($H_{0}$). The precision of the $H_{0}$ measurement depends on the uncertainties in the time-delay measurements, the mass distribution of the main deflector(s), and the mass distribution along the line of sight. Tie and Kochanek (2018) have proposed a new microlensing effect on time delays based on differential magnification of the coherent accretion disc variability of the lensed quasar. If real, this effect could significantly broaden the uncertainty on the time delay measurements by up to $30\%$ for lens systems such as PG1115+080, which have relatively short time delays and monitoring over several different epochs. In this paper we develop a new technique that uses the time-delay ratios and simulated microlensing maps within a Bayesian framework in order to limit the allowed combinations of microlensing delays and thus to lessen the uncertainties due to the proposed effect. We show that, under the assumption of Tie and Kochanek (2018), the uncertainty on the time-delay distance ($D_{\Delta t}$, which is proportional to 1/$H_{0}$) of short time-delay ($\sim18$ days) lens, PG1115+080, increases from $\sim7\%$ to $\sim10\%$ by simultaneously fitting the three time-delay measurements from the three different datasets across twenty years, while in the case of long time-delay ($\sim90$ days) lens, the microlensing effect on time delays is negligible as the uncertainty on $D_{\Delta t}$ of RXJ1131-1231 only increases from $\sim2.5\%$ to $\sim2.6\%$.

We present time-delay estimates for the quadruply imaged quasar PG 1115+080. Our resuls are based on almost daily observations for seven months at the ESO MPIA 2.2m telescope at La Silla Observatory, reaching a signal-to-noise ratio of about 1000 per quasar image. In addition, we re-analyse existing light curves from the literature that we complete with an additional three seasons of monitoring with the Mercator telescope at La Palma Observatory. When exploring the possible source of bias we consider the so-called microlensing time delay, a potential source of systematic error so far never directly accounted for in previous time-delay publications. In fifteen years of data on PG 1115+080, we find no strong evidence of microlensing time delay. Therefore not accounting for this effect, our time-delay estimates on the individual data sets are in good agreement with each other and with the literature. Combining the data sets, we obtain the most precise time-delay estimates to date on PG 1115+080, with Dt(AB) = 8.3+1.5-1.6 days (18.7% precision), Dt(AC) = 9.9+1.1-1.1 days (11.1%) and Dt(BC) = 18.8+1.6-1.6 days (8.5%). Turning these time delays into cosmological constraints is done in a companion paper that makes use of ground-based Adaptive Optics (AO) with the Keck telescope.

Wave dark matter ($\psi$DM) predicts a compact soliton core and a granular halo in every galaxy. This work presents the first simulation study of an elliptical galaxy by including both stars and $\psi$DM, focusing on the systematic changes of the central soliton and halo granules. With the addition of stars in the inner halo, we find the soliton core consistently becomes more prominent by absorbing mass from the host halo than that without stars, and the halo granules become "non-isothermal", "hotter" in the inner halo and "cooler" in the outer halo, as opposed to the isothermal halo in pure $\psi$DM cosmological simulations. Moreover, the composite (star+$\psi$DM) mass density is found to follow a $r^{-2}$ isothermal profile near the half-light radius in most cases. Most striking is the velocity dispersion of halo stars that increases rapidly toward the galactic center by a factor of at least 2 inside the half-light radius caused by the deepened soliton gravitational potential, a result that compares favorably with observations of elliptical galaxies and bulges in spiral galaxies. However in some rare situations we find a phase segregation turning a compact distribution of stars into two distinct populations with high and very low velocity dispersions; while the high-velocity component mostly resides in the halo, the very low-velocity component is bound to the interior of the soliton core, resembling stars in faint dwarf spheroidal galaxies.

Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

Using integral field spectroscopy we investigate the kinematic properties of 35 massive centrally-dense and compact star-forming galaxies (${\log{\overline{M}_*}}=11.1$, $\log{(\Sigma_\mathrm{1kpc})}>9.5$, $\log{(M_\ast/r_e^{1.5})}>10.3$) at $z\sim0.7-3.7$ within the KMOS$^\mathrm{3D}$survey. We spatially resolve 23 compact star-forming galaxies (SFGs) and find that the majority are dominated by rotational motions with velocities ranging from {$95-500$ km s$^{-1}$}. The range of rotation velocities is reflected in a similar range of integrated H$\alpha$ linewidths, $75-400$ km s$^{-1}$, consistent with the kinematic properties of mass-matched extended galaxies from the full KMOS$^\mathrm{3D}$ sample. The fraction of compact SFGs that are classified as `rotation-dominated' or `disk-like' also mirrors the fractions of the full KMOS$^\mathrm{3D}$ sample. We show that integrated line-of-sight gas velocity dispersions from KMOS$^\mathrm{3D}$ are best approximated by a linear combination of their rotation and turbulent velocities with a lesser but still significant contribution from galactic scale winds. The H$\alpha$ exponential disk sizes of compact SFGs are on average $2.5\pm0.2$ kpc, $1-2\times$ the continuum sizes, in agreement with previous work. The compact SFGs have a $1.4\times$ higher AGN incidence than the full KMOS$^\mathrm{3D}$ sample at fixed stellar mass with average AGN fraction of 76\%. Given their high and centrally concentrated stellar masses as well as stellar to dynamical mass ratios close to unity, the compact SFGs are likely to have low molecular gas fractions and to quench on a short time scale unless replenished with inflowing gas. The rotation in these compact systems suggests that their direct descendants are rotating passive galaxies.

Using a sample of 4 galaxy clusters at $1.35 < z < 1.65$ and 10 galaxy clusters at $0.85 < z < 1.35$, we measure the environmental quenching timescale, $t_Q$, corresponding to the time required after a galaxy is accreted by a cluster for it to fully cease star formation. Cluster members are selected by a photometric-redshift criterion, and categorized as star-forming, quiescent, or intermediate according to their dust-corrected rest-frame colors and magnitudes. We employ a "delayed-then-rapid" quenching model that relates a simulated cluster mass accretion rate to the observed numbers of each type of galaxy in the cluster to constrain $t_Q$. For galaxies of mass $M_* \gtrsim 10^{10.5}~ \mathrm{M}_\odot$, we find a quenching timescale of $t_Q=$ 1.24 Gyr in the $z\sim1.5$ cluster sample, and $t_Q=$ 1.50 Gyr at $z\sim1$. Using values drawn from the literature, we compare the redshift evolution of $t_Q$ to timescales predicted for different physical quenching mechanisms. We find $t_Q$ to depend on host halo mass such that quenching occurs over faster timescales in clusters relative to groups, suggesting that properties of the host halo are responsible for quenching high-mass galaxies. Between $z=0$ and $z=1.5$, we find that $t_Q$ evolves faster than the molecular gas depletion timescale and slower than an SFR-outflow timescale, but is consistent with the evolution of the dynamical time. This suggests that environmental quenching in these galaxies is driven by the motion of satellites relative to the cluster environment, although due to uncertainties in the atomic gas budget at high redshift, we cannot rule out quenching due to simple gas depletion.

We present results on the structural properties of massive passive galaxies in three clusters at $1.39<z<1.61$ from the KMOS Cluster Survey. We measure light-weighted and mass-weighted sizes from optical and near-infrared Hubble Space Telescope imaging and spatially resolved stellar mass maps. The rest-frame $R$-band sizes of these galaxies are a factor of $\sim2-3$ smaller than their local counterparts. The slopes of the relation between the stellar mass and the light-weighted size are consistent with recent studies in clusters and the field. Their mass-weighted sizes are smaller than the rest frame $R$-band sizes, with an average mass-weighted to light-weighted size ratio that varies between $\sim0.45$ and $0.8$ among the clusters. We find that the median light-weighted size of the passive galaxies in the two more evolved clusters is $\sim24\%$ larger than for field galaxies, independent of the use of circularized effective radii or semi-major axes. These two clusters also show a smaller size ratio than the less evolved cluster, which we investigate using color gradients to probe the underlying $M_{*}/L_{H_{160}}$ gradients. The median color gradients are $\nabla{z-H} \sim-0.4$ mag dex$^{-1}$, twice the local value. Using stellar populations models, these gradients are best reproduced by a combination of age and metallicity gradients. Our results favor the minor merger scenario as the dominant process responsible for the observed galaxy properties and the environmental differences at this redshift. The environmental differences support that clusters experience accelerated structural evolution compared to the field, likely via an epoch of enhanced minor merger activity during cluster assembly.

The Hyper Suprime-Cam Subaru Strategic Program (HSC SSP) is an excellent survey for the search for strong lenses, thanks to its area, image quality and depth. We use three different methods to look for lenses among 43,000 luminous red galaxies from the Baryon Oscillation Spectroscopic Survey (BOSS) sample with photometry from the S16A internal data release of the HSC SSP. The first method is a newly developed algorithm, named YATTALENS, which looks for arc-like features around massive galaxies and then estimates the likelihood of an object being a lens by performing a lens model fit. The second method, CHITAH, is a modeling-based algorithm originally developed to look for lensed quasars. The third method makes use of spectroscopic data to look for emission lines from objects at a different redshift from that of the main galaxy. We find 15 definite lenses, 36 highly probable lenses and 282 possible lenses. Among the three methods, YATTALENS, which was developed specifically for this problem, performs best in terms of both completeness and purity. Nevertheless five highly probable lenses were missed by YATTALENS but found by the other two methods, indicating that the three methods are highly complementary. Based on these numbers we expect to find $\sim$300 definite or probable lenses by the end of the HSC SSP.

We present data for 16 galaxies in the overdensity JKCS 041 at $z \simeq 1.80$ as part of the K-band Multi-Object Spectrograph (KMOS) Cluster Survey (KCS). With 20-hour integrations, we have obtained deep absorption-line spectra from which we derived velocity dispersions for seven quiescent galaxies. We combined photometric parameters derived from Hubble Space Telescope images with the dispersions to construct a fundamental plane (FP) for quiescent galaxies in JKCS 041. From the zero-point evolution of the FP, we derived a formation redshift for the galaxies of $z_{form} = 3.0\pm0.3$, corresponding to a mean age of $1.4\pm0.2$ Gyrs. We tested the effect of structural and velocity dispersion evolution on our FP zero point and found a negligible contribution when using dynamical mass-normalized parameters ($\sim 3\%$), but a significant contribution from stellar-mass-normalized parameters ($\sim 42 \%$). From the relative velocities of the galaxies, we probed the three-dimensional structure of these 16 confirmed members of JKCS 041, and found that a group of galaxies in the south west of the overdensity had systematically higher velocities. We derived ages for the galaxies in the different groups from the FP. We found the east-extending group had typically older galaxies ($2.1\substack{+0.3\\-0.2}$ Gyrs), than those in the south-west group ($0.3\pm0.2$ Gyrs). Although based on small numbers, the overdensity dynamics, morphology, and age results could indicate that JKCS 041 is in formation and may comprise two merging groups of galaxies. The result could link large-scale structure to ages of galaxies for the first time at this redshift.

We present the analysis of the fundamental plane (FP) for a sample of 19 massive red-sequence galaxies ($M_{\star} >4\times10^{10} M_{\odot}$) in 3 known overdensities at $1.39<z<1.61$ from the KMOS Cluster Survey, a guaranteed time program with spectroscopy from the K-band Multi-Object Spectrograph (KMOS) at the VLT and imaging from the Hubble Space Telescope. As expected, we find that the FP zero-point in $B$ band evolves with redshift, from the value 0.443 of Coma to $-0.10\pm0.09$, $-0.19\pm0.05$, $-0.29\pm0.12$ for our clusters at $z=1.39$, $z=1.46$, and $z=1.61$, respectively. For the most massive galaxies ($\log M_{\star}/M_{\odot}>11$) in our sample, we translate the FP zero-point evolution into a mass-to-light-ratio $M/L$ evolution finding $\Delta \log M/L_{B}=(-0.46\pm0.10)z$, $\Delta \log M/L_{B}=(-0.52\pm0.07)z$, to $\Delta \log M/L_{B}=(-0.55\pm0.10)z$, respectively. We assess the potential contribution of the galaxies structural and stellar velocity dispersion evolution to the evolution of the FP zero-point and find it to be $\sim$6-35 % of the FP zero-point evolution. The rate of $M/L$ evolution is consistent with galaxies evolving passively. By using single stellar population models, we find an average age of $2.33^{+0.86}_{-0.51}$ Gyr for the $\log M_{\star}/M_{\odot}>11$ galaxies in our massive and virialized cluster at $z=1.39$, $1.59^{+1.40}_{-0.62}$ Gyr in a massive but not virialized cluster at $z=1.46$, and $1.20^{+1.03}_{-0.47}$ Gyr in a protocluster at $z=1.61$. After accounting for the difference in the age of the Universe between redshifts, the ages of the galaxies in the three overdensities are consistent within the errors, with possibly a weak suggestion that galaxies in the most evolved structure are older.

We describe a new Large Program in progress on the Gemini North and South telescopes: Gemini Observations of Galaxies in Rich Early Environments (GOGREEN). This is an imaging and deep spectroscopic survey of 21 galaxy systems at $1<z<1.5$, selected to span a factor $>10$ in halo mass. The scientific objectives include measuring the role of environment in the evolution of low-mass galaxies, and measuring the dynamics and stellar contents of their host haloes. The targets are selected from the SpARCS, SPT, COSMOS and SXDS surveys, to be the evolutionary counterparts of today's clusters and groups. The new red-sensitive Hamamatsu detectors on GMOS, coupled with the nod-and-shuffle sky subtraction, allow simultaneous wavelength coverage over $\lambda\sim 0.6$--$1.05\mu$m, and this enables a homogeneous and statistically complete redshift survey of galaxies of all types. The spectroscopic sample targets galaxies with AB magnitudes $z^{\prime}<24.25$ and [3.6]$\mu$m$<22.5$, and is therefore statistically complete for stellar masses $M_\ast\gtrsim10^{10.3}M_\odot$, for all galaxy types and over the entire redshift range. Deep, multiwavelength imaging has been acquired over larger fields for most systems, spanning $u$ through $K$, in addition to deep IRAC imaging at 3.6$\mu$m. The spectroscopy is $\sim 50$ per cent complete as of semester 17A, and we anticipate a final sample of $\sim 500$ new cluster members. Combined with existing spectroscopy on the brighter galaxies from GCLASS, SPT and other sources, GOGREEN will be a large legacy cluster and field galaxy sample at this redshift that spectroscopically covers a wide range in stellar mass, halo mass, and clustercentric radius.

We investigate the stellar mass and baryonic mass Tully-Fisher relations (TFRs) of massive star-forming disk galaxies at redshift z~2.3 and z~0.9 as part of the KMOS^3D integral field spectroscopy survey. Our spatially resolved data allow reliable modelling of individual galaxies, including the effect of pressure support on the inferred gravitational potential. At fixed circular velocity, we find higher baryonic masses and similar stellar masses at z~2.3 as compared to z~0.9. Together with the decreasing gas-to-stellar mass ratios with decreasing redshift, this implies that the contribution of dark matter to the dynamical mass at the galaxy scale increases towards lower redshift. A comparison to local relations reveals a negative evolution of the stellar and baryonic TFR zero-points from z=0 to z~0.9, no evolution of the stellar TFR zero-point from z~0.9 to z~2.3, and a positive evolution of the baryonic TFR zero-point from z~0.9 to z~2.3. We discuss a toy model of disk galaxy evolution to explain the observed, non-monotonic TFR evolution, taking into account the empirically motivated redshift dependencies of galactic gas fractions, and of the relative amount of baryons to dark matter on galaxy and halo scales.

We explore the H-alpha emission in the massive quiescent galaxies observed by the KMOS-3D survey at 0.7 < z < 2.7. The H-alpha line is robustly detected in 20 out of 120 UVJ-selected quiescent galaxies, and we classify the emission mechanism using the H-alpha line width and the [NII]/H-alpha line ratio. We find that AGN are likely to be responsible for the line emission in more than half of the cases. We also find robust evidence for star formation activity in nine quiescent galaxies, which we explore in detail. The H-alpha kinematics reveal rotating disks in five of the nine galaxies. The dust-corrected H-alpha star formation rates are low (0.2 - 7 Msun/yr), and place these systems significantly below the main sequence. The 24micron-based infrared luminosities, instead, overestimate the star formation rates. These galaxies present a lower gas-phase metallicity compared to star-forming objects with similar stellar mass, and many of them have close companions. We therefore conclude that the low-level star formation activity in these nine quiescent galaxies is likely to be fueled by inflowing gas or minor mergers, and could be a sign of rejuvenation events.

We present a list of galaxy-scale lens candidates including a highly probable interacting galaxy-scale lens in the Hyper Suprime-Cam (HSC) imaging survey. We combine HSC imaging with the blended-spectra catalog from the Galaxy And Mass Assembly (GAMA) survey to identify lens candidates, and use lens mass modeling to confirm the candidates. We find 46 matches between the HSC S14A_0b imaging data release and the GAMA catalog. Ten of them are probable lens systems according to their morphology and redshifts. There is one system with an interacting galaxy pair, HSC J084928+000949, that has a valid mass model. We predict the total mass enclosed by the Einstein radius of $\sim0.72$" ($\sim1.65$kpc) for this new expected lens system to be $\sim10^{10.59}M_{\odot}$. Using the photometry in the {\it grizy} bands of the HSC survey and stellar population synthesis modeling with a Salpeter stellar initial mass function, we estimate the stellar mass within the Einstein radius to be $\sim10^{10.46}\,M_{\odot}$. We thus find a dark matter mass fraction within the Einstein radius of $\sim25\%$. Further spectroscopy or high-resolution imaging would allow confirmation of the nature of these lens candidates. The particular system with the interacting galaxy pair, if confirmed, would provide an opportunity to study the interplay between dark matter and stars as galaxies build up through hierarchical mergers.

We exploit the deep resolved Halpha kinematic data from the KMOS^3D and SINS/zC-SINF surveys to examine the largely unexplored outer disk kinematics of star-forming galaxies (SFGs) out to the peak of cosmic star formation. Our sample contains 101 SFGs representative of the more massive (9.3 < log(M*/Msun) < 11.5) main sequence population at 0.6<z<2.6. Through a novel stacking approach we are able to constrain a representative rotation curve extending out to ~4 effective radii. This average rotation curve exhibits a significant drop in rotation velocity beyond the turnover, with a slope of Delta(V)/Delta(R) = $-0.26^{+0.10}_{-0.09}$ in units of normalized coordinates V/V_max and R/R_turn. This result confirms that the fall-off seen previously in some individual galaxies is a common feature of our sample of high-z disks. We show that this outer fall-off strikingly deviates from the flat or mildly rising rotation curves of local spiral galaxies of similar masses. We furthermore compare our data with models including baryons and dark matter demonstrating that the falling stacked rotation curve can be explained by a high mass fraction of baryons relative to the total dark matter halo (m_d>~0.05) in combination with a sizeable level of pressure support in the outer disk. These findings are in agreement with recent studies demonstrating that star-forming disks at high redshift are strongly baryon dominated within the disk scale, and furthermore suggest that pressure gradients caused by large turbulent gas motions are present even in their outer disks. We demonstrate that these results are largely independent of our model assumptions such as the presence of a central stellar bulge, the effect of adiabatic contraction at fixed m_d, and variations in the concentration parameter.

In cold dark matter cosmology, the baryonic components of galaxies are thought to be mixed with and embedded in non-baryonic and non-relativistic dark matter, which dominates the total mass of the galaxy and its dark matter halo. In the local Universe, the mass of dark matter within a galactic disk increases with disk radius, becoming appreciable and then dominant in the outer, baryonic regions of the disks of star-forming galaxies. This results in rotation velocities of the visible matter within the disk that are constant or increasing with disk radius. Comparison between the dynamical mass and the sum of stellar and cold gas mass at the peak epoch of galaxy formation, inferred from ancillary data, suggest high baryon factions in the inner, star-forming regions of the disks. Although this implied baryon fraction may be larger than in the local Universe, the systematic uncertainties (stellar initial mass function, calibration of gas masses) render such comparisons inconclusive in terms of the mass of dark matter. Here we report rotation curves for the outer disks of six massive star-forming galaxies, and find that the rotation velocities are not constant, but decrease with radius. We propose that this trend arises because of two main factors: first, a large fraction of the massive, high-redshift galaxy population was strongly baryon dominated, with dark matter playing a smaller part than in the local Universe; and second, the large velocity dispersion in high-redshift disks introduces a substantial pressure term that leads to a decrease in rotation velocity with increasing radius. The effect of both factors appears to increase with redshift. Qualitatively, the observations suggest that baryons in the early Universe efficiently condensed at the centres of dark matter halos when gas fractions were high, and dark matter was less concentrated. [Abridged]

We report the serendipitous discovery of a quadruply lensed source at $z_{\rm s}=3.76$, HSC~J115252+004733, from the Hyper Suprime-Cam (HSC) Survey. The source is lensed by an early-type galaxy at $z_{\rm l}=0.466$ and a satellite galaxy. Here, we investigate the properties of the source by studying its size and luminosity from the imaging and the luminosity and velocity width of the Ly-$\alpha$ line from the spectrum. Our analyses suggest that the source is most probably a low-luminosity active galactic nucleus (LLAGN) but the possibility of it being a compact bright galaxy (e.g., a Lyman-$\alpha$ emitter or Lyman Break Galaxy) cannot be excluded. The brighter pair of lensed images appears point-like except in the HSC $i$-band (with a seeing $\sim0.5"$). The extended emission in the $i$-band image could be due to the host galaxy underneath the AGN, or alternatively, due to a highly compact lensed galaxy (without AGN) which appears point-like in all bands except in $i$-band. We also find that the flux ratio of the brighter pair of images is different in the Ks-band compared to optical wavelengths. Phenomena such as differential extinction and intrinsic variability cannot explain this chromatic variation. While microlensing from stars in the foreground galaxy is less likely to be the cause, it cannot be ruled out completely. If the galaxy hosts an AGN, then this represents the highest redshift quadruply imaged AGN known to date, enabling study of a distant LLAGN. Discovery of this unusually compact and faint source demonstrates the potential of the HSC survey.

We present gravitational lens models of the multiply imaged quasar DES J0408-5354, recently discovered in the Dark Energy Survey (DES) footprint, with the aim of interpreting its remarkable quad-like configuration. We first model the DES single-epoch $grizY$ images as a superposition of a lens galaxy and four point-like objects, obtaining spectral energy distributions (SEDs) and relative positions for the objects. Three of the point sources (A,B,D) have SEDs compatible with the discovery quasar spectra, while the faintest point-like image (G2/C) shows significant reddening and a `grey' dimming of $\approx0.8$mag. In order to understand the lens configuration, we fit different models to the relative positions of A,B,D. Models with just a single deflector predict a fourth image at the location of G2/C but considerably brighter and bluer. The addition of a small satellite galaxy ($R_{\rm E}\approx0.2$") in the lens plane near the position of G2/C suppresses the flux of the fourth image and can explain both the reddening and grey dimming. All models predict a main deflector with Einstein radius between $1.7"$ and $2.0",$ velocity dispersion $267-280$km/s and enclosed mass $\approx 6\times10^{11}M_{\odot},$ even though higher resolution imaging data are needed to break residual degeneracies in model parameters. The longest time-delay (B-A) is estimated as $\approx 85$ (resp. $\approx125$) days by models with (resp. without) a perturber near G2/C. The configuration and predicted time-delays of J0408-5354 make it an excellent target for follow-up aimed at understanding the source quasar host galaxy and substructure in the lens, and measuring cosmological parameters. We also discuss some lessons learnt from J0408-5354 on lensed quasar finding strategies, due to its chromaticity and morphology.

We report the discovery and spectroscopic confirmation of the quad-like lensed quasar system DES J0408-5354 found in the Dark Energy Survey (DES) Year 1 (Y1) data. This system was discovered during a search for DES Y1 strong lensing systems using a method that identified candidates as red galaxies with multiple blue neighbors. DES J0408-5354 consists of a central red galaxy surrounded by three bright (i < 20) blue objects and a fourth red object. Subsequent spectroscopic observations using the Gemini South telescope confirmed that the three blue objects are indeed the lensed images of a quasar with redshift z = 2.375, and that the central red object is an early-type lensing galaxy with redshift z = 0.597. DES J0408-5354 is the first quad lensed quasar system to be found in DES and begins to demonstrate the potential of DES to discover and dramatically increase the sample size of these very rare objects.

Curvelets are efficient to represent highly anisotropic signal content, such as a local linear and curvilinear structure. First-generation curvelets on the sphere, however, suffered from blocking artefacts. We present a new second-generation curvelet transform, where scale-discretised curvelets are constructed directly on the sphere. Scale-discretised curvelets exhibit a parabolic scaling relation, are well-localised in both spatial and harmonic domains, support the exact analysis and synthesis of both scalar and spin signals, and are free of blocking artefacts. We present fast algorithms to compute the exact curvelet transform, reducing computational complexity from $\mathcal{O}(L^5)$ to $\mathcal{O}(L^3\log_{2}{L})$ for signals band-limited at $L$. The implementation of these algorithms is made publicly available. Finally, we present an illustrative application demonstrating the effectiveness of curvelets for representing directional curve-like features in natural spherical images.

We make publicly available a catalog of calibrated environmental measures for galaxies in the five 3D-HST/CANDELS deep fields. Leveraging the spectroscopic and grism redshifts from the 3D-HST survey, multi wavelength photometry from CANDELS, and wider field public data for edge corrections, we derive densities in fixed apertures to characterize the environment of galaxies brighter than $JH_{140} < 24$ mag in the redshift range $0.5<z<3.0$. By linking observed galaxies to a mock sample, selected to reproduce the 3D-HST sample selection and redshift accuracy, each 3D-HST galaxy is assigned a probability density function of the host halo mass, and a probability that is a central or a satellite galaxy. The same procedure is applied to a $z=0$ sample selected from SDSS. We compute the fraction of passive central and satellite galaxies as a function of stellar and halo mass, and redshift, and then derive the fraction of galaxies that were quenched by environment specific processes. Using the mock sample, we estimate that the timescale for satellite quenching is $t_{\rm quench} \sim 2-5$ Gyr; longer at lower stellar mass or lower redshift, but remarkably independent of halo mass. This indicates that, in the range of environments commonly found within the 3D-HST sample, satellites are quenched by exhaustion of their gas reservoir in absence of cosmological accretion. We find that the quenching times can be separated into a delay phase during which satellite galaxies behave similarly to centrals at fixed stellar mass, and a phase where the star formation rate drops rapidly ($\sim 0.4-0.6$ Gyr), as shown previously at $z=0$. We conclude that this scenario requires satellite galaxies to retain a large reservoir of multi-phase gas upon accretion, even at high redshift, and that this gas sustains star formation for the long quenching times observed.

We exploit deep integral-field spectroscopic observations with KMOS/VLT of 240 star-forming disks at 0.6 < z < 2.6 to dynamically constrain their mass budget. Our sample consists of massive ($\gtrsim 10^{9.8} M_\odot$) galaxies with sizes $R_e \gtrsim 2$ kpc. By contrasting the observed velocity and dispersion profiles to dynamical models, we find that on average the stellar content contributes $32^{+8}_{-7}\%$ of the total dynamical mass, with a significant spread among galaxies (68th percentile range f_star ~ 18 - 62%). Including molecular gas as inferred from CO- and dust-based scaling relations, the estimated baryonic mass adds up to $56^{+17}_{-12}\%$ of total for the typical galaxy in our sample, reaching ~ 90% at z > 2. We conclude that baryons make up most of the mass within the disk regions of high-redshift star-forming disk galaxies, with typical disks at z > 2 being strongly baryon-dominated within $R_e$. Substantial object-to-object variations in both stellar and baryonic mass fractions are observed among the galaxies in our sample, larger than what can be accounted for by the formal uncertainties in their respective measurements. In both cases, the mass fractions correlate most strongly with measures of surface density. High $\Sigma_{star}$ galaxies feature stellar mass fractions closer to unity, and systems with high inferred gas or baryonic surface densities leave less room for additional mass components other than stars and molecular gas. Our findings can be interpreted as more extended disks probing further (and more compact disks probing less far) into the dark matter halos that host them.

We analyze the angular momenta of massive star forming galaxies (SFGs) at the peak of the cosmic star formation epoch (z~0.8-2.6). Our sample of ~360 log(M*/Msun) ~ 9.3-11.8 SFGs is mainly based on the KMOS3D and SINS/zC-SINF surveys of H$\alpha$ kinematics, and collectively provides a representative subset of the massive star forming population. The inferred halo scale angular momentum distribution is broadly consistent with that theoretically predicted for their dark matter halos, in terms of mean spin parameter <$\lambda$> ~ 0.037 and its dispersion ($\sigma_{log(\lambda)}$~0.2). Spin parameters correlate with the disk radial scale, and with their stellar surface density, but do not depend significantly on halo mass, stellar mass, or redshift. Our data thus support the long-standing assumption that on average, even at high redshifts, the specific angular momentum of disk galaxies reflects that of their dark matter halos (j_d = j_DM). The lack of correlation between $\lambda$ x (j_d/j_DM) and the nuclear stellar density $\Sigma_{*}$(1kpc) favors a scenario where disk-internal angular momentum redistribution leads to "compaction" inside massive high-redshift disks. For our sample, the inferred average stellar-to-dark matter mass ratio is ~2%, consistent with abundance matching results. Including the molecular gas, the total baryonic disk-to-dark matter mass ratio is ~5% for halos near $10^{12}$ Msun, which corresponds to 31% of the cosmologically available baryons, implying that high-redshift disks are strongly baryon dominated.

We present measurements of the [NII]/Ha ratio as a probe of gas-phase oxygen abundance for a sample of 419 star-forming galaxies at z=0.6-2.7 from the KMOS3D near-IR multi-IFU survey. The mass-metallicity relation (MZR) is determined consistently with the same sample selection, metallicity tracer, and methodology over the wide redshift range probed by the survey. We find good agreement with long-slit surveys in the literature, except for the low-mass slope of the relation at z~2.3, where this sample is less biased than previous samples based on optical spectroscopic redshifts. In this regime we measure a steeper slope than some literature results. Excluding the AGN contribution from the MZR reduces sensitivity at the high mass end, but produces otherwise consistent results. There is no significant dependence of the [NII]/Ha ratio on SFR or environment at fixed redshift and stellar mass. The IFU data allow spatially resolved measurements of [NII]/Ha, from which we can infer abundance gradients for 180 galaxies, thus tripling the current sample in the literature. The observed gradients are on average flat, with only 15 gradients statistically offset from zero at >3sigma. We have modelled the effect of beam-smearing, assuming a smooth intrinsic radial gradient and known seeing, inclination and effective radius for each galaxy. Our seeing-limited observations can recover up to 70% of the intrinsic gradient for the largest, face-on disks, but only 30% for the smaller, more inclined galaxies. We do not find significant trends between observed or corrected gradients and any stellar population, dynamical or structural galaxy parameters, mostly in agreement with existing studies with much smaller sample sizes. In cosmological simulations, strong feedback is generally required to produce flat gradients at high redshift.

We analyse the sizes, colour gradients, and resolved stellar mass distributions for 36 massive and passive galaxies in the cluster XMMUJ2235-2557 at z=1.39 using optical and near-infrared Hubble Space Telescope imaging. We derive light-weighted S\'ersic fits in five HST bands ($i_{775},z_{850},Y_{105},J_{125},H_{160}$), and find that the size decreases by ~20% going from $i_{775}$ to $H_{160}$ band, consistent with recent studies. We then generate spatially resolved stellar mass maps using an empirical relationship between $M_{*}/L_{H_{160}}$ and $(z_{850}-H_{160})$ and use these to derive mass-weighted S\'ersic fits: the mass-weighted sizes are ~41% smaller than their rest-frame $r$-band counterparts compared with an average of ~12% at z~0. We attribute this evolution to the evolution in the $M_{*}/L_{H_{160}}$ and colour gradient. Indeed, as expected, the ratio of mass-weighted to light-weighted size is correlated with the $M_{*}/L$ gradient, but is also mildly correlated with the mass surface density and mass-weighted size. The colour gradients $(\nabla_{z-H})$ are mostly negative, with a median value of $\sim0.45$ mag dex$^{-1}$, twice the local value. The evolution is caused by an evolution in age gradients along the semi-major axis ($a$), with $\nabla_{age} = d \log(age) / d \log(a)$ $\sim-0.33$, while the survival of weaker colour gradients in old, local galaxies implies that metallicity gradients are also required, with $\nabla_{Z} = d \log(Z) / d \log(a)$ $\sim-0.2$. This is consistent with recent observational evidence for the inside-out growth of passive galaxies at high redshift, and favours a gradual mass growth mechanism, such as minor mergers.

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

We investigate the stellar populations of 25 massive, galaxies ($\log[M_\ast/M_\odot] \geq 10.9$) at $1.5 < z < 2$ using data obtained with the K-band Multi-Object Spectrograph (KMOS) on the ESO VLT. Targets were selected to be quiescent based on their broadband colors and redshifts using data from the 3D-HST grism survey. The mean redshift of our sample is $\bar{z} = 1.75$, where KMOS YJ-band data probe age- and metallicity-sensitive absorption features in the rest-frame optical, including the $G$ band, Fe I, and high-order Balmer lines. Fitting simple stellar population models to a stack of our KMOS spectra, we derive a mean age of $1.03^{+0.13}_{-0.08}$ Gyr. We confirm previous results suggesting a correlation between color and age for quiescent galaxies, finding mean ages of $1.22^{+0.56}_{-0.19}$ Gyr and $0.85^{+0.08}_{-0.05}$ Gyr for the reddest and bluest galaxies in our sample. Combining our KMOS measurements with those obtained from previous studies at $0.2 < z < 2$ we find evidence for a $2-3$ Gyr spread in the formation epoch of massive galaxies. At $z < 1$ the measured stellar ages are consistent with passive evolution, while at $1 < z \lesssim2$ they appear to saturate at $\sim$1 Gyr, which likely reflects changing demographics of the (mean) progenitor population. By comparing to star-formation histories inferred for "normal" star-forming galaxies, we show that the timescales required to form massive galaxies at $z \gtrsim 1.5$ are consistent with the enhanced $\alpha$-element abundances found in massive local early-type galaxies.

We present the KMOS^3D survey, a new integral field survey of over 600 galaxies at 0.7<z<2.7 using KMOS at the Very Large Telescope (VLT). The KMOS^3D survey utilizes synergies with multi-wavelength ground and space-based surveys to trace the evolution of spatially-resolved kinematics and star formation from a homogeneous sample over 5 Gyrs of cosmic history. Targets, drawn from a mass-selected parent sample from the 3D-HST survey, cover the star formation-stellar mass ($M_*$) and rest-frame $(U-V)-M_*$ planes uniformly. We describe the selection of targets, the observations, and the data reduction. In the first year of data we detect Halpha emission in 191 $M_*=3\times10^{9}-7\times10^{11}$ Msun galaxies at z=0.7-1.1 and z=1.9-2.7. In the current sample 83% of the resolved galaxies are rotation-dominated, determined from a continuous velocity gradient and $v_{rot}/\sigma>1$, implying that the star-forming 'main sequence' (MS) is primarily composed of rotating galaxies at both redshift regimes. When considering additional stricter criteria, the Halpha kinematic maps indicate at least ~70% of the resolved galaxies are disk-like systems. Our high-quality KMOS data confirm the elevated velocity dispersions reported in previous IFS studies at z>0.7. For rotation-dominated disks, the average intrinsic velocity dispersion decreases by a factor of two from 50 km/s at z~2.3 to 25 km/s at z~0.9 while the rotational velocities at the two redshifts are comparable. Combined with existing results spanning z~0-3, disk velocity dispersions follow an approximate (1+z) evolution that is consistent with the dependence of velocity dispersion on gas fractions predicted by marginally-stable disk theory.

In this paper we follow up on our previous detection of nuclear ionized outflows in the most massive (log(M*/Msun) >= 10.9) z~1-3 star-forming galaxies (Forster Schreiber et al.), by increasing the sample size by a factor of six (to 44 galaxies above log(M*/Msun) >= 10.9) from a combination of the SINS/zC-SINF, LUCI, GNIRS, and KMOS^3D spectroscopic surveys. We find a fairly sharp onset of the incidence of broad nuclear emission (FWHM in the Ha, [NII], and [SII] lines ~ 450-5300 km/s), with large [NII]/Ha ratios, above log(M*/Msun) ~ 10.9, with about two thirds of the galaxies in this mass range exhibiting this component. Broad nuclear components near and above the Schechter mass are similarly prevalent above and below the main sequence of star-forming galaxies, and at z~1 and ~2. The line ratios of the nuclear component are fit by excitation from active galactic nuclei (AGN), or by a combination of shocks and photoionization. The incidence of the most massive galaxies with broad nuclear components is at least as large as that of AGNs identified by X-ray, optical, infrared or radio indicators. The mass loading of the nuclear outflows is near unity. Our findings provide compelling evidence for powerful, high-duty cycle, AGN-driven outflows near the Schechter mass, and acting across the peak of cosmic galaxy formation.

A new approximation method for inverting the Poisson's equation is presented for a continuously distributed and finite-sized source in an unbound domain. The advantage of this image multipole method arises from its ability to place the computational error close to the computational domain boundary, making the source region almost error free. It is contrasted to the modified Green's function method that has small but finite errors in the source region. Moreover, this approximation method also has a systematic way to greatly reduce the errors at the expense of somewhat greater computational efforts. Numerical examples of three-dimensional and two-dimensional cases are given to illustrate the advantage of the new method.

We present the correlations between stellar mass, star formation rate (SFR) and [NII]/Ha flux ratio as indicator of gas-phase metallicity for a sample of 222 galaxies at 0.8 < z < 2.6 and log(M*/Msun)=9.0-11.5 from the LUCI, SINS/zC-SINF and KMOS3D surveys. This sample provides a unique analysis of the mass-metallicity relation (MZR) over an extended redshift range using consistent data analysis techniques and strong-line metallicity indicator. We find a constant slope at the low-mass end of the relation and can fully describe its redshift evolution through the evolution of the characteristic turnover mass where the relation begins to flatten at the asymptotic metallicity. At fixed mass and redshift, our data do not show a correlation between the [NII]/Ha ratio and SFR, which disagrees with the 0.2-0.3dex offset in [NII]/Ha predicted by the "fundamental relation" between stellar mass, SFR and metallicity discussed in recent literature. However, the overall evolution towards lower [NII]/Ha at earlier times does broadly agree with these predictions.

The SAGE-Spec Spitzer Legacy program is a spectroscopic follow-up to the SAGE-LMC photometric survey of the Large Magellanic Cloud carried out with the Spitzer Space Telescope. We present an overview of SAGE-Spec and some of its first results. The SAGE-Spec program aims to study the life cycle of gas and dust in the Large Magellanic Cloud, and to provide information essential to the classification of the point sources observed in the earlier SAGE-LMC photometric survey. We acquired 224.6 hours of observations using the InfraRed Spectrograph and the SED mode of the Multiband Imaging Photometer for Spitzer. The SAGE-Spec data, along with archival Spitzer spectroscopy of objects in the Large Magellanic Cloud, are reduced and delivered to the community. We discuss the observing strategy, the specific data reduction pipelines applied and the dissemination of data products to the scientific community. Initial science results include the first detection of an extragalactic "21 um" feature towards an evolved star and elucidation of the nature of disks around RV Tauri stars in the Large Magellanic Cloud. Towards some young stars, ice features are observed in absorption. We also serendipitously observed a background quasar, at a redshift of z~0.14, which appears to be host-less.