Search

We conducted an investigation of the Coma cluster of galaxies by running a series of constrained hydrodynamic simulations with GIZMO-SIMBA and GADGET-3, based on initial conditions reconstructed from the SDSS survey volume in the ELUCID project. We compared simulation predictions and observations for galaxies, ICM and IGM in and around the Coma cluster to constrain galaxy formation physics. Our results demonstrate that this type of constrained investigation allows us to probe in more detail the implemented physical processes, because the comparison between simulations and observations is free of cosmic variance and hence can be conducted in a ''one-to-one'' manner. We found that an increase in the earlier star formation rate and the supernova feedback of the original GIZMO-SIMBA model is needed to match observational data on stellar, ISM and ICM metallicity. The simulations without AGN feedback can well reproduce the observational ICM electron density, temperature, and entropy profiles, ICM substructures, and the IGM temperature-density relation, while the ones with AGN feedback usually fail. However, one requires something like AGN feedback to reproduce a sufficiently large population of quiescent galaxies, particularly in low-density regions. The constrained simulations of the Coma cluster thus provide a test bed to understand processes that drive galaxy formation and evolution., Comment: 21 pages, 9 figures

We define a physically-motivated measure for galactic bar length, called the dynamical length. The dynamical length of the bar corresponds to the radial extent of the orbits that are the backbone supporting the bar feature. We propose a direct observational technique using integral field unit spectroscopy to measure it. Identifying these orbits and using the dynamical length is a more faithful tracer of the secular evolution and influence of the bar. We demonstrate the success of the metric for recovering the maximal bar-parenting orbit in a range of simulations, and to show its promise we perform its measurement on a real galaxy. We also study the difference between traditionally used ellipse fit approaches to determine bar length and the dynamical length proposed here in a wide range of bar-forming N-body simulations of a stellar disc and dark matter halo. We find that ellipse fitting may severely overestimate measurements of the bar length by a factor of 1.5-2.5 relative to the extent of the orbits that are trapped and actually comprise the bar. This bias leads to overestimates of both bar mass and the ratio of corotation radius to bar length, i.e. the bar speed, affecting inferences about the evolution of bars in the real universe., Comment: Submitted to MNRAS, comments welcome. arXiv admin note: text overlap with arXiv:1903.08203

Using reconstructed initial conditions in the SDSS survey volume, we carry out constrained hydrodynamic simulations in three regions representing different types of the cosmic web: the Coma cluster of galaxies; the SDSS great wall; and a large low-density region at $z\sim 0.05$. These simulations, which include star formation and stellar feedback but no AGN formation and feedback, are used to investigate the properties and evolution of intergalactic and intra-cluster media. About half of the warm-hot intergalactic gas is associated with filaments in the local cosmic web. Gas in the outskirts of massive filaments and halos can be heated significantly by accretion shocks generated by mergers of filaments and halos, respectively, and there is a tight correlation between gas temperature and the strength of the local tidal field. The simulations also predict some discontinuities associated with shock fronts and contact edges, which can be tested using observations of the thermal SZ effect and X-rays. A large fraction of the sky is covered by Ly$\alpha$ and OVI absorption systems, and most of the OVI systems and low-column density HI systems are associated with filaments in the cosmic web. The constrained simulations, which follow the formation and heating history of the observed cosmic web, provide an important avenue to interpret observational data. With full information about the origin and location of the cosmic gas to be observed, such simulations can also be used to develop observational strategies., Comment: 21 pages, 15 figures, 1 table

The bimodal absorption system imaging campaign (BASIC) aims to characterize the galaxy environments of a sample of 36 HI-selected partial Lyman limit systems (pLLSs) and Lyman limit systems (LLSs) in 23 QSO fields at $z \lesssim 1$. These pLLSs/LLSs provide a unique sample of absorbers with unbiased and well-constrained metallicities, allowing us to explore the origins of metal-rich and low-metallicity circumgalactic medium (CGM) at $z<1$. Here we present Keck/KCWI and VLT/MUSE observations of 11 of these QSO fields (19 pLLSs) that we combine with HST/ACS imaging to identify and characterize the absorber-associated galaxies. We find 23 unique absorber-associated galaxies, with an average of one associated galaxy per absorber. For seven absorbers, all with $<10\%$ solar metallicities, we find no associated galaxies with $\log M_\star \gtrsim 9.0$ within $\rho/R_{vir}$ and $|\Delta v|/v_{esc} \le$ 1.5 with respect to the absorber. We do not find any strong correlations between the metallicities or HI column densities of the gas and most of the galaxy properties, except for the stellar mass of the galaxies: the low-metallicity ([X/H] $\le -1.4$) systems have a probability of $0.39^{+0.16}_{-0.15}$ for having a host galaxy with $\log M_\star \ge 9.0$ within $\rho/R_{vir} \le 1.5$, while the higher metallicity absorbers have a probability of $0.78^{+0.10}_{-0.13}$. This implies metal-enriched pLLSs/LLSs at $z<1$ are typically associated with the CGM of galaxies with $\log M_\star > 9.0$, whereas low-metallicity pLLSs/LLSs are found in more diverse locations, with one population arising in the CGM of galaxies and another more broadly distributed in overdense regions of the universe. Using absorbers not associated with galaxies, we estimate the unweighted geometric mean metallicity of the intergalactic medium to be [X/H] $\lesssim -2.1$ at $z<1$, which is lower than previously estimated., Comment: Key figures: 13 and 18. Modified abstract. Accepted for publication in ApJ. Updated after referee comments, expanded subsection 4.4

