We present an analysis of the physical and dynamical states of two sets of EAGLE zoom simulations of galaxy haloes, one at high redshift ($z=2-3$) and the other at low redshift ($z=0$), with masses of $\approx 10^{12} M_{\odot}$. Our focus is how the circumgalactic medium (CGM) of these $L^*$ star-forming galaxies change over the last 10 Gyr. We find that the high-$z$ CGM is almost equally divided between the "cool" ($T<10^5$ K) and "hot" ($T\geq 10^5$ K) phases, while the low-$z$ hot CGM phase contains $5\times$ more mass. The high-$z$ hot CGM contains 60% more metals than the cool CGM, while the low-$z$ cool CGM contains 35% more metals than the hot CGM content. The metals are evenly distributed radially between the hot and cool phases throughout the high-$z$ CGM. At high $z$, the CGM volume is dominated by hot outflows, cool gas is mainly inflowing, but cool metals are flowing outward. At low $z$, the cool metals dominate the interior and the hot metals are more prevalent at larger radii. The low-$z$ cool CGM has tangential motions consistent with rotational support out to $0.2 R_{200}$, often exhibiting $r \approx 40$ kpc disc-like structures. The low-$z$ hot CGM has several times greater angular momentum than the cool CGM, and a more flattened radial density profile than the high-$z$ hot CGM. This study verifies that, just as galaxies demonstrate significant evolutionary stages over cosmic time, the gaseous haloes surrounding them also undergo considerable changes of their own both in physical characteristics of density, temperature and metallicity, and dynamic properties of velocity and angular momentum.

We use the EAGLE cosmological simulation to study the distribution of baryons, and FUV (O VI), EUV (Ne VIII) and X-ray (O VII, O VIII, Ne IX, and Fe XVII) line absorbers, around low-redshift galaxies and haloes of different masses. We work with column densities, but quantify their relation with equivalent widths by analysing virtual spectra. Halo gas dominates the highest column density absorption for X-ray lines, but lower-density gas contributes to strong UV absorption lines from O VI and Ne VIII. Of the O VII (O VIII) absorbers detectable in an Athena X-IFU blind survey, we find that 41 (56) per cent arise from haloes with $\mathrm{M}_{\mathrm{200c}} = 10^{12.0 - 13.5} \, \mathrm{M}_{\odot}$. We predict that the Athena X-IFU will be able to observe O VII (O VIII) covering fractions within 100 pkpc of 77 (46) per cent around $\mathrm{M}_{\star} = 10^{10.5 - 11.0} \, \mathrm{M}_{\odot}$ galaxies (59 (82) per cent for $\mathrm{M}_{\star} >10^{11.0} \, \mathrm{M}_{\odot}$), comparable to the covering fractions observable for O VI and Ne VIII by current UV spectrographs. The circumgalactic medium (CGM) contains more metals than the interstellar medium across halo masses. The ions we study here trace the warm-hot, volume-filling phase of the CGM, but are biased towards gas at temperatures corresponding to the collisional ionization peak for each ion, and towards high-metallicity gas. Gas well within the virial radius is mostly collisionally ionized, but around and beyond this radius, and for O VI, photo-ionization becomes significant. Virial-temperature collisional ionization equilibrium ion fractions are a good predictor of column density trends with halo mass, but haloes contain a larger diversity of ions than this single-temperature model predicts.

We simulate stacked observations of nearby hot X-ray coronae associated with galaxies in the EAGLE and Illustris-TNG hydrodynamic simulations. A forward modeling pipeline is developed to predict 4-year eROSITA observations and stacked image analysis, including the effects of instrumental and astrophysical backgrounds. We propose an experiment to stack z~0.01 galaxies separated by specific star-formation rate (sSFR) to examine how the hot (T>=10^6 K) circumgalactic medium (CGM) differs for high- and low-sSFR galaxies. The simulations indicate that the hot CGM of low-mass (M_*~10^{10.5} Msol), high-sSFR (defined as the top one-third ranked by sSFR) central galaxies will be detectable to a galactocentric radius r~30-50 kpc. Both simulations predict lower luminosities at fixed stellar mass for the low-sSFR galaxies (the lower third of sSFR) with Illustris-TNG predicting 3x brighter coronae around high-sSFR galaxies than EAGLE. Both simulations predict detectable emission out to r~150-200 kpc for stacks centered on high-mass (M_*~10^{11.0} Msol) galaxies, with EAGLE predicting brighter X-ray halos. The extended soft X-ray luminosity correlates strongly and positively with the mass of circumgalactic gas within the virial radius (f_{CGM}). Prior analyses of both simulations have established that f_{CGM} is reduced by expulsive feedback driven mainly by black hole growth, which quenches galaxy growth by inhibiting replenishment of the ISM. Both simulations predict that eROSITA stacks should not only conclusively detect and resolve the hot CGM around L^* galaxies for the first time, but provide a powerful probe of how the baryon cycle operates, for which there remains an absence of consensus between state-of-the-art simulations.

CMB-HD is a proposed ultra-deep (0.5 uk-arcmin), high-resolution (15 arcseconds) millimeter-wave survey over half the sky that would answer many outstanding questions in both fundamental physics of the Universe and astrophysics. This survey would be delivered in 7.5 years of observing 20,000 square degrees, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors.

We present a study of the galaxy environment of 9 strong HI+CIV absorption line systems ($16.2<{\rm log}(N({\rm HI}))<21.2$) spanning a wide range in metallicity at $z\sim4-5$, using MUSE integral field and X-Shooter spectroscopic data collected in a $z\approx 5.26$ quasar field. We identify galaxies within a 250 kpc and $\pm1000$ km s$^{-1}$ window for 6 out of the 9 absorption systems, with 2 of the absorption line systems showing multiple associated galaxies within the MUSE field of view. The space density of Ly$\alpha$ emitting galaxies (LAEs) around the HI and CIV systems is $\approx10-20$ times the average sky density of LAEs given the flux limit of our survey, showing a clear correlation between the absorption and galaxy populations. Further, we find that the strongest CIV systems in our sample are those that are most closely aligned with galaxies in velocity space, i.e. within velocities of $\pm500$ km s$^{-1}$. The two most metal poor systems lie in the most dense galaxy environments, implying we are potentially tracing gas that is infalling for the first time into star-forming groups at high redshift. Finally, we detect an extended Ly$\alpha$ nebula around the $z\approx 5.26$ quasar, which extends up to $\approx50$ kpc at the surface brightness limit of $3.8 \times 10^{-18}$ erg s$^{-1}$ cm$^{-2}$ arcsec$^{-2}$. After scaling for surface brightness dimming, we find that this nebula is centrally brighter, having a steeper radial profile than the average for nebulae studied at $z\sim3$ and is consistent with the mild redshift evolution seen from $z\approx 2$.

Many phenomenologically successful cosmological galaxy formation simulations employ kinetic winds to model galactic outflows, a crucial ingredient in obtaining predictions that agree with various observations. Yet systematic studies of how variations in kinetic wind scalings might alter observable galaxy properties are rare. Here we employ GADGET-3 simulations to study how the baryon cycle, stellar mass function, and other galaxy and CGM predictions vary as a function of the assumed outflow speed $v_w$ and the scaling of the mass loading factor $\eta$ with velocity dispersion $\sigma$. We design our fiducial model to reproduce the measured wind properties at 25% of the virial radius from the Feedback In Realistic Environments (FIRE) simulations. We find that a strong dependence of $\eta \sim \sigma^5$ in low mass haloes with $\sigma < 106\ \mathrm{km\ s^{-1}}$ is required to match the faint end of the stellar mass functions at $z > 1$. The wind speed also has a major impact, with faster winds significantly reducing wind recycling and heating more halo gas. Both effects result in less stellar mass growth in massive haloes and impact high ionization absorption in halo gas. We cannot simultaneously match the stellar content at $z=2$ and $z=0$ within a single model, suggesting that an additional feedback source such as AGN might be required in massive galaxies at lower redshifts, but the amount needed depends strongly on assumptions regarding the outflow properties. We run a 50 $\mathrm{Mpc/h}$, $2\times576^3$ simulation with our fiducial parameters and show that it matches a range of star-forming galaxy properties at $z\sim0-2$. In closing, the results from simulations of galaxy formation are much more sensitive to small changes in the feedback implementation than to the hydrodynamic technique.

We test and improve the numerical schemes in our smoothed particle hydrodynamics (SPH) code for cosmological simulations, including the pressure-entropy formulation (PESPH), a time-dependent artificial viscosity, a refined timestep criterion, and metal-line cooling that accounts for photoionisation in the presence of a recently refined Haardt \& Madau (2012) model of the ionising background. The PESPH algorithm effectively removes the artificial surface tension present in the traditional SPH formulation, and in our test simulations it produces better qualitative agreement with mesh-code results for Kelvin-Helmholtz instability and cold cloud disruption. Using a set of cosmological simulations, we examine many of the quantities we have studied in previous work. Results for galaxy stellar and HI mass functions, star formation histories, galaxy scaling relations, and statistics of the Ly$\alpha$ forest are robust to the changes in numerics and microphysics. As in our previous simulations, cold gas accretion dominates the growth of high-redshift galaxies and of low mass galaxies at low redshift, and recycling of winds dominates the growth of massive galaxies at low redshift. However, the PESPH simulation removes spurious cold clumps seen in our earlier simulations, and the accretion rate of hot gas increases by up to an order of magnitude at some redshifts. The new numerical model also influences the distribution of metals among gas phases, leading to considerable differences in the statistics of some metal absorption lines, most notably NeVIII.

We examine the connection between the properties of the circumgalactic medium (CGM) and the quenching and morphological evolution of central galaxies in the EAGLE and IllustrisTNG simulations. The simulations yield very different median CGM mass fractions, $f_{\rm CGM}$, as a function of halo mass, $M_{200}$, with low-mass haloes being significantly more gas-rich in IllustrisTNG than in EAGLE. Nonetheless, in both cases scatter in $f_{\rm CGM}$ at fixed $M_{200}$ is strongly correlated with the specific star formation rate and the kinematic morphology of central galaxies. The correlations are strongest for $\sim L^\star$ galaxies, corresponding to the mass scale at which AGN feedback becomes efficient. This feedback elevates the CGM cooling time, preventing gas from accreting onto the galaxy to fuel star formation, and thus establishing a preference for quenched, spheroidal galaxies to be hosted by haloes with low $f_{\rm CGM}$ for their mass. In both simulations, $f_{\rm CGM}$ correlates negatively with the host halo's intrinsic concentration, and hence with its binding energy and formation redshift, primarily because early halo formation fosters the rapid early growth of the central black hole (BH). This leads to a lower $f_{\rm CGM}$ at fixed $M_{200}$ in EAGLE because the BH reaches high accretion rates sooner, whilst in IllustrisTNG it occurs because the central BH reaches the mass threshold at which AGN feedback is assumed to switch from thermal to kinetic injection earlier. Despite these differences, there is consensus from these state-of-the-art simulations that the expulsion of efficiently-cooling gas from the CGM is a crucial step in the quenching and morphological evolution of central galaxies.

We use high-resolution hydrodynamical simulations run with the EAGLE model of galaxy formation to study the differences between the properties of - and subsequently the lensing signal from - subhaloes of massive elliptical galaxies at redshift 0.2, in Cold and Sterile Neutrino (SN) Dark matter models. We focus on the two 7 keV SN models that bracket the range of matter power spectra compatible with resonantly-produced SN as the source of the observed 3.5 keV line. We derive an accurate parametrisation for the subhalo mass function in these two SN models relative to CDM, as well as the subhalo spatial distribution, density profile, and projected number density and the dark matter fraction in subhaloes. We create mock lensing maps from the simulated haloes to study the differences in the lensing signal in the framework of subhalo detection. We find that subhalo convergence is well described by a log-normal distribution and that signal of subhaloes in the power spectrum is lower in SN models with respect to CDM, at a level of 10 to 80 per cent, depending on the scale. However, the scatter between different projections is large and might make the use of power-spectrum studies on the typical scales of current lensing images very difficult. Moreover, in the framework of individual detections through gravitational imaging a sample of ~30 lenses with an average sensitivity of M_sub=5x10^7M_sun would be required to discriminate between CDM and the considered sterile neutrino models.

We report a tentative detection of the circumgalactic medium (CGM) of WLM, an isolated, low-mass (log$M_*/M_\odot\approx7.6$), dwarf irregular galaxy in the Local Group (LG). We analyze an HST/COS archival spectrum of a quasar sightline (PHL2525), which is 45 kpc (0.5 virial radius) from WLM and close to the Magellanic Stream (MS). Along this sightline, two ion absorbers are detected in Si II, Si III, Si IV, C II, and C IV at velocities of $\sim$-220 km s$^{-1}$ (Component v-220) and $\sim$-150 km s$^{-1}$ (Component v-150). To identify their origins, we study the position-velocity alignment of the components with WLM and the nearby MS. Near the Magellanic longitude of PHL2525, the MS-related neutral and ionized gas moves at $\lesssim-190$ km s$^{-1}$, suggesting an MS origin for Component v-220, but not for Component v-150. Because PHL2525 passes near WLM and Component v-150 is close to WLM's systemic velocity ($\sim$-132 km s$^{-1}$), it is likely that Component v-150 arises from the galaxy's CGM. This results in a total Si mass in WLM's CGM of $M_{\rm Si}^{\rm CGM}\sim(0.2-1.0)\times10^5~M_\odot$ using assumption from other COS dwarf studies. Comparing $M_{\rm Si}^{\rm CGM}$ to the total Si mass synthesized in WLM over its lifetime ($\sim$1.3$\times10^5~M_\odot$), we find $\sim$3% is locked in stars, $\sim$6% in the ISM, $\sim$15%-77% in the CGM, and the rest ($\sim$14%-76%) is likely lost beyond the virial radius. Our finding resonates with other COS dwarf galaxy studies and theoretical predictions that low-mass galaxies can easily lose metals into their CGM due to stellar feedback and shallow gravitational potential.

We characterize the metallicities and physical properties of cool, photoionized gas in a sample of 152 z<1 strong Lya forest systems (SLFSs, absorbers with 15<log N(HI)<16.2). The sample is drawn from our COS circumgalactic medium (CGM) compendium (CCC), an ultraviolet survey of HI-selected circumgalactic gas around z<1 galaxies that targets 262 absorbers with 15<log N(HI)<19. We show that the metallicity probability distribution function of the SLFSs at z<1 is unimodal, skewed to low metallicities with a mean and median of [X/H]=-1.47$ and -1.18 dex. Very metal-poor gas with [X/H]<-1.4 represents about half of the population of absorbers with 15<log N(HI)<18. Thus, there are important reservoirs of primitive (though not pristine) gas around z<1 galaxies. The photoionized gas around z<1 galaxies is highly inhomogeneous based on the wide range of metallicities observed (-3<[X/H]<+0.4) and that there are large metallicity variations (factors of 2 to 25) for most of the closely-spaced absorbers (Dv<300 km/s) along the same sightlines. These absorbers show a complex evolution with redshift and HI column density, and we identify subtle cosmic evolution effects that affect the interpretation of metallicity distributions and comparison with other of absorbers samples. We discuss the physical conditions and cosmic baryon and metal budgets of the CCC absorbers. Finally, we compare the CCC results to recent cosmological zoom simulations and explore the origins of the 15<log N(HI)<19 absorbers within the EAGLE high-resolution simulations.

A millimeter-wave survey over half the sky, that spans frequencies in the range of 30 to 350 GHz, and that is both an order of magnitude deeper and of higher-resolution than currently funded surveys would yield an enormous gain in understanding of both fundamental physics and astrophysics. By providing such a deep, high-resolution millimeter-wave survey (about 0.5 uK-arcmin noise and 15 arcsecond resolution at 150 GHz), CMB-HD will enable major advances. It will allow 1) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10/hMpc), which probes dark matter particle properties. It will also allow 2) measurements of the thermal and kinetic Sunyaev-Zel'dovich effects on small scales to map the gas density and gas pressure profiles of halos over a wide field, which probes galaxy evolution and cluster astrophysics. In addition, CMB-HD would allow us to cross critical thresholds in fundamental physics: 3) ruling out or detecting any new, light (< 0.1eV), thermal particles, which could potentially be the dark matter, and 4) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe. Such a survey would also 5) monitor the transient sky by mapping the full observing region every few days, which opens a new window on gamma-ray bursts, novae, fast radio bursts, and variable active galactic nuclei. Moreover, CMB-HD would 6) provide a census of planets, dwarf planets, and asteroids in the outer Solar System, and 7) enable the detection of exo-Oort clouds around other solar systems, shedding light on planet formation. CMB-HD will deliver this survey in 5 years of observing half the sky, using two new 30-meter-class off-axis cross-Dragone telescopes to be located at Cerro Toco in the Atacama Desert. The telescopes will field about 2.4 million detectors (600,000 pixels) in total.

We use the EAGLE cosmological, hydrodynamical simulations to predict the column density and equivalent width distributions of intergalactic O VII ($E=574$ eV) and O VIII ($E=654$ eV) absorbers at low redshift. These two ions are predicted to account for 40% of the gas-phase oxygen, which implies that they are key tracers of cosmic metals. We find that their column density distributions evolve little at observable column densities from redshift 1 to 0, and that they are sensitive to AGN feedback, which strongly reduces the number of strong (column density $N \gtrsim 10^{16} \, \mathrm{cm}^{-2})$ absorbers. The distributions have a break at $N \sim 10^{16} \, \mathrm{cm}^{-2}$, corresponding to overdensities of $\sim 10^{2}$, likely caused by the transition from sheet/filament to halo gas. Absorption systems with $N \gtrsim 10^{16} \mathrm{cm}^{-2}$ are dominated by collisionally ionized O VII and O VIII, while the ionization state of oxygen at lower column densities is also influenced by photoionization. At these high column densities, O VII and O VIII arising in the same structures probe systematically different gas temperatures, meaning their line ratio does not translate into a simple estimate of temperature. While O VII and O VIII column densities and covering fractions correlate poorly with the H I column density at $N_{\mathrm{H \, I}} \gtrsim 10^{15} \, \mathrm{cm}^{-2}$, O VII and O VIII column densities are higher in this regime than at the more common, lower H I column densities. The column densities of O VI and especially Ne VIII, which have strong absorption lines in the UV, are good predictors of the strengths of O VII and O VIII absorption and can hence aid in the detection of the X-ray lines.

We examine the origin of scatter in the relationship between the gas fraction and mass of dark matter haloes hosting present-day $\sim L^\star$ central galaxies in the EAGLE simulations. The scatter is uncorrelated with the accretion rate of the central galaxy's black hole (BH), but correlates strongly and negatively with the BH's mass, implicating differences in the expulsion of gas by active galactic nucleus feedback, throughout the assembly of the halo, as the main cause of scatter. Haloes whose central galaxies host undermassive BHs also tend to retain a higher gas fraction, and exhibit elevated star formation rates (SFRs). Diversity in the mass of central BHs stems primarily from diversity in the dark matter halo binding energy, as these quantities are strongly and positively correlated at fixed halo mass, such that $\sim L^\star$ galaxies hosted by haloes that are more (less) tightly-bound develop central BHs that are more (less) massive than is typical for their halo mass. Variations in the halo gas fraction at fixed halo mass are reflected in both the soft X-ray luminosity and thermal Sunyaev-Zel'dovich flux, suggesting that the prediction of a strong coupling between the properties of galaxies and their halo gas fractions can be tested with measurements of these diagnostics for galaxies with diverse SFRs but similar halo masses.

