The James Webb Space Telescope (JWST) is expected to observe galaxies at $z>10$ that are presently inaccessible. Here, we use a self-consistent empirical model, the UniverseMachine, to generate mock galaxy catalogues and lightcones over the redshift range $z=0-15$. These data include realistic galaxy properties (stellar masses, star formation rates, and UV luminosities), galaxy-halo relationships, and galaxy-galaxy clustering. Mock observables are also provided for different model parameters spanning observational uncertainties at $z<10$. We predict that Cycle 1 JWST surveys will very likely detect galaxies with $M_*>10^7 M_\odot$ and/or $M_{1500}<-17$ out to at least $z\sim 13.5$. Number density uncertainties at $z>12$ expand dramatically, so efforts to detect $z>12$ galaxies will provide the most valuable constraints on galaxy formation models. The faint-end slopes of the stellar mass/luminosity functions at a given mass/luminosity threshold steepen as redshift increases. This is because observable galaxies are hosted by haloes in the exponentially falling regime of the halo mass function at high redshifts. Hence, these faint-end slopes are robustly predicted to become shallower below current observable limits ($M_\ast < 10^7M_\odot$ or $M_\mathrm{1500}>-17$). For reionization models, extrapolating luminosity functions with a constant faint-end slope from $M_{1500}=-17$ down to $M_{1500}=-12$ gives the most reasonable upper limit for the total UV luminosity and cosmic star formation rate up to $z\sim 12$. We compare to three other empirical models and one semi-analytic model, showing that the range of predicted observables from our approach encompasses predictions from other techniques. Public catalogues and lightcones for common fields are available online.

LCDM cosmological models with Early Dark Energy (EDE) have been proposed to resolve tensions between the Hubble constant H0 = 100h km/s/Mpc measured locally, giving h ~ 0.73, and H0 deduced from Planck cosmic microwave background (CMB) and other early universe measurements plus LCDM, giving h ~ 0.67. EDE models do this by adding a scalar field that temporarily adds dark energy equal to about 10% of the cosmological energy density at the end of the radiation-dominated era at redshift z ~ 3500. Here we compare linear and nonlinear predictions of a Planck-normalized LCDM model including EDE giving h = 0.728 with those of standard Planck-normalized LCDM with h = 0.678. We find that nonlinear evolution reduces the differences between power spectra of fluctuations at low redshifts. As a result, at z = 0 the halo mass functions on galactic scales are nearly the same, with differences only 1-2%. However, the differences dramatically increase at high redshifts. The EDE model predicts 50% more massive clusters at z = 1 and twice more galaxy-mass halos at z = 4. Even greater increases in abundances of galaxy-mass halos at higher redshifts may make it easier to reionize the universe with EDE. Predicted galaxy abundances and clustering will soon be tested by JWST observations. Positions of baryonic acoustic oscillations (BAOs) and correlation functions differ by about 2% between the models -- an effect that is not washed out by nonlinearities. Both standard LCDM and the EDE model studied here agree well with presently available acoustic-scale observations, but DESI and Euclid measurements will provide stringent new tests.

Galaxies forming in low-mass halos are thought to be primarily responsible for reionizing the Universe during the first billion years after the Big Bang. Yet, these halos are extremely inefficient at forming stars in the nearby Universe. In this work, we address this apparent tension, and ask whether a physically motivated model of galaxy formation that reproduces the observed abundance of faint galaxies in the nearby Universe is also consistent with available observational constraints on the reionization history. By interfacing the Santa Cruz semi-analytic model for galaxy formation with an analytic reionization model, we constructed a computationally efficient pipeline that connects `ground-level' galaxy formation physics to `top-level' cosmological-scale observables. Based on photometric properties of the galaxy populations predicted up to $z=15$, we compute the reionization history of intergalactic hydrogen. We quantify the three degenerate quantities that influence the total ionizing photon budget, including the abundance of galaxies, the intrinsic production rate of ionizing photons, and the LyC escape fraction. We explore covariances between these quantities using a Markov chain Monte Carlo method. We find that our locally calibrated model is consistent with all currently available constraints on the reionization history, under reasonable assumptions about the LyC escape fraction. We quantify the fraction of ionizing photons contributed by galaxies of different luminosities and find that the galaxies expected to be detected in JWST NIRCam wide and deep surveys are responsible for producing $\sim 40$-$80\%$ of ionizing photons throughout the EoR. All results presented in this work are available at https://www.simonsfoundation.org/semi-analytic-forecasts-for-jwst/.

The correlation between fresh gas accretion onto haloes and galaxy star formation is critical to understanding galaxy formation. Different theoretical models have predicted different correlation strengths between halo accretion rates and galaxy star formation rates, ranging from strong positive correlations to little or no correlation. Here, we present a technique to observationally measure this correlation strength for isolated Milky Way-mass galaxies with $z < 0.123$. This technique is based on correlations between dark matter accretion rates and the projected density profile of neighbouring galaxies; these correlations also underlie past work with splashback radii. We apply our technique to both observed galaxies in the Sloan Digital Sky Survey as well as simulated galaxies in the UniverseMachine where we can test any desired correlation strength. We find that positive correlations between dark matter accretion and recent star formation activity are ruled out with $\gtrsim 85\%$ confidence. Our results suggest that star formation activity may not be correlated with fresh accretion for isolated Milky Way-mass galaxies at $z=0$ and that other processes, such as gas recycling, dominate further galaxy growth.

General-education college astronomy courses offer instructors both a unique audience and a unique challenge. For many students, such a course may be their first time encountering a standalone astronomy class, and it is also likely one of the last science courses they will take. Thus, in a single semester, primary course goals often include both imparting knowledge about the Universe and giving students some familiarity with the processes of science. In traditional course environments, students often compartmentalize information into separate "life files" and "course files" rather than integrating information into a coherent framework. The astronomy course created through this project, taught at the University of Arizona in Spring 2019, was designed around inclusivity-driven guiding principles that help students engage with course content in ways that are meaningful, relevant, and accessible. Our course bridges the gap between students' "life" and "course files", encourages and respects diverse points of view, and empowers students to connect course content with their personal lives and identities. In this paper, we provide insight into the guiding principles that informed our course design and share research results on the effectiveness of the instructional strategies and assessment techniques implemented in the course.

We continue our empirical study of the emission line flux originating in the cool ($T\sim10^4$ K) gas that populates the halos of galaxies and their environments. Specifically, we present results obtained for a sample of galaxy pairs with a range of projected separations, {\bf $10 < {S_p/\rm kpc} < 200$}, and mass ratios $<$ 1:5, intersected by 5,443 SDSS lines of sight at projected radii of 10 to 50 kpc from either or both of the two galaxies. We find significant enhancement in H$\alpha$ emission and a moderate enhancement in [N {\small II}]6583 emission for low mass pairs (mean stellar mass per galaxy, $\overline{\rm M}_*, <10^{10.4} {\rm M}_\odot$) relative to the results from a control sample. This enhanced H$\alpha$ emission comes almost entirely from sight lines located between the galaxies, consistent with a short-term, interaction-driven origin for the enhancement. We find no enhancement in H$\alpha$ emission, but significant enhancement in [N {\small II}]6583 emission for high mass ($\overline{\rm M}_* >10^{10.4}{\rm M}_\odot$) pairs. Furthermore, we find a dependence of the emission line properties on the galaxy pair mass ratio such that those with a mass ratio below 1:2.5 have enhanced [N {\small II}]6583 and those with a mass ratio between 1:2.5 and 1:5 do not. In all cases, departures from the control sample are only detected for close pairs ($S_p <$ 100 kpc). Attributing an elevated [N {\small II}]6583/H$\alpha$ ratio to shocks, we infer that shocks play a role in determining the CGM properties for close pairs that are among the more massive and have mass ratios closer to 1:1.

We present a new semi-empirical model for the dust continuum number counts of galaxies at 1.1 millimeter and 850 \micron. Our approach couples an observationally motivated model for the stellar mass and SFR distribution of galaxies with empirical scaling relations to predict the dust continuum flux density of these galaxies. Without a need to tweak the IMF, the model reproduces the currently available observations of the 1.1 millimeter and 850 \micron number counts, including the observed flattening in the 1.1 millimeter number counts below 0.3 mJy \citep{Gonzalez2019numbercounts} and the number counts in discrete bins of different galaxy properties. Predictions of our work include : (1) the galaxies that dominate the number counts at flux densities below 1 mJy (3 mJy) at 1.1 millimeter (850 $\mu$m) have redshifts between $z=1$ and $z=2$, stellar masses of $\sim 5\times10^{10}~\rm{M}_\odot$, and dust masses of $\sim 10^{8}~\rm{M}_\odot$; (2) the flattening in the observed 1.1 millimeter number counts corresponds to the knee of the 1.1 millimeter luminosity function. A similar flattening is predicted for the number counts at 850 $\mu$m; (3) the model reproduces the redshift distribution of current 1.1 millimeter detections; (4) to efficiently detect large numbers of galaxies through their dust continuum, future surveys should scan large areas once reaching a 1.1 millimeter flux density of 0.1 mJy rather than integrating to fainter fluxes. Our modeling framework also suggests that the amount of information on galaxy physics that can be extracted from the 1.1 millimeter and 850 $\mu$m number counts is almost exhausted.

We present a search for companion [CII] emitters to known luminous sources at $6<$ z $<6.5$ in deep, archival ALMA observations. The observations are deep enough to detect sources with L$_{\rm [CII]} \sim 10^8$ at z $\sim6$. We identify four robust line detections from a blind search of five deep fields centered on ultra-luminous infrared galaxies and QSOs, over an order of magnitude more than expected based on current observations and predictions, suggesting that these objects may be highly biased tracers of mass in the early Universe. We find these companion lines to have comparable properties to other known galaxies at the same epoch. All companions lie less than 650 km s$^{-1}$ and between 20 -- 70 kpc (projected) from their central source, providing a constraint on their halo masses of the central galaxies ranging from 2.5$\times$10$^{12}$ M$_\odot$ to 4$\times$10$^{13}$ M$_\odot$. To place these discoveries in context, we employ a mock galaxy catalog to estimate the luminosity function for [CII] during reionization and compare to our observations. The simulations support this result by showing a similar level of elevated counts found around such luminous sources.

Observational systematics complicate comparisons with theoretical models limiting understanding of galaxy evolution. In particular, different empirical determinations of the stellar mass function imply distinct mappings between the galaxy and halo masses, leading to diverse galaxy evolutionary tracks. Using our state-of-the-art STatistical sEmi-Empirical modeL, STEEL, we show fully self-consistent models capable of generating galaxy growth histories that simultaneously and closely agree with the latest data on satellite richness and star-formation rates at multiple redshifts and environments. Central galaxy histories are generated using the central halo mass tracks from state-of-the-art statistical dark matter accretion histories coupled to abundance matching routines. We show that too flat high-mass slopes in the input stellar-mass-halo-mass relations as predicted by previous works, imply non-physical stellar mass growth histories weaker than those implied by satellite accretion alone. Our best-fit models reproduce the satellite distributions at the largest masses and highest redshifts probed, the latest data on star formation rates and its bi-modality in the local Universe, and the correct fraction of ellipticals. Our results are important to predict robust and self-consistent stellar-mass-halo-mass relations and to generate reliable galaxy mock catalogues for the next generations of extra-galactic surveys such as Euclid and LSST.

We continue our empirical study of the emission line flux originating in the cool ($T\sim10^4$ K) gas that populates the halos of galaxies and their environments. Specifically, we present results obtained for a sample of nearly half a million individual galaxies, groups, and clusters of galaxies, intersected by more than two million SDSS lines of sight at projected separations of up to a quarter of the virial radius. Adopting simple power law relationships between the circumgalactic (CGM) cool gas fraction and either the halo or stellar mass, we present expressions for the CGM cool gas fraction as a function of either halo or stellar mass, $f_{\rm cool}(M_h) = (0.25^{+0.07}_{-0.06}) \times (M_h/10^{12}M_\odot)^{(-0.39^{+0.06}_{-0.07})}$ or $f_{\rm cool}(M_{*}) = (0.28^{+0.07}_{-0.06}) \times (M_{\rm *}/10^{10.0}M_\odot)^{(-0.33\pm 0.06)}$. Where we can compare, our results are consistent with previous constraints from absorption line studies, our own previous emission line work, and simulations. The cool gas can be the dominant baryonic CGM component, comprising a fraction as high as $> 90\%$ of halo gaseous baryons, in low mass halos, $M_h\sim$ $10^{10.5} M_\odot$, and a minor fraction, $<$ 5\%, in groups and clusters, $M_h > 10^{14} M_\odot$.

We study the effects of the local environmental density and the cosmic web environment (filaments, walls, and voids) on key properties of dark matter halos using the Bolshoi-Planck LCDM cosmological simulation. The z = 0 simulation is analysed into filaments, walls, and voids using the SpineWeb method and also the VIDE package of tools, both of which use the watershed transform. The key halo properties that we study are the specific mass accretion rate, spin parameter, concentration, prolateness, scale factor of the last major merger, and scale factor when the halo had half of its z = 0 mass. For all these properties, we find that there is no discernible difference between the halo properties in filaments, walls, or voids when compared at the same environmental density. As a result, we conclude that environmental density is the core attribute that affects these properties. This conclusion is in line with recent findings that properties of galaxies in redshift surveys are independent of their cosmic web environment at the same environmental density at z ~ 0. We also find that the local web environment of the Milky Way and the Andromeda galaxies near the centre of a cosmic wall does not appear to have any effect on the properties of these galaxies' dark matter halos except for their orientation, although we find that it is rather rare to have such massive halos near the centre of a relatively small cosmic wall.

We present the measurements and modelling of the small-to-intermediate scale (0.1--25 Mpc/h) projected and three-dimensional (3D) redshift-space two-point correlation functions (2PCFs) of local galaxies in the Sloan Digital Sky Survey (SDSS) Data Release 7. We find a clear dependence of galaxy clustering on luminosity in both projected and redshift spaces, generally being stronger for more luminous samples. The measurements are successfully interpreted within the halo occupation distribution (HOD) framework with central and satellite velocity bias parameters to describe galaxy kinematics inside haloes and to model redshift-space distortion (RSD) effects. In agreement with previous studies, we find that more luminous galaxies reside in more massive haloes. Including the redshift-space 2PCFs helps tighten the HOD constraints. Moreover, we find that luminous central galaxies are not at rest at the halo centres, with the velocity dispersion about 30% that of the dark matter. Such a relative motion may reflect the consequence of galaxy and halo mergers, and we find that central galaxies in lower mass haloes tend to be more relaxed with respect to their host haloes. The motion of satellite galaxies in luminous samples is consistent with their following that of the dark matter. For faint samples, satellites tends to have slower motion, with velocity dispersion inside haloes about 85% that of the dark matter. We discuss possible applications of the velocity bias constraints on studying galaxy evolution and cosmology. In the appendix, we characterize the distribution of galaxy redshift measurement errors, which is well described by a Gaussian-convolved double exponential distribution.

ATLAS (Astrophysics Telescope for Large Area Spectroscopy) is a concept for a NASA probe-class space mission. It is the spectroscopic follow-up mission to WFIRST, boosting its scientific return by obtaining deep NIR & MIR slit spectroscopy for most of the galaxies imaged by the WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for ~200M galaxies out to z=7 and beyond, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS and WFIRST together will produce a definitive 3D map of the Universe over 2000 sq deg. ATLAS Science Goals are: (1) Discover how galaxies have evolved in the cosmic web of dark matter from cosmic dawn through the peak era of galaxy assembly. (2) Discover the nature of cosmic acceleration. (3) Probe the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy. (4) Discover the bulk compositional building blocks of planetesimals formed in the outer Solar System. These flow down to the ATLAS Scientific Objectives: (1A) Trace the relation between galaxies and dark matter with less than 10% shot noise on relevant scales at 1<z<7. (1B) Probe the physics of galaxy evolution at 1<z<7. (2) Obtain definitive measurements of dark energy and tests of General Relativity. (3) Measure the 3D structure and stellar content of the inner Milky Way to a distance of 25 kpc. (4) Detect and quantify the composition of 3,000 planetesimals in the outer Solar System. ATLAS is a 1.5m telescope with a FoV of 0.4 sq deg, and uses Digital Micro-mirror Devices (DMDs) as slit selectors. It has a spectroscopic resolution of R = 1000, and a wavelength range of 1-4 microns. ATLAS has an unprecedented spectroscopic capability based on DMDs, with a spectroscopic multiplex factor ~6,000. ATLAS is designed to fit within the NASA probe-class space mission cost envelope.

