We present the first stellar velocity dispersion measurement of a massive quenching galaxy at z=4.01. The galaxy is first identified as a massive z>~4 galaxy with suppressed star formation from photometric redshifts based on deep multi-band data in the UKIDSS Ultra Deep Survey field. A follow-up spectroscopic observation with MOSFIRE on Keck revealed strong multiple absorption features, which are identified as Balmer absorption lines, giving a secure redshift of z=4.01. Thanks to the high S/N of the spectrum, we are able to estimate the stellar velocity dispersion, sigma=268+/-59 km/s. This velocity dispersion is consistent with that of massive galaxies today, implying no significant evolution in stellar velocity dispersion over the last 12 Gyr. Based on an upper limit on its physical size from deep optical images (r_eff<1.3 kpc), we find that its dynamical mass is consistent with the stellar mass inferred from photometry. Furthermore, the galaxy is located on the mass fundamental plane extrapolated from lower redshift galaxies. Combining all these results, we find that the velocity dispersion does not significantly evolve with redshift, although the size and mass of massive quenched galaxies do. This suggests that the mass in the core of massive galaxies does not evolve significantly, while most of the mass growth occurs in the outskirts of the galaxies, which also increases the size. This picture is consistent with a two-phase formation scenario in which mass and size growth is due to accretion in the outskirts of galaxies via mergers.

The duration of a gamma-ray burst (GRB) is a key indicator of its physics origin, with long bursts perhaps associated with the collapse of massive stars and short bursts with mergers of neutron stars.However, there is substantial overlap in the properties of both short and long GRBs and neither duration nor any other parameter so far considered completely separates the two groups. Here we unambiguously classify every GRB using a machine-learning, dimensionality-reduction algorithm, t-distributed stochastic neighborhood embedding (t-SNE), providing a catalog separating all Swift GRBs into two groups. Although the classification takes place only using prompt emission light curves,every burst with an associated supernova is found in the longer group and bursts with kilonovae in the short, suggesting along with the duration distributions that these two groups are truly long and short GRBs. Two bursts with a clear absence of a supernova belong to the longer class, indicating that these might have been direct-collapse black holes, a proposed phenomenon that may occur in the deaths of more massive stars.

Redshifts used in current cosmological supernova samples are measured using two primary techniques, one based on well-measured host galaxy spectral lines and the other based on supernova-dominated spectra. Here, we construct an updated Pantheon catalog with revised redshifts, redshift sources and estimated uncertainties for the entire sample to investigate whether these two techniques yield consistent results. The best-fit cosmological parameters using these two measurement techniques disagree, and we explore several possible sources of bias which could result from using the lower-precision supernova-dominated redshifts. In a pilot study, we show that using a host redshift-only subsample of the Pantheon catalog will result in lower $\Omega_m$ and matter density $\Omega_m h^2$ and slightly higher $H_0$ than previous analysis which, among other possibilities, could result in supernova and CMB measurements agreeing on $\Omega_m h^2$ despite tension in $H_0$. To obtain rigorous results, though, the Pantheon catalog should be improved by obtaining host spectra for supernova that have faded and future surveys should be designed to use host galaxy redshifts rather than lower-precision methods.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

Large photometric surveys provide a rich source of observations of quiescent galaxies, including a surprisingly large population at z>1. However, identifying large, but clean, samples of quiescent galaxies has proven difficult because of their near-degeneracy with interlopers such as dusty, star-forming galaxies. We describe a new technique for selecting quiescent galaxies based upon t-distributed stochastic neighbor embedding (t-SNE), an unsupervised machine learning algorithm for dimensionality reduction. This t-SNE selection provides an improvement both over UVJ, removing interlopers which otherwise would pass color selection, and over photometric template fitting, more strongly towards high redshift. Due to the similarity between the colors of high- and low-redshift quiescent galaxies, under our assumptions t-SNE outperforms template fitting in 63% of trials at redshifts where a large training sample already exists. It also may be able to select quiescent galaxies more efficiently at higher redshifts than the training sample.

We report two secure ($z=3.775, 4.012$) and one tentative ($z\approx3.767$) spectroscopic confirmations of massive and quiescent galaxies through $K$-band observations with Keck/MOSFIRE and VLT/X-Shooter. The stellar continuum emission, the absence of strong nebular emission lines and the lack of significant far-infrared detections confirm the passive nature of these objects, disfavoring the alternative solution of low-redshift dusty star-forming interlopers. We derive stellar masses of $\mathrm{log}(M_{\star}/M_\odot)\sim11$ and ongoing star formation rates placing these galaxies $\gtrsim 1-2$ dex below the main sequence at their redshifts. The adopted parametrization of the star formation history suggests that these sources experienced a strong ($\langle \rm SFR \rangle \sim 1200-3500\,M_\odot\,\mathrm{yr}^{-1}$) and short ($\sim 50$ Myr) burst of star formation, peaking $\sim 150-500$ Myr before the time of observation, all properties reminiscent of the characteristics of sub-millimeter galaxies (SMGs) at $z>4$. We investigate this connection by comparing the comoving number densities and the properties of these two populations. We find a fair agreement only with the deepest sub-mm surveys detecting not only the most extreme starbursts, but also more normal galaxies. We support these findings by further exploring the Illustris-TNG cosmological simulation, retrieving populations of both fully quenched massive galaxies at $z\sim3-4$ and SMGs at $z\sim4-5$, with number densities and properties in agreement with the observations at $z\sim3$, but in increasing tension at higher redshift. Nevertheless, as suggested by the observations, not all the progenitors of quiescent galaxies at these redshifts shine as bright SMGs in their past and, similarly, not all bright SMGs quench by $z\sim3$, both fractions depending on the threshold assumed to define the SMGs themselves.