Davies et al. (2019) established that for L^* galaxies the fraction of baryons in the circumgalactic medium (CGM) is inversely correlated with the mass of their central supermassive black holes (BHs) in the EAGLE hydrodynamic simulation. The interpretation is that, over time, a more massive BH has provided more energy to transport baryons beyond the virial radius, which additionally reduces gas accretion and star formation. We continue this research by focusing on the relationship between the 1) BH masses, 2) physical and observational properties of the CGM, and 3) galaxy colours for Milky Way-mass systems. The ratio of the cumulative BH feedback energy over the gaseous halo binding energy is a strong predictor of the CGM gas content, with BHs injecting >~10x the binding energy resulting in gas-poor haloes. Observable tracers of the CGM, including CIV, OVI, and HI absorption line measurements, are found to be effective tracers of the total z~0 CGM halo mass. We use high-cadence simulation outputs to demonstrate that BH feedback pushes baryons beyond the virial radius within 100 Myr timescales, but that CGM metal tracers take longer (0.5-2.5 Gyr) to respond. Secular evolution of galaxies results in blue, star-forming or red, passive populations depending on the cumulative feedback from BHs. The reddest quartile of galaxies with M_*=10^{10.2-10.7} M_solar (median u-r = 2.28) has a CGM mass that is 2.5x lower than the bluest quartile (u-r=1.59). We propose strategies for observing the predicted lower CGM column densities and covering fractions around galaxies hosting more massive BHs using the Cosmic Origins Spectrograph on Hubble.

A large fraction of the baryons and most of the metals in the Universe are unaccounted for. They likely lie in extended galaxy halos, galaxy groups, and the cosmic web, and measuring their nature is essential to understanding galaxy formation. These environments have virial temperatures >10^5.5 K, so the gas should be visible in X-rays. Here we show the breakthrough capabilities of grating spectrometers to 1) detect these reservoirs of hidden metals and mass, and 2) quantify hot gas flows, turbulence, and rotation around the Milky Way and external galaxies. Grating spectrometers are essential instruments for future X-ray missions, and existing technologies provide 50-1500-fold higher throughput compared to current orbiting instruments.

The majority of baryons reside beyond the optical extent of a galaxy in the circumgalactic and intergalactic media (CGM/IGM). Gaseous halos are inextricably linked to the appearance of their host galaxies through a complex story of accretion, feedback, and continual recycling. The energetic processes, which define the state of gas in the CGM, are the same ones that 1) regulate stellar growth so that it is not over-efficient, and 2) create the diversity of today's galaxy colors, SFRs, and morphologies spanning Hubble's Tuning Fork Diagram. They work in concert to set the speed of growth on the star-forming Main Sequence, transform a galaxy across the Green Valley, and maintain a galaxy's quenched appearance on the Red Sequence. Most baryons in halos more massive than 10^12 Msolar along with their high-energy physics and dynamics remain invisible because that gas is heated above the UV ionization states. We argue that information on many of the essential drivers of galaxy evolution is primarily contained in this "missing" hot gas phase. Completing the picture of galaxy formation requires uncovering the physical mechanisms behind stellar and SMBH feedback driving mass, metals, and energy into the CGM. By opening galactic hot halos to new wavebands, we not only obtain fossil imprints of >13 Gyrs of evolution, but observe on-going hot-mode accretion, the deposition of superwind outflows into the CGM, and the re-arrangement of baryons by SMBH feedback. A description of the flows of mass, metals, and energy will only be complete by observing the thermodynamic states, chemical compositions, structure, and dynamics of T>=10^6 K halos. These measurements are uniquely possible with a next-generation X-ray observatory if it provides the sensitivity to detect faint CGM emission, spectroscopic power to measure absorption lines and gas motions, and high spatial resolution to resolve structures.

The past decade has seen an explosion of discoveries and new insights into the diffuse gas within galaxies, galaxy clusters, and the filaments composing the Cosmic Web. A new decade will bring fresh opportunities to further this progress towards developing a comprehensive view of the composition, thermal state, and physical processes of diffuse gas in the Universe. Ultraviolet (UV) spectroscopy, probing diffuse 10^4-10^6 K gas at high spectral resolution, is uniquely poised to (1) witness environmental galaxy quenching processes in action, such as strangulation and tidal- and ram-pressure stripping, (2) directly account for the baryon content of galaxy clusters in the cold-warm (T<10^6 K) gas, (3) determine the phase structure and kinematics of gas participating in the equilibrium-regulating exchange of energy at the cores of galaxy clusters, and (4) map cold streams and filaments of the Cosmic Web that feed galaxies and clusters. With a substantial UV undertaking beyond the Hubble Space Telescope, all of the above would be achievable over the entire epoch of galaxy cluster formation. Such capabilities, coupled with already-planned advancements at other wavelengths, will transform extragalactic astronomy by revealing the dominant formation and growth mechanisms of gaseous halos over the mass spectrum, settling the debate between early- and late-time metal enrichment scenarios, and revealing how the ecosystems in which galaxies reside ultimately facilitate their demise.

The past decade began with the first light of ALMA and will end at the start of the new era of hyperdimensional astrophysics. Our community-wide movement toward highly multiwavelength and multidimensional datasets has enabled immense progress in each science frontier identified by the 2010 Decadal Survey, particularly with regard to black hole feedback and the cycle of baryons in galaxies. Facilities like ALMA and the next generation of integral field unit (IFU) spectrographs together enable mapping the physical conditions and kinematics of warm ionized and cold molecular gas in galaxies in unprecedented detail (Fig. 1). JWST's launch at the start of the coming decade will push this capability to the rest-frame UV at redshifts z > 6, mapping the birth of stars in the first galaxies at cosmic dawn. Understanding of their subsequent evolution, however, now awaits an ability to map the processes that transform galaxies directly, rather than the consequences of those processes in isolation. In this paper, we argue that doing so requires an equivalent revolution in spatially resolved spectroscopy for the hot plasma that pervades galaxies, the atmospheres in which they reside, and the winds that are the engines of their evolution.

A small survey of the UV-absorbing gas in 12 low-$z$ galaxy groups has been conducted using the Cosmic Origins Spectrograph (COS) on-board the Hubble Space Telescope (HST). Targets were selected from a large, homogeneously-selected sample of groups found in the Sloan Digital Sky Survey (SDSS). A critical selection criterion excluded sight lines that pass close ($<1.5$ virial radii) to a group galaxy, to ensure absorber association with the group as a whole. Deeper galaxy redshift observations are used both to search for closer galaxies and also to characterize these $10^{13.5}$ to $10^{14.5} M_{\odot}$ groups, the most massive of which are highly-virialized with numerous early-type galaxies (ETGs). This sample also includes two spiral-rich groups, not yet fully-virialized. At group-centric impact parameters of 0.3-2 Mpc, these $\mathrm{S/N}=15$-30 spectra detected HI absorption in 7 of 12 groups; high (OVI) and low (SiIII) ion metal lines are present in 2/3 of the absorption components. None of the three most highly-virialized, ETG-dominated groups are detected in absorption. Covering fractions $\gtrsim50$% are seen at all impact parameters probed, but do not require large filling factors despite an enormous extent. Unlike halo clouds in individual galaxies, group absorbers have radial velocities which are too low to escape the group potential well without doubt. This suggests that these groups are "closed boxes" for galactic evolution in the current epoch. Evidence is presented that the cool and warm group absorbers are not a pervasive intra-group medium (IGrM), requiring a hotter ($T\sim10^6$ to $10^7$ K) IGrM to be present to close the baryon accounting.

We present the results from our COS circumgalactic medium (CGM) compendium (CCC), a survey of the CGM at z<1 using HI-selected absorbers with 15<log N(HI) <19. We focus here on 82 partial Lyman limit systems (pLLSs, 16.2<log N(HI) <17.2) and 29 LLSs (17.2<log N(HI) <19). Using Bayesian techniques and Markov-chain Monte Carlo sampling of a grid of photoionization models, we derive the posterior probability distribution functions (PDFs) for the metallicity of each absorber in CCC. We show that the combined pLLS metallicity PDF at z<1 has two main peaks at [X/H]=-1.7 and -0.4, with a strong dip at [X/H]=-1. The metallicity PDF of the LLSs might be more complicated than an unimodal or bimodal distribution. The pLLSs and LLSs probe a similar range of metallicities -3<[X/H]<+0.4, but the fraction of very metal-poor absorbers with [X/H]<-1.4 is much larger for the pLLSs than the LLSs. In contrast, absorbers with log N(HI)>19 have mostly -1<[X/H]<0 at z<1. The metal-enriched gas probed by pLLSs and LLSs confirms that galaxies that have been enriching their CGM over billions of years. Surprisingly, despite this enrichment, there is also abundant metal-poor CGM gas (41-59% of the pLLSs have [X/H]<-1.4), representing a reservoir of near-pristine gas around z<1 galaxies. We compare our empirical results to recent cosmological zoom simulations, finding some discrepancies, including an overabundance of metal-enriched CGM gas in simulations.

We present a neutral hydrogen-selected absorption-line survey of gas with HI column densities 15<log N(HI)<19 at z<1 using the Cosmic Origins Spectrograph on the Hubble Space Telescope. Our main aim is to determine the metallicity distribution of these absorbers. Our sample consists of 224 absorbers selected on the basis of their HI absorption strength. Here we discuss the properties of our survey and the immediate empirical results. We find singly and doubly ionized metal species and HI typically have similar velocity profiles, implying they probe gas in the same or similar environments. The column density ionic ratios (e.g., CII/CIII, OI/CII) indicate the gas in these absorbers is largely ionized, and the ionization conditions are quite comparable across the sampled N(HI) range. The Doppler parameters of the HI imply T<50,000 K on average, consistent with the gas being photoionized. The MgII column densities span >2 orders of magnitude at any given N(HI), indicating a wide range of metallicities (from solar to <1/100 solar). In the range 16.2<log N(HI)<17, there is a gap in the N(MgII) distribution corresponding to gas with ~10% solar metallicity, consistent with the gap seen in the previously identified bimodal metallicity distribution in this column density regime. Less than 3% of the absorbers in our sample show no detectable metal absorption, implying truly-pristine gas at z<1 is uncommon. We find [FeII/MgII] = -0.4+/-0.3, and since alpha-enhancement can affect this ratio, dust depletion is extremely mild.

Hydrostatic equilibrium (HSE), where the thermal pressure gradient balances the force of gravity, is tested across a range of simulated EAGLE haloes from Milky Way L* haloes (M_200~10^12 Msol) to cluster scales. Clusters (M_200>=10^14 Msol) reproduce previous results with thermal pressure responsible for ~90% of the support against gravity, but this fraction drops for group-sized haloes (M_200~10^13 Msol) and is even lower (40-70%) for L* haloes between 0.1-0.3 R_200. Energy from feedback grows relative to the binding energy of a halo toward lower mass resulting in greater deviations from HSE. Tangential motions comprise the largest deviation from HSE in L* haloes indicating that the hot circumgalactic medium (CGM) has significant sub- centrifugal rotation and angular momentum spin parameters 2-3x higher than the dark matter spin parameters. Thermal feedback can buoyantly rise to the outer CGM of M_200<10^12 Msol haloes, both moving baryons beyond R_200 and feeding uncorrelated tangential motions. The resulting hot halo density and rotation profiles show promising agreement with X-ray observations of the inner Milky Way halo, and we discuss future observational prospects to detect spinning hot haloes around other galaxies. Acceleration and radial streaming motions also comprise significant deviations from HSE, especially net outward accelerations seen in L* and group haloes indicating active feedback. Black hole feedback acts in a preventative manner during the later growth of group haloes, applying significant accelerations via shocks that do not feed tangential motions. We argue that HSE is a poor assumption for the CGM, especially in the inner regions, and rotating baryonic hot haloes are a critical consideration for analytic models of the CGM.

We present the Technicolor Dawn simulations, a suite of cosmological radiation-hydrodynamic simulations of the first 1.2 billion years. By modeling a spatially-inhomogeneous UVB on-the-fly with 24 frequencies and resolving dark matter halos down to $10^8 M_\odot$ within 12 $h^{-1}$ Mpc volumes, our simulations unify observations of the intergalactic and circumgalactic media, galaxies, and reionization into a common framework. The only empirically-tuned parameter, the fraction $f_{\mathrm{esc,gal}}(z)$ of ionizing photons that escape the interstellar medium, is adjusted to match observations of the Lyman-$\alpha$ forest and the cosmic microwave background. With this single calibration, our simulations reproduce the history of reionization; the stellar mass-star formation rate relation of galaxies; the number density and metallicity of damped Lyman-$\alpha$ absorbers (DLAs) at $z\sim5$; the abundance of weak metal absorbers; the ultraviolet background (UVB) amplitude; and the Lyman-$\alpha$ flux power spectrum at $z=5.4$. The galaxy stellar mass and UV luminosity functions are underproduced by $\leq2\times$, suggesting an overly vigorous feedback model. The mean transmission in the Lyman-$\alpha$ forest is underproduced at $z<6$, indicating tension between measurements of the UVB amplitude and Lyman-$\alpha$ transmission. The observed SiIV column density distribution is reasonably well-reproduced ($\sim 1\sigma$ low). By contrast, CIV remains significantly underproduced despite being boosted by an intense $>4$ Ryd UVB. Solving this problem by increasing metal yields would overproduce both weak absorbers and DLA metallicities. Instead, the observed strength of high-ionization emission from high-redshift galaxies and absorption from their environments suggest that the ionizing flux from conventional stellar population models is too soft.

Active galactic nuclei (AGN) are thought to play a critical role in shaping galaxies, but their effect on the circumgalactic medium (CGM) is not well studied. We present results from the COS-AGN survey: 19 quasar sightlines that probe the CGM of 20 optically-selected AGN host galaxies with impact parameters $80 < \rho_{imp} < 300$ kpc. Absorption lines from a variety of species are measured and compared to a stellar mass and impact parameter matched sample of sightlines through non-AGN galaxies. Amongst the observed species in the COS-AGN sample (HI, CII, SiII, SiIII, CIV, SiIV, NV), only Ly$\alpha$ shows a high covering fraction ($94^{+6}_{-23}$% for rest-frame equivalent widths EW $> 124$ m\AA) whilst many of the metal ions are not detected in individual sightlines. A sightline-by-sightline comparison between COS-AGN and the control sample yields no significant difference in EW distribution. However, stacked spectra of the COS-AGN and control samples show significant (> 3 sigma) enhancements in the EW of both Ly$\alpha$ and SiIII at impact parameters $> 164$ kpc by a factor of $+0.45\pm0.05$ dex and $> +0.75$ dex respectively. The lack of detections of both high-ionization species near the AGN and strong kinematic offsets between the absorption systemic galaxy redshifts indicates that neither the AGN's ionization nor its outflows are the origin of these differences. Instead, we suggest the observed differences could result from either AGN hosts residing in haloes with intrinsically distinct gas properties, or that their CGM has been affected by a previous event, such as a starburst, which may also have fuelled the nuclear activity.

The metallicity of strong HI systems, spanning from damped Lyman-alpha absorbers (DLAs) to Lyman-limit systems (LLSs) is explored between z = 5->0 using the EAGLE high-resolution cosmological hydrodynamic simulation of galaxy formation. The metallicities of LLSs and DLAs steadily increase with time in agreement with observations. DLAs are more metal rich than LLSs, although the metallicities in the LLS column density range (NHI = 10^17 -10^20 cm^-2) are relatively flat, evolving from a median HI-weighted metallicity of Z<10^-2 Zsol at z = 3 to ~10^-0.5 Zsol by z = 0. The metal content of HI systems tracks the increasing stellar content of the Universe, holding ~5% of the integrated total metals released from stars at z = 0. We also consider partial LLS (pLLS, NHI = 10^16-10^17 cm^-2) metallicities, and find good agreement with Wotta et al. (2016) for the fraction of systems above (37%) and below (63%) 0.1 Zsol. We also find a large dispersion of pLLS metallicities, although we do not reproduce the observed metallicity bimodality and instead we make the prediction that a larger sample will yield more pLLSs around 0.1 Zsol. We under-predict the median metallicity of strong LLSs, and predict a population of Z < 10^-3 Zsol DLAs at z > 3 that are not observed, which may indicate more widespread early enrichment in the real Universe compared to EAGLE.

Proximity zone fossils (PZFs) are ionization signatures around recently active galactic nuclei (AGN) where metal species in the circumgalactic medium remain over-ionized after the AGN has shut-off due to their long recombination timescales. We explore cosmological zoom hydrodynamic simulations using the EAGLE model paired with a non-equilibrium ionization and cooling module including time-variable AGN radiation to model PZFs around star-forming, disk galaxies in the z~0.2 Universe. Previous simulations typically under-estimated the O VI content of galactic haloes, but we show that plausible PZF models increase O VI column densities by 2-3x to achieve the levels observed around COS-Halos star-forming galaxies out to 150 kpc. Models with AGN bolometric luminosities >~10^43.6 erg s^-1, duty cycle fractions <~10%, and AGN lifetimes <~10^6 yr are the most promising, because their super-massive black holes grow at the cosmologically expected rate and they mostly appear as inactive AGN, consistent with COS-Halos. The central requirement is that the typical star-forming galaxy hosted an active AGN within a timescale comparable to the recombination time of a high metal ion, which for circumgalactic O VI is 10^7 years. H I, by contrast, returns to equilibrium much more rapidly due to its low neutral fraction and does not show a significant PZF effect. O VI absorption features originating from PZFs appear narrow, indicating photo-ionization, and are often well-aligned with lower metal ion species. PZFs are highly likely to affect the physical interpretation of circumgalactic high ionization metal lines if, as expected, normal galaxies host flickering AGN.

We use observations of spatially-aligned C II, C IV, Si II, Si IV, and O I absorbers to probe the slope and intensity of the ultraviolet background (UVB) at $z \sim 6$. We accomplish this by comparing observations with predictions from a cosmological hydrodynamic simulation using three trial UVBs applied in post-processing: a spectrally soft, fluctuating UVB calculated using multi-frequency radiative transfer; a soft, spatially-uniform UVB; and a hard, spatially-uniform "quasars-only" model. When considering our paired high-ionization absorbers (C IV/Si IV), the observed statistics strongly prefer the hard, spatially-uniform UVB. This echoes recent findings that cosmological simulations generically underproduce strong C IV absorbers at $z>5$. A single low/high ionization pair (Si II/Si IV), by contrast, shows a preference for the HM12 UVB, while two more (C II/C IV and O I/C IV) show no preference for any of the three UVBs. Despite this, future observations of specific absorbers, particularly Si IV/C IV, with next-generation telescopes probing to lower column densities should yield tighter constraints on the UVB.

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$.

We explore the circumgalactic metal content traced by commonly observed low ion absorbers, including C II, Si II, Si III, Si IV, and Mg II. We use a set of cosmological hydrodynamical zoom simulations run with the EAGLE model and including a non-equilibrium ionization and cooling module that follows 136 ions. The simulations of z~0.2 L* (M_200=10^11.7-10^12.3 Msol) haloes hosting star-forming galaxies and group-sized (M_200=10^12.7-10^13.3 Msol) haloes hosting mainly passive galaxies reproduce key trends observed by the COS-Halos survey-- low ion column densities show 1) little dependence on galaxy specific star formation rate, 2) a patchy covering fraction indicative of 10^4 K clouds with a small volume filling factor, and 3) a declining covering fraction as impact parameter increases from 20-160 kpc. Simulated Si II, Si III, Si IV, C II, and C III column densities show good agreement with observations, while Mg II is under-predicted. Low ions trace a significant metal reservoir, ~10^8 Msol, residing primarily at 10-100 kpc from star-forming and passive central galaxies. These clouds tend to flow inwards and most will accrete onto the central galaxy within the next several Gyr, while a small fraction are entrained in strong outflows. A two-phase structure describes the inner CGM (<0.5 R_200) with low-ion metal clouds surrounded by a hot, ambient medium. This cool phase is separate from the O VI observed by COS-Halos, which arises from the outer CGM (>0.5 R_200) tracing virial temperature gas around L* galaxies. Physical parameters derived from standard photo-ionization modelling of observed column densities (e.g. aligned Si II/Si III absorbers) are validated against our simulations. Our simulations therefore support previous ionization models indicating that cloud covering factors decline while densities and pressures show little variation with increasing impact parameter.

