The Subaru Strategic Program (SSP) is an ambitious multi-band survey using the Hyper Suprime-Cam (HSC) on the Subaru telescope. The Wide layer of the SSP is both wide and deep, reaching a detection limit of i~26.0 mag. At these depths, it is challenging to achieve accurate, unbiased, and consistent photometry across all five bands. The HSC data are reduced using a pipeline that builds on the prototype pipeline for the Large Synoptic Survey Telescope. We have developed a Python-based, flexible framework to inject synthetic galaxies into real HSC images called SynPipe. Here we explain the design and implementation of SynPipe and generate a sample of synthetic galaxies to examine the photometric performance of the HSC pipeline. For stars, we achieve 1% photometric precision at i~19.0 mag and 6% precision at i~25.0 in the i-band. For synthetic galaxies with single-Sersic profiles, forced CModel photometry achieves 13% photometric precision at i~20.0 mag and 18% precision at i~25.0 in the i-band. We show that both forced PSF and CModel photometry yield unbiased color estimates that are robust to seeing conditions. We identify several caveats that apply to the version of HSC pipeline used for the first public HSC data release (DR1) that need to be taking into consideration. First, the degree to which an object is blended with other objects impacts the overall photometric performance. This is especially true for point sources. Highly blended objects tend to have larger photometric uncertainties, systematically underestimated fluxes and slightly biased colors. Second, >20% of stars at 22.5< i < 25.0 mag can be misclassified as extended objects. Third, the current CModel algorithm tends to strongly underestimate the half-light radius and ellipticity of galaxy with i>21.5 mag.

We investigate two-dimensional image decomposition of nearby, morphologically selected early-type galaxies (ETGs). We are motivated by recent observational evidence of significant size growth of quiescent galaxies and theoretical development advocating a two-phase formation scenario for ETGs. We find that a significant fraction of nearby ETGs show changes in isophotal shape that require multi-component models. The characteristic sizes of the inner and outer component are $\sim 3$ and $\sim 15$ kpc. The inner component lies on the mass-size relation of ETGs at $z \sim 0.25-0.75$, while the outer component tends to be more elliptical and hints at a stochastic buildup process. We find real physical differences between the single- and double-component ETGs, with the double-component galaxies being younger and more metal-rich. The fraction of double component ETGs increases with increasing $\sigma$ and decreases in denser environments. We hypothesize that double-component systems were able to accrete gas and small galaxies until later times, boosting their central densities, building up their outer parts, and lowering their typical central ages. In contrast, the oldest galaxies, perhaps due to residing in richer environments, have no remaining hints of their last accretion episode.

Luminous Red Galaxies (LRG) from the Sloan Digital Sky Survey are considered among the best understood samples of galaxies, and they are employed in a broad range of cosmological studies. Because they form a relatively homogeneous population, with high stellar masses and red colors, they are expected to occupy halos in a relatively simple way. In this paper, we study how LRGs occupy massive halos via direct counts in clusters and we reveal several unexpected trends suggesting that the connection between LRGs and dark matter halos may not be straightforward. Using the redMaPPer cluster catalog, we derive the central occupation of LRGs as a function richness, Ncen({\lambda}). Assuming no correlation between cluster mass and central galaxy luminosity at fixed richness, we show that clusters contain a significantly lower fraction of central LRGs than predicted from the two-point correlation function. At halo masses of 10^14.5 Msun, we find Ncen=0.73, compared to Ncen of 0.89 from correlation studies. Our central occupation function for LRGs converges to 0.95 at large halo masses. A strong anti-correlation between central luminosity and cluster mass at fixed richness is required to reconcile our results with those based on clustering studies. We also derive P_BNC, the probability that the brightest cluster member is not the central galaxy. We find P_BNC ~ 20-30% which is a factor of ~2 lower than the value found by Skibba et al. 2011. Finally, we study the radial offsets of bright non-central LRGs from cluster centers and show that bright non-central LRGs follow a different radial distribution compared to red cluster members, which follow a Navarro-Frank-White profile. This work demonstrates that even the most massive clusters do not always have an LRG at the center, and that the brightest galaxy in a cluster is not always the central galaxy.

