Showing 1–50 of 58 results for author: Schneider, E

Dust Survival in Galactic Winds

Authors: Helena M. Richie, Evan E. Schneider, Matthew W. Abruzzo, Paul Torrey

Abstract: We present a suite of high-resolution numerical simulations to study the evolution and survival of dust in hot galactic winds. We implement a novel dust framework in the Cholla hydrodynamics code and use wind tunnel simulations of cool, dusty clouds to understand how thermal sputtering affects the dust content of galactic winds. Our simulations illustrate how various regimes of cloud evolution imp… ▽ More We present a suite of high-resolution numerical simulations to study the evolution and survival of dust in hot galactic winds. We implement a novel dust framework in the Cholla hydrodynamics code and use wind tunnel simulations of cool, dusty clouds to understand how thermal sputtering affects the dust content of galactic winds. Our simulations illustrate how various regimes of cloud evolution impact dust survival, dependent on cloud size, wind properties, and dust grain size. We find that significant amounts of dust can survive in winds in all scenarios, even without shielding from the cool phase of outflows. We present an analytic framework that explains this result, along with an analysis of the impact of cloud evolution on the total fraction of dust survival. Using these results, we estimate that 60 percent of dust that enters a starburst-driven wind could survive to populate the halo, based on a simulated distribution of cloud properties. We also investigate how these conclusions depend on grain size, exploring grains from 0.1 micron to 10 Angstrom. Under most circumstances, grains smaller than 0.01 micron cannot withstand hot-phase exposure, suggesting that the small grains observed in the CGM are either formed in situ due to shattering of larger grains, or must be carried there in cool phase outflows. Finally, we show that the dust-to-gas ratio of clouds declines as a function of distance from the galaxy due to cloud-wind mixing and condensation. These results provide an explanation for the vast amounts of dust observed in the CGMs of galaxies and beyond. △ Less

Submitted 6 March, 2024; originally announced March 2024.

Comments: Submitted to ApJ. 23 pages, 16 figures, 2 tables

CGOLS V: Disk-wide Stellar Feedback and Observational Implications of the Cholla Galactic Wind Model

Authors: Evan E. Schneider, S. Alwin Mao

Abstract: We present the fifth simulation in the CGOLS project -- a set of isolated starburst galaxy simulations modeled over large scales ($10\kpc$) at uniformly high resolution ($Δx \approx 5\pc$). Supernova feedback in this simulation is implemented as a disk-wide distribution of clusters, and we assess the impact of this geometry on several features of the resulting outflow, including radial profiles of… ▽ More We present the fifth simulation in the CGOLS project -- a set of isolated starburst galaxy simulations modeled over large scales ($10\kpc$) at uniformly high resolution ($Δx \approx 5\pc$). Supernova feedback in this simulation is implemented as a disk-wide distribution of clusters, and we assess the impact of this geometry on several features of the resulting outflow, including radial profiles of various phases; mass, momentum, and energy outflow rates; covering fraction of cool gas; mock absorption-line spectra; and X-ray surface brightness. In general, we find that the outflow generated by this model is cooler, slower, and contains more mass in the cool phase than a more centrally concentrated outflow driven by a similar number of supernovae. In addition, the energy loading factors in the hot phase are an order-of-magnitude lower, indicating much larger losses due to radiative cooling in the outflow. However, coupling between the hot and cool phases is more efficient than in the nuclear burst case, with almost 50\% of the total outflowing energy flux carried by the cool phase at a radial distance of 5 kpc. These physical differences have corresponding signatures in observable quantities: the covering fraction of cool gas is much larger, and there is greater evidence of absorption in low and intermediate ionization-energy lines. Taken together, our simulations indicate that centrally-concentrated starbursts are more effective at driving hot, low-density outflows that will expand far into the halo, while galaxy-wide bursts may be more effective at removing cool gas from the disk. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Comments: 22 pages, 13 figures, accepted in ApJ

Cholla-MHD: An Exascale-Capable Magnetohydrodynamic Extension to the Cholla Astrophysical Simulation Code

Authors: Robert V. Caddy, Evan E. Schneider

Abstract: We present an extension of the massively parallel, GPU native, astrophysical hydrodynamics code Cholla to magnetohydrodynamics (MHD). Cholla solves the ideal MHD equations in their Eulerian form on a static Cartesian mesh utilizing the Van Leer + Constrained Transport integrator, the HLLD Riemann solver, and reconstruction methods at second and third order. Cholla's MHD module can perform… ▽ More We present an extension of the massively parallel, GPU native, astrophysical hydrodynamics code Cholla to magnetohydrodynamics (MHD). Cholla solves the ideal MHD equations in their Eulerian form on a static Cartesian mesh utilizing the Van Leer + Constrained Transport integrator, the HLLD Riemann solver, and reconstruction methods at second and third order. Cholla's MHD module can perform $\approx260$ million cell updates per GPU-second on an NVIDIA A100 while using the HLLD Riemann solver and second order reconstruction. The inherently parallel nature of GPUs combined with increased memory in new hardware allows Cholla's MHD module to perform simulations with resolutions $\sim500^3$ cells on a single high end GPU (e.g. an NVIDIA A100 with 80GB of memory). We employ GPU direct MPI to attain excellent weak scaling on the exascale supercomputer \textit{Frontier}, while using 74,088 GPUs and simulating a total grid size of over 7.2 trillion cells. A suite of test problems highlights the accuracy of Cholla's MHD module and demonstrates that zero magnetic divergence in solutions is maintained to round off error. We also present new testing and continuous integration tools using GoogleTest, GitHub Actions, and Jenkins that have made development more robust and accurate and ensure reliability in the future. △ Less

Submitted 1 May, 2024; v1 submitted 7 February, 2024; originally announced February 2024.

Comments: 36 pages, 14 figures

CloudFlex: A Flexible Parametric Model for the Small-Scale Structure of the Circumgalactic Medium

Authors: Cameron B. Hummels, Kate H. R. Rubin, Evan E. Schneider, Drummond B. Fielding

Abstract: We present CloudFlex, a new open-source tool for predicting the absorption-line signatures of cool gas in galaxy halos with complex small-scale structure. Motivated by analyses of cool material in hydrodynamical simulations of turbulent, multiphase media, we model individual cool gas structures as assemblies of cloudlets with a power-law distribution of cloudlet mass $\propto m_{\rm cl}^{-α}$ and… ▽ More We present CloudFlex, a new open-source tool for predicting the absorption-line signatures of cool gas in galaxy halos with complex small-scale structure. Motivated by analyses of cool material in hydrodynamical simulations of turbulent, multiphase media, we model individual cool gas structures as assemblies of cloudlets with a power-law distribution of cloudlet mass $\propto m_{\rm cl}^{-α}$ and relative velocities drawn from a turbulent velocity field. The user may specify $α$, the lower limit of the cloudlet mass distribution ($m_{\rm cl,min}$), and several other parameters that set the total mass, size, and velocity distribution of the complex. We then calculate the MgII 2796 absorption profiles induced by the cloudlets along pencil-beam lines of sight. We demonstrate that at fixed metallicity, the covering fraction of sightlines with equivalent widths $W_{2796} < 0.3$ Ang increases significantly with decreasing $m_{\rm cl,min}$, cool cloudlet number density ($n_{\rm cl}$), and cloudlet complex size. We then present a first application, using this framework to predict the projected $W_{2796}$ distribution around ${\sim}L^*$ galaxies. We show that the observed incidences of $W_{2796}>0.3$ Ang sightlines within 10 kpc < $R_{\perp}$ < 50 kpc are consistent with our model over much of parameter space. However, they are underpredicted by models with $m_{\rm cl,min}\ge100M_{\odot}$ and $n_{\rm cl}\ge0.03$ $\rm cm^{-3}$, in keeping with a picture in which the inner cool circumgalactic medium (CGM) is dominated by numerous low-mass cloudlets ($m_{\rm cl}\lesssim100M_{\odot}$) with a volume filling factor ${\lesssim}1\%$. When used to simultaneously model absorption-line datasets built from multi-sightline and/or spatially-extended background probes, CloudFlex will enable detailed constraints on the size and velocity distributions of structures comprising the photoionized CGM. △ Less

Submitted 9 November, 2023; originally announced November 2023.

Comments: 22 pages, 7 figures. Submitted to AAS Journals, with minor modifications. Comments welcome. (1) Co-first authors who made equal contributions to this work

The Large Magellanic Cloud's $\sim30$ Kiloparsec Bow Shock and its Impact on the Circumgalactic Medium

Authors: David J. Setton, Gurtina Besla, Ekta Patel, Cameron Hummels, Yong Zheng, Evan Schneider

Abstract: The interaction between the supersonic motion of the Large Magellanic Cloud (LMC) and the Circumgalactic Medium (CGM) is expected to result in a bow shock that leads the LMC's gaseous disk. In this letter, we use hydrodynamic simulations of the LMC's recent infall to predict the extent of this shock and its effect on the Milky Way's (MW) CGM. The simulations clearly predict the existence of an asy… ▽ More The interaction between the supersonic motion of the Large Magellanic Cloud (LMC) and the Circumgalactic Medium (CGM) is expected to result in a bow shock that leads the LMC's gaseous disk. In this letter, we use hydrodynamic simulations of the LMC's recent infall to predict the extent of this shock and its effect on the Milky Way's (MW) CGM. The simulations clearly predict the existence of an asymmetric shock with a present day stand-off radius of $\sim6.7$ kpc and a transverse diameter of $\sim30$ kpc. Over the past 500 Myr, $\sim8\%$ of the MW's CGM in the southern hemisphere should have interacted with the shock front. This interaction may have had the effect of smoothing over inhomogeneities and increasing mixing in the MW CGM. We find observational evidence of the existence of the bow shock in recent $Hα$ maps of the LMC, providing a potential explanation for the envelope of ionized gas surrounding the LMC. Furthermore, the interaction of the bow shock with the MW CGM may also explain observations of ionized gas surrounding the Magellanic Stream. Using recent orbital histories of MW satellites, we find that many satellites have likely interacted with the LMC shock. Additionally, the dwarf galaxy Ret2 is currently sitting inside the shock, which may impact the interpretation of reported gamma ray excess in Ret2. This work highlights bow shocks associated with infalling satellites are an under-explored, yet potentially very important dynamical mixing process in the circumgalactic and intracluster media. △ Less

Submitted 17 November, 2023; v1 submitted 21 August, 2023; originally announced August 2023.

Comments: Re-uploaded after acceptance to ApJ Letters. 5 figures and 1 table. Comments welcome!

Highly-mass-loaded hot galactic winds are unstable to cool filament formation

Authors: Dustin D. Nguyen, Todd A. Thompson, Evan E. Schneider, Ashley P. Tarrant

Abstract: When cool clouds are ram-pressure accelerated by a hot supersonic galactic wind, some of the clouds may be shredded by hydrodynamical instabilities and incorporated into the hot flow. Recent one-dimensional steady-state calculations show how cool cloud entrainment directly affects the bulk thermodynamics, kinematics, and observational characteristics of the hot gas. In particular, mass-loading dec… ▽ More When cool clouds are ram-pressure accelerated by a hot supersonic galactic wind, some of the clouds may be shredded by hydrodynamical instabilities and incorporated into the hot flow. Recent one-dimensional steady-state calculations show how cool cloud entrainment directly affects the bulk thermodynamics, kinematics, and observational characteristics of the hot gas. In particular, mass-loading decelerates the hot flow and changes its entropy. Here, we investigate the stability of planar and spherical mass-loaded hot supersonic flows using both perturbation analysis and three-dimensional time-dependent radiative hydrodynamical simulations. We show that mass-loading is stable over a broad range of parameters and that the 1D time-steady analytic solutions exactly reproduce the 3D time-dependent calculations, provided that the flow does not decelerate sufficiently to become subsonic. For higher values of the mass-loading, the flow develops a sonic point and becomes thermally unstable, rapidly cooling and forming elongated dense cometary filaments. We explore the mass-loading parameters required to reach a sonic point and the radiative formation of these filaments. For certain approximations, we can derive simple analytic criteria. In general a mass-loading rate similar to the initial mass outflow rate is required. In this sense, the destruction of small cool clouds by a hot flow may ultimately spontaneously generate fast cool filaments, as observed in starburst superwinds. Lastly, we find that the kinematics of filaments is sensitive to the slope of the mass-loading function. Filaments move faster than the surrounding wind if mass-loading is over long distances whereas filaments move slower than their surroundings if mass-loading is abrupt. △ Less

