Showing papers by "Joel R. Primack published in 2013"

Gravitationally Consistent HALO Catalogs and Merger Trees for Precision Cosmology

Abstract: We present a new algorithm for generating merger trees and halo catalogs which explicitly ensures consistency of halo properties (mass, position, and velocity) across time steps. Our algorithm has demonstrated the ability to improve both the completeness (through detecting and inserting otherwise missing halos) and purity (through detecting and removing spurious objects) of both merger trees and halo catalogs. In addition, our method is able to robustly measure the self-consistency of halo finders; it is the first to directly measure the uncertainties in halo positions, halo velocities, and the halo mass function for a given halo finder based on consistency between snapshots in cosmological simulations. We use this algorithm to generate merger trees for two large simulations (Bolshoi and Consuelo) and evaluate two halo finders (ROCKSTAR and BDM). We find that both the ROCKSTAR and BDM halo finders track halos extremely well; in both, the number of halos which do not have physically consistent progenitors is at the 1%-2% level across all halo masses. Our code is publicly available at http://code.google.com/p/consistent-trees. Our trees and catalogs are publicly available at http://hipacc.ucsc.edu/Bolshoi/.

CANDELS: The progenitors of compact quiescent galaxies at z ~ 2

Abstract: We combine high-resolution Hubble Space Telescope/WFC3 images with multi-wavelength photometry to track the evolution of structure and activity of massive (M_*> 10^10 M_☉) galaxies at redshifts z = 1.4-3 in two fields of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey. We detect compact, star-forming galaxies (cSFGs) whose number densities, masses, sizes, and star formation rates (SFRs) qualify them as likely progenitors of compact, quiescent, massive galaxies (cQGs) at z = 1.5-3. At z≲2, cSFGs present SFR = 100-200 M_☉ yr^–1, yet their specific star formation rates (sSFR ~ 10^–9 yr^–1) are typically half that of other massive SFGs at the same epoch, and host X-ray luminous active galactic nuclei (AGNs) 30 times (~30%) more frequently. These properties suggest that cSFGs are formed by gas-rich processes (mergers or disk-instabilities) that induce a compact starburst and feed an AGN, which, in turn, quench the star formation on dynamical timescales (few 10^8 yr). The cSFGs are continuously being formed at z = 2-3 and fade to cQGs down to z ~ 1.5. After this epoch, cSFGs are rare, thereby truncating the formation of new cQGs. Meanwhile, down to z = 1, existing cQGs continue to enlarge to match local QGs in size, while less-gas-rich mergers and other secular mechanisms shepherd (larger) SFGs as later arrivals to the red sequence. In summary, we propose two evolutionary tracks of QG formation: an early (z≲2), formation path of rapidly quenched cSFGs fading into cQGs that later enlarge within the quiescent phase, and a late-arrival (z≳2) path in which larger SFGs form extended QGs without passing through a compact state.

Toy models for galaxy formation versus simulations

Abstract: We describe simple useful toy models for key processes of galaxy formation in its most active phase, at z > 1, and test the approximate expressions against the typical behaviour in a suite of high-resolution hydro-cosmological simulations of massive galaxies at z = 4-1. We address in particular the evolution of (a) the total mass inflow rate from the cosmic web into galactic haloes based on the EPS approximation, (b) the penetration of baryonic streams into the inner galaxy, (c) the disc size, (d) the implied steady-state gas content and star-formation rate (SFR) in the galaxy subject to mass conservation and a universal star-formation law, (e) the inflow rate within the disc to a central bulge and black hole as derived using energy conservation and self-regulated Q ~ 1 violent disc instability (VDI), and (f) the implied steady state in the disc and bulge. The toy models provide useful approximations for the behaviour of the simulated galaxies. We find that (a) the inflow rate is proportional to mass and to (1+z)^5/2, (b) the penetration to the inner halo is ~50% at z = 4-2, (c) the disc radius is ~5% of the virial radius, (d) the galaxies reach a steady state with the SFR following the accretion rate into the galaxy, (e) there is an intense gas inflow through the disc, comparable to the SFR, following the predictions of VDI, and (f) the galaxies approach a steady state with the bulge mass comparable to the disc mass, where the draining of gas by SFR, outflows and disc inflows is replenished by fresh accretion. Given the agreement with simulations, these toy models are useful for understanding the complex phenomena in simple terms and for back-of-the-envelope predictions.