Upcoming weak lensing surveys will survey large cosmological volumes to measure the growth of cosmological structure with time and thereby constrain dark energy. One major systematic uncertainty in this process is the calibration of the weak lensing shape distortions, or shears. Most upcoming surveys plan to test several aspects of their shear estimation algorithms using sophisticated image simulations that include realistic galaxy populations based on high-resolution data from the Hubble Space Telescope (HST). However, existing datasets from the (HST) cover very small cosmological volumes, so cosmic variance could cause the galaxy populations in them to be atypical. A narrow redshift slice from such surveys could be dominated by a single large overdensity or underdensity. In that case, the morphology-density relation could alter the local galaxy populations and yield an incorrect calibration of shear estimates as a function of redshift. We directly test this scenario using the COSMOS survey, the largest-area (HST) survey to date, and show how the statistical distributions of galaxy shapes and morphological parameters (e.g., S\'{e}rsic $n$) are influenced by redshift-dependent cosmic variance. The typical variation in RMS ellipticity due to environmental effects is 5 per cent (absolute, not relative) for redshift bins of width $\Delta z=0.05$, which could result in uncertain shear calibration at the 1 per cent level. We conclude that the cosmic variance effects are large enough to exceed the systematic error budget of future surveys, but can be mitigated with careful choice of training dataset and sufficiently large redshift binning.

Accretion is thought to primarily contribute to the mass accumulation history of supermassive black holes throughout cosmic time. While this may be true at high redshifts, at lower redshifts and for the most massive black holes mergers themselves might add significantly to the mass budget. We evolve SMBHs from $4 > z > 0$ using merger trees derived from hydrodynamical cosmological simulations of a cluster and void region, scaled to the observed value of the stellar mass fraction to account for overcooling. Mass gains from gas accretion proportional to bulge growth and BH-BH mergers are tracked, as are black holes that remain "orbiting" due to insufficient dynamical friction in a merger remnant, as well as those that are ejected due to gravitational recoil. We find that gas accretion remains the dominant source of mass accumulation in almost all SMBHs; mergers contribute $2.5\pm0.1\%$ for all SMBHs in the cluster and $1.0\pm0.1\%$ in the void since $z = 4$. However, mergers are significant for massive SMBHs. The fraction of mass accumulated from mergers for central BHs generally increases for larger values of the host bulge mass: in the void, the fraction is $2\%$ at $M_{*, bul} = 10^{10} M_{\odot}$, increasing to $4\%$ at $M_{*, bul} \gtrsim 10^{11} M_{\odot}$, and in the cluster it is $4\%$ at $M_{*, bul} = 10^{10} M_{\odot}$ and $23\%$ at $10^{12} M_{\odot}$. We find that $40\%$ of SMBHs and $\approx 8\%$ of the total SMBH mass is found orbiting in the cluster region at $z = 0$. The existence of orbiting and ejected SMBHs requires modification of the Soltan argument. We estimate this correction to the integrated accreted mass density of SMBHs to be in the range $6-21\%$, with a mean value of $11\pm3\%$. We also calculate the total energy output and strain from gravitational waves emitted by merging SMBHs, and obtain a signal potentially detectable by pulsar timing arrays.

Central Galaxies (CGs) in massive halos live in unique environments with formation histories closely linked to that of the host halo. In local clusters they have larger sizes ($R_e$) and lower velocity dispersions (sigma) at fixed stellar mass M_star, and much larger R_e at a fixed $\sigma$ than field and satellite galaxies (non-CGs). Using spectroscopic observations of group galaxies selected from the COSMOS survey, we compare the dynamical scaling relations of early-type CGs and non-CGs at z~0.6, to distinguish possible mechanisms that produce the required evolution. CGs are systematically offset towards larger R_e at fixed $\sigma$ compared to non-CGs with similar M_star. The CG R_e-M_star relation also shows differences, primarily driven by a sub-population (~15%) of galaxies with large $R_e$, while the M_star-sigma relations are indistinguishable. These results are accentuated when double Sersic profiles, which better fit light in the outer regions of galaxies, are adopted. They suggest that even group-scale CGs can develop extended components by these redshifts that can increase total $R_e$ and M_star estimates by factors of ~2. To probe the evolutionary link between our sample and cluster CGs, we also analyze two cluster samples at z~0.6 and z~0. We find similar results for the more massive halos at comparable z, but much more distinct CG scaling relations at low-z. Thus, the rapid, late-time accretion of outer components, perhaps via the stripping and accretion of satellites, would appear to be a key feature that distinguishes the evolutionary history of CGs.

The role of major mergers in galaxy and black hole formation is not well constrained. To help address this, we develop an automated method to identify late-stage galaxy mergers before coalescence of the galactic cores. The resulting sample of mergers is distinct from those obtained using pair-finding and morphological indicators. Our method relies on median-filtering of high-resolution images in order to distinguish two concentrated galaxy nuclei at small separations. Using mock images, we derive statistical contamination and incompleteness corrections for the fraction of late-stage mergers. We apply our new method to a magnitude-limited (I < 23) sample of 44,164 galaxies from the COSMOS HST/ACS catalog. Using a mass-complete sample with $\log M_*/M_\odot > 10.6$ and $0.25 < z \leq 1.00$, we find ~5% of systems are late-stage mergers with separations between 2.2 and 8 kpc. Correcting for incompleteness and contamination, the fractional merger rate increases strongly with redshift as $(1+z)^{3.8\pm0.9}$, in agreement with earlier studies and with dark matter halo merger rates. Separating the sample into star-forming and quiescent galaxies shows that the merger rate for star-forming galaxies increases strongly redshift, $(1+z)^{4.5\pm1.3}$, while the merger rate for quiescent galaxies is consistent with no evolution, $(1+z)^{1.1\pm1.2}$. Limiting our sample to galaxies with spectroscopic redshifts from zCOSMOS, we find that the star formation rates and X-ray selected AGN activity in likely late-stage mergers are enhanced by factors of ~2 relative to a control sample. Combining our sample with more widely separated pairs, we find that $8\pm5\%$ of star formation and $20\pm8\%$ of AGN activity is triggered by close encounters (<143 kpc) or mergers, once more suggesting that major mergers are not the only channels for star formation and black hole growth. (abridged)

The GRavitational lEnsing Accuracy Testing 3 (GREAT3) challenge is the third in a series of image analysis challenges, with a goal of testing and facilitating the development of methods for analyzing astronomical images that will be used to measure weak gravitational lensing. This measurement requires extremely precise estimation of very small galaxy shape distortions, in the presence of far larger intrinsic galaxy shapes and distortions due to the blurring kernel caused by the atmosphere, telescope optics, and instrumental effects. The GREAT3 challenge is posed to the astronomy, machine learning, and statistics communities, and includes tests of three specific effects that are of immediate relevance to upcoming weak lensing surveys, two of which have never been tested in a community challenge before. These effects include realistically complex galaxy models based on high-resolution imaging from space; spatially varying, physically-motivated blurring kernel; and combination of multiple different exposures. To facilitate entry by people new to the field, and for use as a diagnostic tool, the simulation software for the challenge is publicly available, though the exact parameters used for the challenge are blinded. Sample scripts to analyze the challenge data using existing methods will also be provided. See http://great3challenge.info and http://great3.projects.phys.ucl.ac.uk/leaderboard/ for more information.

We examine the changes in the properties of galactic bulges and discs with environment for a volume-limited sample of 12500 nearby galaxies from SDSS. We focus on galaxies with classical bulges. Classical bulges seem to have the same formation history as ellipticals of the same mass, and we test if environment determines whether or not a classical bulge possesses a disc. Using the projected fifth nearest neighbour density as a measure of local environment, we look for correlations with environment at fixed bulge stellar mass. In groups with fewer than 20 members, we find no evidence for changes in disc morphology with local density. At fixed bulge mass, disc mass and disc scale length are independent of local density. However, disc colour does increase (Delta(g - r) ~ 0.05 mag) as a function of local density in relatively poor groups. Therefore, the colour--density relation for classical bulge+disc galaxies in the field and in poor groups is due solely to changes in disc colour with density. In contrast, we find no correlations between disc colour and local density for classical bulge+disc galaxies in large, relaxed groups and clusters. However, there is a weak correlation between disc mass and group crossing time, suggesting morphological transformation takes places in rich groups. Our results add to the evidence that star formation is quenched in group environments, instead of clusters, and that star formation quenching and morphological transformation are separate processes. Overall, we show that environment has two effects on galactic discs: relatively low density environments can quench star formation in discs, while processes occurring in higher density environments contribute to the morphological transformation from disc-dominated systems to bulge-dominated systems.

We examine galaxy formation in a cosmological AMR simulation, which includes two high resolution boxes, one centered on a 3 \times 10^14 M\odot cluster, and one centered on a void. We examine the evolution of 611 massive (M\ast > 10^10M\odot) galaxies. We find that the fraction of the final stellar mass which is accreted from other galaxies is between 15 and 40% and increases with stellar mass. The accreted fraction does not depend strongly on environment at a given stellar mass, but the galaxies in groups and cluster environments are older and underwent mergers earlier than galaxies in lower density environments. On average, the accreted stars are ~2.5 Gyrs older, and ~0.15 dex more metal poor than the stars formed in-situ. Accreted stellar material typically lies on the outskirts of galaxies; the average half-light radius of the accreted stars is 2.6 times larger than that of the in-situ stars. This leads to radial gradients in age and metallicity for massive galaxies, in qualitative agreement with observations. Massive galaxies grow by mergers at a rate of approximately 2.6% per Gyr. These mergers have a median (mass-weighted) mass ratio less than 0.26 \pm 0.21, with an absolute lower limit of 0.20, for galaxies with M\ast ~ 10^12 M\odot. This suggests that major mergers do not dominate in the accretion history of massive galaxies. All of these results agree qualitatively with results from SPH simulations by Oser et al. (2010, 2012).

We present a set of bulge-disk decompositions for a sample of 71,825 SDSS main-sample galaxies in the redshift range 0.003<z<0.05. We have fit each galaxy with either a de Vaucouleurs ('classical') or an exponential ('pseudo-') bulge and an exponential disk. Two dimensional Sersic fits are performed when the 2-component fits are not statistically significant or when the fits are poor, even in the presence of high signal-to-noise. We study the robustness of our 2-component fits by studying a bright subsample of galaxies and we study the systematics of these fits with decreasing resolution and S/N. Only 30% of our sample have been fit with two-component fits in which both components are non-zero. The g-r and g-i colours of each component for the two-component models are determined using linear templates derived from the r-band model. We attempt a physical classification of types of fits into disk galaxies, pseudo-bulges, classical bulges, and ellipticals. Our classification of galaxies agrees well with previous large B+D decomposed samples. Using our galaxy classifications, we find that Petrosian concentration is a good indicator of B/T, while overall Sersic index is not. Additionally, we find that the majority of green valley galaxies are bulge+disk galaxies. Furthermore, in the transition from green to red B+D galaxies, the total galaxy colour is most strongly correlated with the disk colour.

In this paper, we derive scaling relations between photometric observable quantities and disk galaxy rotation velocity V_rot, or Tully-Fisher relations (TFRs). Our methodology is dictated by our purpose of obtaining purely photometric, minimal-scatter estimators of V_rot applicable to large galaxy samples from imaging surveys. To achieve this goal, we have constructed a sample of 189 disk galaxies at redshifts z<0.1 with long-slit H-alpha spectroscopy from Pizagno et al. (2007) and new observations. By construction, this sample is a fair subsample of a large, well-defined parent disk sample of ~170 000 galaxies selected from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7). The optimal photometric estimator of V_rot we find is stellar mass M_* from Bell et al. (2003), based on the linear combination of a luminosity and a colour. Assuming a Kroupa IMF, we find: log [V_{80}/(km s^-1)] = (2.142 +/- 0.004)+(0.278 +/- 0.010)[log (M_*/M_sun)-10.10], where V_{80} is the rotation velocity measured at the radius R_{80} containing 80 per cent of the i-band galaxy light. This relation has an intrinsic Gaussian scatter of 0.036 +/- 0.005 dex and a measured scatter of 0.056 dex in log V_{80}. For a fixed IMF, we find that the dynamical-to-stellar mass ratios within R_{80}, (M_dyn/M_*)(R_{80}), decrease from approximately 10 to 3, as stellar mass increases from M_* ~ 10^9 to 10^{11} M_sun. At a fixed stellar mass, (M_dyn/M_*)(R_{80}) increases with disk size, so that it correlates more tightly with stellar surface density than with stellar mass or disk size alone. In future work, we will use these results to study disk galaxy formation and evolution, and perform a fair, statistical analysis of the dynamics and masses of a photometrically-selected sample of disk galaxies. [Abridged]