Submitted 21 July, 2023; originally announced July 2023.

Comments: 12 pages, 15 figures, submitted to MNRAS (21 July 2023)

Searching in HI for Massive Low Surface Brightness Galaxies: Samples from HyperLeda and the UGC

Authors: K. O'Neil, Stephan E. Schneider, W. van Driel, G. Liu, T. Joseph, A. C. Schwortz, Z. Butcher

Abstract: A search has been made for 21 cm HI line emission in a total of 350 unique galaxies from two samples whose optical properties indicate they may be massive The first consists of 241 low surface brightness (LSB) galaxies of morphological type Sb and later selected from the HyperLeda database and the the second consists of 119 LSB galaxies from the UGC with morphological types Sd-m and later. Of the… ▽ More A search has been made for 21 cm HI line emission in a total of 350 unique galaxies from two samples whose optical properties indicate they may be massive The first consists of 241 low surface brightness (LSB) galaxies of morphological type Sb and later selected from the HyperLeda database and the the second consists of 119 LSB galaxies from the UGC with morphological types Sd-m and later. Of the 350 unique galaxies, 239 were observed at the Nancay Radio Telescope, 161 at the Green Bank Telescope, and 66 at the Arecibo telescope. A total of 295 (84.3%) were detected, of which 253 (72.3%) appear to be uncontaminated by any other galaxies within the telescope beam. Finally, of the total detected, uncontaminated galaxies, at least 31 appear to be massive LSB galaxies, with a total HI mass $\ge$ 10$^{10}$ M$_{sol}$, for H$_0$ = 70 km/s/Mpc. If we expand the definition to also include galaxies with significant total (rather than just gas) mass, i.e., those with inclination-corrected HI line width W$_{50}$,cor > 500 km/s, this bring the total number of massive LSB galaxies to 41. There are no obvious trends between the various measured global galaxy properties, particularly between mean surface brightness and galaxy mass. △ Less

Submitted 18 July, 2023; originally announced July 2023.

Comments: 71 pages, including all tables and figures; Accepted by AJ

Dynamics of hot galactic winds launched from spherically-stratified starburst cores

Authors: Dustin D. Nguyen, Todd A. Thompson, Evan E. Schneider, Sebastian Lopez, Laura A. Lopez

Abstract: The analytic galactic wind model derived by Chevalier and Clegg in 1985 (CC85) assumes $\textit{uniform}$ energy and mass-injection within the starburst galaxy nucleus. However, the structure of nuclear star clusters, bulges, and star-forming knots are non-uniform. We generalize to cases with spherically-symmetric energy/mass injection that scale as $r^{-Δ}$ within the starburst volume $R$, provid… ▽ More The analytic galactic wind model derived by Chevalier and Clegg in 1985 (CC85) assumes $\textit{uniform}$ energy and mass-injection within the starburst galaxy nucleus. However, the structure of nuclear star clusters, bulges, and star-forming knots are non-uniform. We generalize to cases with spherically-symmetric energy/mass injection that scale as $r^{-Δ}$ within the starburst volume $R$, providing solutions for $Δ= 0$, 1/2, 1, 3/2, and 2. In marked contrast with the CC85 model ($Δ=0$), which predicts zero velocity at the center, for a singular isothermal sphere profile ($Δ=2$), we find that the flow maintains a $\textit{constant}$ Mach number of $\mathcal{M}=\sqrt{3/5} \simeq 0.77$ throughout the volume. The fast interior flow can be written as $v_{r < R} = (\dot{E}_T/3\dot{M}_T)^{1/2} \simeq 0.41 \, v_\infty$, where $v_\infty$ is the asymptotic velocity, and $\dot{E}_T$ and $\dot{M}_T$ are the total energy and mass injection rates. For $v_\infty \simeq 2000 \, \mathrm{km \, s^{-1}}$, $v_{r<R} \simeq 820 \, \mathrm{km\, s^{-1}}$ throughout the wind-driving region. The temperature and density profiles of the non-uniform models may be important for interpreting spatially-resolved maps of starburst nuclei. We compute velocity resolved spectra to contrast the $Δ=0$ (CC85) and $Δ=2$ models. Next generation X-ray space telescopes such as XRISM may assess these kinematic predictions. △ Less

Submitted 8 November, 2022; v1 submitted 13 October, 2022; originally announced October 2022.

Comments: submitted (13 October 2022), accepted MNRAS Letters (08 November 2022)

New Constraints on Warm Dark Matter from the Lyman-$α$ Forest Power Spectrum

Authors: Bruno Villasenor, Brant Robertson, Piero Madau, Evan Schneider

Abstract: The forest of Lyman-$α$ absorption lines detected in the spectra of distant quasars encodes information on the nature and properties of dark matter and the thermodynamics of diffuse baryonic material. Its main observable -- the 1D flux power spectrum (FPS) -- should exhibit a suppression on small scales and an enhancement on large scales in warm dark matter (WDM) cosmologies compared to standard… ▽ More The forest of Lyman-$α$ absorption lines detected in the spectra of distant quasars encodes information on the nature and properties of dark matter and the thermodynamics of diffuse baryonic material. Its main observable -- the 1D flux power spectrum (FPS) -- should exhibit a suppression on small scales and an enhancement on large scales in warm dark matter (WDM) cosmologies compared to standard $Λ$CDM. Here, we present an unprecedented suite of 1080 high-resolution cosmological hydrodynamical simulations run with the Graphics Processing Unit-accelerated code Cholla to study the evolution of the Lyman-$α$ forest under a wide range of physically-motivated gas thermal histories along with different free-streaming lengths of WDM thermal relics in the early Universe. A statistical comparison of synthetic data with the forest FPS measured down to the smallest velocity scales ever probed at redshifts $4.0\lesssim z\lesssim 5.2$ (Boera et al. 2019) yields a lower limit $m_{\rm WDM}>3.1$ keV (95 percent CL) for the WDM particle mass and constrains the amplitude and spectrum of the photoheating and photoionizing background produced by star-forming galaxies and active galactic nuclei at these redshifts. Interestingly, our Bayesian inference analysis appears to weakly favor WDM models with a peak likelihood value at the thermal relic mass of $m_{\rm WDM}=4.5$ keV. We find that the suppression of the FPS from free-streaming saturates at $k\gtrsim 0.1\,\mathrm{s}\,\mathrm{km}^{-1}$ because of peculiar velocity smearing, and this saturated suppression combined with a slightly lower gas temperature provides a moderately better fit to the observed small-scale FPS for WDM cosmologies. △ Less

Submitted 9 June, 2023; v1 submitted 28 September, 2022; originally announced September 2022.

Comments: Version accepted to Phys. Rev. D

Revealing the Galaxy-Halo Connection Through Machine Learning

Authors: Ryan Hausen, Brant E. Robertson, Hanjue Zhu, Nickolay Y. Gnedin, Piero Madau, Evan E. Schneider, Bruno Villasenor, Nicole E. Drakos

Abstract: Understanding the connections between galaxy stellar mass, star formation rate, and dark matter halo mass represents a key goal of the theory of galaxy formation. Cosmological simulations that include hydrodynamics, physical treatments of star formation, feedback from supernovae, and the radiative transfer of ionizing photons can capture the processes relevant for establishing these connections. T… ▽ More Understanding the connections between galaxy stellar mass, star formation rate, and dark matter halo mass represents a key goal of the theory of galaxy formation. Cosmological simulations that include hydrodynamics, physical treatments of star formation, feedback from supernovae, and the radiative transfer of ionizing photons can capture the processes relevant for establishing these connections. The complexity of these physics can prove difficult to disentangle and obfuscate how mass-dependent trends in the galaxy population originate. Here, we train a machine learning method called Explainable Boosting Machines (EBMs) to infer how the stellar mass and star formation rate of nearly 6 million galaxies simulated by the Cosmic Reionization on Computers (CROC) project depend on the physical properties of halo mass, the peak circular velocity of the galaxy during its formation history $v_\mathrm{peak}$, cosmic environment, and redshift. The resulting EBM models reveal the relative importance of these properties in setting galaxy stellar mass and star formation rate, with $v_\mathrm{peak}$ providing the most dominant contribution. Environmental properties provide substantial improvements for modeling the stellar mass and star formation rate in only $\lesssim10\%$ of the simulated galaxies. We also provide alternative formulations of EBM models that enable low-resolution simulations, which cannot track the interior structure of dark matter halos, to predict the stellar mass and star formation rate of galaxies computed by high-resolution simulations with detailed baryonic physics. △ Less

Submitted 21 April, 2022; originally announced April 2022.

Comments: 27 pages, 18 figures, Submitted to AAS Journals

Inferring the Thermal History of the Intergalactic Medium from the Properties of the Hydrogen and Helium Lyman-alpha Forest

Authors: Bruno Villasenor, Brant Robertson, Piero Madau, Evan Schneider

Abstract: The filamentary network of intergalactic medium (IGM) gas that gives origin to the Lyman-alpha forest in the spectra of distant quasars encodes information on the physics of structure formation and the early thermodynamics of diffuse baryonic material. Here, we use a massive suite of more than 400 high-resolution cosmological hydrodynamical simulations run with the Graphics Processing Unit-acceler… ▽ More The filamentary network of intergalactic medium (IGM) gas that gives origin to the Lyman-alpha forest in the spectra of distant quasars encodes information on the physics of structure formation and the early thermodynamics of diffuse baryonic material. Here, we use a massive suite of more than 400 high-resolution cosmological hydrodynamical simulations run with the Graphics Processing Unit-accelerated code Cholla to study the IGM at high spatial resolution maintained over the entire computational volume. The simulations capture a wide range of possible thermal histories of intergalactic gas by varying the amplitude and timing of the photoheating and photoionizing background produced by star-forming galaxies and active galactic nuclei. A statistical comparison of synthetic spectra with the observed 1D flux power spectra of hydrogen in 14 redshift bins over the full range 2.2 <= z <= 5.0 and with the Lyman-alpha opacity of helium in 5 redshift bins over the range 2.4 < z < 2.9 tightly constrains the photoionization and photoheating history of the IGM. By leveraging the constraining power of the available Lyman-alpha forest data to break model degeneracies, we find that the IGM experienced two main reheating events from the non-equilibrium ionization of hydrogen and helium over 1.2 Gyr of cosmic time. For our best-fit model, hydrogen reionization completes by z~6.0 with a first IGM temperature peak T_0 ~ 1.3x10^4 K, and is followed by the reionization of HeII that completes by z~3.0 and yields a second temperature peak of T_0 ~ 1.4x10^4 K. We discuss how our results can be used to obtain information on the timing and the sources of hydrogen and helium reionization. △ Less

Submitted 5 June, 2022; v1 submitted 29 October, 2021; originally announced November 2021.

Comments: 33 pages, 22 figures. Accepted for publication in AAS Journal. Revised version (V3)

Simulated Observations of Multiphase Galactic Winds

Authors: Lita M. de la Cruz, Evan E. Schneider, Eve C. Ostriker