In the local Universe, OH megamasers (OHMs) are detected almost exclusively in infrared-luminous galaxies, with a prevalence that increases with IR luminosity, suggesting that they trace gas-rich galaxy mergers. Given the proximity of the rest frequencies of OH and the hyperfine transition of neutral atomic hydrogen (HI), radio surveys to probe the cosmic evolution of HI in galaxies also offer exciting prospects for exploiting OHMs to probe the cosmic history of gas-rich mergers. Using observations for the Looking At the Distant Universe with the MeerKAT Array (LADUMA) deep HI survey, we report the first untargeted detection of an OHM at $z > 0.5$, LADUMA J033046.20$-$275518.1 (nicknamed "Nkalakatha"). The host system, WISEA J033046.26$-$275518.3, is an infrared-luminous radio galaxy whose optical redshift $z \approx 0.52$ confirms the MeerKAT emission line detection as OH at a redshift $z_{\rm OH} = 0.5225 \pm 0.0001$ rather than HI at lower redshift. The detected spectral line has 18.4$\sigma$ peak significance, a width of $459 \pm 59\,{\rm km\,s^{-1}}$, and an integrated luminosity of $(6.31 \pm 0.18\,{\rm [statistical]}\,\pm 0.31\,{\rm [systematic]}) \times 10^3\,L_\odot$, placing it among the most luminous OHMs known. The galaxy's far-infrared luminosity $L_{\rm FIR} = (1.576 \pm 0.013) \times 10^{12}\,L_\odot$ marks it as an ultra-luminous infrared galaxy; its ratio of OH and infrared luminosities is similar to those for lower-redshift OHMs. A comparison between optical and OH redshifts offers a slight indication of an OH outflow. This detection represents the first step towards a systematic exploitation of OHMs as a tracer of galaxy growth at high redshifts., Comment: 10 pages, 4 figures. Accepted to ApJ Letters

Observations of cold molecular gas reservoirs are critical for understanding the shutdown of star formation in massive galaxies. While dust continuum is an efficient and affordable tracer, this method relies upon the assumption of a "normal" molecular-gas to dust mass ratio, $\delta_{\mathrm{GDR}}$, typically of order one hundred. Recent null detections of quiescent galaxies in deep dust continuum observations support a picture where the cold gas and dust has been rapidly depleted or expelled. In this work, we present another viable explanation: a significant fraction of galaxies with low star formation per unit stellar mass are predicted to have extreme $\delta_{\mathrm{GDR}}$ ratios. We show that simulated massive quiescent galaxies at $0 < z < 3$ in the \textsc{simba} cosmological simulations have $\delta_{\mathrm{GDR}}$ values that extend $>$4 orders of magnitude. The dust in most simulated quiescent galaxies is destroyed significantly more rapidly than the molecular gas depletes, and cannot be replenished. The transition from star-forming to quiescent halts dust formation via star formation processes, with dust subsequently destroyed by supernova shocks and thermal sputtering of dust grains embedded in hot plasma. After this point, the dust growth rate in the models is not sufficient to overcome the loss of $>$3 orders of magnitude in dust mass to return to normal values of $\delta_{\mathrm{GDR}}$ despite having high metallicity. Our results indicate that it is not straight forward to use a single observational indicator to robustly pre-select exotic versus normal ratios. These simulations make strong predictions that can be tested with millimeter facilities., Comment: 9 pages, 5 figures, version accepted for publication in the Astrophysical Journal Letters

We develop a model to establish the interconnection between galaxies and their dark matter halos. We use Principal Component Analysis (PCA) to reduce the dimensionality of both the mass assembly histories of halos/subhalos and the star formation histories of galaxies, and Gradient Boosted Decision Trees (GBDT) to transform halo/subhalo properties into galaxy properties. We use two sets of hydrodynamic simulations to motivate our model architecture and to train the transformation. We then apply the two sets of trained models to dark matter only (DMO) simulations to show that the transformation is reliable and statistically accurate. The model trained by a high-resolution hydrodynamic simulation, or by a set of such simulations implementing the same physics of galaxy formation, can thus be applied to large DMO simulations to make "mock" copies of the hydrodynamic simulation. The model is both flexible and interpretable, which paves the way for future applications in which we will constrain the model using observations at different redshifts simultaneously and explore how galaxies form and evolve in dark matter halos empirically., Comment: 21 pages, 16 figures, 2 tables, MNRAS accepted

