We present the COSMOS2015 catalog which contains precise photometric redshifts and stellar masses for more than half a million objects over the 2deg$^{2}$ COSMOS field. Including new $YJHK_{\rm s}$ images from the UltraVISTA-DR2 survey, $Y$-band from Subaru/Hyper-Suprime-Cam and infrared data from the Spitzer Large Area Survey with the Hyper-Suprime-Cam Spitzer legacy program, this near-infrared-selected catalog is highly optimized for the study of galaxy evolution and environments in the early Universe. To maximise catalog completeness for bluer objects and at higher redshifts, objects have been detected on a $\chi^{2}$ sum of the $YJHK_{\rm s}$ and $z^{++}$ images. The catalog contains $\sim 6\times 10^5$ objects in the 1.5 deg$^{2}$ UltraVISTA-DR2 region, and $\sim 1.5\times 10^5$ objects are detected in the "ultra-deep stripes" (0.62 deg$^{2}$) at $K_{\rm s}\leq 24.7$ (3$\sigma$, 3", AB magnitude). Through a comparison with the zCOSMOS-bright spectroscopic redshifts, we measure a photometric redshift precision of $\sigma_{\Delta z/(1+z_s)}$ = 0.007 and a catastrophic failure fraction of $\eta=0.5$%. At $3

We present the X-ray spectral analysis of the 1855 extragalactic sources in the Chandra COSMOS-Legacy survey catalog having more than 30 net counts in the 0.5-7 keV band. 38% of the sources are optically classified Type 1 active galactic nuclei (AGN), 60% are Type 2 AGN and 2% are passive, low-redshift galaxies. We study the distribution of AGN photon index and of the intrinsic absorption N(H,z) based on the sources optical classification: Type 1 have a slightly steeper mean photon index than Type 2 AGN, which on the other hand have average intrinsic absorption ~3 times higher than Type 1 AGN. We find that ~15% of Type 1 AGN have N(H,z)>1E22 cm^(-2), i.e., are obscured according to the X-ray spectral fitting; the vast majority of these sources have L(2-10keV)>$1E44 erg/s. The existence of these objects suggests that optical and X-ray obscuration can be caused by different phenomena, the X-ray obscuration being for example caused by dust-free material surrounding the inner part of the nuclei. ~18% of Type 2 AGN have N(H,z)<1E22 cm^(-2), and most of these sources have low X-ray luminosities (L(2-10keV)<$1E43 erg/s). We expect a part of these sources to be low-accretion, unobscured AGN lacking of broad emission lines. Finally, we also find a direct proportional trend between N(H,z) and host galaxy mass and star formation rate, although part of this trend is due to a redshift selection effect.

We present the largest, most homogeneous catalogue of merging galaxies in the nearby universe obtained through the Galaxy Zoo project - an interface on the world-wide web enabling large-scale morphological classification of galaxies through visual inspection of images from the Sloan Digital Sky Survey (SDSS). The method converts a set of visually-inspected classifications for each galaxy into a single parameter (the `weighted-merger-vote fraction,' $f_m$) which describes our confidence that the system is part of an ongoing merger. We describe how $f_m$ is used to create a catalogue of 3003 visually-selected pairs of merging galaxies from the SDSS in the redshift range 0.005 &lt; z &lt;0.1. We use our merger sample and values of $f_m$ applied to the SDSS Main Galaxy Spectral sample (MGS) to estimate that the fraction of volume-limited (M_r &lt; -20.55) major mergers (1/3 &lt; {M}^*_1/{M}^*_2 &lt; 3) in the nearby universe is $1 - 3 \times C$ where $C \sim 1.5$ is a correction factor for spectroscopic incompleteness. Having visually classified the morphologies of the constituent galaxies in our mergers, we find that the spiral-to-elliptical ratio of galaxies in mergers is higher by a factor $\sim 2$ relative to the global population. In a companion paper, we examine the internal properties of these merging galaxies and conclude that this high spiral-to-elliptical ratio in mergers is due to a longer time-scale over which mergers with spirals are detectable compared to mergers with ellipticals.

We have utilized the largest sample of $\gamma$-ray selected Fermi flat-spectrum radio quasars (FSRQs) ever used (519 sources) to construct the luminosity function and its evolution through the cosmic history. In addition to spanning large redshift ($0

The relationship between galaxy size and environment has remained enigmatic, with over a decade of conflicting results. We present one of the first comprehensive studies of the variation of galaxy radius with environment beyond the local Universe and demonstrate that large-scale environmental density is correlated with galaxy radius independent of stellar mass and galaxy morphology. We confirm with &gt;5\sigma confidence that galaxies in denser environments are up to $\sim25\%$ larger than their equally massive counterparts with similar morphology in less dense regions of the Universe. We achieve this result by correlating projected two-dimensional densities over $\sim360$ deg$^2$ with the structural parameters of $\sim3$ million Hyper Suprime-Cam galaxies at 0.3 \leq z &lt; 0.7 with $\log M/M_{\odot} \geq 8.9$. Compared to most previous studies, this sample is $\sim100-10,000$ times larger and goes $\sim1$ dex deeper in mass-completeness. We demonstrate that past conflicting results have been driven by small sample sizes and a lack of robust measurement uncertainties. We verify the presence of the above correlation separately for disk-dominated, bulge-dominated, star-forming, and quiescent subpopulations. We find the strength of the correlation to be dependent on redshift, stellar mass, and morphology. The correlation is strongest at lower redshifts and systematically weakens or disappears beyond $z \geq 0.5$. At $z\geq0.5$, more massive galaxies still display a statistically significant correlation. Although some existing theoretical frameworks can be selectively invoked to explain some of the observed correlations, our work demonstrates the need for more comprehensive theoretical investigations of the correlation between galaxy size and environment.

We present a study of the connection between black hole accretion, star formation, and galaxy morphology at z<=2.5. We focus on active galactic nuclei (AGNs) selected by their mid-IR power-law emission. By fitting optical to far-IR photometry with state-of-the-art spectral energy distribution (SED) techniques, we derive stellar masses, star formation rates, dust properties, and AGN contributions in galaxies over the whole COSMOS field. We find that obscured AGNs lie within or slightly above the star-forming sequence. We confirm our previous finding about compact host galaxies of obscured AGNs at z~1, and find that galaxies with 20-50% AGN contributions tend to have smaller sizes, by ~25-50%, compared to galaxies without AGNs. Furthermore, we find that a high merger fraction of up to 0.5 is appropriate for the most luminous (log (LIR/Lsun) ~ 12.5) AGN hosts and non-AGN galaxies, but not for the whole obscured AGN sample. Moreover, merger fraction depends on the total and star-forming infrared luminosity, rather than the decomposed AGN infrared luminosity. Our results suggest that major mergers are not the main driver of AGN activity, and therefore obscured AGNs might be triggered by internal mechanisms, such as secular processes, disk instabilities, and compaction in a particular evolutionary stage. We make the SED modeling results publicly available.

Exploring the ICM power spectrum can help us to probe the physics of galaxy clusters. Using high-resolution 3D plasma simulations, we study the statistics of the velocity field and its relation with the thermodynamic perturbations. The normalization of the ICM spectrum (density, entropy, or pressure) is linearly tied to the level of large-scale motions, which excite both gravity and sound waves due to stratification. For low 3D Mach number M~0.25, gravity waves mainly drive entropy perturbations, traced by preferentially tangential turbulence. For M>0.5, sound waves start to significantly contribute, passing the leading role to compressive pressure fluctuations, associated with isotropic (or slightly radial) turbulence. Density and temperature fluctuations are then characterized by the dominant process: isobaric (low M), adiabatic (high M), or isothermal (strong conduction). Most clusters reside in the intermediate regime, showing a mixture of gravity and sound waves, hence drifting towards isotropic velocities. Remarkably, regardless of the regime, the variance of density perturbations is comparable to the 1D Mach number. This linear relation allows to easily convert between gas motions and ICM perturbations, which can be exploited by Chandra, XMM data and by the forthcoming Astro-H. At intermediate and small scales (10-100 kpc), the turbulent velocities develop a Kolmogorov cascade. The thermodynamic perturbations act as effective tracers of the velocity field, broadly consistent with the Kolmogorov-Obukhov-Corrsin advection theory. Thermal conduction acts to damp the gas fluctuations, washing out the filamentary structures and steepening the spectrum, while leaving unaltered the velocity cascade. The ratio of the velocity and density spectrum thus inverts the downtrend shown by the non-diffusive models, allowing to probe the presence of significant conductivity in the ICM.

We present deep 3 GHz VLA observations of the potentially recoiling black hole CID-42 in the COSMOS field. This galaxy shows two optical nuclei in the HST/ACS image and a large velocity offset of ~ 1300 km/s between the broad and narrow H beta emission line although the spectrum is not spacially resolved (Civano et al. 2010). The new 3 GHz VLA data has a bandwidth of 2 GHz and to correctly interpret the flux densities imaging was done with two different methods: multi-scale multi-frequency synthesis and spectral windows stacking. The final resolutions and sensitivities of these maps are 0.7" with rms = 4.6 muJy/beam and 0.9" with rms = 4.8 muJy/beam respectively. With a 7 sigma detection we find that the entire observed 3 GHz radio emission can be associated with the South-Eastern component of CID-42, coincident with the detected X-ray emission. We use our 3 GHz data combined with other radio data from the literature ranging from 320 MHz to 9 GHz, which include the VLA, VLBA and GMRT data, to construct a radio synchrotron spectrum of CID-42. The radio spectrum suggests a type I unobscured radio-quiet flat-spectrum AGN in the South-Eastern component which may be surrounded by a more extended region of old synchrotron electron population or shocks generated by the outflow from the supermassive black hole. Our data are consistent with the recoiling black hole picture but cannot rule out the presence of an obscured and radio-quiet SMBH in the North-Western component.

We obtained follow-up HST observations of the seven low surface brightness galaxies discovered with the Dragonfly Telephoto Array in the field of the massive spiral galaxy M101. Out of the seven galaxies, only three were resolved into stars and are potentially associated with the M101 group at $D=7\text{ Mpc}$. Based on HST ACS photometry in the broad F606W and F814W filters, we use a maximum likelihood algorithm to locate the Tip of the Red Giant Branch (TRGB) in galaxy color-magnitude diagrams. Distances are $6.38^{+0.35}_{-0.35}, 6.87^{+0.21}_{-0.30}$ and $6.52^{+0.25}_{-0.27} \text{ Mpc}$ and we confirm that they are members of the M101 group. Combining the three confirmed low luminosity satellites with previous results for brighter group members, we find the M101 galaxy group to be a sparsely populated galaxy group consisting of seven group members, down to $M_V = -9.2 \text{ mag}$. We compare the M101 cumulative luminosity function to that of the Milky Way and M31. We find that they are remarkably similar; In fact, the cumulative luminosity function of the M101 group gets even flatter for fainter magnitudes, and we show that the M101 group might exhibit the two known small-scale flaws in the $\Lambda\textrm{CDM}$ model, namely `the missing satellite' problem and the `too big to fail' problem. Kinematic measurements of M101$'$s satellite galaxies are required to determine whether the `too big to fail' problem does in fact exist in the M101 group.

We have imaged the Sunyaev-Zel'dovich (SZ) effect signals at 140 and 270 GHz towards ten galaxy clusters with Bolocam and AzTEC/ASTE. We also used Planck data to constrain the signal at large angular scales, Herschel-SPIRE images to subtract the brightest galaxies that comprise the cosmic infrared background (CIB), Chandra imaging to map the electron temperature $T_e$ of the intra-cluster medium (ICM), and HST imaging to derive models of each galaxy cluster's mass density. The galaxy clusters gravitationally lens the background CIB, which produced an on-average reduction in brightness towards the galaxy clusters' centers after the brightest galaxies were subtracted. We corrected for this deficit, which was between 5-25% of the 270 GHz SZ effect signal within $R_{2500}$. Using the SZ effect measurements, along with the X-ray constraint on $T_e$, we measured each galaxy cluster's average line of sight (LOS) velocity $v_z$ within $R_{2500}$, with a median per-cluster uncertainty of +-700 km/s. We found an ensemble-mean <$v_z$> of 430+-210 km/s, and an intrinsic cluster-to-cluster scatter $\sigma_{int}$ of 470+-340 km/s. We also obtained maps of $v_z$ over each galaxy cluster's face with an angular resolution of 70". All four galaxy clusters previously identified as having a merger oriented along the LOS showed an excess variance in these maps at a significance of 2-4$\sigma$, indicating an internal $v_z$ rms of $\gtrsim$1000 km/s. None of the six galaxy clusters previously identified as relaxed or plane of sky mergers showed any such excess variance.

We report total abundances and related parameters for the full sample of the FUSE survey of molecular hydrogen in 38 translucent lines of sight. New results are presented for the "second half" of the survey involving 15 lines of sight to supplement data for the first 23 lines of sight already published. We assess the correlations between molecular hydrogen and various extinction parameters in the full sample, which covers a broader range of conditions than the initial sample. In particular, we are now able to confirm that many, but not all, lines of sight with shallow far-UV extinction curves and large values of the total-to-selective extinction ratio, $R_V$ = $A_V$ / $E(B-V)$ -- characteristic of larger than average dust grains -- are associated with particularly low hydrogen molecular fractions ($f_{\rm H2}$). In the lines of sight with large $R_V$, there is in fact a wide range in molecular fractions, despite the expectation that the larger grains should lead to less H$_2$ formation. However, we see specific evidence that the molecular fractions in this sub-sample are inversely related to the estimated strength of the UV radiation field and thus the latter factor is more important in this regime. We have provided an update to previous values of the gas-to-dust ratio, $N$(H$_{\rm tot}$)/$E(B-V)$, based on direct measurements of $N$(H$_2$) and $N$(H I). Although our value is nearly identical to that found with Copernicus data, it extends the relationship by a factor of 2 in reddening. Finally, as the new lines of sight generally show low to moderate molecular fractions, we still find little evidence for single monolithic "translucent clouds" with $f_{\rm H2}$ $\sim$ 1.

We identify an abundant population of extreme emission line galaxies (EELGs) at redshift z~1.7 in the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS) imaging from Hubble Space Telescope/Wide Field Camera 3 (HST/WFC3). 69 EELG candidates are selected by the large contribution of exceptionally bright emission lines to their near-infrared broad-band magnitudes. Supported by spectroscopic confirmation of strong [OIII] emission lines -- with rest-frame equivalent widths ~1000\AA -- in the four candidates that have HST/WFC3 grism observations, we conclude that these objects are galaxies with 10^8 Msol in stellar mass, undergoing an enormous starburst phase with M_*/(dM_*/dt) of only ~15 Myr. These bursts may cause outflows that are strong enough to produce cored dark matter profiles in low-mass galaxies. The individual star formation rates and the co-moving number density (3.7x10^-4 Mpc^-3) can produce in ~4 Gyr much of the stellar mass density that is presently contained in 10^8-10^9 Msol dwarf galaxies. Therefore, our observations provide a strong indication that many or even most of the stars in present-day dwarf galaxies formed in strong, short-lived bursts, mostly at z>1.

Star-forming galaxies are considered to be the leading candidate sources that dominate the cosmic reionization at z>7, and the search for analogs at moderate redshift showing Lyman continuum (LyC) leakage is currently a active line of research. We have observed a star-forming galaxy at z=3.2 with Hubble/WFC3 in the F336W filter, corresponding to the 730-890A rest-frame, and detect LyC emission. This galaxy is very compact and also has large Oxygen ratio [OIII]5007/[OII]3727 (>=10). No nuclear activity is revealed from optical/near-infrared spectroscopy and deep multi-band photometry (including the 6Ms X-ray, Chandra). The measured escape fraction of ionizing radiation spans the range 50-100\%, depending on the IGM attenuation. The LyC emission is detected at S/N=10 with m(F336W)=27.57+/-0.11 and it is spatially unresolved, with effective radius R_e<200pc. Predictions from photoionization and radiative transfer models are in line with the properties reported here, indicating that stellar winds and supernova explosions in a nucleated star-forming region can blow cavities generating density-bounded conditions compatible with optically thin media. Irrespective to the nature of the ionizing radiation, spectral signatures of these sources over the entire electromagnetic spectrum are of central importance for their identification during the epoch of reionization, when the LyC is unobservable. Intriguingly, the Spitzer/IRAC photometric signature of intense rest-frame optical emissions ([OIII]+Hbeta) observed recently at z~7.5-8.5 is similar to what is observed in this galaxy. Only the James Webb Space Telescope will measure optical line ratios at z>7 allowing a direct comparison with lower redshift LyC emitters, as reported here.

Using the phat-ELVIS suite of Milky Way-size halo simulations, we show that subhalo orbital pericenters, $r_{\rm peri}$, correlate with their dark matter halo structural properties. Specifically, at fixed maximum circular velocity, $V_{\rm max}$, subhalos with smaller $r_{\rm peri}$ are more concentrated (have smaller $r_{\rm max}$ values) and have lost more mass, with larger peak circular velocities, $V_{\rm peak}$, prior to infall. These trends provide information that can tighten constraints on the inferred $V_{\rm max}$ and $V_{\rm peak}$ values for known Milky Way satellites. We illustrate this using published pericenter estimates enabled by Gaia for the nine classical Milky Way dwarf spheroidal satellites. The two densest dSph satellites (Draco and Ursa Minor) have relatively small pericenters, and this pushes their inferred $r_{\rm max}$ and $V_{\rm max}$ values lower than they would have been without pericenter information. For Draco, we infer $V_{\rm max} = 23.5 \, \pm 3.3$ km s$^{-1}$ (compared to $27.3 \, \pm 7.1$ km s$^{-1}$ without pericenter information). Such a shift exacerbates the traditional Too Big to Fail problem. Draco's peak circular velocity range prior to infall narrows from $V_{\rm peak} = 21 - 49$km s$^{-1}$without pericenter information to$V_{\rm peak} = 25-37$ km s$^{-1}$ with the constraint. Over the full population of classical dwarf spheroidals, we find no correlation between $V_{\rm peak}$ and stellar mass today, indicative of a high level of stochasticity in galaxy formation at stellar masses below $\sim 10^7$ M$_\odot$. As proper motion measurements for dwarf satellites become more precise, they should enable useful priors on the expected structure and evolution of their host dark matter subhalos.

We investigate the burstiness of star formation and the ionizing efficiency of a large sample of galaxies at 0.7 &lt; z &lt; 1.5 using HST grism spectroscopy and deep ultraviolet (UV) imaging in the GOODS-N and GOODS-S fields. The star formation history (SFH) in these strong emission line low-mass galaxies indicates an elevated star formation rate (SFR) based on the H$\alpha$ emission line at a given stellar mass when compared to the standard main sequence. Moreover, when comparing the H$\alpha$ and UV SFR indicators, we find that an excess in SFR(H$\alpha$) compared to SFR(UV) is preferentially observed in lower-mass galaxies below $10^{9}$ M$\odot$, which are also the highest-EW galaxies. These findings suggest that the burstiness parameters of these strong emission line galaxies may differ from those inferred from hydrodynamical simulations and previous observations. For instance, a larger burstiness duty cycle would explain the observed SFR(H$\alpha$) excess. We also estimate the ionizing photon production efficiency $\xi_{ion}$, finding a median value of Log($\xi_{ion}$/erg$^{-1}$ Hz)$=24.80 \pm 0.26$ when adopting a Galactic dust correction for H$\alpha$ and an SMC one for the stellar component. We observe an increase of $\xi_{ion}$ with redshift, further confirming similar results at higher redshifts. We also find that $\xi_{ion}$ is strongly correlated with EW(H$\alpha$), which provides an approach for deriving $\xi_{ion}$ in early galaxies. Lower-mass, lower-luminosity galaxies have a higher $\xi_{ion}$. Overall, these results provide further support for faint galaxies playing a major role in the reionization of the Universe.

We introduce a novel machine learning framework for estimating the Bayesian posteriors of morphological parameters for arbitrarily large numbers of galaxies. The Galaxy Morphology Posterior Estimation Network (GaMPEN) estimates values and uncertainties for a galaxy's bulge-to-total light ratio ($L_B/L_T$), effective radius ($R_e$), and flux ($F$). To estimate posteriors, GaMPEN uses the Monte Carlo Dropout technique and incorporates the full covariance matrix between the output parameters in its loss function. GaMPEN also uses a Spatial Transformer Network (STN) to automatically crop input galaxy frames to an optimal size before determining their morphology. This will allow it to be applied to new data without prior knowledge of galaxy size. Training and testing GaMPEN on galaxies simulated to match z &lt; 0.25 galaxies in Hyper Suprime-Cam Wide $g$-band images, we demonstrate that GaMPEN achieves typical errors of $0.1$ in $L_B/L_T$, $0.17$ arcsec ($\sim 7\%$) in $R_e$, and $6.3\times10^4$ nJy ($\sim 1\%$) in $F$. GaMPEN's predicted uncertainties are well-calibrated and accurate (&lt;5\% deviation) -- for regions of the parameter space with high residuals, GaMPEN correctly predicts correspondingly large uncertainties. We also demonstrate that we can apply categorical labels (i.e., classifications such as "highly bulge-dominated") to predictions in regions with high residuals and verify that those labels are $\gtrsim 97\%$ accurate. To the best of our knowledge, GaMPEN is the first machine learning framework for determining joint posterior distributions of multiple morphological parameters and is also the first application of an STN to optical imaging in astronomy.

We present results from a large volume simulation of Hydrogen reionization. We combine 3d radiative transfer calculations and an N-body simulation, describing structure formation in the intergalactic medium, to detail the growth of HII regions around high redshift galaxies. Our simulation tracks 1024^3 dark matter particles, in a box of co-moving side length 65.6 Mpc/h. This large volume allows us to accurately characterize the size distribution of HII regions throughout most of the reionization process. At the same time, our simulation resolves many of the small galaxies likely responsible for reionization. It confirms a picture anticipated by analytic models: HII regions grow collectively around highly-clustered sources, and have a well-defined characteristic size, which evolves from a sub-Mpc scale at the beginning of reionization to R > 10 Mpc towards the end. We present a detailed statistical description of our results, and compare them with a numerical scheme based on the analytic model by Furlanetto, Zaldarriaga, and Hernquist. We find that the analytic calculation reproduces the size distribution of HII regions and the 21 cm power spectrum of the radiative transfer simulation remarkably well. The ionization field from the simulation, however, has more small scale structure than the analytic calculation, owing to Poisson scatter in the simulated abundance of galaxies on small scales. We propose and validate a simple scheme to incorporate this scatter into our calculations. Our results suggest that analytic calculations are sufficiently accurate to aid in predicting and interpreting the results of future 21 cm surveys. In particular, our fast numerical scheme is useful for forecasting constraints from future 21 cm surveys, and in constructing mock surveys to test data analysis procedures.

We constrain the number density and evolution of Compton-thick Active Galactic Nuclei (AGN). In the local Universe we use the wide area surveys from the Swift and INTEGRAL satellites, while for high redshifts we explore candidate selections based on a combination of X-ray and mid-IR parameters. We find a significantly lower space density of Compton-thick AGN in the local Universe than expected from published AGN population synthesis models to explain the X-ray background. This can be explained by the numerous degeneracies in the parameters of those models; we use the high-energy surveys described here to remove those degeneracies. We show that only direct observations of CT AGN can currently constrain the number of heavily-obscured supermassive black holes. At high redshift, the inclusion of IR-selected Compton-thick AGN candidates leads to a much higher space density, implying (a) a different (steeper) evolution for these sources compared to less-obscured AGN, (b) that the IR selection includes a large number of interlopers, and/or (c) that there is a large number of reflection-dominated AGN missed in the INTEGRAL and Swift observations. The contribution of CT AGN to the X-ray background is small, ~9%, with a comparable contribution to the total cosmic accretion, unless reflection-dominated CT AGN significantly outnumber transmission-dominated CT AGN, in which case their contribution can be much higher. Using estimates derived here for the accretion luminosity over cosmic time we estimate the local mass density in supermassive black holes and find a good agreement with available constraints for an accretion efficiency of ~10%. Transmission-dominated CT AGN contribute only ~8% to total black hole growth.