Abstract: Supernova-driven galactic winds are multiphase streams of gas that are often observed flowing at a range of velocities out of star-forming regions in galaxies. In this study, we use high resolution 3D simulations of multiphase galactic winds modeled with the hydrodynamics code Cholla to investigate the connection between numerical studies and observations. Using a simulated interaction between a h… ▽ More Supernova-driven galactic winds are multiphase streams of gas that are often observed flowing at a range of velocities out of star-forming regions in galaxies. In this study, we use high resolution 3D simulations of multiphase galactic winds modeled with the hydrodynamics code Cholla to investigate the connection between numerical studies and observations. Using a simulated interaction between a hot $T\sim 10^{7}\,\rm K$ supernova-driven wind and a cool $T\sim 10^{4}\,\rm K$ cloud of interstellar material, we create mock observables, including the optical depth and covering fraction of six commonly observed ions (Si II, C II, Si IV, C IV, N V, and O VI) as a function of gas velocity. We compare our mock observables to surveys of galactic winds in the literature, finding good agreement with velocities and profiles of the low ions. We then compute "empirical" values for the optical depth and covering fraction following observational techniques, and compare them to the values calculated directly from the simulation data. We find that the empirically computed values tend to underestimate the "true" value of $τ$ for ions with high optical depth, and overestimate the "true" value of $τ$ for ions with low optical depth, relative to the simulated data. △ Less

Submitted 3 December, 2020; originally announced December 2020.

The Characteristic Momentum of Radiatively Cooling Energy-Driven Galactic Winds

Authors: Cassandra Lochhaas, Todd A. Thompson, Evan E. Schneider

Abstract: Energy injection by supernovae may drive hot supersonic galactic winds in rapidly star-forming galaxies, driving metal-enriched gas into the circumgalactic medium and potentially accelerating cool gas. If sufficiently mass-loaded, such flows become radiative within the wind-driving region, reducing the overall mass outflow rate from the host galaxy. We show that this sets a maximum on the total ou… ▽ More Energy injection by supernovae may drive hot supersonic galactic winds in rapidly star-forming galaxies, driving metal-enriched gas into the circumgalactic medium and potentially accelerating cool gas. If sufficiently mass-loaded, such flows become radiative within the wind-driving region, reducing the overall mass outflow rate from the host galaxy. We show that this sets a maximum on the total outflow momentum for hot energy-driven winds. For a spherical wind of Solar metallicity driven by continuous star formation, $\dot p_\rm{max} \simeq 1.9\times10^4\ M_\odot/\rm{yr\ km/s}\ (α/0.9)^{0.86}(R_\star/300\ \rm{pc})^{0.14}(\dot M_\star/20\ M_\odot/\rm{yr})^{0.86}$, where $α$ is the fraction of supernova energy that thermalizes the wind, and $\dot M_\star$ and $R_\star$ are the star formation rate and radius of the wind-driving region. This maximum momentum for hot winds can also apply to cool, ionized outflows that are typically observed in starburst galaxies, if the hot wind undergoes bulk radiative cooling or if the hot wind transfers mass and momentum to cool clouds within the flow. We show that requiring the hot wind to undergo single-phase cooling on large scales sets a minimum on the total outflow momentum rate. These maximum and minimum outflow momenta have similar values, setting a characteristic momentum rate of hot galactic winds that can become radiative on large scales. We find that most observations of photoionized outflow wind momentum fall below the theoretical maximum and thus may be signatures of cooling hot flows. On the other hand, many systems fall below the minimum momentum required for bulk cooling, indicating that perhaps the cool material observed has instead been entrained in or mixed with the hot flow. △ Less

Submitted 15 April, 2021; v1 submitted 11 November, 2020; originally announced November 2020.

Comments: 14 pages, 7 figures, accepted by MNRAS

Effects of Photoionization and Photoheating on Lyman-alpha Forest Properties from Cholla Cosmological Simulations

Authors: Bruno Villasenor, Brant Robertson, Piero Madau, Evan Schneider

Abstract: The density and temperature properties of the intergalactic medium (IGM) reflect the heating and ionization history during cosmological structure formation, and are primarily probed by the Lyman-alpha forest of neutral hydrogen absorption features in the observed spectra of background sources (Gunn & Peterson 1965). We present the methodology and initial results from the Cholla IGM Photoheating Si… ▽ More The density and temperature properties of the intergalactic medium (IGM) reflect the heating and ionization history during cosmological structure formation, and are primarily probed by the Lyman-alpha forest of neutral hydrogen absorption features in the observed spectra of background sources (Gunn & Peterson 1965). We present the methodology and initial results from the Cholla IGM Photoheating Simulation (CHIPS) suite performed with the Graphics Process Unit-accelerated Cholla code to study the IGM at high, uniform spatial resolution maintained over large volumes. In this first paper, we examine the IGM structure in CHIPS cosmological simulations that include IGM uniform photoheating and photoionization models where hydrogen reionization completes early (Haardt & Madau 2012) or by redshift z~6 (Puchwein et al. 2019). Comparing with observations of the large- and small-scale Lyman-alpha transmitted flux power spectra P(k) at redshifts 2 <~ z <~ 5.5, the relative agreement of the models depends on scale, with the self-consistent Puchwein et al. (2019) IGM photoheating and photoionization model in good agreement with the flux P(k) at k >~ 0.01 s/km at redshifts 2 <~ z <~ 3.5. On larger scales the P(k) measurements increase in amplitude from z~4.6 to z~2.2 faster than the models, and lie in between the model predictions at 2.2 <~ z <~ 4.6 for k~= 0.002-0.01 s/km. We argue the models could improve by changing the HeII photoheating rate associated with active galactic nuclei to reduce the IGM temperature at z~3. At higher redshifts z>~4.5 the observed flux P(k) amplitude increases at a rate intermediate between the models, and we argue that for models where hydrogen reionization completes late (z~5.5 - 6) resolving this disagreement will require inhomogeneous or 'patchy' reionization. (Abridged) △ Less

Submitted 24 March, 2021; v1 submitted 14 September, 2020; originally announced September 2020.

Comments: 30 pages, 19 figures, 2 tables. Accepted for publication at ApJ

First results from SMAUG: Uncovering the Origin of the Multiphase Circumgalactic Medium with a Comparative Analysis of Idealized and Cosmological Simulations

Authors: Drummond B. Fielding, Stephanie Tonnesen, Daniel DeFelippis, Miao Li, Kung-Yi Su, Greg L. Bryan, Chang-Goo Kim, John C. Forbes, Rachel S. Somerville, Nicholas Battaglia, Evan E. Schneider, Yuan Li, Ena Choi, Christopher C. Hayward, Lars Hernquist

Abstract: We examine the properties of the circumgalactic medium (CGM) at low redshift in a range of simulated Milky Way mass halos. The sample is comprised of seven idealized simulations, an adaptive mesh refinement cosmological zoom-in simulation, and two groups of 50 halos with star forming or quiescent galaxies taken from the IllustrisTNG100 simulation. The simulations have very different setups, resolu… ▽ More We examine the properties of the circumgalactic medium (CGM) at low redshift in a range of simulated Milky Way mass halos. The sample is comprised of seven idealized simulations, an adaptive mesh refinement cosmological zoom-in simulation, and two groups of 50 halos with star forming or quiescent galaxies taken from the IllustrisTNG100 simulation. The simulations have very different setups, resolution, and feedback models, but are analyzed in a uniform manner. By comparing median radial profiles and mass distributions of CGM properties, we isolate key similarities and differences. In doing so, we advance the efforts of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project that aims to understand the inherently multiscale galaxy formation process. In the cosmological simulations, the CGM exhibits nearly flat temperature distributions, and broad pressure and radial velocity distributions. In the idealized simulations, similar distributions are found in the inner CGM ($\lesssim 0.5 \, r_{\rm 200c}$) when strong galactic feedback models are employed, but the outer CGM ($\gtrsim 0.5 \, r_{\rm 200c}$) has a much less prominent cold phase, and narrower pressure and velocity distributions even in models with strong feedback. This comparative analysis demonstrates the dominant role feedback plays in shaping the inner CGM and the increased importance of cosmological effects, such as nonspherical accretion and satellite galaxies, in the outer CGM. Furthermore, our findings highlight that while cosmological simulations are required to capture the multiphase structure of the CGM at large radii, idealized simulations provide a robust framework to study how galactic feedback interacts with the inner CGM and thereby provide a reliable avenue to constrain feedback prescriptions. △ Less

Submitted 29 June, 2020; originally announced June 2020.

Comments: More information on the SMAUG project and the other first result papers can be found here: https://www.simonsfoundation.org/flatiron/center-for-computational-astrophysics/galaxy-formation/smaug/papersplash1

The Physical Nature of Starburst-Driven Galactic Outflows

Authors: Evan E. Schneider, Eve C. Ostriker, Brant E. Robertson, Todd A. Thompson

Abstract: We present the fourth of the Cholla Galactic OutfLow Simulations suite (CGOLS). Using a physically-motivated prescription for clustered supernova feedback, we successfully drive a multiphase outflow from a disk galaxy. The high resolution ($< 5\,\mathrm{pc}$) across a relatively large domain ($20\,\mathrm{kpc}$) allows us to capture the hydrodynamic mixing and dynamical interactions between the ho… ▽ More We present the fourth of the Cholla Galactic OutfLow Simulations suite (CGOLS). Using a physically-motivated prescription for clustered supernova feedback, we successfully drive a multiphase outflow from a disk galaxy. The high resolution ($< 5\,\mathrm{pc}$) across a relatively large domain ($20\,\mathrm{kpc}$) allows us to capture the hydrodynamic mixing and dynamical interactions between the hot and cool ($T \sim 10^4\,\mathrm{K}$) phases in the outflow, which in turn leads to direct evidence of a qualitatively new mechanism for cool gas acceleration in galactic winds. We show that mixing of momentum from the hot phase to the cool phase accelerates the cool gas to $800\,\mathrm{km}\,\mathrm{s}^{-1}$ on kpc scales, with properties inconsistent with the physical models of ram pressure acceleration or with bulk cooling from the hot phase. The mixing process also affects the hot phase, modifying its radial profiles of temperature, density, and velocity from the expectations of radial supersonic flow. This mechanism provides a physical explanation for the high velocity, blue shifted, low ionization absorption lines often observed in the spectra of starburst and high redshift galaxies. △ Less

Submitted 24 February, 2020; originally announced February 2020.

Comments: 29 pages, 21 figures. Submitted to ApJ. Comments welcome

Astro2020: Training the Future Generation of Computational Researchers

Authors: Gurtina Besla, Daniela Huppenkothen, Nicole Lloyd-Ronning, Evan Schneider, Peter Behroozi, Blakesley Burkhart, C. K. Chan, Seth A. Jacobson, Sarah Morrison, Hai Ah Nam, Smadar Naoz, Annika Peter, Enrico Ramirez-Ruiz

Abstract: The current disparity in computational knowledge is a critical hindrance to the diversity and success of the field. Recommendations are outlined for policies and funding models to enable the growth and retention of a new generation of computational researchers that reflect the demographics of the undergraduate population in Astronomy and Physics. The current disparity in computational knowledge is a critical hindrance to the diversity and success of the field. Recommendations are outlined for policies and funding models to enable the growth and retention of a new generation of computational researchers that reflect the demographics of the undergraduate population in Astronomy and Physics. △ Less

Submitted 9 July, 2019; originally announced July 2019.

Comments: Astro2020 APC White Paper: State of the Profession Consideration

Infrared Galaxies in the Field of the Massive Cluster Abell S1063: Discovery of a Luminous Kiloparsec-Sized HII Region in a Gravitationally Lensed IR-Luminous Galaxy at $z=0.6$

Authors: Gregory L. Walth, Eiichi Egami, Benjamin Clément, Timothy D. Rawle, Marie Rex, Johan Richard, Pablo Pérez-González, Frédéric Boone, Miroslava Dessauges-Zavadsky, Jeff Portouw, Benjamin Weiner, Ian McGreer, Evan Schneider