Although galactic winds play a critical role in regulating galaxy formation, hydrodynamic cosmological simulations do not resolve the scales that govern the interaction between winds and the ambient circumgalactic medium (CGM). We implement the Physically Evolved Wind (PhEW) model of Huang et al. (2020) in the GIZMO hydrodynamics code and perform test cosmological simulations with different choices of model parameters and numerical resolution. PhEW adopts an explicit subgrid model that treats each wind particle as a collection of clouds that exchange mass, metals, and momentum with their surroundings and evaporate by conduction and hydrodynamic instabilities as calibrated on much higher resolution cloud scale simulations. In contrast to a conventional wind algorithm, we find that PhEW results are robust to numerical resolution and implementation details because the small scale interactions are defined by the model itself. Compared to conventional wind simulations with the same resolution, our PhEW simulations produce similar galaxy stellar mass functions at $z\geq 1$ but are in better agreement with low-redshift observations at $M_* < 10^{11}M_\odot$ because PhEW particles shed mass to the CGM before escaping low mass halos. PhEW radically alters the CGM metal distribution because PhEW particles disperse metals to the ambient medium as their clouds dissipate, producing a CGM metallicity distribution that is skewed but unimodal and is similar between cold and hot gas. While the temperature distributions and radial profiles of gaseous halos are similar in simulations with PhEW and conventional winds, these changes in metal distribution will affect their predicted UV/X-ray properties in absorption and emission., Comment: 23 pages, 17 figures, MNRAS accepted

We present the N-body simulation techniques in EXP. EXP uses empirically-chosen basis functions to expand the potential field of an ensemble of particles. Unlike other basis function expansions, the derived basis functions are adapted to an input mass distribution, enabling accurate expansion of highly non-spherical objects, such as galactic discs. We measure the force accuracy in three models, one based on a spherical or aspherical halo, one based on an exponential disc, and one based on a bar-based disc model. We find that EXP is as accurate as a direct-summation or tree-based calculation, and in some ways is better, while being considerably less computationally intensive. We discuss optimising the computation of the basis function representation. We also detail numerical improvements for performing orbit integrations, including timesteps., Comment: Submitted to MNRAS, comments welcome

Quasar absorption systems encode a wealth of information about the abundances, ionization structure, and physical conditions in intergalactic and circumgalactic media. Simple (often single-phase) photoionization models are frequently used to decode such data. Using five discrete absorbers from the COS Absorption Survey of Baryon Harbors (CASBaH) that exhibit a wide range of detected ions (e.g., Mg II, S II--S VI, O II--O VI, Ne VIII), we show several examples where single-phase ionization models cannot reproduce the full set of measured column densities. To explore models that can self-consistently explain the measurements and kinematic alignment of disparate ions, we develop a Bayesian multiphase ionization modeling framework that characterizes discrete phases by their unique physical conditions and also investigates variations in the shape of the UV flux field, metallicity, and relative abundances. Our models require at least two (but favor three) distinct ionization phases ranging from $T \approx 10^{4}$ K photoionized gas to warm-hot phases at $T \lesssim 10^{5.8}$ K. For some ions, an apparently single absorption "component" includes contributions from more than one phase, and up to 30% of the H I is not from the lowest ionization phase. If we assume that all of the phases are photoionized, we cannot find solutions in thermal pressure equilibrium. By introducing hotter, collisionally ionized phases, however, we can achieve balanced pressures. The best models indicate moderate metallicities, often with sub-solar N/$\alpha$, and, in two cases, ionizing flux fields that are softer and brighter than the fiducial Haardt & Madau UV background model., Comment: Accepted for publication in MNRAS

The propagation and evolution of cold galactic winds in galactic haloes is crucial to galaxy formation models. However, modelling of this process in hydrodynamic simulations of galaxy formation is over-simplified owing to a lack of numerical resolution and often neglects critical physical processes such as hydrodynamic instabilities and thermal conduction. We propose an analytic model, Physically Evolved Winds (PhEW), that calculates the evolution of individual clouds moving supersonically through a uniform ambient medium. Our model reproduces predictions from very high resolution cloud-crushing simulations that include isotropic thermal conduction over a wide range of physical conditions. We discuss the implementation of this model into cosmological hydrodynamic simulations of galaxy formation as a sub-grid prescription to model galactic winds more robustly both physically and numerically., Comment: 21 pages, 9 figures, accepted for publication in MNRAS

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., Comment: 31 pages, 19 figures, accepted by MNRAS

We study the mechanisms and evolutionary phases of bar formation in n-body simulations of a stellar disc and dark matter halo system using harmonic basis function expansion analysis to characterize the dynamical mechanisms in bar evolution. We correlate orbit families with phases of bar evolution by using empirical orthogonal functions that act as a spatial filter and form the gravitational potential basis. In both models we find evidence for three phases in evolution with unique harmonic signatures. We recover known analytic results, such as bar slowdown owing to angular momentum transfer. We also find new dynamical mechanisms for bar evolution: a steady-state equilibrium configuration and harmonic interaction resulting in harmonic mode locking, both of which may be observable. Additionally, we find that ellipse fitting may severely overestimate measurements of bar length by a factor of two relative to the measurements based on orbits that comprise the true backbone supporting the bar feature. The bias will lead to overestimates of both bar mass and bar pattern speed, affecting inferences about the evolution of bars in the real universe, such as the fraction of bars with fast pattern speeds. We propose a direct observational technique to compute the radial extent of trapped orbits and determine a dynamical length for the bar., Comment: Submitted to MNRAS, comments welcome