We study the effect of a fluctuating active galactic nucleus (AGN) on the abundance of circumgalactic OVI in galaxies selected from the EAGLE simulations. We follow the time-variable OVI abundance in post-processing around four galaxies - two at $z=0.1$ with stellar masses of $M_{\ast} \sim 10^{10}$ M$_{\odot}$ and $M_{\ast} \sim 10^{11}$ M$_{\odot}$, and two at $z=3$ with similar stellar masses - out to impact parameters of twice their virial radii, implementing a fluctuating central source of ionizing radiation. Due to delayed recombination, the AGN leave significant `AGN proximity zone fossils' around all four galaxies, where OVI and other metal ions are out of ionization equilibrium for several megayears after the AGN fade. The column density of OVI is typically enhanced by $\approx 0.3-1.0$ dex at impact parameters within $0.3R_{\rm vir}$, and by $\approx 0.06-0.2$ dex at $2R_{\rm vir}$, thereby also enhancing the covering fraction of OVI above a given column density threshold. The fossil effect tends to increase with increasing AGN luminosity, and towards shorter AGN lifetimes and larger AGN duty cycle fractions. In the limit of short AGN lifetimes, the effect converges to that of a continuous AGN with a luminosity of $(f_{\rm duty}/100\%)$ times the AGN luminosity. We also find significant fossil effects for other metal ions, where low-ionization state ions are decreased (SiIV, CIV at $z=3$) and high-ionization state ions are increased (CIV at $z=0.1$, NeVIII, MgX). Using observationally motivated AGN parameters, we predict AGN proximity zone fossils to be ubiquitous around $M_{\ast} \sim 10^{10-11}$ M$_{\odot}$ galaxies, and to affect observations of metals in the circumgalactic medium at both low and high redshifts.

In preparation for a Hubble Space Telescope (HST) observing project using the Cosmic Origins Spectrograph (COS), the positions of all AGN targets having high-S/N far-UV G130M spectra were cross-correlated with a large catalog of low-redshift galaxy groups homogeneously selected from the spectroscopic sample of the Sloan Digital Sky Survey (SDSS). Searching for targets behind only those groups at z = 0.1-0.2 (which places the OVI doublet in the wavelength region of peak COS sensitivity) we identified only one potential S/N = 15-20 target, FBQS 1010+3003. An OVI-only absorber was found in its G130M spectrum at z = 0.11326, close to the redshift of a foreground small group of luminous galaxies at z = 0.11685. Because there is no associated Lyalpha absorption, any characterization of this absorber is necessarily minimal; however, the OVI detection likely traces "warm" gas in collisional ionization equilibrium at T ~ 300,000 K. While this discovery is consistent with being interface gas between cooler, photoionized clouds and a hotter intra-group medium, it could also be warm, interface gas associated with the circum-galactic medium (CGM) of the single closest galaxy. In this case a detailed analysis of the galaxy distribution (complete to 0.2 L*) strongly favors the individual galaxy association. This analysis highlights the necessity of both high-S/N > 20 COS data and a deep galaxy redshift survey of the region in order to test more rigorously the association of OVI-absorbing gas with a galaxy group. A Cycle 23 HST/COS program currently is targeting 10 UV-bright AGN behind 12 low-redshift galaxy groups to test the warm, group gas hypothesis.

The total contribution of diffuse halo gas to the galaxy baryon budget strongly depends on its dominant ionization state. In this paper, we address the physical conditions in the highly ionized circumgalactic medium (CGM) traced by OVI absorption lines observed in COS-Halos spectra. We analyze the observed ionic column densities, absorption-line widths and relative velocities, along with the ratios of NV/OVI for 39 fitted Voigt profile components of OVI. We compare these quantities with the predictions given by a wide range of ionization models. Photoionization models that include only extragalactic UV background radiation are ruled out; conservatively, the upper limits to NV/OVI and measurements of N$_{\rm OVI}$ imply unphysically large path lengths > 100 kpc. Furthermore, very broad OVI absorption (b > 40 km s$^{-1}$) is a defining characteristic of the CGM of star-forming L* galaxies. We highlight two possible origins for the bulk of the observed OVI: (1) highly structured gas clouds photoionized primarily by local high energy sources or (2) gas radiatively cooling on large scales behind a supersonic wind. Approximately 20% of circumgalactic OVI does not align with any low-ionization state gas within $\pm$50 km s$^{-1}$ and is found only in halos with M$_{\rm halo}$ < 10$^{12}$ M$_{\odot}$. We suggest that this type of unmatched OVI absorption traces the hot corona itself at a characteristic temperature of 10$^{5.5}$ K. We discuss the implications of these very distinct physical origins for the dynamical state, gas cooling rates, and total baryonic content of L* gaseous halos.

Recent analysis of strongly-lensed sources in the Hubble Frontier Fields indicates that the rest-frame UV luminosity function of galaxies at $z=$6--8 rises as a power law down to $M_\mathrm{UV}=-15$, and possibly as faint as -12.5. We use predictions from a cosmological radiation hydrodynamic simulation to map these luminosities onto physical space, constraining the minimum dark matter halo mass and stellar mass that the Frontier Fields probe. While previously-published theoretical studies have suggested or assumed that early star formation was suppressed in halos less massive than $10^9$--$10^{11} M_\odot$, we find that recent observations demand vigorous star formation in halos at least as massive as (3.1, 5.6, 10.5)$\times10^9 M_\odot$ at $z=(6,7,8)$. Likewise, we find that Frontier Fields observations probe down to stellar masses of (8.1, 18, 32)$\times10^6 M_\odot$; that is, they are observing the likely progenitors of analogues to Local Group dwarfs such as Pegasus and M32. Our simulations yield somewhat different constraints than two complementary models that have been invoked in similar analyses, emphasizing the need for further observational constraints on the galaxy-halo connection.

We introduce a series of 20 cosmological hydrodynamical simulations of Lstar (M_200 =10^11.7 - 10^12.3 Msol) and group-sized (M_200 = 10^12.7 - 10^13.3 Msol) haloes run with the model used for the EAGLE project, which additionally includes a non-equilibrium ionization and cooling module that follows 136 ions. The simulations reproduce the observed correlation, revealed by COS-Halos at z~0.2, between O VI column density at impact parameters b < 150 kpc and the specific star formation rate (sSFR=SFR/Mstar) of the central galaxy at z~0.2. We find that the column density of circumgalactic O VI is maximal in the haloes associated with Lstar galaxies, because their virial temperatures are close to the temperature at which the ionization fraction of O VI peaks (T~10^5.5 K). The higher virial temperature of group haloes (> 10^6 K) promotes oxygen to higher ionization states, suppressing the O VI column density. The observed NO VI-sSFR correlation therefore does not imply a causal link, but reflects the changing characteristic ionization state of oxygen as halo mass is increased. In spite of the mass-dependence of the oxygen ionization state, the most abundant circumgalactic oxygen ion in both Lstar and group haloes is O VII; O VI accounts for only 0.1% of the oxygen in group haloes and 0.9-1.3% with Lstar haloes. Nonetheless, the metals traced by O VI absorbers represent a fossil record of the feedback history of galaxies over a Hubble time; their characteristic epoch of ejection corresponds to z > 1 and much of the ejected metal mass resides beyond the virial radius of galaxies. For both Lstar and group galaxies, more of the oxygen produced and released by stars resides in the circumgalactic medium (within twice the virial radius) than in the stars and ISM of the galaxy.

The sources that drove cosmological reionization left clues regarding their identity in the slope and inhomogeneity of the ultraviolet ionizing background (UVB): Bright quasars (QSOs) generate a hard UVB with predominantly large-scale fluctuations while Population II stars generate a softer one with smaller-scale fluctuations. Metal absorbers probe the UVB's slope because different ions are sensitive to different energies. Likewise, they probe spatial fluctuations because they originate in regions where a galaxy-driven UVB is harder and more intense. We take a first step towards studying the reionization-epoch UVB's slope and inhomogeneity by comparing observations of 12 metal absorbers at $z\sim6$ versus predictions from a cosmological hydrodynamic simulation using three different UVBs: a soft, spatially-inhomogeneous "galaxies+QSOs" UVB; a homogeneous "galaxies+QSOs" UVB (Haardt & Madau 2012); and a QSOs-only model. All UVBs reproduce the observed column density distributions of CII, SiIV, and CIV reasonably well although high-column, high-ionization absorbers are underproduced, reflecting numerical limitations. With upper limits treated as detections, only a soft, fluctuating UVB reproduces both the observed SiIV/CIV and CII/CIV distributions. The QSOs-only UVB overpredicts both CIV/CII and CIV/SiIV, indicating that it is too hard. The Haardt & Madau (2012) UVB underpredicts CIV/SiIV, suggesting that it lacks amplifications near galaxies. Hence current observations prefer a soft, fluctuating UVB as expected from a predominantly Population II background although they cannot rule out a harder one. Future observations probing a factor of two deeper in metal column density will distinguish between the soft, fluctuating and QSOs-only UVBs.

We study the distribution and evolution of highly ionised intergalactic metals in the Evolution and Assembly of Galaxies and their Environment (EAGLE) cosmological, hydrodynamical simulations. EAGLE has been shown to reproduce a wide range of galaxy properties while its subgrid feedback was calibrated without considering gas properties. We compare the predictions for the column density distribution functions (CDDFs) and cosmic densities of SiIV, CIV, NV, OVI and NeVIII absorbers with observations at redshift z = 0 to ~ 6 and find reasonable agreement, although there are some differences. We show that the typical physical densities of the absorbing gas increase with column density and redshift, but decrease with the ionization energy of the absorbing ion. The typical metallicity increases with both column density and time. The fraction of collisionally ionized metal absorbers increases with time and ionization energy. While our results show little sensitivity to the presence or absence of AGN feedback, increasing/decreasing the efficiency of stellar feedback by a factor of two substantially decreases/increases the CDDFs and the cosmic densities of the metal ions. We show that the impact of the efficiency of stellar feedback on the CDDFs and cosmic densities is largely due to its effect on the metal production rate. However, the temperatures of the metal absorbers, particularly those of strong OVI, are directly sensitive to the strength of the feedback.

The thermal and chemical properties of the hot diffuse intragroup medium (IGrM) provide important constraints on the feedback processes associated with massive galaxy formation and evolution. Here we explore these constraints via a detailed analysis of the global properties of simulated z<3 galaxy groups from a cosmological simulation that includes a well-constrained prescription for stellar/supernovae-powered galactic outflows but no AGN feedback. Our aims are to (a) establish a baseline against which we will compare future models; (b) identify model successes due to stellar/supernovae-powered outflows; and (c) pinpoint features that signal the need for, and constrain the nature of, AGN feedback. Our simulation successfully reproduces key observed z=0 group IGrM properties, including the various X-ray Lx - Tx - entropy scaling relations, for all but the most massive groups. The z<1 redshift evolution of these also agree with observations. Contrary to expectations, the simulated groups' IGrM does not suffer catastrophic cooling. Yet, the z=0 group stellar mass is ~ 2X too large. This is due to the build-up of cold gas in the massive galaxies before they are incorporated inside groups. This in turn indicates that other feedback mechanisms must activate in real galaxies once their stellar masses grow to a few X 10^{10} M_sun. We show that these must be powerful enough to expel a significant fraction of the gas from the galactic halos. Gentle maintenance-mode (quenching) AGN feedback, as seen in galaxy clusters, will not do. Just as importantly, we find that the stellar/supernovae-powered winds are essential for understanding the IGrM metal abundances. Our simulation is able to reproduce the observed relationship between the global IGrM iron and silicon abundance and the group X-ray temperature, and these results ought to be relatively insensitive to the addition of AGN feedback.

We quantitatively examine the effects of accretion and environment on the HI content of galaxies within a cosmological hydrodynamic simulation that reproduces basic observed trends of HI in galaxies. We show that our model broadly reproduces the observed scatter in HI at a given stellar mass as quantified by the HI mass function in bins of stellar mass, as well as the HI richness versus local galaxy density. This shows that the predicted HI fluctuations and environmental effects are roughly consistent with data with few minor discrepancies. For satellite galaxies in >= 10^12M_* halos, the HI richness distribution is bimodal and drops towards the largest halo masses. The depletion rate of HI once a galaxy enters a more massive halo is more rapid at higher halo mass, in contrast to the specific star formation rate which shows much less variation in the attenuation rate versus halo mass. This suggests that, up to halo mass scales probed here (<= 10^14M_*), star formation is mainly attenuated by starvation, but HI is additionally removed by stripping once a hot gaseous halo is present. In low mass halos, the HI richness of satellites is independent of radius, while in high mass halos they become gas-poor towards the center. We show that the gas fraction of satellite and central galaxies decreases from z=5 to z=0,tracking each other until z~1 after which the satellites' HI content drops much more quickly, particularly for the highest halo masses. Mergers somewhat increase the HI richness and its scatter about the mean relation, tracking the metallicity in a way consistent with arising from inflow fluctuations, while star formation is significantly boosted relative to HI.

Simulations of the formation of large-scale structure predict that dark matter, low density highly ionized gas, and galaxies form 10 40 Mpc scale filaments. These structure are easily recognized in the distribution of galaxies, but have not been directly observed in the distribution of the gas. We use Ly-alpha absorption lines in the spectra of 24 AGN to present a new way to probe these filaments. We use a new catalogue of nearby (cz<10,000 km/s) galaxies, complete down to a luminosity of about 0.05 L* for the region of space analyzed here. Using HST spectra of 24 AGN we sample the gas associated with a 30x5 Mpc galaxy filament at cz~3500 km/s. All of our sightlines pass outside the virial radius of any known filament galaxy. Within 500 kpc of the filament axis the detection rate is ~80%, while no detections are seen more than 2.1 Mpc from the filament. The width of the Lya lines correlates with filament impact parameter and the four BLAs in our sample all occur within 400 kpc of the filament axis, indicating increased temperature and/or turbulence. Comparing to simulations, we find that the recent Haardt & Madau (2012) extragalactic ionizing background predicts a factor 3-5 too few ionizing photons. Using a more intense radiation field matches the hydrogen density profile within 2.1 Mpc of the filament axis, but the simulations still overpredict the detection rate between 2.1 and 5 Mpc from the axis.

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.

We combine cosmological hydrodynamic simulations with analytic models to evaluate the role of galaxy-scale gravitational torques on the evolution of massive black holes at the centers of star-forming galaxies. We confirm and extend our earlier results to show that torque-limited growth yields black holes and host galaxies evolving on average along the Mbh-Mbulge relation from early times down to z = 0 and that convergence onto the scaling relation occurs independent of the initial conditions and with no need for mass averaging through mergers or additional self-regulation processes. Smooth accretion dominates the long-term evolution, with black hole mergers with mass ratios >1:5 representing typically a small fraction of the total growth. Winds from the accretion disk are required to eject significant mass to suppress black hole growth, but there is no need for coupling this wind to galactic-scale gas to regulate black holes in a non-linear feedback loop. Torque-limited growth yields a close-to-linear relation for the star formation rate and the black hole accretion rate averaged over galaxy evolution time scales. However, the SFR-AGN connection has significant scatter owing to strong variability of black hole accretion at all resolved time scales. Eddington ratios can be described by a broad lognormal distribution with median value evolving roughly as (1 + z)^1.9, suggesting a main sequence for black hole growth similar to the cosmic evolution of specific SFRs. Our results offer an attractive scenario consistent with available observations in which cosmological gas infall and transport of angular momentum in the galaxy by gravitational instabilities regulate the long-term co-evolution of black holes and star-forming galaxies.

Observations suggest that CII was more abundant than CIV in the intergalactic medium towards the end of the hydrogen reionization epoch. This transition provides a unique opportunity to study the enrichment history of intergalactic gas and the growth of the ionizing background (UVB) at early times. We study how carbon absorption evolves from z=10-5 using a cosmological hydrodynamic simulation that includes a self-consistent multifrequency UVB as well as a well-constrained model for galactic outflows to disperse metals. Our predicted UVB is within 2-4 times that of Haardt & Madau (2012), which is fair agreement given the uncertainties. Nonetheless, we use a calibration in post-processing to account for Lyman-alpha forest measurements while preserving the predicted spectral slope and inhomogeneity. The UVB fluctuates spatially in such a way that it always exceeds the volume average in regions where metals are found. This implies both that a spatially-uniform UVB is a poor approximation and that metal absorption is not sensitive to the epoch when HII regions overlap globally even at column densites of 10^{12} cm^{-2}. We find, consistent with observations, that the CII mass fraction drops to low redshift while CIV rises owing the combined effects of a growing UVB and continued addition of carbon in low-density regions. This is mimicked in absorption statistics, which broadly agree with observations at z=6-3 while predicting that the absorber column density distributions rise steeply to the lowest observable columns. Our model reproduces the large observed scatter in the number of low-ionization absorbers per sightline, implying that the scatter does not indicate a partially-neutral Universe at z=6.

We report new observations of circumgalactic gas from the COS-Dwarfs survey, a systematic investigation of the gaseous halos around 43 low-mass z $\leq$ 0.1 galaxies using background QSOs observed with the Cosmic Origins Spectrograph. From the projected 1D and 2D distribution of C IV absorption, we find that C IV absorption is detected out to ~ 0.5 R$_{vir}$ of the host galaxies. The C IV absorption strength falls off radially as a power law and beyond 0.5 R$_{vir}$, no C IV absorption is detected above our sensitivity limit of ~ 50-100 m$\AA$. We find a tentative correlation between detected C IV absorption strength and star formation, paralleling the strong correlation seen in highly ionized oxygen for L~L* galaxies by the COS-Halos survey. The data imply a large carbon reservoir in the CGM of these galaxies, corresponding to a minimum carbon mass of $\gtrsim$ 1.2$\times 10^6$ $M_\odot$ out to ~ 110 kpc. This mass is comparable to the carbon mass in the ISM and more than the carbon mass currently in stars of these galaxies. The C IV absorption seen around these sub-L* galaxies can account for almost two-thirds of all $W_r$> 100 m$\AA$ C IV absorption detected at low z. Comparing the C IV covering fraction with hydrodynamical simulations, we find that an energy-driven wind model is consistent with the observations whereas a wind model of constant velocity fails to reproduce the CGM or the galaxy properties.

We analyze the physical conditions of the cool, photoionized (T $\sim 10^4$ K) circumgalactic medium (CGM) using the COS-Halos suite of gas column density measurements for 44 gaseous halos within 160 kpc of $L \sim L^*$ galaxies at $z \sim 0.2$. These data are well described by simple photoionization models, with the gas highly ionized (n$_{\rm HII}$/n$_{\rm H} \gtrsim 99\%$) by the extragalactic ultraviolet background (EUVB). Scaling by estimates for the virial radius, R$_{\rm vir}$, we show that the ionization state (tracked by the dimensionless ionization parameter, U) increases with distance from the host galaxy. The ionization parameters imply a decreasing volume density profile n$_{\rm H}$ = (10$^{-4.2 \pm 0.25}$)(R/R$_{\rm vir})^{-0.8\pm0.3}$. Our derived gas volume densities are several orders of magnitude lower than predictions from standard two-phase models with a cool medium in pressure equilibrium with a hot, coronal medium expected in virialized halos at this mass scale. Applying the ionization corrections to the HI column densities, we estimate a lower limit to the cool gas mass M$_{\rm CGM}^{\rm cool} > 6.5 \times 10^{10}$ M$_{\odot}$ for the volume within R $<$ R$_{\rm vir}$. Allowing for an additional warm-hot, OVI-traced phase, the CGM accounts for at least half of the baryons purported to be missing from dark matter halos at the 10$^{12}$ M$_{\odot}$ scale.

Using the recently installed Gemini Planet Imager (GPI), we have taken the first H-band spectrum of the planetary companion to the nearby young star beta Pictoris. GPI is designed to image and provide low-resolution spectra of Jupiter sized, self-luminous planetary companions around young nearby stars. These observations were taken covering the H-band (1.65 microns). The spectrum has a resolving power of $\sim$ 45 and demonstrates the distinctive triangular shape of a cool substellar object with low surface gravity. Using atmospheric models, we find an effective temperature of $1650 \pm 50$ K and a surface gravity of $\log(g) = 4.0 \pm 0.25$ (cgs units). These values agree well with predictions from planetary evolution models for a gas giant with mass between 10 and 12 $M_{\rm Jup}$ and age between 10 and 20 Myrs.

We present the first results from the polarimetry mode of the Gemini Planet Imager (GPI), which uses a new integral field polarimetry architecture to provide high contrast linear polarimetry with minimal systematic biases between the orthogonal polarizations. We describe the design, data reduction methods, and performance of polarimetry with GPI. Point spread function subtraction via differential polarimetry suppresses unpolarized starlight by a factor of over 100, and provides sensitivity to circumstellar dust reaching the photon noise limit for these observations. In the case of the circumstellar disk around HR 4796A, GPI's advanced adaptive optics system reveals the disk clearly even prior to PSF subtraction. In polarized light, the disk is seen all the way in to its semi-minor axis for the first time. The disk exhibits surprisingly strong asymmetry in polarized intensity, with the west side >9 times brighter than the east side despite the fact that the east side is slightly brighter in total intensity. Based on a synthesis of the total and polarized intensities, we now believe that the west side is closer to us, contrary to most prior interpretations. Forward scattering by relatively large silicate dust particles leads to the strong polarized intensity on the west side, and the ring must be slightly optically thick in order to explain the lower brightness in total intensity there. These findings suggest that the ring is geometrically narrow and dynamically cold, perhaps shepherded by larger bodies in the same manner as Saturn's F ring.

The Gemini Planet Imager (GPI) combines extreme adaptive optics, an integral field spectrograph, and a high performance coronagraph to directly image extrasolar planets in the near-infrared. Because the coronagraph blocks most of the light from the star, it prevents the properties of the host star from being measured directly. Instead, satellite spots, which are created by diffraction from a square grid in the pupil plane, can be used to locate the star and extract its spectrum. We describe the techniques implemented into the GPI Data Reduction Pipeline to measure the properties of the satellite spots and discuss the precision of the reconstructed astrometry and spectrophotometry of the occulted star. We find the astrometric precision of the satellite spots in an $H$-band datacube to be $0.05$ pixels and is best when individual satellite spots have a signal to noise ratio (SNR) of $> 20$. In regards to satellite spot spectrophotometry, we find that the total flux from the satellite spots is stable to $\sim 7\%$ and scales linearly with central star brightness and that the shape of the satellite spot spectrum varies on the $2\%$ level.

We present the results of both laboratory and on sky astrometric characterization of the Gemini Planet Imager (GPI). This characterization includes measurement of the pixel scale of the integral field spectrograph (IFS), the position of the detector with respect to north, and optical distortion. Two of these three quantities (pixel scale and distortion) were measured in the laboratory using two transparent grids of spots, one with a square pattern and the other with a random pattern. The pixel scale in the laboratory was also estimate using small movements of the artificial star unit (ASU) in the GPI adaptive optics system. On sky, the pixel scale and the north angle are determined using a number of known binary or multiple systems and Solar System objects, a subsample of which had concurrent measurements at Keck Observatory. Our current estimate of the GPI pixel scale is 14.14 $\pm$ 0.01 millarcseconds/pixel, and the north angle is -1.00 $\pm$ 0.03$\deg$. Distortion is shown to be small, with an average positional residual of 0.26 pixels over the field of view, and is corrected using a 5th order polynomial. We also present results from Monte Carlo simulations of the GPI Exoplanet Survey (GPIES) assuming GPI achieves ~1 milliarcsecond relative astrometric precision. We find that with this precision, we will be able to constrain the eccentricities of all detected planets, and possibly determine the underlying eccentricity distribution of widely separated Jovians.

We extend the non-equilibrium model for the chemical and thermal evolution of diffuse interstellar gas presented in Richings et al. (2014) to account for shielding from the UV radiation field. We attenuate the photochemical rates by dust and by gas, including absorption by HI, H2, HeI, HeII and CO where appropriate. We then use this model to investigate the dominant cooling and heating processes in interstellar gas as it becomes shielded from the UV radiation. We consider a one-dimensional plane-parallel slab of gas irradiated by the interstellar radiation field, either at constant density and temperature or in thermal and pressure equilibrium. The dominant thermal processes tend to form three distinct regions in the clouds. At low column densities cooling is dominated by ionised metals such as SiII, FeII, FeIII and CII, which are balanced by photoheating, primarily from HI. Once the hydrogen-ionising radiation becomes attenuated by neutral hydrogen, photoelectric dust heating dominates, while CII becomes dominant for cooling. Finally, dust shielding triggers the formation of CO and suppresses photoelectric heating. The dominant coolants in this fully shielded region are H2 and CO. The column density of the HI-H2 transition predicted by our model is lower at higher density (or at higher pressure for gas clouds in pressure equilibrium) and at higher metallicity, in agreement with previous PDR models. We also compare the HI-H2 transition in our model to two prescriptions for molecular hydrogen formation that have been implemented in hydrodynamic simulations.

The discovery of the first methane brown dwarf provides a framework for describing the important advances in both fundamental physics and astrophysics that are due to the study of companions of stars. I present a few highlights of the history of this subject along with details of the discovery of the brown dwarf Gliese 229B. The nature of companions of stars is discussed with an attempt to avoid biases induced by anthropocentric nomenclature. With the newer types of remote reconnaissance of nearby stars and their systems of companions, an exciting and perhaps even more profound set of contributions to science is within reach in the near future. This includes an exploration of the diversity of planets in the universe and perhaps soon the first solid evidence for biological activity outside our Solar System.

We examine the statistics of the low-redshift Lyman-alpha forest from smoothed particle hydrodynamic simulations in light of recent improvements in the estimated evolution of the cosmic ultraviolet background (UVB) and recent observations from the Cosmic Origins Spectrograph (COS). We find that the value of the metagalactic photoionization rate required by our simulations to match the observed properties of the low-redshift Lyman-alpha forest is a factor of 5 larger than the value predicted by state-of-the art models for the evolution of this quantity. This mismatch results in the mean flux decrement of the Lyman-alpha forest being underpredicted by at least a factor of 2 (a 10-sigma discrepancy with observations) and a column density distribution of Lyman-alpha forest absorbers systematically and significantly elevated compared to observations over nearly two decades in column density. We examine potential resolutions to this mismatch and find that either conventional sources of ionizing photons (galaxies and quasars) must be significantly elevated relative to current observational estimates or our theoretical understanding of the low-redshift universe is in need of substantial revision.

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.

The Gemini Planet Imager (GPI) is a dedicated facility for directly imaging and spectroscopically characterizing extrasolar planets. It combines a very high-order adaptive optics system, a diffraction-suppressing coronagraph, and an integral field spectrograph with low spectral resolution but high spatial resolution. Every aspect of GPI has been tuned for maximum sensitivity to faint planets near bright stars. During first light observations, we achieved an estimated H band Strehl ratio of 0.89 and a 5-sigma contrast of $10^6$ at 0.75 arcseconds and $10^5$ at 0.35 arcseconds. Observations of Beta Pictoris clearly detect the planet, Beta Pictoris b, in a single 60-second exposure with minimal post-processing. Beta Pictoris b is observed at a separation of $434 \pm 6$ milli-arcseconds and position angle $211.8 \pm 0.5$ deg. Fitting the Keplerian orbit of Beta Pic b using the new position together with previous astrometry gives a factor of three improvement in most parameters over previous solutions. The planet orbits at a semi-major axis of $9.0^{+0.8}_{-0.4}$ AU near the 3:2 resonance with the previously-known 6 AU asteroidal belt and is aligned with the inner warped disk. The observations give a 4% posterior probability of a transit of the planet in late 2017.

An accurate treatment of the multiphase interstellar medium (ISM) in hydrodynamic galaxy simulations requires that we follow not only the thermal evolution of the gas, but also the evolution of its chemical state, including its molecular chemistry, without assuming chemical (including ionisation) equilibrium. We present a reaction network that can be used to solve for this thermo-chemical evolution. Our model follows the evolution of all ionisation states of the 11 elements that dominate the cooling rate, along with important molecules such as H2 and CO, and the intermediate molecular species that are involved in their formation (20 molecules in total). We include chemical reactions on dust grains, thermal processes involving dust, cosmic ray ionisation and heating and photochemical reactions. We focus on conditions typical for the diffuse ISM, with densities of 10^-2 cm^-3 < nH < 10^4 cm^-3 and temperatures of 10^2 K < T < 10^4 K, and we consider a range of radiation fields, including no UV radiation. In this paper we consider only gas that is optically thin, while paper II considers gas that becomes shielded from the radiation field. We verify the accuracy of our model by comparing chemical abundances and cooling functions in chemical equilibrium with the photoionisation code Cloudy. We identify the major coolants in diffuse interstellar gas to be CII, SiII and FeII, along with OI and H2 at densities nH > 10^2 cm^-3. Finally, we investigate the impact of non-equilibrium chemistry on the cooling functions of isochorically or isobarically cooling gas. We find that, at T < 10^4 K, recombination lags increase the electron abundance above its equilibrium value at a given temperature, which can enhance the cooling rate by up to two orders of magnitude. The cooling gas also shows lower H2 abundances than in equilibrium, by up to an order of magnitude.

We present a budget and accounting of metals in and around star-forming galaxies at $z\sim 0$. We combine empirically derived star formation histories with updated supernova and AGB yields and rates to estimate the total mass of metals produced by galaxies with present-day stellar mass of $10^{9.3}$--$10^{11.6} M_{\odot}$. On the accounting side of the ledger, we show that a surprisingly constant 20--25% mass fraction of produced metals remain in galaxies' stars, interstellar gas and interstellar dust, with little dependence of this fraction on the galaxy stellar mass (omitting those metals immediately locked up in remnants). Thus, the bulk of metals are outside of galaxies, produced in the progenitors of today's $L^*$ galaxies. The COS-Halos survey is uniquely able to measure the mass of metals in the circumgalactic medium (to impact parameters of $< 150$ kpc) of low-redshift $\sim L^*$ galaxies. Using these data, we map the distribution of CGM metals as traced by both the highly ionized OVI ion and a suite of low-ionization species; combined with constraints on circumgalactic dust and hotter X-ray emitting gas out to similar impact parameters, we show that $\sim 40$% of metals produced by $M_{\star}\sim 10^{10}M_{\odot}$ galaxies can be easily accounted for out to 150 kpc. With the current data, we cannot rule out a constant mass of metals within this fixed physical radius. This census provides a crucial boundary condition for the eventual fate of metals in galaxy evolution models.

We use high-resolution cosmological zoom simulations with ~200 pc resolution at z = 2 and various prescriptions for galactic outflows in order to explore the impact of winds on the morphological, dynamical, and structural properties of eight individual galaxies with halo masses ~ 10^11--2x10^12 Msun at z = 2. We present a detailed comparison to spatially and spectrally resolved H{\alpha} and other observations of z ~ 2 galaxies. We find that simulations without winds produce massive, compact galaxies with low gas fractions, super-solar metallicities, high bulge fractions, and much of the star formation concentrated within the inner kpc. Strong winds are required to maintain high gas fractions, redistribute star-forming gas over larger scales, and increase the velocity dispersion of simulated galaxies, more in agreement with the large, extended, turbulent disks typical of high-redshift star-forming galaxies. Winds also suppress early star formation to produce high-redshift cosmic star formation efficiencies in better agreement with observations. Sizes, rotation velocities, and velocity dispersions all scale with stellar mass in accord with observations. Our simulations produce a diversity of morphological characteristics - among our three most massive galaxies, we find a quiescent grand-design spiral, a very compact star-forming galaxy, and a clumpy disk undergoing a minor merger; the clumps are evident in H{\alpha} but not in the stars. Rotation curves are generally slowly rising, particularly when calculated using azimuthal velocities rather than enclosed mass. Our results are broadly resolution-converged. These results show that cosmological simulations including outflows can produce disk galaxies similar to those observed during the peak epoch of cosmic galaxy growth.

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.

We study the physical conditions of the circum-galactic medium (CGM) around z=0.25 galaxies as traced by HI and metal line absorption, using cosmological hydrodynamic simulations that include galactic outflows. Using lines of sight targeted at impact parameters from 10 kpc to 1 Mpc around galaxies with halo masses from 10^11-10^13 M_solar, we study the physical conditions and their variation with impact parameter b and line-of-sight velocity delta v in the CGM as traced by HI, MgII, SiIV, CIV, OVI, and NeVIII absorbers. All ions show a strong excess of absorption near galaxies compared to random lines of sight. The excess continues beyond 1 Mpc, reflecting the correlation of metal absorption with large-scale structure. Absorption is particularly enhanced within about v<300 km/sec and roughly 300 kpc of galaxies (with distances somewhat larger for the highest ion), approximately delineating the CGM; this range contains the majority of global metal absorption. Low ions like MgII and SiIV predominantly arise in denser gas closer to galaxies and drop more rapidly with b, while high ions OVI and NeVIII trace more diffusely distributed gas with a comparatively flat radial profile; CIV is intermediate. All ions predominantly trace T~10^4-4.5 K photo-ionised gas at all b, but when hot CGM gas is present (mostly in larger halos), we see strong collisionally-ionised OVI and NeVIII at b <= 100 kpc. Larger halo masses generally produce more absorption, though overall the trends are not as strong as that with impact parameter. These findings arise using our favoured outflow scalings as expected for momentum-driven winds; with no winds, the CGM gas remains mostly unenriched, while our outflow model with a constant velocity and mass loading factor produce hotter, more widely dispersed metals.

We examine how HI and metal absorption lines within low-redshift galaxy halos trace the dynamical state of circumgalactic gas, using cosmological hydrodynamic simulations that include a well-vetted heuristic model for galactic outflows. We categorize inflowing, outflowing, and ambient gas based on its history and fate as tracked in our simulation. Following our earlier work showing that the ionisation level of absorbers was a primary factor in determining the physical conditions of absorbing gas, we show here that it is also a governing factor for its dynamical state. Low-ionisation metal absorbers (e.g. MgII) tend to arise in gas that will fall onto galaxies within several Gyr, while high-ionisation metal absorbers (e.g. OVI) generally trace material that was deposited by outflows many Gyr ago. Inflowing gas is dominated by enriched material that was previously ejected in an outflow, hence accretion at low redshifts is typically substantially enriched. Recycling wind material is preferentially found closer to galaxies, and is more dominant in lower-mass halos since high-mass halos have more hot gas that is able to support itself against infall. Low-mass halos also tend to re-eject more of their accreted material, owing to our outflow prescription that employs higher mass loading factors for lower-mass galaxies. Typical HI absorbers trace unenriched ambient material that is not participating in the baryon cycle, but stronger HI absorbers arise in cool, enriched inflowing gas. Instantaneous radial velocity measures of absorbers are generally poor at distinguishing between inflowing and outflowing gas, except in the case of very recent outflows. These results suggest that probing halo gas using a range of absorbers can provide detailed information about the amount and physical conditions of material that is participating in the baryon cycle.

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.

The Project 1640 instrument on the 200-inch Hale telescope at Palomar Observatory is a coronagraphic instrument with an integral field spectrograph at the back end, designed to find young, self-luminous planets around nearby stars. To reach the necessary contrast for this, the PALM-3000 adaptive optics system corrects for fast atmospheric speckles, while CAL, a phase-shifting interferometer in a Mach-Zehnder configuration, measures the quasistatic components of the complex electric field in the pupil plane following the coronagraphic stop. Two additional sensors measure and control low-order modes. These field measurements may then be combined with a system model and data taken separately using a white-light source internal to the AO system to correct for both phase and amplitude aberrations. Here, we discuss and demonstrate the procedure to maintain a half-plane dark hole in the image plane while the spectrograph is taking data, including initial on-sky performance.

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 model the time-dependent evolution of metal-enriched intergalactic and circumgalactic gas exposed to the fluctuating radiation field from an active galactic nucleus (AGN). We consider diffuse gas densities (n_H=10^-5-10^-3 cm^-3) exposed to the extra-galactic background (EGB) and initially in thermal equilibrium (T \sim 10^4-10^4.5 K). Once the proximate AGN field turns on, additional photo-ionisation rapidly ionises the HI and metals. The enhanced AGN radiation field turns off after a typical AGN lifetime (tau_AGN=1-20 Myr) and the field returns to the EGB intensity, but the metals remain out of ionisation equilibrium for timescales that can significantly exceed tau_AGN. We define this phase as the AGN proximity zone "fossil" phase and show that high ionisation stages (e.g. OVI, NeVIII, MgX) are in general enhanced, while the abundances of low ions are reduced. In contrast, HI re-equilibrates rapidly (<<tau_AGN) owing to its low neutral fraction at diffuse densities. We demonstrate that metal column densities of intervening gas observed in absorption in quasar sight lines are significantly affected by delayed recombination for a wide range of densities, metallicities, and AGN strengths, lifetimes, and duty cycles. We model the exceptionally strong z=0.9 NeVIII absorbers observed by Tripp et al. (2011) as arising in a possible fossil zone or near a recently turned-on AGN and we demonstrate that at low redshift even moderate strength AGN could significantly enhance the high-ion metal columns in the circumgalactic media of galaxies observed without active AGN. Fossil proximity zones may be very important during the quasar era, z \sim 2-5. AGN proximity zone fossils allow a whole new class of non-equilibrium solutions that may be applicable to a large fraction of observed metal absorbers and which could potentially change the inferred physical conditions and masses of diffuse gases.

Simulations of the formation of galaxies, as well as ionisation models used to interpret observations of quasar absorption lines, generally either assume ionisation equilibrium or ignore the presence of the extra-galactic background (EGB) radiation. We introduce a method to compute the non-equilibrium ionisation and cooling of diffuse gas exposed to the EGB. Our method iterates the ionisation states of the 11 elements that dominate the cooling and uses tabulated ion-by-ion cooling and photo-heating efficiencies to update the temperature of the gas. Our reaction network includes radiative and di-electric recombination, collisional ionisation, photo-ionisation, Auger ionisation, and charge transfer. We verify that our method reproduces published results for collisional equilibrium, collisional non-equilibrium, and photo-ionised equilibrium. Non-equilibrium effects can become very important in cooling gas, particularly below 10^6 K. Photo-ionisation and non-equilibrium effects both tend to boost the degree of ionisation and to reduce cooling efficiencies. The effect of the EGB is larger for lower densities. Hence, photo-ionisation affects cooling more under isochoric than under isobaric conditions. Non-equilibrium effects are smaller in the presence of the EGB and are thus overestimated when using collisional-only processes. The inclusion of the EGB alters the observational diagnostics of diffuse, metal-enriched gas (e.g. metal absorption lines) even more significantly than the cooling efficiencies. We argue that the cooling efficiency should be considered if ionisation models are used to infer physical conditions from observed line ratios, as the a priori probability of observing gas is lower if its cooling time is shorter. We provide on-line tables of ionisation fractions and cooling efficiencies for equilibrium and non-equilibrium scenarios. (abridged)

Kappa Andromedae is a B9IVn star at 52 pc for which a faint substellar companion separated by 55 AU was recently announced. In this work, we present the first spectrum of the companion, "kappa And B," using the Project 1640 high-contrast imaging platform. Comparison of our low-resolution YJH-band spectra to empirical brown dwarf spectra suggests an early-L spectral type. Fitting synthetic spectra from PHOENIX model atmospheres to our observed spectrum allows us to constrain the effective temperature to ~2000K, as well as place constraints on the companion surface gravity. Further, we use previously reported log(g) and effective temperature measurements of the host star to argue that the kappa And system has an isochronal age of 220 +/- 100 Myr, older than the 30 Myr age reported previously. This interpretation of an older age is corroborated by the photometric properties of kappa And B, which appear to be marginally inconsistent with other 10-100 Myr low-gravity L-dwarfs for the spectral type range we derive. In addition, we use Keck aperture masking interferometry combined with published radial velocity measurements to rule out the existence of any tight stellar companions to kappa And A that might be responsible for the system's overluminosity. Further, we show that luminosity enhancements due to a nearly "pole-on" viewing angle coupled with extremely rapid rotation is unlikely. Kappa And A is thus consistent with its slightly evolved luminosity class (IV) and we propose here that kappa And, with a revised age of 220 +/- 100 Myr, is an interloper to the 30 Myr Columba association with which it was previously associated. The photometric and spectroscopic evidence for kappa And B combined with our re-assesment of the system age implies a substellar companion mass of 50^{+16}_{-13} Jupiter Masses, consistent with a brown dwarf rather than a planetary mass companion.

We use a radiation hydrodynamic simulation of the hydrogen reionization epoch to study OI absorbers at z~6. The intergalactic medium (IGM) is reionized before it is enriched, hence OI absorption originates within dark matter halos. The predicted abundance of OI absorbers is in reasonable agreement with observations. At z=10, roughly 70% of sightlines through atomically-cooled halos encounter a visible (N_OI > 10^14 cm^-2) column. Reionization ionizes and removes gas from halos less massive than 10^8.4 M_0, but 20% of sightlines through more massive halos encounter visible columns even at z=5. The mass scale of absorber host halos is 10-100 times smaller than the halos of Lyman break galaxies and Lyman-alpha emitters, hence absorption probes the dominant ionizing sources more directly. OI absorbers have neutral hydrogen columns of 10^19-10^21 cm^-2, suggesting a close resemblance between objects selected in OI and HI absorption. Finally, the absorption in the foreground of the z=7.085 quasar ULASJ1120+0641 cannot originate in a dark matter halo because halo gas at the observed HI column density is enriched enough to violate the upper limits on the OI column. By contrast, gas at less than one third the cosmic mean density satisfies the constraints. Hence the foreground absorption likely originates in the IGM.

We examine the global HI properties of galaxies in quarter-billion particle cosmological simulations using Gadget-2, focusing on how galactic outflows impact HI content. We consider four outflow models, including a new one (ezw) motivated by recent interstellar medium simulations in which the wind speed and mass loading factor scale as expected for momentum-driven outflows for larger galaxies and energy-driven outflows for dwarfs (sigma<75 km/s). To obtain predicted HI masses, we employ a simple but effective local correction for particle self-shielding, and an observationally-constrained transition from neutral to molecular hydrogen. Our ezw simulation produces an HI mass function whose faint-end slope of -1.3 agrees well with observations from the ALFALFA survey; other models agree less well. Satellite galaxies have a bimodal distribution in HI fraction versus halo mass, with smaller satellites and/or those in larger halos more often being HI-deficient. At a given stellar mass, HI content correlates with star formation rate and inversely correlates with metallicity, as expected if driven by stochasticity in the accretion rate. To higher redshifts, massive HI galaxies disappear and the mass function steepens. The global cosmic HI density conspires to remain fairly constant from z~5-0, but the relative contribution from smaller galaxies increases with redshift.

We obtained spectra, in the wavelength range \lambda = 995 - 1769 nm, of all four known planets orbiting the star HR 8799. Using the suite of instrumentation known as Project 1640 on the Palomar 5-m Hale Telescope, we acquired data at two epochs. This allowed for multiple imaging detections of the companions and multiple extractions of low-resolution (R ~ 35) spectra. Data reduction employed two different methods of speckle suppression and spectrum extraction, both yielding results that agree. The spectra do not directly correspond to those of any known objects, although similarities with L and T-dwarfs are present, as well as some characteristics similar to planets such as Saturn. We tentatively identify the presence of CH_4 along with NH_3 and/or C_2H_2, and possibly CO_2 or HCN in varying amounts in each component of the system. Other studies suggested red colors for these faint companions, and our data confirm those observations. Cloudy models, based on previous photometric observations, may provide the best explanation for the new data presented here. Notable in our data is that these presumably co-eval objects of similar luminosity have significantly different spectra; the diversity of planets may be greater than previously thought. The techniques and methods employed in this paper represent a new capability to observe and rapidly characterize exoplanetary systems in a routine manner over a broad range of planet masses and separations. These are the first simultaneous spectroscopic observations of multiple planets in a planetary system other than our own.

We present adaptive optics photometry and spectra in the JHKL-bands along with high spectral resolution K-band spectroscopy for each component of the Z Canis Majoris system. Our high angular resolution photometry of this very young (<1 Myr) binary, comprised of an FU Ori object and a Herbig Ae/Be star, were gathered shortly after the 2008 outburst while our high resolution spectroscopy was gathered during a quiescent phase. Our photometry conclusively determine that the outburst was due solely to the embedded Herbig Ae/Be member, supporting results from earlier works, and that the optically visible FU Ori component decreased slightly (~30%) in luminosity during the same period, consistent with previous works on the variability of FU Ori type systems. Further, our high-resolution K-band spectra definitively demonstrate that the 2.294 micron CO absorption feature seen in composite spectra of the system is due solely to the FU Ori component, while a prominent CO emission feature at the same wavelength, long suspected to be associated with the innermost regions of a circumstellar accretion disk, can be assigned to the Herbig Ae/Be member. These findings are in contrast to previous analyses (e.g. Malbet et al 2010, Benisty et al. 2010) of this complex system which assigned the CO emission to the FU Ori component.

We report low resolution near infrared spectroscopic observations of the eruptive star FU Orionis using the Integral Field Spectrograph Project 1640 installed at the Palomar Hale telescope. This work focuses on elucidating the nature of the faint source, located 0.5" south of FU Ori, and identified in 2003 as FU Ori S. We first use our observations in conjunction with published data to demonstrate that the two stars are indeed physically associated and form a true binary pair. We then proceed to extract J and H band spectro-photometry using the damped LOCI algorithm, a reduction method tailored for high contrast science with IFS. This is the first communication reporting the high accuracy of this technique, pioneered by the Project 1640 team, on a faint astronomical source. We use our low resolution near infrared spectrum in conjunction with 10.2 micron interferometric data to constrain the infrared excess of FU Ori S. We then focus on estimating the bulk physical properties of FU Ori S. Our models lead to estimates of an object heavily reddened, A_V =8-12, with an effective temperature of ~ 4000-6500 K . Finally we put these results in the context of the FU Ori N-S system and argue that our analysis provides evidence that FU Ori S might be the more massive component of this binary system

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.

Recent studies have found a dramatic difference between the observed number density evolution of low mass galaxies and that predicted by semi-analytic models. While models accurately reproduce the z=0 number density, they require that the evolution occurs rapidly at early times, which is incompatible with the strong late evolution found observationally. We report here the same discrepancy in two state-of-the-art cosmological hydrodynamical simulations, which is evidence that the problem is fundamental. We search for the underlying cause of this problem using two complementary methods. Firstly, we look for evidence of a different history of today's low mass galaxies in models and observations and we find that the models yield too few young, strongly star-forming galaxies. Secondly, we construct a toy model to link the observed evolution of specific star formation rates (sSFR) with the evolution of the galaxy stellar mass function. We infer from this model that a key problem in both semi-analytic and hydrodynamical models is the presence of a positive instead of a negative correlation between sSFR and stellar mass. A similar positive correlation is found between the specific dark matter halo accretion rate and the halo mass, indicating that model galaxies are growing in a way that follows the growth of their host haloes too closely. It therefore appears necessary to find a mechanism that decouples the growth of low mass galaxies, which occurs primarily at late times, from the growth of their host haloes, which occurs primarily at early times. We argue that the current form of star-formation driven feedback implemented in most galaxy formation models is unlikely to achieve this goal, owing to its fundamental dependence on host halo mass and time. [Abridged]

Small-angle coronagraphy is technically and scientifically appealing because it enables the use of smaller telescopes, allows covering wider wavelength ranges, and potentially increases the yield and completeness of circumstellar environment - exoplanets and disks - detection and characterization campaigns. However, opening up this new parameter space is challenging. Here we will review the four posts of high contrast imaging and their intricate interactions at very small angles (within the first 4 resolution elements from the star). The four posts are: choice of coronagraph, optimized wavefront control, observing strategy, and post-processing methods. After detailing each of the four foundations, we will present the lessons learned from the 10+ years of operations of zeroth and first-generation adaptive optics systems. We will then tentatively show how informative the current integration of second-generation adaptive optics system is, and which lessons can already be drawn from this fresh experience. Then, we will review the current state of the art, by presenting world record contrasts obtained in the framework of technological demonstrations for space-based exoplanet imaging and characterization mission concepts. Finally, we will conclude by emphasizing the importance of the cross-breeding between techniques developed for both ground-based and space-based projects, which is relevant for future high contrast imaging instruments and facilities in space or on the ground.

We used the Project 1640 near-infrared coronagraph and integral field spectrograph to observe 19 young solar type stars. Five of these stars are known binary stars and we detected the late-type secondaries and were able to measure their JH spectra with a resolution of R\sim30. The reduced, extracted, and calibrated spectra were compared to template spectra from the IRTF spectral library. With this comparison we test the accuracy and consistency of spectral type determination with the low-resolution near-infrared spectra from P1640. Additionally, we determine effective temperature and surface gravity of the companions by fitting synthetic spectra calculated with the PHOENIX model atmosphere code. We also present several new epochs of astrometry of each of the systems. Together these data increase our knowledge and understanding of the stellar make up of these systems. In addition to the astronomical results, the analysis presented helps validate the Project 1640 data reduction and spectral extraction processes and the utility of low-resolution, near-infrared spectra for characterizing late-type companions in multiple systems.

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.

High-contrast imaging instruments are now being equipped with integral field spectrographs (IFS) to facilitate the detection and characterization of faint substellar companions. Algorithms currently envisioned to handle IFS data, such as the Locally Optimized Combination of Images (LOCI) algorithm, rely upon aggressive point-spread-function (PSF) subtraction, which is ideal for initially identifying companions but results in significantly biased photometry and spectroscopy due to unwanted mixing with residual starlight. This spectro-photometric issue is further complicated by the fact that algorithmic color response is a function of the companion's spectrum, making it difficult to calibrate the effects of the reduction without using iterations involving a series of injected synthetic companions. In this paper, we introduce a new PSF calibration method, which we call "damped LOCI", that seeks to alleviate these concerns. By modifying the cost function that determines the weighting coefficients used to construct PSF reference images, and also forcing those coefficients to be positive, it is possible to extract companion spectra with a precision that is set by calibration of the instrument response and transmission of the atmosphere, and not by post-processing. We demonstrate the utility of this approach using on-sky data obtained with the Project 1640 IFS at Palomar. Damped-LOCI does not require any iterations on the underlying spectral type of the companion, nor does it rely upon priors involving the chromatic and statistical properties of speckles. It is a general technique that can readily be applied to other current and planned instruments that employ IFS's.

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.

We present an analytic formalism that describes the evolution of the stellar, gas, and metal content of galaxies. It is based on the idea, inspired by hydrodynamic simulations, that galaxies live in a slowly-evolving equilibrium between inflow, outflow, and star formation. We argue that this formalism broadly captures the behavior of galaxy properties evolving in simulations. The resulting equilibrium equations for the star formation rate, gas fraction, and metallicity depend on three key free parameters that represent ejective feedback, preventive feedback, and re-accretion of ejected material. We schematically describe how these parameters are constrained by models and observations. Galaxies perturbed off the equilibrium relations owing to inflow stochasticity tend to be driven back towards equilibrium, such that deviations in star formation rate at a given mass are correlated with gas fraction and anti-correlated with metallicity. After an early gas accumulation epoch, quiescently star-forming galaxies are expected to be in equilibrium over most of cosmic time. The equilibrium model provides a simple intuitive framework for understanding the cosmic evolution of galaxy properties, and centrally features the cycle of baryons between galaxies and surrounding gas as the driver of galaxy growth.

We investigate the metallicity evolution and content of the intergalactic medium (IGM) and galactic halo gas from z=2->0 using 110-million particle cosmological hydrodynamic simulations. We focus on the detectability and physical properties of UV resonance metal-line absorbers observable with Hubble's Cosmic Origins Spectrograph (COS). We confirm that galactic superwind outflows are required to enrich the IGM to observed levels down to z=0 using three wind prescriptions contrasted to a no-wind simulation. Our favoured momentum-conserved wind prescription deposits metals closer to galaxies owing to its moderate energy input, while the more energetic constant wind model enriches the warm-hot IGM 6.4x more. Despite these significant differences, all wind models produce metal-line statistics within a factor of two of existing observations. This is because OVI, CIV, SiIV, and NeVIII absorbers primarily arise from T<10^5 K, photo-ionised gas that is enriched to similar levels in the three feedback schemes. OVI absorbers trace the diffuse phase with delta<100, which is enriched to ~1/50 Zsol at z=0, although the absorbers themselves usually exceed 0.3 Zsol and arise from inhomogeneously distributed, un-mixed winds. CIV and SiIV absorbers trace primarily T~10^4 K gas inside haloes. We predict COS will observe a population of NeVIII photo-ionised absorbers tracing T<10^5 K, delta~10 gas. MgX and SiXII are rarely detected in COS S/N=30 simulated sight lines although detections trace T=10^(6-7) K halo gas. In general, the IGM is enriched in an outside-in manner, where wind-blown metals released at higher redshift reach lower overdensities, resulting in higher ionisation species tracing lower-density, older metals. At z=0, the 90% of baryons outside of galaxies are enriched to 0.096 Zsol, but the 65% of unbound baryons in the IGM have 0.018 Zsol and contain only 4% of all metals.

Many observed massive star-forming z\approx2 galaxies are large disks that exhibit irregular morphologies, with \sim1kpc, \sim10^(8-10)Msun clumps. We present the largest sample to date of high-resolution cosmological SPH simulations that zoom-in on the formation of individual M*\sim10^(10.5)Msun galaxies in \sim10^(12)Msun halos at z\approx2. Our code includes strong stellar feedback parameterized as momentum-driven galactic winds. This model reproduces many characteristic features of this observed class of galaxies, such as their clumpy morphologies, smooth and monotonic velocity gradients, high gas fractions (f_g\sim50%) and high specific star-formation rates (\gtrsim1Gyr^(-1)). In accord with recent models, giant clumps (Mclump\sim(5x10^8-10^9)Msun) form in-situ via gravitational instabilities. However, the galactic winds are critical for their subsequent evolution. The giant clumps we obtain are short-lived and are disrupted by wind-driven mass loss. They do not virialise or migrate to the galaxy centers as suggested in recent work neglecting strong winds. By phenomenologically implementing the winds that are observed from high-redshift galaxies and in particular from individual clumps, our simulations reproduce well new observational constraints on clump kinematics and clump ages. In particular, the observation that older clumps appear closer to their galaxy centers is reproduced in our simulations, as a result of inside-out formation of the disks rather than inward clump migration.

Massive galaxies today typically are not forming stars despite being surrounded by hot gaseous halos with short central cooling times. This likely owes to some form of "quenching feedback" such as merger-driven quasar activity or radio jets emerging from central black holes. Here we implement heuristic prescriptions for these phenomena on-the-fly within cosmological hydrodynamic simulations. We constrain them by comparing to observed luminosity functions and color-magnitude diagrams from SDSS. We find that quenching from mergers alone does not produce a realistic red sequence, because 1 - 2 Gyr after a merger the remnant accretes new fuel and star formation reignites. In contrast, quenching by continuously adding thermal energy to hot gaseous halos quantitatively matches the red galaxy luminosity function and produces a reasonable red sequence. Small discrepancies remain - a shallow red sequence slope suggests that our models underestimate metal production or retention in massive red galaxies, while a deficit of massive blue galaxies may reflect the fact that observed heating is intermittent rather than continuous. Overall, injection of energy into hot halo gas appears to be a necessary and sufficient condition to broadly produce red and dead massive galaxies as observed.

We use cosmological hydrodynamic simulations to investigate how inflows, star formation, and outflows govern the the gaseous and metal content of galaxies. In our simulations, galaxy metallicities are established by a balance between inflows and outflows as governed by the mass outflow rate, implying that the mass-metallicity relation reflects how the outflow rate varies with stellar mass (M*). Gas content is set by a competition between inflow into and gas consumption within the ISM, the latter being governed by the SF law, while the former is impacted by both wind recycling and preventive feedback. Stochastic variations in the inflow rate move galaxies off the equilibrium M*-Z and Z*-fgas relations in a manner correlated with star formation rate, and the scatter is set by the timescale to re-equilibrate. The evolution of both relations from z=3-0 is slow, as individual galaxies tend to evolve mostly along the relations. Gas fractions at a given M* slowly decrease with time because the cosmic inflow rate diminishes faster than the consumption rate, while metallicities slowly increase as infalling gas becomes more enriched. Observations from z~3-0 are better matched by simulations employing momentum-driven wind scalings rather than constant wind speeds, but all models predict too low gas fractions at low masses and too high metallicities at high M*. All our models reproduce observed second-parameter trends of the mass-metallicity relation with star formation rate and environment, indicating that these are a consequence of equilibrium and not feedback. Overall, the analytical framework of our equilibrium scenario broadly captures the relevant physics establishing the galaxy gas and metal content in simulations, which suggests that the cycle of baryonic inflows and outflows centrally governs the cosmic evolution of these properties in typical star-forming galaxies.

Project 1640 is a high contrast near-infrared instrument probing the vicinities of nearby stars through the unique combination of an integral field spectrograph with a Lyot coronagraph and a high-order adaptive optics system. The extraordinary data reduction demands, similar those which several new exoplanet imaging instruments will face in the near future, have been met by the novel software algorithms described herein. The Project 1640 Data Cube Extraction Pipeline (PCXP) automates the translation of 3.8*10^4 closely packed, coarsely sampled spectra to a data cube. We implement a robust empirical model of the spectrograph focal plane geometry to register the detector image at sub-pixel precision, and map the cube extraction. We demonstrate our ability to accurately retrieve source spectra based on an observation of Saturn's moon Titan.

We examine the growth of the stellar content of galaxies from z=3-0 in cosmological hydrodynamic simulations incorporating parameterised galactic outflows. Without outflows, galaxies overproduce stellar masses (M*) and star formation rates (SFRs) compared to observations. Winds introduce a three-tier form for the galaxy stellar mass and star formation rate functions, where the middle tier depends on differential (i.e. mass-dependent) recycling of ejected wind material back into galaxies. A tight M*-SFR relation is a generic outcome of all these simulations, and its evolution is well-described as being powered by cold accretion, although current observations at z>2 suggest that star formation in small early galaxies must be highly suppressed. Roughly one-third of z=0 galaxies at masses below M^* are satellites, and star formation in satellites is not much burstier than in centrals. All models fail to suppress star formation and stellar mass growth in massive galaxies at z<2, indicating the need for an external quenching mechanism such as black hole feedback. All models also fail to produce dwarfs as young and rapidly star-forming as observed. An outflow model following scalings expected for momentum-driven winds broadly matches observed galaxy evolution around M^* from z=0-3, which is a significant success since these galaxies dominate cosmic star formation, but the failures at higher and lower masses highlight the challenges still faced by this class of models. We argue that central star-forming galaxies are well-described as living in a slowly-evolving equilibrium between inflows from gravity and recycled winds, star formation, and strong and ubiquitous outflows that regulate how much inflow forms into stars. Star-forming galaxy evolution is thus primarily governed by the continual cycling of baryons between galaxies and intergalactic gas.

Recently Menard et al. detected a subtle but systematic change in the mean color of quasars as a function of their projected separation from foreground galaxies, extending to comoving separations of ~10Mpc/h, which they interpret as a signature of reddening by intergalactic dust. We present theoretical models of this remarkable observation, using SPH cosmological simulations of a (50Mpc/h)^3 volume. Our primary model uses a simulation with galactic winds and assumes that dust traces the intergalactic metals. The predicted galaxy-dust correlation function is similar in form to the galaxy-mass correlation function, and reproducing the MSFR data requires a dust-to-metal mass ratio of 0.24, about half the value in the Galactic ISM. Roughly half of the reddening arises in dust that is more than 100Kpc/h from the nearest massive galaxy. We also examine a simulation with no galactic winds, which predicts a much smaller fraction of intergalactic metals (3% vs. 35%) and therefore requires an unphysical dust-to-metal ratio of 2.18 to reproduce the MSFR data. In both models, the signal is dominated by sightlines with E(g-i)=0.001-0.1. The no-wind simulation can be reconciled with the data if we also allow reddening to arise in galaxies up to several x 10^10 Msun. The wind model predicts a mean visual extinction of A_V ~0.0133 mag out to z=0.5, with a sightline-to-sightline dispersion similar to the mean, which could be significant for future supernova cosmology studies. Reproducing the MSFR results in these simulations requires that a large fraction of ISM dust survive its expulsion from galaxies and its residence in the intergalactic medium. Future observational studies that provide higher precision and measure the dependence on galaxy type and environment will allow detailed tests for models of enriched galactic outflows and the survival of IG dust.

Project 1640 is a high-contrast imaging instrument recently commissioned at Palomar observatory. A combination of a coronagraph with an integral field spectrograph (IFS), Project 1640 is designed to detect and characterize extrasolar planets, brown dwarfs, and circumstellar material orbiting nearby stars. In this paper, we present our data processing techniques for improving upon instrument raw sensitivity via the removal of quasi-static speckles. Our approach utilizes the chromatic image diversity provided by the IFS in combination with the locally-optimized combination of images (LOCI) algorithm to suppress the intensity of residual contaminating light in close angular proximity to target stars. We describe the Project 1640 speckle suppression pipeline (PSSP) and demonstrate the ability to detect companions with brightness comparable to and below that of initial speckle intensities using on-sky commissioning data. Our preliminary results indicate that suppression factors of at least one order of magnitude are consistently possible, reaching $5\sigma$ contrast levels of $2.1\times10^{-5}$ at $1\arcsec$ in the H-band in 20 minutes of on-source integration time when non-common-path errors are reasonably well-calibrated. These results suggest that near-infrared contrast levels of order $\approx10^{-7}$ at subarcsecond separations will soon be possible for Project 1640 and similarly designed instruments that receive a diffraction-limited beam corrected by adaptive optics (AO) systems employing deformable mirrors with high actuator-density.

We describe a new instrument that forms the core of a long-term high contrast imaging program at the 200-inch Hale Telescope at Palomar Observatory. The primary scientific thrust is to obtain images and low-resolution spectroscopy of brown dwarfs and young Jovian mass exoplanets in the vicinity of stars within 50 parsecs of the Sun. The instrument is a microlens-based integral field spectrograph integrated with a diffraction limited, apodized-pupil Lyot coronagraph, mounted behind the Palomar adaptive optics system. The spectrograph obtains imaging in 23 channels across the J and H bands (1.06 - 1.78 microns). In addition to obtaining spectra, this wavelength resolution allows suppression of the chromatically dependent speckle noise, which we describe. We have recently installed a novel internal wave front calibration system that will provide continuous updates to the AO system every 0.5 - 1.0 minutes by sensing the wave front within the coronagraph. The Palomar AO system is undergoing an upgrade to a much higher-order AO system ("PALM-3000"): a 3388-actuator tweeter deformable mirror working together with the existing 241-actuator mirror. This system will allow correction with subapertures as small as 8cm at the telescope pupil using natural guide stars. The coronagraph alone has achieved an initial dynamic range in the H-band of 2 X 10^-4 at 1 arcsecond, without speckle noise suppression. We demonstrate that spectral speckle suppression is providing a factor of 10-20 improvement over this bringing our current contrast at an arcsecond to ~2 X 10^-5. This system is the first of a new generation of apodized pupil coronagraphs combined with high-order adaptive optics and integral field spectrographs (e.g. GPI, SPHERE, HiCIAO), and we anticipate this instrument will make a lasting contribution to high contrast imaging in the Northern Hemisphere for years.

The nearby star Alpha Oph (Ras Alhague) is a rapidly rotating A5IV star spinning at ~89% of its breakup velocity. This system has been imaged extensively by interferometric techniques, giving a precise geometric model of the star's oblateness and the resulting temperature variation on the stellar surface. Fortuitously, Alpha Oph has a previously known stellar companion, and characterization of the orbit provides an independent, dynamically-based check of both the host star and the companion mass. Such measurements are crucial to constrain models of such rapidly rotating stars. In this study, we combine eight years of Adaptive Optics imaging data from the Palomar, AEOS, and CFHT telescopes to derive an improved, astrometric characterization of the companion orbit. We also use photometry from these observations to derive a model-based estimate of the companion mass. A fit was performed on the photocenter motion of this system to extract a component mass ratio. We find masses of 2.40^{0.23}_{0.37} solar masses and 0.85^{0.06}_{0.04} solar masses for Alpha Oph A and Alpha Oph B, respectively. Previous orbital studies of this system found a mass too high for this system, inconsistent with stellar evolutionary calculations. Our measurements of the host star mass are more consistent with these evolutionary calculations, but with slightly higher uncertainties. In addition to the dynamically-derived masses, we use IJHK photometry to derive a model-based mass for Alpha Oph B, of 0.77 +/- 0.05 solar masses marginally consistent with the dynamical masses derived from our orbit. Our model fits predict a periastron passage on 2012 April 19, with the two components having a ~50 milliarcsec separation from March to May 2012. A modest amount of interferometric and radial velocity data during this period could provide a mass determination of this star at the few percent level.

The kinematics of damped Lyman alpha absorbers (DLAs) are difficult to reproduce in hierarchical galaxy formation models, particularly the preponderance of wide systems. We investigate DLA kinematics at z=3 using high-resolution cosmological hydrodynamical simulations that include a heuristic model for galactic outflows. Without outflows, our simulations fail to yield enough wide DLAs, as in previous studies. With outflows, predicted DLA kinematics are in much better agreement with observations. Comparing two outflow models, we find that a model based on momentum-driven wind scalings provides the best match to the observed DLA kinematic statistics of Prochaska & Wolfe. In this model, DLAs typically arise a few kpc away from galaxies that would be identified in emission. Narrow DLAs can arise from any halo and galaxy mass, but wide ones only arise in halos with mass >10^11 Mo, from either large central or small satellite galaxies. This implies that the success of this outflow model originates from being most efficient at pushing gas out from small satellite galaxies living in larger halos. This increases the cross-section for large halos relative to smaller ones, thereby yielding wider kinematics. Our simulations do not include radiative transfer effects or detailed metal tracking, and outflows are modeled heuristically, but they strongly suggest that galactic outflows are central to understanding DLA kinematics. An interesting consequence is that DLA kinematics may place constraints on the nature and efficiency of gas ejection from high-z galaxies.

Cosmological hydrodynamic simulations robustly predict that high-redshift galaxy star formation histories (SFHs) are smoothly-rising and vary with mass only by a scale factor. We use our latest simulations to test whether this scenario can account for recent observations at z>=6 from WFC3/IR, NICMOS, and IRAC. Our simulations broadly reproduce the observed ultraviolet (UV) luminosity functions and stellar mass densities and their evolution at z=6-8, all of which are nontrivial tests of the mean SFH. In agreement with observations, simulated galaxies possess blue UV continua owing to young ages (50-150 Myr), low metallicities (0.1-0.5 Zsun), and low dust columns (E(B-V) <= 0.05). Observations imply a near-unity slope in the stellar mass--star formation rate relation at all z=6-8, confirming the prediction that SFH shapes are invariant. Current surveys detect the majority of galaxies with stellar masses exceeding 10^9 Msun and few galaxies less massive than 10^{8.5} Msun, implying that they probe no more than the brightest 30% of the complete star formation and stellar mass densities at z>=6. Finally, we demonstrate that there is no conflict between smoothly-rising SFHs and recent clustering observations. This is because momentum-driven outflows suppress star formation in low-mass halos, leading to overall occupancies of 0.2-0.4 even though the star formation duty cycle is one. This leads to many interesting predictions at z>=4, among them that (1) optically-selected and UV-selected samples largely overlap; (2) few galaxies exhibit significantly suppressed specific star formation rates; and (3) occupancy is constant or increasing with decreasing luminosity. These predictions are in tentative agreement with current observations, but further analysis of existing and upcoming data sets is required in order to test them more thoroughly. (abridged)

The intergalactic medium (IGM) is the dominant reservoir of baryons at all cosmic epochs. We investigate the evolution of the IGM from z=2-0 in 48 Mpc/h, 110-million particle cosmological hydrodynamic simulations using three prescriptions for galactic outflows. We focus on the evolution of IGM physical properties, and how such properties are traced by Ly-alpha absorption as detectable using HST/COS. Our results broadly confirm the canonical picture that most Ly-alpha absorbers arise from highly ionized gas tracing filamentary large-scale structure. Growth of structure causes gas to move from the diffuse photoionized IGM into other cosmic phases, namely stars, cold and hot gas within galaxy halos, and the unbound and shock-heated warm-hot intergalactic medium (WHIM). By today, baryons are roughly equally divided between bound phases (35%), the diffuse IGM (41%), and the WHIM (24%). Here we (re)define the WHIM as gas with overdensities lower than that in halos and temperatures >10^5 K, in order to more closely align it with "missing baryons". When we tune our photoionizing background to match the observed evolution of the Ly-alpha mean flux decrement, we obtain a line count evolution that broadly agrees with available data. We predict a column density distribution slope of -1.70 for our favored momentum-driven wind model, in agreement with recent observations, and it becomes shallower with redshift. With improved statistics, the frequency of strong lines can be a valuable diagnostic of outflows, and our favored wind model matches existing data best among our models. The relationship between column density and physical density is fairly tight from z=2-0, and evolves as rho N_HI^0.74 10^(-0.37z) for diffuse absorbers. Linewidths only loosely reflect the temperature of the absorbing gas, which will hamper attempts to quantify the WHIM using broad Ly-alpha absorbers. [Abridged]

The Lyot project used an optimized Lyot coronagraph with Extreme Adaptive Optics at the 3.63m Advanced Electro-Optical System telescope (AEOS) to observe 86 stars from 2004 to 2007. In this paper we give an overview of the survey results and a statistical analysis of the observed nondetections around 58 of our targets to place constraints on the population of substellar companions to nearby stars. The observations did not detect any companion in the substellar regime. Since null results can be as important as detections, we analyzed each observation to determine the characteristics of the companions that can be ruled out. For this purpose we use a Monte Carlo approach to produce artificial companions, and determine their detectability by comparison with the sensitivity curve for each star. All the non-detection results are combined using a Bayesian approach and we provide upper limits on the population of giant exoplanets and brown dwarfs for this sample of stars. Our nondetections confirm the rarity of brown dwarfs around solar-like stars and we constrain the frequency of massive substellar companions (M>40Mjup) at orbital separation between and 10 and 50 AU to be <20%.

The bimodality in observed present-day galaxy colours has long been a challenge for hierarchical galaxy formation models, as it requires some physical process to quench (and keep quenched) star formation in massive galaxies. Here we examine phenomenological models of quenching by post-processing the star formation histories of galaxies from cosmological hydrodynamic simulations that reproduce observations of star-forming galaxies reasonably well. We consider recipes for quenching based on major mergers, halo mass thresholds, gas temperature thresholds, and variants thereof. We compare the resulting simulated star formation histories to observed g-r colour-magnitude diagrams and red and blue luminosity functions from SDSS. The merger and halo mass quenching scenarios each yield a distinct red sequence and blue cloud of galaxies that are in broad agreement with data, albeit only under rather extreme assumptions. In detail, however, the simulated red sequence slope and amplitude in each scenario is somewhat discrepant, perhaps traceable to low metallicities in simulated galaxies. Merger quenching produces more massive blue galaxies, earlier quenching, and more frosting of young stars; comparing to relevant data tends to favor merger over halo mass quenching. Although physically-motivated quenching models can produce a red sequence, interesting generic discrepancies remain that indicate that additional physics is required to reproduce the star formation and enrichment histories of red and dead galaxies.

We study the nature of rapidly star-forming galaxies at z=2 in cosmological hydrodynamic simulations, and compare their properties to observations of sub-millimetre galaxies (SMGs). We identify simulated SMGs as the most rapidly star-forming systems that match the observed number density of SMGs. In our models, SMGs are massive galaxies sitting at the centres of large potential wells, being fed by smooth infall and gas-rich satellites at rates comparable to their star formation rates (SFR). They are not typically undergoing major mergers that significantly boost their quiescent SFR, but they still often show complex gas morphologies and kinematics. Our simulated SMGs have stellar masses of log M*/Mo~11-11.7, SFRs of ~180-500 Mo/yr, a clustering length of 10 Mpc/h, and solar metallicities. The SFRs are lower than those inferred from far-IR data by a factor of 3, which we suggest may owe to one or more systematic effects in the SFR calibrations. SMGs at z=2 live in ~10^13 Mo halos, and by z=0 they mostly end up as brightest group galaxies in ~10^14 Mo halos. We predict that higher-M* SMGs should have on average lower specific SFRs, less disturbed morphologies, and higher clustering. We also predict that deeper far-IR surveys will smoothly join SMGs onto the massive end of the SFR-M* relationship defined by lower-mass z=2 galaxies. Overall, our simulated rapid star-formers provide as good a match to available SMG data as merger-based scenarios, offering an alternative scenario that emerges naturally from cosmological simulations.

We analyse cosmological hydrodynamic simulations that include observationally-constrained prescriptions for galactic outflows. If these simulated winds accurately represent winds in the real Universe, then material previously ejected in winds provides the dominant source of gas infall for new star formation at redshifts z<1. This recycled wind accretion, or wind mode, provides a third physically distinct accretion channel in addition to the "hot" and "cold" modes emphasised in recent theoretical studies. Because of the interaction between outflows and gas in and around halos, the recycling timescale of wind material (t_rec) is shorter in higher-mass systems, which reside in denser gaseous environments. In these simulations, this differential recycling plays a central role in shaping the present-day galaxy stellar mass function (GSMF). If we remove all particles that were ever ejected in a wind, then the predicted GSMFs are much steeper than observed; galaxy masses are suppressed both by the direct removal of gas and by the hydrodynamic heating of their surroundings, which reduces subsequent infall. With wind recycling included, the simulation that incorporates our favoured momentum-driven wind scalings reproduces the observed GSMF for stellar masses 10^9 < M < 5x10^10 Msolar. At higher masses, wind recycling leads to excessive galaxy masses and excessive star formation rates relative to observations. In these massive systems, some quenching mechanism must suppress the re-accretion of gas ejected from star-forming galaxies. In short, as has long been anticipated, the form of the GSMF is governed by outflows; the unexpected twist here for our simulated winds is that it is not primarily the ejection of material but how the ejected material is re-accreted that governs the GSMF.

We use 2167 Sloan Digital Sky Survey (SDSS) quasar spectra to search for low-density oxygen in the Intergalactic Medium (IGM). Oxygen absorption is detected on a pixel-by-pixel basis by its correlation with Lyman-alpha forest absorption. We have developed a novel Locally Calibrated Pixel (LCP) search method that uses adjacent regions of the spectrum to calibrate interlopers and spectral artifacts, which would otherwise limit the measurement of OVI absorption. Despite the challenges presented by searching for weak OVI within the Lyman-alpha forest in spectra of moderate resolution and signal-to-noise, we find a highly significant detection of absorption by oxygen at 2.7 < z < 3.2 (the null hypothesis has a chi^2=80 for 9 data points). We interpret our results using synthetic spectra generated from a lognormal density field assuming a mixed quasar-galaxy photoionizing background (Haardt & Madau 2001) and that it dominates the ionization fraction of detected OVI. The LCP search data can be fit by a constant metallicity model with [O/H] = -2.15_(-0.09)^(+0.07), but also by models in which low-density regions are unenriched and higher density regions have a higher metallicity. The density-dependent enrichment model by Aguirre et al. (2008) is also an acceptable fit. All our successful models have similar mass-weighted oxygen abundance, corresponding to [<O/H>_MW] = -2.45+-0.06. This result can be used to find the cosmic oxygen density in the Lyman-alpha forest, Omega_(Oxy, IGM) = 1.4(+-0.2)x10^(-6) = 3x10^(-4) Omega_b. This is the tightest constraint on the mass-weighted mean oxygen abundance and the cosmic oxygen density in the Lyman-alpha forest to date and indicates that it contains approximately 16% of metals produced by star formation (Bouch\'e et al. 2008) up to z = 3.

Through the combination of high-order Adaptive Optics and coronagraphy, we report the discovery of a faint stellar companion to the A3V star zeta Virginis. This companion is ~7 magnitudes fainter than its host star in the H-band, and infrared imaging spanning 4.75 years over five epochs indicates this companion has common proper motion with its host star. Using evolutionary models, we estimate its mass to be 0.168+/-.016 solar masses, giving a mass ratio for this system q = 0.082. Assuming the two objects are coeval, this mass suggests a M4V-M7V spectral type for the companion, which is confirmed through integral field spectroscopic measurements. We see clear evidence for orbital motion from this companion and are able to constrain the semi-major axis to be greater than 24.9 AU, the period > 124$ yrs, and eccentricity > 0.16. Multiplicity studies of higher mass stars are relatively rare, and binary companions such as this one at the extreme low end of the mass ratio distribution are useful additions to surveys incomplete at such a low mass ratio. Moreover, the frequency of binary companions can help to discriminate between binary formation scenarios that predict an abundance of low-mass companions forming from the early fragmentation of a massive circumstellar disk. A system such as this may provide insight into the anomalous X-ray emission from A stars, hypothesized to be from unseen late-type stellar companions. Indeed, we calculate that the presence of this M-dwarf companion easily accounts for the X-ray emission from this star detected by ROSAT.

High contrast imaging can find and characterize gas giant planets around nearby young stars and the closest M stars, complementing radial velocity and astrometric searches by exploring orbital separations inaccessible to indirect methods. Ground-based coronagraphs are already probing within 25 AU of nearby young stars to find objects as small as ~ 3 Jupiter masses. This paper compares near-term and future ground-based capabilities with high contrast imaging modes of the James Webb Space Telescope (JWST). Monte Carlo modeling reveals that JWST can detect planets with masses as small as 0.2 MJup across a broad range of orbital separations. We present new calculations for planet brightness as a function of mass and age for specific JWST filters and extending to 0.1 MJup.

The A5V star Alcor has an M3-M4 dwarf companion, as evidenced by a novel astrometric technique. Imaging spectroscopy combined with adaptive optics coronagraphy allowed for the detection and spectrophotometric characterization of the point source at a contrast of ~6 J- and H-band magnitudes and separation of 1" from the primary star. The use of an astrometric pupil plane grid allowed us to determine the projected separations between the companion and the coronagraphically occulted primary star to <=3 milliarcsecond precision at two observation epochs. Our measurements demonstrate common parallactic and proper motion over the course of 103 days, significantly shorter than the period of time needed for most companion confirmations through proper motion measurements alone. This common parallax method is potentially more rigorous than common proper motion, ensuring that the neighboring bodies lie at the same distance, rather than relying on the statistical improbability that two objects in close proximity to each other on the sky move in the same direction. The discovery of a low-mass (~0.25M_sun) companion around a bright (V = 4.0), nearby (d = 25 pc) star highlights a region of binary star parameter space that to date has not been fully probed.

The Terrestrial Planet Finder Coronagraph (TPF-C) mission presented here is an existence proof for a flagship-class internal coronagraph space mission capable of detecting and characterizing Earth-like planets and planetary systems at visible wavelengths around nearby stars, using an existing launch vehicle. TPF-C will use spectroscopy to measure key properties of exoplanets including the presence of atmospheric water or oxygen, powerful signatures in the search for habitable worlds.

Over 300 extrasolar planets (exoplanets) have been detected orbiting nearby stars. We now hope to conduct a census of all planets around nearby stars and to characterize their atmospheres and surfaces with spectroscopy. Rocky planets within their star's habitable zones have the highest priority, as these have the potential to harbor life. Our science goal is to find and characterize all nearby exoplanets; this requires that we measure the mass, orbit, and spectroscopic signature of each one at visible and infrared wavelengths. The techniques for doing this are at hand today. Within the decade we could answer long-standing questions about the evolution and nature of other planetary systems, and we could search for clues as to whether life exists elsewhere in our galactic neighborhood.

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.

We demonstrate the versatility of a dual imaging polarimeter working in tandem with a Lyot coronagraph and Adaptive Optics to suppress the highly static speckle noise pattern--the greatest hindrance to ground-based direct imaging of planets and disks around nearby stars. Using a double difference technique with the polarimetric data, we quantify the level of speckle suppression, and hence improved sensitivity, by placing an ensemble of artificial faint companions into real data, with given total brightness and polarization. For highly polarized sources within 0.5 arcsec, we show that we achieve 3 to 4 magnitudes greater sensitivity through polarimetric speckle suppression than simply using a coronagraph coupled to a high-order Adaptive Optics system. Using such a polarimeter with a classical Lyot coronagraph at the 3.63m AEOS telescope, we have obtained a 6.5 sigma detection in the H-band of the 76 AU diameter circumstellar debris disk around the star HR 4796A. Our data represent the first definitive, ground-based, near-IR polarimetric image of the HR 4796A debris disk and clearly show the two outer ansae of the disk, evident in Hubble Space Telescope NICMOS/STIS imaging. We derive a lower limit to the fractional linear polarization of 29% caused by dust grains in the disk. In addition, we fit simple morphological models of optically thin disks to our data allowing us to constrain the dust disk scale height to 2.5{+5.0}_{-1.3} AU and scattering asymmetry parameter (g=0.20^{+.07}_{-.10}). These values are consistent with several lines of evidence suggesting that the HR 4796A disk is dominated by a micron-sized dust population, and are indeed typical of disks in transition between those surrounding the Herbig Ae stars to those associated with Vega-like stars.

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)

High-contrast observations in optical and infrared astronomy are defined as any observation requiring a technique to reveal a celestial object of interest that is in such close angular proximity to another source brighter by a factor of at least 10^5 that optical effects hinder or prevent the collection of photons directly from the target of observation. This is a relatively new type of observation that enables research on previously obscured parts of the Universe. In particular, it is most applicable to Comparative Planetary Science, a field that directly attacks such questions as "how common are planetary systems? What types of planets exist, and are there planets other than Earth that are capable of supporting life as we know it?" We survey the scientific motivations for high-contrast observations, provide an overview of the techniques currently being used or developed, and discuss some ideas and studies for future prospects.

This White Paper, submitted to the National Academy of Sciences' Astro2010 Decadal Review Committee, focuses on 2 central themes in the study of young brown dwarfs -- their formation mechanism and disk characteristics -- which are of direct relevance to fundamental questions of stellar and planetary origins and properties.

IGM metal absorption lines observed in z>6 spectra offer the opportunity to probe early feedback processes, the nature of enriching sources, and the topology of reionization. We run high-resolution cosmological simulations including galactic outflows to study the observability and physical properties of 5 ions (C II, C IV, O I, Si II, Si IV) in absorption between z=8->5. We apply three cases for ionization conditions: Fully neutral, fully reionized, and a patchy model based on the flux from the nearest galaxy. We find that our simulations broadly fit available z~5-6 IGM metal-line data, although all observations cannot be accommodated with a single ionization condition. Variations in O I absorbers among sight lines seen by Becker et al. (2006) suggest significant neutral IGM patches down to z~6. Strong C IV absorbers at z~6 may be the result of ionization by their parent galaxy. Our outflows have typical speeds of ~200 km/s and mass loading factors of ~6. Such high mass loading is critical for enriching the IGM to the observed levels while curtailing star formation to match the observed z~6 rest-frame UV luminosity function. The volume filling factor of metals increases during this epoch, but only reaches ~1% for Z>10^(-3) Zsolar by z=5. C IV is an ideal tracer of IGM metals at z~5-6, with dropping global ionization fractions to either higher or lower redshifts. This results in a strongly increasing global Omega(C IV) from z=8->5, in contrast to its relative constancy from z=5->2. Our simulations do not support widespread early IGM enrichment from e.g. Pop III stars. High-z absorbers arise from metals on their first outward journey from galaxies, at distances less than 50 kpc. The galaxies responsible for early IGM enrichment have typical M*=10^(7.0-8.5) Msolar.

The O VI ion observed in quasar absorption line spectra is the most accessible tracer of the cosmic metal distribution in the low redshift (z<0.5) intergalactic medium (IGM). We explore the nature and origin of O VI absorbers using cosmological hydrodynamic simulations including galactic outflows. We consider the effects of ionization background variations, non-equilibrium ionization and cooling, uniform metallicity, and small-scale (sub-resolution) turbulence. Our main results are 1) IGM O VI is predominantly photo-ionized with T= 10^(4.2+/-0.2) K. A key reason for this is that O VI absorbers preferentially trace over-enriched regions of the IGM at a given density, which enhances metal-line cooling such that absorbers can cool within a Hubble time. As such, O VI is not a good tracer of the WHIM. 2) The predicted O VI properties fit observables only if sub-resolution turbulence is added. The required turbulence increases with O VI absorber strength such that stronger absorbers arise from more recent outflows with turbulence dissipating on the order of a Hubble time. The amount of turbulence is consistent with other examples of turbulence observed in the IGM and galactic halos. 3) Metals traced by O VI and H I do not trace exactly the same baryons, but reside in the same large-scale structure. Observed alignment statistics are reproduced in our simulations. 4) Photo-ionized O VI traces gas in a variety of environments, and is not directly associated with the nearest galaxy, though is typically nearest to ~0.1L* galaxies. Weaker O VI components trace some of the oldest cosmic metals. 5) Very strong absorbers are more likely to be collisionally ionized, tracing more recent enrichment (<2 Gyr) within or near galactic halos.

This White Paper to the National Academy of Sciences Astro2010 Decadal Review Committee highlights cross-disciplinary science opportunities over the next decade with cold brown dwarfs, sources defined here as having photospheric temperatures less than ~1000 K.

Despite the revolution in our knowledge resulting from the detection of planets around mature stars, we know almost nothing about planets orbiting young stars because rapid rotation and active photospheres preclude detection by radial velocities or transits and because direct imaging has barely penetrated the requisite range of high contrast and angular resolution. Of the techniques presently under consideration for the coming decade, only space-based astrometry offers the prospect of discovering gas giants (100 to >> 300 Mearth), lower mass systems such as icy giants (10 to 100 Mearth), and even a few rocky, super-Earths 300 Mearth) orbiting stars ranging in age from 1 to 100 Myr. Astrometry will complement high contrast imaging which should be able to detect gas giants (1~10 MJup) in orbits from a few to a few hundred AU. An astrometric survey in combination with imaging data for a subsample of objects will allow a detailed physical understanding of the formation and evolution of young gas giant planets impossible to achieve by any one technique.

We examine metal and entropy content in galaxy groups having T_X~0.5-2 keV in cosmological hydrodynamic simulations. Our simulations include a well-constrained prescription for galactic outflows following momentum-driven wind scalings, and a sophisticated chemical evolution model. Our simulation with no outflows reproduces observed iron abundances in X-ray emitting gas, but the oxygen abundance is too low; including outflows yields iron and oxygen abundances in good agreement with data. X-ray measures of [O/Fe] primarily reflect metal distribution mechanisms into hot gas, not the ratio of Type Ia to Type II supernovae within the group. Iron abundance increases by x2 from z=1-0 independent of group size, consistent with that seen in clusters, while [O/Fe] drops by ~30%. Core entropy versus temperature is elevated over self-similar predictions regardless of outflows due to radiative cooling removing low-entropy gas, but outflows provide an additional entropy contribution below 1 keV. This results in a noticeable break in the L_X-T_X relation below 1 keV, as observed. Importantly, outflows serve to reduce the stellar content of groups to observed levels. Radial profiles from simulations are in broad agreement with observations, but there remain non-trivial discrepancies that may reflect an excess of late-time star formation in central group galaxies in our simulations. Our model with outflows suggests a connection between physical processes of galaxy formation and both pre-heating and enrichment in intragroup gas, though more definitive conclusions must await a model that simultaneously suppresses cooling flows as observed.

Using Gadget-2 cosmological hydrodynamic simulations including an observationally-constrained model for galactic outflows, we investigate how feedback from star formation distributes mass, metals, and energy on cosmic scales from z=6->0. We include instantaneous enrichment from Type II SNe, delayed enrichment from Type Ia SNe and stellar (AGB) mass loss, and we individually track C, O, Si, and Fe. Following on the successes of the momentum-driven wind scalings, we improve our implementation with an on-the-fly galaxy finder to derive wind properties based on host galaxy masses. By tracking wind particles in a suite of simulations, we find: (1) Wind material reaccretes onto a galaxy on a recycling timescale that varies inversely with galaxy mass. Hence metals driven into the IGM by galactic superwinds cannot be assumed to leave their galaxy forever. Wind material is typically recycled several times; the median number of ejections for a given wind particle is 3, so by z=0 the total mass ejected in winds exceeds 0.5\Omega_b. (2) The physical distance winds travel is fairly independent of redshift and galaxy mass (60-100 physical kpc). For sizable galaxies at later epochs, winds typically do not escape the galaxy halo, and rain back down in a halo fountain. High-z galaxies enrich a significantly larger comoving volume of the IGM, with metals migrating back into galaxies to lower z. (3) The energy imparted into winds scales with M_{gal}^{1/3}, and energy from another source besides SNe (such as photons from young stars) may be required to distributed cosmic metals as observed. (4) The production of all 4 metals tracked is globally dominated by Type II SNe at all epochs. However, intracluster gas iron content triples as a result of non-Type II sources, and the low-z IGM carbon content is boosted by AGB feedback.

The young star AB Aurigae is surrounded by a complex combination of gas-rich and dust dominated structures. The inner disk which has not been studied previously at sufficient resolution and imaging dynamic range seems to contain very little gas inside a radius of least 130 astronomical units (AU) from the star. Using adaptive-optics coronagraphy and polarimetry we have imaged the dust in an annulus between 43 and 302 AU from the star, a region never seen before. An azimuthal gap in an annulus of dust at a radius of 102 AU, along with a clearing at closer radii inside this annulus, suggests the formation of at least one small body at an orbital distance of about 100 AU. This structure seems consistent with crude models of mean motion resonances, or accumulation of material at two of the Lagrange points relative to the putative object and the star. We also report a low significance detection of a point source in this outer annulus of dust. This source may be an overdensity in the disk due to dust accreting onto an unseen companion. An alternate interpretation suggests that the object's mass is between 5 and 37 times the mass of Jupiter. The results have implications for circumstellar disk dynamics and planet formation.

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 present the results of the Gemini Deep Planet Survey, a near-infrared adaptive optics search for giant planets and brown dwarfs around nearby young stars. The observations were obtained with the Altair adaptive optics system at the Gemini North telescope and angular differential imaging was used to suppress the speckle noise of the central star. Detection limits for the 85 stars observed are presented, along with a list of all faint point sources detected around them. Typically, the observations are sensitive to angular separations beyond 0.5" with 5-sigma contrast sensitivities in magnitude difference at 1.6 micron of 9.5 at 0.5", 12.9 at 1", 15.0 at 2", and 16.5 at 5". For the typical target of the survey, a 100 Myr old K0 star located 22 pc from the Sun, the observations are sensitive enough to detect planets more massive than 2 Mjup with a projected separation in the range 40-200 AU. Second epoch observations of 48 stars with candidates (out of 54) have confirmed that all candidates are unrelated background stars. A detailed statistical analysis of the survey results, yielding upper limits on the fractions of stars with giant planet or low mass brown dwarf companions, is presented. Assuming a planet mass distribution dN/dm m^{-1.2} and a semi-major axis distribution dN/da a^{-1}, the 95% credible upper limits on the fraction of stars with at least one planet of mass 0.5-13 Mjup are 0.28 for the range 10-25 AU, 0.13 for 25-50 AU, and 0.093 for 50-250 AU. The 95% credible interval for the fraction of stars with at least one brown dwarf companion having a semi-major axis in the range 25-250 AU is 0.019 (-0.015/+0.083), irrespective of any assumption on the mass and semi-major axis distributions. The stars HD 14802, HD 166181, and HD 213845 have been resolved into binaries for the first time.

The Gemini Planet (GPI) imager is an "extreme" adaptive optics system being designed and built for the Gemini Observatory. GPI combines precise and accurate wavefront control, diffraction suppression, and a speckle-suppressing science camera with integral field and polarimetry capabilities. GPI's primary science goal is the direct detection and characterization of young, Jovian-mass exoplanets. For systems younger than 2 Gyr exoplanets more massive than 6 MJ and semimajor axes beyond 10 AU are detected with completeness greater than 50%. GPI will also discover faint debris disks, explore icy moons and minor planets in the solar system, reveal high dynamic range main-sequence binaries, and study mass loss from evolved stars. This white paper explains the role of GPI in exoplanet discovery and characterization and summarizes our recommendations to the NSF-NASA-DOE Astronomy and Astrophysics Advisory Committee ExoPlanet Task Force.

The status of 38 halo white dwarf candidates identified by Oppenheimer et al. (2001) has been intensively discussed by various authors. In analyses undertaken to date, trigonometric parallaxes are crucial missing data. Distance measurements are mandatory to kinematically segregate halo object from disk objects and hence enable a more reliable estimate of the local density of halo dark matter residing in such objects. We present trigonometric parallax measurements for 15 candidate halo white dwarfs (WDs) selected from the Oppenheimer et al. (2001) list. We observed the stars using the ESO 1.56-m Danish Telescope and ESO 2.2-m telescope from August 2001 to July 2004. Parallaxes with accuracies of 1--2 mas were determined yielding relative errors on distances of $\sim5$% for 6 objects, $\sim12$% for 3 objects, and $\sim20$% for two more objects. Four stars appear to be too distant (probably farther than 100 pc) to have measurable parallaxes in our observations. Distances, absolute magnitudes and revised space velocities were derived for the 15 halo WDs from the Oppenheimer et al. (2001) list. Halo membership is confirmed unambiguously for 6 objects while 5 objects may be thick disk members and 4 objects are too distant to draw any conclusion based solely on kinematics. Comparing our trigonometric parallaxes with photometric parallaxes used in previous work reveals an overestimation of distance as derived from photometric techniques. This new data set can be used to revise the halo white dwarf space density, and that analysis will be presented in a subsequent publication.

We conduct a detailed comparison of broad-band spectral energy distributions of six z >= 5.5 galaxies against galaxies drawn from cosmological hydrodynamic simulations. We employ a new tool called SPOC, which constrains the physical properties of observed galaxies through a Bayesian likelihood comparison with model galaxies. For five out of six observed z>=5.5 objects, our simulated galaxies match the observations at least as well as simple star formation histories such as tau-models, with similar favored values obtained for the intrinsic physical parameters such as stellar mass and star formation rate, but with substantially smaller uncertainties. Our results are broadly insensitive to simulation choices for galactic outflows and dust reddening. Hence the existence of early galaxies as observed is broadly consistent with current hierarchical structure formation models. However, one of the six objects has photometry that is best fit by a bursty SFH unlike anything produced in our simulations, driven primarily by a high K-band flux. These findings illustrate how SPOC provides a robust tool for optimally utilizing hydrodynamic simulations (or any model that predicts galaxy SFHs) to constrain the physical properties of individual galaxies having only photometric data, as well as identify objects that challenge current models. (abridged)

We use cosmological hydrodynamic simulations including galactic feedback based on observations of local starbursts to find a self-consistent evolutionary model capable of fitting the observations of the intergalactic metallicity history as traced by C IV between z=6.0->1.5. Our main finding is that despite the relative invariance in the measurement of Omega(C IV) as well as the column density and linewidth distributions over this range, continual feedback from star formation-driven winds are able to reproduce the observations, while an early enrichment scenario where a majority of the metals are injected into the IGM at z>6 is disfavored. The constancy of the C IV observations results from a rising IGM metallicity content balanced by a declining C IV ionization fraction due to a 1) decreasing physical densities, 2) increasing ionization background strength, and 3) metals becoming more shock-heated at lower redshift. Our models predict that ~20x more metals are injected into the IGM between z=6->2 than at z>6. We show that the median C IV absorber at z=2 traces metals injected 1 Gyr earlier indicating that the typical metals traced by C IV are neither from very early times nor from very recent feedback.

We present predictions for the cosmic metal budget in various phases of baryons from redshift z=6-0, taken from a cosmological hydrodynamic simulation that includes a well-constrained model for enriched galactic outflows. We find that substantial amounts of metals are found in every baryonic phase at all epochs, with diffuse intergalactic gas dominating the metal budget at early epochs and stars and halo gas dominating at recent epochs. We provide a full accounting of metals in the context of the missing metals problem at z~2.5, showing that ~40% of the metals are in galaxies, and the remainder is divided between diffuse IGM gas and shocked gas in halos and filamentary structures. Comparisons with available observations of metallicity and metal mass fraction evolution show broad agreement. We predict stars have a mean metallicity of one-tenth solar already at z=6, which increases slowly to one-half solar today, while stars just forming today have typically solar metallicity. Our HI column density-weighted mean metallicity (comparable to Damped Ly-alpha system metallicities) slowly increases from one-tenth to one-third solar from z=6-1, then falls to one-quarter solar at z=0. The global mean metallicity of the universe tracks ~50% higher than that of the diffuse phase down to z~1, and by z=0 it has a value around one-tenth solar. Metals move towards higher densities and temperatures with time, peaking around the mean cosmic density at z=2 and an overdensity of 100 at z=0. We study how carbon and oxygen ions trace the path of metals in phase space, and show that OIII-OVII lines provide the most practical option for constraining intergalactic medium metals at z<2.

The major obstacle to the direct detection of companions to nearby stars is the overwhelming brightness of the host star. Current instruments employing the combination of adaptive optics (AO) and coronagraphy can typically detect objects within 2'' of the star that are 10^{4-5} times fainter. Correlated speckle noise is one of the biggest obstacles limiting such high-contrast imaging. We have obtained a series of 284 8 s, AO-corrected, coronagraphically occulted H-band images of the star Vega at the 3.63 m AEOS telescope located on Haleakala, Hawaii. This dataset is unique for studying the temporal behavior of speckle noise, and represents the first time such a study on highly corrected _coronagraphic_ AO images has been carried out in a quantitative way. We find the speckle pattern to be highly stable in both position and time in our data. This is due to the fact that the AO system corrects disturbances to the stellar wave front at the level where the instrumental wave front errors dominate. Because of this, we find that our detection limit is not significantly improved simply with increased exposure time alone. However, we are able to improve our dynamic range by 1.5-2 magnitudes through subtraction of static/quasi-static speckles in two rotating frames: the telescope pupil frame and the deformable mirror frame. Furthermore, from our data, we are able to constrain the mass of any purported companion to Vega to be less than 45 M_J at 8 AU and less than 30 M_J at 16 AU, radii not previously probed at these sensitivities.

We investigate models of self-consistent chemical enrichment of the intergalactic medium (IGM) from z=6 to 1.5, based on hydrodynamic simulations of structure formation that incorporate galactic outflows. Our main result is that outflow parameterizations based on momentum-driven winds as seen in local starburst galaxies provide the best agreement with observations of CIV absorption at z~2-5. Such models sufficiently enrich the high-z IGM to produce a global mass density of CIV absorbers that is relatively invariant from z=5.5 to 1.5, in agreement with observations. This occurs despite an increase in the volume-averaged metallicity by x5-10 over this redshift range, because energy input from outflows causes a drop in the global ionization fraction of CIV. Comparisons to observed CIV column density and linewidth distributions and CIV-based pixel optical depth ratios provide significant constraints on wind models. Our best-fitting models show mean IGM temperatures only slightly above our no-outflow case, metal filling factors of just a few % with volume-weighted metallicities around 0.001 at z~3, significant amounts of collisionally-ionized CIV absorption, and a metallicity-density relationship that rises rapidly at low overdensities and flattens at higher ones. In general, we find that outflow speeds must be high enough to enrich the low-density IGM at early times but low enough not to overheat it, and concurrently must significantly suppress early star formation while still producing enough early metals. It is therefore non-trivial that locally-calibrated momentum-driven wind scenarios naturally yield the desired strength and evolution of outflows, and suggest that such models represent a significant step towards understanding the impact of galactic outflows on galaxies and the IGM across cosmic time.

We have analyzed HST/NICMOS2 F110W, F160W, F165M, and F207M band images covering the central 1'x1' of the cluster associated with Mon R2 in order to constrain the Initial Mass Function (IMF) down to 20 Mjup. The flux ratio between the F165M and F160W bands was used to measure the strength of the water band absorption feature and select a sample of 12 out of the total sample of 181 objects that have effective temperatures between 2700 K and 3300 K. These objects are placed in the HR diagram together with sources observed by Carpenter et al. (1997) to estimate an age of ~1 Myr for the low mass cluster population. By constructing extinction limited samples, we are able to constrain the IMF and the fraction of stars with a circumstellar disk in a sample that is 90% complete for both high and low mass objects. For stars with estimated masses between 0.1 Msun and 1.0 Msun for a 1 Myr population with Av < 19 mag, we find that 27+-9% have a near-infrared excess indicative of a circumstellar disk. The derived fraction is similar to, or slightly lower than, the fraction found in other star forming regions of comparable age. We constrain the number of stars in the mass interval 0.08-1.0Msun to the number of objects in the mass interval 0.02-0.08 Msun by forming the ratio, R**=N(0.08-1Msun)/N(0.02-0.08Msun) for objects in an extinction limited sample complete for Av < 7 mag. The ratio is found to be R^**=2.2+-1.3 assuming an age of 1 Myr, consistent with the similar ratio predicted by the system IMF proposed by Chabrier (2003). The ratio is similar to the ratios observed towards the Orion Nebula Cluster and IC 348 as well as the ratio derived in the 28 square degree survey of Taurus by Guieu et al. (2006).

We use cosmological hydrodynamic simulations with enriched galactic outflows to compare predictions for the galaxy mass-metallicity (M*-Z) relation versus observations at z~2 from Erb et al. (2006). With no outflows, galaxies are over-enriched, indicating that outflows are required not only to suppress star formation and enrich the IGM but also to lower galaxy metal content. The observed M*-Z slope is matched both in our model without winds as well as our favored outflow model where the outflow velocity scales as the escape velocity, but is too steep in a model with constant outflow speeds. If outflows are too widespread at early times, the IGM out of which smaller galaxies form can become pre-polluted, resulting in a low-mass flattening of the M*-Z relation that is inconsistent with data. Remarkably, the same momentum-driven wind model that provides the best agreement with IGM enrichment data also yields the best agreement with the z~2 M*-Z relation, showing the proper outflow scaling and strength to match the observed slope and amplitude. In this model, the M*-Z relation evolves slowly from z=6-2; an (admittedly uncertain) extrapolation to z=0 broadly matches local M*-Z observations. Overall, the M*-Z relation provides critical constraints on galactic outflow processes during the heydey of star formation in the Universe.

We propose a solution to the problem of astrometric and photometric calibration of coronagraphic images with a simple optical device which, in theory, is easy to use. Our design uses the Fraunhofer approximation of Fourier optics. Placing a periodic grid of wires (we use a square grid) with known width and spacing in a pupil plane in front of the occulting coronagraphic focal plane mask produces fiducial images of the obscured star at known locations relative to the star. We also derive the intensity of these fiducial images in the coronagraphic image. These calibrator images can be used for precise relative astrometry, to establish companionship of other objects in the field of view through measurement of common proper motion or common parallax, to determine orbits, and to observe disk structure around the star quantitatively. The calibrator spots also have known brightness, selectable by the coronagraph designer, permitting accurate relative photometry in the coronagraphic image. This technique, which enables precision exoplanetary science, is relevant to future coronagraphic instruments, and is particularly useful for `extreme' adaptive optics and space-based coronagraphy.

