Using the catalogues of galaxy clusters from The Three Hundred project, modelled with both hydrodynamic simulations, (Gadget-X and Gadget-MUSIC), and semi-analytic models (SAMs), we study the scatter and self-similarity of the profiles and distributions of the baryonic components of the clusters: the stellar and gas mass, metallicity, the stellar age, gas temperature, and the (specific) star formation rate. Through comparisons with observational results, we find that the shape and the scatter of the gas density profiles matches well the observed trends including the reduced scatter at large radii which is a signature of self-similarity suggested in previous studies. One of our simulated sets, Gadget-X, reproduces well the shape of the observed temperature profile, while Gadget-MUSIC has a higher and flatter profile in the cluster centre and a lower and steeper profile at large radii. The gas metallicity profiles from both simulation sets, despite following the observed trend, have a relatively lower normalisation. The cumulative stellar density profiles from SAMs are in better agreement with the observed result than both hydrodynamic simulations which show relatively higher profiles. The scatter in these physical profiles, especially in the cluster centre region, shows a dependence on the cluster dynamical state and on the cool-core/non-cool-core dichotomy. The stellar age, metallicity and (s)SFR show very large scatter, which are then presented in 2D maps. We also do not find any clear radial dependence of these properties. However, the brightest central galaxies have distinguishable features compared to the properties of the satellite galaxies.

We present the first predictions for the $L_{\rm [CII]}$ - SFR relation and [CII] luminosity function (LF) in the Epoch of Reionization (EoR) based on cosmological hydrodynamics simulations using the SIMBA suite plus radiative transfer calculations via S\'IGAME. The sample consists of 11,137 galaxies covering halo mass $\log M_{\rm halo}\in$[9, 12.4] $M_\odot$, star formation rate SFR$\in$[0.01, 330] $M_\odot$ yr$^{-1}$, and metallicity $<Z_{\rm gas}>_{\rm SFR}\in$[0.1, 0.9] $Z_\odot$. The simulated $L_{\rm [CII]}$-SFR relation is consistent with the range observed, but with a spread of $\simeq$0.3 dex at the high end of SFR ($>$100 $M_\odot$ yr$^{-1}$) and $\simeq$0.6 dex at the lower end, and there is tension between our predictions and the values of $L_{\rm [CII]}$ above 10$^{8.5}$ $L_\odot$ observed in some galaxies reported in the literature. The scatter in the $L_{\rm [CII]}$-SFR relation is mostly driven by galaxy properties, such that at a given SFR, galaxies with higher molecular gas mass and metallicity have higher $L_{\rm [CII]}$. The [CII] LF predicted by SIMBA is consistent with the upper limits placed by the only existing untargeted flux-limited [CII] survey at the EoR (ASPECS) and those predicted by semi-analytic models. We compare our results with existing models and discuss differences responsible for the discrepant slopes in the $L_{\rm [CII]}$-SFR relatiion.

High-redshift Lyman-alpha blobs (LABs) are an enigmatic class of objects that have been the subject of numerous observational and theoretical investigations. It is of particular interest to determine the dominant power sources for the copious luminosity, as direct emission from HII regions, cooling gas, and fluorescence due to the presence of active galactic nuclei (AGN) can all contribute significantly. In this paper, we present the first theoretical model to consider all of these physical processes in an attempt to develop an evolutionary model for the origin of high-z LABs. This is achieved by combining a series of high-resolution cosmological zoom-in simulations with ionization and Lyman-alpha (Lya) radiative transfer models. We find that massive galaxies display a range of Lya luminosities and spatial extents (which strongly depend on the limiting surface brightness used) over the course of their lives, though regularly exhibit luminosities and sizes consistent with observed LABs. The model LABs are typically powered from a combination of recombination in star-forming galaxies, as well as cooling emission from gas associated with accretion. When AGN are included in the model, the fluorescence caused by AGN-driven ionization can be a significant contributor to the total Lya luminosity as well. We propose that the presence of an AGN may be predicted from the Gini coefficient of the blob's surface brightness. Within our modeled mass range, there are no obvious threshold physical properties that predict appearance of LABs, and only weak correlations of the luminosity with the physical properties of the host galaxy. This is because the emergent Lya luminosity from a system is a complex function of the gas temperature, ionization state, and Lya escape fraction.

The past two decades have seen a tremendous investment in observational facilities that promise to reveal new and unprecedented discoveries about the universe. In comparison, the investment in theoretical work is completely dwarfed, even though theory plays a crucial role in the interpretation of these observations, predicting new types of phenomena, and informing observing strategies. In this white paper, we argue that in order to reach the promised critical breakthroughs in astrophysics over the next decade and well beyond, the national agencies must take a serious approach to investment in theoretical astrophysics research. We discuss the role of theory in shaping our understanding of the universe, and then we provide a multi-level strategy, from the grassroots to the national, to address the current underinvestment in theory relative to observational work.

Galaxy cluster outskirts mark the transition region from the mildly non-linear cosmic web to the highly non-linear, virialised, cluster interior. It is in this transition region that the intra-cluster medium (ICM) begins to influence the properties of accreting galaxies and groups, as ram pressure impacts a galaxy's cold gas content and subsequent star formation rate. Conversely, the thermodynamical properties of the ICM in this transition region should also feel the influence of accreting substructure (i.e. galaxies and groups), whose passage can drive shocks. In this paper, we use a suite of cosmological hydrodynamical zoom simulations of a single galaxy cluster, drawn from the nIFTy comparison project, to study how the dynamics of substructure accreted from the cosmic web influences the thermodynamical properties of the ICM in the cluster's outskirts. We demonstrate how features evident in radial profiles of the ICM (e.g. gas density and temperature) can be linked to strong shocks, transient and short-lived in nature, driven by the passage of substructure. The range of astrophysical codes and galaxy formation models in our comparison are broadly consistent in their predictions (e.g. agreeing when and where shocks occur, but differing in how strong shocks will be); this is as we would expect of a process driven by large-scale gravitational dynamics and strong, inefficiently radiating, shocks. This suggests that mapping such shock structures in the ICM in a cluster's outskirts (via e.g. radio synchrotron emission) could provide a complementary measure of its recent merger and accretion history.

We present a framework for characterizing the large scale movement of baryons relative to dark matter in cosmological simulations, requiring only the initial conditions and final state of the simulation. This is performed using the spread metric which quantifies the distance in the final conditions between initially neighbouring particles, and by analysing the baryonic content of final haloes relative to that of the initial Lagrangian regions defined by their dark matter component. Applying this framework to the SIMBA cosmological simulations, we show that 40% (10%) of cosmological baryons have moved $> 1h^{-1}~ {\rm Mpc}^{-1}$ ($3h^{-1}~ {\rm Mpc}^{-1}$) by $z=0$, due primarily to entrainment of gas by jets powered by AGN, with baryons moving up to $12h^{-1}~ {\rm Mpc}^{-1}$ away in extreme cases. Baryons decouple from the dynamics of the dark matter component due to hydrodynamic forces, radiative cooling, and feedback processes. As a result, only 60% of the gas content in a given halo at $z=0$ originates from its Lagrangian region, roughly independent of halo mass. A typical halo in the mass range $M_{\rm vir} = 10^{12}$--$10^{13}{\rm M}_\odot$ only retains 20% of the gas originally contained in its Lagrangian region. We show that up to 20% of the gas content in a typical Milky Way mass halo may originate in the region defined by the dark matter of another halo. This inter-Lagrangian baryon transfer may have important implications for the origin of gas and metals in the circumgalactic medium of galaxies, as well as for semi-analytic models of galaxy formation and "zoom-in" simulations.

We describe the survey for galaxies in the fields surrounding 9 sightlines to far-UV bright, z~1 quasars that define the COS Absorption Survey of Baryon Harbors (CASBaH) program. The photometry and spectroscopy that comprise the dataset come from a mixture of public surveys (SDSS, DECaLS) and our dedicated efforts on private facilities (Keck, MMT, LBT). We report the redshifts and stellar masses for 5902 galaxies within ~10 comoving-Mpc (cMpc) of the sightlines with a median of z=0.28 and M_* ~ 10^(10.1) Msun. This dataset, publicly available as the CASBaH specDB, forms the basis of several recent and ongoing CASBaH analyses. Here, we perform a clustering analysis of the galaxy sample with itself (auto-correlation) and against the set of OVI absorption systems (cross-correlation) discovered in the CASBaH quasar spectra with column densities N(O^+5) >= 10^(13.5)/cm^2. For each, we describe the measured clustering signal with a power-law correlation function xi(r) = (r/r_0)^(-gamma) and find that (r_0,gamma) = (5.48 +/- 0.07 h_100^-1 Mpc, 1.33 +/- 0.04) for the auto-correlation and (6.00 +/- 1 h^-1 Mpc, 1.25 +/- 0.18) for galaxy-OVI cross-correlation. We further estimate a bias factor of b_gg = 1.3 +/- 0.1 from the galaxy-galaxy auto-correlation indicating the galaxies are hosted by halos with mass M_halo ~ 10^(12.1 +/- 0.05) Msun. Finally, we estimate an OVI-galaxy bias factor b_OVI = 1.0 +/- 0.1 from the cross-correlation which is consistent with OVI absorbers being hosted by dark matter halos with typical mass M_halo ~ 10^(11) Msun. Future works with upcoming datasets (e.g., CGM^2) will improve upon these results and will assess whether any of the detected OVI arises in the intergalactic medium.

We study the spin alignment of galaxies and halos with respect to filaments and walls of the cosmic web, identified with DisPerSE, using the SIMBA simulation from z=0-2. Massive halos' spins are oriented perpendicularly to their closest filament's axis and walls, while low mass halos tend to have their spins parallel to filaments and in the plane of walls. A similar mass-dependent spin flip is found for galaxies, albeit with a weaker signal particularly at low mass and low-z, suggesting that galaxies' spins retain memory of their larger-scale environment. Low-z star-forming and rotation-dominated galaxies tend to have spins parallel to nearby filaments, while quiescent and dispersion-dominated galaxies show preferentially perpendicular orientation; the star formation trend can be fully explained by the stellar mass correlation, but the morphology trend cannot. There is a strong dependence on HI mass, such that high-HI galaxies tend to have parallel spins while low-HI galaxies are perpendicular, which persists even when matching samples in stellar mass, suggesting that HI content traces anisotropic infall more faithfully than the stellar component. Finally, at fixed stellar mass, the strength of spin alignments correlates with the filament's density, with parallel alignment for galaxies in high density environments. These findings are consistent with conditional tidal torque theory, and highlight a significant correlation between galactic spin and the larger scale tides that are important e.g. for interpreting weak lensing studies. SIMBA allows us to rule out numerical grid locking as the cause of previously-seen low mass alignment.

Following Cui et al. 2018 (hereafter Paper I) on the classification of large-scale environments (LSE) at z = 0, we push our analysis to higher redshifts and study the evolution of LSE and the baryon distributions in them. Our aim is to investigate how baryons affect the LSE as a function of redshift. In agreement with Paper I, the baryon models have negligible effect on the LSE over all investigated redshifts. We further validate the conclusion obtained in Paper I that the gas web is an unbiased tracer of total matter -- even better at high redshifts. By separating the gas mainly by temperature, we find that about 40 per cent of gas is in the so-called warm-hot intergalactic medium (WHIM). This fraction of gas mass in the WHIM decreases with redshift, especially from z = 1 (29 per cent) to z = 2.1 (10 per cent). By separating the whole WHIM gas mass into the four large-scale environments (i.e. voids, sheets, filaments, and knots), we find that about half of the WHIM gas is located in filaments. Although the total gas mass in WHIM decreases with redshift, the WHIM mass fractions in the different LSE seem unchanged.

The Sloan Digital Sky Survey (SDSS) was revolutionary because of the extraordinary breadth and ambition of its optical imaging and spectroscopy. We argue that a 'sub-millimeter SDSS' - a sensitive large-area imaging+spectroscopic survey in the sub-mm window - will revolutionize our understanding of galaxy evolution in the early Universe. By detecting the thermal dust continuum emission and atomic and molecular line emission of galaxies out to z~10 it will be possible to measure the redshifts, star formation rates, dust and gas content of hundreds of thousands of high-z galaxies down to ~L*. Many of these galaxies will have counterparts visible in the deep optical imaging of the Large Synoptic Survey Telescope. This 3D map of galaxy evolution will span the peak epoch of galaxy formation all the way back to cosmic dawn, measuring the co-evolution of the star formation rate density and molecular gas content of galaxies, tracking the production of metals and charting the growth of large-scale structure.

Recent numerical studies of the dark matter density profiles of massive galaxy clusters ($M_{\rm halo} > 10^{15}$M$_{\odot}$) show that their median radial mass density profile remains unchanged up to $z > 1$, displaying a highly self-similar evolution. We verify this by using the data set of the THE THREE HUNDRED project, i.e. 324 cluster-sized haloes as found in full physics hydrodynamical simulations. We track the progenitors of the mass-complete sample of clusters at $z=0$, and find that their median shape is already in place by $z=2.5$. However, selecting a dynamically relaxed subsample ($\sim16$ per cent of the clusters), we observe a shift of the scale radius $r_s$ towards larger values at earlier times. Classifying the whole sample by formation time, this evolution is understood as a result of a two-phase halo mass accretion process. Early-forming clusters -- identified as relaxed today -- have already entered their slow accretion phase, hence their mass growth occurs mostly at the outskirts. Late-forming clusters -- which are still unrelaxed today -- are in their fast accretion phase, thus the central region of the clusters is still growing. We conclude that the density profile of galaxy clusters shows a profound self-similarity out to redshifts $z\sim2.5$. This result holds for both gas and total density profiles when including baryonic physics, as reported here for two rather distinct sub-grid models.