Abstract: Using the Spitzer Space Telescope and Herschel Space Observatory, we have conducted a survey of infrared galaxies in the field of the galaxy cluster Abell S1063 (AS1063) at $z=0.347$, which is one of the most massive clusters known and a target of the HST CLASH and Frontier-Field surveys. The Spitzer/MIPS 24 $μ$m and Herschel/PACS & SPIRE images revealed that the core of AS1063 is surprisingly dev… ▽ More Using the Spitzer Space Telescope and Herschel Space Observatory, we have conducted a survey of infrared galaxies in the field of the galaxy cluster Abell S1063 (AS1063) at $z=0.347$, which is one of the most massive clusters known and a target of the HST CLASH and Frontier-Field surveys. The Spitzer/MIPS 24 $μ$m and Herschel/PACS & SPIRE images revealed that the core of AS1063 is surprisingly devoid of infrared sources, showing only a few detectable sources within the central r$\sim1^{\prime}$. There is, however, one particularly bright source (2.3 mJy at 24 $μ$m; 106 mJy at 160 $μ$m), which corresponds to a background galaxy at $z=0.61$. The modest magnification factor (4.0$\times$) implies that this galaxy is intrinsically IR-luminous (L$_{\rm IR}=3.1\times10^{11}\ \rm L_{\odot}$). What is particularly interesting about this galaxy is that HST optical/near-infrared images show a remarkably bright and large (1 kpc) clump at one edge of the disk. Our follow-up optical/near-infrared spectroscopy shows Balmer (H$α$-H8) and forbidden emission from this clump ([OII] $λ$3727, [OIII] $λλ$4959,5007, [NII] $λλ$6548,6583), indicating that it is a HII region. The HII region appears to have formed in-situ, as kinematically it is part of a rotating disk, and there is no evidence of nearby interacting galaxies. With an extinction correction of A$_{\rm V}=1.5$ mag, the star formation rate of this giant HII region is $\sim$10 M$_{\odot}$ yr$^{-1}$, which is exceptionally large, even for high redshift HII regions. Such a large and luminous HII region is often seen at $z\sim2$ but quite rare in the nearby Universe. △ Less

Submitted 15 April, 2019; originally announced April 2019.

Comments: Accepted for publication in the Astrophysical Journal. 29 pages, 15 figures

Astro2020 Science White Paper: Cold Gas Outflows, Feedback, and the Shaping of Galaxies

Authors: Alberto D. Bolatto, Lee Armus, Sylvain Veilleux, Adam K. Leroy, Fabian Walter, Richard Mushotzky, Karin M. Sandstrom, Paul Martini, Evan E. Schneider, Tony Wong, Roberto Decarli, Caitlin Casey, Dominik Riechers, David Meier, Desika Narayana

Abstract: There is wide consensus that galaxy outflows are one of the most important processes determining the evolution of galaxies through cosmic time, for example playing a key role in shaping the galaxy mass function. Our understanding of outflows and their drivers, however, is in its infancy --- this is particularly true for the cold (neutral atomic and molecular) phases of outflows, which present obse… ▽ More There is wide consensus that galaxy outflows are one of the most important processes determining the evolution of galaxies through cosmic time, for example playing a key role in shaping the galaxy mass function. Our understanding of outflows and their drivers, however, is in its infancy --- this is particularly true for the cold (neutral atomic and molecular) phases of outflows, which present observational and modeling challenges. Here we outline several key open questions, briefly discussing the requirements of the observations necessary to make progress, and the relevance of several existing and planned facilities. It is clear that galaxy outflows, and particularly cold outflows, will remain a topic of active research for the next decade and beyond. △ Less

Submitted 3 April, 2019; originally announced April 2019.

Comments: Submitted as a Science White Paper to the Astro2020 Decadal Survey. arXiv admin note: substantial text overlap with arXiv:1810.06737

Galaxy Winds in the Age of Hyperdimensional Astrophysics

Authors: Grant R. Tremblay, Evan E. Schneider, Alexey Vikhlinin, Lars Hernquist, Mateusz Ruszkowski, Benjamin D. Oppenheimer, Ralph P. Kraft, John ZuHone, Michael A. McDonald, Massimo Gaspari, Megan Donahue, G. Mark Voit

Abstract: 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. Facilit… ▽ More 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. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Comments: A Science White Paper submitted to the Astro2020 Decadal Survey (7 papers, 4 figures)

Introducing CGOLS: The Cholla Galactic OutfLow Simulation Suite

Authors: Evan E. Schneider, Brant E. Robertson

Abstract: We present the Cholla Galactic OutfLow Simulations (CGOLS) suite, a set of extremely high resolution global simulations of isolated disk galaxies designed to clarify the nature of multiphase structure in galactic winds. Using the GPU-based code Cholla, we achieve unprecedented resolution in these simulations, modeling galaxies over a 20 kpc region at a constant resolution of 5 pc. The simulations… ▽ More We present the Cholla Galactic OutfLow Simulations (CGOLS) suite, a set of extremely high resolution global simulations of isolated disk galaxies designed to clarify the nature of multiphase structure in galactic winds. Using the GPU-based code Cholla, we achieve unprecedented resolution in these simulations, modeling galaxies over a 20 kpc region at a constant resolution of 5 pc. The simulations include a feedback model designed to test the effects of different mass- and energy-loading factors on galactic outflows over kiloparsec scales. In addition to describing the simulation methodology in detail, we also present the results from an adiabatic simulation that tests the frequently-adopted analytic galactic wind model of Chevalier & Clegg. Our results indicate that the Chevalier & Clegg model is a good fit to nuclear starburst winds in the non-radiative region of parameter space. Lastly, we investigate the role of resolution and convergence in large-scale simulations of multiphase galactic winds. While our largest-scale simulations do show convergence of observable features like soft X-ray emission, our tests demonstrate that simulations of this kind with resolutions greater than 10 pc are not yet converged, confirming the need for extreme resolution in order to study the structure of winds and their effects on the circumgalactic medium. △ Less

Submitted 25 June, 2018; v1 submitted 2 March, 2018; originally announced March 2018.

Comments: 18 pages, 15 figures; version accepted in ApJ

Production of Cool Gas in Thermally-Driven Outflows

Authors: Evan E. Schneider, Brant E. Robertson, Todd A. Thompson

Abstract: Galactic outflows commonly contain multiphase gas, and its physical origin requires explanation. Using the CGOLS (Cholla Galactic OutfLow Simulations) suite of high-resolution isolated galaxy models, we demonstrate the viability of rapid radiative cooling as a source of fast-moving ($v \sim 1000$ km/s), cool ($10^4$ K) gas observed in absorption line studies of outflows around some star-forming ga… ▽ More Galactic outflows commonly contain multiphase gas, and its physical origin requires explanation. Using the CGOLS (Cholla Galactic OutfLow Simulations) suite of high-resolution isolated galaxy models, we demonstrate the viability of rapid radiative cooling as a source of fast-moving ($v \sim 1000$ km/s), cool ($10^4$ K) gas observed in absorption line studies of outflows around some star-forming galaxies. By varying the mass-loading and geometry of the simulated winds, we identify a region of parameter space that leads to cool gas in outflows. In particular, when using an analytically-motivated central feedback model, we find that cooling flows can be produced with reasonable mass-loading rates ($\dot{M}_{wind} / \dot{M}_{SFR} \sim 0.5$), provided the star formation rate surface density is high. When a more realistic clustered feedback model is applied, destruction of high density clouds near the disk and interactions between different outflow regions indicate that lower mass-loading rates of the hot gas within the feedback region may still produce multiphase outflows. These results suggest an origin for fast-moving cool gas in outflows that does not rely on directly accelerating cool gas from the interstellar medium. These cooling flows may additionally provide an explanation for the multiphase gas ubiquitously observed in the halos of star-forming galaxies at low redshift. △ Less

Submitted 7 August, 2018; v1 submitted 2 March, 2018; originally announced March 2018.

Comments: 19 pages, 13 figures; version accepted to ApJ

Hydrodynamical Coupling of Mass and Momentum in Multiphase Galactic Winds

Authors: Evan E. Schneider, Brant E. Robertson

Abstract: Using a set of high resolution hydrodynamical simulations run with the Cholla code, we investigate how mass and momentum couple to the multiphase components of galactic winds. The simulations model the interaction between a hot wind driven by supernova explosions and a cooler, denser cloud of interstellar or circumgalactic media. By resolving scales of $Δx<0.1$ pc over $>100$ pc distances our calc… ▽ More Using a set of high resolution hydrodynamical simulations run with the Cholla code, we investigate how mass and momentum couple to the multiphase components of galactic winds. The simulations model the interaction between a hot wind driven by supernova explosions and a cooler, denser cloud of interstellar or circumgalactic media. By resolving scales of $Δx<0.1$ pc over $>100$ pc distances our calculations capture how the cloud disruption leads to a distribution of densities and temperatures in the resulting multiphase outflow, and quantify the mass and momentum associated with each phase. We find the multiphase wind contains comparable mass and momenta in phases over a wide range of densities and temperatures extending from the hot wind ($n \approx 10^{-2.5}$ $\mathrm{cm}^{-3}$, $T \approx 10^{6.5}$ K) to the coldest components ($n \approx 10^2$ $\mathrm{cm}^{-3}$, $T \approx 10^2$ K). We further find that the momentum distributes roughly in proportion to the mass in each phase, and the mass-loading of the hot phase by the destruction of cold, dense material is an efficient process. These results provide new insight into the physical origin of observed multiphase galactic outflows, and inform galaxy formation models that include coarser treatments of galactic winds. Our results confirm that cool gas observed in outflows at large distances from the galaxy ($\gtrsim1$kpc) likely does not originate through the entrainment of cold material near the central starburst. △ Less

Submitted 22 September, 2016; v1 submitted 6 July, 2016; originally announced July 2016.

Comments: 28 pages, 20 figures. Version accepted by ApJ

Cholla : A New Massively-Parallel Hydrodynamics Code For Astrophysical Simulation

Authors: Evan E. Schneider, Brant E. Robertson

Abstract: We present Cholla (Computational Hydrodynamics On ParaLLel Architectures), a new three-dimensional hydrodynamics code that harnesses the power of graphics processing units (GPUs) to accelerate astrophysical simulations. Cholla models the Euler equations on a static mesh using state-of-the-art techniques, including the unsplit Corner Transport Upwind (CTU) algorithm, a variety of exact and approxim… ▽ More We present Cholla (Computational Hydrodynamics On ParaLLel Architectures), a new three-dimensional hydrodynamics code that harnesses the power of graphics processing units (GPUs) to accelerate astrophysical simulations. Cholla models the Euler equations on a static mesh using state-of-the-art techniques, including the unsplit Corner Transport Upwind (CTU) algorithm, a variety of exact and approximate Riemann solvers, and multiple spatial reconstruction techniques including the piecewise parabolic method (PPM). Using GPUs, Cholla evolves the fluid properties of thousands of cells simultaneously and can update over ten million cells per GPU-second while using an exact Riemann solver and PPM reconstruction. Owing to the massively-parallel architecture of GPUs and the design of the Cholla code, astrophysical simulations with physically interesting grid resolutions (> 256^3) can easily be computed on a single device. We use the Message Passing Interface library to extend calculations onto multiple devices and demonstrate nearly ideal scaling beyond 64 GPUs. A suite of test problems highlights the physical accuracy of our modeling and provides a useful comparison to other codes. We then use Cholla to simulate the interaction of a shock wave with a gas cloud in the interstellar medium, showing that the evolution of the cloud is highly dependent on its density structure. We reconcile the computed mixing time of a turbulent cloud with a realistic density distribution destroyed by a strong shock with the existing analytic theory for spherical cloud destruction by describing the system in terms of its median gas density. △ Less

Submitted 6 March, 2015; v1 submitted 15 October, 2014; originally announced October 2014.

Comments: 38 pages, 20 figures. Version accepted to ApJS

Unravelling the Mysteries of the Leo Ring: An Absorption Line Study of an Unusual Gas Cloud

Authors: J. L. Rosenberg, Karl Haislmaier, M. L. Giroux, B. A. Keeney, S. E. Schneider