In the context of the BUFFALO (Beyond Ultra-deep Frontier Fields And Legacy Observations) survey, we present a new analysis of the merging galaxy cluster MACS\,J0416.1-2403 ($z = 0.397$) and its parallel field using the data collected by the Hubble Frontier Fields (HFF) campaign. In this work, we measure the surface mass density from a weak-lensing analysis, and characterise the overall matter distribution in both the cluster and parallel fields. The surface mass distribution derived for the parallel field shows clumpy overdensities connected by filament-like structures elongated in the direction of the cluster core. We also characterise the X-ray emission of the cluster, and compare it with the lensing mass distribution. We identify five substructures at the $>5\sigma$ level over the two fields, four of them being in the cluster one. Furthermore, three of them are located close to the edges of the field of view, and border issues can significantly hamper the determination of their physical parameters. Finally, we compare our results with the predicted subhalo distribution of one of the Hydrangea/C-EAGLE simulated cluster. Significant differences are obtained suggesting the simulated cluster is at a more advanced evolutionary state than MACS\,J0416.1-2403. Our results anticipate the upcoming BUFFALO observations that will link the two HFF fields, extending further the \emph{HST} coverage, and thus allowing a better characterisation of the reported substructures.

We present a detailed analysis of a large sample of spectroscopically confirmed ultra-massive quiescent galaxies (${\rm{log}}(M_{\ast}/M_{\odot})\sim11.5$) at $z\gtrsim2$. This sample comprises 15 galaxies selected in the COSMOS and UDS fields by their bright K-band magnitudes and followed up with VLT/X-shooter spectroscopy and HST/WFC3 $H_{F160W}$ imaging. These observations allow us to unambiguously confirm their redshifts ascertain their quiescent nature and stellar ages, and to reliably assess their internal kinematics and effective radii. We find that these galaxies are compact, consistent with the high mass end of the mass-size relation for quiescent galaxies at $z=2$. Moreover, the distribution of the measured stellar velocity dispersions of the sample is consistent with the most massive local early-type galaxies from the MASSIVE Survey showing that evolution in these galaxies, is dominated by changes in size. The HST images reveal, as surprisingly high, that $40\ \%$ of the sample have tidal features suggestive of mergers and companions in close proximity, including three galaxies experiencing ongoing major mergers. The absence of velocity dispersion evolution from $z=2$ to $0$, coupled with a doubling of the stellar mass, with a factor of four size increase and the observed disturbed stellar morphologies support dry minor mergers as the primary drivers of the evolution of the massive quiescent galaxies over the last 10 billion years.

We argue that the interplay between cosmic rays, the initial mass function, and star formation plays a crucial role in regulating the star-forming "main sequence". To explore these phenomena we develop a toy model for galaxy evolution in which star formation is regulated by a combination of a temperature-dependent initial mass function and heating due to starlight, cosmic rays and, at very high redshift, the cosmic microwave background. This produces an attractor, near-equilibrium solution which is consistent with observations of the star-forming main sequence over a broad redshift range. Additional solutions to the same equations may correspond to other observed phases of galaxy evolution including quiescent galaxies. This model makes several falsifiable predictions, including higher metallicities and dust masses than anticipated at high redshift and isotopic abundances in the Milky Way. It also predicts that stellar mass-to-light ratios are lower than produced using a Milky Way-derived IMF, so that inferences of stellar masses and star formation rates for high redshift galaxies are overestimated. In some cases, this may also transform inferred dark matter profiles from core-like to cusp-like.

Studying the formation and evolution of galaxies at the earliest cosmic times, and their role in reionization, requires the deepest imaging possible. Ultra-deep surveys like the HUDF and HFF have pushed to mag \mAB$\,\sim\,$30, revealing galaxies at the faint end of the LF to $z$$\,\sim\,$9$\,-\,$11 and constraining their role in reionization. However, a key limitation of these fields is their size, only a few arcminutes (less than a Mpc at these redshifts), too small to probe large-scale environments or clustering properties of these galaxies, crucial for advancing our understanding of reionization. Achieving HUDF-quality depth over areas $\sim$100 times larger becomes possible with a mission like the Wide Field Infrared Survey Telescope (WFIRST), a 2.4-m telescope with similar optical properties to HST, with a field of view of $\sim$1000 arcmin$^2$, $\sim$100$\times$ the area of the HST/ACS HUDF. This whitepaper motivates an Ultra-Deep Field survey with WFIRST, covering $\sim$100$\,-\,$300$\times$ the area of the HUDF, or up to $\sim$1 deg$^2$, to \mAB$\,\sim\,$30, potentially revealing thousands of galaxies and AGN at the faint end of the LF, at or beyond $z$\,$\sim$\,9$\,-\,$10 in the epoch of reionization, and tracing their LSS environments, dramatically increasing the discovery potential at these redshifts. (Note: This paper is a somewhat expanded version of one that was submitted as input to the Astro2020 Decadal Survey, with this version including an Appendix (which exceeded the Astro2020 page limits), describing how the science drivers for a WFIRST Ultra Deep Field might map into a notional observing program, including the filters used and exposure times needed to achieve these depths.)

We argue that an increased temperature in star-forming clouds alters the stellar initial mass function to be more bottom-light than in the Milky Way. At redshifts $z \gtrsim 6$, heating from the cosmic microwave background radiation produces this effect in all galaxies, and it is also present at lower redshifts in galaxies with very high star formation rates (SFRs). A failure to account for it means that at present, photometric template fitting likely overestimates stellar masses and star formation rates for the highest-redshift and highest-SFR galaxies. In addition this may resolve several outstanding problems in the chemical evolution of galactic halos.