We track the angular momentum transfer in n-body simulations of barred galaxies by measuring torques to understand the dynamical mechanisms responsible for the evolution of the bar-disc-dark matter halo system. We find evidence for three distinct phases of barred galaxy evolution: assembly, secular growth, and steady-state equilibrium. Using a decomposition of the disc into orbital families, we track bar mass and angular momentum through time and correlate the quantities with the phases of evolution. We follow the angular momentum transfer between particles and identify the dominant torque channels. We find that the halo model mediates the assembly and growth of the bar for a high central density halo, and the outer disc mediates the assembly and growth of the bar in a low central density halo model. Both galaxies exhibit a steady-state equilibrium phase where the bar is neither lengthening nor slowing. The steady-state equilibrium results from the balance of torque between particles that are gaining and losing angular momentum. We propose observational metrics for barred galaxies that can be used to help determine the evolutionary phase of a barred galaxy, and discuss the implications of the phases for galaxy evolution as a whole., Comment: Accepted for publication in MNRAS

We interpret simulations of secularly-evolving disc galaxies through orbit morphology. Using a new algorithm that measures the volume of orbits in real space using a tessellation, we rapidly isolate commensurate (resonant) orbits. We identify phase-space regions occupied by different orbital families. Compared to spectral methods, the tessellation algorithm can identify resonant orbits within a few dynamical periods, crucial for understanding an evolving galaxy model. The flexible methodology accepts arbitrary potentials, enabling detailed descriptions of the orbital families. We apply the machinery to four different potential models, including two barred models, and fully characterise the orbital membership. We identify key differences in the content of orbit families, emphasising the presence of orbit families indicative of the bar evolutionary state and the shape of the dark matter halo. We use the characterisation of orbits to investigate the shortcomings of analytic and self-consistent studies, comparing our findings to the evolutionary epochs in self-consistent barred galaxy simulations. Using insight from our orbit analysis, we present a new observational metric that uses spatial and kinematic information from integral field spectrometers that may reveal signatures of commensurabilities and allow for a differentiation between models.

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., Comment: 29 pages, 25 figures, accepted by MNRAS

We survey the highly ionized circumgalactic media (CGM) of 29 blindly selected galaxies at 0.49 < z_(gal) < 1.44 based on high-S/N ultraviolet spectra of z > 1 QSOs and the galaxy database from the COS Absorption Survey of Baryon Harbors (CASBaH). We detect the Ne VIII doublet in nine of the galaxies, and for gas with N(Ne VIII) > 10^13.3 cm^-2 (> 10^13.5 cm^-2), we derive a Ne VIII covering fraction f_c = 75 +15/-25% (44 +22/-20%) within impact parameter (rho) < 200 kpc of M_* = 10^(9.5-11.5) Msol galaxies and f_c = 70 +16/-22% (f_c = 42 +20/-17%) within rho < 1.5 virial radii. We estimate the mass in Ne VIII-traced gas to be M_gas(Ne VIII) > 10^9.5 Msol (Z/Zsol)^-1, or 6-20% of the expected baryonic mass if the Ne VIII absorbers have solar metallicity. Ionizing Ne VII to Ne VIII requires 207 eV, and photons with this energy are scarce in the CGM. However, for the median halo mass and redshift of our sample, the virial temperature is close to the peak temperature for the Ne VIII ion, and the Ne VIII-bearing gas is plausibly collisionally ionized near this temperature. Moreover, we find that photoionized Ne VIII requires cool and low-density clouds that would be highly underpressured (by approximately two orders of magnitude) relative to the putative, ambient virialized medium, complicating scenarios where such clouds could survive. Thus, more complex (e.g., non-equilibrium) models may be required; this first statistical sample of Ne VIII absorber/galaxy systems will provide stringent constraints for future CGM studies., Comment: Published in ApJL, Volume 877, Issue 2, Article L20

We explore the application of Bayesian image analysis to infer the properties of an SDSS early-type galaxy sample including AGN. We use GALPHAT (Yoon et al. 2010) with a Bayes-factor model comparison to photometrically infer an AGN population and verify this using spectroscopic signatures. Our combined posterior sample for the SDSS sample reveals distinct low and high concentration modes after the point-source flux is modeled. This suggests that ETG parameters are intrinsically bimodal. The bimodal signature was weak when analyzed by GALFIT (Peng et al. 2002, 2010). This led us to create several ensembles of synthetic images to investigate the bias of inferred structural parameters and compare with GALFIT. GALPHAT inferences are less biased, especially for high-concentration profiles: GALPHAT S\'ersic index $n$, $r_{e}$ and MAG deviate from the true values by $6\%$, $7.6\%$ and $-0.03 \,\mathrm{mag}$, respectively, while GALFIT deviates by $15\%$, $22\%$ and $-0.09$\, mag, respectively. In addition, we explore the reliability for the photometric detection of AGN using Bayes factors. For our SDSS sample with $r_{e}\ge 7.92\,$arcsec, we correctly identify central point sources with $\mathrm{Mag_{PS}}-\mathrm{Mag_{Sersic}}\le 5$ for $n\le6$ and $\mathrm{Mag_{PS}}-\mathrm{Mag_{Sersic}}\le 3$ for $n>6$. The magnitude range increases and classification error decreases with increasing resolution, suggesting that this approach will excel for upcoming high-resolution surveys. Future work will extend this to models that test hypotheses of galaxy evolution through the cosmic time.