Abstract: Since the 1980's discovery of the large (2x10^9 Msun) intergalactic cloud known as the Leo Ring, this object has been the center of a lively debate about its origin. Determining the origin of this object is still important as we develop a deeper understanding of the accretion and feedback processes that shape galaxy evolution. We present HST/COS observations of three sightlines near the Ring, two… ▽ More Since the 1980's discovery of the large (2x10^9 Msun) intergalactic cloud known as the Leo Ring, this object has been the center of a lively debate about its origin. Determining the origin of this object is still important as we develop a deeper understanding of the accretion and feedback processes that shape galaxy evolution. We present HST/COS observations of three sightlines near the Ring, two of which penetrate the high column density neutral hydrogen gas visible in 21 cm observations of the object. These observations provide the first direct measurement of the metallicity of the gas in the Ring, an important clue to its origins. Our best estimate of the metallicity of the ring is ~10% Zsun, higher than expected for primordial gas but lower than expected from an interaction. We discuss possible modifications to the interaction and primordial gas scenarios that would be consistent with this metallicity measurement. △ Less

Submitted 6 June, 2014; originally announced June 2014.

Comments: 11 pages, 7 figures, accepted ApJ

New Constraints on Cosmic Reionization from the 2012 Hubble Ultra Deep Field Campaign

Authors: Brant E. Robertson, Steven R. Furlanetto, Evan Schneider, Stephane Charlot, Richard S. Ellis, Daniel P. Stark, Ross J. McLure, James S. Dunlop, Anton Koekemoer, Matthew A. Schenker, Masami Ouchi, Yoshiaki Ono, Emma Curtis-Lake, Alexander B. Rogers, Rebecca A. A. Bowler, Michele Cirasuolo

Abstract: Understanding cosmic reionization requires the identification and characterization of early sources of hydrogen-ionizing photons. The 2012 Hubble Ultra Deep Field (UDF12) campaign has acquired the deepest infrared images with the Wide Field Camera 3 aboard Hubble Space Telescope and, for the first time, systematically explored the galaxy population deep into the era when cosmic microwave backgroun… ▽ More Understanding cosmic reionization requires the identification and characterization of early sources of hydrogen-ionizing photons. The 2012 Hubble Ultra Deep Field (UDF12) campaign has acquired the deepest infrared images with the Wide Field Camera 3 aboard Hubble Space Telescope and, for the first time, systematically explored the galaxy population deep into the era when cosmic microwave background (CMB) data indicates reionization was underway. The UDF12 campaign thus provides the best constraints to date on the abundance, luminosity distribution, and spectral properties of early star-forming galaxies. We synthesize the new UDF12 results with the most recent constraints from CMB observations to infer redshift-dependent ultraviolet (UV) luminosity densities, reionization histories, and electron scattering optical depth evolution consistent with the available data. Under reasonable assumptions about the escape fraction of hydrogen ionizing photons and the intergalactic medium clumping factor, we find that to fully reionize the universe by redshift z~6 the population of star-forming galaxies at redshifts z~7-9 likely must extend in luminosity below the UDF12 limits to absolute UV magnitudes of M_UV\sim -13 or fainter. Moreover, low levels of star formation extending to redshifts z~15-25, as suggested by the normal UV colors of z\simeq7-8 galaxies and the smooth decline in abundance with redshift observed by UDF12 to z\simeq10, are additionally likely required to reproduce the optical depth to electron scattering inferred from CMB observations. △ Less

Submitted 20 March, 2013; v1 submitted 7 January, 2013; originally announced January 2013.

Comments: Version accepted by ApJ (originally submitted Jan 5, 2013). The UDF12 website can be found at http://udf12.arizona.edu

The UV Luminosity Function of star-forming galaxies via dropout selection at redshifts z ~ 7 and 8 from the 2012 Ultra Deep Field campaign

Authors: Matthew A. Schenker, Brant E. Robertson, Richard S. Ellis, Yoshiaki Ono, Ross J. McLure, James S. Dunlop, Anton Koekemoer, Rebecca A. A. Bowler, Masami Ouchi, Emma Curtis-Lake, Alexander B. Rogers, Evan Schneider, Stephane Charlot, Daniel P. Stark, Steven R. Furlanetto, Michele Cirasuolo

Abstract: We present a catalog of high redshift star-forming galaxies selected to lie within the redshift range z ~ 7-8 using the Ultra Deep Field 2012 (UDF12), the deepest near-infrared (near-IR) exposures yet taken with the Hubble Space Telescope. As a result of the increased near-infrared exposure time compared to previous HST imaging in this field, we probe 0.65 (0.25) mag fainter in absolute UV magnitu… ▽ More We present a catalog of high redshift star-forming galaxies selected to lie within the redshift range z ~ 7-8 using the Ultra Deep Field 2012 (UDF12), the deepest near-infrared (near-IR) exposures yet taken with the Hubble Space Telescope. As a result of the increased near-infrared exposure time compared to previous HST imaging in this field, we probe 0.65 (0.25) mag fainter in absolute UV magnitude, at z ~ 7 (8), which increases confidence in a measurement of the faint end slope of the galaxy luminosity function. Through a 0.7 mag deeper limit in the key F105W filter that encompasses or lies just longward of the Lyman break, we also achieve a much-refined color-color selection that balances high redshift completeness and a low expected contamination fraction. We improve the number of drop-out selected UDF sources to 47 at z ~ 7 and 27 at z ~ 8. Incorporating brighter archival and ground-based samples, we measure the z ~ 7 UV luminosity function to an absolute magnitude limit of M_UV = -17 and find a faint end Schechter slope of α= -1.87+/- 0.18. Using a similar color-color selection at z ~ 8 that takes account of our newly-added imaging in the F140W filter, and incorporating archival data from the HIPPIES and BoRG campaigns, we provide a robust estimate of the faint end slope at z ~ 8, α= -1.94 +/- 0.23. We briefly discuss our results in the context of earlier work and that derived using the same UDF12 data but with an independent photometric redshift technique (McLure et al 2012). △ Less

Submitted 20 December, 2012; v1 submitted 19 December, 2012; originally announced December 2012.

Comments: 15 pages, 5 figures, To be submitted to ApJ

Evolution of the Sizes of Galaxies over 7<z<12 Revealed by the 2012 Hubble Ultra Deep Field Campaign

Authors: Yoshiaki Ono, Masami Ouchi, Emma Curtis-Lake, Matthew A. Schenker, Richard S. Ellis, Ross J. McLure, James S. Dunlop, Brant E. Robertson, Anton M. Koekemoer, Rebecca A. A. Bowler, Alexander B. Rogers, Evan Schneider, Stephane Charlot, Daniel P. Stark, Kazuhiro Shimasaku, Steven R. Furlanetto, Michele Cirasuolo

Abstract: We analyze the redshift- and luminosity-dependent sizes of dropout galaxy candidates in the redshift range z~7-12 using deep images from the UDF12 campaign, data which offers two distinct advantages over that used in earlier work. Firstly, we utilize the increased S/N ratio offered by the UDF12 imaging to provide improved size measurements for known galaxies at z=6.5-8 in the HUDF. Specifically, w… ▽ More We analyze the redshift- and luminosity-dependent sizes of dropout galaxy candidates in the redshift range z~7-12 using deep images from the UDF12 campaign, data which offers two distinct advantages over that used in earlier work. Firstly, we utilize the increased S/N ratio offered by the UDF12 imaging to provide improved size measurements for known galaxies at z=6.5-8 in the HUDF. Specifically, we stack the new deep F140W image with the existing F125W data in order to provide improved measurements of the half-light radii of z-dropouts. Similarly we stack this image with the new deep UDF12 F160W image to obtain new size measurements for a sample of Y-dropouts. Secondly, because the UDF12 data have allowed the construction of the first robust galaxy sample in the HUDF at z>8, we have been able to extend the measurement of average galaxy size out to significantly higher redshifts. Restricting our size measurements to sources which are now detected at >15sigma, we confirm earlier indications that the average half-light radii of z~7-12 galaxies are extremely small, 0.3-0.4 kpc, comparable to the sizes of giant molecular associations in local star-forming galaxies. We also confirm that there is a clear trend of decreasing half-light radius with increasing redshift, and provide the first evidence that this trend continues beyond z~8. Modeling the evolution of the average half-light radius as a power-law (1+z)^s, we obtain a best-fit index of s=-1.28+/-0.13 over the redshift range z~4-12, mid-way between the physically expected evolution for baryons embedded in dark halos of constant mass (s=-1) and constant velocity (s=-1.5). A clear size-luminosity relation, such as that found at lower redshift, is also evident in both our z- and Y-dropout sample. This relation can be interpreted in terms of a constant surface density of star formation over a range in luminosity of 0.05-1.0L*_z=3.(abridged) △ Less

Submitted 16 December, 2012; originally announced December 2012.

Comments: 15 pages, 13 figures, submitted to ApJ

The 2012 Hubble Ultra Deep Field (UDF12): Observational Overview

Authors: Anton M. Koekemoer, Richard S Ellis, Ross J. McLure, James S. Dunlop, Brant E Robertson, Yoshiaki Ono, Matthew A. Schenker, Masami Ouchi, Rebecca A. A. Bowler, Alexander B. Rogers, Emma Curtis-Lake, Evan Schneider, Stephane Charlot, Daniel P. Stark, Steven R. Furlanetto, Michele Cirasuolo, V. Wild, T. Targett

Abstract: We present the 2012 Hubble Ultra Deep Field campaign (UDF12), a large 128-orbit Cycle 19 \HST\ program aimed at extending previous WFC3/IR observations of the UDF by quadrupling the exposure time in the F105W filter, imaging in an additional F140W filter, and extending the F160W exposure time by 50%. The principal scientific goal of this project is to determine whether galaxies reionized the unive… ▽ More We present the 2012 Hubble Ultra Deep Field campaign (UDF12), a large 128-orbit Cycle 19 \HST\ program aimed at extending previous WFC3/IR observations of the UDF by quadrupling the exposure time in the F105W filter, imaging in an additional F140W filter, and extending the F160W exposure time by 50%. The principal scientific goal of this project is to determine whether galaxies reionized the universe; our observations are designed to provide a robust determination of the star formation density at $z$$\,\gtrsim\,$8, improve measurements of the ultraviolet continuum slope at $z$$\,\sim\,7\,-\,$8, facilitate the construction of new samples of $z$$\,\sim\,9\,-\,$10 candidates, and enable the detection of sources up to $z$$\,\sim\,$12. For this project we committed to combining these and other WFC3/IR imaging observations of the UDF area into a single homogeneous dataset, to provide the deepest near-infrared observations of the sky currently achievable. In this paper we present the observational overview of the project, motivated by its scientific goals, and describe the procedures used in reducing the data as well as the final products that are produced. We have used the most up up-to-date methods for calibrating and combining the images, in particular paying attention to correcting several instrumental effects. We release the full combined mosaics, comprising a single, unified set of mosaics of the UDF, providing the deepest near-infrared blank-field view of the universe obtained to date, reaching magnitudes as deep as AB$\,\sim\,$30 in the near-infrared, and yielding a legacy dataset on this field of lasting scientific value to the community. △ Less

Submitted 6 December, 2012; originally announced December 2012.

Comments: 16 pages, 9 figures, submitted to ApJS

The UV continua and inferred stellar populations of galaxies at z ~ 7 - 9 revealed by the Hubble Ultra Deep Field 2012 campaign

Authors: J. S. Dunlop, A. B. Rogers, R. J. McLure, R. S. Ellis, B. E. Robertson, A. Koekemoer, P. Dayal, E. Curtis-Lake, V. Wild, S. Charlot, R. A. A. Bowler, M. A. Schenker, M. Ouchi, Y. Ono, M. Cirasuolo, S. R. Furlanetto, D. P. Stark, T. A. Targett, E. Schneider