We present a catalog of 10718 objects in the COSMOS field observed through multi-slit spectroscopy with the Deep Imaging Multi-Object Spectrograph (DEIMOS) on the Keck II telescope in the wavelength range ~5500-9800A. The catalog contains 6617 objects with high-quality spectra (two or more spectral features), and 1798 objects with a single spectroscopic feature confirmed by the photometric redshift. For 2024 typically faint objects we could not obtain reliable redshifts. The objects have been selected from a variety of input catalogs based on multi-wavelength observations in the field, and thus have a diverse selection function, which enables the study of the diversity in the galaxy population. The magnitude distribution of our objects is peaked at I_AB~23 and K_AB~21, with a secondary peak at K_AB~24. We sample a broad redshift distribution in the range 0<z<6, with one peak at z~1, and another one around z~4. We have identified 13 redshift spikes at z>0.65 with chance probabilities <4xE-4$, some of which are clearly related to protocluster structures of sizes >10 Mpc. An object-to-object comparison with a multitude of other spectroscopic samples in the same field shows that our DEIMOS sample is among the best in terms of fraction of spectroscopic failures and relative redshift accuracy. We have determined the fraction of spectroscopic blends to about 0.8% in our sample. This is likely a lower limit and at any rate well below the most pessimistic expectations. Interestingly, we find evidence for strong lensing of Ly-alpha background emitters within the slits of 12 of our target galaxies, increasing their apparent density by about a factor of 4.

Dust-enshrouded, starbursting, submillimeter galaxies (SMGs) at $z \geq 3$ have been proposed as progenitors of $z \geq 2$ compact quiescent galaxies (cQGs). To test this connection, we present a detailed spatially resolved study of the stars, dust and stellar mass in a sample of six submillimeter-bright starburst galaxies at $z \sim 4.5$. The stellar UV emission probed by HST is extended, irregular and shows evidence of multiple components. Informed by HST, we deblend Spitzer/IRAC data at rest-frame optical finding that the systems are undergoing minor mergers, with a typical stellar mass ratio of 1:6.5. The FIR dust continuum emission traced by ALMA locates the bulk of star formation in extremely compact regions (median $r_{\rm{e}} = 0.70 \pm 0.29$ kpc) and it is in all cases associated with the most massive component of the mergers (median $\log (M_{*}/M_{\odot}) = 10.49 \pm 0.32$). We compare spatially resolved UV slope ($\beta$) maps with the FIR dust continuum to study the infrared excess ($\rm{IRX} = L_{\rm{IR}}/L_{\rm{UV}}$)-$\beta$ relation. The SMGs display systematically higher $\rm{IRX}$ values than expected from the nominal trend, demonstrating that the FIR and UV emissions are spatially disconnected. Finally, we show that the SMGs fall on the mass-size plane at smaller stellar masses and sizes than cQGs at $z = 2$. Taking into account the expected evolution in stellar mass and size between $z = 4.5$ and $z = 2$ due to the ongoing starburst and mergers with minor companions, this is in agreement with a direct evolutionary connection between the two populations.

Telescopes are much more expensive than astronomers, so it is essential to minimize required sample sizes by using the most data-efficient statistical methods possible. However, the most commonly used model-independent techniques for finding the relationship between two variables in astronomy are flawed. In the worst case they can lead without warning to subtly yet catastrophically wrong results, and even in the best case they require more data than necessary. Unfortunately, there is no single best technique for nonparametric regression. Instead, we provide a guide for how astronomers can choose the best method for their specific problem and provide a python library with both wrappers for the most useful existing algorithms and implementations of two new algorithms developed here.

At redshift z = 2, when the Universe was just three billion years old, half of the most massive galaxies were extremely compact and had already exhausted their fuel for star formation(1-4). It is believed that they were formed in intense nuclear starbursts and that they ultimately grew into the most massive local elliptical galaxies seen today, through mergers with minor companions(5,6), but validating this picture requires higher-resolution observations of their centres than is currently possible. Magnification from gravitational lensing offers an opportunity to resolve the inner regions of galaxies(7). Here we report an analysis of the stellar populations and kinematics of a lensed z = 2.1478 compact galaxy, which surprisingly turns out to be a fast-spinning, rotationally supported disk galaxy. Its stars must have formed in a disk, rather than in a merger-driven nuclear starburst(8). The galaxy was probably fed by streams of cold gas, which were able to penetrate the hot halo gas until they were cut off by shock heating from the dark matter halo(9). This result confirms previous indirect indications(10-13) that the first galaxies to cease star formation must have gone through major changes not just in their structure, but also in their kinematics, to evolve into present-day elliptical galaxies.

We measure the stellar mass function (SMF) of galaxies in the COSMOS field up to $z\sim6$. We select them in the near-IR bands of the COSMOS2015 catalogue, which includes ultra-deep photometry from UltraVISTA-DR2, SPLASH, and Subaru/Hyper-SuprimeCam. At $z>2.5$ we use new precise photometric redshifts with error $\sigma_z=0.03(1+z)$ and an outlier fraction of $12\%$, estimated by means of the unique spectroscopic sample of COSMOS. The increased exposure time in the DR2, along with our panchromatic detection strategy, allow us to improve the stellar mass completeness at high $z$ with respect to previous UltraVISTA catalogues. We also identify passive galaxies through a robust colour-colour selection, extending their SMF estimate up to $z=4$. Our work provides a comprehensive view of galaxy stellar mass assembly between $z=0.1$ and 6, for the first time using consistent estimates across the entire redshift range. We fit these measurements with a Schechter function, correcting for Eddington bias. We compare the SMF fit with the halo mass function predicted from $\Lambda$CDM simulations. We find that at $z>3$ both functions decline with a similar slope in the high-mass end. This feature could be explained assuming that the mechanisms that quench star formation in massive haloes become less effective at high redshift; however further work needs to be done to confirm this scenario. Concerning the SMF low-mass end, it shows a progressive steepening as moving towards higher redshifts, with $\alpha$ decreasing from $-1.47_{-0.02}^{+0.02}$ at $z\simeq0.1$ to $-2.11_{-0.13}^{+0.30}$ at $z\simeq5$. This slope depends on the characterisation of the observational uncertainties, which is crucial to properly remove the Eddington bias. We show that there is currently no consensus on the method to quantify such errors: different error models result in different best-fit Schechter parameters. [Abridged]