We examine the properties of the galaxies and dark matter haloes residing in the cluster infall region surrounding the simulated $\Lambda$CDM galaxy cluster studied by Elahi et al. (2016) at z=0. The $1.1\times10^{15}h^{-1}\text{M}_{\odot}$ galaxy cluster has been simulated with eight different hydrodynamical codes containing a variety of hydrodynamic solvers and subgrid schemes. All models completed a dark-matter only, non-radiative and full-physics run from the same initial conditions. The simulations contain dark matter and gas with mass resolution $m_{\text{DM}}=9.01\times 10^8h^{-1}\text{M}_{\odot}$ and $m_{\text{gas}}=1.9\times 10^8h^{-1}\text{M}_{\odot}$ respectively. We find that the synthetic cluster is surrounded by clear filamentary structures that contain ~60% of haloes in the infall region with mass ~$10^{12.5} - 10^{14} h^{-1}\text{M}_{\odot}$, including 2-3 group-sized haloes ($> 10^{13}h^{-1}\text{M}_{\odot}$). However, we find that only ~10% of objects in the infall region are subhaloes residing in haloes, which may suggest that there is not much ongoing preprocessing occurring in the infall region at z=0. By examining the baryonic content contained within the haloes, we also show that the code-to-code scatter in stellar fraction across all halo masses is typically ~2 orders of magnitude between the two most extreme cases, and this is predominantly due to the differences in subgrid schemes and calibration procedures that each model uses. Models that do not include AGN feedback typically produce too high stellar fractions compared to observations by at least ~1 order of magnitude., Comment: 14 pages, 9 figures, 1 table; accepted for publication in MNRAS

Over a handful of rotation periods, dynamical processes in barred galaxies induce non-axisymmetric structure in dark matter halos. Using n-body simulations of a Milky Way-like barred galaxy, we identify both a trapped dark-matter component, a shadow bar, and a strong response wake in the dark-matter distribution that affects the predicted dark-matter detection rates for current experiments. The presence of a baryonic disk together with well-known dynamical processes (e.g. spiral structure and bar instabilities) increase the dark matter density in the disk plane. We find that the magnitude of the combined stellar and shadow bar evolution, when isolated from the effect of the axisymmetric gravitational potential of the disk, accounts for >30% of this overall increase in disk-plane density. This is significantly larger that of previously claimed deviations from the standard halo model. The dark-matter density and kinematic wakes driven by the Milky Way bar increase the detectability of dark matter overall, especially for the experiments with higher $v_{min}$. These astrophysical features increase the detection rate by more than a factor of two when compared to the standard halo model and by a factor of ten for experiments with high minimum recoil energy thresholds. These same features increase (decrease) the annual modulation for low (high) minimum recoil energy experiments. We present physical arguments for why these dynamics are generic for barred galaxies such as the Milky Way rather than contingent on a specific galaxy model., Comment: 16 pages, 7 figures, submitted to Phys Rev D

Building on the initial results of the nIFTy simulated galaxy cluster comparison, we compare and contrast the impact of baryonic physics with a single massive galaxy cluster, run with 11 state-of-the-art codes, spanning adaptive mesh, moving mesh, classic and modern SPH approaches. For each code represented we have a dark matter only (DM) and non-radiative (NR) version of the cluster, as well as a full physics (FP) version for a subset of the codes. We compare both radial mass and kinematic profiles, as well as global measures of the cluster (e.g. concentration, spin, shape), in the NR and FP runs with that in the DM runs. Our analysis reveals good consistency (<= 20%) between global properties of the cluster predicted by different codes when integrated quantities are measured within the virial radius R200. However, we see larger differences for quantities within R2500, especially in the FP runs. The radial profiles reveal a diversity, especially in the cluster centre, between the NR runs, which can be understood straightforwardly from the division of codes into classic SPH and non-classic SPH (including the modern SPH, adaptive and moving mesh codes); and between the FP runs, which can also be understood broadly from the division of codes into those that include AGN feedback and those that do not. The variation with respect to the median is much larger in the FP runs with different baryonic physics prescriptions than in the NR runs with different hydrodynamics solvers., Comment: 24 pages, 13 figures, MNRAS submitted

We investigate the complex interactions between the stellar disc and the dark-matter halo during bar formation and evolution using N-body simulations with fine temporal resolution and optimally chosen spatial resolution. We find that the forming stellar bar traps dark matter in the vicinity of the stellar bar into bar-supporting orbits. We call this feature the shadow bar. The shadow bar modifies both the location and magnitude of the angular momentum transfer between the disc and dark matter halo and adds 10 per cent to the mass of the stellar bar over 4 Gyr. The shadow bar is potentially observable by its density and velocity signature in spheroid stars and by direct dark matter detection experiments. Numerical tests demonstrate that the shadow bar can diminish the rate of angular momentum transport from the bar to the dark matter halo by more than a factor of three over the rate predicted by dynamical friction with an untrapped dark halo, and thus provides a possible physical explanation for the observed prevalence of fast bars in nature., Comment: 16 pages, 12 figures. MNRAS accepted

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., Comment: 26 pages, 16 figures, accepted by MNRAS