Abstract: We use the new ultra-deep, near-infrared imaging of the Hubble Ultra-Deep Field (HUDF) provided by our UDF12 HST WFC3/IR campaign to explore the rest-frame UV properties of galaxies at redshifts z > 6.5. We present the first unbiased measurement of the average UV power-law index, beta, for faint galaxies at z ~ 7, the first meaningful measurements of beta at z ~ 8, and tentative estimates for a ne… ▽ More We use the new ultra-deep, near-infrared imaging of the Hubble Ultra-Deep Field (HUDF) provided by our UDF12 HST WFC3/IR campaign to explore the rest-frame UV properties of galaxies at redshifts z > 6.5. We present the first unbiased measurement of the average UV power-law index, beta, for faint galaxies at z ~ 7, the first meaningful measurements of beta at z ~ 8, and tentative estimates for a new sample of galaxies at z ~ 9. Utilising galaxy selection in the new F140W imaging to minimize colour bias, and applying both colour and power-law estimators of beta, we find beta = -2.1 (+/-0.2) at z ~ 7 for galaxies with M_UV ~ -18. This means that the faintest galaxies uncovered at this epoch have, on average, UV colours no more extreme than those displayed by the bluest star-forming galaxies at low redshift. At z ~ 8 we find a similar value, beta = -1.9 (+/-0.3). At z ~ 9, we find beta = -1.8 (+/-0.6), essentially unchanged from z ~ 6 - 7 (albeit highly uncertain). Finally, we show that there is as yet no evidence for a significant intrinsic scatter in beta within our new, robust z ~ 7 galaxy sample. Our results are most easily explained by a population of steadily star-forming galaxies with either ~ solar metallicity and zero dust, or moderately sub-solar (~ 10-20%) metallicity with modest dust obscuration (A_V ~ 0.1-0.2). This latter interpretation is consistent with the predictions of a state-of-the-art galaxy-formation simulation, which also suggests that a significant population of very-low metallicity, dust-free galaxies with beta ~ -2.5 may not emerge until M_UV > -16, a regime likely to remain inaccessible until the James Webb Space Telescope. △ Less

Submitted 22 April, 2013; v1 submitted 4 December, 2012; originally announced December 2012.

Comments: 15 pages, 7 figures, updated to match accepted version, now including an appendix tabulating the individual beta values for the objects used in this study

The Abundance of Star-Forming Galaxies in the Redshift Range 8.5 to 12: New Results from the 2012 Hubble Ultra Deep Field Campaign

Authors: Richard S Ellis, Ross J McLure, James S Dunlop, Brant E Robertson, Yoshiaki Ono, Matthew A Schenker, Anton Koekemoer, Rebecca A A Bowler, Masami Ouchi, Alexander B Rogers, Emma Curtis-Lake, Evan Schneider, Stephane Charlot, Daniel P Stark, Steven R Furlanetto, Michele Cirasuolo

Abstract: We present the results of the deepest search to date for star-forming galaxies beyond a redshift z~8.5 utilizing a new sequence of near-infrared Wide Field Camera 3 images of the Hubble Ultra Deep Field. This `UDF12' campaign completed in September 2012 doubles the earlier exposures with WFC3/IR in this field and quadruples the exposure in the key F105W filter used to locate such distant galaxies.… ▽ More We present the results of the deepest search to date for star-forming galaxies beyond a redshift z~8.5 utilizing a new sequence of near-infrared Wide Field Camera 3 images of the Hubble Ultra Deep Field. This `UDF12' campaign completed in September 2012 doubles the earlier exposures with WFC3/IR in this field and quadruples the exposure in the key F105W filter used to locate such distant galaxies. Combined with additional imaging in the F140W filter, the fidelity of high redshift candidates is greatly improved. Using spectral energy distribution fitting techniques on objects selected from a deep multi-band near-infrared stack we find 7 promising z>8.5 candidates. As none of the previously claimed UDF candidates with 8.5<z<10 is confirmed by our deeper multi-band imaging, our campaign has transformed the measured abundance of galaxies in this redshift range. Although we recover the candidate UDFj-39546284 (previously proposed at z=10.3), it is undetected in the newly added F140W image, implying it lies at z=11.9 or is an intense emission line galaxy at z~2.4. Although no physically-plausible model can explain the required line intensity given the lack of Lyman alpha or broad-band UV signal, without an infrared spectrum we cannot rule out an exotic interloper. Regardless, our robust z ~ 8.5 - 10 sample demonstrates a luminosity density that continues the smooth decline observed over 6 < z < 8. Such continuity has important implications for models of cosmic reionization and future searches for z>10 galaxies with JWST. △ Less

Submitted 28 November, 2012; originally announced November 2012.

Comments: Submitted to Astrophys. J. Lett. on Nov. 7, 2012. Resubmitted after responding to referee report. UDF12 public webpage available at http://udf12.arizona.edu

TW Hya: Spectral Variability, X-Rays, and Accretion Diagnostics

Authors: A. K. Dupree, N. S. Brickhouse, S. R. Cranmer, G. J. M. Luna, E. E. Schneider, M. S. Bessell, A. Bonanos, L. A. Crause, W. A. Lawson, S. V. Mallik, S. C. Schuler

Abstract: The nearest accreting T Tauri star, TW Hya was observed with spectroscopic and photometric measurements simultaneous with a long se gmented exposure using the CHANDRA satellite. Contemporaneous optical photometry from WASP-S indicates a 4.74 day period was present during this time. Absence of a similar periodicity in the H-alpha flux and the total X-ray flux points to a different source of photome… ▽ More The nearest accreting T Tauri star, TW Hya was observed with spectroscopic and photometric measurements simultaneous with a long se gmented exposure using the CHANDRA satellite. Contemporaneous optical photometry from WASP-S indicates a 4.74 day period was present during this time. Absence of a similar periodicity in the H-alpha flux and the total X-ray flux points to a different source of photometric variations. The H-alpha emission line appears intrinsically broad and symmetric, and both the profile and its variability suggest an origin in the post-shock cooling region. An accretion event, signaled by soft X-rays, is traced spectroscopically for the first time through the optical emission line profiles. After the accretion event, downflowing turbulent material observed in the H-alpha and H-beta lines is followed by He I (5876A) broadening. Optical veiling increases with a delay of about 2 hours after the X-ray accretion event. The response of the stellar coronal emission to an increase in the veiling follows about 2.4 hours later, giving direct evidence that the stellar corona is heated in part by accretion. Subsequently, the stellar wind becomes re-established. We suggest a model that incorporates this sequential series of events: an accretion shock, a cooling downflow in a supersonically turbulent region, followed by photospheric and later, coronal heating. This model naturally explains the presence of broad optical and ultraviolet lines, and affects the mass accretion rates determined from emission line profiles. △ Less

Submitted 28 February, 2012; originally announced February 2012.

Comments: 61 pages; 22 figures; to appear in The Astrophysical Journal

HI line observations of 2MASS galaxies in the Zone of Avoidance

Authors: W. van Driel, S. E. Schneider, R. C. Kraan-Korteweg, D. Monnier Ragaigne

Abstract: A pilot survey has been made to obtain 21cm HI emission line profiles for 197 objects in the Zone of Avoidance (ZoA) that were classified as galaxies in the 2MASS all-sky near-infrared Extended Source Catalog (2MASX). 116 of the sources were observed using the Nancay radio telescope, in the 325 to 11,825 km/s range, and the other 81 sources were observed at Arecibo in the -500 to 11,000 km/s ran… ▽ More A pilot survey has been made to obtain 21cm HI emission line profiles for 197 objects in the Zone of Avoidance (ZoA) that were classified as galaxies in the 2MASS all-sky near-infrared Extended Source Catalog (2MASX). 116 of the sources were observed using the Nancay radio telescope, in the 325 to 11,825 km/s range, and the other 81 sources were observed at Arecibo in the -500 to 11,000 km/s range, and for 9 also in the 9,500 to 21,000 km/s range. Global HI line parameters are presented for the 22 and 29 2MASX objects that were detected at Nancay and Arecibo, respectively, as well as upper limits for the undetected 2MASX objects. Whereas object 2MASX J08170147-3410277 appears to be a very massive galaxy with an HI mass of 4.6 x 10E10 Msun, it is clear that only radio synthesis HI imaging observations will allow a firm conclusion on this. Overall, the global properties of the detected galaxies match those of other ZoA HI surveys. Although the detections are as yet too sparse to give further insight into suspected or unknown large-scale structures in the ZoA, they already indicate that an extension of the present pilot survey is bound to quantify filaments, clusters, and voids behind this part of the Milky Way. It is shown that the number of candidate 2MASS-selected ZoA galaxies to be observed in HI could have been reduced by about 15% through examination of composite near-infrared images and the application of extinction-corrected near-infrared colour limits. Present results confirm that the Galactic extinction values from Schlegel et al. (1998) are valid for latitudes above |b|=5 degrees, but increasingly less so for lower latitudes. △ Less

Submitted 11 May, 2009; originally announced May 2009.

Comments: 22 pages, 10 figures

The Arecibo Galaxy Environments Survey - Description and Early Results

Authors: R. F. Minchin, R. Auld, J. I. Davies, B. Catinella, L. Cortese, S. Linder, E. Momjian, E. Muller, K. O'Neil, J. Rosenberg, S Sabatini, S. E. Schneider, M. Stage, W. van Driel

Abstract: The Arecibo Galaxy Environments Survey (AGES) is a 2000-hour neutral hydrogen (HI) survey using the new Arecibo L-band Feed Array (ALFA) multibeam instrument at Arecibo Observatory. It will cover 200 square degrees of sky, sampling a range of environments from the Local Void through to the Virgo Cluster with higher sensitivity, spatial resolution and velocity resolution than previous neutral hyd… ▽ More The Arecibo Galaxy Environments Survey (AGES) is a 2000-hour neutral hydrogen (HI) survey using the new Arecibo L-band Feed Array (ALFA) multibeam instrument at Arecibo Observatory. It will cover 200 square degrees of sky, sampling a range of environments from the Local Void through to the Virgo Cluster with higher sensitivity, spatial resolution and velocity resolution than previous neutral hydrogen surveys. △ Less

Submitted 20 October, 2006; originally announced October 2006.

Comments: IAU Symposium 235 - Galaxy Evolution across the Hubble Time. 3 pages (only 1 st to appear in printed proceedings), 3 figures. Uses iaus.cls. Full resolution PDF is available from http://www2.naic.edu/alfa/ealfa/AGES-IAUS235-Minchin.pdf and a PDF of the poster is available from http://www2.naic.edu/alfa/ealfa/AGES-IAU-Minchin.pdf

Galaxy transmutations: The double ringed galaxy ESO 474-G26

Authors: V. Reshetnikov, F. Bournaud, F. Combes, M. Faundez-Abans, M. de Oliveira-Abans, W. van Driel, S. E. Schneider

Abstract: Surface photometry and a 21cm HI line spectrum of the giant double-ringed galaxy ESO 474-G26 are presented. The morphology of this system is unique among the 30,000 galaxies with >B15. Two almost orthogonal optical rings with diameters of 60 and 40 kpc surround the central body (assuming H0=70 km/s/Mpc). The outer one is an equatorial ring, while the inner ring lies in a nearly polar plane. The… ▽ More Surface photometry and a 21cm HI line spectrum of the giant double-ringed galaxy ESO 474-G26 are presented. The morphology of this system is unique among the 30,000 galaxies with >B15. Two almost orthogonal optical rings with diameters of 60 and 40 kpc surround the central body (assuming H0=70 km/s/Mpc). The outer one is an equatorial ring, while the inner ring lies in a nearly polar plane. The rings have blue optical colors typical of late-type spirals. Both appear to be rotating around the central galaxy, so that this system can be considered as a kinematically confirmed polar ring galaxy. Its observational characteristics are typical of galaxy merger remnants. Although the central object has a surface brightness distribution typical of elliptical galaxies, it has a higher surface brightness for its effective radius than ordinary ellipticals. Possible origins of this galaxy are discussed and numerical simulations are presented that illustrate the formation of the two rings in the merging process of two spiral galaxies, in which the observed appearance of ESO 474-G26 appears to be a transient stage. △ Less

Submitted 1 December, 2004; originally announced December 2004.

Comments: Accepted for publication in A&A

Gas and Stars in an HI-Selected Galaxy Sample

Authors: J. L. Rosenberg, S. E. Schneider, J. Posson-Brown