We present predictions drawn from cosmological hydrodynamic simulations for the physical, photometric and emission line properties of galaxies present during the latter stages of reionization from z=9-6. We find significant numbers of galaxies that have stellar masses exceeding 10^8 Mo during this epoch, with metallicities exceeding one-thirtieth solar. Far from primeval "first-star" objects,these objects are likely to have reionized their infall regions prior to z=9, are dominated by atomic rather than molecular cooling, and are not expected to be forming Population III stars. Galaxies exhibit a slowly evolving comoving autocorrelation length from z=9-6, continuing a trend seen at lower redshifts in which the rapidly dropping bias counteracts the rapidly increasing matter clustering. These sources can be marginally detected using current instruments, but modest increases in sensitivity or survey area would yield significantly increased samples. We compare to current observations of the z~6 rest-UV and Ly-alpha line luminosity functions, and find good agreement. We also compare with the z~7 object observed by Egami et al., and find that such systems are ubiquitous in our simulations. The intrinsic Ly-alpha luminosity function evolves slowly from z=9-6, meaning that it should also be possible to detect these objects with upcoming narrow band surveys such as DAzLE. We make predictions for near-IR surveys with JWST, and show that while a high density of sources will be found, Population III objects may remain elusive. We present and compare simulations with several recipes for superwind feedback, and show that while our broad conclusions are insensitive to this choice, a feedback model based on momentum-driven winds is favored in comparisons with available data. (abridged)