We examine two extreme models for the build-up of the stellar component of luminous elliptical galaxies. In one case, we assume the build-up of stars is dissipational, with centrally accreted gas radiating away its orbital and thermal energy; the dark matter halo will undergo adiabatic contraction and the central dark matter density profile will steepen. For the second model, we assume the central galaxy is assembled by a series of dissipationless mergers of stellar clumps that have formed far from the nascent galaxy. In order to be accreted, these clumps lose their orbital energy to the dark matter halo via dynamical friction, thereby heating the central dark matter and smoothing the dark matter density cusp. The central dark matter density profiles differ drastically between these models. For the isolated elliptical galaxy, NGC 4494, the central dark matter densities follow the power-laws r^(-0.2) and r^(-1.7) for the dissipational and dissipationless models, respectively. By matching the dissipational and dissipationless models to observations of the stellar component of elliptical galaxies, we examine the relative contributions of dissipational and dissipationless mergers to the formation of elliptical galaxies and look for observational tests that will distinguish between these models. Comparisons to strong lensing brightest cluster galaxies yield median M*/L_B ratios of 2.1+/-0.8 and 5.2+/-1.7 at z=0.39 for the dissipational and dissipationless models, respectively. For NGC 4494, the best-fit dissipational and dissipationless models have M*/L_B=2.97 and 3.96. Comparisons to expected stellar mass-to-light ratios from passive evolution and population syntheses appear to rule out a purely dissipational formation mechanism for the central stellar regions of giant elliptical galaxies.

Tidal dissipation may be important for the internal evolution as well as the orbits of short-period massive planets--hot Jupiters. We revisit a mechanism proposed by Ogilvie and Lin for tidal forcing of inertial waves, which are short-wavelength, low-frequency disturbances restored primarily by Coriolis rather than buoyancy forces. This mechanism is of particular interest for hot Jupiters because it relies upon a rocky core, and because these bodies are otherwise largely convective. Compared to waves excited at the base of the stratified, externally heated atmosphere, waves excited at the core are more likely to deposit heat in the convective region and thereby affect the planetary radius. However, Ogilvie and Lin's results were numerical, and the manner of the wave excitation was not clear. Using WKB methods, we demonstrate the production of short waves by scattering of the equilibrium tide off the core at critical latitudes. The tidal dissipation rate associated with these waves scales as the fifth power of the core radius, and the implied tidal $Q$ is of order ten million for nominal values of the planet's mass, radius, orbital period, and core size. We comment upon an alternative proposal by Wu for exciting inertial waves in an unstratified fluid body by means of compressibility rather than a core. We also find that even a core of rock is unlikely to be rigid. But Ogilvie and Lin's mechanism should still operate if the core is substantially denser than its immediate surroundings.