We present a method to flexibly and self-consistently determine individual galaxies' star formation rates (SFRs) from their host haloes' potential well depths, assembly histories, and redshifts. The method is constrained by galaxies' observed stellar mass functions, SFRs (specific and cosmic), quenched fractions, UV luminosity functions, UV-SM relations, IRX-UV relations, auto- and cross-correlation functions (including quenched and star-forming subsamples), and quenching dependence on environment; each observable is reproduced over the full redshift range available, up to 0<z<10. Key findings include: galaxy assembly correlates strongly with halo assembly; quenching at z>1 correlates strongly with halo mass; quenched fractions at fixed halo mass decrease with increasing redshift; massive quenched galaxies reside in higher-mass haloes than star-forming galaxies at fixed galaxy mass; star-forming and quenched galaxies' star formation histories at fixed mass differ most at z<0.5; satellites have large scatter in quenching timescales after infall, and have modestly higher quenched fractions than central galaxies; Planck cosmologies result in up to 0.3 dex lower stellar mass-halo mass ratios at early times; and, nonetheless, stellar mass-halo mass ratios rise at z>5. Also presented are revised stellar mass-halo mass relations for all, quenched, star-forming, central, and satellite galaxies; the dependence of star formation histories on halo mass, stellar mass, and galaxy SSFR; quenched fractions and quenching timescale distributions for satellites; and predictions for higher-redshift galaxy correlation functions and weak lensing surface densities. The public data release (DR1) includes the massively parallel (>10^5 cores) implementation (the UniverseMachine), the newly compiled and remeasured observational data, derived galaxy formation constraints, and mock catalogues including lightcones.

Most galaxies are hosted by massive, invisible dark matter haloes, yet little is known about the scatter in the stellar mass--halo mass relation for galaxies with host halo masses $M_{h}\le 10^{11}M_{\odot}$. Using mock catalogues based on dark matter simulations, we find that two observable signatures are sensitive to scatter in the stellar mass--halo mass relation even at these mass scales; i.e., conditional stellar mass functions and velocity distribution functions for neighbouring galaxies. We compute these observables for 179,373 galaxies in the Sloan Digital Sky Survey (SDSS) with stellar masses $M_{\ast} > 10^9 M_{\odot}$ and redshifts 0.01 $< z <$ 0.307. We then compare to mock observations generated from the $\textit{Bolshoi-Planck}$ dark matter simulation for stellar mass--halo mass scatters ranging from 0 to 0.6 dex. The observed results are consistent with simulated results for most values of scatter ($<$0.6 dex), and SDSS statistics are insufficient to provide firm constraints. However, this method could provide much tighter constraints on stellar mass--halo mass scatter in the future if applied to larger data sets, especially the anticipated Dark Energy Spectroscopic Instrument Bright Galaxy Survey. Constraining the value of scatter could have important implications for galaxy formation and evolution.

We present Searches After Gravitational-waves Using ARizona Observatories (SAGUARO), a comprehensive effort dedicated to the discovery and characterization of optical counterparts to gravitational wave (GW) events. SAGUARO utilizes ground-based facilities ranging from 1.5m to 10m in diameter, located primarily in the Northern Hemisphere. We provide an overview of SAGUARO's telescopic resources, pipeline for transient detection, and database for candidate visualization. We describe SAGUARO's discovery component, which utilizes the $5$~deg$^2$ field-of-view optical imager on the Mt. Lemmon 1.5m telescope, reaching limits of $\approx 21.3$~AB mag while rapidly tiling large areas. We also describe the follow-up component of SAGUARO, used for rapid vetting and monitoring of optical candidates. With the onset of Advanced LIGO/Virgo's third observing run, we present results from the first three SAGUARO searches following the GW events S190408an, S190425z and S190426c, which serve as a valuable proof-of-concept of SAGUARO. We triggered and searched 15, 60 and 60 deg$^{2}$ respectively, 17.6, 1.4 and 41.8 hrs after the initial GW alerts. We covered 7.8, 3.0 and 5.1\% of the total probability within the GW event localizations, reaching 3$\sigma$ limits of 19.8, 21.3 and 20.8 AB mag, respectively. Although no viable counterparts associated with these events were found, we recovered 6 known transients and ruled out 5 potential candidates. We also present Large Binocular Telescope spectroscopy of PS19eq/SN2019ebq, a promising kilonova candidate that was later determined to be a supernova. With the ability to tile large areas and conduct detailed follow-up, SAGUARO represents a significant addition to GW counterpart searches.

Hubble Space Telescope observations show that low-mass ($M_*=10^9-10^{10}M_{\odot}$) galaxies at high redshift ($z=1.0-2.5$) tend to be elongated (prolate) rather than disky (oblate) or spheroidal. This is explained in zoom-in cosmological hydrodynamical simulations by the fact that these galaxies are forming in cosmic web filaments where accretion happens preferentially along the direction of elongation. We ask whether the elongated morphology of these galaxies allows them to be used as effective tracers of cosmic web filaments at high redshift via their intrinsic alignments. Using mock lightcones and spectroscopically-confirmed galaxy pairs from the CANDELS survey, we test two types of alignments: (1) between the galaxy major axis and the direction to nearby galaxies of any mass, and (2) between the major axes of nearby pairs of low-mass, likely prolate, galaxies. The mock lightcones predict strong signals in 3D real space, 3D redshift space, and 2D projected redshift space for both types of alignments (assuming prolate galaxy orientations are the same as those of their host prolate halos), but we do not detect significant alignment signals in CANDELS observations. However, we show that spectroscopic redshifts have been obtained for only a small fraction of highly elongated galaxies, and accounting for spectroscopic incompleteness and redshift errors significantly degrades the 2D mock signal. This may partly explain the alignment discrepancy and highlights one of several avenues for future work.

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.

Galaxy clustering measurements can be used to constrain many aspects of galaxy evolution, including galaxy host halo masses, satellite quenching efficiencies, and merger rates. We simulate JWST galaxy clustering measurements at z$\sim$4$-$10 by utilizing mock galaxy samples produced by an empirical model, the UniverseMachine. We also adopt the survey footprints and typical depths of the planned joint NIRCam and NIRSpec Guaranteed Time Observation program planned for Cycle 1 to generate realistic JWST survey realizations and to model high-redshift galaxy selection completeness. We find that galaxy clustering will be measured with $\gtrsim$5$\sigma$ significance at z$\sim$4$-$10. Halo mass precisions resulting from Cycle 1 angular clustering measurements will be $\sim$0.2 dex for faint (-18 $\gtrsim$ $\mathit{M}_{\mathrm{UV}}^{ }$ $\gtrsim$ -19) galaxies at z$\sim$4$-$10 as well as $\sim$0.3 dex for bright ($\mathit{M}_{\mathrm{UV}}^{ }$ $\sim$ -20) galaxies at z$\sim$4$-$7. Dedicated spectroscopic follow-up over $\sim$150 arcmin$^2$ would improve these precisions by $\sim$0.1 dex by removing chance projections and low-redshift contaminants. Future JWST observations will therefore provide the first constraints on the stellar-halo mass relation in the epoch of reionization and substantially clarify how this relation evolves at z$>$4. We also find that $\sim$1000 individual satellites will be identifiable at z$\sim$4$-$8 with JWST, enabling strong tests of satellite quenching evolution beyond currently available data (z$\lesssim$2). Finally, we find that JWST observations can measure the evolution of galaxy major merger pair fractions at z$\sim$4$-$8 with $\sim$0.1$-$0.2 dex uncertainties. Such measurements would help determine the relative role of mergers to the build-up of stellar mass into the epoch of reionization.

We measure differences in the emission line flux from the circumgalactic medium, CGM, of galaxies in different environments. Such differences could be a critical clue in explaining a range of galaxy properties that depend on environment. Using large samples of stacked archival spectra from the Sloan Digital Sky Survey, we find that the H$\alpha$ + [N {\small II}] emission line flux from the CGM within 50 kpc of $\sim$ L$^*$ galaxies is lower both for galaxies that lie within a projected distance of $\sim$ 500 kpc from a massive ($M_*>10^{11} M_\odot$) galaxy and for galaxies in richer/denser environments. The environmental differences are statistically significant even after we control for galaxy mass and morphology. We interpret these observations as a direct signature of environmentally-caused strangulation. We present a simple, heuristic model for the effect of a massive parent galaxy. In this model, the CGM cool gas fraction within 50 kpc is significantly decreased for galaxies that lie within 700 kpc of a massive galaxy, with about 80\% of the cool gas removed even when the galaxy is at a distance of 500 kpc from its massive parent. However, we discuss alternative physical causes for the observed behavior and discuss ways forward in addressing open questions.

We use the UniverseMachine to analyze the source of scatter between the central galaxy mass, the total stellar mass in the halo, and the dark matter halo mass. We also propose a new halo mass estimator, the cen+N mass: the sum of the stellar mass of the central and the N most massive satellites. We show that, when real space positions are perfectly known, the cen+N mass has scatter competitive with that of richness-based estimators. However, in redshift space, the cen+N mass suffers less from projection effects in the UniverseMachine model. The cen+N mass is therefore a viable low scatter halo mass estimator, and should be considered an important tool to constrain cosmology with upcoming spectroscopic data from DESI. We analyze the scatter in stellar mass at fixed halo mass and show that the total stellar mass in a halo is uncorrelated with secondary halo properties, but that the central stellar mass is a function of both halo mass and halo age. This is because central galaxies in older halos have had more time to grow via accretion. If the UniverseMachine model is correct, accurate galaxy-halo modeling of mass selected samples therefore needs to consider halo age in addition to mass.

We explore scenarios for reionizing the intergalactic medium with low galaxy ionizing photon escape fractions. We combine simulation-based halo-mass dependent escape fractions with an extrapolation of the observed galaxy rest-ultraviolet luminosity functions to solve for the reionization history from z=20 to z=4. We explore the posterior distributions for key unknown quantities, including the limiting halo mass for star-formation, the ionizing photon production efficiency, and a potential contribution from active galactic nuclei (AGN). We marginalize over the allowable parameter space using a Markov Chain Monte Carlo method, finding a solution which satisfies the most model-independent constraints on reionization. Our fiducial model can match observational constraints with an average escape fraction of <5% throughout the bulk of the epoch of reionization if: i) galaxies form stars down to the atomic cooling limit before reionization and a photosuppression mass of log(M_h/Msol)~9 during/after reionization (-13<M_UV,lim<-11); ii) galaxies become more efficient producers of ionizing photons at higher redshifts and fainter magnitudes, and iii) there is a significant, but sub-dominant, contribution by AGN at z < 7. In this model the faintest galaxies (M_UV>-15) dominate the ionizing emissivity, leading to an earlier start to reionization and a smoother evolution of the ionized volume filling fraction than models which assume a single escape fraction at all redshifts and luminosities. The ionizing emissivity from this model is consistent with observations at z=4-5 (and below, when extrapolated), in contrast to some models which assume a single escape fraction. Our predicted ionized volume filling fraction at z=7 of Q_HII=78% (+\- 8%) is in ~1-2 sigma tension with observations of Lya emitters at z~7 and the damping wing analyses of the two known z>7 quasars, which prefer Q_HII,z=7~40-50%.

The circumgalactic medium (CGM), which harbors > 50% of all the baryons in a galaxy, is both the reservoir of gas for subsequent star formation and the depository of chemically processed gas, energy, and angular momentum from feedback. As such, the CGM obviously plays a critical role in galaxy evolution. We discuss the opportunity to image this component using recombination line emission, beginning with the early results coming from recent statistical detection of this emission to the final goal of realizing spectral-line images of the CGM in individual nearby galaxies. Such work will happen in the next decade and provide new insights on the galactic baryon cycle.

Galaxies form and evolve in the context of their local and large-scale environments. Their baryonic content that we observe with imaging and spectroscopy is intimately connected to the properties of their dark matter halos, and to their location in the "cosmic web" of large-scale structure. Very large spectroscopic surveys of the local universe (e.g., SDSS and GAMA) measure galaxy positions (location within large-scale structure), statistical clustering (a direct constraint on dark matter halo masses), and spectral features (measuring physical conditions of the gas and stars within the galaxies, as well as internal velocities). Deep surveys with the James Webb Space Telescope (JWST) will revolutionize spectroscopic measurements of redshifts and spectral properties for galaxies out to the epoch of reionization, but with numerical statistics and over cosmic volumes that are too small to map large-scale structure and to constrain halo properties via clustering. Here, we consider advances in understanding galaxy evolution that would be enabled by very large spectroscopic surveys at high redshifts: very large numbers of galaxies (outstanding statistics) over large co-moving volumes (large-scale structure on all scales) over broad redshift ranges (evolution over most of cosmic history). The required observational facility can be established as part of the probe portfolio by NASA within the next decade.

Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution.

The nature of dark energy remains a profound mystery 20 years after the discovery of cosmic acceleration. A very high number density galaxy redshift survey over a wide area (HD GRS Wide) spanning the redshift range of 0.5<z<4 using the same tracer, carried out using massively parallel wide field multi-object slit spectroscopy from space, will provide definitive dark energy measurements with minimal observational systematics by design. The HD GRS Wide will illuminate the nature of dark energy, and lead to revolutionary advances in particle physics and cosmology. It will also trace the cosmic web of dark matter and provide key insight into large-scale structure in the Universe. The required observational facility can be established as part of the probe portfolio by NASA within the next decade.

Over the past decade, empirical constraints on the galaxy-dark matter halo connection have significantly advanced our understanding of galaxy evolution. Past techniques have focused on connections between halo properties and galaxy stellar mass and/or star formation rates. Empirical techniques in the next decade will link halo assembly histories with galaxies' circumgalactic media, supermassive black holes, morphologies, kinematics, sizes, colors, metallicities, and transient rates. Uncovering these links will resolve many critical uncertainties in galaxy formation and will enable much higher-fidelity mock catalogs essential for interpreting observations. Achieving these results will require broader and deeper spectroscopic coverage of galaxies and their circumgalactic media; survey teams will also need to meet several criteria (cross-comparisons, public access, and covariance matrices) to facilitate combining data across different surveys. Acting on these recommendations will continue enabling dramatic progress in both empirical modeling and galaxy evolution for the next decade.

ATLAS (Astrophysics Telescope for Large Area Spectroscopy) Probe is a concept for a NASA probe-class space mission. It is the follow-up space mission to WFIRST, boosting its scientific return by obtaining deep IR slit spectroscopy for 70% of all galaxies imaged by a 2000 sq deg WFIRST High Latitude Survey at z>0.5. ATLAS will measure accurate and precise redshifts for 200M galaxies out to z < 7, and deliver spectra that enable a wide range of diagnostic studies of the physical properties of galaxies over most of cosmic history. ATLAS Probe science spans four broad categories: (1) Revolutionizing galaxy evolution studies by tracing the relation between galaxies and dark matter from galaxy groups to cosmic voids and filaments, from the epoch of reionization through the peak era of galaxy assembly; (2) Opening a new window into the dark Universe by weighing the dark matter filaments using 3D weak lensing with spectroscopic redshifts, and obtaining definitive measurements of dark energy and modification of General Relativity using galaxy clustering; (3) Probing the Milky Way's dust-enshrouded regions, reaching the far side of our Galaxy; and (4) Exploring the formation history of the outer Solar System by characterizing Kuiper Belt Objects. ATLAS Probe is a 1.5m telescope with a field of view of 0.4 sq deg, and uses Digital Micro-mirror Devices (DMDs) as slit selectors. It has a spectroscopic resolution of R = 1000 over 1-4 microns, and a spectroscopic multiplex factor >5,000. ATLAS is designed to fit within the NASA probe-class space mission cost envelope; it has a single instrument, a telescope aperture that allows for a lighter launch vehicle, and mature technology. ATLAS Probe will lead to transformative science over the entire range of astrophysics: from galaxy evolution to the dark Universe, from Solar System objects to the dusty regions of the Milky Way.