Using Hubble Space Telescope Cosmic Origins Spectrograph observations of 89 QSO sightlines through the Sloan Digital Sky Survey footprint, we study the relationships between C IV absorption systems and the properties of nearby galaxies as well as large-scale environment. To maintain sensitivity to very faint galaxies, we restrict our sample to 0.0015 < z < 0.015, which defines a complete galaxy survey to L > 0.01 L* or stellar mass log M_* > 8 Msun. We report two principal findings. First, for galaxies with impact parameter rho < 1 rvir, C IV detection strongly depends on the luminosity/stellar mass of the nearby galaxy. C IV is preferentially associated with galaxies with log M_* > 9.5 Msun; lower mass galaxies rarely exhibit significant C IV absorption (covering fraction f = 9 +12-6% for 11 galaxies with log M_* < 9.5 Msun). Second, C IV detection within the log M_* > 9.5 Msun population depends on environment. Using a fixed-aperture environmental density metric for galaxies with rho < 160 kpc at z < 0.055, we find that 57+/-12% (8/14) of galaxies in low-density regions (regions with fewer than seven L > 0.15 L* galaxies within 1.5 Mpc) have affiliated C IV absorption; however, none (0/7) of the galaxies in denser regions show C IV. Similarly, the C IV detection rate is lower for galaxies residing in groups with dark-matter halo masses of log Mhalo > 12.5 Msun. In contrast to C IV, H I is pervasive in the CGM without regard to mass or environment. These results indicate that C IV absorbers with log N(C IV) > 13.5 cm^-2 trace the halos of log M_* > 9.5 Msun galaxies but also reflect larger scale environmental conditions., Comment: 26 pages, 13 figures. ApJ, in press

We examine subhaloes and galaxies residing in a simulated LCDM galaxy cluster ($M^{\rm crit}_{200}=1.1\times10^{15}M_\odot/h$) produced by hydrodynamical codes ranging from classic Smooth Particle Hydrodynamics (SPH), newer SPH codes, adaptive and moving mesh codes. These codes use subgrid models to capture galaxy formation physics. We compare how well these codes reproduce the same subhaloes/galaxies in gravity only, non-radiative hydrodynamics and full feedback physics runs by looking at the overall subhalo/galaxy distribution and on an individual objects basis. We find the subhalo population is reproduced to within $\lesssim10\%$ for both dark matter only and non-radiative runs, with individual objects showing code-to-code scatter of $\lesssim0.1$ dex, although the gas in non-radiative simulations shows significant scatter. Including feedback physics significantly increases the diversity. Subhalo mass and $V_{max}$ distributions vary by $\approx20\%$. The galaxy populations also show striking code-to-code variations. Although the Tully-Fisher relation is similar in almost all codes, the number of galaxies with $10^{9}M_\odot/h\lesssim M_*\lesssim 10^{12}M_\odot/h$ can differ by a factor of 4. Individual galaxies show code-to-code scatter of $\sim0.5$ dex in stellar mass. Moreover, strong systematic differences exist, with some codes producing galaxies $70\%$ smaller than others. The diversity partially arises from the inclusion/absence of AGN feedback. Our results combined with our companion papers demonstrate that subgrid physics is not just subject to fine-tuning, but the complexity of building galaxies in all environments remains a challenge. We argue even basic galaxy properties, such as the stellar mass to halo mass, should be treated with errors bars of $\sim0.2-0.4$ dex., Comment: 17 pages (+4 page appendix), 16 figures, 2 tables; accepted for publication in MNRAS

We have simulated the formation of a massive galaxy cluster (M$_{200}^{\rm crit}$ = 1.1$\times$10$^{15}h^{-1}M_{\odot}$) in a $\Lambda$CDM universe using 10 different codes (RAMSES, 2 incarnations of AREPO and 7 of GADGET), modeling hydrodynamics with full radiative subgrid physics. These codes include Smoothed-Particle Hydrodynamics (SPH), spanning traditional and advanced SPH schemes, adaptive mesh and moving mesh codes. Our goal is to study the consistency between simulated clusters modeled with different radiative physical implementations - such as cooling, star formation and AGN feedback. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. We find that, with respect to non-radiative simulations, dark matter is more centrally concentrated, the extent not simply depending on the presence/absence of AGN feedback. The scatter in global quantities is substantially higher than for non-radiative runs. Intriguingly, adding radiative physics seems to have washed away the marked code-based differences present in the entropy profile seen for non-radiative simulations in Sembolini et al. (2015): radiative physics + classic SPH can produce entropy cores. Furthermore, the inclusion/absence of AGN feedback is not the dividing line -as in the case of describing the stellar content- for whether a code produces an unrealistic temperature inversion and a falling central entropy profile. However, AGN feedback does strongly affect the overall stellar distribution, limiting the effect of overcooling and reducing sensibly the stellar fraction., Comment: 20 pages, 13 figures, submitted to MNRAS