Abstract: We present the results of a J-band study of the HI-selected Arecibo Dual-Beam Survey and Arecibo Slice Survey galaxy samples using the 2 Micron All-Sky Survey data. We find that these galaxies span a wide range of stellar and gas properties. However, despite the diversity within the samples, we find a very tight correlation between luminosity and size in the J-band, similar to that we previously… ▽ More We present the results of a J-band study of the HI-selected Arecibo Dual-Beam Survey and Arecibo Slice Survey galaxy samples using the 2 Micron All-Sky Survey data. We find that these galaxies span a wide range of stellar and gas properties. However, despite the diversity within the samples, we find a very tight correlation between luminosity and size in the J-band, similar to that we previously found (Rosenberg & Schneider 2003) between the HI mass and size. We also find that the correlation between the baryonic mass and the J-band diameter is even tighter than between the baryonic mass and the rotational velocity. △ Less

Submitted 17 November, 2004; originally announced November 2004.

Comments: AJ in press, 17 pages (including tables and figures) + 6 additional jpg figures

The Molecular Gas Distribution and Schmidt Law in M33

Authors: M. H. Heyer, E. Corbelli, S. E. Schneider, J. S. Young

Abstract: The relationship between the star formation rate and surface density of neutral gas within the disk of M33 is examined with new imaging observations of CO J=1-0 emission gathered with the FCRAO 14m telescope and IRAS HiRes images of the 60 micron and 100 micron emission. The Schmidt law, Sigma_SFR ~ Sigma_gas^n, is constructed using radial profiles of the HI 21cm, CO, and far infrared emission.… ▽ More The relationship between the star formation rate and surface density of neutral gas within the disk of M33 is examined with new imaging observations of CO J=1-0 emission gathered with the FCRAO 14m telescope and IRAS HiRes images of the 60 micron and 100 micron emission. The Schmidt law, Sigma_SFR ~ Sigma_gas^n, is constructed using radial profiles of the HI 21cm, CO, and far infrared emission. A strong correlation is identified between the star formation rate and molecular gas surface density. This suggests that the condensation of giant molecular clouds is the limiting step to star formation within the M33 disk. The corresponding molecular Schmidt index, n_{mol}, is 1.36 +/- 0.08. The star formation rate has a steep dependence on total mass gas surface density, (Sigma_{HI}+Sigma_{H_2}), owing to the shallow radial profile of the atomic gas which dominates the total gas surface density for most radii. The disk pressure of the gas is shown to play a prominent role in regulating the molecular gas fraction in M33. △ Less

Submitted 10 November, 2003; originally announced November 2003.

Comments: 19 pages + 5 figures. Accepted for publication in ApJ

Journal ref: Astrophys.J. 602 (2004) 723-729

A search for Low Surface Brightness galaxies in the near-infrared III. Nancay HI line observation

Authors: D. Monnier Ragaigne, W. van Driel, S. E. Schneider, C. Balkowski, T. H. Jarrett

Abstract: A total of 334 Low Surface Brightness galaxies detected in the 2MASS all-sky near-infrared survey have been observed in the 21 cm HI line using the Nancay telescope. All have a K_s-band mean central surface brightness, measured within a 5 arcsec radius, fainter than 18 mag/arcsec^2 and a K_s-band isophotal radius at the 20 mag/arcsec^2 level larger than 20 arcsec. We present global HI line param… ▽ More A total of 334 Low Surface Brightness galaxies detected in the 2MASS all-sky near-infrared survey have been observed in the 21 cm HI line using the Nancay telescope. All have a K_s-band mean central surface brightness, measured within a 5 arcsec radius, fainter than 18 mag/arcsec^2 and a K_s-band isophotal radius at the 20 mag/arcsec^2 level larger than 20 arcsec. We present global HI line parameters for the 171 clearly detected objects and the 23 marginal detections, as well as upper limits for the undetected objects. The 171 clear detections comprise 50 previously uncatalogued objects and 41 objects with a PGC entry only. △ Less

Submitted 17 May, 2003; originally announced May 2003.

Comments: Accepted for publication by Astronomy and Astrophysics on 15 May 2003, 41 pages

Journal ref: Astron.Astrophys. 408 (2003) 465-478

A search for Low Surface Brightness galaxies in the near-infrared II. Arecibo HI line observations

Authors: D. Monnier Ragaigne, W. van Driel, K. O'Neil, S. E. Schneider, C. Balkowski, T. H. Jarrett

Abstract: A total of 367 Low Surface Brightness galaxies detected in the 2MASS all-sky near-infrared survey have been observed in the 21 cm HI line using the Arecibo telescope. All have a K_s-band mean central surface brightness, measured within a 5 arcsec radius, fainter than 18 mag/arcsec^(2). We present global HI line parameters for the 107 clearly detected objects and the 21 marginal detections, as we… ▽ More A total of 367 Low Surface Brightness galaxies detected in the 2MASS all-sky near-infrared survey have been observed in the 21 cm HI line using the Arecibo telescope. All have a K_s-band mean central surface brightness, measured within a 5 arcsec radius, fainter than 18 mag/arcsec^(2). We present global HI line parameters for the 107 clearly detected objects and the 21 marginal detections, as well as upper limits for the undetected objects. The 107 clear detections comprise 15 previously uncatalogued objects and 36 with a PGC entry only. △ Less

Submitted 30 April, 2003; originally announced April 2003.

Comments: Accepted for publication by Astronomy and Astrophysics on 30 April 2003, 34 pages

Journal ref: Astron.Astrophys. 408 (2003) 67-78

A search for Low Surface Brightness galaxies in the near-infrared I. Selection of the sample

Authors: D. Monnier Ragaigne, W. van Driel, S. E. Schneider, T. H. Jarrett, C. Balkowski