Using deep images from the Hyper Suprime-Cam (HSC) survey and taking advantage of its unprecedented weak lensing capabilities, we reveal a remarkably tight connection between the stellar mass distribution of massive central galaxies and their host dark matter halo mass. Massive galaxies with more extended stellar mass distributions tend to live in more massive dark matter haloes. We explain this connection with a phenomenological model that assumes, (1) a tight relation between the halo mass and the total stellar content in the halo, (2) that the fraction of in-situ and ex-situ mass at $r<10$ kpc depends on halo mass. This model provides an excellent description of the stellar mass functions (SMF) of total stellar mass ($M_{\star}^{\rm Max}$) and stellar mass within inner 10 kpc ($M_{\star}^{10}$) and also reproduces the HSC weak lensing signals of massive galaxies with different stellar mass distributions. The best-fit model shows that halo mass varies significantly at fixed total stellar mass (as much as 0.4 dex) with a clear dependence on $M_{\star}^{10}$. Our two-parameter $M_{\star}^{\rm Max}$-$M_{\star}^{10}$ description provides a more accurate picture of the galaxy-halo connection at the high-mass end than the simple stellar-halo mass relation (SHMR) and opens a new window to connect the assembly history of halos with those of central galaxies. The model also predicts that the ex-situ component dominates the mass profiles of galaxies at $r< 10$ kpc for $\log M_{\star} \ge 11.7$). The code used for this paper is available online: https://github.com/dr-guangtou/asap

We find significantly different diagnostic emission line ratios for the circumgalactic gas associated with galaxies of stellar masses above and below $10^{10.4}$ M$_\odot$ using SDSS spectroscopy. Specifically, in a sample of 17,393 galaxies, intersected by 18,535 lines of sight at projected radii between 10 and 50 kpc, we stack measured fluxes for nebular strong emission lines, [O {\small III}] $\lambda$5007, H$\alpha$ and [N {\small II}] $\lambda6583$, and find that the gas surrounding the lower mass galaxies exhibits similar line ratios to those of gas ionized by star formation and that surrounding the higher mass galaxies similar to those of gas ionized by AGN or shocks. This finding highlights yet another characteristic of galaxies that is distinctly different above and below this stellar mass threshold, but one that is more closely connected to the gas accretion behavior hypothesized to be responsible for this dichotomy.

The relation between galaxies and dark matter halos reflects the combined effects of many distinct physical processes. Observations indicate that the $z=0$ stellar mass-halo mass (SMHM) relation has remarkably small scatter in stellar mass at fixed halo mass ($\lesssim$ 0.2 dex) with little dependence on halo mass. We investigate the origins of this scatter by combining N-body simulations with observational constraints on the SMHM relation. We find that at the group and cluster scale ($M_{\rm vir}>10^{14}{M_\odot}$) the scatter due purely to hierarchical assembly is $\approx0.16$ dex, which is comparable to recent direct observational estimates. At lower masses, mass buildup since $z\approx2$ is driven largely by in-situ growth. We include a model for the in-situ buildup of stellar mass and find that an intrinsic scatter in this growth channel of $0.2$ dex produces a relation between scatter and halo mass that is consistent with observations from $10^{12}{M_\odot}<M_{\rm vir}<10^{14.75}{M_\odot}$. The approximately constant scatter across a wide range of halo masses at $z=0$ thus appears to be a coincidence as it is determined largely by in-situ growth at low masses and by hierarchical assembly at high masses. These results indicate that the scatter in the SMHM relation can provide unique insight into the regularity of the galaxy formation process.

We simulate recoiling black hole trajectories from $z=20$ to $z=0$ in dark matter halos, quantifying how parameter choices affect escape velocities. These choices include the strength of dynamical friction, the presence of stars and gas, the accelerating expansion of the universe (Hubble acceleration), host halo accretion and motion, and seed black hole mass. $\Lambda$CDM halo accretion increases escape velocities by up to 0.6 dex and significantly shortens return timescales compared to non-accreting cases. Other parameters change orbit damping rates but have subdominant effects on escape velocities; dynamical friction is weak at halo escape velocities, even for extreme parameter values. We present formulae for black hole escape velocities as a function of host halo mass and redshift. Finally, we discuss how these findings affect black hole mass assembly as well as minimum stellar and halo masses necessary to retain supermassive black holes.

Given a galaxy's stellar mass, its host halo mass has a lower limit from the cosmic baryon fraction and known baryonic physics. At z>4, galaxy stellar mass functions place lower limits on halo number densities that approach expected $\Lambda$CDM halo mass functions. High-redshift galaxy stellar mass functions can thus place interesting limits on number densities of massive haloes, which are otherwise very difficult to measure. Although halo mass functions at z<8 are consistent with observed galaxy stellar masses if galaxy baryonic conversion efficiencies increase with redshift, JWST and WFIRST will more than double the redshift range over which useful constraints are available. We calculate maximum galaxy stellar masses as a function of redshift given expected halo number densities from $\Lambda$CDM. We apply similar arguments to black holes. If their virial mass estimates are accurate, number density constraints alone suggest that the quasars SDSS J1044-0125 and SDSS J010013.02+280225.8 likely have black hole mass -- stellar mass ratios higher than the median z=0 relation, confirming the expectation from Lauer bias. Finally, we present a public code to evaluate the probability of an apparently $\Lambda$CDM-inconsistent high-mass halo being detected given the combined effects of multiple surveys and observational errors.

Using a sample of nearly half a million galaxies, intersected by over 8 million lines of sight from the Sloan Digital Sky Survey Data Release 12, we extend our previous study of the recombination radiation emitted by the gaseous halos of nearby galaxies. We identify an inflection in the radial profile of the H$\alpha$+N[{\small II}] radial emission profile at a projected radius of $\sim 50$ kpc and suggest that beyond this radius the emission from ionized gas in spatially correlated halos dominates the profile. We confirm that this is a viable hypothesis using results from a highly simplified theoretical treatment in which the dark matter halo distribution from cosmological simulations is straightforwardly populated with gas. Whether we fit the fraction of halo gas in a cooler (T $= 12,000$ K), smooth ($c = 1$) component (0.26 for galaxies with M$_* = 10^{10.88}$ M$_\odot$ and 0.34 for those with M$_* = 10^{10.18}$ M$_\odot$) or take independent values of this fraction from published hydrodynamical simulations (0.19 and 0.38, respectively), this model successfully reproduces the radial location and amplitude of the observed inflection. We also observe that the physical nature of the gaseous halo connects to primary galaxy morphology beyond any relationship to the galaxy's stellar mass and star formation rate. We explore whether the model reproduces behavior related to the central galaxy's stellar mass, star formation rate, and morphology. We find that it is unsuccessful in reproducing the observations at this level of detail and discuss various shortcomings of our simple model that may be responsible.

This is the first in a series of papers examining the demographics of star-forming galaxies at $0.2<z<2.5$ in CANDELS. We study 9,100 galaxies from GOODS-S and UDS having published values of redshifts, masses, star-formation rates (SFRs), and dust attenuation ($A_V$) derived from UV-optical SED fitting. In agreement with previous works, we find that the $UVJ$ colors of a galaxy are closely correlated with its specific star-formation rate (SSFR) and $A_V$. We define rotated $UVJ$ coordinate axes, termed $S_\mathrm{SED}$ and $C_\mathrm{SED}$, that are parallel and perpendicular to the star-forming sequence and derive a quantitative calibration that predicts SSFR from $C_\mathrm{SED}$ with an accuracy of ~0.2 dex. SFRs from UV-optical fitting and from UV+IR values based on Spitzer/MIPS 24 $\mu\mathrm{m}$ agree well overall, but systematic differences of order 0.2 dex exist at high and low redshifts. A novel plotting scheme conveys the evolution of multiple galaxy properties simultaneously, and dust growth, as well as star-formation decline and quenching, exhibit "mass-accelerated evolution" ("downsizing"). A population of transition galaxies below the star-forming main sequence is identified. These objects are located between star-forming and quiescent galaxies in $UVJ$ space and have lower $A_V$ and smaller radii than galaxies on the main sequence. Their properties are consistent with their being in transit between the two regions. The relative numbers of quenched, transition, and star-forming galaxies are given as a function of mass and redshift.

We explore the phenomenon commonly known as halo assembly bias, whereby dark matter halos of the same mass are found to be more or less clustered when a second halo property is considered, for halos in the mass range $3.7 \times 10^{11} \; h^{-1} \mathrm{M_{\odot}} - 5.0 \times 10^{13} \; h^{-1} \mathrm{M_{\odot}}$. Using the Large Suite of Dark Matter Simulations (LasDamas) we consider nine commonly used halo properties and find that a clustering bias exists if halos are binned by mass or by any other halo property. This secondary bias implies that no single halo property encompasses all the spatial clustering information of the halo population. The mean values of some halo properties depend on their halo's distance to a more massive neighbor. Halo samples selected by having high values of one of these properties therefore inherit a neighbor bias such that they are much more likely to be close to a much more massive neighbor. This neighbor bias largely accounts for the secondary bias seen in halos binned by mass and split by concentration or age. However, halos binned by other mass-like properties still show a secondary bias even when the neighbor bias is removed. The secondary bias of halos selected by their spin behaves differently than that for other halo properties, suggesting that the origin of the spin bias is different than of other secondary biases.

We investigate the evolution of dust content in galaxies from redshifts z=0 to z=9.5. Using empirically motivated prescriptions, we model galactic-scale properties -- including halo mass, stellar mass, star formation rate, gas mass, and metallicity -- to make predictions for the galactic evolution of dust mass and dust temperature in main sequence galaxies. Our simple analytic model, which predicts that galaxies in the early Universe had greater quantities of dust than their low-redshift counterparts, does a good job at reproducing observed trends between galaxy dust and stellar mass out to z~6. We find that for fixed galaxy stellar mass, the dust temperature increases from z=0 to z=6. Our model forecasts a population of low-mass, high-redshift galaxies with interstellar dust as hot as, or hotter than, their more massive counterparts; but this prediction needs to be constrained by observations. Finally, we make predictions for observing 1.1-mm flux density arising from interstellar dust emission with the Atacama Large Millimeter Array.

We empirically constrain how galaxy size relates to halo virial radius using new measurements of the size- and stellar mass-dependent clustering of galaxies in the Sloan Digital Sky Survey. We find that small galaxies cluster much more strongly than large galaxies of the same stellar mass. The magnitude of this clustering difference increases on small scales, and decreases with increasing stellar mass. Using Halotools to forward model the observations, we test an empirical model in which present-day galaxy size is proportional to the size of the virial radius at the time the halo reached its maximum mass. This simple model reproduces the observed size-dependence of galaxy clustering in striking detail. The success of this model provides strong support for the conclusion that satellite galaxies have smaller sizes relative to central galaxies of the same halo mass. Our findings indicate that satellite size is set prior to the time of infall, and that a remarkably simple, linear size--virial radius relation emerges from the complex physics regulating galaxy size. We make quantitative predictions for future measurements of galaxy-galaxy lensing, including dependence upon size, scale, and stellar mass, and provide a scaling relation of the ratio of mean sizes of satellites and central galaxies as a function of their halo mass that can be used to calibrate hydrodynamical simulations and semi-analytic models.

We study the properties of distinct dark matter halos (i.e., those that are not subhalos) that have a final virial mass $M_{\mathrm{vir}}$ at $z = 0$ less than their peak mass ($M_{\mathrm{peak}}$) in the Bolshoi-Planck cosmological simulation. We identify two primary causes of halo mass loss: relaxation after a major merger and tidal stripping by a massive neighbouring halo. Major mergers initially boost $M_{\mathrm{vir}}$ and typically cause the final halo to become more prolate and less relaxed and to have higher spin and lower NFW concentration. As the halo relaxes, high energy material from the recent merger gradually escapes beyond the virial radius, temporarily resulting in a net negative accretion rate that reduces the halo mass by $5-15\%$ on average. Halos that experience a major merger around $z = 0.4$ typically reach a minimum mass near $z = 0$. Tidal stripping mainly occurs in dense regions, and it causes halos to become less prolate and have lower spins and higher NFW concentrations. Tidally stripped halos often lose a large fraction of their peak mass ($> 20\%$) and most never recover (or even reattain a positive accretion rate). Low mass halos can be strongly affected by both post-merger mass loss and tidal stripping, while high mass halos are predominantly influenced by post-merger mass loss and show few signs of significant tidal stripping.

We present the public release of the MultiDark-Galaxies: three distinct galaxy catalogues derived from one of the Planck cosmology MultiDark simulations (i.e. MDPL2, with a volume of (1 Gpc/$h$)$^{3}$ and mass resolution of $1.5 \times 10^{9} M_{\odot}/h$) by applying the semi-analytic models GALACTICUS, SAG, and SAGE to it. We compare the three models and their conformity with observational data for a selection of fundamental properties of galaxies like stellar mass function, star formation rate, cold gas fractions, and metallicities - noting that they sometimes perform differently reflecting model designs and calibrations. We have further selected galaxy subsamples of the catalogues by number densities in stellar mass, cold gas mass, and star formation rate in order to study the clustering statistics of galaxies. We show that despite different treatment of orphan galaxies, i.e. galaxies that lost their dark-matter host halo due to the finite mass resolution of the N-body simulation or tidal stripping, the clustering signal is comparable, and reproduces the observations in all three models - in particular when selecting samples based upon stellar mass. Our catalogues provide a powerful tool to study galaxy formation within a volume comparable to those probed by on-going and future photometric and redshift surveys. All model data consisting of a range of galaxy properties - including broad-band SDSS magnitudes - are publicly available.

We present the first stable release of Halotools (v0.2), a community-driven Python package designed to build and test models of the galaxy-halo connection. Halotools provides a modular platform for creating mock universes of galaxies starting from a catalog of dark matter halos obtained from a cosmological simulation. The package supports many of the common forms used to describe galaxy-halo models: the halo occupation distribution (HOD), the conditional luminosity function (CLF), abundance matching, and alternatives to these models that include effects such as environmental quenching or variable galaxy assembly bias. Satellite galaxies can be modeled to live in subhalos, or to follow custom number density profiles within their halos, including spatial and/or velocity bias with respect to the dark matter profile. The package has an optimized toolkit to make mock observations on a synthetic galaxy population, including galaxy clustering, galaxy-galaxy lensing, galaxy group identification, RSD multipoles, void statistics, pairwise velocities and others, allowing direct comparison to observations. Halotools is object-oriented, enabling complex models to be built from a set of simple, interchangeable components, including those of your own creation. Halotools has an automated testing suite and is exhaustively documented on http://halotools.readthedocs.io, which includes quickstart guides, source code notes and a large collection of tutorials. The documentation is effectively an online textbook on how to build and study empirical models of galaxy formation with Python.