To investigate the evolution of metal-enriched gas over recent cosmic epochs as well as to characterize the diffuse, ionized, metal-enriched circumgalactic medium (CGM), we have conducted a blind survey for C IV absorption systems in 89 QSO sightlines observed with the Hubble Space Telescope (HST) Cosmic Origins Spectrograph (COS). We have identified 42 absorbers at z < 0.16, comprising the largest uniform blind sample size to date in this redshift range. Our measurements indicate an increasing C IV absorber number density per comoving path length (dN/dX = 7.5 +/- 1.1) and modestly increasing mass density relative to the critical density of the Universe (Omega(C IV) = 10.0 +/- 1.5 x 10^-8 ) from z ~ 1.5 to the present epoch, consistent with predictions from cosmological hydrodynamical simulations. Furthermore, the data support a functional form for the column density distribution function that deviates from a single power-law, also consistent with independent theoretical predictions. As the data also probe heavy element ions in addition to C IV at the same redshifts, we identify, measure, and search for correlations between column densities of these species where components appear aligned in velocity. Among these ion-ion correlations, we find evidence for tight correlations between C II and Si II, C II and Si III, and C IV and Si IV, suggesting that these pairs of species arise in similar ionization conditions. However, the evidence for correlations decreases as the difference in ionization potential increases. Finally, when controlling for observational bias, we find only marginal evidence for a correlation (86.8% likelihood) between the Doppler line width b(C IV) and column density N(C IV)., Comment: 32 pages, 18 figures, 4 tables. Resubmitted to ApJ after referee review

We investigate the time evolution of non-parametric morphological quantities and their relationship to major mergers between $4\geq z \geq 2$ in high-resolution cosmological zoom simulations of disk galaxies that implement kinetic wind feedback, $H_2$-based star formation, and minimal ISM pressurisation. We show that the resulting galaxies broadly match basic observed physical properties of $z\sim 2$ objects. We measure the galaxies' concentrations ($C$), asymmetries ($A$), and $Gini$ ($G$) and $M_{20}$ coefficients, and correlate these with major merger events identified from the mass growth history. We find that high values of asymmetry provide the best indicator for identifying major mergers of $>1:4$ mass ratio within our sample, with $Gini$-$M_{20}\,$ merger classification only as effective for face-on systems and much less effective for edge-on or randomly-oriented galaxies. The canonical asymmetry cut of $A\geq0.35$, however, is only able to correctly identify major mergers $\sim 10\%$ of the time, while a higher cut of $A\geq 0.8$ more efficiently picks out mergers at this epoch. We further examine the temporal correlation between morphological statistics and mergers, and show that for randomly-oriented galaxies, half the galaxies with $A\geq0.8$ undergo a merger within $\pm0.2\,{\rm Gyr}$, whereas $Gini$-$M_{20}\,$ identification only identifies about a third correctly. The fraction improves further using $A\geq 1.5$, but about the half the mergers are missed by this stringent cut., Comment: 19 pages, 19 figures, submitted

We have simulated the formation of a galaxy cluster in a $\Lambda$CDM universe using twelve different codes modeling only gravity and non-radiative hydrodynamics (\art, \arepo, \hydra\ and 9 incarnations of GADGET). This range of codes includes particle based, moving and fixed mesh codes as well as both Eulerian and Lagrangian fluid schemes. The various GADGET implementations span traditional and advanced smoothed-particle hydrodynamics (SPH) schemes. The goal of this comparison is to assess the reliability of cosmological hydrodynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. The underlying gravitational framework can be aligned very accurately for all the codes allowing a detailed investigation of the differences that develop due to the various gas physics implementations employed. As expected, the mesh-based codes ART and AREPO form extended entropy cores in the gas with rising central gas temperatures. Those codes employing traditional SPH schemes show falling entropy profiles all the way into the very centre with correspondingly rising density profiles and central temperature inversions. We show that methods with modern SPH schemes that allow entropy mixing span the range between these two extremes and the latest SPH variants produce gas entropy profiles that are essentially indistinguishable from those obtained with grid based methods., Comment: 21 pages, 13 figures, 4 tables - submitted to MNRAS

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., Comment: 21 pages, 11 figures, submitted to MNRAS

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., Comment: 18 Pages, 16 Figures

The empirical model of Lu et al. 2014 is updated with recent data and used to study galaxy star formation and assembly histories. At $z > 2$, the predicted galaxy stellar mass functions are steep, and a significant amount of star formation is hosted by low-mass haloes that may be missed in current observations. Most of the stars in cluster centrals formed earlier than $z\approx 2$ but have been assembled much later. Milky Way mass galaxies have had on-going star formation without significant mergers since $z\approx 2$, and are thus free of significant (classic) bulges produced by major mergers. In massive clusters, stars bound in galaxies and scattered in the halo form a homogeneous population that is old and with solar metallicity. In contrast, in Milky Way mass systems the two components form two distinct populations, with halo stars being older and poorer in metals by a factor of $\approx 3$. Dwarf galaxies in haloes with $M_{\rm h} < 10^{11}h^{-1}M_{\odot}$ have experienced a star formation burst accompanied by major mergers at $z > 2$, followed by a nearly constant star formation rate after $z = 1$. The early burst leaves a significant old stellar population that is distributed in spheroids., Comment: 17 pages, 17 figures

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., Comment: 18 Pages, 11 Figures, ApJ 796 136

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., Comment: Submitted to ApJ Letters; 6 pages including 3 figures

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., Comment: 28 pages, 18 figures

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., Comment: 19 pages, 12 Figures, and a 37-page Appendix with 36 additional figures. Accepted to ApJ June 21 2014