The $\rm H_2$ mass of molecular clouds has traditionally been traced by the CO(J=1-0) rotational transition line. This said, CO is relatively easily photodissociated, and can also be destroyed by cosmic rays, thus rendering some fraction of molecular gas to be "CO-dark". We investigate the amount and physical properties of CO-dark gas in two $z \sim 0$ disc galaxies, and develop predictions for the expected intensities of promising alternative tracers ([CI 609 $\mu$m and [CII] 158 $\mu$m emission). We do this by combining cosmological zoom simulations of disc galaxies with thermal-radiative-chemical equilibrium interstellar medium (ISM) calculations to model the predicted H~\textsc{i} and $\rm H_2$ abundances and CO(J=1-0), [CI] 609 $\mu$m and [CII] 158 $\mu$m emission properties. Our model treats the ISM as a collection of radially stratified clouds whose properties are dictated by their volume and column densities, the gas-phase metallicity, and the interstellar radiation field and cosmic ray ionization rates. Our main results follow. Adopting an observationally motivated definition of CO-dark gas, i.e. $\rm H_2$ gas with $W_{\rm CO} < 0.1 $ K-km/s, we find that a significant amount ($> 50\%$) of the total $\rm H_2$ mass lies in CO-dark gas, most of which is diffuse gas, poorly shielded due to low dust column density. The CO-dark molecular gas tends to be dominated by [CII], though [CI] also serves as a bright tracer of the dark gas in many instances. At the same time, [CII] also tends to trace neutral atomic gas. As a result, when we quantify the conversion factors for the three carbon-based tracers of molecular gas, we find that [CI] suffers the least contamination from diffuse atomic gas, and is relatively insensitive to secondary parameters.

The next generation Very Large Array (ngVLA) will revolutionize our understanding of the distant Universe via the detection of cold molecular gas in the first galaxies. Its impact on studies of galaxy characterization via detailed gas dynamics will provide crucial insight on dominant physical drivers for star-formation in high redshift galaxies, including the exchange of gas from scales of the circumgalactic medium down to resolved clouds on mass scales of $\sim10^{5}\,M_\odot$. In this study, we employ a series of high-resolution, cosmological, hydrodynamic zoom simulations from the MUFASA simulation suite and a CASA simulator to generate mock ngVLA observations of a $z\sim4.5$ gas rich star-forming galaxy. Using the DESPOTIC radiative transfer code that encompasses simultaneous thermal, chemical, and statistical equilibrium in calculating the molecular and atomic level transitions of CO from ALMA for comparison. We find that observations of CO(1-0) are especially important for tracing the systemic redshift of the galaxy and the total mass of the well-shielded molecular gas reservoir, while even CO(2-1) can predominantly trace denser gas regions distinct from CO(1-0). The factor of 100 times improvement in mapping speed for the ngVLA beyond the Jansky VLA and the proposed ALMA Band 1 will make these detailed, high-resolution imaging and kinematic studies of CO(1-0) routine at $z\sim2-5$.

We introduce the THE THREE HUNDRED project, an endeavour to model 324 large galaxy clusters with full-physics hydrodynamical re-simulations. Here we present the data set and study the differences to observations for fundamental galaxy cluster properties and scaling relations. We find that the modelled galaxy clusters are generally in reasonable agreement with observations with respect to baryonic fractions and gas scaling relations at redshift z = 0. However, there are still some (model-dependent) differences, such as central galaxies being too massive, and galaxy colours (g - r) being bluer (about 0.2 dex lower at the peak position) than in observations. The agreement in gas scaling relations down to 10^{13} h^{-1} M_{\odot} between the simulations indicates that particulars of the sub-grid modelling of the baryonic physics only has a weak influence on these relations. We also include - where appropriate - a comparison to three semi-analytical galaxy formation models as applied to the same underlying dark-matter-only simulation. All simulations and derived data products are publicly available.

We analyze the metal accumulation in dwarf and spiral galaxies by following the history of metal enrichment and outflows in a suite of twenty high-resolution simulated galaxies. These simulations agree with the observed stellar and gas-phase mass-metallicity relation, an agreement that relies on large fractions of the produced metals escaping into the CGM. For instance, in galaxies with Mvir ~ 1e9.5 -- 1e10 solar masses, we find that about ~ 85% of the available metals are outside of the galactic disk at z = 0, although the fraction decreases to a little less than half in Milky Way-mass galaxies. In many cases, these metals are spread far beyond the virial radius. We analyze the metal deficit within the ISM and stars in the context of previous work tracking the inflow and outflow of baryons. Outflows are prevalent across the entire mass range, as is reaccretion. We find that between 40 and 80% of all metals removed from the galactic disk are later reaccreted. The outflows themselves are metal enriched relative to the ISM by a factor of 0.2 dex because of the correspondence between sites of metal enrichment and outflows. As a result, the metal mass loading factor scales as eta_metals \propto v_circ^-0.91, a somewhat shallower scaling than the total mass loading factor. We analyze the simulated galaxies within the context of analytic chemical evolution models by determining their net metal expulsion efficiencies, which encapsulate the rates of metal loss and reaccretion. We discuss these results in light of the inflow and outflow properties necessary for reproducing the mass-metallicity relation.

In this paper, we provide a physical model for the origin of variations in the shapes and bump strengths of dust attenuation laws in galaxies by combining a large suite of cosmological "zoom-in" galaxy formation simulations with 3D Monte Carlo dust radiative transfer calculations. We model galaxies over 3 orders of magnitude in stellar mass, ranging from Milky Way like systems through massive galaxies at high-redshift. Critically, for these calculations we employ a constant underlying dust extinction law in all cases, and examine how the role of geometry and radiative transfer effects impact the resultant attenuation curves. Our main results follow. Despite our usage of a constant dust extinction curve, we find dramatic variations in the derived attenuation laws. The slopes of normalized attenuation laws depend primarily on the complexities of star-dust geometry. Increasing fractions of unobscured young stars flatten normalized curves, while increasing fractions of unobscured old stars steepen curves. Similar to the slopes of our model attenuation laws, we find dramatic variation in the 2175 Angstrom ultraviolet (UV) bump strength, including a subset of curves with little to no bump. These bump strengths are primarily influenced by the fraction of unobscured O and B stars in our model, with the impact of scattered light having only a secondary effect. Taken together, these results lead to a natural relationship between the attenuation curve slope and 2175 Angstrom bump strength. Finally, we apply these results to a 25 Mpc/h box cosmological hydrodynamic simulation in order to model the expected dispersion in attenuation laws at integer redshifts from z=0-6. A significant dispersion is expected at low redshifts, and decreases toward z=6. We provide tabulated results for the best fit median attenuation curve at all redshifts.

We test the galactic outflow model by probing associated galaxies of four strong intergalactic CIV absorbers at $z=5$--6 using the Hubble Space Telescope (HST) ACS ramp narrowband filters. The four strong CIV absorbers reside at $z=5.74$, $5.52$, $4.95$, and $4.87$, with column densities ranging from $N_{\rm{CIV}}=10^{13.8}$ cm$^{-2}$ to $10^{14.8}$ cm$^{-2}$. At $z=5.74$, we detect an i-dropout Ly$\alpha$ emitter (LAE) candidate with a projected impact parameter of 42 physical kpc from the CIV absorber. This LAE candidate has a Ly$\alpha$-based star formation rate (SFR$_{\rm{Ly\alpha}}$) of 2 $M_\odot$ yr$^{-1}$ and a UV-based SFR of 4 $M_\odot$ yr$^{-1}$. Although we cannot completely rule out that this $i$-dropout emitter may be an [OII] interloper, its measured properties are consistent with the CIV powering galaxy at $z=5.74$. For CIV absorbers at $z=4.95$ and $z=4.87$, although we detect two LAE candidates with impact parameters of 160 kpc and 200 kpc, such distances are larger than that predicted from the simulations. Therefore we treat them as non-detections. For the system at $z=5.52$, we do not detect LAE candidates, placing a 3-$\sigma$ upper limit of SFR$_{\rm{Ly\alpha}}\approx 1.5\ M_\odot$ yr$^{-1}$. In summary, in these four cases, we only detect one plausible CIV source at $z=5.74$. Combining the modest SFR of the one detection and the three non-detections, our HST observations strongly support that smaller galaxies (SFR$_{\rm{Ly\alpha}} \lesssim 2\ M_\odot$ yr$^{-1}$) are main sources of intergalactic CIV absorbers, and such small galaxies play a major role in the metal enrichment of the intergalactic medium at $z\gtrsim5$.

The MIGHTEE large survey project will survey four of the most well-studied extragalactic deep fields, totalling 20 square degrees to $\mu$Jy sensitivity at Giga-Hertz frequencies, as well as an ultra-deep image of a single ~1 square degree MeerKAT pointing. The observations will provide radio continuum, spectral line and polarisation information. As such, MIGHTEE, along with the excellent multi-wavelength data already available in these deep fields, will allow a range of science to be achieved. Specifically, MIGHTEE is designed to significantly enhance our understanding of, (i) the evolution of AGN and star-formation activity over cosmic time, as a function of stellar mass and environment, free of dust obscuration; (ii) the evolution of neutral hydrogen in the Universe and how this neutral gas eventually turns into stars after moving through the molecular phase, and how efficiently this can fuel AGN activity; (iii) the properties of cosmic magnetic fields and how they evolve in clusters, filaments and galaxies. MIGHTEE will reach similar depth to the planned SKA all-sky survey, and thus will provide a pilot to the cosmology experiments that will be carried out by the SKA over a much larger survey volume.

We utilise a series of high-resolution cosmological zoom simulations of galaxy formation to investigate the relationship between the ultraviolet (UV) slope, beta, and the ratio of the infrared luminosity to UV luminosity (IRX) in the spectral energy distributions (SEDs) of galaxies. We employ dust radiative transfer calculations in which the SEDs of the stars in galaxies propagate through the dusty interstellar medium. Our main goals are to understand the origin of, and scatter in the IRX-beta relation; to assess the efficacy of simplified stellar population synthesis screen models in capturing the essential physics in the IRX-beta relation; and to understand systematic deviations from the canonical local IRX-beta relations in particular populations of high-redshift galaxies. Our main results follow. Galaxies that have young stellar populations with relatively cospatial UV and IR emitting regions and a Milky Way-like extinction curve fall on or near the standard Meurer relation. This behaviour is well captured by simplified screen models. Scatter in the IRX-beta relation is dominated by three major effects: (i) older stellar populations drive galaxies below the relations defined for local starbursts due to a reddening of their intrinsic UV SEDs; (ii) complex geometries in high-z heavily star forming galaxies drive galaxies toward blue UV slopes owing to optically thin UV sightlines; (iii) shallow extinction curves drive galaxies downward in the IRX-beta plane due to lowered NUV/FUV extinction ratios. We use these features of the UV slopes of galaxies to derive a fitting relation that reasonably collapses the scatter back toward the canonical local relation. Finally, we use these results to develop an understanding for the location of two particularly enigmatic populations of galaxies in the IRX-beta plane: z~2-4 dusty star forming galaxies, and z>5 star forming galaxies.

The recent low value of Planck (2016) integrated optical depth to Thomson scattering suggests that the reionization occurred fairly suddenly, disfavoring extended reionization scenarios. This will have a significant impact on the 21cm power spectrum. Using a semi-numerical framework, we improve our model from Hassan et al. (2016) to include time-integrated ionisation and recombination effects, and find that this leads to more sudden reionisation. It also yields larger HII bubbles which leads to an order of magnitude more 21cm power on large scales, while suppressing the small scale ionization power. Local fluctuations in the neutral hydrogen density play the dominant role in boosting the 21cm power spectrum on large scales, while recombinations are subdominant. We use a Monte Carlo Markov Chain approach to constrain our model to observations of the star formation rate functions at z = 6,7,8 from Bouwens et al. (2015), the Planck (2016) optical depth measurements, and the Becker & Bolton (2013) ionising emissivity data at z~5. We then use this constrained model to perform 21cm forecasting for LOFAR, HERA, and SKA in order to determine how well such data can characterise the sources driving reionisation. We find that the 21cm power spectrum alone can somewhat constrain the halo mass dependence of ionising sources, the photon escape fraction and ionising amplitude, but combining the 21cm data with other current observations enables us to separately constrain all these parameters. Our framework illustrates how 21cm data can play a key role in understanding the sources and topology of reionisation as observations improve.

We present Herschel PACS observations of the [CII] 158 micron emission line in a sample of 24 intermediate mass (9<logM$_\ast$/M$_\odot$<10) and low metallicity (0.4< Z/Z$_\odot$<1.0) galaxies from the xCOLD GASS survey. Combining them with IRAM CO(1-0) measurements, we establish scaling relations between integrated and molecular region [CII]/CO(1-0) luminosity ratios as a function of integrated galaxy properties. A Bayesian analysis reveals that only two parameters, metallicity and offset from the star formation main sequence, $\Delta$MS, are needed to quantify variations in the luminosity ratio; metallicity describes the total dust content available to shield CO from UV radiation, while $\Delta$MS describes the strength of this radiation field. We connect the [CII]/CO luminosity ratio to the CO-to-H$_2$ conversion factor and find a multivariate conversion function $\alpha_{CO}$, which can be used up to z~2.5. This function depends primarily on metallicity, with a second order dependence on $\Delta$MS. We apply this to the full xCOLD GASS and PHIBSS1 surveys and investigate molecular gas scaling relations. We find a flattening of the relation between gas mass fraction and stellar mass at logM$_\ast$/M$_\odot$<10. While the molecular gas depletion time varies with sSFR, it is mostly independent of mass, indicating that the low L$_{CO}$/SFR ratios long observed in low mass galaxies are entirely due to photodissociation of CO, and not to an enhanced star formation efficiency.

In this paper we investigate environment driven gas depletion in satellite galaxies, taking full advantage of the atomic hydrogen (HI) spectral stacking technique to quantify the gas content for the entire gas-poor to -rich regime. We do so using a multi-wavelength sample of 10,600 satellite galaxies, selected according to stellar mass (log M$_{\star}$/M$_{\odot}$ $\geq$ 9) and redshift (0.02 $\leq$ z $\leq$ 0.05) from the Sloan Digital Sky Survey, with HI data from the Arecibo Legacy Fast ALFA (ALFALFA) survey. Using key HI-to-stellar mass scaling relations, we present evidence that the gas content of satellite galaxies is, to a significant extent, dependent on the environment in which a galaxy resides. For the first time, we demonstrate that systematic environmental suppression of gas content at both fixed stellar mass and fixed specific star formation rate (sSFR) in satellite galaxies begins in halo masses typical of the group regime (log M$_{h}$/M$_{\odot}$ < 13.5), well before galaxies reach the cluster environment. We also show that environment driven gas depletion is more closely associated to halo mass than local density. Our results are then compared with state-of-the-art semi-analytic models and hydrodynamical simulations and discussed within this framework, showing that more work is needed if models are to reproduce the observations. We conclude that the observed decrease of gas content in the group and cluster environments cannot be reproduced by starvation of the gas supply alone and invoke fast acting processes such as ram-pressure stripping of cold gas to explain this.

We use the combined data from the COS-GASS and COS-Halos surveys to characterize the Circum-Galactic Medium (CGM) surrounding typical low-redshift galaxies in the mass range $\rm~M_*\sim~10^{9.5-11.5}~M_{\odot} $, and over a range of impact parameters extending to just beyond the halo virial radius ($\rm~R_{vir}$). We find the radial scale length of the distributions of the equivalent widths of the Lyman~$\alpha$ and Si III absorbers to be 0.9 and 0.4 $\rm~R_{vir}$, respectively. The radial distribution of equivalent widths is relatively uniform for the blue galaxies, but highly patchy (low covering fraction) for the red galaxies. We also find that the Lyman~$\alpha$ and Si III equivalent widths show significant positive correlations with the specific star-formation rate (sSFR) of the galaxy. We find a surprising lack of correlations between the halo mass (virial velocity) and either the velocity dispersions or velocity offsets of the Lyman~$\alpha$ lines. The ratio of the velocity offset to the velocity dispersion for the Lyman~$\alpha$ absorbers has a mean value of $\sim$ 4, suggesting that a given the line-of-sight is intersecting a dynamically coherent structure in the CGM rather than a sea of orbiting clouds. The kinematic properties of the CGM are similar in the blue and red galaxies, although we find that a significantly larger fraction of the blue galaxies have large Lyman~$\alpha$ velocity offsets (>200 km s$^{-1}$). We show that - if the CGM clouds represent future fuel for star-formation - our new results could imply a large drop in the specific star-formation rate across the galaxy mass-range we probe.

Although there has been much progress in understanding how galaxies evolve, we still do not understand how and when they stop forming stars and become quiescent. We address this by applying our galaxy spectral energy distribution models, which incorporate physically motivated star formation histories (SFHs) from cosmological simulations, to a sample of quiescent galaxies at $0.2<z<2.1$. A total of 845 quiescent galaxies with multi-band photometry spanning rest-frame ultraviolet through near-infrared wavelengths are selected from the CANDELS dataset. We compute median SFHs of these galaxies in bins of stellar mass and redshift. At all redshifts and stellar masses, the median SFHs rise, reach a peak, and then decline to reach quiescence. At high redshift, we find that the rise and decline are fast, as expected because the Universe is young. At low redshift, the duration of these phases depends strongly on stellar mass. Low-mass galaxies ($\log(M_{\ast}/M_{\odot})\sim9.5$) grow on average slowly, take a long time to reach their peak of star formation ($\gtrsim 4$ Gyr), and the declining phase is fast ($\lesssim 2$ Gyr). Conversely, high-mass galaxies ($\log(M_{\ast}/M_{\odot})\sim11$) grow on average fast ($\lesssim 2$ Gyr), and, after reaching their peak, decrease the star formation slowly ($\gtrsim 3$ Gyr). These findings are consistent with galaxy stellar mass being a driving factor in determining how evolved galaxies are, with high-mass galaxies being the most evolved at any time (i.e., downsizing). The different durations we observe in the declining phases also suggest that low- and high-mass galaxies experience different quenching mechanisms that operate on different timescales.

The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a new 400-800MHz radio interferometer under development for deployment in South Africa. HIRAX will comprise 1024 six meter parabolic dishes on a compact grid and will map most of the southern sky over the course of four years. HIRAX has two primary science goals: to constrain Dark Energy and measure structure at high redshift, and to study radio transients and pulsars. HIRAX will observe unresolved sources of neutral hydrogen via their redshifted 21-cm emission line (`hydrogen intensity mapping'). The resulting maps of large-scale structure at redshifts 0.8-2.5 will be used to measure Baryon Acoustic Oscillations (BAO). HIRAX will improve upon current BAO measurements from galaxy surveys by observing a larger cosmological volume (larger in both survey area and redshift range) and by measuring BAO at higher redshift when the expansion of the universe transitioned to Dark Energy domination. HIRAX will complement CHIME, a hydrogen intensity mapping experiment in the Northern Hemisphere, by completing the sky coverage in the same redshift range. HIRAX's location in the Southern Hemisphere also allows a variety of cross-correlation measurements with large-scale structure surveys at many wavelengths. Daily maps of a few thousand square degrees of the Southern Hemisphere, encompassing much of the Milky Way galaxy, will also open new opportunities for discovering and monitoring radio transients. The HIRAX correlator will have the ability to rapidly and eXperimentciently detect transient events. This new data will shed light on the poorly understood nature of fast radio bursts (FRBs), enable pulsar monitoring to enhance long-wavelength gravitational wave searches, and provide a rich data set for new radio transient phenomena searches. This paper discusses the HIRAX instrument, science goals, and current status.

We present a study of stellar populations in a sample of spectroscopically-confirmed Lyman-break galaxies (LBGs) and Ly$\alpha$ emitters (LAEs) at $5.7<z<7$. These galaxies have deep optical and infrared images from Subaru, $HST$, and $Spitzer$/IRAC. We focus on a subset of 27 galaxies with IRAC detections, and characterize their stellar populations utilizing galaxy synthesis models based on the multi-band data and secure redshifts. By incorporating nebular emission estimated from the observed Ly$\alpha$ flux, we are able to break the strong degeneracy of model spectra between young galaxies with prominent nebular emission and older galaxies with strong Balmer breaks. The results show that our galaxies cover a wide range of ages from several to a few hundred million years (Myr), and a wide range of stellar masses from $\sim10^8$ to $\sim10^{11}$ $M_{\odot}$. These galaxies can be roughly divided into an `old' subsample and a `young' subsample. The `old' subsample consists of galaxies older than 100 Myr, with stellar masses higher than $10^9$ $M_{\odot}$. The galaxies in the `young' subsample are younger than $\sim$30 Myr, with masses ranging between $\sim10^8$ and $\sim3\times10^9$ $M_{\odot}$. Both subsamples display a correlation between stellar mass and star-formation rate (SFR), but with very different normalizations. The average specific SFR (sSFR) of the `old' subsample is 3--4 Gyr$^{-1}$, consistent with previous studies of `normal' star-forming galaxies at $z\ge6$. The average sSFR of the `young' subsample is an order of magnitude higher, likely due to starburst activity. Our results also indicate little or no dust extinction in the majority of the galaxies, as already suggested by their steep rest-frame UV slopes. Finally, LAEs and LBGs with strong Ly$\alpha$ emission are indistinguishable in terms of age, stellar mass, and SFR.

We present a study exploring the nature and properties of the Circum-Galactic Medium (CGM) and its connection to the atomic gas content in the interstellar medium (ISM) of galaxies as traced by the HI 21cm line. Our sample includes 45 low-z (0.026-0.049) galaxies from the GALEX Arecibo SDSS Survey. Their CGM was probed via absorption in the spectra of background Quasi-Stellar Objects at impact parameters of 63 to 231kpc. The spectra were obtained with the Cosmic Origins Spectrograph aboard the Hubble Space Telescope. We detected neutral hydrogen (Ly$\alpha$ absorption-lines) in the CGM of 92% of the galaxies. We find the radial profile of the CGM as traced by the Ly$\alpha$ equivalent width can be fit as an exponential with a scale length of roughly the virial radius of the dark matter halo. We found no correlation between the orientation of sightline relative to the galaxy major axis and the Ly$\alpha$ equivalent width. The velocity spread of the circumgalactic gas is consistent with that seen in the atomic gas in the interstellar medium. We find a strong correlation (99.8% confidence) between the gas fraction (M(HI)/M*) and the impact-parameter-corrected Ly$\alpha$ equivalent width. This is stronger than the analogous correlation between corrected Ly$\alpha$ equivalent width and SFR/M* (97.5% confidence). These results imply a physical connection between the HI disk and the CGM, which is on scales an order-of-magnitude larger. This is consistent with the picture in which the HI disk is nourished by accretion of gas from the CGM.

In this paper, we present the FATS (Feature Analysis for Time Series) library. FATS is a Python library which facilitates and standardizes feature extraction for time series data. In particular, we focus on one application: feature extraction for astronomical light curve data, although the library is generalizable for other uses. We detail the methods and features implemented for light curve analysis, and present examples for its usage.

In this paper we present the MOSFIRE Deep Evolution Field (MOSDEF) survey. The MOSDEF survey aims to obtain moderate-resolution (R=3000-3650) rest-frame optical spectra (~3700-7000 Angstrom) for ~1500 galaxies at 1.37<z<3.80 in three well-studied CANDELS fields: AEGIS, COSMOS, and GOODS-N. Targets are selected in three redshift intervals: 1.37<z<1.70, 2.09<z<2.61, and 2.95<z<3.80, down to fixed H_AB (F160W) magnitudes of 24.0, 24.5 and 25.0, respectively, using the photometric and spectroscopic catalogs from the 3D-HST survey. We target both strong nebular emission lines (e.g., [OII], Hbeta, [OIII], 5008, Halpha, [NII], and [SII]) and stellar continuum and absorption features (e.g., Balmer lines, Ca-II H and K, Mgb, 4000 Angstrom break). Here we present an overview of our survey, the observational strategy, the data reduction and analysis, and the sample characteristics based on spectra obtained during the first 24 nights. To date, we have completed 21 masks, obtaining spectra for 591 galaxies. For ~80% of the targets we derive a robust redshift from either emission or absorption lines. In addition, we confirm 55 additional galaxies, which were serendipitously detected. The MOSDEF galaxy sample includes unobscured star-forming, dusty star-forming, and quiescent galaxies and spans a wide range in stellar mass (~10^9-10^11.5 Msol) and star formation rate (~10^0-10^3 Msol/yr). The spectroscopically confirmed sample is roughly representative of an H-band limited galaxy sample at these redshifts. With its large sample size, broad diversity in galaxy properties, and wealth of available ancillary data, MOSDEF will transform our understanding of the stellar, gaseous, metal, dust, and black hole content of galaxies during the time when the universe was most active.

Four years after the last LISA meeting, the NASA Astrophysics Data System (ADS) finds itself in the middle of major changes to the infrastructure and contents of its database. In this paper we highlight a number of features of great importance to librarians and discuss the additional functionality that we are currently developing. Starting in 2011, the ADS started to systematically collect, parse and index full-text documents for all the major publications in Physics and Astronomy as well as many smaller Astronomy journals and arXiv e-prints, for a total of over 3.5 million papers. Our citation coverage has doubled since 2010 and now consists of over 70 million citations. We are normalizing the affiliation information in our records and, in collaboration with the CfA library and NASA, we have started collecting and linking funding sources with papers in our system. At the same time, we are undergoing major technology changes in the ADS platform which affect all aspects of the system and its operations. We have rolled out and are now enhancing a new high-performance search engine capable of performing full-text as well as metadata searches using an intuitive query language which supports fielded, unfielded and functional searches. We are currently able to index acknowledgments, affiliations, citations, funding sources, and to the extent that these metadata are available to us they are now searchable under our new platform. The ADS private library system is being enhanced to support reading groups, collaborative editing of lists of papers, tagging, and a variety of privacy settings when managing one's paper collection. While this effort is still ongoing, some of its benefits are already available through the ADS Labs user interface and API at http://adslabs.org/adsabs/

We analyze the low-redshift (z~0.2) circumgalactic medium by comparing absorption-line data from the COS-Halos Survey to absorption around a matched galaxy sample from two cosmological hydrodynamic simulations. The models include different prescriptions for galactic outflows, namely hybrid energy/momentum driven wind (ezw), and constant winds (cw). We extract for comparison direct observables including equivalent widths, covering factors, ion ratios, and kinematics. Both wind models are generally in good agreement with these observations for HI and certain low ionization metal lines, but show poorer agreement with higher ionization metal lines including SiIII and OVI that are well-observed by COS-Halos. These discrepancies suggest that both wind models predict too much cool, metal-enriched gas and not enough hot gas, and/or that the metals are not sufficiently well-mixed. This may reflect our model assumption of ejecting outflows as cool and unmixing gas. Our ezw simulation includes a heuristic prescription to quench massive galaxies by super-heating its ISM gas, which we show yields sufficient low ionisation absorption to be broadly consistent with observations, but also substantial OVI absorption that is inconsistent with data, suggesting that gas around quenched galaxies in the real Universe does not cool. At impact parameters of <50 kpc, recycling winds dominate the absorption of low ions and even HI, while OVI almost always arises from metals ejected longer than 1 Gyr ago. The similarity between the wind models is surprising, since we show that they differ substantially in their predicted amount and phase distribution of halo gas. We show that this similarity owes mainly to our comparison here at fixed stellar mass rather than at fixed halo mass in our previous works, which suggests that CGM properties are more closely tied to the stellar mass of galaxies rather than halo mass.

Numerous theoretical models have long proposed that a strong He II 1640 emission line is the most prominent and unique feature of massive Population III (Pop III) stars in high redshift galaxies. The He II 1640 line strength can constrain the mass and IMF of Pop III stars. We use F132N narrowband filter on the Hubble Space Telescope's (HST) Wide Field Camera 3 (WFC3) to look for strong He II lambda 1640 emission in the galaxy BDF-521 at z=7.01, one of the most distant spectroscopically-confirmed galaxies to date. Using deep F132N narrowband imaging, together with our broadband imaging with F125W and F160W filters, we do not detect He II emission from this galaxy, but place a 2-sigma upper limit on the flux of 5.3x10^-19 ergs s^-1 cm^-2. This measurement corresponds to a 2-sigma upper limit on the Pop III star formation rate (SFR_PopIII) of ~ 0.2 M_solar yr^-1, assuming a Salpeter IMF with 50< M/M_solar < 1000. From the high signal-to-noise broadband measurements in F125W and F160W, we fit the UV continuum for BDF-521. The spectral flux density is ~ 3.6x 10^-11 lambda^-2.32 ergs s^-1 cm^-2 A^-1, which corresponds to an overall unobscured SFR of ~ 5 M_solar yr^-1. Our upper limit on SFR_PopIII suggests that massive Pop III stars represent < 4% of the total star formation. Further, the HST high resolution imaging suggests that BDF-521 is an extremely compact galaxy, with a half-light radius of 0.6 kpc.

The Herschel very wide-field surveys have charted hundreds of square degrees in multiple far-IR (FIR) bands. While the Sloan Digital Sky Survey (SDSS) is currently the best resource for optical counterpart identifications over such wide areas, it does not detect a large number of Herschel FIR sources and leaves their nature undetermined. As a test case, we studied seven "SDSS-invisible", very bright 250um sources (S_{250} > 55 mJy) in the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) fields where we have a rich multi-wavelength data set. We took a new approach to decompose the FIR sources, using the near-IR or the optical images directly for position priors. This is an improvement over the previous decomposition efforts where the priors are from mid-IR data that still suffer from the source blending problem in the first place. We found that in most cases the single Herschel sources are made of multiple components that are not necessarily at the same redshifts. Our decomposition succeeded in identifying and extracting their major contributors. We show that these are all ULIRGs at z ~ 1--2 whose high L_IR is mainly due to dust-obscured star formation. Most of them would not be selected as sub-mm galaxies. They all have complicated morphologies indicative of merger or violent instability, and their stellar populations are heterogeneous in terms of stellar masses, ages and formation histories. Their current ULIRG phases are of various degrees of importance in their stellar mass assembly. Our practice provides a promising starting point to develop an automatic routine to reliably study bright Herschel sources.

We examine the star formation histories (SFHs) of galaxies in smoothed particle hydrodynamics (SPH) simulations, compare them to parametric models that are commonly used in fitting observed galaxy spectral energy distributions, and examine the efficacy of these parametric models as practical tools for recovering the physical parameters of galaxies. The commonly used tau-model, with SFR ~ exp(-t/tau), provides a poor match to the SFH of our SPH galaxies, with a mismatch between early and late star formation that leads to systematic errors in predicting colours and stellar mass-to-light ratios. A one-parameter lin-exp model, with SFR ~ t*exp(-t/tau), is much more successful on average, but it fails to match the late-time behavior of the bluest, most actively star-forming galaxies and the passive, "red and dead" galaxies. We introduce a 4-parameter model, which transitions from lin-exp to a linear ramp after a transition time, which describes our simulated galaxies very well. We test the ability of these parametrised models to recover (at z=0, 0.5, and 1) the stellar mass-to-light ratios, specific star formation rates, and stellar population ages from the galaxy colours, computed from the full SPH star formation histories using the FSPS code of Conroy et al. (2009). Fits with tau-models systematically overestimate M/L by ~ 0.2 dex, overestimate population ages by ~ 1-2 Gyr, and underestimate sSFR by ~ 0.05 dex. Fits with lin-exp are less biased on average, but the 4-parameter model yields the best results for the full range of galaxies. Marginalizing over the free parameters of the 4-parameter model leads to slightly larger statistical errors than 1-parameter fits but essentially removes all systematic biases, so this is our recommended procedure for fitting real galaxies.

We derive physical properties of 10 submillimeter galaxies located in the CANDELS coverage of the GOODS-S field. The galaxies were first identified as submillimeter sources with the LABOCA bolometer and subsequently targeted for 870um continuum observation with ALMA. The high angular resolution of the ALMA imaging allows secure counterparts to be identified in the CANDELS multiband dataset. The CANDELS data provide deep photometric data from UV through near-infrared wavelengths. Using synthetic spectral energy distributions, we derive photometric redshifts, stellar masses, extinction, ages, and the star formation history. The redshift range is z=1.65-4.76, with two of the galaxies located at z>4. Two SMG counterparts have stellar masses 2-3 orders of magnitude lower than the rest. The remaining SMG counterparts have stellar masses around 1x10^11 Msun. The stellar population in the SMGs is typically older than the expected duration of the submillimeter phase, suggesting that the star formation history of submillimeter galaxies is more complex than a single burst. Non-parametric morphology indices suggest that the SMG counterparts are among the most asymmetric systems compared with galaxies of the same stellar mass and redshift. The HST images shows that 3 of the SMGs are associated with on-going mergers. The remaining counterparts are isolated. Estimating the dust and molecular gas mass from the submm fluxes, and comparing with our stellar masses shows that the molecular gas mass fraction of SMGs is ~28% and that the final stellar mass is likely to be (1-2)x10^11 Msun.

We study the high-ionization phase and kinematics of the circumgalactic medium around low-redshift galaxies using a sample of 23 Lyman Limit Systems (LLSs) at 0.08<z<0.93 observed with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope. In Lehner et al. (2013), we recently showed that low-z LLSs have a bimodal metallicity distribution. Here we extend that analysis to search for differences between the high-ion and kinematic properties of the metal-poor and metal-rich branches. We find that metal-rich LLSs tend to show higher O VI columns and broader O VI profiles than metal-poor LLSs. The total H I line width (dv90 statistic) in LLSs is not correlated with metallicity, indicating that the H I kinematics alone cannot be used to distinguish inflow from outflow and gas recycling. Among the 17 LLSs with O VI detections, all but two show evidence of kinematic sub-structure, in the form of O VI-H I centroid offsets, multiple components, or both. Using various scenarios for how the metallicity in the high-ion and low-ion phases of each LLS compare, we constrain the ionized hydrogen column in the O VI phase to lie in the range log N(H II)~17.6-20. The O VI phase of LLSs is a substantial baryon reservoir, with M(high-ion)~10^{8.5-10.9}(r/150 kpc)^2 solar masses, similar to the mass in the low-ion phase. Accounting for the O VI phase approximately doubles the contribution of low-z LLSs to the cosmic baryon budget.

Galactic winds are a key physical mechanism for understanding galaxy formation and evolution, yet empirical and theoretical constraints for the character of winds are limited and discrepant. Recent empirical models find that local star-forming galaxies have a deficit of oxygen that scales with galaxy stellar mass. The oxygen deficit provides unique empirical constraints on the magnitude of mass loss, composition of outflowing material and metal reaccretion onto galaxies. We formulate the oxygen deficit constraints so they may be easily implemented into theoretical models of galaxy evolution. We parameterize an effective metal loading factor which combines the uncertainties of metal outflows and metal reaccretion into a single function of galaxy virial velocity. We determine the effective metal loading factor by forward-fitting the oxygen deficit. The effective metal loading factor we derive has important implications for the implementation of mass loss in models of galaxy evolution.

We present the design and methods of the COS-Halos survey, a systematic investigation of the gaseous halos of 44 z = 0.15-0.35 galaxies using background QSOs observed with the Cosmic Origins Spectrograph aboard the Hubble Space Telescope. This survey has yielded 39 spectra of z_em ~ 0.5 QSOs with S/N ~ 10-15 per resolution element. The QSO sightlines pass within 150 physical kpc of the galaxies, which span early and late types over stellar mass log M* / Msun= 9.5 - 11.5. We find that the CGM exhibits strong HI, averaging 1 Ang in Lya equivalent width out to 150 kpc, with 100% covering fraction for star-forming galaxies and 75% covering for passive galaxies. We find good agreement in column densities between this survey and previous studies over similar range of impact parameter. There is weak evidence for a difference between early- and late-type galaxies in the strength and distribution of HI. Kinematics indicate that the detected material is bound to the host galaxy, such that >~90% of the detected column density is confined within +/-200 km s^-1 of the galaxies. This material generally exists well below the halo virial temperatures at T<~ 10^5 K. We evaluate a number of possible origin scenarios for the detected material, and in the end favor a simple model in which the bulk of the detected HI arises in a bound, cool, low-density photoionized diffuse medium that is generic to all L* galaxies and may harbor a total gaseous mass comparable to galactic stellar masses.

We present a new optical spectroscopic survey of 1777 'star-forming' ('SF') and 366 'non-star-forming' ('non-SF') galaxies at redshifts z < 1 (2143 in total), 22 AGN and 423 stars, observed by instruments such as DEIMOS, VIMOS and GMOS, in 3 fields containing 5 quasi-stellar objects (QSOs) with HST UV spectroscopy. We also present a new spectroscopic survey of 165 'strong' (10^14 < NHI < 10^17 cm^-2), and 489 'weak' (10^13 < NHI < 10^14 cm^-2) intervening HI absorption line systems at z < 1 (654 in total), observed in the spectra of 8 QSOs by COS and FOS on the HST. Combining these new data with previously published galaxy catalogs such as VVDS and GDDS, we have gathered a sample of 654 HI absorption systems and 17509 galaxies at transverse scales < 50 Mpc. We present observational results on the HI-galaxy and galaxy-galaxy correlations at transverse scales r < 10 Mpc, and the HI-HI auto-correlation at transverse scales r < 2 Mpc. The two-point correlation functions are measured both along and transverse to the line-of-sight. We constrain the HI-galaxy statistical connection, as a function of both HI column density and galaxy star-forming activity. Our results are consistent with the following conclusions: (1) the bulk of HI systems on Mpc scales have little velocity dispersion (<120 km/s) with respect to the bulk of galaxies; (2) the vast majority of strong HI systems and SF galaxies are distributed in the same locations, together with 75+-15% of non-SF galaxies, all of which typically reside in dark matter haloes of similar masses; (3) 25+-15% of non-SF galaxies reside in galaxy clusters and are not correlated with strong HI systems at scales < 2 Mpc; and (4) 50% of weak HI systems reside within galaxy voids (hence not correlated with galaxies), and are confined in dark matter haloes of masses smaller than those hosting... [abridged]

We investigate the differential effects of metal cooling and galactic stellar winds on the cosmological formation of individual galaxies with three sets of cosmological, hydrodynamical zoom simulations of 45 halos in the mass range 10^11<M_halo<10^13M_sun. Models including both galactic winds and metal cooling (i) suppress early star formation at z>1 and predict reasonable star formation histories, (ii) produce galaxies with high cold gas fractions (30-60 per cent) at high redshift, (iii) significantly reduce the galaxy formation efficiencies for halos (M_halo<10^12M_sun) at all redshifts in agreement with observational and abundance matching constraints, (iv) result in high-redshift galaxies with reduced circular velocities matching the observed Tully-Fisher relation at z~2, and (v) significantly increase the sizes of low-mass galaxies (M_stellar<3x10^10M_sun) at high redshift resulting in a weak size evolution - a trend in agreement with observations. However, the low redshift (z<0.5) star formation rates of massive galaxies are higher than observed (up to ten times). No tested model predicts the observed size evolution for low-mass and high-mass galaxies simultaneously. Due to the delayed onset of star formation in the wind models, the metal enrichment of gas and stars is delayed and agrees well with observational constraints. Metal cooling and stellar winds are both found to increase the ratio of in situ formed to accreted stars - the relative importance of dissipative vs. dissipationless assembly. For halo masses below ~10^12M_sun, this is mainly caused by less stellar accretion and compares well to predictions from semi-analytical models but still differs from abundance matching models. For higher masses, the fraction of in situ stars is over-predicted due to the unrealistically high star formation rates at low redshifts.

We develop a dust efflux model of radiation pressure acting on dust grains which successfully reproduces the relation between stellar mass, dust opacity and star formation rate observed in local star-forming galaxies. The dust content of local star-forming galaxies is set by the competition between the physical processes of dust production and dust loss in our model. The dust loss rate is proportional to the dust opacity and star formation rate. Observations of the relation between stellar mass and star formation rate at several epochs imply that the majority of local star-forming galaxies are best characterized as having continuous star formation histories. Dust loss is a consequence of sustained interaction of dust with the radiation field generated by continuous star formation. Dust efflux driven by radiation pressure rather than dust destruction offers a more consistent physical interpretation of the dust loss mechanism. By comparing our model results with the observed relation between stellar mass, dust extinction and star formation rate in local star-forming galaxies we are able to constrain the timescale and magnitude of dust loss. The timescale of dust loss is long and therefore dust is effluxed in a "Slow Flow". Dust loss is modest in low mass galaxies but massive galaxies may lose up to 70~80% of their dust over their lifetime. Our Slow Flow model shows that mass loss driven by dust opacity and star formation may be an important physical process for understanding normal star-forming galaxy evolution.

We present a deep LBT/LBC U-band imaging survey (9 deg2) covering the NOAO Bootes field. A total of 14,485 Lyman Break Galaxies (LBGs) at z~3 are selected, which are used to measure the rest-frame UV luminosity function (LF). The large sample size and survey area reduce the LF uncertainties due to Poisson statistics and cosmic variance by >3 compared to previous studies. At the bright end, the LF shows excess power compared to the best-fit Schechter function, which can be attributed to the contribution of $z\sim3$ quasars. We compute the rest-frame near-infrared LF and stellar mass function (SMF) of z~3 LBGs based on the R-band and IRAC [4.5 micro m]-band flux relation. We investigate the evolution of the UV LFs and SMFs between z~7 and z~3, which supports a rising star formation history in the LBGs. We study the spatial correlation function of two bright LBG samples and estimate their average host halo mass. We find a tight relation between the host halo mass and the galaxy star formation rate (SFR),which follows the trend predicted by the baryonic accretion rate onto the halo, suggesting that the star formation in LBGs is fueled by baryonic accretion through the cosmic web. By comparing the SFRs with the total baryonic accretion rates, we find that cosmic star formation efficiency is about 5%-20% and it does not evolve significantly with redshift, halo mass, or galaxy luminosity.

We present a detailed structural and morphological study of a large sample of spectroscopically-confirmed galaxies at z >= 6, using deep HST near-IR broad-band images and Subaru optical narrow-band images. The galaxy sample consists of 51 Lyman-alpha emitters (LAEs) at z ~ 5.7, 6.5, and 7.0, and 16 Lyman-break galaxies (LBGs) at 5.9 < z < 6.5. These galaxies exhibit a wide range of rest-frame UV continuum morphology in the HST images, from compact features to multiple component systems. The fraction of merging/interacting galaxies reaches 40% ~ 50% at the brightest end of M_1500 <= -20.5 mag. The intrinsic half-light radii r_{hl,in}, after correction for PSF broadening, are roughly between r_{hl,in} ~ 0.05" (0.3 kpc) and 0.3" (1.7 kpc) at M_1500 <= -19.5 mag. The median r_{hl,in} value is 0.16" (~0.9 kpc). This is consistent with the sizes of bright LAEs and LBGs at z >= 6 in previous studies. In addition, more luminous galaxies tend to have larger sizes, exhibiting a weak size-luminosity relation r_{hl,in} \propto L^{0.14} at M_1500 <= -19.5 mag. The slope of 0.14 is significantly flatter than those in fainter LBG samples. We discuss the morphology of z >= 6 galaxies with nonparametric methods, including the CAS system and the Gini and M_20 parameters, and demonstrate their validity through simulations. We search for extended Lyman-alpha emission halos around LAEs at z ~ 5.7 and 6.5, by stacking a number of narrow-band images. We do not find evidence of extended halos predicted by cosmological simulations. Such Lyman-alpha halos, if they exist, could be weaker than predicted. Finally, we investigate any positional misalignment between UV continuum and Lyman-alpha emission in LAEs. While the two positions are generally consistent, several merging galaxies show significant positional differences. This is likely caused by a disturbed ISM distribution due to merging activity.

We use the chemical evolution predictions of cosmological hydrodynamic simulations with our latest theoretical stellar population synthesis, photoionization and shock models to predict the strong line evolution of ensembles of galaxies from z=3 to the present day. In this paper, we focus on the brightest optical emission-line ratios, [NII]/H-alpha and [OIII]/H-beta. We use the optical diagnostic Baldwin-Phillips-Terlevich (BPT) diagram as a tool for investigating the spectral properties of ensembles of active galaxies. We use four redshift windows chosen to exploit new near-infrared multi-object spectrographs. We predict how the BPT diagram will appear in these four redshift windows given different sets of assumptions. We show that the position of star-forming galaxies on the BPT diagram traces the ISM conditions and radiation field in galaxies at a given redshift. Galaxies containing AGN form a mixing sequence with purely star-forming galaxies. This mixing sequence may change dramatically with cosmic time, due to the metallicity sensitivity of the optical emission-lines. Furthermore, the position of the mixing sequence may probe metallicity gradients in galaxies as a function of redshift, depending on the size of the AGN narrow line region. We apply our latest slow shock models for gas shocked by galactic-scale winds. We show that at high redshift, galactic wind shocks are clearly separated from AGN in line ratio space. Instead, shocks from galactic winds mimic high metallicity starburst galaxies. We discuss our models in the context of future large near-infrared spectroscopic surveys.

We present deep HST near-IR and Spitzer mid-IR observations of a large sample of spectroscopically-confirmed galaxies at z >= 6. The sample consists of 51 Lyman-alpha emitters (LAEs) at z ~ 5.7, 6.5, and 7.0, and 16 Lyman-break galaxies (LBGs) at 5.9 < z < 6.5. The near-IR images were mostly obtained with WFC3 in the F125W and F160W bands, and the mid-IR images were obtained with IRAC in the 3.6um and 4.5um bands. Our galaxies also have deep optical imaging data from Subaru Suprime-Cam. We utilize the multi-band data and secure redshifts to derive their rest-frame UV properties. These galaxies have steep UV continuum slopes roughly between beta ~ -1.5 and -3.5, with an average value of beta ~ -2.3, slightly steeper than the slopes of LBGs in previous studies. The slope shows little dependence on UV continuum luminosity except for a few of the brightest galaxies. We find a statistically significant excess of galaxies with slopes around beta ~ -3, suggesting the existence of very young stellar populations with extremely low metallicity and dust content. Our galaxies have moderately strong rest-frame Lyman-alpha equivalent width (EW) in a range of ~10 to ~200 \AA. The star-formation rates are also moderate, from a few to a few tens solar masses per year. The LAEs and LBGs in this sample share many common properties, implying that LAEs represent a subset of LBGs with strong Lyman-alpha emission. Finally, the comparison of the UV luminosity functions between LAEs and LBGs suggests that there exists a substantial population of faint galaxies with weak Lyman-alpha emission (EW < 20 \AA) that could be the dominant contribution to the total ionizing flux at z >= 6.

In this paper we present a detailed study of the structures and morphologies of a sample of 1188 massive galaxies with Mstar>10^10Msun between redshifts z=1-3 within the Ultra Deep Survey (UDS) region of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) field. Using this sample we determine how galaxy structure and morphology evolve with time. We visually classify our sample into disks, ellipticals and peculiar systems and correct for redshift effects on our classifications through simulations. We find evolution in the fractions of galaxies at a given visual classification as a function of redshift. The peculiar population is dominant at z>2 with a substantial spheroid population, and a negligible disk population. We compute the transition redshift, ztrans, where the combined fraction of spheroids and disks is equal to that of peculiar galaxies, as ztrans=1.86+/-0.62 for galaxies in our stellar mass range. We find that this changes as a function of stellar mass, with Hubble-type systems becoming dominant at higher redshifts for higher mass galaxies (ztrans=2.22+/-0.82), than for the lower mass galaxies (ztrans=1.73+/-0.57). Higher mass galaxies become morphologically settled before their lower mass counterparts, a form of morphological downsizing. We furthermore compare our visual classifications with Sersic index, the concentration, asymmetry and clumpiness (CAS) parameters, star formation rate and rest frame U-B colour. We find links between the colour of a galaxy, its star formation rate and how extended or peculiar it appears. Finally, we discuss the negligible z>2 disk fraction based on visual morphologies and speculate that this is an effect of forming disks appearing peculiar through processes such as violent disk instabilities or mergers. We conclude that to properly define high redshift morphology a new and more exact classification scheme is needed.