We study the way Lyot coronagraphs with unapodized entrance pupils respond to small, low order phase aberrations. This study is applicable to ground-based adaptive optics coronagraphs operating at 90% and higher Strehl ratios, as well as to some space-based coronagraphs with intrinsically higher Strehl ratio imaging. We utilize a second order expansion of the monochromatic point-spread function (written as a power spectrum of a power series in the phase aberration over clear aperture) to derive analytical expressions for the response of a `band-limited' Lyot coronagraph (BLC) to small, low order, phase aberrations. The BLC possesses a focal plane mask with an occulting spot whose opacity profile is a spatially band-limited function rather than a hard-edged, opaque disk. The BLC is, to first order, insensitive to tilt and astigmatism. Undersizing the stop in the re-imaged pupil plane (the Lyot plane) following the focal plane mask can alleviate second order effects of astigmatism, at the expense of system throughput and angular resolution. The optimal degree of such undersizing depends on individual instrument designs and goals. Our analytical work engenders physical insight, and complements existing numerical work on this subject. Our methods can be extended to treat the passage of higher order aberrations through band-limited Lyot coronagraphs, by using our polynomial decomposition or an analogous Fourier approach.

Any predictive theory of star formation must explain observed variations (or lack thereof) in the initial mass function. Recent work suggests that we might expect quantitative variations in the IMF as a function of metallicity (Larson 2005) or magnetic field strength (Shu et al. 2004). We summarize results from several on-going studies attempting to constrain the ratio of high to low mass stars, as well as stars to sub- stellar objects, in a variety of different environments, all containing high mass stars. First, we examine the ratio of stars to sub--stellar objects in the nearby Mon R2 region utilizing NICMOS/HST data. We compare our results to the IMF by Kroupa (2002)]} and to the observed ratios for IC 348 and Orion. Second, we present preliminary results for the ratio of high to low mass stars in W51, the most luminous HII region in the galaxy. Based on ground--based multi--colour images of the cluster obtained with the MMT adaptive optics system, we derive a lower limit to the ratio of high-mass to low-mass stars and compare it to the ratios for nearby clusters. Finally, we present the derived IMF for the R136 region in the LMC where the metallicity is 1/4 solar using HST/NICMOS data. We find that the IMF is consistent with that characterizing the field (Chabrier 2003), as well as nearby star--forming regions, down to 1.0 Msun outside 2 pc. Whereas the results for both Mon R2 and R136 are consistent with the nearby clusters, the ratio of high to low mass stars in W51 tentatively indicates a lack of low--mass objects.

Adaptive optics (AO) systems have significantly improved astronomical imaging capabilities over the last decade, and are revolutionizing the kinds of science possible with 4-5m class ground-based telescopes. A thorough understanding of AO system performance at the telescope can enable new frontiers of science as observations push AO systems to their performance limits. We look at recent advances with wave front reconstruction (WFR) on the Advanced Electro-Optical System (AEOS) 3.6 m telescope to show how progress made in improving WFR can be measured directly in improved science images. We describe how a "waffle mode" wave front error (which is not sensed by a Fried geometry Shack-Hartmann wave front sensor) affects the AO point-spread function (PSF). We model details of AEOS AO to simulate a PSF which matches the actual AO PSF in the I-band, and show that while the older observed AEOS PSF contained several times more waffle error than expected, improved WFR techniques noticeably improve AEOS AO performance. We estimate the impact of these improved WFRs on H-band imaging at AEOS, chosen based on the optimization of the Lyot Project near-infrared coronagraph at this bandpass.

We present high resolution (~0.1"), very high Strehl ratio (0.97+-0.03) mid-infrared (IR) adaptive optics (AO) images of the AGB star RV Boo utilizing the MMT adaptive secondary AO system. RV Boo was observed at a number of wavelengths over two epochs (9.8 um in May 2003, 8.8, 9.8 and 11.7 um in February 2004) and appeared slightly extended at all wavelengths. While the extension is very slight at 8.8 and 11.7 um data, the extension is somewhat more pronounced at 9.8 um. With such high Strehls we can achieve super-resolutions of 0.1" by deconvolving RV Boo with a point-spread function (PSF) derived from an unresolved star. We tentatively resolve RV Boo into a 0.16" FWHM extension at a position angle of 120 degrees. At a distance of 390(+250)(-100) pc, this corresponds to a FWHM of 60(+40)(-15) AU. We measure a total flux at 9.8 um of 145+-24 Jy for the disk and star. Based on a dust thermal emission model for the observed IR spectral energy distribution and the 9.8 um AO image, we derive a disk dust mass of 1.6x10^-6 Msun and an inclination of 30 to 45 degrees from edge-on. We discuss whether the dust disk observed around RV Boo is an example of the early stages in the formation of asymmetric structure in planetary nebula.

In this paper we present new and improved data on 38 cool white dwarfs identified by Oppenheimer et al. 2001 (OHDHS) as candidate dark halo objects. Using the high-res spectra obtained with LRIS, we measure radial velocities for 13 WDs that show an H alpha line. We show that the knowledge of RVs decreases the UV-plane velocities by only 6%. The radial velocity sample has a W-velocity dispersion of sig_W = 59 km/s--in between the values associated with the thick disk and the stellar halo. We also see indications for the presence of two populations by analyzing the velocities in the UV plane. In addition, we present CCD photometry for half of the sample, and with it recalibrate the photographic photometry of the remaining WDs. Using the new photometry in standard bands, and by applying the appropriate color-magnitude relations for H and He atmospheres, we obtain new distance estimates. New distances of the WDs that were not originally selected as halo candidates yield 13 new candidates. On average, new distances produce velocities in the UV plane that are larger by 10%, with already fast objects gaining more. Using the new data, while applying the same UV-velocity cut (94 km/s) as in OHDHS, we find a density of cool WDs of 1.7e-4 pc^-3, confirming the value of OHDHS. In addition, we derive the density as a function of the UV-velocity cutoff. The density (corrected for losses due to higher UV cuts) starts to flatten out at 150 km/s (0.4e-4 pc^-3), and is minimized (thus minimizing a possible non-halo contamination) at 190 km/s (0.3e-4 pc^-3). These densities are in a rough agreement with the estimates for the stellar halo WDs, corresponding to a factor of 1.9 and 1.4 higher values.

The properties of metal thin films have been largely overlooked in discussions of the technical limitations and problems that arise in the field of direct detection of exoplanets. Here, polarization properties and anisotropy properties of highly reflecting thin metal films are examined within the context of the requirements for the ultra-low scattered-light system performance of coronagraphs applied to space and ground-based high-contrast, white-light astronomy. Wavelength-dependent optical constants for highly reflecting thin metal films, taken from the literature are used to calculate the polarization-dependent transmissivity of a typical coronagraph. The effects of degraded performance on the astronomical science are examined. Suggestions are made for future work.

We describe the symmetries present in the point-spread function (PSF) of an optical system either located in space or corrected by an adaptive o to Strehl ratios of about 70% and higher. We present a formalism for expanding the PSF to arbitrary order in terms of powers of the Fourier transform of the residual phase error, over an arbitrarily shaped and apodized entrance aperture. For traditional unapodized apertures at high Strehl ratios, bright speckles pinned to the bright Airy rings are part of an antisymmetric perturbation of the perfect PSF, arising from the term that is first order in the residual phase error. There are two symmetric second degree terms. One is negative at the center, and, like the first order term, is modulated by the perfect image's field strength -- it reduces to the Marechal approximation at the center of the PSF. The other is non-negative everywhere, zero at the image center, and can be responsible for an extended halo -- which limits the dynamic range of faint companion detection in the darkest portions of the image. In regimes where one or the other term dominates the speckles in an image, the symmetry of the dominant term can be exploited to reduce the effect of those speckles, potentially by an order of magnitude or more. We demonstrate the effects of both secondary obscuration and pupil apodization on the structure of residual speckles, and discuss how these symmetries can be exploited by appropriate telescope and instrument design, observing strategies, and filter bandwidths to improve the dynamic range of high dynamic range AO and space-based observations. Finally, we show that our analysis is relevant to high dynamic range coronagraphy.

We report differential I-band photometry of four early L-dwarfs obtained to study variability. We detect variability on the timescale of hours in two objects, 2M0746425+200032 (at a level of 0.007 mag -- 6.5 sigma) and 2M1108307+683017 (0.012 mag -- 5 sigma). We also place upper limits of 0.02 mag (1 sigma) on the variability of two others.

The low infrared background and high atmospheric transparency are the principal advantages of Antarctic Plateau sites for astronomy. However, the poor seeing (between one and three arcseconds) negates much of the sensitivity improvements that the Antarctic atmosphere offers, compared to mid-latitude sites such as Mauna Kea or Cerro Paranal. The seeing at mid-latitude sites, though smaller in amplitude, is dominated by turbulence at altitudes of 10 to 20 km. Over the Antarctic plateau, virtually no high altitude turbulence is present in the winter. The mean square error for an astrometric measurement with a dual-beam, differential astrometric interferometer in the very narrow angle regime is proportional to the integral of $h^2 C_N^2(h)$. Therefore, sites at which the turbulence occurs only at low altitudes offer large gains in astrometric precision. We show that a modest interferometer at the South Pole can achieve 10 $\mu$as differential astrometry 300 times faster than a comparable interferometer at a good mid-latitude site, in median conditions. Science programs that would benefit from such an instrument include planet detection and orbit determination and astrometric observation of stars microlensed by dark matter candidates.

We reply to the criticism of Gibson & Flynn and of Graff regarding "Direct detection of galactic halo dark matter" by Oppenheimer et al. (2001) (Science, Vol. 292, p. 698), which reported on the discovery of a significant population of halo white dwarfs found in our South Galactic Cap survey. None of their arguments against our claims withstand close scrutiny.

The Milky Way Galaxy contains a large, spherical component which is believed to harbor a substantial amount of unseen matter. Recent observations indirectly suggest that as much as half of this ``dark matter'' may be in the form of old, very cool white dwarfs, the remnants of an ancient population of stars as old as the Galaxy itself. We conducted a survey to find faint, cool white dwarfs with large space velocities, indicative of their membership in the Galaxy's spherical halo component. The survey reveals a substantial, directly observed population of old white dwarfs, too faint to be seen in previous surveys. This newly discovered population accounts for at least 2% of the halo dark matter. It provides a natural explanation for the indirect observations, and represents a direct detection of Galactic halo dark matter.

We present the technique and results of a survey of stars within 8 pc of the Sun with declinations $\bf\delta > -35\arcdeg$ (J2000.00). The survey, designed to find without color bias faint companions, consists of optical coronagraphic images of the 1\arcmin field of view centered on each star and infrared direct images with a 32\asc field of view. The images were obtained through the optical Gunn $r$ and $z$ filters and the infrared J and K filters. The survey achieves sensitivities up to four absolute magnitudes fainter than the prototype brown dwarf, Gliese 229B. However, this sensitivity varies with the seeing conditions, the intrinsic brightness of the star observed and the angular distance from the star. As a result we tabulate sensitivity limits for each star in the survey. We used the criterion of common proper motion to distinguish companions and to determine their luminosities. In addition to the brown dwarf Gliese 229B, we have identified 6 new stellar companions of the sample stars. Since the survey began, accurate trigonometric parallax measurements for most of the stars have become available. As a result some of the stars we originally included should no longer be included in the 8 pc sample. In addition, the 8 pc sample is incomplete at the faint end of the main sequence complicating our calculation of the binary fraction of brown dwarfs. We assess the sensitivity of the survey to stellar companions and to brown dwarf companions of various masses and ages.

We present new spectroscopic and photometric measurements of the white dwarfs LHS 3250 and WD 0346+246. Along with F351-50, these white dwarfs are the coolest ones known, all with effective temperatures below 4000 K. Their membership in the Galactic halo population is discussed, and detailed comparisons of all three objects with new atmosphere models are presented. The new models consider the effects of mixed H/He atmospheres and indicate that WD 0346+246 and F351-50 have predominantly helium atmospheres with only traces of hydrogen. LHS 3250 may be a double degenerate whose average radiative temperature is between 2000 and 4000 K, but the new models fail to explain this object.

We have measured the rest-frame B,V, and I-band light curves of a high-redshift type Ia supernova (SN Ia), SN 1999Q (z=0.46), using HST and ground-based near-infrared detectors. A goal of this study is the measurement of the color excess, E_{B-I}, which is a sensitive indicator of interstellar or intergalactic dust which could affect recent cosmological measurements from high-redshift SNe Ia. Our observations disfavor a 30% opacity of SN Ia visual light by dust as an alternative to an accelerating Universe. This statement applies to both Galactic-type dust (rejected at the 3.4 sigma confidence level) and greyer dust (grain size > 0.1 microns; rejected at the 2.3 to 2.6 sigma confidence level) as proposed by Aguirre (1999). The rest-frame $I$-band light cur ve shows the secondary maximum a month after B maximum typical of nearby SNe Ia of normal luminosi ty, providing no indication of evolution as a function of redshift out to z~0.5. A n expanded set of similar observations could improve the constraints on any contribution of extragalactic dust to the dimming of high-redshift SNe Ia.

After a discussion of the physical processes in brown dwarfs, we present a complete, precise definition of brown dwarfs and of planets inspired by the internal physics of objects between 0.1 and 0.001 M_sun. We discuss observational techniques for characterizing low-luminosity objects as brown dwarfs, including the use of the lithium test and cooling curves. A brief history of the search for brown dwarfs leads to a detailed review of known isolated brown dwarfs with emphasis on those in the Pleiades star cluster. We also discuss brown dwarf companions to nearby stars, paying particular attention to Gliese 229B, the only known cool brown dwarf.

We present a spectrum of the cool (T_eff = 900 K) brown dwarf Gliese 229B. This spectrum, with a relatively high signal-to-noise ratio per spectral resolution element (> 30), spans the wavelength range from 0.837 microns to 5.0 microns. We identify a total of four different major methane absorption features, including the fundamental band at 3.3 microns, at least four steam bands, and two neutral cesium features. We confirm the recent detection of carbon monoxide (CO) in excess of what is predicted by thermochemical equilibrium calculations. Carbon is primarily involved in a chemical balance between methane and CO at the temperatures and pressures present in the outer parts of a brown dwarf. At lower temperatures, the balance favors methane, while in the deeper, hotter regions, the reaction reverses to convert methane into CO. The presence of CO in the observable part of the atmosphere is therefore a sensitive indicator of vertical flows. The high signal-to-noise ratio in the 1 to 2.5 microns region permits us to place constraints on the quantity of dust in the atmosphere of the brown dwarf. We are unable to reconcile the observed spectrum with synthetic spectra that include the presences of dust. The presence of CO but lack of dust may be a clue to the location of the boundaries of the outer convective region of the atmosphere: The lack of dust may mean that it is not being conveyed into the photosphere by convection, or that it exists in patchy clouds. If the dust is not in clouds, but rather sits below the outer convective region, we estimate that the boundary between outer convective and inner radiative layers is between 1250 K and 1600 K, in agreement with recent models.