We present clustering properties from 579,492 Lyman break galaxies (LBGs) at z~4-6 over the 100 deg^2 sky (corresponding to a 1.4 Gpc^3 volume) identified in early data of the Hyper Suprime-Cam (HSC) Subaru strategic program survey. We derive angular correlation functions (ACFs) of the HSC LBGs with unprecedentedly high statistical accuracies at z~4-6, and compare them with the halo occupation distribution (HOD) models. We clearly identify significant ACF excesses in 10"<$\theta$<90", the transition scale between 1- and 2-halo terms, suggestive of the existence of the non-linear halo bias effect. Combining the HOD models and previous clustering measurements of faint LBGs at z~4-7, we investigate dark-matter halo mass (Mh) of the z~4-7 LBGs and its correlation with various physical properties including the star-formation rate (SFR), the stellar-to-halo mass ratio (SHMR), and the dark matter accretion rate (dotMh) over a wide-mass range of Mh/M$_\odot$=4x10^10-4x10^12. We find that the SHMR increases from z~4 to 7 by a factor of ~4 at Mh~1x10^11 M$_\odot$, while the SHMR shows no strong evolution in the similar redshift range at Mh~1x10^12 M$_\odot$. Interestingly, we identify a tight relation of SFR/dotMh-Mh showing no significant evolution beyond 0.15 dex in this wide-mass range over z~4-7. This weak evolution suggests that the SFR/dotMh-Mh relation is a fundamental relation in high-redshift galaxy formation whose star formation activities are regulated by the dark matter mass assembly. Assuming this fundamental relation, we calculate the cosmic SFR densities (SFRDs) over z=0-10 (a.k.a. Madau-Lilly plot). The cosmic SFRD evolution based on the fundamental relation agrees with the one obtained by observations, suggesting that the cosmic SFRD increase from z~10 to 4-2 (decrease from z~4-2 to 0) is mainly driven by the increase of the halo abundance (the decrease of the accretion rate).

We explore observational and theoretical constraints on how galaxies might transition between the "star-forming main sequence" (SFMS) and varying "degrees of quiescence" out to $z=3$. Our analysis is focused on galaxies with stellar mass $M_*>10^{10}M_{\odot}$, and is enabled by GAMA and CANDELS observations, a semi-analytic model (SAM) of galaxy formation, and a cosmological hydrodynamical "zoom in" simulation with momentum-driven AGN feedback. In both the observations and the SAM, transition galaxies tend to have intermediate S\'ersic indices, half-light radii, and surface stellar mass densities compared to star-forming and quiescent galaxies out to $z=3$. We place an observational upper limit on the average population transition timescale as a function of redshift, finding that the average high-redshift galaxy is on a "fast track" for quenching whereas the average low-redshift galaxy is on a "slow track" for quenching. We qualitatively identify four physical origin scenarios for transition galaxies in the SAM: oscillations on the SFMS, slow quenching, fast quenching, and rejuvenation. Quenching timescales in both the SAM and the hydrodynamical simulation are not fast enough to reproduce the quiescent population that we observe at $z\sim3$. In the SAM, we do not find a clear-cut morphological dependence of quenching timescales, but we do predict that the mean stellar ages, cold gas fractions, SMBH masses, and halo masses of transition galaxies tend to be intermediate relative to those of star-forming and quiescent galaxies at $z<3$.

We explore empirical constraints on the statistical relationship between the radial size of galaxies and the radius of their host dark matter halos from $z\sim 0.1$--3 using the GAMA and CANDELS surveys. We map dark matter halo mass to galaxy stellar mass using relationships from abundance matching, applied to the Bolshoi-Planck dissipationless N-body simulation. We define SRHR$\equiv r_e/R_h$ as the ratio of galaxy radius to halo virial radius, and SRHR$\lambda \equiv r_e/(\lambda R_h)$ as the ratio of galaxy radius to halo spin parameter times halo radius. At $z\sim 0.1$, we find an average value of SRHR $\simeq 0.018$ and SRHR$\lambda \simeq 0.5$ with very little dependence on stellar mass. SRHR and SRHR$\lambda$ have a weak dependence on cosmic time since $z\sim 3$. SRHR shows a mild decrease over cosmic time for low mass galaxies, but increases slightly or does not evolve for more massive galaxies. We find hints that at high redshift ($z\sim 2$--3), SRHR$\lambda$ is lower for more massive galaxies, while it shows no significant dependence on stellar mass at $z\lesssim 0.5$. We find that for both the GAMA and CANDELS samples, at all redshifts from $z\sim 0.1$--3, the observed conditional size distribution in stellar mass bins is remarkably similar to the conditional distribution of $\lambda R_h$. We discuss the physical interpretation and implications of these results.

We report and provide fitting functions for the abundance of dark matter halos and subhalos as a function of mass, circular velocity, and redshift from the new Bolshoi-Planck and MultiDark-Planck $\Lambda$CDM cosmological simulations, based on the Planck cosmological parameters. We also report the halo mass accretion rates, which may be connected with galaxy star formation rates. We show that the higher cosmological matter density of the Planck parameters compared with the WMAP parameters leads to higher abundance of massive halos at high redshifts. We find that the median halo spin parameter $\lambda_{\rm B} = J(2M_{\rm vir}R_{\rm vir}V_{\rm vir})^{-1}$ is nearly independent of redshift, leading to predicted evolution of galaxy sizes that is consistent with observations, while the significant decrease with redshift in median $\lambda_{\rm P} = J|E|^{-1/2}G^{-1}M^{-5/2}$ predicts more decrease in galaxy sizes than is observed. Using the Tully-Fisher and Faber-Jackson relations between galaxy velocity and mass, we show that a simple model of how galaxy velocity is related to halo maximum circular velocity leads to increasing overprediction of cosmic stellar mass density as redshift increases beyond redshifts $z\sim1$, implying that such velocity-mass relations must change at redshifts $z>1$. By making a realistic model of how observed galaxy velocities are related to halo circular velocity, we show that recent optical and radio observations of the abundance of galaxies are in good agreement with our $\Lambda$CDM simulations. Our halo demographics are based on updated versions of the \rockstar\ and \ctrees\ codes, and this paper includes appendices explaining all of their outputs. This paper is an introduction to a series of related papers presenting other analyses of the Bolshoi-Planck and MultiDark-Planck simulations.

In current and future surveys, quasars play a key role. The new data will extend our knowledge of the Universe as it will be used to better constrain the cosmological model at redshift $z>1$ via baryon acoustic oscillation and redshift space distortion measurements. Here, we present the first clustering study of quasars observed by the extended Baryon Oscillation Spectroscopic Survey. We measure the clustering of $\sim 70,000$ quasars located in the redshift range $0.9<z<2.2$ that cover 1,168 deg$^2$. We model the clustering and produce high-fidelity quasar mock catalogues based on the BigMultiDark Planck simulation. Thus, we use a modified (Sub)Halo Abundance Matching model to account for the specificities of the halo population hosting quasars. We find that quasars are hosted by halos with masses $\sim10^{12.7}M_\odot$ and their bias evolves from 1.54 ($z=1.06$) to 3.15 ($z=1.98$). Using the current eBOSS data, we cannot distinguish between models with different fractions of satellites. The high-fidelity mock light-cones, including properties of halos hosting quasars, are made publicly available.

We study how properties of discrete dark matter halos depend on halo environment, characterized by the mass density around the halos on scales from 0.5 to 16 $h^{-1}{\rm Mpc}$. We find that low mass halos (those less massive than the characteristic mass $M_{\rm C}$ of halos collapsing at a given epoch) in high-density environments have lower accretion rates, lower spins, higher concentrations, and rounder shapes than halos in median density environments. Halos in median and low-density environments have similar accretion rates and concentrations, but halos in low density environments have lower spins and are more elongated. Halos of a given mass in high-density regions accrete material earlier than halos of the same mass in lower-density regions. All but the most massive halos in high-density regions are losing mass (i.e., being stripped) at low redshifts, which causes artificially lowered NFW scale radii and increased concentrations. Tidal effects are also responsible for the decreasing spins of low mass halos in high density regions at low redshifts $z < 1$, by preferentially removing higher angular momentum material from halos. Halos in low-density regions have lower than average spins because they lack nearby halos whose tidal fields can spin them up. We also show that the simulation density distribution is well fit by an Extreme Value Distribution, and that the density distribution becomes broader with cosmic time.

We present high signal-to-noise galaxy-galaxy lensing measurements of the BOSS CMASS sample using 250 square degrees of weak lensing data from CFHTLenS and CS82. We compare this signal with predictions from mock catalogs trained to match observables including the stellar mass function and the projected and two dimensional clustering of CMASS. We show that the clustering of CMASS, together with standard models of the galaxy-halo connection, robustly predicts a lensing signal that is 20-40% larger than observed. Detailed tests show that our results are robust to a variety of systematic effects. Lowering the value of $S_{\rm 8}=\sigma_{\rm 8} \sqrt{\Omega_{\rm m}/0.3}$ compared to Planck2015 reconciles the lensing with clustering. However, given the scale of our measurement ($r<10$ $h^{-1}$ Mpc), other effects may also be at play and need to be taken into consideration. We explore the impact of baryon physics, assembly bias, massive neutrinos, and modifications to general relativity on $\Delta\Sigma$ and show that several of these effects may be non-negligible given the precision of our measurement. Disentangling cosmological effects from the details of the galaxy-halo connection, the effects of baryons, and massive neutrinos, is the next challenge facing joint lensing and clustering analyses. This is especially true in the context of large galaxy samples from Baryon Acoustic Oscillation surveys with precise measurements but complex selection functions.

We test for galactic conformity at $0.2<z<1.0$ to a projected distance of 5 Mpc using spectroscopic redshifts from the PRism MUlti-object Survey (PRIMUS). Our sample consists of $\sim60,000$ galaxies in five separate fields covering a total of $\sim5.5$ square degrees, which allows us to account for cosmic variance. We identify star-forming and quiescent "isolated primary" (i.e., central) galaxies using isolation criteria and cuts in specific star formation rate. We match the redshift and stellar mass distributions of these samples, to control for correlations between quiescent fraction and redshift and stellar mass. We detect a significant $(>3\sigma)$ one-halo conformity signal, or an excess of star-forming neighbors around star-forming central galaxies, of $\sim5$% on scales of 0-1 Mpc and a $2.5\sigma$ two-halo signal of $\sim1$% on scales of 1-3 Mpc. These signals are weaker than those detected in SDSS and are consistent with galactic conformity being the result of large-scale tidal fields and reflecting assembly bias. We also measure the star-forming fraction of central galaxies at fixed stellar mass as a function of large-scale environment and find that central galaxies are more likely to be quenched in overdense environments, independent of stellar mass. However, we find that environment does not affect the star formation efficiency of central galaxies, as long as they are forming stars. We test for redshift and stellar mass dependence of the conformity signal within our sample and show that large volumes and multiple fields are required at intermediate redshift to adequately account for cosmic variance.

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

In the next few years, intensity-mapping surveys that target lines such as CO, Ly$\alpha$, and CII stand to provide powerful probes of high-redshift astrophysics. However, these line emissions are highly non-Gaussian, and so the typical power-spectrum methods used to study these maps will leave out a significant amount of information. We propose a new statistic, the probability distribution of voxel intensities, which can access this extra information. Using a model of a CO intensity map at $z\sim3$ as an example, we demonstrate that this voxel intensity distribution (VID) provides substantial constraining power beyond what is obtainable from the power spectrum alone. We find that a future survey similar to the planned COMAP Full experiment could constrain the CO luminosity function to order $\sim10\%$. We also explore the effects of contamination from continuum emission, interloper lines, and gravitational lensing on our constraints and find that the VID statistic retains significant constraining power even in pessimistic scenarios.

We present the measurements of the luminosity-dependent redshift-space three-point correlation functions (3PCFs) for the Sloan Digital Sky Survey (SDSS) DR7 Main galaxy sample. We compare the 3PCF measurements to the predictions from three different halo and subhalo models. One is the halo occupation distribution (HOD) model and the other two are extensions of the subhalo abundance matching (SHAM) model by allowing the central and satellite galaxies to have different occupation distributions in the host halos and subhalos. Parameters in all the models are chosen to best describe the projected and redshift-space two-point correlation functions (2PCFs) of the same set of galaxies. All three model predictions agree well with the 3PCF measurements for the most luminous galaxy sample, while the HOD model better performs in matching the 3PCFs of fainter samples (with luminosity threshold below $L^*$), which is similar in trend to the case of fitting the 2PCFs. The decomposition of the model 3PCFs into contributions from different types of galaxy triplets shows that on small scales the dependence of the 3PCFs on triangle shape is driven by nonlinear redshift-space distortion (and not by the intrinsic halo shape) while on large scales it reflects the filamentary structure. The decomposition also reveals more detailed differences in the three models, which are related to the radial distribution, the mean occupation function, and the velocity distribution of satellite galaxies inside halos. The results suggest that galaxy 3PCFs can further help constrain the above galaxy-halo relation and test theoretical models.

We use subhalo abundance matching (SHAM) to model the stellar mass function (SMF) and clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) "CMASS" sample at $z\sim0.5$. We introduce a novel method which accounts for the stellar mass incompleteness of CMASS as a function of redshift, and produce CMASS mock catalogs which include selection effects, reproduce the overall SMF, the projected two-point correlation function $w_{\rm p}$, the CMASS $dn/dz$, and are made publicly available. We study the effects of assembly bias above collapse mass in the context of "age matching" and show that these effects are markedly different compared to the ones explored by Hearin et al. (2013) at lower stellar masses. We construct two models, one in which galaxy color is stochastic ("AbM" model) as well as a model which contains assembly bias effects ("AgM" model). By confronting the redshift dependent clustering of CMASS with the predictions from our model, we argue that that galaxy colors are not a stochastic process in high-mass halos. Our results suggest that the colors of galaxies in high-mass halos are determined by other halo properties besides halo peak velocity and that assembly bias effects play an important role in determining the clustering properties of this sample.

We present a study of the clustering and halo occupation distribution of BOSS CMASS galaxies in the redshift range 0.43 < z < 0.7 drawn from the Final SDSS-III Data Release. We compare the BOSS results with the predictions of a Halo Abundance Matching (HAM) clustering model that assigns galaxies to dark matter halos selected from the large BigMultiDark $N$-body simulation of a flat $\Lambda$CDM Planck cosmology. We compare the observational data with the simulated ones on a light-cone constructed from 20 subsequent outputs of the simulation. Observational effects such as incompleteness, geometry, veto masks and fiber collisions are included in the model, which reproduces within 1-$\sigma$ errors the observed monopole of the 2-point correlation function at all relevant scales: from the smallest scales, 0.5 $h^{-1}$ Mpc, up to scales beyond the Baryonic Acoustic Oscillation feature. This model also agrees remarkably well with the BOSS galaxy power spectrum (up to $k\sim1$ $h$ Mpc$^{-1}$), and the Three-point correlation function. The quadrupole of the correlation function presents some tensions with observations. We discuss possible causes that can explain this disagreement, including target selection effects. Overall, the standard HAM model describes remarkably well the clustering statistics of the CMASS sample. We compare the stellar to halo mass relation for the CMASS sample measured using weak lensing in the CFHT Stripe 82 Survey with the prediction of our clustering model, and find a good agreement within 1-$\sigma$. The BigMD-BOSS light-cone including properties of BOSS galaxies and halo properties is made publicly available.

We model the luminosity-dependent projected and redshift-space two-point correlation functions (2PCFs) of the Sloan Digital Sky Survey (SDSS) DR7 Main galaxy sample, using the halo occupation distribution (HOD) model and the subhalo abundance matching (SHAM) model and its extension. All the models are built on the same high-resolution $N$-body simulations. We find that the HOD model generally provides the best performance in reproducing the clustering measurements in both projected and redshift spaces. The SHAM model with the same halo-galaxy relation for central and satellite galaxies (or distinct haloes and subhaloes), when including scatters, has a best-fitting $\chi^2/\rm{dof}$ around $2$--$3$. We therefore extend the SHAM model to the subhalo clustering and abundance matching (SCAM) by allowing the central and satellite galaxies to have different galaxy--halo relations. We infer the corresponding halo/subhalo parameters by jointly fitting the galaxy 2PCFs and abundances and consider subhaloes selected based on three properties, the mass $M_{\rm acc}$ at the time of accretion, the maximum circular velocity $V_{\rm acc}$ at the time of accretion, and the peak maximum circular velocity $V_{\rm peak}$ over the history of the subhaloes. The three subhalo models work well for luminous galaxy samples (with luminosity above $L_*$). For low-luminosity samples, the $V_{\rm acc}$ model stands out in reproducing the data, with the $V_{\rm peak}$ model slightly worse, while the $M_{\rm acc}$ model fails to fit the data. We discuss the implications of the modeling results.