One key goal of the Hubble Space Telescope Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey is to track galaxy evolution back to z ~ 8. Its two-tiered "wide and deep" strategy bridges significant gaps in existing near-infrared surveys. Here we report on z ~ 8 galaxy candidates selected as F105W-band dropouts in one of its deep fields, which covers 50.1 square arcmin to 4 ks depth in each of three near-infrared bands in the Great Observatories Origins Deep Survey southern field. Two of our candidates have J<26.2 mag, and are > 1 mag brighter than any previously known F105W-dropouts. We derive constraints on the bright-end of the rest-frame ultraviolet luminosity function of galaxies at z ~ 8, and show that the number density of such very bright objects is higher than expected from the previous Schechter luminosity function estimates at this redshift. Another two candidates are securely detected in Spitzer Infrared Array Camera images, which are the first such individual detections at z ~ 8. Their derived stellar masses are on the order of a few x 10^9 M_sun, from which we obtain the first measurement of the high-mass end of the galaxy stellar mass function at z ~ 8. The high number density of very luminous and very massive galaxies at z ~ 8, if real, could imply a large stellar-to-halo mass ratio and an efficient conversion of baryons to stars at such an early time.

We report new observations of circumgalactic gas in the halos of early type galaxies obtained by the COS-Halos Survey with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope. We find that detections of HI surrounding early type galaxies are typically as common and strong as around star-forming galaxies, implying that the total mass of circumgalactic material is comparable in the two populations. For early type galaxies, the covering fraction for HI absorption above 10^16 cm^2 is ~40-50% within ~150 kpc. Line widths and kinematics of the detected material show it to be cold (T ~< 10^5 K) in comparison to the virial temperature of the host halos. The implied masses of cool, photoionized CGM baryons may be up to 10^9 --- 10^11 Msun. Contrary to some theoretical expectations, strong halo HI absorbers do not disappear as part of the quenching of star-formation. Even passive galaxies retain significant reservoirs of halo baryons which could replenish the interstellar gas reservoir and eventually form stars. This halo gas may feed the diffuse and molecular gas that is frequently observed inside ETGs.