The AGORA high-resolution galaxy simulations comparison project

Abstract: We introduce the Assembling Galaxies Of Resolved Anatomy (AGORA) project, a comprehensive numerical study of well-resolved galaxies within the ΛCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of eight galaxies with halo masses M_(vir) ≃ 10^(10), 10^(11), 10^(12), and 10^(13) M_☉ at z = 0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit and validated against observations to verify that the solutions are robust—i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The initial conditions for the AGORA galaxies as well as simulation outputs at various epochs will be made publicly available to the community. The proof-of-concept dark-matter-only test of the formation of a galactic halo with a z = 0 mass of M_(vir) ≃ 1.7 × 10^(11) M_☉ by nine different versions of the participating codes is also presented to validate the infrastructure of the project.

The AGORA High-Resolution Galaxy Simulations Comparison Project

Abstract: We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project.

The MultiDark Database: release of the Bolshoi and MultiDark cosmological simulations

Abstract: We present the online MultiDark Database – a Virtual Observatory-oriented, relational database for hosting various cos-mological simulations. The data is accessible via an SQL (Structured Query Language) query interface, which also allows users to directly pose scientific questions, as shown in a number of examples in this paper. Further examples for the usage of the database are given in its extensive online documentation. The database is based on the same technology as the Millennium Database, a fact that will greatly facilitate the usage of both suites of cosmological simulations. The first release of the MultiDark Database hosts two 8.6 billion particle cosmological N-body simulations: the Bolshoi (250 h–1 Mpc simulation box, 1 h–1 kpc resolution) and MultiDark Run1 simulation (MDR1, or BigBolshoi, 1000 h–1 Mpc simulation box, 7 h–1 kpc resolution). The extraction methods for halos/subhalos from the raw simulation data, and how this data is structured in the database are explained in this paper. With the first data release, users get full access to halo/subhalo catalogs, various profiles of the halos at redshifts z = 0–15, and raw dark matter data for one time-step of the Bolshoi and four time-steps of the MultiDark simulation. Later releases will also include galaxy mock catalogs and additional merger trees for both simulations as well as new large volume simulations with high resolution. This project is further proof of the viability to store and present complex data using relational database technology. We encourage other simulators to publish their results in a similar manner. (© 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

Detection of the Cosmic γ-Ray Horizon from Multiwavelength Observations of Blazars

Abstract: The first statistically significant detection of the cosmic γ-ray horizon (CGRH) that is independent of any extragalactic background light (EBL) model is presented. The CGRH is a fundamental quantity in cosmology. It gives an estimate of the opacity of the universe to very high energy (VHE) γ-ray photons due to photon-photon pair production with the EBL. The only estimations of the CGRH to date are predictions from EBL models and lower limits from γ-ray observations of cosmological blazars and γ-ray bursts. Here, we present homogeneous synchrotron/synchrotron self-Compton (SSC) models of the spectral energy distributions of 15 blazars based on (almost) simultaneous observations from radio up to the highest energy γ-rays taken with the Fermi satellite. These synchrotron/SSC models predict the unattenuated VHE fluxes, which are compared with the observations by imaging atmospheric Cherenkov telescopes. This comparison provides an estimate of the optical depth of the EBL, which allows us a derivation of the CGRH through a maximum likelihood analysis that is EBL-model independent. We find that the observed CGRH is compatible with the current knowledge of the EBL.

A Comparison between Semi-Analytic Model Predictions for the CANDELS Survey

Abstract: We compare the predictions of three independently developed semi-analytic galaxy formation models that are being used to aid in the interpretation of results from the CANDELS survey. These models are each applied to the same set of halo merger trees extracted from the "Bolshoi" simulation and are carefully tuned to match the local galaxy stellar mass function using the powerful method of Bayesian Inference coupled with MCMC or by hand. The comparisons reveal that in spite of the significantly different parameterizations for star formation and feedback processes, the three models yield qualitatively similar predictions for the assembly histories of galaxy stellar mass and star formation over cosmic time. We show that the SAMs generally require strong outflows to suppress star formation in low-mass halos to match the present day stellar mass function. However, all of the models considered produce predictions for the star formation rates and metallicities of low-mass galaxies that are inconsistent with existing data and diverge between the models. We suggest that large differences in the metallicity relations and small differences in the stellar mass assembly histories of model galaxies stem from different assumptions for the outflow mass-loading factor. Importantly, while more accurate observational measurements for stellar mass, SFR and metallicity of galaxies at 1

Detection of the cosmic \gamma-ray horizon from multiwavelength observations of blazars

Abstract: The first statistically significant detection of the cosmic \gamma-ray horizon (CGRH) that is independent of any extragalactic background light (EBL) model is presented. The CGRH is a fundamental quantity in cosmology. It gives an estimate of the opacity of the Universe to very high energy (VHE) \gamma-ray photons due to photon-photon pair production with the EBL. The only estimations of the CGRH to date are predictions from EBL models and lower limits from \gamma-ray observations of cosmological blazars and \gamma-ray bursts. Here, we present homogeneous synchrotron/synchrotron self-Compton (SSC) models of the spectral energy distributions of 15 blazars based on (almost) simultaneous observations from radio up to the highest energy \gamma-rays taken with the Fermi satellite. These synchrotron/SSC models predict the unattenuated VHE fluxes, which are compared with the observations by imaging atmospheric Cherenkov telescopes. This comparison provides an estimate of the optical depth of the EBL, which allows a derivation of the CGRH through a maximum likelihood analysis that is EBL-model independent. We find that the observed CGRH is compatible with the current knowledge of the EBL.

Confronting simulations of optically thick gas in massive halos with observations at z = 2-3

Abstract: Cosmological hydrodynamic simulations predict the physical state of baryons in the circumgalactic medium (CGM), which can be directly tested via quasar absorption line observations. We use high-resolution "zoom-in" simulations of 21 galaxies to characterize the distribution of neutral hydrogen around halos in the mass range M vir ~ 2 × 1011 to 4 × 1012 M ☉ at z ~ 2. We find that both the mass fraction of cool (T ≤ 3 × 104 K) gas and the covering fraction of optically thick Lyman limit systems (LLSs) depend only weakly on halo mass, even around the critical value for the formation of stable virial shocks. The covering fraction of LLSs interior to the virial radius varies between f c ~ 0.05-0.2, with significant scatter among halos. Our simulations of massive halos (M vir ≥ 1012 M ☉) underpredict the covering fraction of optically thick gas observed in the quasar CGM by a large factor. The reason for this discrepancy is unclear, but several possibilities are discussed. In the lower mass halos (M vir ≥ 5 × 1011 M ☉) hosting star-forming galaxies, the predicted covering factor agrees with observations; however, current samples of quasar-galaxy pairs are too small for a conclusive comparison. To overcome this limitation, we propose a new observable: the small-scale autocorrelation function of optically thick absorbers detected in the foreground of close quasar pairs. We show that this new observable can constrain the underlying dark halos hosting LLSs at z ~ 2-3, as well as the characteristic size and covering factor of the CGM.

Spatially unassociated galaxies contribute significantly to the blended submillimetre galaxy population: predictions for follow-up observations of ALMA sources

Abstract: There is anecdotal evidence that spatially and physically unassociated galaxies blended into a single submillimetre (submm) source contribute to the submm galaxy (SMG) population. This work is the first to theoretically predict the number counts of such sources. We generate mock SMG catalogues using lightcones derived from the Bolshoi cosmological simulation; to assign submm flux densities to the mock galaxies, we use a fitting function previously derived from the results of dust radiative transfer performed on hydrodynamical simulations of isolated disc and merging galaxies. We then calculate submm number counts for different beam sizes and without blending. We predict that > ~50 per cent of blended SMGs have at least one spatially unassociated component with S_850 > 1 mJy. For a 15-arcsec beam, blends of >2 galaxies in which at least one component is spatially unassociated dominate the blended sources with total S_850 > ~3 mJy. The distribution of the redshift separations amongst the components is strongly bimodal. The typical redshift separation of spatially unassociated blended sources is ~1. Our predictions for the contributions of spatially unassociated components and the distribution of redshift separations are not testable with currently available data, but they will be easily tested once sufficiently accurate redshifts for the individual subcomponents (resolved by, e.g., ALMA) of a sufficient number of single-dish-detected blended SMGs are available.

THE FIRM REDSHIFT LOWER LIMIT OF THE MOST DISTANT TeV-DETECTED BLAZAR PKS 1424+240

Abstract: We present the redshift lower limit of z ≥ 0.6035 for the very high energy (VHE; E ≥ 100 GeV) emitting blazar PKS 1424+240 (PG 1424+240). This limit is inferred from Lyβ and Lyγ absorption observed in the far-ultraviolet spectra from the Hubble Space Telescope/Cosmic Origins Spectrograph. No VHE-detected blazar has shown solid spectroscopic evidence of being more distant. At this distance, VHE observations by VERITAS are shown to sample historically large gamma-ray opacity values at 500 GeV, extending beyond τ = 4 for low-level models of the extragalactic background light (EBL) and beyond τ = 5 for high levels. The majority of the z = 0.6035 absorption-corrected VHE spectrum appears to exhibit a lower flux than an extrapolation of the contemporaneous Large Area Telescope power-law fit beyond 100 GeV. However, the highest energy VERITAS point is the only point showing agreement with this extrapolation, possibly implying the overestimation of the gamma-ray opacity or the onset of an unexpected VHE spectral feature. A curved log parabola is favored when fitting the full range of gamma-ray data (0.5-500 GeV). While fitting the absorption-corrected VHE data alone results in a harder differential power law than that from the full range, the indices derived using three EBL models are consistent with the physically motivated limit set by Fermi acceleration processes.

The Firm Redshift Lower Limit of the Most Distant TeV-Detected Blazar PKS 1424+240

Abstract: We present the redshift lower limit of z>0.6035 for the very-high-energy (VHE; E>100 GeV) emitting blazar PKS 1424+240 (PG 1424+240). This limit is inferred from Lyman beta and gamma absorption observed in the far-ultraviolet spectra from the Hubble Space Telescope/Cosmic Origins Spectrograph. No VHE-detected blazar has shown solid spectroscopic evidence of being more distant. At this distance, VHE observations by VERITAS are shown to sample historically large gamma-ray opacity values at 500 GeV, extending beyond tau=4 for low-level models of the extragalactic background light (EBL) and beyond tau=5 for high-levels. The majority of the z=0.6035 absorption-corrected VHE spectrum appears to exhibit a lower flux than an extrapolation of the contemporaneous LAT power-law fit beyond 100 GeV. However, the highest energy VERITAS point is the only point showing agreement with this extrapolation, possibly implying the overestimation of the gamma-ray opacity or the onset of an unexpected VHE spectral feature. A curved log parabola is favored when fitting the full range of gamma-ray data (0.5 to 500 GeV). While fitting the absorption-corrected VHE data alone results in a harder differential power law than that from the full range, the indices derived using three EBL models are consistent with the physically motivated limit set by Fermi acceleration processes.

IACT observations of gamma-ray bursts: prospects for the Cherenkov Telescope Array

Abstract: Gamma rays at rest frame energies as high as 90 GeV have been reported from gamma-ray bursts (GRBs) by the Fermi Large Area Telescope (LAT). There is considerable hope that a confirmed GRB detection will be possible with the upcoming Cherenkov Telescope Array (CTA), which will have a larger effective area and better low-energy sensitivity than current-generation imaging atmospheric Cherenkov telescopes (IACTs). To estimate the likelihood of such a detection, we have developed a phenomenological model for GRB emission between 1 GeV and 1 TeV that is motivated by the high-energy GRB detections of Fermi-LAT, and allows us to extrapolate the statistics of GRBs seen by lower energy instruments such as the Swift-BAT and BATSE on the Compton Gamma-ray Observatory. We show a number of statistics for detected GRBs, and describe how the detectability of GRBs with CTA could vary based on a number of parameters, such as the typical observation delay between the burst onset and the start of ground observations. We also consider the possibility of using GBM on Fermi as a finder of GRBs for rapid ground follow-up. While the uncertainty of GBM localization is problematic, the small field-of-view for IACTs can potentially be overcome by scanning over the GBM error region. Overall, our results indicate that CTA should be able to detect one GRB every 20–30 months with our baseline instrument model, assuming consistently rapid pursuit of GRB alerts, and provided that spectral breaks below ~100 GeV are not a common feature of the bright GRB population. With a more optimistic instrument model, the detection rate can be as high as 1 to 2 GRBs per year.

Confronting Simulations of Optically Thick Gas in Massive Halos with Observations at z=2-3

Abstract: Cosmological hydrodynamic simulations predict the physical state of baryons in the circumgalactic medium (CGM), which can be directly tested via quasar absorption line observations. We use high resolution "zoom-in" simulations of 21 galaxies to characterize the distribution of neutral hydrogen around halos in the mass range M_vir~ 2x10^11 - 4x10^12 M_sun at z~2. We find that both the mass fraction of cool (T = 10^12 M_sun) underpredict the covering fraction of optically-thick gas observed in the quasar CGM by a large factor. The reason for this discrepancy is unclear, but several possibilities are discussed. In the lower mass halos (M_vir >= 5x10^11 M_sun) hosting star-forming galaxies, the predicted covering factor agrees with observations, however current samples of quasar-galaxy pairs are too small for a conclusive comparison. To overcome this limitation, we propose a new observable: the small-scale auto-correlation function of optically-thick absorbers detected in the foreground of close quasar pairs. We show that this new observable can constrain the underlying dark halos hosting LLSs at z~2-3, as well as the characteristic size and covering factor of the CGM.

Mergers and Mass Accretion for Infalling Halos Both End Well Outside Cluster Virial Radii

Abstract: We find that infalling dark matter halos (i.e., the progenitors of satellite halos) begin losing mass well outside the virial radius of their eventual host halos. The peak mass occurs at a range of clustercentric distances, with median and 68th percentile range of 1.8 +2.3/-1.0 R_(vir,host) for progenitors of z=0 satellites. The peak circular velocity for infalling halos occurs at significantly larger distances (3.7 +3.3/-2.2 R_(vir,host) at z=0). This difference arises because different physical processes set peak circular velocity (typically, ~1:5 and larger mergers which cause transient circular velocity spikes) and peak mass (typically, smooth accretion) for infalling halos. We find that infalling halos also stop having significant mergers well before they enter the virial radius of their eventual hosts. Mergers larger than a 1:40 ratio in halo mass end for infalling halos at similar clustercentric distances (~ 1.9 R_(vir,host)) as the end of overall mass accretion. However, mergers larger than 1:3 typically end for infalling halos at more than 4 virial radial away from their eventual hosts. This limits the ability of mergers to affect quenching and morphology changes in clusters. We also note that the transient spikes which set peak circular velocity may lead to issues with abundance matching on that parameter, including unphysical galaxy stellar mass growth profiles near clusters; we propose a simple observational test to check if a better halo proxy for galaxy stellar mass exists.