Merger trees are routinely used to follow the growth and merging history of dark matter haloes and subhaloes in simulations of cosmic structure formation. Srisawat et al. (2013) compared a wide range of merger-tree-building codes. Here we test the influence of output strategies and mass resolution on tree-building. We find that, somewhat surprisingly, building the tree from more snapshots does not generally produce more complete trees; instead, it tends to short- en them. Significant improvements are seen for patching schemes which attempt to bridge over occasional dropouts in the underlying halo catalogues or schemes which combine the halo-finding and tree-building steps seamlessly. The adopted output strategy does not affec- t the average number of branches (bushiness) of the resultant merger trees. However, mass resolution has an influence on both main branch length and the bushiness. As the resolution increases, a halo with the same mass can be traced back further in time and will encounter more small progenitors during its evolutionary history. Given these results, we recommend that, for simulations intended as precursors for galaxy formation models where of order 100 or more snapshots are analysed, the tree-building routine should be integrated with the halo finder, or at the very least be able to patch over multiple adjacent snapshots.

We present post-cryogenic Spitzer imaging at 3.6 and 4.5 micron with the Infrared Array Camera (IRAC) of the Spitzer/HETDEX Exploratory Large-Area (SHELA) survey. SHELA covers $\sim$deg$^2$ of the Sloan Digital Sky Survey "Stripe 82" region, and falls within the footprints of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) and the Dark Energy Survey. The HETDEX blind R $\sim$ 800 spectroscopy will produce $\sim$ 200,000 redshifts from the Lyman-$\alpha$ emission for galaxies in the range 1.9 < z < 3.5, and an additional $\sim$200,000 redshifts from the [OII] emission for galaxies at z < 0.5. When combined with deep ugriz images from the Dark Energy Camera, K-band images from NEWFIRM, and other ancillary data, the IRAC photometry from Spitzer will enable a broad range of scientific studies of the relationship between structure formation, galaxy stellar mass, halo mass, AGN, and environment over a co-moving volume of $\sim$0.5 Gpc$^3$ at 1.9 < z < 3.5. Here, we discuss the properties of the SHELA IRAC dataset, including the data acquisition, reduction, validation, and source catalogs. Our tests show the images and catalogs are 80% (50%) complete to limiting magnitudes of 22.0 (22.6) AB mag in the detection image, which is constructed from the weighted sum of the IRAC 3.6 and 4.5 micron images. The catalogs reach limiting sensitivities of 1.1 $\mu$Jy at both 3.6 and 4.5 micron (1$\sigma$, for R=2 arcsec circular apertures). As a demonstration of science, we present IRAC number counts, examples of highly temporally variable sources, and galaxy surface density profiles of rich galaxy clusters. In the spirit of Spitzer Exploratory programs we provide all images and catalogs as part of the publication.

Recent observations have shown that the characteristic luminosity of the rest-frame ultraviolet (UV) luminosity function does not significantly evolve at 4 < z < 7 and is approximately M*_UV ~ -21. We investigate this apparent non-evolution by examining a sample of 178 bright, M_UV < -21 galaxies at z=4 to 7, analyzing their stellar populations and host halo masses. Including deep Spitzer/IRAC imaging to constrain the rest-frame optical light, we find that M*_UV galaxies at z=4-7 have similar stellar masses of log(M/Msol)=9.6-9.9 and are thus relatively massive for these high redshifts. However, bright galaxies at z=4-7 are less massive and have younger inferred ages than similarly bright galaxies at z=2-3, even though the two populations have similar star formation rates and levels of dust attenuation. We match the abundances of these bright z=4-7 galaxies to halo mass functions from the Bolshoi Lambda-CDM simulation to estimate the halo masses. We find that the typical halo masses in ~M*_UV galaxies decrease from log(M_h/Msol)=11.9 at z=4 to log(M_h/Msol)=11.4 at z=7. Thus, although we are studying galaxies at a similar mass across multiple redshifts, these galaxies live in lower mass halos at higher redshift. The stellar baryon fraction in units of the cosmic mean Omega_b/Omega_m rises from 5.1% at z=4 to 11.7% at z=7; this evolution is significant at the ~3-sigma level. This rise does not agree with simple expectations of how galaxies grow, and implies that some effect, perhaps a diminishing efficiency of feedback, is allowing a higher fraction of available baryons to be converted into stars at high redshifts.

We measure the evolution of the quiescent fraction and quenching efficiency of satellites around star-forming and quiescent central galaxies with stellar mass $\log(M_{\mathrm{cen}}/M_{\odot})>10.5$ at $0.3<z<2.5$. We combine imaging from three deep near-infrared-selected surveys (ZFOURGE/CANDELS, UDS, and UltraVISTA), which allows us to select a stellar-mass complete sample of satellites with $\log(M_{\mathrm{sat}}/M_{\odot})>9.3$. Satellites for both star-forming and quiescent central galaxies have higher quiescent fractions compared to field galaxies matched in stellar mass at all redshifts. We also observe "galactic conformity": satellites around quiescent centrals are more likely to be quenched compared to the satellites around star-forming centrals. In our sample, this conformity signal is significant at $\gtrsim3\sigma$ for $0.6<z<1.6$, whereas it is only weakly significant at $0.3<z<0.6$ and $1.6<z<2.5$. Therefore, conformity (and therefore satellite quenching) has been present for a significant fraction of the age of the universe. The satellite quenching efficiency increases with increasing stellar mass of the central, but does not appear to depend on the stellar mass of the satellite to the mass limit of our sample. When we compare the satellite quenching efficiency of star-forming centrals with stellar masses 0.2 dex higher than quiescent centrals (which should account for any difference in halo mass), the conformity signal decreases, but remains statistically significant at $0.6<z<0.9$. This is evidence that satellite quenching is connected to the star-formation properties of the central as well as to the mass of the halo. We discuss physical effects that may contribute to galactic conformity, and emphasize that they must allow for continued star-formation in the central galaxy even as the satellites are quenched.

The galaxy stellar-to-halo mass relation (SHMR) is nearly time-independent for $z < 4$. We therefore construct a time-independent SHMR model for central galaxies, wherein the in-situ star formation rate (SFR) is determined by the halo mass accretion rate (MAR), which we call Stellar-Halo Accretion Rate Coevolution (SHARC). We show that the $\sim0.3$ dex dispersion of the halo MAR matches the observed dispersion of the SFR on the star-formation main sequence (MS). In the context of "bathtub"-type models of galaxy formation, SHARC leads to mass-dependent constraints on the relation between SFR and MAR. Despite its simplicity and the simplified treatment of mass growth from mergers, the SHARC model is likely to be a good approximation for central galaxies with $M_*=10^9- 10^{10.5}M_\odot$ that are on the MS, representing most of the star formation in the Universe. SHARC predictions agree with observed SFRs for galaxies on the MS at low redshifts, agree fairly well at $z\sim4$, but exceed observations at $z>4$. Assuming that the interstellar gas mass is constant for each galaxy (the "equilibrium condition" in bathtub models), the SHARC model allows calculation of net mass loading factors for inflowing and outflowing gas. With assumptions about preventive feedback based on simulations, SHARC allows calculation of galaxy metallicity evolution. If galaxy SFRs indeed track halo MARs, especially at low redshifts, that may help explain the success of models linking galaxy properties to halos (including age-matching) and the similarities between two-halo galaxy conformity and halo mass accretion conformity.

We present the first comprehensive analysis of the segregation of dark matter subhaloes in their host haloes. Using numerical simulations, we examine the segregation of twelve different subhalo properties with respect to both orbital energy and halo-centric radius (in real space as well as in projection). Subhaloes are strongly segregated by accretion redshift, which is an outcome of the inside-out assembly of their host haloes. Since subhaloes that were accreted earlier have experienced more tidal stripping, subhaloes that have lost a larger fraction of their mass at infall are on more bound orbits. Subhaloes are also strongly segregated in their masses and maximum circular velocities at accretion. We demonstrate that part of this segregation is already imprinted in the infall conditions. For massive subhaloes it is subsequently boosted by dynamical friction, but only during their first radial orbit. The impact of these two effects is counterbalanced, though, by the fact that subhaloes with larger accretion masses are accreted later. Because of tidal stripping, subhaloes reveal little to no segregation by present-day mass or maximum circular velocity, while the corresponding torques cause subhaloes on more bound orbits to have smaller spin. There is a weak tendency for subhaloes that formed earlier to be segregated towards the center of their host halo, which is an indirect consequence of the fact that (sub)halo formation time is correlated with other, strongly segregated properties. We discuss the implications of our results for the segregation of satellite galaxies in galaxy groups and clusters.

We combine constraints on galaxy formation histories with planet formation models, yielding the Earth-like and giant planet formation histories of the Milky Way and the Universe as a whole. In the Hubble Volume (10^13 Mpc^3), we expect there to be ~10^20 Earth-like and ~10^20 giant planets; our own galaxy is expected to host ~10^9 and ~10^10 Earth-like and giant planets, respectively. Proposed metallicity thresholds for planet formation do not significantly affect these numbers. However, the metallicity dependence for giant planets results in later typical formation times and larger host galaxies than for Earth-like planets. The Solar System formed at the median age for existing giant planets in the Milky Way, and consistent with past estimates, formed after 80% of Earth-like planets. However, if existing gas within virialised dark matter haloes continues to collapse and form stars and planets, the Universe will form over 10 times more planets than currently exist. We show that this would imply at least a 92% chance that we are not the only civilisation the Universe will ever have, independent of arguments involving the Drake Equation.

It is often claimed that overdensities of (or even individual bright) submillimetre-selected galaxies (SMGs) trace the assembly of the most-massive dark matter structures in the Universe. We test this claim by performing a counts-in-cells analysis of mock SMG catalogues derived from the Bolshoi cosmological simulation to investigate how well SMG associations trace the underlying dark matter structure. We find that SMGs exhibit a relatively complex bias: some regions of high SMG overdensity are underdense in terms of dark matter mass, and some regions of high dark matter overdensity contain no SMGs. Because of their rarity, Poisson noise causes scatter in the SMG overdensity at fixed dark matter overdensity. Consequently, rich associations of less-luminous, more-abundant galaxies (i.e. Lyman-break galaxy analogues) trace the highest dark matter overdensities much better than SMGs. Even on average, SMG associations are relatively poor tracers of the most significant dark matter overdensities because of 'downsizing': at z < ~2.5, the most-massive galaxies that reside in the highest dark matter overdensities have already had their star formation quenched and are thus no longer SMGs. At a given redshift, of the 10 per cent most-massive overdensities, only ~25 per cent contain at least one SMG, and less than a few per cent contain more than one SMG.

We present a robust measurement and analysis of the rest-frame ultraviolet (UV) luminosity function at z=4-8. We use deep Hubble Space Telescope imaging over the CANDELS/GOODS fields, the Hubble Ultra Deep Field and the Year 1 Hubble Frontier Field deep parallel observations. These surveys provides an effective volume of 0.6-1.2 x 10^6 Mpc^3 over this epoch, allowing us to perform a robust search for faint (M_UV=-18) and bright (M_UV < -21) galaxies. We select candidate galaxies using a well-tested photometric redshift technique with careful screening of contaminants, finding a sample of 7446 galaxies at 3.5<z<8.5, with >1000 galaxies at z~6-8. We measure the luminosity function using a Markov Chain Monte Carlo analysis to measure robust uncertainties. At the faint end our results agree with previous studies, yet we find a higher abundance of UV-bright galaxies at z>6, with M* ~ -21 at z>5, different than that inferred based on previous trends at lower redshift. At z=8, a single power-law provides an equally good fit to the UV luminosity function, while at z=6 and 7, an exponential cutoff at the bright-end is moderately preferred. We compare to semi-analytical models, and find that the lack of evolution in M* is consistent with models where the impact of dust attenuation on the bright-end of the luminosity function decreases at higher redshift. We measure the evolution of the cosmic star-formation rate density, correcting for dust attenuation, and find that it declines as (1+z)^(-4.3 +/- 0.5) at z>4, consistent with observations at z>9. Our observations are consistent with a reionization history that starts at z>10, completes at z>6, and reaches a midpoint (x_HII = 0.5) at 6.7<z<9.4. Finally, our observations predict that the abundance of bright z=9 galaxies is likely higher than previous constraints, though consistent with recent estimates of bright z~10 galaxies. [abridged]

Currently-proposed galaxy quenching mechanisms predict very different behaviours during major halo mergers, ranging from significant quenching enhancement (e.g., clump-induced gravitational heating models) to significant star formation enhancement (e.g., gas starvation models). To test real galaxies' behaviour, we present an observational galaxy pair method for selecting galaxies whose host haloes are preferentially undergoing major mergers. Applying the method to central L* (10^10 Msun < M_* < 10^10.5 Msun) galaxies in the Sloan Digital Sky Survey (SDSS) at z<0.06, we find that major halo mergers can at most modestly reduce the star-forming fraction, from 59% to 47%. Consistent with past research, however, mergers accompany enhanced specific star formation rates for star-forming L* centrals: ~10% when a paired galaxy is within 200 kpc (approximately the host halo's virial radius), climbing to ~70% when a paired galaxy is within 30 kpc. No evidence is seen for even extremely close pairs (<30 kpc separation) rejuvenating star formation in quenched galaxies. For galaxy formation models, our results suggest: (1) quenching in L* galaxies likely begins due to decoupling of the galaxy from existing hot and cold gas reservoirs, rather than a lack of available gas or gravitational heating from infalling clumps, (2) state-of-the-art semi-analytic models currently over-predict the effect of major halo mergers on quenching, and (3) major halo mergers can trigger enhanced star formation in non-quenched central galaxies.

Merging haloes with similar masses (i.e., major mergers) pose significant challenges for halo finders. We compare five halo finding algorithms' (AHF, HBT, Rockstar, SubFind, and VELOCIraptor) recovery of halo properties for both isolated and cosmological major mergers. We find that halo positions and velocities are often robust, but mass biases exist for every technique. The algorithms also show strong disagreement in the prevalence and duration of major mergers, especially at high redshifts (z>1). This raises significant uncertainties for theoretical models that require major mergers for, e.g., galaxy morphology changes, size changes, or black hole growth, as well as for finding Bullet Cluster analogues. All finders not using temporal information also show host halo and subhalo relationship swaps over successive timesteps, requiring careful merger tree construction to avoid problematic mass accretion histories. We suggest that future algorithms should combine phase-space and temporal information to avoid the issues presented.

Correlations between the star formation rates (SFRs) of nearby galaxies (so-called galactic conformity) have been observed for projected separations up to 4 Mpc, an effect not predicted by current semi-analytic models. We investigate correlations between the mass accretion rates (dMvir/dt) of nearby halos as a potential physical origin for this effect. We find that pairs of host halos "know about" each others' assembly histories even when their present-day separation is greater than thirty times the virial radius of either halo. These distances are far too large for direct interaction between the halos to explain the correlation in their dMvir/dt. Instead, halo pairs at these distances reside in the same large-scale tidal environment, which regulates dMvir/dt for both halos. Larger halos are less affected by external forces, which naturally gives rise to a mass dependence of the halo conformity signal. SDSS measurements of galactic conformity exhibit a qualitatively similar dependence on stellar mass, including how the signal varies with distance. Based on the expectation that halo accretion and galaxy SFR are correlated, we predict the scale-, mass- and redshift-dependence of large-scale galactic conformity, finding that the signal should drop to undetectable levels by z > 1. These predictions are testable with current surveys to z ~ 1; confirmation would establish a strong correlation between dark matter halo accretion rate and central galaxy SFR.