We combine star formation along the `main sequence', quiescence, and clustering and merging to produce an empirical model for the evolution of individual galaxies. Main sequence star formation alone would significantly steepen the stellar mass function towards low redshift, in sharp conflict with observation. However, a combination of star formation and merging produces a consistent result for correct choice of the merger rate function. As a result, we are motivated to propose a model in which hierarchical merging is disconnected from environmentally-independent star formation. This model can be tested via correlation functions and would produce new constraints on clustering and merging.

The current hierarchical merging paradigm and $\Lambda$CDM predict that the $z \sim 4-8$ universe should be a time in which the most massive galaxies are transitioning from their initial halo assembly to the later baryonic evolution seen in star-forming galaxies and quasars. However, no evidence of this transition has been found in many high redshift galaxy surveys including CFHTLS, CANDELS and SPLASH, the first studies to probe the high-mass end at these redshifts. Indeed, if halo mass to stellar mass ratios estimated at lower-redshift continue to $z \sim 6-8$, CANDELS and SPLASH report several orders of magnitude more $M \sim 10^{12-13} M_\odot$ halos than are possible to have formed by those redshifts, implying these massive galaxies formed impossibly early. We consider various systematics in the stellar synthesis models used to estimate physical parameters and possible galaxy formation scenarios in an effort to reconcile observation with theory. Although known uncertainties can greatly reduce the disparity between recent observations and cold dark matter merger simulations, even taking the most conservative view of the observations, there remains considerable tension with current theory.

We measure the H{\alpha} and [OIII] emission line properties as well as specific star-formation rates (sSFR) of spectroscopically confirmed 3<z<6 galaxies in COSMOS from their observed colors vs. redshift evolution. Our model describes consistently the ensemble of galaxies including intrinsic properties (age, metallicity, star-formation history), dust-attenuation, and optical emission lines. We forward-model the measured H{\alpha} equivalent-widths (EW) to obtain the sSFR out to z~6 without stellar mass fitting. We find a strongly increasing rest-frame H{\alpha} EW that is flattening off above z~2.5 with average EWs of 300-600A at z~6. The sSFR is increasing proportional to (1+z)^2.4 at z<2.2 and (1+z)^1.5 at higher redshifts, indicative of a fast mass build-up in high-z galaxies within e-folding times of 100-200Myr at z~6. The redshift evolution at z>3 cannot be fully explained in a picture of cold accretion driven growth. We find a progressively increasing [OIII]{\lambda}5007/H{\beta} ratio out to z~6, consistent with the ratios in local galaxies selected by increasing H{\alpha} EW (i.e., sSFR). This demonstrates the potential of using "local high-z analogs" to investigate the spectroscopic properties and relations of galaxies in the re-ionization epoch.

We investigate the Eddington ratio distribution of X-ray selected broad-line active galactic nuclei (AGN) in the redshift range 1.0<z<2.2, where the number density of AGNs peaks. Combining the optical and Subaru/FMOS near-infrared spectroscopy, we estimate black hole masses for broad-line AGNs in the Chandra Deep Field-South (CDF-S), Extended Chandra Deep Field-South (E-CDF-S), and the XMM-Newton Lockman Hole (XMM-LH) surveys. AGNs with similar black hole masses show a broad range of AGN bolometric luminosities, which are calculated from X-ray luminosities, indicating that the accretion rate of black holes is widely distributed. We find that a substantial fraction of massive black holes accreting significantly below the Eddington limit at z~2, in contrast to what is generally found for luminous AGNs at high redshift. Our analysis of observational selection biases indicates that the "AGN cosmic downsizing" phenomenon can be simply explained by the strong evolution of the co-moving number density at the bright end of the AGN luminosity function, together with the corresponding selection effects. However, it might need to consider a correlation between the AGN luminosity and the accretion rate of black holes that luminous AGNs have higher Eddington ratios than low-luminosity AGNs in order to understand the relatively small fraction of low-luminosity AGNs with high accretion rates in this epoch. Therefore, the observed downsizing trend could be interpreted as massive black holes with low accretion rates, which are relatively fainter than less massive black holes with efficient accretion.

Calibrating the photometric redshifts of >10^9 galaxies for upcoming weak lensing cosmology experiments is a major challenge for the astrophysics community. The path to obtaining the required spectroscopic redshifts for training and calibration is daunting, given the anticipated depths of the surveys and the difficulty in obtaining secure redshifts for some faint galaxy populations. Here we present an analysis of the problem based on the self-organizing map, a method of mapping the distribution of data in a high-dimensional space and projecting it onto a lower-dimensional representation. We apply this method to existing photometric data from the COSMOS survey selected to approximate the anticipated Euclid weak lensing sample, enabling us to robustly map the empirical distribution of galaxies in the multidimensional color space defined by the expected Euclid filters. Mapping this multicolor distribution lets us determine where - in galaxy color space - redshifts from current spectroscopic surveys exist and where they are systematically missing. Crucially, the method lets us determine whether a spectroscopic training sample is representative of the full photometric space occupied by the galaxies in a survey. We explore optimal sampling techniques and estimate the additional spectroscopy needed to map out the color-redshift relation, finding that sampling the galaxy distribution in color space in a systematic way can efficiently meet the calibration requirements. While the analysis presented here focuses on the Euclid survey, similar analysis can be applied to other surveys facing the same calibration challenge, such as DES, LSST, and WFIRST.

Using a grid of $\sim 2$ million elements ($\Delta z = 0.005$) adapted from COSMOS photometric redshift (photo-z) searches, we investigate the general properties of template-based photo-z likelihood surfaces. We find these surfaces are filled with numerous local minima and large degeneracies that generally confound rapid but "greedy" optimization schemes, even with additional stochastic sampling methods. In order to robustly and efficiently explore these surfaces, we develop BAD-Z [Brisk Annealing-Driven Redshifts (Z)], which combines ensemble Markov Chain Monte Carlo (MCMC) sampling with simulated annealing to sample arbitrarily large, pre-generated grids in approximately constant time. Using a mock catalog of 384,662 objects, we show BAD-Z samples $\sim 40$ times more efficiently compared to a brute-force counterpart while maintaining similar levels of accuracy. Our results represent first steps toward designing template-fitting photo-z approaches limited mainly by memory constraints rather than computation time.

We use observed optical to near infrared spectral energy distributions (SEDs) of 266 galaxies in the COSMOS survey to derive the wavelength dependence of the dust attenuation at high redshift. All of the galaxies have spectroscopic redshifts in the range z = 2 to 6.5. The presence of the CIV absorption feature, indicating that the rest-frame UV-optical SED is dominated by OB stars, is used to select objects for which the intrinsic, unattenuated spectrum has a well-established shape. Comparison of this intrinsic spectrum with the observed broadband photometric SED then permits derivation of the wavelength dependence of the dust attenuation. The derived dust attenuation curve is similar in overall shape to the Calzetti curve for local starburst galaxies. We also see the 2175 \AA~bump feature which is present in the Milky Way and LMC extinction curves but not seen in the Calzetti curve. The bump feature is commonly attributed to graphite or PAHs. No significant dependence is seen with redshift between sub-samples at z = 2 - 4 and z = 4 - 6.5. The 'extinction' curve obtained here provides a firm basis for color and extinction corrections of high redshift galaxy photometry.

We study how AGN activity changes across environments from galaxy pairs to clusters using $143\, 843$ galaxies with $z<0.2$ from the Sloan Digital Sky Survey (SDSS). Using a refined technique, we apply a continuous measure of AGN activity, characteristic of the ionization state of the narrow-line emitting gas. Changes in key emission-line ratios ([NII]$\lambda6548$/H$\alpha$, [OIII]$\lambda5007$/H$\beta$) between different samples allow us to disentangle different environmental effects while removing contamination. We confirm that galaxy interactions enhance AGN activity. However, conditions in the central regions of clusters are inhospitable for AGN activity even if galaxies are in pairs. These results can be explained through models of gas dynamics in which pair interactions stimulate the transfer of gas to the nucleus and clusters suppress gas availability for accretion onto the central black hole.

Recent observations indicate a remarkable similarity in the properties of evolving galaxies at fixed mass and redshift, prompting us to consider the possibility that most galaxies may evolve with a common history encompassing star formation, quasar accretion, and eventual quiescence. We quantify this by defining a "synchronization timescale" for galaxies as a function of mass and redshift that characterizes the extent to which different galaxies of a common mass are evolving in the same manner at various cosmic epochs. We measure this synchronization timescale using 9 different star-forming galaxy observations from the literature and SDSS quasar observations spanning $0 < z \lesssim 6$. Surprisingly, this synchronization timescale is a constant, approximately 1.5 Gyr for all combinations of mass and time. We also find that the ratio between the stellar mass of galaxies turning off star formation and black hole mass of turnoff quasars is approximately 30:1, much lower than the 500:1 for quiescent galaxies at low redshift. As a result, we propose a model in which the star-forming "main sequence", analogous quasar behavior, and other observations form a galactic evolution "main sequence", in which star formation occurs earliest, followed by supermassive black hole accretion, and feedback between the two are dominated by deterministic rather than stochastic processes.

Using the first 50% of data collected for the Spitzer Large Area Survey with Hyper-Suprime-Cam (SPLASH) observations on the 1.8 deg$^2$ Cosmological Evolution Survey (COSMOS) we estimate the masses and star formation rates of 3398 $M_*>10^{10}M_\odot $ star-forming galaxies at $4 < z < 6$ with a substantial population up to $M_* \gtrsim 10^{11.5} M_\odot$. We find that the strong correlation between stellar mass and star formation rate seen at lower redshift (the "main sequence" of star-forming galaxies) extends to $z\sim6$. The observed relation and scatter is consistent with a continued increase in star formation rate at fixed mass in line with extrapolations from lower-redshift observations. It is difficult to explain this continued correlation, especially for the most massive systems, unless the most massive galaxies are forming stars near their Eddington-limited rate from their first collapse. Furthermore, we find no evidence for moderate quenching at higher masses, indicating quenching either has not occurred prior to $z \sim 6$ or else occurs rapidly, so that few galaxies are visible in transition between star-forming and quenched.

Using a compilation of 25 studies from the literature, we investigate the evolution of the star-forming galaxy (SFG) Main Sequence (MS) in stellar mass and star formation rate (SFR) out to $z \sim 6$. After converting all observations to a common set of calibrations, we find a remarkable consensus among MS observations ($\sim 0.1$ dex 1$\sigma$ interpublication scatter). By fitting for time evolution of the MS in bins of constant mass, we deconvolve the observed scatter about the MS within each observed redshift bins. After accounting for observed scatter between different SFR indicators, we find the width of the MS distribution is $\sim 0.2$ dex and remains constant over cosmic time. Our best fits indicate the slope of the MS is likely time-dependent, with our best fit $\log\textrm{SFR}(M_*,t) = \left(0.84 \pm 0.02 - 0.026 \pm 0.003 \times t\right) \log M_* - \left(6.51 \pm 0.24 - 0.11 \pm 0.03 \times t\right)$, with $t$ the age of the Universe in Gyr. We use our fits to create empirical evolutionary tracks in order to constrain MS galaxy star formation histories (SFHs), finding that (1) the most accurate representations of MS SFHs are given by delayed-$\tau$ models, (2) the decline in fractional stellar mass growth for a "typical" MS galaxy today is approximately linear for most of its lifetime, and (3) scatter about the MS can be generated by galaxies evolving along identical evolutionary tracks assuming an initial $1\sigma$ spread in formation times of $\sim 1.4$ Gyr.

We took advantage of the observations carried out by XMM in the COSMOS field during 3.5 years, to study the long term variability of a large sample of AGN (638 sources), in a wide range of redshift (0.1<z<3.5) and X-ray luminosity ($10^{41}<$L(2-10)$<10^{45.5}$). Both a simple statistical method to asses the significance of variability, and the Normalized Excess Variance ($\sigma^{2}_{rms}$) parameter, where used to obtain a quantitative measurement of the variability. Variability is found to be prevalent in most AGN, whenever we have good statistic to measure it, and no significant differences between type-1 and type-2 AGN were found. A flat (slope -0.23+/-0.03) anti-correlation between $\sigma^{2}_{rms}$ and X-ray luminosity is found, when significantly variable sources are considered all together. When divided in three redshift bins, the anti-correlation becomes stronger and evolving with z, with higher redshift AGN being more variable. We prove however that this effect is due to the pre-selection of variable sources: considering all the sources with available $\sigma^{2}_{rms}$ measurement, the evolution in redshift disappears. For the first time we were also able to study the long term X-ray variability as a function of $M_{\rm BH}$ and Eddington ratio, for a large sample of AGN spanning a wide range of redshift. An anti-correlation between $\sigma^{2}_{rms}$ and $M_{\rm BH}$ is found, with the same slope of the anti-correlation between $\sigma^{2}_{rms}$ and X-ray luminosity, suggesting that the latter can be a byproduct of the former one. No clear correlation is found between $\sigma^{2}_{rms}$ and the Eddington ratio in our sample. Finally, no correlation is found between the X-ray $\sigma^{2}_{rms}$ and the optical variability.

We present an investigation into how well the properties of the accretion flow onto a supermassive black hole may be coupled to those of the overlying hot corona. To do so, we specifically measure the characteristic spectral index, Gamma, of a power-law energy distribution, over an energy range of 2 to 10 keV, for X-ray selected, broad-lined radio-quiet AGN up to z~2 in COSMOS and E-CDF-S. We test the previously reported dependence between Gamma and black hole mass, FWHM and Eddington ratio using a sample of AGN covering a broad range in these parameters based on both the Mg ii and H-alpha emission lines with the later afforded by recent near infrared spectroscopic observations using Subaru/FMOS. We calculate the Eddington ratios, lambda_Edd, for sources where a bolometric luminosity (L_Bol) has been presented in the literature, based on SED fitting, or, for sources where these data do not exist, we calculate L_Bol using a bolometric correction to the LX, derived from a relationship between the bolometric correction, and LX/L3000. From a sample of 69 X-ray bright sources (>250 counts), where Gamma can be measured with greatest precision, with an estimate of L_Bol, we find a statistically significant correlation between Gamma and lambda_Edd, which is highly significant with a chance probability of 6.59x10^-8. A statistically significant correlation between Gamma and the FWHM of the optical lines is confirmed, but at lower significance than with lambda_Edd indicating that lambda_Edd is the key parameter driving conditions in the corona. Linear regression analysis reveals that Gamma=(0.32+/-0.05)log10 lambda_Edd+(2.27+/-0.06) and Gamma=(-0.69+/-0.11)log10(FWHM/km/s)+(4.44+/-0.42). Our results on Gamma-lambda_Edd are in very good agreement with previous results. (ABRIDGED)

The H\beta\ emission line in a typical Type I quasar is composed of a broad base and a narrow core, with the core velocity typical of narrow-line region emission, and line-fitting routines typically assume this picture. We test the effects of removing this constraint, and find a substantial group of Type I quasars in the Sloan Digital Sky Survey catalog with H\beta\ emission line cores broader than 1200 km/s, above the velocity believed possible for gas in the quasar narrow-line region. We identify this group of "anomalous H\beta\ quasars" (AHQs) as a distinct population because of a variety of spectral and photometric signatures common to these AHQs but atypical of other quasars. These features are similar to some aspects of narrow-line Seyfert 1s and correlations identified by Eigenvector 1, but also contain distinct features that make AHQs difficult to classify. We demonstrate that AHQs comprise at least 11% and most likely approximately one quarter of the SDSS Type I quasar population at 0.2 < z < 0.8. For AHQs, the [O III]\lambda 4959,5007 profile is often better fit by de-linking it from the H\beta\ core, while a more standard linked fit produces a tight correlation between narrow- and broad-line velocities. We find that [O III] in AHQs sometimes has a standard narrow-line profile and other times matches the H\beta core, but is rarely in between the two, implying that the broadened core emission arises from a distinct physical region. Another feature of AHQs is a diminished [O II] line, which might indicate a connection between AHQs and the interstellar mediums of their host galaxies, through reduced photoionization or star formation. We find that it is difficult to produce AHQs using the current quasar standard model.

We present an analysis of broad emission lines observed in moderate-luminosity active galactic nuclei (AGNs), typical of those found in X-ray surveys of deep fields, with the aim to test the validity of single-epoch virial black hole mass estimates. We have acquired near-infrared (NIR) spectra of AGNs up to z ~ 1.8 in the COSMOS and Extended Chandra Deep Field-South Survey, with the Fiber Multi-Object Spectrograph (FMOS) mounted on the Subaru Telescope. These low-resolution NIR spectra provide a significant detection of the broad Halpha line that has been shown to be a reliable probe of black hole mass at low redshift. Our sample has existing optical spectroscopy which provides a detection of MgII, a broad emission line typically used for black hole mass estimation at z > 1. We carry out a spectral-line fitting procedure using both Halpha and MgII to determine the virial velocity of gas in the broad line region, the monochromatic continuum luminosity at 3000 A, and the total Halpha line luminosity. With a sample of 43 AGNs spanning a range of two decades in luminosity (i.e., L ~ 10^44-46 ergs/s), we find a tight correlation between the continuum and line luminosity with a distribution characterized by <log(L_3000/L_Halpha)> = 1.52 and a dispersion sigma = 0.16. There is also a close one-to-one relationship between the FWHM of Halpha and of MgII up to 10000 km/s with a dispersion of 0.14 in the distribution of the logarithm of their ratios. Both of these then lead to there being very good agreement between Halpha- and MgII-based masses over a wide range in black hole mass (i.e., M_BH ~ 10^7-9 M_sun). We do find a small offset in MgII-based masses, relative to those based on Halpha, of +0.17 dex and a dispersion sigma = 0.32. In general, these results demonstrate that local scaling relations, using MgII or Halpha, are applicable for AGN at moderate luminosities and up to z ~ 2.

We report on the discovery of a Type 1 quasar, SDSS 0956+5128, with a surprising combination of extreme velocity offsets. SDSS 0956+5128 is a broad-lined quasar exhibiting emission lines at three substantially different redshifts: a systemic redshift of z ~ 0.714 based on narrow emission lines, a broad MgII emission line centered 1200 km/s bluer than the systemic velocity, at z ~ 0.707, and broad H\alpha and H\beta emission lines centered at z ~ 0.690. The Balmer line peaks are 4100 km/s bluer than the systemic redshift. There are no previously known objects with such an extreme difference between broad MgII and broad Balmer emission. The two most promising explanations are either an extreme disk emitter or a high-velocity black hole recoil. However, neither explanation appears able to explain all of the observed features of SDSS 0956+5128, so the object may provide a challenge to our general understanding of quasar physics.

We consider the formation of dwarf carbon stars via accretion from a carbon AGB companion in light of the new 107 object sample of Downes et al. (2004). This sample is now large enough to allow good mass determination via comparison of a composite spectrum to theoretical atmospheric models. Carbon dwarfs of spectral type M are indeed main sequence M dwarfs with enhanced metallicity and carbon abundance. We also calculate the predicted abundance of both M and of F/G carbon dwarfs, and show that the latter should be falsifiable in the near future.

Anomalous Narrow-Line Quasars (ANLs) are a population of quasars with narrow H\beta, and sometimes [O III] broader than ~1000 km/s, in total comprising \sim 10-30% (most likely ~25%) of Type I quasars at 0.2 < z < 0.8. We find that virial masses using the H\beta and Mg II lines systematically differ for ANLs by an average of as much as 0.5 dex. Because the broad H\beta component width increases in ANLs but Mg II does not, we might suspect H\beta-based virial masses for ANLs are wrong but Mg II masses are correct. If this is due to an outflow reaching the lower-ionization potential H\beta line, C IV masses will be similarly flawed. However, we cannot be certain of this explanation without followup work, and may be unable to identify which quasars are ANLs at z > 0.8. Therefore, it is essential that ANLs be well-understood and well-modeled in order to allow the use of virial mass estimators on large optical spectroscopic catalogs, particularly at z < 0.4 and z > 2.0 where only one broad line is available for use in mass estimation.

Recent work using virial mass estimates and the quasar mass-luminosity plane has yielded several new puzzles regarding quasar accretion, including a sub-Eddington boundary on most quasar accretion, near-independence of the accretion rate from properties of the host galaxy, and a cosmic synchronization of accretion among black holes of a common mass. We consider how these puzzles might change if virial mass estimation turns out to have a systematic bias. As examples, we consider two recent claims of mass-dependent biases in MgII masses. Under any such correction, the surprising cosmic synchronization of quasar accretion rates and independence from the host galaxy remain. The slope and location of the sub-Eddington boundary are very sensitive to biases in virial mass estimation, and various mass calibrations appear to favor different possible physical explanations for feedback between the central black hole and its environment. The alternative mass estimators considered do not simply remove puzzling quasar behavior, but rather replace it with new puzzles that may be more difficult to solve than those using current virial mass estimators and the Shen et al. (2008) catalog.

Previous work on the quasar mass-luminosity plane indicates the possibility that quasars of the same central black hole mass might follow a common evolutionary track, independent of the properties of the host galaxy. We consider two simple models for the evolution of individual quasars. Requiring these tracks to lie within the observed quasar locus at all redshifts strongly constrains the model parameters, but does allow some solutions. These solutions include a family of tracks with similar shape but different initial masses that might match the observed quasar distributions at all redshifts z < 2.0. This family of solutions is characterized by short (1-2 Gyr) lifetimes, a duty cycle in which the quasar is on at least 25% of the time, and a rapid decline in Eddington ratio, perhaps with L/L_Edd ~ t^-6 or steeper.

We use 62185 quasars from the Sloan Digital Sky Survey Data Release 5 sample to explore the relationship between black hole mass and luminosity. Black hole masses were estimated based on the widths of their H{\beta}, MgII and CIV lines and adjacent continuum luminosities using standard virial mass estimate scaling laws. We find that, over the range 0.2 < z < 4.0, the most luminous low-mass quasars are at their Eddington luminosity, but the most luminous high-mass quasars in each redshift bin fall short of their Eddington luminosities, with the shortfall of the order of 10 or more at 0.2 < z < 0.6. We examine several potential sources of measurement uncertainty or bias and show that none of them can account for this effect. We also show the statistical uncertainty in virial mass estimation to have an upper bound of ~0.15 dex, smaller than the 0.4 dex previously reported. We also examine the highest mass quasars in every redshift bin in an effort to learn more about quasars that are about to cease their luminous accretion. We conclude that the quasar mass-luminosity locus contains a number of new puzzles that must be explained theoretically.

We use 62185 quasars from the Sloan Digital Sky Survey (SDSS) DR5 sample to explore the quasar mass-luminosity plane view of virial mass estimation. Previous work shows deviations of ~0.4 dex between virial and reverberation masses. The decline in quasar number density for the highest Eddington ratio quasars at each redshift provides an upper bound of between 0.13 and 0.29 dex for virial mass estimate statistical uncertainties. Across different redshift bins, the maximum possible MgII mass uncertainties average 0.15 dex, while H{\beta} uncertainties average 0.21 dex and CIV uncertainties average 0.27 dex. Any physical spread near the high-Eddington-ratio boundary will produce a more restrictive bound. A comparison of the sub-Eddington boundary slope using H{\beta} and MgII masses finds better agreement with uncorrected MgII masses than with recently proposed corrections. The best agreement for these bright objects is produced by a multiplicative correction by a factor of 1.19, smaller than the factor of 1.8 previously reported as producing the best agreement for the entire SDSS sample.

Comment on "Biases in the Quasar Mass-Luminosity Plane"

We use 62,185 quasars from the Sloan Digital Sky Survey DR5 sample and standard virial mass scaling laws based on the widths of H beta, Mg II, and C IV lines and adjacent continuum luminosities to explore the maximum mass of quasars as a function of redshift, which we find to be sharp and evolving. This evolution is in the sense that high-mass black holes cease their luminous accretion at higher redshift than lower-mass black holes. Further, turnoff for quasars at any given mass is more highly synchronized than would be expected given the dynamics of their host galaxies. We investigate potential signatures of the quasar turnoff mechanism, including a dearth of high-mass quasars at low Eddington ratio. These new results allow a closer examination of several common assumptions used in modeling quasar accretion and turnoff.

We discuss two new observational techniques that use observations of the Cosmic Microwave Background (CMB) to place constraints upon the mass, distance, and size distribution of small objects in the Kuiper Belt and inner Oort Cloud, collectively known as Trans-Neptunian Objects (TNOs). The first new technique considers the spectral distortion of the isotropic, or monopole, CMB by TNOs that have been heated by solar radiation to temperatures above that of the CMB. We apply this technique to the spectral measurements of the CMB by the Far Infrared Absolute Spectrophotometer (FIRAS) on the Cosmic Background Explorer (COBE). The second technique utilizes the change in amplitude of the TNO signal due to the orbital motion of the observer to separate the TNO signal from the invariant extra-galactic CMB and construct a map of the mass distribution in the outer Solar System. We estimate the ability of future CMB experiments to create such a map.

Recent theoretical ideas and observational claims suggest that the fine structure constant alpha may be variable. We examine a spectrum of models in which alpha is a function of a scalar field. Specifically, we consider three scenarios: oscillating alpha, monotonic time variation of alpha, and time-independent alpha that is spatially varying. We examine the constraints imposed upon these theories by cosmological observations, particle detector experiments, and ``fifth force'' experiments. These constraints are very strong on models involving oscillation, but cannot compete with bounds from the Oklo sub-nuclear reactor on models with monotonic time-like variation of alpha. One particular model with spatial variation is consistent with all current experimental and observational measurements, including those from two seemingly conflicting measurements of the fine structure constant using the many multiplet method on absorption lines.

We use the strong nebular lines of O III, 5007 A and 4959 A, to set a robust upper limit on the time dependence of the fine structure constant. We find |alpha^{-1} d alpha(t)/dt | < 2*10^{-13} yr, corresponding to Delta alpha/alpha(0) = (0.7 +- 1.4) * 10^{-4} for quasars with 0.16 < z < 0.80 obtained from the SDSS Early Data Release. Using a blind analysis, we show that the upper limit given here is invariant with respect to 17 different ways of selecting the sample and analyzing the data. As a by-product, we show that the ratio of transition probabilities for the 5007 A and the 4959 A lines is2.99 +- 0.02, in good agreement with (but more accurate than) theoretical estimates. We compare and contrast the O III emission line method used here with the Many-Multiplet method that has been used recently to suggest evidence for a time-dependent alpha. In an Appendix, we have analyzed the larger quasar sample from the SDSS Data Release One Sample and find Delta alpha/alpha(0) = (1.2 +- 0.7) * 10^{-4}.