Abstract: A sample of about 3,800 Low Surface Brightness (LSB) galaxies was selected using the all-sky near-infrared (J, H and K_s-band) 2MASS survey. The selected objects have a mean central surface brightness within a 5 arcsec radius around their centre fainter than 18 mag/sq.arcsec in the K_s band, making them the lowest surface brightness galaxies detected by 2MASS. A description is given of the relev… ▽ More A sample of about 3,800 Low Surface Brightness (LSB) galaxies was selected using the all-sky near-infrared (J, H and K_s-band) 2MASS survey. The selected objects have a mean central surface brightness within a 5 arcsec radius around their centre fainter than 18 mag/sq.arcsec in the K_s band, making them the lowest surface brightness galaxies detected by 2MASS. A description is given of the relevant properties of the 2MASS survey and the LSB galaxy selection procedure, as well as of basic photometric properties of the selected objects. The latter properties are compared to those of other samples of galaxies, of both LSBs and `classical' high surface brightness (HSB) objects, which were selected in the optical. The 2MASS LSBs have a (B_T_c)-(K_T) colour which is on average 0.9 mag bluer than that of HSBs from the NGC. The 2MASS sample does not appear to contain a significant population of red objects. △ Less

Submitted 22 April, 2003; originally announced April 2003.

Comments: accepted for publication in Astronomy and Astrophysics on 24/2/2003; 62 pages

Journal ref: Astron.Astrophys. 405 (2003) 99-110

Clusters of galaxies in the local universe

Authors: C. S. Kochanek, Martin White, J. Huchra, L. Macri, T. H. Jarrett, S. E. Schneider, J. Mader

Abstract: We use a matched filter algorithm to find and study clusters in both N-body simulations artificially populated with galaxies and the 2MASS survey. In addition to numerous checks of the matched filter algorithm, we present results on the halo multiplicity function and the cluster number function. For a subset of our identified clusters we have information on X-ray temperatures and luminosities wh… ▽ More We use a matched filter algorithm to find and study clusters in both N-body simulations artificially populated with galaxies and the 2MASS survey. In addition to numerous checks of the matched filter algorithm, we present results on the halo multiplicity function and the cluster number function. For a subset of our identified clusters we have information on X-ray temperatures and luminosities which we cross-correlate with optical richness and galaxy velocity dispersions. With all quantities normalized by the spherical radius corresponding to a mass overdensity of Delta_M=200 or the equivalent galaxy number overdensity of Delta_N=200Omega^{-1}~666, we find that the number of L>L_* galaxies in a cluster of mass M_{200} is log N_{*666} = (1.44 +/- 0.17)+(1.10 +/- 0.09)log(M_{200}h/10^{15}Msun) where the uncertainties are dominated by the scatter created by three choices for relating the observed quantities to the cluster mass. The region inside the virial radius has a K-band cluster mass-to-light ratio of (M/L)_K=(116 +/- 46)h which is essentially independent of cluster mass. Integrating over all clusters more massive than M_{200}=10^{14}Msun/h, the virialized regions of clusters contain ~7% of the local stellar luminosity, quite comparable to the mass fraction in such objects in currently popular LambdaCDM models. △ Less

Submitted 7 August, 2002; originally announced August 2002.

Comments: 22 pages, 15 figures

Journal ref: Astrophys.J.585:161-181,2003

The Contribution of HI-rich Galaxies to the Damped Lyman-alpha Absorber Population at z=0

Authors: Jessica L. Rosenberg, Stephen E. Schneider

Abstract: We present a study of the expected properties of the low redshift damped Lyman-alpha absorber population determined from a sample of HI-selected galaxies in the local universe. Because of a tight correlation between the HI mass and HI cross-section, which we demonstrate spans all galaxy types, we can use our HI-selected sample to predict the properties of the absorption line systems. We use meas… ▽ More We present a study of the expected properties of the low redshift damped Lyman-alpha absorber population determined from a sample of HI-selected galaxies in the local universe. Because of a tight correlation between the HI mass and HI cross-section, which we demonstrate spans all galaxy types, we can use our HI-selected sample to predict the properties of the absorption line systems. We use measurements of the number density and HI cross-section of galaxies to show that the total HI cross-section at column densities sufficient to produce damped Lyman-alpha absorption is consistent with no evolution of the absorber population. We also find that the dN/dz distribution is dominated by galaxies with HI masses near 10^9 Msolar. However, because of the large dispersion in the correlation between HI mass and stellar luminosity, we find that the distribution of dN/dz as a function of L_J is fairly flat. Additionally, we examine the line widths of the HI-selected galaxies and show that there may be evolution in the kinematics of HI-rich galaxies, but it is not necessary for the higher redshift population to contain a greater proportion of high mass galaxies than we find locally. △ Less

Submitted 4 November, 2002; v1 submitted 11 February, 2002; originally announced February 2002.

Comments: 14 pages, accepted ApJ - Includes editorial changes and a new appendix on the HI mass/HI size correlation

Journal ref: Astrophys.J. 585 (2003) 256-267

The Arecibo Dual-Beam Survey: The HI Mass Function of Galaxies

Authors: Jessica L. Rosenberg, Stephen E. Schneider

Abstract: We use the HI-selected galaxy sample from the Arecibo Dual-Beam Survey (Rosenberg & Schneider 2000) to determine the shape of the HI mass function of galaxies in the local universe using both the step-wise maximum likelihood and the 1/V_tot methods. Our survey region spanned all 24 hours of right ascension at selected declinations between 8 and 29 degrees covering ~430 deg^2 of sky in the main b… ▽ More We use the HI-selected galaxy sample from the Arecibo Dual-Beam Survey (Rosenberg & Schneider 2000) to determine the shape of the HI mass function of galaxies in the local universe using both the step-wise maximum likelihood and the 1/V_tot methods. Our survey region spanned all 24 hours of right ascension at selected declinations between 8 and 29 degrees covering ~430 deg^2 of sky in the main beam. The survey is not as deep as some previous Arecibo surveys, but it has a larger total search volume and samples a much larger area of the sky. We conducted extensive tests on all aspects of the galaxy detection process, allowing us to empirically correct for our sensitivity limits, unlike the previous surveys. The mass function for the entire sample is quite steep, with a power-law slope of α~ -1.5. We find indications that the slope of the HI mass function is flatter near the Virgo cluster, suggesting that evolutionary effects in high density environments may alter the shape of the HI mass function. These evolutionary effects may help to explain differences in the HI mass function derived by different groups. We are sensitive to the most massive sources (log M > 5x10^10 M\solar) over most of the declination range, \~1 sr, and do not detect any massive low surface brightness galaxies. These statistics restrict the population of Malin 1-like galaxies to <5.5x10^-6 Mpc^-3. △ Less

Submitted 31 October, 2001; originally announced November 2001.

Comments: ApJ accepted, 12 pages

The Contribution of HI-Rich Galaxies to the Damped Absorber Population at z=0

Authors: Jessica L. Rosenberg, Stephen E. Schneider

Abstract: We present a study of HI-rich galaxies in the local universe selected from blind emission-line surveys. These galaxies represent the emission-line counterparts of local damped Lyman-alpha systems. We find that the HI cross-section of galaxies is drawn from a large range of galaxy masses below M_star, 66% of the area comes from galaxies in the range 8.5 < Log M_star < 9.7. Both because of the low… ▽ More We present a study of HI-rich galaxies in the local universe selected from blind emission-line surveys. These galaxies represent the emission-line counterparts of local damped Lyman-alpha systems. We find that the HI cross-section of galaxies is drawn from a large range of galaxy masses below M_star, 66% of the area comes from galaxies in the range 8.5 < Log M_star < 9.7. Both because of the low mass galaxy contribution, and because of the range of galaxy types and luminosities at any given HI mass, the galaxies contributing to the HI cross-section are not exclusively L_star spirals, as is often expected. The optical and near infrared counterparts of these galaxies cover a range of types (from spirals to irregulars), luminosities (from L_star to <0.01 L_star), and surface brightnesses. The range of optical and near infrared properties as well as the kinematics for this population are consistent with the properties for the low-z damped Lyman-alpha absorbers. We also show that the number of HI-rich galaxies in the local universe does not preclude evolution of the low-z damped absorber population, but it is consistent with no evolution. △ Less

Submitted 26 July, 2001; originally announced July 2001.

Comments: 10 pages, 7 figures. To appear in "Extragalactic Gas at Low Redshift" (ASP Conf. Series, Weymann Conf.)

The K-Band Galaxy Luminosity Function

Authors: C. S. Kochanek, M. A. Pahre, E. E. Falco, J. P. Huchra, J. Mader, T. H Jarrett, T. Chester, R. Cutri, S. E. Schneider

Abstract: We measured the K-band luminosity function using a complete sample of 4192 morphologically-typed 2MASS galaxies with 7 < K < 11.25 mag spread over 2.12 str. Early-type (T < -0.5) and late-type (T > -0.5) galaxies have similarly shaped luminosity functions, alpha_e=-0.92+/-0.10 and alpha_l=-0.87+/-0.09. The early-type galaxies are brighter, M_*e=-23.53+/-0.06 mag compared to M_*l=-22.98\pm0.06 ma… ▽ More We measured the K-band luminosity function using a complete sample of 4192 morphologically-typed 2MASS galaxies with 7 < K < 11.25 mag spread over 2.12 str. Early-type (T < -0.5) and late-type (T > -0.5) galaxies have similarly shaped luminosity functions, alpha_e=-0.92+/-0.10 and alpha_l=-0.87+/-0.09. The early-type galaxies are brighter, M_*e=-23.53+/-0.06 mag compared to M_*l=-22.98\pm0.06 mag, but less numerous, n_*e=(0.0045+/-0.0006)h^3/Mpc^3 compared to n_*l=(0.0101+/-0.0013)h^3/Mpc^3 for H_0=100h km/s Mpc, such that the late-type galaxies slightly dominate the K-band luminosity density, j_late/j_early=1.17+/-0.12. Our morphological classifications are internally consistent, consistent with previous classifications and lead to luminosity functions unaffected by the estimated uncertainties in the classifications. These luminosity functions accurately predict the K-band number counts and redshift distributions for K < 18 mag, beyond which the results depend on galaxy evolution and merger histories. △ Less

Submitted 24 November, 2000; originally announced November 2000.

Comments: submitted to ApJ, 25 pages, 6 figures, complete redshift survey. Table 1 included in source

Journal ref: Astrophys.J. 560 (2001) 566-579

The Arecibo Dual-Beam Survey: Filling the Gap in the Extragalactic Census

Authors: J. L. Rosenberg, S. E. Schneider

Abstract: We present the Arecibo Dual-Beam Survey (ADBS), a "blind" 21cm survey for galaxies. The ADBS is the largest to date with the Arecibo radio telescope. We were able to achieve an rms sensitivity of 3-4 mJy in only 7 seconds. The ADBS covers substantially more volume than the Spitzak & Schneider (1998) and the Zwaan et al. (1997) surveys. It is surpassed by the Parkes HIPASS survey (Kilborn et al.… ▽ More We present the Arecibo Dual-Beam Survey (ADBS), a "blind" 21cm survey for galaxies. The ADBS is the largest to date with the Arecibo radio telescope. We were able to achieve an rms sensitivity of 3-4 mJy in only 7 seconds. The ADBS covers substantially more volume than the Spitzak & Schneider (1998) and the Zwaan et al. (1997) surveys. It is surpassed by the Parkes HIPASS survey (Kilborn et al. 1999) but has much better angular resolution and sensitivity. We will discuss the HI sample and our efforts to quantify the completeness and reliability. This is the first survey to use synthetic sources to establish the relation between the flux, line width, and completeness of the survey. An empirical understanding of these quantities is vital to the determination of the mass function. Even if we do not understand all of the reasons why we do not detect galaxies, the empirical understanding allows us to account for them properly. In addition to the more general discussion of the ADBS survey, we will discuss galaxies which have large ratios of HI mass to optical area and probably have extended HI disks. We detect a much higher percentage of these high HI mass to optical area galaxies than are seen in optical surveys indicating that optical surveys miss these galaxies. Among these galaxies are also 7 galaxies with HI masses < 10^8 M_sun. Based on the mass function of Schneider & Rosenberg (this volume), the low mass galaxies are an important part of the general galaxy population. △ Less

Submitted 18 October, 2000; originally announced October 2000.

Comments: 7 pages, presented at Mapping the Hidden Universe: The Universe in HI. eds raan-Korteweg, Henning, Andernach

The Tail of the HI Mass Function

Authors: S. E. Schneider, J. L. Rosenberg

Abstract: The contribution of extragalactic objects with HI masses below $10^8 M_\odot$ to the HI mass function remains uncertain. Several aspects of the detection of low-mass sources in HI surveys are not always considered, and as a result different analysis techniques yield widely different estimates for their number density. It is suggested at one extreme that the number density of galaxies follows a s… ▽ More The contribution of extragalactic objects with HI masses below $10^8 M_\odot$ to the HI mass function remains uncertain. Several aspects of the detection of low-mass sources in HI surveys are not always considered, and as a result different analysis techniques yield widely different estimates for their number density. It is suggested at one extreme that the number density of galaxies follows a shallow Schechter power-law slope, and at the other extreme that it follows a steep faint-end rise like that found for field optical sources. Here we examine a variety of selection effects, issues of completeness, and consequences of LSS. We derive results for the large Arecibo Dual Beam Survey which indicate that the field mass function does rise steeply, while within the Virgo Cluster environs, the slope appears to be much shallower. Dependence on the local density of galaxies may partially explain differences between surveys. △ Less

Submitted 18 October, 2000; originally announced October 2000.

Comments: 8 pages, presented at Mapping the Hidden Universe: The Universe in HI. eds Kraan-Korteweg, Henning, Andernach

Serendipitous 2MASS Discoveries Near the Galactic Plane: A Spiral Galaxy and Two Globular Clusters

Authors: Robert L. Hurt, Tom H. Jarrett, J. Davy Kirkpatrick, Roc M. Cutri, Stephen E. Schneider, Mike Skrutskie, Willem van Driel

Abstract: We present the basic properties of three objects near the Galactic Plane - a large galaxy and two candidate globular clusters - discovered in the Two Micron All Sky Survey (2MASS) dataset. All were noted during spot-checks of the data during 2MASS quality assurance reviews. The galaxy is a late-type spiral galaxy (Sc-Sd), ~11 Mpc distant. We estimate an extinction-corrected optical diameter of ~… ▽ More We present the basic properties of three objects near the Galactic Plane - a large galaxy and two candidate globular clusters - discovered in the Two Micron All Sky Survey (2MASS) dataset. All were noted during spot-checks of the data during 2MASS quality assurance reviews. The galaxy is a late-type spiral galaxy (Sc-Sd), ~11 Mpc distant. We estimate an extinction-corrected optical diameter of ~9.5', making it larger than most Messier galaxies. The candidate globular clusters are ~2-3' in extent and are hidden optically behind foreground extinctions of Av ~18-21 mag at l ~ 10 deg, b ~ 0 deg. These chance discoveries were not the result of any kind of systematic search but they do hint at the wealth of obscured sources of all kinds, many previously unknown, that are in the 2MASS dataset. △ Less

Submitted 19 June, 2000; originally announced June 2000.

Comments: Accepted 6/15/00 for publication in the Astronomical Journal

IR-TF Relation In The Zone Of Avoidance With 2MASS

Authors: N. Bouche, S. E. Schneider

Abstract: Using the Tully-Fisher relation, one can map the peculiar velocity field and estimate the mass in regions such as the Great Attractor. 2MASS is surveying the full sky in J, H and K bands and has the great advantage that it allows us to detect galaxies more uniformly to lower Galactic latitude than optical surveys. This study shows the feasibility of using 2MASS and the TF relation to probe large… ▽ More Using the Tully-Fisher relation, one can map the peculiar velocity field and estimate the mass in regions such as the Great Attractor. 2MASS is surveying the full sky in J, H and K bands and has the great advantage that it allows us to detect galaxies more uniformly to lower Galactic latitude than optical surveys. This study shows the feasibility of using 2MASS and the TF relation to probe large scale structure. We found that (i) we can use axis ratio $b/a$ up to 0.5; (ii) intrinsic extinction is present (up to 0.5mag at J, 0.1 mag at K); (iii) the zero-point of the TF relation is independent of the 2MASS magnitude measurement and is consistent with the HST Key-Project value; (iv) the 2MASS K-band 20th mag/arcsec$^2$ isophotal aperture magnitude produces the best TF relation; and (v) there is no type dependence of the residuals △ Less

Submitted 21 November, 2000; v1 submitted 10 May, 2000; originally announced May 2000.

Comments: 8 pages, Latex using newpasp.sty. Matches version to appear in Mapping the Hidden Universe: The Universe Behind the Milky Way - The Universe in HI, ASP Conference Series 2000, eds R.C.Kraan-Korteweg, P.A. Henning and H. Andernach (matched to published version; abstract bug in v2 fixed)

The Arecibo Dual-Beam Survey: Arecibo and VLA Observations

Authors: Jessica L. Rosenberg, Stephen E. Schneider

Abstract: The Arecibo Dual-Beam Survey is a "blind" 21 cm search for galaxies covering \~430 deg^2 of sky. We present the data from the detection survey as well as from the follow-up observations to confirm detections and improve positions and flux measurements. We find 265 galaxies, many of which are extremely low surface brightness. Some of these previously uncataloged galaxies lie within the zone of av… ▽ More The Arecibo Dual-Beam Survey is a "blind" 21 cm search for galaxies covering \~430 deg^2 of sky. We present the data from the detection survey as well as from the follow-up observations to confirm detections and improve positions and flux measurements. We find 265 galaxies, many of which are extremely low surface brightness. Some of these previously uncataloged galaxies lie within the zone of avoidance where they are obscured by the gas and dust in our Galaxy. 81 of these sources are not previously cataloged optically and there are 11 galaxies that have no associated optical counterpart or are only tentatively associated with faint wisps of nebulosity on the Digitized Sky Survey images. We discuss the properties of the survey and in particular we make direct determinations of the completeness and reliability of the sample. The behavior of the completeness and its dependencies is essential for determining the HI mass function. We leave the discussion of the mass function for a later paper, but do note that we find many low surface brightness galaxies and 7 sources with M_HI < 10^8 Msolar. △ Less

Submitted 16 April, 2000; v1 submitted 14 April, 2000; originally announced April 2000.

Comments: 23 pages, 20 figures, accepted ApJS. For tables 2 and 3 only the first page has been included. ASCII tables are provided separately

Journal ref: Astrophys.J.Suppl. 130 (2000) 177-200