This article is the second in a series in which we perform an extensive comparison of various galaxy-based cluster mass estimation techniques that utilise the positions, velocities and colours of galaxies. Our aim is to quantify the scatter, systematic bias and completeness of cluster masses derived from a diverse set of 25 galaxy-based methods using two contrasting mock galaxy catalogues based on a sophisticated halo occupation model and a semi-analytic model. Analysing 968 clusters, we find a wide range in the RMS errors in log M200c delivered by the different methods (0.18 to 1.08 dex, i.e., a factor of ~1.5 to 12), with abundance matching and richness methods providing the best results, irrespective of the input model assumptions. In addition, certain methods produce a significant number of catastrophic cases where the mass is under- or over-estimated by a factor greater than 10. Given the steeply falling high-mass end of the cluster mass function, we recommend that richness or abundance matching-based methods are used in conjunction with these methods as a sanity check for studies selecting high mass clusters. We see a stronger correlation of the recovered to input number of galaxies for both catalogues in comparison with the group/cluster mass, however, this does not guarantee that the correct member galaxies are being selected. We do not observe significantly higher scatter for either mock galaxy catalogues. Our results have implications for cosmological analyses that utilise the masses, richnesses, or abundances of clusters, which have different uncertainties when different methods are used.

We present a semi-empirical model to infer the atomic and molecular hydrogen content of galaxies as a function of halo mass and time. Our model combines the SFR-halo mass-redshift relation (constrained by galaxy abundances) with inverted SFR-surface density relations to infer galaxy H I and H2 masses. We present gas scaling relations, gas fractions, and mass functions from z = 0 to z = 3 and the gas properties of galaxies as a function of their host halo masses. Predictions of our work include: 1) there is a ~ 0.2 dex decrease in the H I mass of galaxies as a function of their stellar mass since z = 1.5, whereas the H2 mass of galaxies decreases by > 1 dex over the same period. 2) galaxy cold gas fractions and H2 fractions decrease with increasing stellar mass and time. Galaxies with M* > 10^10 Msun are dominated by their stellar content at z < 1, whereas less-massive galaxies only reach these gas fractions at z = 0. We find the strongest evolution in relative gas content at z < 1.5. 3) the SFR to gas mass ratio decreases by an order of magnitude from z = 3 to z = 0. This is consistent with lower H2 fractions; these lower fractions in combination with smaller gas reservoirs correspond to decreased present-day galaxy SFRs. 4) an H2-based star- formation relation can simultaneously fuel the evolution of the cosmic star-formation and reproduce the observed weak evolution in the cosmic HI density. 5) galaxies residing in haloes with masses near 10^12 Msun are most efficient at obtaining large gas reservoirs and forming H2 at all redshifts. These two effects lie at the origin of the high star-formation efficiencies in haloes with the same mass.

We show that the ratio of galaxies' specific star formation rates (SSFRs) to their host halos' specific mass accretion rates (SMARs) strongly constrains how the galaxies' stellar masses, specific star formation rates, and host halo masses evolve over cosmic time. This evolutionary constraint provides a simple way to probe z>8 galaxy populations without direct observations. Tests of the method with galaxy properties at z=4 successfully reproduce the known evolution of the stellar mass--halo mass (SMHM) relation, galaxy SSFRs, and the cosmic star formation rate (CSFR) for 5<z<8. We then predict the continued evolution of these properties for 8<z<15. In contrast to the non-evolution in the SMHM relation at z<4, the median galaxy mass at fixed halo mass increases strongly at z>4. We show that this result is closely linked to the flattening in galaxy SSFRs at z>2 compared to halo specific mass accretion rates; we expect that average galaxy SSFRs at fixed stellar mass will continue their mild evolution to z~15. The expected CSFR shows no breaks or features at z>8.5; this constrains both reionization and the possibility of a steep falloff in the CSFR at z=9-10. Finally, we make predictions for stellar mass and luminosity functions for the James Webb Space Telescope (JWST), which should be able to observe one galaxy with M* > ~10^8 Msun per 10^3 Mpc^3 at z=9.6 and one such galaxy per 10^4 Mpc^3 at z=15.

Distant star-forming galaxies show a correlation between their star formation rates (SFR) and stellar masses, and this has deep implications for galaxy formation. Here, we present a study on the evolution of the slope and scatter of the SFR-stellar mass relation for galaxies at $3.5\leq z\leq 6.5$ using multi-wavelength photometry in GOODS-S from the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (CANDELS) and Spitzer Extended Deep Survey. We describe an updated, Bayesian spectral-energy distribution fitting method that incorporates effects of nebular line emission, star formation histories that are constant or rising with time, and different dust attenuation prescriptions (starburst and Small Magellanic Cloud). From $z$=6.5 to $z$=3.5 star-forming galaxies in CANDELS follow a nearly unevolving correlation between stellar mass and SFR that follows SFR $\sim$ $M_\star^a$ with $a = 0.54 \pm 0.16$ at $z\sim 6$ and $0.70 \pm 0.21$ at $z\sim 4$. This evolution requires a star formation history that increases with decreasing redshift (on average, the SFRs of individual galaxies rise with time). The observed scatter in the SFR-stellar mass relation is tight, $\sigma(\log \mathrm{SFR}/\mathrm{M}_\odot$ yr$^{-1})< 0.3\ - $ 0.4 dex, for galaxies with $\log M_\star/\mathrm{M}_\odot > 9$ dex. Assuming that the SFR is tied to the net gas inflow rate (SFR $\sim$ $\dot{M}_\mathrm{gas}$), then the scatter in the gas inflow rate is also smaller than 0.3$-$0.4 dex for star-forming galaxies in these stellar mass and redshift ranges, at least when averaged over the timescale of star formation. We further show that the implied star formation history of objects selected on the basis of their co-moving number densities is consistent with the evolution in the SFR-stellar mass relation.

Galaxies with stellar masses near M* contain the majority of stellar mass in the universe, and are therefore of special interest in the study of galaxy evolution. The Milky Way (MW) and Andromeda (M31) have present day stellar masses near M*, at 5x10^10 Msol (MW-mass) and 10^11 Msol (M31-mass). We study the typical progenitors of these galaxies using ZFOURGE, a deep medium-band near-IR imaging survey, which is sensitive to the progenitors of these galaxies out to z~3. We use abundance-matching techniques to identify the main progenitors of these galaxies at higher redshifts. We measure the evolution in the stellar mass, rest-frame colors, morphologies, far-IR luminosities, and star-formation rates combining our deep multiwavelength imaging with near-IR HST imaging from CANDELS, and far-IR imaging from GOODS-H and CANDELS-H. The typical MW-mass and M31-mass progenitors passed through the same evolution stages, evolving from blue, star-forming disk galaxies at the earliest stages, to redder dust-obscured IR-luminous galaxies in intermediate stages, and to red, more quiescent galaxies at their latest stages. The progenitors of the MW-mass galaxies reached each evolutionary stage at later times (lower redshifts) and with stellar masses that are a factor of 2-3 lower than the progenitors of the M31-mass galaxies. The process driving this evolution, including the suppression of star-formation in present-day M* galaxies requires an evolving stellar-mass/halo-mass ratio and/or evolving halo-mass threshold for quiescent galaxies. The effective size and star-formation rates imply that the baryonic cold-gas fractions drop as galaxies evolve from high redshift to z~0 and are strongly anticorrelated with an increase in the S\'ersic index. Therefore, the growth of galaxy bulges in M* galaxies corresponds to a rapid decline in the galaxy gas fractions and/or a decrease in the star-formation efficiency.

The galaxy stellar mass function (GSMF) at high-z provides key information on star-formation history and mass assembly in the young Universe. We aimed to use the unique combination of deep optical/NIR/MIR imaging provided by HST, Spitzer and the VLT in the CANDELS-UDS, GOODS-South, and HUDF fields to determine the GSMF over the redshift range 3.5<z<7.5. We utilised the HST WFC3/IR NIR imaging from CANDELS and HUDF09, reaching H~27-28.5 over a total area of 369 arcmin2, in combination with associated deep HST ACS optical data, deep Spitzer IRAC imaging from the SEDS programme, and deep Y and K-band VLT Hawk-I images from the HUGS programme, to select a galaxy sample with high-quality photometric redshifts. These have been calibrated with more than 150 spectroscopic redshifts in the range 3.5<z<7.5, resulting in an overall precision of sigma_z/(1+z)~0.037. We have determined the low-mass end of the high-z GSMF with unprecedented precision, reaching down to masses as low as M*~10^9 Msun at z=4 and ~6x10^9 Msun at z=7. We find that the GSMF at 3.5<z<7.5 depends only slightly on the recipes adopted to measure the stellar masses, namely the photo-z, the SFHs, the nebular contribution or the presence of AGN on the parent sample. The low-mass end of the GSMF is steeper than has been found at lower redshifts, but appears to be unchanged over the redshift range probed here. Our results are very different from previous GSMF estimates based on converting UV galaxy luminosity functions into mass functions via tight M/L relations. Integrating our evolving GSMF over mass, we find that the growth of stellar mass density is barely consistent with the time-integral of the SFR density over cosmic time at z>4. These results confirm the unique synergy of the CANDELS+HUDF, HUGS, and SEDS surveys for the discovery and study of moderate/low-mass galaxies at high redshifts.