CO measurements of z~1-4 galaxies have found that their baryonic gas fractions are significantly higher than galaxies at z=0, with values ranging from 20-80 %. Here, we suggest that the gas fractions inferred from observations of star-forming galaxies at high-z are overestimated, owing to the adoption of locally-calibrated CO-H2 conversion factors (Xco). Evidence from both observations and numerical models suggest that Xco varies smoothly with the physical properties of galaxies, and that Xco can be parameterised simply as a function of both gas phase metallicity and observed CO surface brightness. When applying this functional form, we find fgas ~10-40 % in galaxies with M*=10^10-10^12 Msun at high-z. Moreover, the scatter in the observed fgas-M* relation is lowered by a factor of two. The lower inferred gas fractions arise physically because the interstellar media of high-z galaxies have higher velocity dispersions and gas temperatures than their local counterparts, which results in an Xco that is lower than the z=0 value for both quiescent discs and starbursts. We further compare these gas fractions to those predicted by cosmological galaxy formation models. We show that while the canonically inferred gas fractions from observations are a factor of 2-3 larger at a given stellar mass than predicted by models, our rederived Xco values for z=1-4 galaxies results in revised gas fractions that agree significantly better with the simulations.

In this contribution we present the first census of oxygen in star-forming galaxies in the local universe. We examine three samples of galaxies with metallicities and star formation rates at z = 0.07, 0.8 and 2.26, including the SDSS and DEEP2 surveys. We infer the total mass of oxygen produced and mass of oxygen found in the gas-phase from our local SDSS sample. The star formation history is determined by requiring that galaxies evolve along the relation between stellar mass and star formation rate observed in our three samples. We show that the observed relation between stellar mass and star formation rate for our three samples is consistent with other samples in the literature. The mass-metallicity relation is well established for our three samples and from this we empirically determine the chemical evolution of star-forming galaxies. Thus, we are able to simultaneously constrain the star formation rates and metallicities of galaxies over cosmic time allowing us to estimate the mass of oxygen locked up in stars. Combining this work with independent measurements reported in the literature we conclude that the loss of oxygen from the interstellar medium of local star-forming galaxies is likely to be a ubiquitous process with the oxygen mass loss scaling (almost) linearly with stellar mass. We estimate the total baryonic mass loss and argue that only a small fraction of the baryons inferred from cosmological observations accrete onto galaxies.