We infer mechanisms of galaxy formation for a broad family of semi-analytic models (SAMs) constrained by the K-band luminosity function and HI mass function of local galaxies using tools of Bayesian analysis. Even with a broad search in parameter space the whole model family fails to match to constraining data. In the best fitting models, the star formation and feedback parameters in low-mass haloes are tightly constrained by the two data sets, and the analysis reveals several generic failures of models that similarly apply to other existing SAMs. First, based on the assumption that baryon accretion follows the dark matter accretion, large mass-loading factors are required for haloes with circular velocities lower than 200 km/s, and most of the wind mass must be expelled from the haloes. Second, assuming that the feedback is powered by Type-II supernovae with a Chabrier IMF, the outflow requires more than 25% of the available SN kinetic energy. Finally, the posterior predictive distributions for the star formation history are dramatically inconsistent with observations for masses similar to or smaller than the Milky-Way mass. The inferences suggest that the current model family is still missing some key physical processes that regulate the gas accretion and star formation in galaxies with masses below that of the Milky Way., Comment: 17 pages, 9 figures, 1 table, accepted for publication in MNRAS

Recent theoretical work has implicated fast bar formation modes and subsequent evolution as the creation mechanism for the observed peanut-shaped bulges in some edge-on disk galaxies. We demonstrate an N-body simulation of a disk undergoing a contrasting slow mode of bar growth, unsubjected to a buckling instability, which nonetheless grows the 4:1 orbit family responsible for a peanut-shaped bulge. We also present a simulation with fast mode bar growth, which exhibits thickening similar to other work. A novel orbit classification method that finds dynamically distinct families is presented, allowing for a detailed analysis of angular momentum transfer channels within the disk., Comment: 4 pages, 3 figures. To appear in the proceedings of the conference "The Structure and Dynamics of Disk Galaxies." M.S. Seigar and P. Treuthardt, eds

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

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., Comment: 19 pages, 12 figures, 2 tables. ApJ, in press

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., Comment: Accepted to ApJ, 87 pages, including 5 tables and 66 figures

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., Comment: 26 pages, 15 figures, replaced by published version

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., Comment: 25 pages, 13 Figures, submitted to MNRAS. Updated citations from original submission

We develop an empirical approach to infer the star formation rate in dark matter halos from the galaxy stellar mass function (SMF) at different redshifts and the local cluster galaxy luminosity function (CGLF), which has a steeper faint end relative to the SMF of local galaxies. As satellites are typically old galaxies which have been accreted earlier, this feature can cast important constraint on the formation of low-mass galaxies at high-redshift. The evolution of the SMFs suggests the star formation in high mass halos ($>10^{12}M_{\odot}h^{-1}$) has to be boosted at high redshift beyond what is expected from a simple scaling of the dynamical time. The faint end of the CGLF implies a characteristic redshift $z_c\approx2$ above which the star formation rate in low mass halos with masses $< 10^{11}M_{\odot}h^{-1}$ must be enhanced relative to that at lower z. This is not directly expected from the standard stellar feedback models. Also, this enhancement leads to some interesting predictions, for instance, a significant old stellar population in present-day dwarf galaxies with $M_* < 10^8 M_{\odot}h^{-2}$ and steep slopes of high redshift stellar mass and star formation rate functions., Comment: 27 pages, 19 figures, accepted by MNRAS

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., Comment: 21 pages, MNRAS, accepted

Weinberg (2012) described a constructive algorithm for computing the marginal likelihood, Z, from a Markov chain simulation of the posterior distribution. Its key point is: the choice of an integration subdomain that eliminates subvolumes with poor sampling owing to low tail-values of posterior probability. Conversely, this same idea may be used to choose the subdomain that optimizes the accuracy of Z. Here, we explore using the simulated distribution to define a small region of high posterior probability, followed by a numerical integration of the sample in the selected region using the volume tessellation algorithm described in Weinberg (2012). Even more promising is the resampling of this small region followed by a naive Monte Carlo integration. The new enhanced algorithm is computationally trivial and leads to a dramatic improvement in accuracy. For example, this application of the new algorithm to a four-component mixture with random locations in 16 dimensions yields accurate evaluation of Z with 5% errors. This enables Bayes-factor model selection for real-world problems that have been infeasible with previous methods., Comment: 14 pages, 3 figures, submitted to Bayesian Analysis

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., Comment: 6 pages, 5 figures. Accepted to ApJL

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., Comment: 26 pages, 15 figures, published in MNRAS. Updates to citations from previous version