We use high resolution Hydro$+$N-Body cosmological simulations to compare the assembly and evolution of a small field dwarf (stellar mass ~ 10$^{6-7}$ M$\odot$, total mass 10$^{10}$ M$\odot$ in $\Lambda$ dominated CDM and 2keV WDM cosmologies. We find that star formation (SF) in the WDM model is reduced and delayed by 1-2 Gyr relative to the CDM model, independently of the details of SF and feedback. Independent of the DM model, but proportionally to the SF efficiency, gas outflows lower the central mass density through `dynamical heating', such that all realizations have circular velocities $<$ 20kms at 500$~$pc, in agreement with local kinematic constraints. As a result of dynamical heating, older stars are less centrally concentrated than younger stars, similar to stellar population gradients observed in nearby dwarf galaxies. Introducing an important diagnostic of SF and feedback models, we translate our simulations into artificial color-magnitude diagrams and star formation histories in order to directly compare to available observations. The simulated galaxies formed most of their stars in many $\sim$10 Myr long bursts. The CDM galaxy has a global SFH, HI abundance and Fe/H and alpha-elements distribution well matched to current observations of dwarf galaxies. These results highlight the importance of directly including `baryon physics' in simulations when 1) comparing predictions of galaxy formation models with the kinematics and number density of local dwarf galaxies and 2) differentiating between CDM and non-standard models with different DM or power spectra.

The stellar mass-halo mass relation is a key constraint in all semi-analytic, numerical, and semi-empirical models of galaxy formation and evolution. However, its exact shape and redshift dependence remain debated. Several recent works support a relation in the local Universe steeper than previously thought. Based on the comparisons with a variety of data on massive central galaxies, we show that this steepening holds up to z~1, for stellar masses Mstar>2e11 Msun. Specifically, we find significant evidence for a high-mass end slope of \beta>0.35-0.70, instead of the usual \beta~0.20-0.30 reported by a number of previous results. When including the independent constraints from the recent BOSS clustering measurements, the data, independent of any systematic errors in stellar masses, tend to favor a model with a very small scatter (< 0.15 dex) in stellar mass at fixed halo mass, in the redshift range z < 0.8 and for Mstar>3e11 Msun, suggesting a close connection between massive galaxies and host halos even at relatively recent epochs. We discuss the implications of our results with respect to the evolution of the most massive galaxies since z~1.

In this paper, we test the age matching hypothesis that the star formation rate (SFR) of a galaxy of fixed stellar mass is determined by its dark matter halo formation history, and as such, that more quiescent galaxies reside in older halos. This simple model has been remarkably successful at predicting color-based galaxy statistics at low redshift as measured in the Sloan Digital Sky Survey (SDSS). To further test this method with observations, we present new SDSS measurements of the galaxy two-point correlation function and galaxy-galaxy lensing as a function of stellar mass and SFR, separated into quenched and star-forming galaxy samples. We find that our age matching model is in excellent agreement with these new measurements. We also employ a galaxy group finder and show that our model is able to predict: (1) the relative SFRs of central and satellite galaxies, (2) the SFR-dependence of the radial distribution of satellite galaxy populations within galaxy groups, rich groups, and clusters and their surrounding larger scale environments, and (3) the interesting feature that the satellite quenched fraction as a function of projected radial distance from the central galaxy exhibits an ~ r^-.15 slope, independent of environment. The accurate prediction for the spatial distribution of satellites is intriguing given the fact that we do not explicitly model satellite-specific processes after infall, and that in our model the virial radius does not mark a special transition region in the evolution of a satellite, contrary to most galaxy evolution models. The success of the model suggests that present-day galaxy SFR is strongly correlated with halo mass assembly history.

A halo merger tree forms the essential backbone of a semi-analytic model for galaxy formation and evolution. Recent studies have pointed out that extracting merger trees from numerical simulations of structure formation is non-trivial; different tree building algorithms can give differing merger histories. These differences should be carefully understood before merger trees are used as input for models of galaxy formation. We investigate the impact of different halo merger trees on a semi-analytic model. We find that the z=0 galaxy properties in our model show differences between trees when using a common parameter set. The star formation history of the Universe and the properties of satellite galaxies can show marked differences between trees with different construction methods. Independently calibrating the semi-analytic model for each tree can reduce the discrepancies between the z=0 global galaxy properties, at the cost of increasing the differences in the evolutionary histories of galaxies. Furthermore, the underlying physics implied can vary, resulting in key quantities such as the supernova feedback efficiency differing by factors of 2. Such a change alters the regimes where star formation is primarily suppressed by supernovae. Therefore, halo merger trees extracted from a common halo catalogue using different, but reliable, algorithms can result in a difference in the semi-analytic model. Given the uncertainties in galaxy formation physics, however, these differences may not necessarily be viewed as significant.

We study the shapes of subhalo distributions from four dark-matter-only simulations of Milky Way type haloes. Comparing the shapes derived from the subhalo distributions at high resolution to those of the underlying dark matter fields we find the former to be more triaxial if theanalysis is restricted to massive subhaloes. For three of the four analysed haloes the increased triaxiality of the distributions of massive subhaloes can be explained by a systematic effect caused by the low number of objects. Subhaloes of the fourth halo show indications for anisotropic accretion via their strong triaxial distribution and orbit alignment with respect to the dark matter field. These results are independent of the employed subhalo finder. Comparing the shape of the observed Milky Way satellite distribution to those of high-resolution subhalo samples from simulations, we find an agreement for samples of bright satellites, but significant deviations if faint satellites are included in the analysis. These deviations might result from observational incompleteness.

Nearly 20% of short gamma-ray bursts (sGRBs) have no observed host galaxies. Combining this finding with constraints on galaxies' dark matter halo potential wells gives strong limits on the natal kick velocity distribution for sGRB progenitors. For the best-fitting velocity distribution, one in five sGRB progenitors receives a natal kick above 150 km/s, consistent with merging neutron star models but not with merging white dwarf binary models. This progenitor model constraint is robust to a wide variety of systematic uncertainties, including the sGRB progenitor time-delay model, the Swift redshift sensitivity, and the shape of the natal kick velocity distribution. We also use constraints on the galaxy-halo connection to determine the host halo and host galaxy demographics for sGRBs, which match extremely well with available data. Most sGRBs are expected to occur in halos near 10^12 Msun and in galaxies near 5x10^10 Msun (L_*); unobserved faint and high-redshift host galaxies contribute a small minority of the observed hostless sGRB fraction. We find that sGRB redshift distributions and host galaxy stellar masses weakly constrain the progenitor time-delay model; the active vs. passive fraction of sGRB host galaxies may offer a stronger constraint. Finally, we discuss how searches for gravitational wave optical counterparts in the local Universe can reduce followup times using these findings.

Spectroscopic + photometric redshifts, stellar mass estimates, and rest-frame colors from the 3D-HST survey are combined with structural parameter measurements from CANDELS imaging to determine the galaxy size-mass distribution over the redshift range 0<z<3. Separating early- and late-type galaxies on the basis of star-formation activity, we confirm that early-type galaxies are on average smaller than late-type galaxies at all redshifts, and find a significantly different rate of average size evolution at fixed galaxy mass, with fast evolution for the early-type population, R_eff ~ (1+z)^-1.48, and moderate evolution for the late-type population, R_eff ~ (1+z)^-0.75. The large sample size and dynamic range in both galaxy mass and redshift, in combination with the high fidelity of our measurements due to the extensive use of spectroscopic data, not only fortify previous results, but also enable us to probe beyond simple average galaxy size measurements. At all redshifts the slope of the size-mass relation is shallow, R_eff ~ M_star^0.22, for late-type galaxies with stellar mass >3x10^9 M_sol, and steep, R_eff M_star^0.75, for early-type galaxies with stellar mass >2x10^10 M_sol. The intrinsic scatter is <~0.2 dex for all galaxy types and redshifts. For late-type galaxies, the logarithmic size distribution is not symmetric, but skewed toward small sizes: at all redshifts and masses a tail of small late-type galaxies exists that overlaps in size with the early-type galaxy population. The number density of massive (~10^11 M_sol), compact (R_eff < 2 kpc) early-type galaxies increases from z=3 to z=1.5-2 and then strongly decreases at later cosmic times.

Merger tree codes are routinely used to follow the growth and merger of dark matter haloes in simulations of cosmic structure formation. Whereas in Srisawat et. al. we compared the trees built using a wide variety of such codes here we study the influence of the underlying halo catalogue upon the resulting trees. We observe that the specifics of halo finding itself greatly influences the constructed merger trees. We find that the choices made to define the halo mass are of prime importance. For instance, amongst many potential options different finders select self-bound objects or spherical regions of defined overdensity, decide whether or not to include substructures within the mass returned and vary in their initial particle selection. The impact of these decisions is seen in tree length (the period of time a particularly halo can be traced back through the simulation), branching ratio (essentially the merger rate of \subhalos) and mass evolution. We therefore conclude that the choice of the underlying halo finder is more relevant to the process of building merger trees than the tree builder itself. We also report on some built-in features of specific merger tree codes that (sometimes) help to improve the quality of the merger trees produced.

We present a study of the substructure finder dependence of subhalo clustering in the Aquarius Simulation. We run 11 different subhalo finders on the haloes of the Aquarius Simulation and we study their differences in the density profile, mass fraction and 2-point correlation function of subhaloes in haloes. We also study the mass and vmax dependence of subhalo clustering. As the Aquarius Simulation has been run at different resolutions, we study the convergence with higher resolutions. We find that the agreement between finders is at around the 10% level inside R200 and at intermediate resolutions when a mass threshold is applied, and better than 5% when vmax is restricted instead of mass. However, some discrepancies appear in the highest resolution, underlined by an observed resolution dependence of subhalo clustering. This dependence is stronger for the smallest subhaloes, which are more clustered in the highest resolution, due to the detection of subhaloes within subhaloes (the sub-subhalo term). This effect modifies the mass dependence of clustering in the highest resolutions. We discuss implications of our results for models of subhalo clustering and their relation with galaxy clustering.

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<z<5 will discriminate between models, the discrepancies between the models and existing data of these observables have already revealed challenging problems in understanding star formation and its feedback in galaxy formation. The three sets of models are being used to construct catalogs of mock galaxies on light cones that have the same geometry as the CANDELS survey, which should be particularly useful for quantifying the biases and uncertainties on measurements and inferences from the real observations. -ABRIDGED

Comparing galaxies across redshifts at fixed cumulative number density is a popular way to estimate the evolution of specific galaxy populations. This method ignores scatter in mass accretion histories and galaxy-galaxy mergers, which can lead to errors when comparing galaxies over large redshift ranges (Delta z > 1). We use abundance matching in the LCDM paradigm to estimate the median change in number density with redshift and provide a simple fit (+0.16 dex per unit Delta z) for progenitors of z = 0 galaxies. We find that galaxy descendants do not evolve in the same way as galaxy progenitors, largely due to scatter in mass accretion histories. We also provide estimates for the 1-sigma range of number densities corresponding to galaxy progenitors and descendants. Finally, we discuss some limits on number density comparisons, which arise due to difficulties measuring physical quantities (e.g., stellar mass) consistently across redshifts. A public tool to calculate number density evolution for galaxies, as well as approximate halo masses, is available online.

Merger trees follow the growth and merger of dark-matter haloes over cosmic history. As well as giving important insights into the growth of cosmic structure in their own right, they provide an essential backbone to semi-analytic models of galaxy formation. This paper is the first in a series to arise from the SUSSING MERGER TREES Workshop in which ten different tree-building algorithms were applied to the same set of halo catalogues and their results compared. Although many of these codes were similar in nature, all algorithms produced distinct results. Our main conclusions are that a useful merger-tree code should possess the following features: (i) the use of particle IDs to match haloes between snapshots; (ii) the ability to skip at least one, and preferably more, snapshots in order to recover subhaloes that are temporarily lost during merging; (iii) the ability to cope with (and ideally smooth out) large, temporary flucuations in halo mass. Finally, to enable different groups to communicate effectively, we defined a common terminology that we used when discussing merger trees and we encourage others to adopt the same language. We also specified a minimal output format to record the results.

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.

We use measurements of the stellar mass function, galaxy clustering, and galaxy-galaxy lensing within the COSMOS survey to constrain the stellar-to-halo mass relation (SHMR) of star forming and quiescent galaxies over the redshift range z=[0.2,1.0]. For massive galaxies, M*>~10^10.6 Msol, our results indicate that star-forming galaxies grow proportionately as fast as their dark matter halos while quiescent galaxies are outpaced by dark matter growth. At lower masses, there is minimal difference in the SHMRs, implying that the majority low-mass quiescent galaxies have only recently been quenched of their star formation. Our analysis also affords a breakdown of all COSMOS galaxies into the relative numbers of central and satellite galaxies for both populations. At z=1, satellite galaxies dominate the red sequence below the knee in the stellar mass function. But the number of quiescent satellites exhibits minimal redshift evolution; all evolution in the red sequence is due to low-mass central galaxies being quenched of their star formation. At M*~10^10 Msol, the fraction of central galaxies on the red sequence increases by a factor of ten over our redshift baseline, while the fraction of quenched satellite galaxies at that mass is constant with redshift. We define a "migration rate" to the red sequence as the time derivative of the passive galaxy abundances. We find that the migration rate of central galaxies to the red sequence increases by nearly an order of magnitude from z=1 to z=0. These results imply that the efficiency of quenching star formation for centrals is increasing with cosmic time, while the mechanisms that quench the star formation of satellite galaxies in groups and clusters is losing efficiency.

The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and Vmax function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -- to a minor extent -- the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general.

With the ever increasing size and complexity of fully self-consistent simulations of galaxy formation within the framework of the cosmic web, the demands upon object finders for these simulations has simultaneously grown. To this extent we initiated the Halo Finder Comparison Project that gathered together all the experts in the field and has so far led to two comparison papers, one for dark matter field haloes (Knebe et al. 2011), and one for dark matter subhaloes (Onions et al. 2012). However, as state-of-the-art simulation codes are perfectly capable of not only following the formation and evolution of dark matter but also account for baryonic physics (e.g. hydrodynamics, star formation, feedback) object finders should also be capable of taking these additional processes into consideration. Here we report on a comparison of codes as applied to the Constrained Local UniversE Simulation (CLUES) of the formation of the Local Group which incorporates much of the physics relevant for galaxy formation. We compare both the properties of the three main galaxies in the simulation (representing the MW, M31, and M33) as well as their satellite populations for a variety of halo finders ranging from phase-space to velocity-space to spherical overdensity based codes, including also a mere baryonic object finder. We obtain agreement amongst codes comparable to (if not better than) our previous comparisons, at least for the total, dark, and stellar components of the objects. However, the diffuse gas content of the haloes shows great disparity, especially for low-mass satellite galaxies. This is primarily due to differences in the treatment of the thermal energy during the unbinding procedure. We acknowledge that the handling of gas in halo finders is something that needs to be dealt with carefully, and the precise treatment may depend sensitively upon the scientific problem being studied.

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.

We investigate unbound dark matter particles in halos by tracing particle trajectories in a simulation run to the far future (a = 100). We find that the traditional sum of kinetic and potential energies is a very poor predictor of which dark matter particles will eventually become unbound from halos. We also study the mass fraction of unbound particles, which increases strongly towards the edges of halos, and decreases significantly at higher redshifts. We discuss implications for dark matter detection experiments, precision calibrations of the halo mass function, the use of baryon fractions to constrain dark energy, and searches for intergalactic supernovae.

We develop empirical methods for modeling the galaxy population and populating cosmological N-body simulations with mock galaxies according to the observed properties of galaxies in survey data. We use these techniques to produce a new set of mock catalogs for the DEEP2 Galaxy Redshift Survey based on the output of the high-resolution Bolshoi simulation, as well as two other simulations with different cosmological parameters, all of which we release for public use. The mock-catalog creation technique uses subhalo abundance matching to assign galaxy luminosities to simulated dark-matter halos. It then adds color information to the resulting mock galaxies in a manner that depends on the local galaxy density, in order to reproduce the measured color-environment relation in the data. In the course of constructing the catalogs, we test various models for including scatter in the relation between halo mass and galaxy luminosity, within the abundance-matching framework. We find that there is no constant-scatter model that can simultaneously reproduce both the luminosity function and the autocorrelation function of DEEP2. This result has implications for galaxy-formation theory, and it restricts the range of contexts in which the mocks can be usefully applied. Nevertheless, careful comparisons show that our new mocks accurately reproduce a wide range of the other properties of the DEEP2 catalog, suggesting that they can be used to gain a detailed understanding of various selection effects in DEEP2.

We present a robust method to constrain average galaxy star formation rates, star formation histories, and the intracluster light as a function of halo mass. Our results are consistent with observed galaxy stellar mass functions, specific star formation rates, and cosmic star formation rates from z=0 to z=8. We consider the effects of a wide range of uncertainties on our results, including those affecting stellar masses, star formation rates, and the halo mass function at the heart of our analysis. As they are relevant to our method, we also present new calibrations of the dark matter halo mass function, halo mass accretion histories, and halo-subhalo merger rates out to z=8. We also provide new compilations of cosmic and specific star formation rates; more recent measurements are now consistent with the buildup of the cosmic stellar mass density at all redshifts. Implications of our work include: halos near 10^12 Msun are the most efficient at forming stars at all redshifts, the baryon conversion efficiency of massive halos drops markedly after z ~ 2.5 (consistent with theories of cold-mode accretion), the ICL for massive galaxies is expected to be significant out to at least z ~ 1-1.5, and dwarf galaxies at low redshifts have higher stellar mass to halo mass ratios than previous expectations and form later than in most theoretical models. Finally, we provide new fitting formulae for star formation histories that are more accurate than the standard declining tau model. Our approach places a wide variety of observations relating to the star formation history of galaxies into a self-consistent framework based on the modern understanding of structure formation in LCDM. Constraints on the stellar mass-halo mass relationship and star formation rates are available for download at http://www.peterbehroozi.com/data.html .

While various codes exist to systematically and robustly find haloes and subhaloes in cosmological simulations (Knebe et al., 2011, Onions et al., 2012), this is the first work to introduce and rigorously test codes that find tidal debris (streams and other unbound substructure) in fully cosmological simulations of structure formation. We use one tracking and three non-tracking codes to identify substructure (bound and unbound) in a Milky Way type simulation from the Aquarius suite (Springel et al., 2008) and post-process their output with a common pipeline to determine the properties of these substructures in a uniform way. By using output from a fully cosmological simulation, we also take a step beyond previous studies of tidal debris that have used simple toy models. We find that both tracking and non-tracking codes agree well on the identification of subhaloes and more importantly, the {\em unbound tidal features} associated with them. The distributions of basic properties of the total substructure distribution (mass, velocity dispersion, position) are recovered with a scatter of $\sim20%$. Using the tracking code as our reference, we show that the non-tracking codes identify complex tidal debris with purities of $\sim40%$. Analysing the results of the substructure finders, we find that the general distribution of {\em substructures} differ significantly from the distribution of bound {\em subhaloes}. Most importantly, both bound and unbound {\em substructures} together constitute $\sim18%$ of the host halo mass, which is a factor of $\sim2$ higher than the fraction in self-bound {\em subhaloes}. However, this result is restricted by the remaining challenge to cleanly define when an unbound structure has become part of the host halo. Nevertheless, the more general substructure distribution provides a more complete picture of a halo's accretion history.

We discuss the properties of subhalos in cluster-size halos, using a high-resolution statistical sample: the Rhapsody simulations introduced in Wu et al. (2012). We demonstrate that the criteria applied to select subhalos have significant impact on the inferred properties of the sample, including the scatter in the number of subhalos, the correlation between the subhalo number and formation time, and the shape of subhalos' spatial distribution and velocity structure. We find that the number of subhalos, when selected using the peak maximum circular velocity in their histories (a property expected to be closely related to the galaxy luminosity), is uncorrelated with the formation time of the main halo. This is in contrast to the previously reported correlation from studies where subhalos are selected by the current maximum circular velocity; we show that this difference is a result of the tidal stripping of the subhalos. We also find that the dominance of the main halo and the subhalo mass fraction are strongly correlated with halo concentration and formation history. These correlations are important to take into account when interpreting results from cluster samples selected with different criteria. Our sample also includes a fossil cluster, which is presented separately and placed in the context of the rest of the sample.

We provide new constraints on the connection between galaxies in the local universe, identified by the Sloan Digital Sky Survey (SDSS), and dark matter halos and their constituent substructures in the $\Lambda$CDM model using WMAP7 cosmological parameters. Predictions for the abundance and clustering properties of dark matter halos, and the relationship between dark matter hosts and substructures, are based on a high-resolution cosmological simulation, the Bolshoi simulation. We associate galaxies with halos and subhalos using subhalo abundance matching, performing a comprehensive analysis which investigates the underlying assumptions of this technique including (a) which halo property is most closely associated with galaxy stellar masses and luminosities, (b) how much scatter is in this relationship, and (c) how much subhalos can be stripped before their galaxies are destroyed. The models are jointly constrained by new measurements of the projected two-point galaxy clustering and the observed conditional stellar mass function of galaxies in groups. The data put tight constraints on the satellite fraction of galaxies as a function of galaxy stellar mass, on the scatter between halo and galaxy properties, and on the underlying conditional stellar mass function. These data rule out several halo properties commonly used in abundance matching, largely because the satellite fractions in the models disagree with those data. We show that an abundance matching model that associates galaxies with the peak circular velocity of their halos is in good agreement with the data, when scatter of $0.20 \pm 0.03$ dex in stellar mass at a given peak velocity is included. This will yield important constraints for galaxy formation models, and also provides encouraging indications that the galaxy--halo connection can be modeled with sufficient fidelity for future precision studies of the dark Universe.

Physical processes regulating star formation in satellite galaxies represent an area of ongoing research, but the projected nature of observed coordinates makes separating different populations of satellites (with different processes at work) difficult. The orbital history of a satellite galaxy leads to its present-day phase space coordinates; we can also work backwards and use these coordinates to statistically infer information about the orbital history. We use merger trees from the MultiDark Run 1 N-body simulation to compile a catalog of the orbits of satellite haloes in cluster environments. We parameterize the orbital history by the time since crossing within 2.5 rvir of the cluster centre and use our catalog to estimate the probability density over a range of this parameter given a set of present-day projected (i.e. observable) phase space coordinates. We show that different populations of satellite haloes, e.g. infalling, backsplash and virialized, occupy distinct regions of phase space, and semi-distinct regions of projected phase space. This will allow us to probabilistically determine the time since infall of a large sample of observed satellite galaxies, and ultimately to study the effect of orbital history on star formation history (the topic of a future paper). We test the accuracy of our method and find that we can reliably recover this time within +/-2.58 Gyr in 68 per cent of cases by using all available phase space coordinate information, compared to +/-2.64 Gyr using only position coordinates and +/-3.10 Gyr guessing 'blindly', i.e. using no coordinate information, but with knowledge of the overall distribution of infall times. In some regions of phase space, the accuracy of the infall time estimate improves to +/-1.85 Gyr. Although we focus on time since infall, our method is easily generalizable to other orbital parameters (e.g. pericentric distance and time).

We present a new algorithm for generating merger trees and halo catalogs which explicitly ensures consistency of halo properties (mass, position, and velocity) across timesteps. 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/ .

We present a new algorithm for identifying dark matter halos, substructure, and tidal features. The approach is based on adaptive hierarchical refinement of friends-of-friends groups in six phase-space dimensions and one time dimension, which allows for robust (grid-independent, shape-independent, and noise-resilient) tracking of substructure; as such, it is named Rockstar (Robust Overdensity Calculation using K-Space Topologically Adaptive Refinement). Our method is massively parallel (up to 10^5 CPUs) and runs on the largest current simulations (>10^10 particles) with high efficiency (10 CPU hours and 60 gigabytes of memory required per billion particles analyzed). A previous paper (Knebe et al 2011) has shown Rockstar to have class-leading recovery of halo properties; we expand on these comparisons with more tests and higher-resolution simulations. We show a significant improvement in substructure recovery as compared to several other halo finders and discuss the theoretical and practical limits of simulations in this regard. Finally, we present results which demonstrate conclusively that dark matter halo cores are not at rest relative to the halo bulk or satellite average velocities and have coherent velocity offsets across a wide range of halo masses and redshifts. For massive clusters, these offsets can be up to 350 km/s at z=0 and even higher at high redshifts. Our implementation is publicly available at http://code.google.com/p/rockstar .

We present the first results from the Rhapsody cluster re-simulation project: a sample of 96 "zoom-in" simulations of dark matter halos of 10^14.8 +- 0.05 Msun/h, selected from a 1 (Gpc/h)^3 volume. This simulation suite is the first to resolve this many halos with ~5x10^6 particles per halo in the cluster-mass regime, allowing us to statistically characterize the distribution of and correlation between halo properties at fixed mass. We focus on the properties of the main halos and how they are affected by formation history, which we track back to z=12, over five decades in mass. We give particular attention to the impact of the formation history on the density profiles of the halos. We find that the deviations from the Navarro-Frenk-White (NFW) model and the Einasto model depend on formation time. Late-forming halos tend to have considerable deviations from both models, partly due to the presence of massive subhalos, while early-forming halos deviate less but still significantly from the NFW model and are better described by the Einasto model. We find that the halo shapes depend only moderately on formation time. Departure from spherical symmetry impacts the density profiles through the anisotropic distribution of massive subhalos. Further evidence of the impact of subhalos is provided by analyzing the phase-space structure. A detailed analysis of the properties of the subhalo population in Rhapsody is presented in a companion paper.

Using reconstructed galaxy star formation histories, we calculate the instantaneous efficiency of galaxy star formation (i.e., the star formation rate divided by the baryon accretion rate) from $z=8$ to the present day. This efficiency exhibits a clear peak near a characteristic halo mass of 10^11.7 Msun, which coincides with longstanding theoretical predictions for the mass scale relevant to virial shock heating of accreted gas. Above the characteristic halo mass, the efficiency falls off as the mass to the minus four-thirds power; below the characteristic mass, the efficiency falls off at an average scaling of mass to the two-thirds power. By comparison, the shape and normalization of the efficiency change very little since z=4. We show that a time-independent star formation efficiency simply explains the shape of the cosmic star formation rate since z=4 in terms of dark matter accretion rates. The rise in the cosmic star formation from early times until z=2 is especially sensitive to galaxy formation efficiency. The mass dependence of the efficiency strongly limits where most star formation occurs, with the result that two-thirds of all star formation has occurred inside halos within a factor of three of the characteristic mass, a range that includes the mass of the Milky Way.

We present a study of a comparison of spin distributions of subhaloes found associated with a host halo. The subhaloes are found within two cosmological simulation families of Milky Way-like galaxies, namely the Aquarius and GHALO simulations. These two simulations use different gravity codes and cosmologies. We employ ten different substructure finders, which span a wide range of methodologies from simple overdensity in configuration space to full 6-d phase space analysis of particles.We subject the results to a common post-processing pipeline to analyse the results in a consistent manner, recovering the dimensionless spin parameter. We find that spin distribution is an excellent indicator of how well the removal of background particles (unbinding) has been carried out. We also find that the spin distribution decreases for substructure the nearer they are to the host halo's, and that the value of the spin parameter rises with enclosed mass towards the edge of the substructure. Finally subhaloes are less rotationally supported than field haloes, with the peak of the spin distribution having a lower spin parameter.

We present a detailed comparison of the substructure properties of a single Milky Way sized dark matter halo from the Aquarius suite at five different resolutions, as identified by a variety of different (sub-)halo finders for simulations of cosmic structure formation. These finders span a wide range of techniques and methodologies to extract and quantify substructures within a larger non-homogeneous background density (e.g. a host halo). This includes real-space, phase-space, velocity-space and time- space based finders, as well as finders employing a Voronoi tessellation, friends-of-friends techniques, or refined meshes as the starting point for locating substructure.A common post-processing pipeline was used to uniformly analyse the particle lists provided by each finder. We extract quantitative and comparable measures for the subhaloes, primarily focusing on mass and the peak of the rotation curve for this particular study. We find that all of the finders agree extremely well on the presence and location of substructure and even for properties relating to the inner part part of the subhalo (e.g. the maximum value of the rotation curve). For properties that rely on particles near the outer edge of the subhalo the agreement is at around the 20 per cent level. We find that basic properties (mass, maximum circular velocity) of a subhalo can be reliably recovered if the subhalo contains more than 100 particles although its presence can be reliably inferred for a lower particle number limit of 20. We finally note that the logarithmic slope of the subhalo cumulative number count is remarkably consistent and <1 for all the finders that reached high resolution. If correct, this would indicate that the larger and more massive, respectively, substructures are the most dynamically interesting and that higher levels of the (sub-)subhalo hierarchy become progressively less important.

We calculate the probability that a Milky-Way-like halo in the standard cosmological model has the observed number of Magellanic Clouds (MCs). The statistics of the number of MCs in the LCDM model are in good agreement with observations of a large sample of SDSS galaxies. Under the sub-halo abundance matching assumption of a relationship with small scatter between galaxy r-band luminosities and halo internal velocities v_max, we make detailed comparisons to similar measurements using SDSS DR7 data by Liu et al. (2010). Models and observational data give very similar probabilities for having zero, one, and two MC-like satellites. In both cases, Milky Way-luminosity hosts have just a \sim 10% chance of hosting two satellites similar to the Magellanic Clouds. In addition, we present a prediction for the probability for a host galaxy to have Nsats satellite galaxies as a function of the magnitudes of both the host and satellite. This probability and its scaling with host properties is significantly different from that of mass-selected objects because of scatter in the mass- luminosity relation and because of variations in the star formation efficiency with halo mass.

Understanding the mechanisms that lead dense environments to host galaxies with redder colors, more spheroidal morphologies, and lower star formation rates than field populations remains an important problem. As most candidate processes ultimately depend on host halo mass, accurate characterizations of the local environment, ideally tied to halo mass estimates and spanning a range in halo mass and redshift are needed. In this work, we present and test a rigorous, probabalistic method for assigning galaxies to groups based on precise photometric redshifts and X-ray selected groups drawn from the COSMOS field. The groups have masses in the range 10^13 < M_200c/M_sun < 10^14 and span redshifts 0<z<1. We characterize our selection algorithm via tests on spectroscopic subsamples, including new data obtained at the VLT, and by applying our method to detailed mock catalogs. We find that our group member galaxy sample has a purity of 84% and completeness of 92% within 0.5 R200c. We measure the impact of uncertainties in redshifts and group centering on the quality of the member selection with simulations based on current data as well as future imaging and spectroscopic surveys. As a first application of our new group member catalog which will be made publicly available, we show that member galaxies exhibit a higher quenched fraction compared to the field at fixed stellar mass out to z~1, indicating a significant relationship between star formation and environment at group scales. We also address the suggestion that dusty star forming galaxies in such groups may impact the high-l power spectrum of the cosmic microwave background and find that such a population cannot explain the low power seen in recent SZ measurements.

Measurements of the total amount of stars locked up in galaxies as a function of host halo mass contain key clues about the efficiency of processes that regulate star formation. We derive the total stellar mass fraction f_star as a function of halo mass M500c from z=0.2 to z=1 using two complementary methods. First, we derive f_star using a statistical Halo Occupation Distribution model jointly constrained by data from lensing, clustering, and the stellar mass function. This method enables us to probe f_star over a much wider halo mass range than with group or cluster catalogs. Second, we derive f_star at group scales using a COSMOS X-ray group catalog and we show that the two methods agree to within 30%. We quantify the systematic uncertainty on f_star using abundance matching methods and we show that the statistical uncertainty on f_star (~10%) is dwarfed by systematic uncertainties associated with stellar mass measurements (~45% excluding IMF uncertainties). Assuming a Chabrier IMF, we find 0.012<f_star<0.025 at M500c=10^13 Msun and 0.0057<f_star<0.015 at M500c=10^14 Msun. These values are significantly lower than previously published estimates. We investigate the cause of this difference and find that previous work has overestimated f_star due to a combination of inaccurate stellar mass estimators and/or because they have assumed that all galaxies in groups are early type galaxies with a constant mass-to-light ratio. Contrary to previous claims, our results suggest that the mean value of f_star is always significantly lower than f_gas for halos above 10^13 Msun. Combining our results with recently published gas mas fractions, we find a shortfall in f_star+f_gas at R500c compared to the cosmic mean. This shortfall varies with halo mass and becomes larger towards lower halos masses.

We develop a theoretical framework that combines measurements of galaxy-galaxy lensing, galaxy clustering, and the galaxy stellar mass function in a self-consistent manner. While considerable effort has been invested in exploring each of these probes individually, attempts to combine them are still in their infancy despite the potential of such combinations to elucidate the galaxy-dark matter connection, to constrain cosmological parameters, and to test the nature of gravity. In this paper, we focus on a theoretical model that describes the galaxy-dark matter connection based on standard halo occupation distribution techniques. Several key modifications enable us to extract additional parameters that determine the stellar-to-halo mass relation and to simultaneously fit data from multiple probes while allowing for independent binning schemes for each probe. In a companion paper, we demonstrate that the model presented here provides an excellent fit to galaxy-galaxy lensing, galaxy clustering, and stellar mass functions measured in the COSMOS survey from z=0.2 to z=1.0. We construct mock catalogs from numerical simulations to investigate the effects of sample variance and covariance on each of the three probes. Finally, we analyze and discuss how trends in each of the three observables impact the derived parameters of the model. In particular, we investigate the various features of the observed galaxy stellar mass function (low-mass slope, plateau, knee, and high-mass cut-off) and show how each feature is related to the underlying relationship between stellar and halo mass. We demonstrate that the observed plateau feature in the stellar mass function at Mstellar~2x10^10 Msun is due to the transition that occurs in the stellar-to-halo mass relation at Mhalo ~ 10^12 Msun from a low-mass power-law regime to a sub-exponential function at higher stellar mass.

We introduce a probabilistic approach to the problem of counting dwarf satellites around host galaxies in databases with limited redshift information. This technique is used to investigate the occurrence of satellites with luminosities similar to the Magellanic Clouds around hosts with properties similar to the Milky Way in the object catalog of the Sloan Digital Sky Survey. Our analysis uses data from SDSS Data Release 7, selecting candidate Milky-Way-like hosts from the spectroscopic catalog and candidate analogs of the Magellanic Clouds from the photometric catalog. Our principal result is the probability for a Milky-Way-like galaxy to host N_{sat} close satellites with luminosities similar to the Magellanic Clouds. We find that 81 percent of galaxies like the Milky Way are have no such satellites within a radius of 150 kpc, 11 percent have one, and only 3.5 percent of hosts have two. The probabilities are robust to changes in host and satellite selection criteria, background-estimation technique, and survey depth. These results demonstrate that the Milky Way has significantly more satellites than a typical galaxy of its luminosity; this fact is useful for understanding the larger cosmological context of our home galaxy.

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).