We present two-dimensional, integral field spectroscopy covering the rest-frame wavelengths of strong optical emission lines in nine sub-mm-luminous galaxies (SMGs) at 2.0<z<2.7. The GEMINI-NIFS and VLT-SINFONI imaging spectroscopy allows the mapping of the gas morphologies and dynamics within the sources, and we measure an average Halpha velocity dispersion of sigma=220+-80km/s and an average half light radius of r=3.7+-0.8kpc. The average dynamical measure, V_obs/2sigma=0.9+-0.1 for the SMGs, is higher than in more quiescent star-forming galaxies at the same redshift, highlighting a difference in the dynamics of the two populations. The SMGs' far-infrared SFRs, measured using Herschel-SPIRE far-infrared photometry, are on average 370+-90Mo/yr which is ~2 times higher than the extinction corrected SFRs of the more quiescent star-forming galaxies. Six of the SMGs in our sample show strong evidence for kinematically distinct multiple components with average velocity offsets of 200+-100km/s and average projected spatial offsets of 8+-2kpc, which we attribute to systems in the early stages of major mergers. Indeed all SMGs are classified as mergers from a kinemetry analysis of the velocity and dispersion field asymmetry. We bring together our sample with seven other SMGs with IFU observations to describe the ionized gas morphologies and kinematics in a sample of 16 SMGs. By comparing the velocity and spatial offsets of the SMG Halpha components with sub-halo offsets in the Millennium simulation database we infer an average halo mass for SMGs of 13<log(M[h^-1Mo])<14. Finally we explore the relationship between the velocity dispersion and star formation intensity within the SMGs, finding the gas motions are consistent with the Kennicutt-Schmidt law and a range of extinction corrections, although might also be driven by the tidal torques from merging or even the star formation itself.

Subhalo abundance matching (SHAM) is a technique for populating simulated dark matter distributions with galaxies, assuming a monotonic relation between a galaxy's stellar mass or luminosity and the mass of its parent dark matter halo or subhalo. We examine the accuracy of SHAM in two cosmological SPH simulations, one of which includes momentum-driven winds. The SPH simulations indeed show a nearly monotonic relation between stellar mass and halo mass provided that, for satellite galaxies, we use the mass of the subhalo at the epoch when it became a satellite. In each simulation, the median relation for central and satellite galaxies is nearly identical, though a somewhat larger fraction of satellites are outliers. SHAM-assigned masses (at z=0-2), luminosities (R-band at z=0), or star formation rates (at z=2) have a 68% scatter of 0.09-0.15 dex relative to the true simulation values. When we apply SHAM to the subhalo population of collisionless N-body simulation with the same initial conditions as the SPH runs, we find generally good agreement for the halo occupation distributions and halo radial profiles of galaxy samples defined by thresholds in stellar mass. However, because a small fraction of SPH galaxies suffer severe stellar mass loss after becoming satellites, SHAM slightly overpopulates high mass halos; this effect is more significant for the wind simulation, which produces galaxies that are less massive and more fragile. SHAM recovers the two-point correlation function of the SPH galaxies in the no-wind simulation to better than 10% at scales 0.1 < r < 10 Mpc/h. For the wind simulation, agreement is better than 15% at r > 2 Mpc/h, but overpopulation of massive halos increases the correlation function by a factor of ~2.5 on small scales.

The circumgalactic medium (CGM) is fed by galaxy outflows and accretion of intergalactic gas, but its mass, heavy element enrichment, and relation to galaxy properties are poorly constrained by observations. In a survey of the outskirts of 42 galaxies with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope, we detected ubiquitous, large (150 kiloparsec) halos of ionized oxygen surrounding star-forming galaxies, but we find much less ionized oxygen around galaxies with little or no star formation. This ionized CGM contains a substantial mass of heavy elements and gas, perhaps far exceeding the reservoirs of gas in the galaxies themselves. It is a basic component of nearly all star-forming galaxies that is removed or transformed during the quenching of star formation and the transition to passive evolution.

The He II \lambda 1640 emission line has been suggested as a direct probe of Population III (Pop III) stars at high-redshift, since it can arise from highly energetic ionizing photons associated with hot, metal free stars. We use the HST WFC3/F130N IR narrowband filter to probe He II \lambda 1640 emission in galaxy IOK-1 at z=6.96. The sensitivity of this measurement is >5x deeper than for previous measurements. From this deep narrowband imaging, combined with broadband observations in the F125W and F160W filters, we find the He II flux to be 1.2+/- 1.0x 10^-18 ergs/s/cm^2, corresponding to a 1\sigma upper limit on the Pop III star formation rate (SFR) of ~ 0.5 M_sun/yr for the case of a Salpeter IMF with 50-500M_sun and mass loss. Given that the broadband measurements can be fit with a UV continuum spectral flux density of ~ 4.85x 10^-10x \lambda ^-2.46 ergs/s/cm^2/A, which corresponds to an overall SFR of ~16+/-2.6 M_sun/yr, massive Pop III stars represent < 6% of the total star formation. This measurement places the strongest limit yet on metal-free star formation at high redshift, although the exact conversion from He II luminosity to Pop III SFR is highly uncertain due to the unknown IMF, stellar evolution, and photoionization effects. Although we have not detected He II \lambda 1640 at more than the 1.2\sigma level, our work suggests that a > 3\sigma level detection is possible with JWST.

In the coming era of data-intensive science, it will be increasingly important to be able to seamlessly move between scientific results, the data analyzed in them, and the processes used to produce them. As observations, derived data products, publications, and object metadata are curated by different projects and archived in different locations, establishing the proper linkages between these resources and describing their relationships becomes an essential activity in their curation and preservation. In this paper we describe initial efforts to create a semantic knowledge base allowing easier integration and linking of the body of heterogeneous astronomical resources which we call the Virtual Observatory (VO). The ultimate goal of this effort is the creation of a semantic layer over existing resources, allowing applications to cross boundaries between archives. The proposed approach follows the current best practices in Semantic Computing and the architecture of the web, allowing the use of off-the-shelf technologies and providing a path for VO resources to become part of the global web of linked data.

The Taiwanese-American Occultation Survey (TAOS) project has collected more than a billion photometric measurements since 2005 January. These sky survey data-covering timescales from a fraction of a second to a few hundred days-are a useful source to study stellar variability. A total of 167 star fields, mostly along the ecliptic plane, have been selected for photometric monitoring with the TAOS telescopes. This paper presents our initial analysis of a search for periodic variable stars from the time-series TAOS data on one particular TAOS field, No. 151 (RA = 17$^{\rm h}30^{\rm m}6\fs$67, Dec = 27\degr17\arcmin 30\arcsec, J2000), which had been observed over 47 epochs in 2005. A total of 81 candidate variables are identified in the 3 square degree field, with magnitudes in the range 8 < R < 16. On the basis of the periodicity and shape of the lightcurves, 29 variables, 15 of which were previously unknown, are classified as RR Lyrae, Cepheid, delta Scuti, SX Phonencis, semi-regular and eclipsing binaries.