Using data from the COSMOS survey, we perform the first joint analysis of galaxy-galaxy weak lensing, galaxy spatial clustering, and galaxy number densities. Carefully accounting for sample variance and for scatter between stellar and halo mass, we model all three observables simultaneously using a novel and self-consistent theoretical framework. Our results provide strong constraints on the shape and redshift evolution of the stellar-to-halo mass relation (SHMR) from z=0.2 to z=1. At low stellar mass, we find that halo mass scales as Mh M*^0.46 and that this scaling does not evolve significantly with redshift to z=1. We show that the dark-to-stellar ratio, Mh/M*, varies from low to high masses, reaching a minimum of Mh/M*~27 at M*=4.5x10^10 Msun and Mh=1.2x10^12 Msun. This minimum is important for models of galaxy formation because it marks the mass at which the accumulated stellar growth of the central galaxy has been the most efficient. We describe the SHMR at this minimum in terms of the "pivot stellar mass", M*piv, the "pivot halo mass", Mhpiv, and the "pivot ratio", (Mh/M*)piv. Thanks to a homogeneous analysis of a single data set, we report the first detection of mass downsizing trends for both Mhpiv and M*piv. The pivot stellar mass decreases from M*piv=5.75+-0.13x10^10 Msun at z=0.88 to M*piv=3.55+-0.17x10^10 Msun at z=0.37. Intriguingly, however, the corresponding evolution of Mhpiv leaves the pivot ratio constant with redshift at (Mh/M*)piv~27. We use simple arguments to show how this result raises the possibility that star formation quenching may ultimately depend on Mh/M* and not simply Mh, as is commonly assumed. We show that simple models with such a dependence naturally lead to downsizing in the sites of star formation. Finally, we discuss the implications of our results in the context of popular quenching models, including disk instabilities and AGN feedback.

We conduct a comprehensive analysis of the relationship between central galaxies and their host dark matter halos, as characterized by the stellar mass-halo mass (SM-HM) relation, with rigorous consideration of uncertainties. Our analysis focuses on results from the abundance matching technique, which assumes that every dark matter halo or subhalo above a specific mass threshold hosts one galaxy. We discuss the quantitative effects of uncertainties in observed galaxy stellar mass functions (GSMFs) (including stellar mass estimates and counting uncertainties), halo mass functions (including cosmology and uncertainties from substructure), and the abundance matching technique used to link galaxies to halos (including scatter in this connection). Our analysis results in a robust estimate of the SM-HM relation and its evolution from z=0 to z=4. The shape and evolution are well constrained for z < 1. The largest uncertainties at these redshifts are due to stellar mass estimates; however, failure to account for scatter in stellar masses at fixed halo mass can lead to errors of similar magnitude in the SM-HM relation for central galaxies in massive halos. We also investigate the SM-HM relation to z=4, although the shape of the relation at higher redshifts remains fairly unconstrained when uncertainties are taken into account. These results will provide a powerful tool to inform galaxy evolution models. [Abridged]