We present a high resolution (down to 0.18"), multi-transition imaging study of the molecular gas in the z = 4.05 submillimeter galaxy GN20. GN20 is one of the most luminous starburst galaxy known at z > 4, and is a member of a rich proto-cluster of galaxies at z = 4.05 in GOODS-North. We have observed the CO 1-0 and 2-1 emission with the VLA, the CO 6-5 emission with the PdBI Interferometer, and the 5-4 emission with CARMA. The H_2 mass derived from the CO 1-0 emission is 1.3 \times 10^{11} (\alpha/0.8) Mo. High resolution imaging of CO 2-1 shows emission distributed over a large area, appearing as partial ring, or disk, of ~ 10kpc diameter. The integrated CO excitation is higher than found in the inner disk of the Milky Way, but lower than that seen in high redshift quasar host galaxies and low redshift starburst nuclei. The VLA CO 2-1 image at 0.2" resolution shows resolved, clumpy structure, with a few brighter clumps with intrinsic sizes ~ 2 kpc. The velocity field determined from the CO 6-5 emission is consistent with a rotating disk with a rotation velocity of ~ 570 km s^{-1} (using an inclination angle of 45^o), from which we derive a dynamical mass of 3 \times 10^{11} \msun within about 4 kpc radius. The star formation distribution, as derived from imaging of the radio synchrotron and dust continuum, is on a similar scale as the molecular gas distribution. The molecular gas and star formation are offset by ~ 1" from the HST I-band emission, implying that the regions of most intense star formation are highly dust-obscured on a scale of ~ 10 kpc. The large spatial extent and ordered rotation of this object suggests that this is not a major merger, but rather a clumpy disk accreting gas rapidly in minor mergers or smoothly from the proto-intracluster medium. ABSTRACT TRUNCATED

We have analyzed the first 3.75 years of data from TAOS, the Taiwanese American Occultation Survey. TAOS monitors bright stars to search for occultations by Kuiper Belt Objects (KBOs). This dataset comprises 5e5 star-hours of multi-telescope photometric data taken at 4 or 5 Hz. No events consistent with KBO occultations were found in this dataset. We compute the number of events expected for the Kuiper Belt formation and evolution models of Pan & Sari (2005), Kenyon & Bromley (2004), Benavidez & Campo Bagatin (2009), and Fraser (2009). A comparison with the upper limits we derive from our data constrains the parameter space of these models. This is the first detailed comparison of models of the KBO size distribution with data from an occultation survey. Our results suggest that the KBO population is comprised of objects with low internal strength and that planetary migration played a role in the shaping of the size distribution.

We analyzed data accumulated during 2005 and 2006 by the Taiwan-American Occultation Survey (TAOS) in order to detect short-period variable stars (periods of <~ 1 hour) such as delta Scuti. TAOS is designed for the detection of stellar occultation by small-size Kuiper Belt Objects (KBOs) and is operating four 50cm telescopes at an effective cadence of 5Hz. The four telescopes simultaneously monitor the same patch of the sky in order to reduce false positives. To detect short-period variables, we used the Fast Fourier Transform algorithm (FFT) inasmuch as the data points in TAOS light-curves are evenly spaced. Using FFT, we found 41 short-period variables with amplitudes smaller than a few hundredths of a magnitude and periods of about an hour, which suggest that they are low-amplitude delta Scuti stars (LADS). The light-curves of TAOS delta Scuti stars are accessible online at the Time Series Center website (http://timemachine.iic.harvard.edu)

We present the results of a search for occultation events by objects at distances between 100 and 1000 AU in lightcurves from the Taiwanese-American Occultation Survey (TAOS). We searched for consecutive, shallow flux reductions in the stellar lightcurves obtained by our survey between 7 February 2005 and 31 December 2006 with a total of $\sim4.5\times10^{9}$ three-telescope simultaneous photometric measurements. No events were detected, allowing us to set upper limits on the number density as a function of size and distance of objects in Sedna-like orbits, using simple models.

We study the topology of reionization using accurate three-dimensional radiative transfer calculations post-processed on outputs from cosmological hydrodynamic simulations. In our simulations, reionization begins in overdense regions and then "leaks" directly into voids, with filaments reionizing last owing to their combination of high recombination rate and low emissivity. This result depends on the uniquely-biased emissivity field predicted by our prescriptions for star formation and feedback, which have previously been shown to account for a wide array of measurements of the post-reionization Universe. It is qualitatively robust to our choice of simulation volume, ionizing escape fraction, and spatial resolution (in fact it grows stronger at higher spatial resolution) even though the exact overlap redshift is sensitive to each of these. However, it weakens slightly as the escape fraction is increased owing to the reduced density contrast at higher redshift. We also explore whether our results are sensitive to commonly-employed approximations such as using optically-thin Eddington tensors or substantially altering the speed of light. Such approximations do not qualitatively change the topology of reionization. However, they can systematically shift the overlap redshift by up to $\Delta z\sim 0.5$, indicating that accurate radiative transfer is essential for computing reionization. Our model cannot simultaneously reproduce the observed optical depth to Thomson scattering and ionization rate per hydrogen atom at $z=6$, which could owe to numerical effects and/or missing early sources of ionization.

The nature of galaxy structures on large scales is a key observational prediction for current models of galaxy formation. The SDSS and 2dF galaxy surveys have revealed a number of structures on 40-150 h^-1 Mpc scales at low redshifts, and some even larger ones. To constrain galaxy number densities, luminosities, and stellar populations in large structures at higher redshift, we have investigated two sheet-like structures of galaxies at z=0.8 and 1.3 spanning 150 h^-1 comoving Mpc embedded in large quasar groups extending over at least 200 h^-1 Mpc. We present first results of an analysis of these sheet--like structures using two contiguous 1deg GALEX fields (FUV and NUV) cross-correlated with optical data from the Sloan Digital Sky Survey (SDSS). We derive a sample of 462 Lyman Break Galaxy (LBG) candidates coincident with the sheets. Using the GALEX and SDSS data, we show that the overall average spectral energy distribution of a LBG galaxy at z~1 is flat (in f_lambda) in the rest frame wavelength range from 1500A, to 4000A, implying evolved populations of stars in the LBGs. From the luminosity functions we get indications for overdensities in the two LQGs compared to their foreground regions. Similar conclusions come from the calculation of the 2-point correlation function, showing a 2sigma overdensity for the LBGs in the z~0.8 LQG on scales of 1.6 to 4.8 Mpc, indicating similar correlation scales for our LBG sample as their z~3 counterparts.

Observations of intergalactic neutral hydrogen can provide a wealth of information about structure and galaxy formation, potentially tracing accretion and feedback processes on Mpc scales. Below a column density of NHI ~ 10^19 cm-2, the "edge" or typical observational limit for HI emission from galaxies, simulations predict a cosmic web of extended emission and filamentary structures. We study the distribution of neutral hydrogen and its 21cm emission properties in a cosmological hydrodynamic simulation, to gain more insights into the distribution of HI below NHI ~ 10^19 cm-2. Such Lyman Limit systems are expected to trace out the cosmic web, and are relatively unexplored. Beginning with a 32 h^-1 Mpc simulation, we extract the neutral hydrogen component by determining the neutral fraction, including a post-processed correction for self-shielding based on the thermal pressure. We take into account molecular hydrogen, assuming an average density ratio Omega_H2 / Omega_HI = 0.3 at z = 0. The statistical properties of the HI emission are compared with observations, to assess the reliability of the simulation. The simulated HI distribution robustly describes the full column density range between NHI ~ 10^14 and NHI ~ 10^21 cm-2 and agrees very well with available measurements from observations. Furthermore there is good correspondence in the statistics when looking at the two-point correlation function and the HI mass function. The reconstructed maps are used to simulate observations of existing and future telescopes by adding noise and taking account of the sensitivity of the telescopes. The general agreement in statistical properties of HI suggests that neutral hydrogen as modeled in this hydrodynamic simulation is a fair representation of that in the Universe. (abridged)

We compare the properties of galaxies that form in a cosmological simulation without strong feedback to observations at z=0. We confirm previous findings that models without strong feedback overproduce the observed galaxy baryonic mass function, especially at the low and high mass extremes. Through post-processing we investigate what kinds of feedback would be required to reproduce observed galaxy masses and star formation rates. To mimic an extreme form of "preventive" feedback (e.g., AGN radio mode) we remove all baryonic mass that was originally accreted via "hot mode" from shock-heated gas. This does not bring the high mass end of the galaxy mass function into agreement with observations because much of the stellar mass in these systems formed at high redshift from baryons that originally accreted via "cold mode" onto lower mass progenitors. An efficient "ejective" feedback mechanism, such as supernova driven winds, must reduce the masses of these progenitors. Feedback must also reduce the masses of lower mass z=0 galaxies, which assemble at lower redshifts and have much lower star formation rates. If we monotonically re-map galaxy masses to reproduce the observed mass function, but retain the simulation's predicted star formation rates, we obtain fairly good agreement with the observed sequence of star-forming galaxies but fail to recover the observed population of passive, low star formation rate galaxies. Suppressing all hot mode accretion improves agreement for high mass galaxies but worsens the agreement at intermediate masses. Reproducing these z=0 observations requires a feedback mechanism that dramatically suppresses star formation in a fraction of galaxies, increasing with mass, while leaving star formation rates of other galaxies essentially unchanged.

The Taiwanese-American Occultation Survey (TAOS) operates four fully automatic telescopes to search for occultations of stars by Kuiper Belt Objects. It is a versatile facility that is also useful for the study of initial optical GRB afterglows. This paper provides a detailed description of the TAOS multi-telescope system, control software, and high-speed imaging.

We study the formation of galaxies in a (50 Mpc/h)^3 cosmological simulation (2x288^3 particles), evolved using the entropy conserving SPH code Gadget-2. Most of the baryonic mass in galaxies of all masses is originally acquired through filamentary "cold mode" accretion of gas that was never shock heated to its halo virial temperature, confirming the key feature of our earlier results obtained with a different SPH code (Keres et al. 2005). Atmospheres of hot, virialized gas develop in halos above ~2.5e11 Msun, a transition mass that is nearly constant from z=3 to z=0. Cold accretion persists in halos above the transition mass, especially at z>=2. It dominates the growth of galaxies in low mass halos at all times, and it is the main driver of the cosmic star formation history. Satellite galaxies have accretion rates similar to central galaxies of the same baryonic mass at high redshifts, but they have less accretion than comparable central galaxies at low redshift. Relative to our earlier results, the Gadget-2 simulations predict much lower rates of "hot mode" accretion from the virialized gas component of massive halos. At z<=1, typical hot accretion rates in halos above 5e12 Msun are below 1 Msun/yr, even though our simulation does not include AGN heating or other forms of "preventive" feedback. The inner density profiles of hot gas in these halos are shallow, with long associated cooling times. The cooling recipes typically used in semi-analytic models can overestimate the accretion rates in these halos by orders of magnitude, so such models may overemphasize the role of preventive feedback in producing observed galaxy masses and colors. A fraction of the massive halos develop cuspy profiles and significant cooling rates between z=1 and z=0, a redshift trend similar to the observed trend in the frequency of cooling flow clusters.

The star BD+29 1748 was resolved to be a close binary from its occultation by the asteroid 87 Sylvia on 2006 December 18 UT. Four telescopes were used to observe this event at two sites separated by some 80 km apart. Two flux drops were observed at one site, whereas only one flux drop was detected at the other. From the long-term variation of Sylvia, we inferred the probable shape of the shadow during the occultation, and this in turn constrains the binary parameters: the two components of BD+29 1748 have a projected separation of 0.097" to 0.110" on the sky with a position angle 104 deg to 107 deg. The asteroid was clearly resolved with a size scale ranging from 130 to 290 km, as projected onto the occultation direction. No occultation was detected for either of the two known moonlets of 87 Sylvia.

We examine the accretion and merger histories of central and satellite galaxies in a smoothed particle hydrodynamics (SPH) cosmological simulation that resolves galaxies down to 7e9 M[Sun]. Most friends-of-friends halos in the simulation have a distinct central galaxy, typically two to five times more massive than the most massive satellite. As expected, satellites have systematically higher assembly redshifts than central galaxies of the same baryonic mass, and satellites in more massive halos form earlier. However, contrary to the simplest expectations, satellite galaxies continue to accrete gas and convert it to stars; the gas accretion declines steadily over a period of 0.5-1 Gyr after the satellite halo merges with a larger parent halo. Satellites in a cluster mass halo eventually begin to lose baryonic mass. Since z=1, 27% of central galaxies (above 3e10 M[Sun]) and 22% of present-day satellite galaxies have merged with a smaller system above a 1:4 mass ratio; about half of the satellite mergers occurred after the galaxy became a satellite and half before. In effect, satellite galaxies can remain "central" objects of halo substructures, with continuing accretion and mergers, making the transition in assembly histories and physical properties a gradual one. Implementing such a gradual transformation in semi-analytic models would improve their agreement with the observed colour distributions of satellite galaxies in groups and with the observed colour dependence of galaxy clustering.

We introduce a new code for computing time-dependent continuum radiative transfer and non-equilibrium ionization states in static density fields with periodic boundaries. Our code solves the moments of the radiative transfer equation, closed by an Eddingtion tensor computed using a long characteristics method. We show that pure (i.e., not source-centered) short characteristics and the optically-thin approximation are inappropriate for computing Eddington factors for the problem of cosmological reionization. We evolve the non-equilibrium ionization field via an efficient and accurate (errors <1%) technique that switches between fully implicit or explicit finite-differencing depending on whether the local timescales are long or short compared to the timestep. We tailor our code for the problem of cosmological reionization. In tests, the code conserves photons, accurately treats cosmological effects, and reproduces analytic Stromgren sphere solutions. Its chief weakness is that the computation time for the long characteristics calculation scales relatively poorly compared to other techniques (t_{LC} \propto N_{cells}^1.5); however, we mitigate this by only recomputing the Eddington tensor when the radiation field changes substantially. Our technique makes almost no physical approximations, so it provides a way to benchmark faster but more approximate techniques. It can readily be extended to evolve multiple frequencies, though we do not do so here. Finally, we note that our method is generally applicable to any problem involving the transfer of continuum radiation through a periodic volume.

Results from the first two years of data from the Taiwanese-American Occultation Survey (TAOS) are presented. Stars have been monitored photometrically at 4 Hz or 5 Hz to search for occultations by small (~3 km) Kuiper Belt Objects (KBOs). No statistically significant events were found, allowing us to present an upper bound to the size distribution of KBOs with diameters 0.5 km < D < 28 km.

This is the first measurement and detection of coherence in the intergalactic medium (IGM) at substantially high redshift (z~3.8) and on large physical scales (~2.5 h^-1 Mpc). We perform the measurement by presenting new observations from Keck LRIS of the high redshift quasar pair PC 1643+4631A, B and their Ly-alpha absorber coincidences. This experiment extends multiple sightline quasar absorber studies to higher redshift, higher opacity, larger transverse separation, and into a regime where coherence across the IGM becomes weak and difficult to detect. We fit 222 discrete Ly-alpha absorbers to sightline A and 211 to sightline B. Relative to a Monte Carlo pairing test (using symmetric, nearest neighbor matching) the data exhibit a 4sigma excess of pairs at low velocity splitting (<150 km/s), thus detecting coherence on transverse scales of ~2.5 h^-1 Mpc. We use spectra extracted from an SPH simulation to analyze symmetric pair matching, transmission distributions as a function of redshift and compute zero-lag cross-correlations to compare with the quasar pair data. The simulations agree with the data with the same strength (~4sigma) at similarly low velocity splitting above random chance pairings. In cross-correlation tests, the simulations agree when the mean flux (as a function of redshift) is assumed to follow the prescription given by Kirkman et al. (2005). While the detection of flux correlation (measured through coincident absorbers and cross-correlation amplitude) is only marginally significant, the agreement between data and simulations is encouraging for future work in which even better quality data will provide the best insight into the overarching structure of the IGM and its understanding as shown by SPH simulations.

We derive a physical model for the observed relations between star formation rate (SFR) and molecular line (CO and HCN) emission in galaxies, and show how these observed relations are reflective of the underlying star formation law. We do this by combining 3D non-LTE radiative transfer calculations with hydrodynamic simulations of isolated disk galaxies and galaxy mergers. We demonstrate that the observed SFR-molecular line relations are driven by the relationship between molecular line emission and gas density, and anchored by the index of the underlying Schmidt law controlling the SFR in the galaxy. Lines with low critical densities (e.g. CO J=1-0) are typically thermalized and trace the gas density faithfully. In these cases, the SFR will be related to line luminosity with an index similar to the Schmidt law index. Lines with high critical densities greater than the mean density of most of the emitting clouds in a galaxy (e.g. CO J=3-2, HCN J=1-0) will have only a small amount of thermalized gas, and consequently a superlinear relationship between molecular line luminosity and mean gas density. This results in a SFR-line luminosity index less than the Schmidt index for high critical density tracers. One observational consequence of this is a significant redistribution of light from the small pockets of dense, thermalized gas to diffuse gas along the line of sight, and prodigious emission from subthermally excited gas. At the highest star formation rates, the SFR-Lmol slope tends to the Schmidt index, regardless of the molecular transition. The fundamental relation is the Kennicutt-Schmidt law, rather than the relation between SFR and molecular line luminosity. We use these results to make imminently testable predictions for the SFR-molecular line relations of unobserved transitions.

The isotropy of the Lyman-alpha forest in real-space uniquely provides a measurement of cosmic geometry at z > 2. The angular diameter distance for which the correlation function along the line of sight and in the transverse direction agree corresponds to the correct cosmological model. However, the Lyman-alpha forest is observed in redshift-space where distortions due to Hubble expansion, bulk flows, and thermal broadening introduce anisotropy. Similarly, a spectrograph's line spread function affects the autocorrelation and cross-correlation differently. In this the second paper of a series on using the Lyman-alpha forest observed in pairs of QSOs for a new application of the Alcock-Paczynski (AP) test, these anisotropies and related sources of potential systematic error are investigated with cosmological hydrodynamic simulations. Three prescriptions for galactic outflow were compared and found to have only a marginal effect on the Lyman-alpha flux correlation (which changed by at most 7% with use of the currently favored variable-momentum wind model vs. no winds at all). An approximate solution for obtaining the zero-lag cross-correlation corresponding to arbitrary spectral resolution directly from the zero-lag cross-correlation computed at full-resolution (good to within 2% at the scales of interest) is presented. Uncertainty in the observationally determined mean flux decrement of the Lyman-alpha forest was found to be the dominant source of systematic error; however, this is reduced significantly when considering correlation ratios. We describe a simple scheme for implementing our results, while mitigating systematic errors, in the context of a future application of the AP test.

We assess the impact of starburst and AGN feedback-driven winds on the CO emission from galaxy mergers, and, in particular, search for signatures of these winds in the simulated CO morphologies and emission line profiles. We do so by combining a 3D non-LTE molecular line radiative transfer code with smoothed particle hydrodynamics (SPH) simulations of galaxy mergers that include prescriptions for star formation, black hole growth, a multiphase interstellar medium (ISM), and the winds associated with star formation and black hole growth. Our main results are: (1) Galactic winds can drive outflows of masses ~10^8-10^9 Msun which may be imaged via CO emission line mapping. (2) AGN feedback-driven winds are able to drive imageable CO outflows for longer periods of time than starburst-driven winds owing to the greater amount of energy imparted to the ISM by AGN feedback compared to star formation. (3) Galactic winds can control the spatial extent of the CO emission in post-merger galaxies, and may serve as a physical motivation for the sub-kiloparsec scale CO emission radii observed in local advanced mergers. (4) Secondary emission peaks at velocities greater than the circular velocity are seen in the CO emission lines in all models. In models with winds, these high velocity peaks are seen to preferentially correspond to outflowing gas entrained in winds, which is not the case in the model without winds. The high velocity peaks seen in models without winds are typically confined to velocity offsets (from the systemic) < 1.7 times the circular velocity, whereas the models with AGN feedback-driven winds can drive high velocity peaks to ~2.5 times the circular velocity.

(Abridged) Using cosmological hydrodynamic simulations in combination with analytic modeling, we show that the galaxy stellar mass-metallicity relation (MZR) provides strong constraints on galactic outflows across cosmic time. We compare three outflow models: No outflows, a "constant wind" (cw) model that emulates the popular Dekel & Silk (1986) scenario, and a ``momentum-driven wind" (vzw) model. We find that only the vzw scaling reproduces the observed z~2 MZR's slope, amplitude, and scatter. Comparing our fully three-dimensional simulations with a simple one-zone chemical evolution model, we find that the MZR can be understood in terms of three parameters: (1) The equilibrium metallicity Z_eq=y*SFR/ACC (where y=true yield), reflecting the enrichment balance between star formation rate SFR and gas accretion rate ACC; (2) the dilution time t_d=M_g/ACC, representing the timescale for a galaxy to return to Z_eq after a metallicity-perturbing interaction; and (3) the blowout mass M_blowout, which is the galaxy stellar mass above which winds can escape its halo. When outflows with mass loading factor MLF are present, galaxies below M_blowout obey Z_eq = y/(1+MLF), while above M_blowout, Z_eq->y. Our cw model has M_blowout ~ 10^10 M_sun, which yields a sharp upturn in the MZR above this scale and a flat MZR with large scatter below it, in strong disagreement with observations. Our vzw model reproduces the observed Z_g \propto M_*^0.3 because Z_eq \propto MLF^-1 \propto M_*^1/3 when MLF >> 1 (i.e. at low masses). The flattening of the MZR at M_* > 10^10.5 M_sun observed by Tremonti et al. (2004) reflects the mass scale where MLF~1, rather than a characteristic wind speed. The tight observed MZR scatter is ensured when t_d<1 dynamical time, which is only satisified at all masses and epochs in our momentum-driven wind model.

We investigate the nature of CO emission from z~6 quasars by combining non-LTE radiative transfer calculations with merger-driven models of z~6 quasar formation that arise naturally in LCDM cosmological simulations. We consider four model quasars formed in 10^12-10^13 M_sun halos from different merging histories. Our main results follow. Owing to massive starbursts and funneling of dense gas into the nuclear regions of merging galaxies, the CO is highly excited and the flux density peaks between J=5-8. The CO morphology of z~6 quasars often exhibits multiple emission peaks which arise from H2 concentrations which have not yet fully coalesced. Quasars at z~6 display a large range of sightline dependent line widths such that the lines are narrowest when the rotating H2 gas associated with the quasar is viewed face-on (when L_B is largest), and broadest when the gas is seen edge-on (when L_B is lowest). Thus for all models selection effects exist such that quasars selected for optical luminosity are preferentially face-on which may result in detected CO line widths narrower than the median. The sightline averaged line width is reflective of the circular velocity (V_c) of the host halo, and ranges from sigma~300-650 km/s. For optically selected QSOs, 10-25% (halo-mass dependant) of sightlines have narrow line widths compatible with the sole CO detection at z>6, J1148+5251. When accounting for both the temporal evolution of CO line widths, as well as the redshift evolution of halo V_c, these models self-consistently account for the CO line widths of both z~2 sub-mm galaxies and QSO's. Finally, the dynamical mass derived from the sightline averaged line widths provides a good estimate of the total mass, and allows for a stellar bulge and SMBH consistent with the local M_BH-M_bulge relation. [abridged]

We measure the 2-point correlation function, xi(AG), between galaxies and quasar absorption line systems at z<1, using the dataset of Morris & Jannuzi (2006; paper I) on 16 lines of sight (LOS) with UV spectroscopy and galaxy multi-object spectroscopy. The measurements are made in 2-D redshift space out to pi=20/h Mpc (comoving) along the LOS and out to 2/h Mpc projected; as a function of HI column density in the range N(HI) = 1E13-1E19 cm^-2, for CIV systems and as a function of galaxy spectral type. This extends the absorber-galaxy pair count analysis of paper I. We find that the peak amplitude of xi(AG) at the smallest separations increases slowly as the lower limit on N(HI) is increased from 1E13 to 1E16 cm^-2, and then jumps sharply (albeit with substantial uncertainty) at N(HI)>1E17 cm^-2. For CIV absorbers, the peak strength of xi(AG) is comparable to that of HI absorbers with N(HI)>1E16.5 cm^-2. We do not reproduce the differences reported by Chen et al. between 1-D xi(AG) measurements using galaxy sub-samples of different spectral types, but the full impact of systematic differences in our samples is hard to quantify. We compare the observations with smoothed particle hydrodynamical (SPH) simulations and discover that in the observations xi(AG) is more concentrated to the smallest separations than in the simulations. The latter also display a `finger of god' elongation of xi(AG) along the LOS in redshift space, which is absent from our data, but similar to that found by Ryan-Weber for the cross-correlation of quasar absorbers and HI-emission-selected galaxies. The physical origin of these `fingers of god' is unclear and we highlight several possible areas for further investigation.

The Taiwan-America Occultation Survey (TAOS) aims to determine the number of small icy bodies in the outer reach of the Solar System by means of stellar occultation. An array of 4 robotic small (D=0.5 m), wide-field (f/1.9) telescopes have been installed at Lulin Observatory in Taiwan to simultaneously monitor some thousand of stars for such rare occultation events. Because a typical occultation event by a TNO a few km across will last for only a fraction of a second, fast photometry is necessary. A special CCD readout scheme has been devised to allow for stellar photometry taken a few times per second. Effective analysis pipelines have been developed to process stellar light curves and to correlate any possible flux changes among all telescopes. A few billion photometric measurements have been collected since the routine survey began in early 2005. Our preliminary result of a very low detection rate suggests a deficit of small TNOs down to a few km size, consistent with the extrapolation of some recent studies of larger (30--100 km) TNOs.

We compare the galaxy population of an SPH simulation to those predicted by the GalICS semi-analytic model and a stripped down version without supernova and AGN feedback. The SPH simulation and the no-feedback GalICS model make similar predictions for the baryonic mass functions of galaxies and for the dependence of these mass functions on environment and redshift. The two methods also make similar predictions for the galaxy content of dark matter haloes as a function of halo mass and for the gas accretion history of galaxies. Both the SPH and no-feedback GalICS models predict a bimodal galaxy population at z=0. The "red'' sequence of gas poor, old galaxies is populated mainly by satellite systems while, contrary to observations, the central galaxies of massive haloes lie on the "blue'' star-forming sequence as a result of continuing hot gas accretion at late times. Furthermore, both models overpredict the observed baryonic mass function, especially at the high mass end. In the full GalICS model, supernova-driven outflows reduce the masses of low and intermediate mass galaxies by about a factor of two. AGN feedback suppresses gas cooling in large haloes, producing a sharp cut-off in the baryonic mass function and moving the central galaxies of these massive haloes to the red sequence. Our results imply that the observational failings of the SPH simulation and the no-feedback GalICS model are a consequence of missing input physics rather than computational inaccuracies, that truncating gas accretion by satellite galaxies automatically produces a bimodal galaxy distribution with a red sequence, but that explaining the red colours of the most massive galaxies requires a mechanism like AGN feedback that suppresses the accretion onto central galaxies in large haloes.

We study the effects of feedback from active galactic nuclei (AGN) on emission from molecular gas in galaxy mergers by combining hydrodynamic simulations which include black holes with a three-dimensional, non-local thermodynamic equilibrium (LTE) radiative transfer code. We find that molecular clouds entrained in AGN winds produce an extended CO morphology with significant off-nuclear emission, which may be detectable via contour mapping. Furthermore, kinematic signatures of these molecular outflows are visible in emission line profiles when the outflow has a large line of sight velocity. Our results can help interpret current and upcoming observations of luminous infrared galaxies, as well as provide a detailed test of subresolution prescriptions for supermassive black hole growth in galaxy-scale hydrodynamic simulations.

We investigate the effect of embedded active galactic nuclei (AGN) in galaxy mergers on the CO molecular line emission by combining non-local thermodynamic equilibrium (LTE) radiative transfer calculations with hydrodynamic simulations. We find that AGN feedback energy in gas rich galaxy mergers can contribute to large molecular outflows which may be detectable via velocity-integrated emission contour maps, as well as through kinematic features in the emission line profiles.

We construct merger trees for galaxies identified in a cosmological hydrodynamical simulation and use them to characterize predicted merger rates as a function of redshift, galaxy mass, and merger mass ratio. At z=0.3, we find a mean rate of 0.054 mergers per galaxy per Gyr above a 1:2 mass ratio threshold for massive galaxies (baryonic mass above 6.4E10 \Msun), but only 0.018 / Gyr for lower mass galaxies. The mass ratio distribution is proportional to R_merg^-1.2 for the massive galaxy sample, so high mass mergers dominate the total merger growth rate. The predicted rates increase rapidly with increasing redshift, and they agree reasonably well with observational estimates. A substantial fraction of galaxies do not experience any resolved mergers during the course of the simulation, and even for the high mass sample only 50% of galaxies experience a greater than 1:4 merger since z=1. Typical galaxies thus have fairly quiescent merger histories. We assign bulge-to-disk ratios to simulated galaxies by assuming that mergers above a mass ratio threshold R_major convert stellar disks into spheroids. With R_major of 1:4, we obtain a fairly good match to the observed dependence of early-type fraction on galaxy mass. However, the predicted fraction of truly bulge-dominated systems (f_bulge} > 0.8) is small, and producing a substantial population of bulge-dominated galaxies may require a mechanism that shuts off gas accretion at late times and/or additional processes (besides major mergers) for producing bulges.

We present STIS 7km/s resolution data of the Ly-a forest toward PKS0405-123 (z=0.574). We make two samples. The strong sample (60 systems) has column density logNHI>13.3 over 0.002<z<0.423, and the weak sample (44) has logNHI>13.1 at 0.020<z<0.234; 7 absorbers show metals, all with OVI, but often offset in velocity from HI. The strong sample Doppler parameters show <b>=47+-22km/s; for the weak, <b>=44+-21km/s. Line blending and s/n effects likely inflate b. The redshift density dN/dz is consistent with previous, lower resolution results for logNHI>14.0. For 13.1<logNHI<14.0, we find a dN/dz overdensity of ~0.2-0.3 dex at 0.127<z<0.234, which we believe arises from cosmic variance. We find Ly-a clustering over Delta v<250km/s for logNHI>13.3, consistent with a numerical model. There is a void in the strong sample at 0.0320<z<0.0814. We detect radial velocity correlations for Delta v<250km/s between Ly-a absorbers and 39 galaxies in a 10' field. The correlation grows with minimum NHI, peaking for logNHI ~> 13.5-14.0. The strength is similar to that of the galaxy-galaxy correlation for our sample, implying that such Ly-a systems have mass log M/Msun= 11.3+1.0-0.6. Including lower column density systems in the sample shows correlations only out to Delta v <125km/s, as would be expected for smaller density perturbations. Galaxy counts and local HI column density appear correlated, peaking at ~5000km/s. We also give column densities for Galactic species. (Abridged)

We analyze the flux power spectrum and its covariance using simulated Lyman alpha forests. We find that pseudo-hydro techniques are good approximations of hydrodynamical simulations at high redshift. However, the pseudo-hydro techniques fail at low redshift because they are insufficient for characterizing some components of the low-redshift intergalactic medium, notably the warm-hot intergalactic medium. Hence, to use the low-redshift Lyman alpha flux power spectrum to constrain cosmology, one would need realistic hydrodynamical simulations. By comparing one-dimensional mass statistics with flux statistics, we show that the nonlinear transform between density and flux quenches the fluctuations so that the flux power spectrum is much less sensitive to cosmological parameters than the one-dimensional mass power spectrum. The covariance of the flux power spectrum is nearly Gaussian. As such, the uncertainties of the underlying mass power spectrum could still be large, even though the flux power spectrum can be precisely determined from a small number of lines of sight.

We study the physical and photometric properties of galaxies at z=4 in cosmological hydrodynamic simulations of a lambda-CDM universe. We focus on galaxies satisfying the GOODS "B-dropout" criteria. Our goals are: (1) to study the nature of high-redshift galaxies; (2) to test the simulations against published measurements of high-redshift galaxies; (3) to find relations between photometric measurements by HST/ACS (0.4 -- 1 micron) and Spitzer/IRAC (3.6 -- 8 micron) and the intrinsic physical properties of GOODS "B-dropouts" such as stellar mass, stellar age, dust reddening, and star-formation rate; and (4) to assess how representative the GOODS survey is at this epoch. Our simulations predict that high-redshift galaxies show strong correlations in star formation rate versus stellar mass, and weaker correlations versus environment and age, such that GOODS galaxies are predicted to be the most massive, most rapidly star-forming galaxies at that epoch, living preferentially in dense regions. The simulated rest-frame UV luminosity function (LF) and integrated luminosity density are in broad agreement with observations at z~4. The predicted rest-frame optical (observed 3.6 micron) LF is similar to the rest-frame UV LF, shifted roughly one magnitude brighter. We predict that GOODS detects less than 50% of the total stellar mass density formed in galaxies more massive than 10^8.7 M_sun by z=4, mainly because of brightness limits in the HST/ACS bands. The most rapidly star forming galaxies in our simulations have rates exceeding 1000 M_sun yr^-1, similar to observed sub-mm galaxies. The star formation rates of these galaxies show at most a mild excess (2--3x) over the rates that would be expected for their stellar mass. Whether these bright galaxies would be observable as LBGs depends on the uncertain effects of dust reddening.

We examine the temperature history of gas accreted by forming galaxies in SPH simulations. About half the gas shock heats to roughly the virial temperature of the galaxy potential well before cooling, condensing, and forming stars, but the other half radiates its acquired gravitational energy at much lower temperatures, typically T<10^5 K, and the histogram of maximum gas temperatures is clearly bimodal. The "cold mode" of gas accretion dominates for low mass galaxies (M_baryon < 10^{10.3}Msun or M_halo < 10^{11.4}Msun), while the conventional "hot mode" dominates the growth of high mass systems. Cold accretion is often directed along filaments, allowing galaxies to efficiently draw gas from large distances, while hot accretion is quasi-spherical. The galaxy and halo mass dependence leads to redshift and environment dependence of cold and hot accretion rates, with cold mode dominating at high redshift and in low density regions today, and hot mode dominating in group and cluster environments at low redshift. Star formation rates closely track accretion rates, and we discuss the physics behind the observed environment and redshift dependence of galactic scale star formation. If we allowed hot accretion to be suppressed by conduction or AGN feedback, then the simulation predictions would change in interesting ways, perhaps resolving conflicts with the colors of ellipticals and the cutoff of the galaxy luminosity function. The transition between cold and hot accretion at M_h ~ 10^{11.4}Msun is similar to that found by Birnboim & Dekel (2003) using 1-d simulations and analytic arguments. The corresponding baryonic mass is tantalizingly close to the scale at which Kauffmann et al. (2003) find a marked shift in galaxy properties. We speculate on connections between these theoretical and observational transitions.

The halo occupation distribution (HOD) describes the relation between galaxies and dark matter at the level of individual dark matter halos. The properties of galaxies residing at the centers of halos differ from those of satellite galaxies because of differences in their formation histories. Using a smoothed particle hydrodynamics (SPH) simulation and a semi-analytic (SA) galaxy formation model, we examine the separate contributions of central and satellite galaxies to the HOD, more specifically to the probability P(N|M) that a halo of virial mass M contains N galaxies of a particular class. In agreement with earlier results for dark matter subhalos, we find that the mean occupation function <N> for galaxies above a baryonic mass threshold can be approximated by a step function for central galaxies plus a power law for satellites, and that the distribution of satellite numbers is close to Poisson at fixed halo mass. For galaxy samples defined by different baryonic mass thresholds, there is a nearly linear relation between the minimum halo mass Mmin required to host a central galaxy and the mass M1 at which an average halo hosts one satellite, with M1 ~ 14 Mmin (SPH) or M1 ~ 18 Mmin (SA). The mean occupation number of young galaxies exhibits a local minimum at M ~ 10 Mmin where halos are too massive to host a young central galaxy but not massive enough to host satellites. We show that the conditional galaxy mass function at fixed halo mass cannot be described by a Schechter function because central galaxies produce a "bump" at high masses. We suggest parameterizations for the HOD and the conditional luminosity function that can be used to model observed galaxy clustering. Many of our predictions are in good agreement with recent results inferred from clustering in the Sloan Digital Sky Survey.

[Abridged] We test the reliability of a method to measure the mean halo mass of Damped Ly-alpha absorbers (DLAs). The method is based on measuring the ratio of the cross-correlation between DLAs and galaxies to the auto-correlation of the galaxies themselves ($w_{\rm dg}/w_{\rm gg}$), which is (in linear theory) the ratio of their bias factor. This is shown to be true irrespective of the galaxy redshift distribution, provided that one uses the same galaxies for the two correlation functions. The method is applicable to all redshifts. Here, we focus on z=3 DLAs and we demonstrate that the method robustly constrains the mean DLA halo mass using smoothed particle hydrodynamics (SPH) cosmological simulations. If we use the bias formalism of Mo & White with the DLA and galaxy mass distributions of these simulations, we predict a bias ratio of 0.771. Direct measurement from the simulations of $w_{\rm dg}/w_{\rm gg}$ st yields a ratio of 0.73+/-0.08, in excellent agreement with that prediction. Equivalently, inverting the measured correlation ratio to infer a mean DLA halo mass yields (log. averaging, in solar units) <log(M_DLA)> =11.13+/-013, in excellent agreement with the true value in the simulations: 11.16. The cross- correlation method thus appears to yield a robust estimate of the average host halo mass even though the DLAs and the galaxies occupy a broad mass spectrum of halos, and massive halos contain multiple galaxies with DLAs. We show that the inferred mean DLA halo mass is independent of the galaxy sub-sample used, i.e. the cross-correlation technique is also reliable. Our results imply that the cross-correlation length between DLAs and LBGs is predicted to be, at most, 2.85 Mpc. Future observations will soon distinguish models in which DLAs are in low mass halos from those in which DLAs are in massive halos.

We present a methodology to discover outliers in catalogs of periodic light-curves. We use cross-correlation as measure of ``similarity'' between two individual light-curves and then classify light-curves with lowest average ``similarity'' as outliers. We performed the analysis on catalogs of variable stars of known type from the MACHO and OGLE projects and established that our method correctly identifies light-curves that do not belong to those catalogs as outliers. We show how our method can scale to large datasets that will be available in the near future such as those anticipated from Pan-STARRS and LSST.

We study the soft X-ray emission (0.1 - 1 keV) from the Warm-Hot Intergalactic Medium (WHIM) in a hydrodynamic simulation of a Cold Dark Matter universe. Our main goal is to investigate how such emission can be explored with a combination of imaging and spectroscopy, and to motivate future X-ray missions. We first present high resolution images of the X-ray emission in several energy bands, in which emission from different ion species dominates. We pick three different areas to study the high resolution spectra of X-rays from the warm-hot IGM: (a) a galaxy group; (b) a filament and (c) an underluminous region. By taking into account the background X-ray emission from AGNs and foreground emission from the Galaxy, we compute composite X-ray spectra of the selected regions. We briefly investigate angular clustering of the soft-X-ray emission, finding a strong signal. Most interestingly, the combination of high spectral resolution and angular information allows us to map the emission from the WHIM in 3 dimensions. We cross-correlate the positions of galaxies in the simulation with this redshift map of emission and detect the presence of six different ion species (Ne IX, Fe XVII, O VII, O VIII, N VII, C VI) in the large-scale structure traced by the galaxies. Finally we show how such emission can be detected and studied with future X-ray satellites, with particular attention to a proposed mission, the Missing Baryon Explorer, or MBE. We present simulated observations of the WHIM gas with MBE.

We present results for QSO PKS 0405-123 (z=0.574, V=14.9), as part of a STIS Investigation Definition Team (IDT) key project to study weak Ly-alpha forest systems at low z. We detect 59 (47) Ly-alpha absorbers at 4.0 sigma significance to an 80% completeness limit of column density log N(HI)=13.3 (13.1) for Doppler parameter V_Dop=40 km/s over 0.002<z<0.423 (0.020<z<0.234). We find 4 intervening O VI systems, useful for studies of hot intergalactic gas. We do not distinguish between metal and Ly-alpha-only systems in the following analysis. The redshift density is consistent with previous measurements for log N(HI)>=14.0, but exhibits twice as many systems at 13.1<log N(HI)<14.0 compared to the mean number density of lines at z<0.07 toward 15 extragalactic objects. The difference possibly arises from cosmic variance. The Doppler parameter distribution has <V_Dop>=48 +- 21 km/s; line blending possibly inflates the value. We find evidence for Ly-alpha-Ly-alpha clustering in our sample on a scale of Delta v<=250 km/s, and there is evidence for a void at 0.032<z<0.081 with probability of occurrance P=0.0005. We find line-of-sight velocity correlations of up to 250 km/s between Ly-alpha absorbers with log N(HI)>=13.1 and 45 galaxies taken from the literature and unpublished data at 0<z<0.47; the transverse distances cover up to 1.5 /h_70 Mpc in the local frame. The Ly-alpha-galaxy clustering is stronger for higher log(N(HI) systems.

Not the way one might have thought. In hydrodynamic simulations of galaxy formation, some gas follows the traditionally envisioned route, shock heating to the halo virial temperature before cooling to the much lower temperature of the neutral ISM. But most gas enters galaxies without ever heating close to the virial temperature, gaining thermal energy from weak shocks and adiabatic compression, and radiating it just as quickly. This ``cold mode'' accretion is channeled along filaments, while the conventional, ``hot mode'' accretion is quasi-spherical. Cold mode accretion dominates high redshift growth by a substantial factor, while at z<1 the overall accretion rate declines and hot mode accretion has greater relative importance. The decline of the cosmic star formation rate at low z is driven largely by geometry, as the typical cross section of filaments begins to exceed that of the galaxies at their intersections.

We analyze the evolution of energy density fluctuations in cosmological scenarios with a mixture of cold dark matter and quintessence, in which the quintessence field is modeled by a constant equation of state. We obtain analytic expressions for the time evolution of the quintessence perturbations in models with light fields. The fluctuations behave analogously to a driven harmonic oscillator, where the driving term arises from the inhomogeneities in the surrounding cosmological fluid. We demonstrate that the homogeneous solution, determined by the initial conditions, is completely sub-dominant to the inhomogeneous solution for physically realistic scenarios. Thus we show that the cosmic microwave background anisotropy predicted for such models is highly insensitive to the initial conditions in the quintessence field.

We present preliminary results on the low-redshift Lyman alpha forest as based on STIS spectra of 3C 273. A total of 121 intergalactic Lyman alpha-absorbing systems were detected, of which 60 are above the 3.5 sigma completness limit, log N(HI)~12.3. The median Doppler parameter, b=27 km/s, is similar to that seen at high redshift. However the distribution of HI column densities (dN/dN(HI) propto N(HI)^-beta) has a steeper slope, beta = 2.02 +- 0.21, than is seen at high redshift. Overall, the observed N(HI)-b distribution is consistent with that derived from a Lambda CDM hydrodynamic simulation.

We use hydrodynamic simulations to examine how the baryonic components of galaxies are assembled, focusing on the relative importance of mergers and smooth accretion in the formation of ~L_* systems. In our primary simulation, which models a (50\hmpc)^3 comoving volume of a Lambda-dominated cold dark matter universe, the space density of objects at our (64-particle) baryon mass resolution threshold, M_c=5.4e10 M_sun, corresponds to that of observed galaxies with L~L_*/4. Galaxies above this threshold gain most of their mass by accretion rather than by mergers. At the redshift of peak mass growth, z~2, accretion dominates over merging by about 4:1. The mean accretion rate per galaxy declines from ~40 M_sun/yr at z=2 to ~10 M_sun/yr at z=0, while the merging rate peaks later (z~1) and declines more slowly, so by z=0 the ratio is about 2:1. We cannot distinguish truly smooth accretion from merging with objects below our mass resolution threshold, but extrapolating our measured mass spectrum of merging objects, dP/dM ~ M^a with a ~ -1, implies that sub-resolution mergers would add relatively little mass. The global star formation history in these simulations tracks the mass accretion rate rather than the merger rate. At low redshift, destruction of galaxies by mergers is approximately balanced by the growth of new systems, so the comoving space density of resolved galaxies stays nearly constant despite significant mass evolution at the galaxy-by-galaxy level. The predicted merger rate at z<~1 agrees with recent estimates from close pairs in the CFRS and CNOC2 redshift surveys.

Spergel and Steinhardt have recently proposed the concept of dark matter with strong self-interactions as a means to address numerous discrepancies between observations of dark matter halos on subgalactic scales and the predictions of the standard collisionless dark matter picture. We review the motivations for this scenario and discuss some recent, successful numerical tests. We also discuss the possibility that the dark matter interacts strongly with ordinary baryonic matter, as well as with itself. We present a new analysis of the experimental constraints and re-evaluate the allowed range of cross-section and mass.

Some form of missing energy may account for the difference between the observed cosmic matter density and the critical density. Two leading candidates are a cosmological constant and quintessence (a time-varying, inhomogenous component with negative pressure). We show that an ideal, full-sky cosmic background anisotropy experiment may not be able to distinguish the two, even when non-linear effects due to gravitational lensing are included. Due to this ambiguity, microwave background experiments alone may not determine the matter density or Hubble constant very precisely. We further show that degeneracy may remain even after considering classical cosmological tests and measurements of large scale structure.

We examine the possibility that a significant component of the energy density of the universe has an equation-of-state different from that of matter, radiation or cosmological constant ($\Lambda$). An example is a cosmic scalar field evolving in a potential, but our treatment is more general. Including this component alters cosmic evolution in a way that fits current observations well. Unlike $\Lambda$, it evolves dynamically and develops fluctuations, leaving a distinctive imprint on the microwave background anisotropy and mass power spectrum.

We study metal line absorption of CIV, CII, SiIV, and NV at redshifts z=3.5 to z=2 within the framework of a cosmological model for the Lyman alpha forest, comparing the results of numerical simulations to recent observations by Songaila & Cowie (1996, SC). In agreement with Rauch, Haehnelt & Steinmetz (1996), we find that the observed mean value of the CIV/HI ratio around z=3 is reproduced if a uniform metallicity of [C/H] = -2.5 is assumed in our model, but that the observed scatter in this ratio is larger than predicted, implying a scatter in the metallicity of the absorbing systems of roughly an order of magnitude. We further argue that absorbers with relatively low column densities (log(N_HI) < 15) likely have a mean metallicity [C/H] less than -2.5, a result which is basically independent of the model considered, and similar to what is predicted in the simulations by Gnedin & Ostriker (1997). Our model predicts no significant evolution in the mean values of metal line column densities between z=3.5 and z=2. Comparison of the predicted and observed numbers of SiIV and NV systems suggests that the photoionizing background radiation field around z=3 is somewhat softer than that proposed by Haardt & Madau (1996). Our model does not account for the increase in the SiIV/CIV ratio at z=3.2 found by SC. While SC suggested that the increase could be explained by a softening of the radiation field due to HeII absorption at z > 3, such a modification does not raise the mean value of SiIV/CIV in our simulation because it shifts numerous weak SiIV systems to just above the detection limit, thus keeping the mean column density of observable SiIV systems low.